05.04.2016 Views

Modern Engineering Thermodynamics

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

632 CHAPTER 15: Chemical <strong>Thermodynamics</strong><br />

EXAMPLE 15.15 (Continued)<br />

a. To determine the required results for part a of this problem, we set p m = p° =0:1 MPa, and then K e can be found in<br />

Table C.17 for various temperatures. Our task is now to pick several reaction temperatures, look up their corresponding<br />

K e values, and solve the previous equation for y.<br />

Squaring both sides of the preceding reversible equilibrium dissociation reaction equation and rearranging gives a<br />

cubic equation in y:<br />

where<br />

y 3 +<br />

3α<br />

1 − α y − 2α<br />

1 − α = 0<br />

<br />

α = Ke<br />

2 p°<br />

p m<br />

A cubic equation in y has a simple algebraic solution. If we let<br />

a =<br />

3α<br />

1 − α<br />

and b = −<br />

2α<br />

1 − α<br />

then the cubic equation becomes y 3 + ay + b = 0, and the three roots of this equation are<br />

where<br />

y 1<br />

y 2<br />

y 3<br />

= A + B<br />

= − A + B + A − B pffiffiffiffiffiffiffiffi<br />

− 3<br />

2 2<br />

= − A + B − A − B pffiffiffiffiffiffiffiffi<br />

− 3<br />

2 2<br />

r ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi! 1/3 r ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi! 1/3<br />

A = − b 2 + b 2<br />

4 + a3 and B = − b 27<br />

2 − b 2<br />

4 + a3<br />

27<br />

Note that, in evaluating these expressions, the cube root of a negative number is computed as the negative of the cube<br />

root of the positive value of the number; that is, (−8) 1/3 = −(8 1/3 ) = −2.<br />

If α is less than 1, then b 2 /4 + a 3 /27 > 0 and there is just one real root, given by y 1 .Butifα is greater than 1, then<br />

b 2 /4 + a 3 /27 < 0 and there are three real roots. The following procedure is a simple way of calculating these roots. Let<br />

J = [(− (a/3) 3 ] 1/2 and K = J 1/3 . Then calculate L = arccos[−b/(2J)], M = cos(L/3), and N = (3 1/2 )sin(L/3). The three roots are<br />

then given by<br />

y 1 = 2K cosðL/3Þ<br />

y 2 = − KðM + NÞ<br />

y 3 = − KðM − NÞ<br />

Note that the only valid root here is the one between 0 and 1. Table 15.8 lists several typical results using these solution<br />

equations with p°=0.1 MPa.<br />

These results can then be plotted either by hand or with the use of computer software. The final plot showing<br />

the variation in the degree of dissociation (y) versus reaction equilibrium temperature is shown in Figure 15.10.<br />

b. With the solution from part a, it is a simple matter to generate reaction equations for various total pressures at T = 3000. K.<br />

The curve in Figure 15.11 shows the results of this effort and provides a graphical representation of the effect of total<br />

pressure on the degree of dissociation (y) at 3000. K.<br />

Table 15.8 Typical Values for the Degree of Dissociation of Carbon Dioxide at Various<br />

Temperatures and Pressures with p°=0.100 MPa<br />

Reaction<br />

Temperature (K)<br />

Reaction Pressure<br />

p m (MPa) log 10 (K e ) K e y<br />

2000. K 0.100 MPa 2.863 0.0014 0.0154<br />

3000. K 0.100 MPa 0.469 0.3396 0.4436<br />

5000. K 0.100 MPa +1.387 24.378 0.9770<br />

3000. K 1.00 MPa 0.469 0.3396 0.2453<br />

3000. K 10.0 MPa 0.469 0.3396 0.1235

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!