Solution Koretsky Completo PDF
Solution Koretsky Completo PDF
Solution Koretsky Completo PDF
Wyatt Tenhaeff
Milo Koretsky
[email protected]
1.2
An approximate solution can be found if we combine Equations 1.4 and 1.5:
1 2
mV ekmolecular
2
3
kT ekmolecular
2
3kT
V
m
The molecules are traveling really, fast (around the length of five football fields every second).
Comment:
We can get a better solution by using the Maxwell-Boltzmann distribution of speeds that is
sketched in Figure 1.4. Looking up the quantitative expression for this expression, we have:
3/ 2
m m 2 2
f (v)dv 4 exp v v dv
2kT 2kT
where f(v) is the fraction of molecules within dv of the speed v. We can find the average speed
by integrating the expression above
f (v)vdv
449 m/s
8kT
V 0
m
f (v)dv
0
2
1.3
Derive the following expressions by combining Equations 1.4 and 1.5:
3kT 3kT
Va2 Vb2
ma mb
Therefore,
Va2 mb
Vb2 ma
Since mb is larger than ma, the molecules of species A move faster on average.
3
1.4
We have the following two points that relate the Reamur temperature scale to the Celsius scale:
At 22 C,
T 17.6 Reamur
4
1.5
(a)
After a short time, the temperature gradient in the copper block is changing (unsteady state), so
the system is not in equilibrium.
(b)
After a long time, the temperature gradient in the copper block will become constant (steady
state), but because the temperature is not uniform everywhere, the system is not in equilibrium.
(c)
After a very long time, the temperature of the reservoirs will equilibrate; The system is then
homogenous in temperature. The system is in thermal equilibrium.
5
1.6
We assume the temperature is constant at 0 C. The molecular weight of air is
0.029 kg/mol9.81 m/s8848 m
P 1 atm exp
J
8.314
273.15 K
mol K
P 0.330 atm 33.4 kPa
Therefore, the liquid boils at 71.4 C. Note: the barometric relationship given assumes that the
temperature remains constant. In reality the temperature decreases with height as we go up the
mountain. However, a solution in which T and P vary with height is not as straight-forward.
6
1.7
To solve these problems, the steam tables were used. The values given for each part constrain
the water to a certain state. In most cases we can look at the saturated table, to determine the
state.
7
1.8
From the steam tables in Appendix B.1:
m3
v critical 0.003155 T 374.15 C, P 22.089 MPa
kg
m3
vlsat 0.001127
kg
m3
vvsat 0.19444
kg
Because the total mass and volume of the closed, rigid system remain constant as the water
condenses, we can develop the following expression:
where x is the quality of the water. Substituting values and solving for the quality, we obtain
x 0.0105 or 1.05 %
8
1.9
The calculation methods will be shown for part (a), but not parts (b) and (c)
(a)
Use the following equation to estimate the specific volume:
v1.9 MPa, 250 C v1.8 MPa, 250 C 0.5v2.0 MPa, 250 C v1.8 MPa, 250 C
Substituting data from the steam tables,
m3
v1.9 MPa, 250 C 0.11821
kg
v v NIST
% Difference 100 % 0.254 %
v NIST
(b)
Linear interpolation:
m3
v1.9 MPa, 300 C 0.13284
kg
NIST website:
m3
v NIST 1.9 MPa, 300 C 0.13249
kg
Therefore,
% Difference 0.264 %
(c)
Linear interpolation:
m3
v1.9 MPa, 270 C 0.12406
kg
9
Note: Double interpolation is required to determine this value. First, find the molar volumes at
270 C and 1.8 MPa and 2.0 MPa using interpolation. Then, interpolate between the results from
the previous step to find the molar volume at 270 C and 1.9 MPa.
NIST website:
m3
v NIST 1.9 MPa, 270 C 0.12389
kg
Therefore,
% Difference 0.137 %
With regards to parts (a), (b), and (c), the values found using interpolation and the NIST website
agree very well. The discrepancies will not significantly affect the accuracy of any subsequent
calculations.
10
1.10
For saturated temperature data at 25 C in the steam tables,
m3 L
vlsat 0.001003 1.003
kg kg
V 1 L
m 0.997 kg
vlsat L
1.003
kg
Because molar volumes of liquids do not depend strongly on pressure, the mass of water at 25 C
and atmospheric pressure in a one liter should approximately be equal to the mass calculated
above unless the pressure is very, very large.
11
1.11
First, find the overall specific volume of the water in the container:
1 m3 m3
v 0.2
5 kg kg
T T sat 120.23 C
Let ml represent the mass of the water in the container that is liquid, and mv represent the mass
of the water in the container that is gas. These two masses are constrained as follows:
ml mv 5 kg
We also have the extensive volume of the system equal the extensive volume of each phase
ml vlsat mv vvsat V
ml 0.001061 m 3 / kg mv 0.8857 m 3 /kg 1 m 3
Solving these two equations simultaneously, we obtain
ml 3.88 kg
mv 1.12 kg
mv 1.12 kg
x 0.224
m 5 kg
The internal energy relative to the reference state in the saturated steam table is
U ml ulsat mv uvsat
12
ulsat 504.47 kJ/kg
uvsat 2529.5 kJ/kg
Therefore,
13
1.12
First, determine the total mass of water in the container. Since we know that 10 % of the mass is
vapor, we can write the following expression
V m 0.9 vlsat 0.1vvsat
From the saturated steam tables for P sat 1 MPa
sat
m3
vl 0.001127
kg
m3
vvsat 0.1944
kg
Therefore,
V 1 m3
m
0.9vlsat 0.1vvsat m3 m3
0.9 0.001127 0.1 0.1944
kg
kg
m 48.9 kg
and
Vv 0.1mvvsat 0.148.9 kg 0.1944 m 3 / kg
Vv 0.950 m 3
14
1.13
In a spreadsheet, make two columns: one for the specific volume of water and another for the
pressure. First, copy the specific volumes of liquid water from the steam tables and the
corresponding saturation pressures. After you have finished tabulating pressure/volume data for
liquid water, list the specific volumes and saturation pressures of water vapor. Every data point
is not required, but be sure to include extra points near the critical values. The data when plotted
on a logarithmic scale should look like the following plot.
100.0000
Critical Point
10.0000
1.0000
0.0001 0.001 0.01 0.1 1 10 100 1000
0.1000
0.0100
0.0010
Specific Volume, v (m 3/kg)
15
1.14
The ideal gas law can be rewritten as
RT
v
P
For each part in the problem, the appropriate values are substituted into this equation where
L bar
R 0.08314
mol K
is used. The values are then converted to the units used in the steam table.
(a)
L
v 30.7
mol
L
30.7
mol 10 3 m 1.70 m
3 3
v
v L kg
MW kg
0.018
mol
For part (a), the calculation of the percent error will be demonstrated. The percent error will be
based on percent error from steam table data, which should be more accurate than using the ideal
gas law. The following equation is used:
v vST
%Error IG 100 %
vST
where v IG is the value calculated using the ideal gas law and vST is the value from the steam
table.
1.70 m 1.6729 m
3 3
kg kg
%Error % 1.62%
m3
1.6729
kg
This result suggests that you would introduce an error of around 1.6% if you characterize boiling
water at atmospheric pressure as an ideal gas.
16
(b)
L
v 64.3
mol
m3
v 3.57
kg
% Error 0.126%
For a given pressure, when the temperature is raised the gas behaves more like an ideal gas.
(c)
L
v 0.643
mol
m3
v 0.0357
kg
%Error = 8.87%
(d)
L
v 1.06
mol
m3
v 0.0588
kg
%Error = 0.823%
17
1.15
The room I am sitting in right now is approximately 16 ft long by 12 feet wide by 10 feet tall
your answer should vary. The volume of the room is
V 12 ft 16 ft 10 ft 1920 ft 3 54.4 m 3
The room is at atmospheric pressure and a temperature of 22 C. Calculate the number of moles
using the ideal gas law:
n
PV
1.01325 10 5 Pa 54.4 m 3
RT J
8.314 295.15 K
mol K
n 2246 mol
18
1.16
First, find the total volume that one mole of gas occupies. Use the ideal gas law:
J
8.314 293.15 K
RT mol K m3
v
P
1105 Pa 0.0244
mol
v occupied
Percentage 100 % 0.0349 %
v
19
1.17
Water condenses on your walls when the water is in liquid-vapor equilibrium. To find the
maximum allowable density at 70 F, we need to find the smallest density of water vapor at 40 F
which satisfies equilibrium conditions. In other words, we must find the saturation density of
water vapor at 40 F. From the steam tables,
m3
vvsat 153.74 (sat. water vapor at 40 F = 4.44 C)
kg
We must correct the specific volume for the day-time temperature of 70 F (21.1 C ) with the
ideal gas law.
vvsat
,277.6 K v
294.3 K
T277.6 T294.3
T
v294.3 K 294.3 vvsat
294.3 K
277.6 K
153 . 74
m3
163
m3
kg
, 277.6 K
T277.6
kg
Therefore,
1 kg
6.13 10 3
vvsat m3
If the density at 70 F is greater than or equal to 6.13 10 3 kg/m 3 , the density at night when the
temperature is 40 F will be greater than the saturation density, so water will condense onto the
wall. However, if the density is less than 6.13 10 3 kg/m 3 , then saturation conditions will not
be obtained, and water will not condense onto the walls.
20
1.18
We consider the air inside the soccer ball as the system. We can answer this question by looking
at the ideal gas law:
Pv RT
If we assume it is a closed system, it will have the same number of moles of air in the winter as
the summer. However, it is colder in the winter (T is lower), so the ideal gas law tells us that Pv
will be lower and the balls will be under inflated.
Alternatively, we can argue that the higher pressure inside the ball causes air to leak out over
time. Thus we have an open system and the number of moles decrease with time leading to the
under inflation.
21
1.19
First, write an equation for the volume of the system in its initial state:
ml vlsat mv vvsat V
Substitute ml 1 x m and mv xm :
m 1 x vlsat xvvsat V
mv critical V
m3
v critical 0.003155 T 374.15 C, P 22.089 MPa
kg
m3
vlsat 0.001043
kg
sat
m3
vv 1.6940
kg
x 0.00125 or 0.125 %
22
1.20
As defined in the problem statement, the relative humidity can be calculated as follows
mass of water
Relative Humidity
mass of water at saturation
We can obtain the saturation pressure at each temperature from the steam tables. At 10 [oC], the
saturation pressure of water is 1.22 [kPa]. This value is proportional to the mass of water in the
vapor at saturation. For 90% relative humidity, the partial pressure of water in the vapor is:
This value represents the vapor pressure of water in the air. At 30 [oC], the saturation pressure of
water is 4.25 [kPa]. For 590% relative humidity, the partial pressure of water in the vapor is:
Since the total pressure in each case is the same (atmospheric), the partial pressure is
proportional to the mass of water in the vapor. We conclude there is about twice the amount of
water in the air in the latter case.
23
1.21
(a)
When you have extensive variables, you do not need to know how many moles of each substance
are present. The volumes can simply be summed.
V1 Va Vb
(b)
The molar volume, v1, is equal to the total volume divided by the total number of moles.
V1 Va Vb
v1
ntot na nb
We can rewrite the above equation to include molar volumes for species a and b.
n v nb vb na nb
v1 a a va vb
na nb na nb na nb
v1 xa va xb vb
(c)
The relationship developed in Part (a) holds true for all extensive variables.
K1 K a K b
(d)
k1 xa k a xb kb
24
Chapter 2 Solutions
Engineering and Chemical Thermodynamics
Wyatt Tenhaeff
Milo Koretsky
[email protected]
2.1
There are many possible solutions to this problem. Assumptions must be made to solve the
problem. One solution is as follows. First, assume that half of a kilogram is absorbed by the
towel when you dry yourself. In other words, let
m H 2 O 0.5 kg
Assume that the pressure is constant at 1.01 bar during the drying process. Performing an energy
balance and neglecting potential and kinetic energy effects reveals
q h
Refer to the development of Equation 2.57 in the text to see how this result is achieved. To find
the minimum energy required for drying the towel, assume that the temperature of the towel
remains constant at T 25 C 298.15 K . In the drying process, the absorbed water is
vaporized into steam. Therefore, the expression for heat is
v l
q hH hH
2O 2O
v
where is hH O
is the specific enthalpy of water vapor at P 1.01 bar and T 298.15 K and
2
h l
H 2O
is the specific enthalpy of liquid water at P 1.01 bar and T 298.15 K . A hypothetical
path must be used to calculate the change in enthalpy. Refer to the diagram below
P
liquid h vapor
1 atm
P = 1 [atm]
h1 h3
sat
P
3.17 kPa
h2
liquid vapor
By adding up each step of the hypothetical path, the expression for heat is
q h1 h2 h3
h l ,sat 25 C h l
H 2O 25 C, 1.01 bar hHv,satO 25 C hHl ,satO 25 C
H 2O
2 2
2
However, the calculation of heat can be simplified by treating the water vapor as an ideal gas,
which is a reasonable assumption at low pressure. The enthalpies of ideal gases depend on
temperature only. Therefore, the enthalpy of the vapor change due to the pressure change is
zero. Furthermore, enthalpy is weakly dependent on pressure in liquids. The leg of the
hypothetical path containing the pressure change of the liquid can be neglected. This leaves
v
q hH
2
O
25 C, 3.17 kPa hHl 2 O 25 C, 3.17 kPa
v, sat kJ
hH O
2547 .2 kg (sat. H2O vapor at 25 C)
2
l , sat kJ
hH O
104.87 kg (sat. H2O liquid at 25 C)
2
kJ
q 2442.3
kg
Therefore,
kJ
Q 0.5 kg 2442.3 1221.2 kJ
kg
To find the efficiency of the drying process, assume the dryer draws 30 A at 208 V and takes 20
minutes (1200 s) to dry the towel. From the definition of electrical work,
Q 1221.2 kJ
100 % 100 % 16.3%
W 7488 kJ
There are a number of ways to improve the drying process. A few are listed below.
Dry the towel outside in the sun.
Use a smaller volume dryer so that less air needs to be heated.
Dry more than one towel at a time since one towel cant absorb all of the available
heat. With more towels, more of the heat will be utilized.
3
2.3
In answering this question, we must distinguish between potential energy and internal energy.
The potential energy of a system is the energy the macroscopic system, as a whole, contains
relative to position. The internal energy represents the energy of the individual atoms and
molecules in the system, which can have contributions from both molecular kinetic energy and
molecular potential energy. Consider the compression of a spring from an initial uncompressed
state as shown below.
uncompressed compressed
x M
State 1 State 2
Since it requires energy to compress the spring, we know that some kind of energy must be
stored within the spring. Since this change in energy can be attributed to a change of the
macroscopic position of the system and is not related to changes on the molecular scale, we
determine the form of energy to be potential energy. In this case, the springs tendency to restore
its original shape is the driving force that is analogous to the gravity for gravitational potential
energy.
This argument can be enhanced by the form of the expression that the increased energy takes. If
we consider the spring as the system, the energy it acquires in a reversible, compression from its
initial uncompressed state may be obtained from an energy balance. Assuming the process is
adiabatic, we obtain:
E Q W W
We have left the energy in terms of the total energy, E. The work can be obtained by integrating
the force over the distance of the compression:
1 2
W F dx kxdx kx
2
Hence:
1 2
E kx
2
We see that the increase in energy depends on macroscopic position through the term x.
It should be noted that there is a school of thought that assigns this increased energy to internal
energy. This approach is all right as long as it is consistently done throughout the energy
balances on systems containing springs.
4
2.4
For the first situation, let the rubber band represent the system. In the second situation, the gas is
the system. If heat transfer, potential and kinetic energy effects are assumed negligible, the
energy balance becomes
U W
Since work must be done on the rubber band to stretch it, the value of the work is positive. From
the energy balance, the change in internal energy is positive, which means that the temperature
of the system rises.
When a gas expands in a piston-cylinder assembly, the system must do work to expand against
the piston and atmosphere. Therefore, the value of work is negative, so the change in internal
energy is negative. Hence, the temperature decreases.
In analogy to the spring in Problem 2.3, it can be argued that some of the work imparted into the
rubber band goes to increase its potential energy; however, a part of it goes into stretching the
polymer molecules which make up the rubber band, and the qualitative argument given above
still is valid.
5
2.5
To explain this phenomenon, you must realize that the water droplet is heated from the bottom.
At sufficiently high temperatures, a portion of the water droplet is instantly vaporized. The
water vapor forms an insulation layer between the skillet and the water droplet. At low
temperatures, the insulating layer of water vapor does not form. The transfer of heat is slower
through a gas than a liquid, so it takes longer for the water to evaporate at higher temperatures.
6
2.6
Apartment Surr.
HOT
System
Fridge
+
W -
If the entire apartment is treated as the system, then only the energy flowing across the apartment
boundaries (apartment walls) is of concern. In other words, the energy flowing into or out of the
refrigerator is not explicitly accounted for in the energy balance because it is within the system.
By neglecting kinetic and potential energy effects, the energy balance becomes
U Q W
The Q term represents the heat from outside passing through the apartments walls. The W term
represents the electrical energy that must be supplied to operate the refrigerator.
To determine whether opening the refrigerator door is a good idea, the energy balance with the
door open should be compared to the energy balance with the door closed. In both situations, Q
is approximately the same. However, the values of W will be different. With the door open,
more electrical energy must be supplied to the refrigerator to compensate for heat loss to the
apartment interior. Therefore,
Wajar Wshut
where the subscript ajar refers the situation where the door is open and the subscript shut
refers to the situation where the door is closed. Since,
Qajar Qshut
U ajar Qajar Wajar U shut Qshut Wshut
Tajar Tshut
7
2.7
The two cases are depicted below.
Surr. Surr.
System System
Q Q
Heater
Heater
W
Lets consider the property changes in your house between the following states. State 1, when
you leave in the morning, and state, the state of your home after you have returned home and
heated it to the same temperature as when you left. Since P and T are identical for states 1 and 2,
the state of the system is the same and U must be zero, so
U Q W 0
or
Q W
where -Q is the total heat that escaped between state 1 and state 2 and W is the total work that
must be delivered to the heater. The case where more heat escapes will require more work and
result in higher energy bills. When the heater is on during the day, the temperature in the system
is greater than when it is left off. Since heat transfer is driven by difference in temperature, the
heat transfer rate is greater, and W will be greater. Hence, it is cheaper to leave the heater off
when you are gone.
8
2.8
The amount of work done at constant pressure can be calculated by applying Equation 2.57
H Q
Hence,
H Q mh
where the specific internal energy is used in anticipation of obtaining data from the steam tables.
The mass can be found from the known volume, as follows:
3
1L 0.001 m
L
1.0 kg
V
m
v m3
0.0010
kg
As in Example 2.2, we use values from the saturated steam tables at the same temperature for
subcooled water at 1 atm. The specific enthalpy is found from values in Appendix B.1:
kJ
kg
kJ
kg
kJ
u ul ,2 at 100 o C ul ,1 at 25 o C 419.02 104.87 314.15
kg
Solve for heat:
kJ
Q mu 1.0 kg 314.05 314.15 kJ
kg
and heat rate:
314.15 kJ
0.52 kW
Q
Q
60 s
t
10 min.
min
9
2.9
(a)
From Steam Tables:
kJ
u1 2967.8 (100 kPa, 400 C)
kg
kJ
u 2 2659.8 (50 kPa, 200 C)
kg
kJ
u u 2 u1 308.0
kg
(b)
From Equations 2.53 and 2.63
T2 T2
u u 2 u1 cv dT c P R dT
T1 T1
c P R( A BT CT 2 DT 2 ET 3 )
T2
u R A BT CT 2 DT 2 ET 3 1 dT
T1
Integrating
B C 1 1 E
u R ( A 1)(T2 T1 ) (T2 2 T12 ) (T23 T13 ) D( ) (T2 4 T14 )
2 3 T2 T1 4
A 3.470
B 1.45 10 3
C 0
D 1.21 10 4
E 0
10
J
R 8.314
mol K
T1 (400 273.15 K) 673.15 K
T2 (200 273.15 K) 473.15 K
provides
J
u 5551
mol
J 1 [mol H 2 O] 1000 g 1 kJ kJ
u 5551
308 . 1 kg
mol 18.0148 [g H 2 O] 1 kg 1000 J
The values in parts (a) and (b) agree very well. The answer from part (a) will serve as the basis
for calculating the percent difference since steam table data should be more accurate.
308 308.1
% Difference 100 % 0.03%
308.0
11
2.10
(a)
Referring to the energy balance for closed systems where kinetic and potential energy are
neglected, Equation 2.30 states
U Q W
(b)
Since internal energy is a function of temperature only for an ideal gas (Equation 2.4) and the
process is isothermal
U 0
P P
W nRT ln 2 n1 RT1 ln 2
P1 P1
n1RT1 P1V1
P
W P1V1 ln 2
P1
5 bar
W 8 105 Pa 2.5 10 3 m 3 ln
8 bar
W 940 J
0 J Q W
Q 940 J
(c)
Since the process is adiabatic
Q0
U W
12
The system must do work on the surroundings to expand. Therefore, the work will be negative
and
U 0
T2
U n cv T 0
T1
T2 T1
T2 will be less than 30 C
13
2.11
(a)
(i).
1 Path B
3
P [bar]
Pa
th
2
A
2
1
(ii).
Since internal energy is a function of temperature only for an ideal gas (Equation 2.4) and the
process is isothermal
u 0
Equation 2.48 states that enthalpy is a function of temperature only for an ideal gas. Therefore,
h 0
Performing an energy balance and neglecting potential and kinetic energy produces
u q w 0
P
W nRT ln 2
P1
or
W P
w RT ln 2
n P1
14
J 1 bar
w 8.314 (88 273.15) K ln
mol K 3 bar
J
w 3299
mol
J
q w 3299
mol
(b)
(i). See path on diagram in part (a)
(ii).
Since the overall process is isothermal and u and h are state functions
u 0
h 0
w PE dv
During the constant volume part of the process, no work is done. The work must be solved for
the constant pressure step. Since it is constant pressure, the above equation simplifies to
w PE dv PE (v2 v1 )
J mol
8.314 (88 273.15) K
RT2 K m3
v2 0 . 030
P2 1 10 5 Pa mol
J mol
8.314 (88 273.15) K
RT1 K m3
v1 0 . 010
P1 3 10 5 Pa mol
Substituting in these values and realizing that PE P1 since the process is isobaric produces
15
m3 m3
w (3 10 5 Pa) 0.030 0.010
mol mol
J
w 6000
mol
Performing an energy balance and neglecting potential and kinetic energy results in
u q w 0
J
q w 6000
mol
16
2.12
First, perform an energy balance. No work is done, and the kinetic and potential energies can be
neglected. The energy balance reduces to
U Q
T2
Q n cv dT
T1
T2
Q n c P dT
T1
since the aluminum is a solid. Using the atomic mass of aluminum we find
5 kg
n 185.3 mol
kg
0.02698
mol
Upon substitution of known values and heat capacity data from Table A.2.3, we get
J
323.15 K
Q 185.3 mol 8.314 2.486 1.49 10 3 T dT
mol K 294.15 K
Q 131.61 kJ
17
2.13
First, start with the energy balance. Potential and kinetic energy effects can be neglected.
Therefore, the energy balance becomes
U Q W
The value of the work will be used to obtain the final temperature. The definition of work
(Equation 2.7) is
V2
W PE dV
V1
Since the piston expands at constant pressure, the above relationship becomes
W PE V2 V1
m3
v1 0.02641 (10 MPa, 400 C)
kg
m3
V1 m1v1 (3 kg) 0.02641
0.07923 m
3
kg
V2 V1
W
PE
0.07923 m3
748740 J
6
2.0 10 Pa
0.4536 m3
V
v2 2
0.4536 m 3
0.1512
m3
m2 3 kg kg
Since v2 and P2 are known, state 2 is constrained. From the steam tables:
3
T2 400 C 20 bar, 0.1512 m
kg
Now U will be evaluated, which is necessary for calculating Q . From the steam tables:
18
kJ 3
u 2 2945.2 20 bar, 0.1512 m
kg
kg
kJ
u1 2832.4 100 bar, 400 C
kg
kJ kJ
U m1 u 2 u1 3 kg 2945.2 2832.4 338.4 kJ
kg kg
Substituting the values of U and W into the energy equation allows calculation of Q
Q U W
Q 338400 [J] 748740 J 1.09 10 6 J
19
2.14
In a reversible process, the system is never out of equilibrium by more than an infinitesimal
amount. In this process the gas is initially at 2 bar, and it expands against a constant pressure of
1 bar. Therefore, a finite mechanical driving force exists, and the process is irreversible.
To solve for the final temperature of the system, the energy balance will be written. The piston-
cylinder assembly is well-insulated, so the process can be assumed adiabatic. Furthermore,
potential and kinetic energy effects can be neglected. The energy balance simplifies to
U W
n1 n2
Let n n1 n2
T2 V2
n cv dT PE dV
T1 V1
ncv T2 T1 PE V2 V1
nRT2
V2
P2
PV
T1 1 1
nR
Substituting these expressions into the energy balance, realizing that PE P2 , and simplifying
the equation gives
5
P2 P1 V1
T2
2
7
nR
2
20
P1 2 bar
P2 1 bar
V1 10 L
n 1.0 mol
L bar
R 0.08314
mol K
results in
T2 206 K
To find the value for work, the energy balance can be used
W U ncv T2 T1
PV
T1 1 1
2 bar 10 L 241 K
1 mol 0.08314 L bar
nR
mol K
W 727 J
21
2.15
The maximum work can be obtained through a reversible expansion of the gas in the piston.
Refer to Section 2.3 for a discussion of reversible processes. The problem states that the piston
assembly is well-insulated, so the heat transfer contribution to the energy balance can be
neglected, in addition to potential and kinetic energy effects. The energy balance reduces to
U W
In this problem, the process is a reversible, adiabatic expansion. For this type of process,
Equation 2.90 states
W
1
P2V2 P1V1
k 1
P1 2 bar
V1 10 L
P2 1 bar
PV k const
1
P k
V2 1 V1k
P2
c 7
Noting that k P and substituting the proper values provides
cv 5
V2 16.4 L
Now all of the needed values are available for calculating the work.
W 9 L bar 900 J
U 900 J
22
The change in internal energy can also be written according to Equation 2.53:
T2
U n cv dT
T1
5
U n R T2 T1
2
T1 240.6 K
and
T2 197.3 K
The temperature is lower because more work is performed during the reversible expansion.
Review the energy balance. As more work is performed, the cooler the gas will become.
23
2.16
Since the vessel is insulated, the rate of heat transfer can be assumed to be negligible.
Furthermore, no work is done on the system and potential and kinetic energy effects can be
neglected. Therefore, the energy balance becomes
u 0
or
u 2 u1
kJ
u1 2619.2 (200 bar, 400 C)
kg
kJ
u 2 2619.2
kg
The values of u 2 and P2 constrain the system. The temperature can be found from the steam
tables using linear interpolation:
kJ
T2 327.5 C 100 bar , u 2 2619.2
kg
m3
v2 0.02012
kg
Therefore,
m3
Vvessel mv2 1.0 kg 0.02012 0.020 m
3
kg
24
2.17
Let the entire tank represent the system. Since no heat or work crosses the system boundaries,
and potential and kinetic energies effects are neglected, the energy balance is
u 0
T2 T1 0
T2 T1 300 K
The final pressure can be found using a combination of the ideal gas law and conservation of
mass.
T1 T
2
P1V1 P2V2
We also know
V2 2V1
Therefore,
P
P2 1 5 bar
2
25
2.18
(a)
First, as always, simplify the energy balance. Potential and kinetic energy effects can be
neglected. Therefore, the energy balance is
U Q W
Since, this system contains water, we can the use the steam tables. Enough thermodynamic
properties are known to constrain the initial state, but only one thermodynamic property is
known for the final state: the pressure. Therefore, the pressure-volume relationship will be used
to find the specific volume of the final state. Since the specific volume is equal to the molar
volume multiplied by the molecular weight and the molecular weight is constant, the given
expression can be written
Pv1.5 const
1
P 1.5
v2 1 v11.5
P2
Using
P1 20 bar
v1
1.0 m 3 m3
0.1
10 kg kg
P2 100 bar
gives
m3
v2 0.0342
kg
Now that the final state is constrained, the steam tables can be used to find the specific internal
energy and temperature.
T2 524.7 K
kJ
u 2 3094.6
kg
26
V2
W PE dV
V1
or
v 2
w PE dv
v1
Since the process is reversible, the external pressure must never differ from the internal pressure
by more than an infinitesimal amount. Therefore, an expression for the pressure must be
developed. From the relationship in the problem statement,
v 2 v 2
P1v11.5 1
w dv P1v11.5 v1.5 dv
v1 v1.5 v1
bar m 3 kJ
w 2.840 284
kg kg
kJ
W 10 kg 284 2840 kJ
kg
2
100
P [bar]
60
1
20 Work
27
To solve for Q , U must first be found, then the energy balance can be used.
kJ kJ
U mu 2 u1 10 kg 3094.6 2602.8 4918 kJ
kg kg
(b)
Since the final state is the same as in Part (a), U remains the same because it is a state function.
The energy balance is also the same, but the calculation of work changes. The pressure from the
weight of the large block and the piston must equal the final pressure of the system since
mechanical equilibrium is reached. The calculation of work becomes:
v 2
W mPE dv
v1
All of the values are known since they are the same as in Part (a), but the following relationship
should be noted
PE P2
W 6580 kJ
2
100
P [bar]
Work
60
1
20
28
Now Q can be solved.
(c)
This part asks us to design a process based on what we learned in Parts (a) and (b). Indeed, as is
characteristic of design problems there are many possible alternative solutions. We first refer to
the energy balance. The value of heat transfer will be zero when
U W
For the same initial and final states as in Parts (a) and (b),
W U 4918 kJ
There are many processed we can construct that give this value of work. We show two
alternatives which we could use:
Design 1:
If the answers to Part (a) and Part (b) are referred to, one can see that two steps can be used: a
reversible compression followed by an irreversible compression. Let the subscript i" represent
the intermediate state where the process switches from a reversible process to an irreversible
process. The equation for the work then becomes
vi
dv 4918000 J
1
W m P2 (v2 vi ) P1v1
1. 5
1 . 5
v1 v
Substituting in known values (be sure to use consistent units) allows calculation of vi :
m3
vi 0.0781
kg
The pressure can be calculated for this state using the expression from part (a) and substituting
the necessary values.
1.5
v
Pi P1 1
vi
Pi 29.0 bar
29
Now that both Pi and vi are known, the process can be plotted on a P-v graph, as follows:
2
100
P [bar]
60 Work
i
1
20
Design 2:
In an alternative design, we can use two irreversible processes. First we drop an intermediate
weight on the piston to compress it to an intermediate state. This step is followed by a step
similar to Part (b) where we drop the remaining mass to lead to 100 bar external pressure. In this
case, we again must find the intermediate state. Writing the equation for work:
v1
1.5
W m P1 (vi v1 ) P2 (v2 vi ) 4918000 J
vi
Solution A:
m3
vi 0.043
kg
30
which gives
Pi 70.8 bar
2
100
P [bar]
i
60
Work 1
20
Solution B
m3
vi 0.0762
kg
which gives
Pi 30.0 bar
2
100
P [bar]
60
Work 1
20
31
2.19
(a)
Force balance to find k:
Fspring=kx Fatm=PatmA
Piston
Fmass=mg
Fgas=PgasA
mg
Pgas Patm A kx
A
since V=Ax
Pgas Patm A k2V
mg
A
2 10 5 Pa 1 10 5 Pa
2040 kg 9.81 m/s 2 k 0.02 m 3
0 .1 m 2 0.1 m 2 2
N
k 5.01 10 4
m
Work can be found graphically (see P-V plot) or analytically as follows: Substituting the
expression in the force balance above:
P
W A k V
mg
atm A A 2 dV
Vf V f
mg kV
W A Patm dV d V
Vi
A
Vi A 2
5
W A 3.00 10 Pa 0.03 0.05 m 2
5.01 10 4 N/m 0.02 m 3 0 2
2
0 .1 m 2
2
2
32
W A 5 103 J
0.5kJ
-W = 10squares square = 5 kJ
3
k V
P [bar]
2
A
2
mg
2 A
Work
1
1
Patm
(b)
You need to find how far the spring extends in the intermediate (int) position. Assume
PVn=const (other assumptions are o.k, such as an isothermal process, and will change the answer
slightly). Since you know P and V for each state in Part (a), you can calculate n.
V n n
V i 10 5 Pa
0.03 m 3
Pi
or 3 n 1.35
Pf f 2 10 5 Pa 0.05 m
mg
Pgas,int Patm A kV
A2
P mg k Vint V i
V n
Pint Pi V i atm A
int A2
This last equality represents 1 equation. and 1 unknown (we know k), which gives
Vint 0.0385 m 3
33
V int Vf
W B Pgas dV PgasdV
Vi Vint
V Vint V V if
W B
int f
5 N
3 10 mN2 5 10 N5 V d V
6 N 5 6
2 10 m2 dV 5 10
m
5 Vd V dV m
Vi V i Vint Vint
Therefore,
WB 3.85 kJ
3
P [bar]
2
mg
2 A
mg
A
1
1
Work
Patm
(c)
The least amount of work is required by adding differential amounts of mass to the piston. This
is a reversible compression. For our assumption that PVn = const, we have the following
expression:
Vf Vf
const
WC , rev Pgas dV V 1.35 dV
Vi Vi
34
const. 1 105 Pa 0.05 m 3 1.35 1.75 103
Therefore,
.03
0 .35 .03
WC,rev 1.75 103 dV 10 V 2800 J
3
1.35 1.750.35
.05 V .05
3
P [bar]
2
2
1
1
Workrev
Patm
35
2.20
Before this problem is solved, a few words must be said about the notation used. The system
was initially broken up into two parts: the constant volume container and the constant pressure
piston-cylinder assembly. The subscript 1 refers to the constant volume container, 2 refers
the piston-cylinder assembly. i" denotes the initial state before the valve is opened, and f
denotes the final state.
To begin the solution, the mass of water present in each part of the system will be calculated.
The mass will be conserved during the expansion process. Since the water in the rigid tank is
saturated and is in equilibrium with the constant temperature surroundings (200 C), the water is
constrained to a specific state. From the steam tables,
kJ
v1l, i 0.001156
kg
kJ
v1v, i 0.12736
kg
kJ
u1l, i 850.64 (Sat. water at 200 C)
kg
kJ
u1v, i 2595.3
kg
P sat 1553.8 kPa
Knowledge of the quality of the water and the overall volume of the rigid container can be used
to calculate the mass present in the container.
V1 0.05m1 v1l, i 0.95m1 v1v,i
Using the values from the steam table and V1 0.5 m 3 provides
m1 4.13 kg
For the piston-cylinder assembly, both P and T are known. From the steam tables
36
m3
v2, i 0.35202
kg (600 kPa, 200 C)
kJ
u 2, i 2638.9
kg
Enough information is available to calculate the mass of water in the piston assembly.
m2 2 0.284 kg
V
v2
Now that the initial state has been characterized, the final state of the system must be determined.
It helps to consider what physically happens when the valve is opened. The initial pressure of
the rigid tank is 1553.8 kPa. When the valve is opened, the water will rush out of the rigid tank
and into the cylinder until equilibrium is reached. Since the pressure of the surroundings is
constant at 600 kPa and the surroundings represent a large temperature bath at 200 C, the final
temperature and pressure of the entire system will match the surroundings. In other words,
kJ
u f u 2, i 2638.9 (600 kPa, 200 C)
kg
U 541.0 kJ
To calculate the work, we realize the gas is expanding against a constant pressure of 600 kPa
(weight of the piston was assumed negligible). From Equation 2.7,
Vf
W PE dV PE (V f Vi )
Vi
where
PE 600000 Pa
V f (m2 m1v, i m1l,i )v2, i 1.55 m 3
Vi 0.1 m 3 0.5 m 3 0.6 m 3
37
Note: v2, i was used to calculate V f because the temperature and pressure are the same
for the final state of the entire system and the initial state of the piston-cylinder assembly.
W 570 kJ
38
2.21
A sketch of the process follows:
system system
Process
A B A B
p 1,A = 10 bar p 1,B = 20 bar p2 p2
boundary T1,A = 700 o C T1,B = 250 o C boundary T2 T2
H2 O H2 O
H2 O H2 O
50 cm 10 cm
The initial states are constrained. Using the steam tables, we get the following:
All the properties in the final state are equal. We need two properties to constrain the system:
We can find the specific volume since we know the total volume and the mass:
v2
V1, A V1, B
0.06 m 3 m3
0.30
m1, A m1, B 0.20 kg kg
We can also find the internal energy of state 2. Since the tank is well insulated, Q=0. Since it is
rigid, W=0. An energy balance gives:
U Q W 0
Thus,
We have constrained the system with u2 and v2, and can find the other properties from the steam
Tables. Very close to
Thus,
39
m3
V2, A v2, A m2, A 0.30 0.09 kg 0.267 m and
kg
40
2.22
We start by defining the system as a bubble of vapor rising through the can. We assume the
initial temperature of the soda is 5 oC. Soda is usually consumed cold; did you use a reasonable
estimate for T1? A schematic of the process gives:
State 2:
T 2= ?
P2= 1.01 bar
State 1:
T 1 = 278 K
P1 = 3 bar
where the initial state is labeled state 1, and the final state is labeled state 2. To find the final
temperature, we perform an energy balance on the system, where the mass of the system (CO2 in
the bubble) remains constant. Assuming the process is adiabatic and potential and kinetic energy
effects are negligible, the energy balance is
u w
cv T2 T1 PE dv PE v 2 v1
where cv = cP R. Since CO2 is assumed an ideal gas, the expression can be rewritten as
T T PT
cv T2 T1 PE R 2 1 R T2 2 1
P2 P1 P1
where the equation was simplified since the final pressure, P2, is equal to the external pressure,
PE. Simplifying, we get:
R RP
T2 1 T1 1 2
cv c v P1
or
RP c
T2 T1 1 2 v 237 K
cv P1 c P
41
2.23
The required amount of work is calculated as follows:
W P V
The initial volume is zero, and the final volume is calculated as follows:
W 1.01325 105 Pa 0.0436 m 3 0 m 3 4417 J
This doesnt account for all of the work because work is required to stretch the rubber that the
balloon is made of.
42
2.24
(a)
Since the water is at its critical point, the system is constrained to a specific temperature,
pressure, and molar volume. From Appendix B.1
m3
vc 0.003155
kg
Therefore,
m
V
0.01 m 3
3.17 kg
vc m3
0.003155
kg
(b)
The quality of the water is defined as the percentage of the water that is vapor. The total volume
of the vessel can be found using specific volumes as follows
V m l v l m v v v 1 x mv l xm v v
where x is the quality of the water. To solve for the quality, realize that starting with saturated
water at a pressure of 1 bar constrains the water. From the steam tables,
m3
v v 1.6940
kg
(sat. H2O at P = 1 bar)
m3
v l 0.001043
kg
x 0.00125
(c)
To determine the required heat input, perform an energy balance. Potential and kinetic energy
effects can be neglected, and no work is done. Therefore,
U Q
43
where
U mu 2 1 x mu1l xm u1v
From the steam tables
kJ
u 2 2029.58 (H2O at its critical point)
kg
l
m3
u1 417.33
kg
(sat. H2O at P=1 bar)
m3
u1v 2506.1
kg
U 5102.6 kJ 5.10 10 6 J
44
2.25
(a)
Consider the air in ChE Hall to be the system. The system is constant volume, and potential and
kinetic energy effects can be neglected. Furthermore, disregard the work. The energy balance is
du
q
dt
since the temperature of the system changes over time. Using the given expression for heat
transfer and the definition of dU , the expression becomes
cv dT
hT Tsurr
dt
We used a negative sign since heat transfer occurs from the system to the surroundings. If cv is
assumed constant, integration provides
cv ln T Tsurr ht C
h
t
T Tsurr C1e cv
where C1 is a constant. Examining this equation reveals that the temperature is an exponential
function of time. Since the temperature is decreasing, we know that the plot of temperature vs.
time shows exponential decay.
T0
Temperature
Tsurr
time
45
(b)
Let time equal zero at 6 PM, when the steam is shut off. At 6 PM, the temperature of the hall is
22 C. Therefore,
h
t
c
T Tsurr C1e v
22 C 2 C C1e (0)
C1 20 C
h
( 4 hr)
12 C 2 C 20 C e
c
v
h
0.173 hr -1
cv
At 6 AM, t 12 hr . Substitution of this value into the expression for temperature results in
T 4 .5 C
46
2.26
The gas leaving the tank does flow work as it exits the valve. This work decreases the internal
energy of the gas lowering the temperature. During this process, water from the atmosphere
will become supersaturated and condense. When the temperature drops below the freezing point
of water, the water forms a solid.
Attractive interactions between the compressed gas molecules can also contribute to this
phenomena, i.e., it takes energy to pull the molecules apart as they escape; we will learn more of
these interactions in Chapter 4.
47
2.27
Mass balance
dm
m in m out m in
dt
dm m in dt
m1 0
or
t
m2 m1 m in dt
0
Energy balance
Since the potential and kinetic energy effects can be neglected, the open system, unsteady state
energy balance is
dU
m out hout m in hin Q W s
dt sys out in
The process is adiabatic and no shaft work is done. Furthermore, there is only one inlet stream
and not outlet stream. Therefore, the energy balance simplifies to
dU
m in hin
dt sys
dU m in hin dt hin m in dt
U1 0 0
U 2 U 1 m2 u 2 m1u1 m2 m1 hin
where the results of the mass balance were used. Both m2 and m1 can be calculated by dividing
the tank volume by the specific volume
48
V
m2
v2
V
m1
v1
u 2 hin u1 hin
0
v2 v1
kJ
u1 2583.6
kg
(sat. H2O vapor at 1 MPa)
m3
v1 0.19444
kg
kJ
hin 3177.2 (6 MPa, 400 C)
kg
There are still two unknowns for this one equation, but the specific volume and internal energy
are coupled to each other. To solve this problem, guess a temperature and then find the
corresponding volume and internal energy values in the steam tables at 6 MPa. The correct
temperature is the one where the above relationship holds.
Interpolation between 500 C and 450 C reveals that the final temperature is
T2 464.4 C
49
2.28
We can pick room temperature to be 295 K
Tin T1 295 K
Mass balance
dn
n in n out n in
dt
dn nin dt
n1 0
or
t
n2 n1 nin dt
0
Energy balance
Neglecting ke and pe, he unsteady energy balance, written in molar units is written as:
dU
nin hin nout hout Q W
dt sys
The terms associated with flow out, heat and work are zero.
dU
nin hin
dt sys
Integrating both sides with respect to time from the initial state where the pressure is 10 bar to
the final state when the tank is at a pressure of 50 bar gives:
U2 t t
since the enthalpy of the inlet stream remains constant throughout the process. Integrating and
using the mass balance above:
50
Now we do some math:
n2 u 2 hin n1 u1 hin
By the definition of h
so
n2 u 2 u1 n2 RT1 n1 u1 u1 n1R T1
n2 cv T2 T1 n2 RT1 n1RT1
3
Since cv c P R R
2
3
n2 T2 T1 n2T1 n1T1
2
or
3n2T2 5n2T1 2n1T1
dividing by n1:
n n
3 2 T2 5 2 T1 2T1
n1 n1
n2 P2T1
n1 P1T2
so
P T P T
3 2 1 T2 5 2 1 T1 2T1
P1T2 P1T2
or
51
P T
5 2 1
T2 P1 434 [K]
P2
3 2
P1
u q w q
u
cv
T2 T1 5R T2 T1 28.9 kJ
MW 2MW kg
kJ
q 28.9
kg
52
2.29
Mass balance
dn
n in n out n in
dt
n2 t
dn nin dt
n1 0 0
or
t
n2 nin dt
0
Energy balance
Neglecting ke and pe, the unsteady energy balance, in molar units, is written as:
dU
nin hin nout hout Q W
dt sys
The terms associated with flow out and heat are zero.
dU
nin hin W
dt sys
Integrating both sides with respect to time from the empty initial state to the final state gives:
U2 t t t
since the enthalpy of the inlet stream remains constant throughout the process. The work is
given by:
53
n2u2 n2 hin Pext v2
Rearranging,
P T
cv T2 Tin RTin R ext 2
P2
so
T2
cv R Tin 333 K
P
cv ext R
P2
54
2.30
valve maintains
pressure in system
constant
v
T1 = 200 oC
x1 = 0.4
V = 0.01 m3
Mass balance
dm
m in m out m out
dt
dm m out dt
m1 0
or
t
m2 m1 m out dt
0
Energy balance
dU
m out hout Q
dt sys
Integrating
m 2 u 2 t t t
55
We can look up property data for state 1 and state 2 from the steam tables:
m3
v1 (1 x)v f xvg 0.6 .001 0.4 0.1274 0.051
kg
m3
v2 0.1274
kg
m1
V1
0.01 m 3
0.196 kg
v1 m3
0.051
kg
m2
V2
0.01 m 3
0.0785 kg
v2 m3
0.1274
kg
m2 m1 0.1175 kg
kJ
u1 (1 x)u f xu g 0.6 850.64 0.4 2597.5 1549
kg
kJ
u 2 2595.3
kg
kJ
hout 2793.2
kg
56
2.31
Consider the tank as the system.
Mass balance
dm
m in m out m in
dt
m2 t
dm m in dt
m1 0
or
t
m2 m1 m in dt
0
Energy balance
Since the potential and kinetic energy effects can be neglected, the open system, unsteady state
energy balance is
dU
m out hout m in hin Q W s
dt sys out in
The process is adiabatic and no shaft work is done. Furthermore, there is only one inlet stream
and not outlet stream. Therefore, the energy balance simplifies to
dU
m in hin
dt sys
U2 t t
dU m
U1 0 0
h dt hin m in dt
in in
0
U 2 m2 u 2 m2 hin
u 2 hin
57
kJ
u 2 3632.5 (9 MPa, 800 C)
kg
so
kJ
hin 3632.5
kg
We can use the value of hin and the fact that the steam in the pipe is at 9 MPa to find the
temperature.
Tin 600 C
58
2.32
(a)
First, the energy balance must be developed. Since the problem asks how much energy is stored
in the battery after 10 hours of operation, the process is not steady-state. Let the battery be the
system. Potential and kinetic energy effects can be neglected. Furthermore, heating of the
battery as it is charged can be ignored. The energy balance is
dU
Q W s
dt sys
No shaft work is performed, but electrical is supplied to the battery, which must be accounted for
in W s . The value of Q is given explicitly in the problem statement. Both of these values remain
constant over time, so integration provides
U Q W s t
From the problem statement
W s 5 kW
Q 1 kW
t 36000 s
Substituting these values allows the calculation of the amount of energy stored:
U 144,000 kJ 144 MJ
(b)
To calculate the velocity of the falling water, an energy balance must be developed with the
water passing through the electricity generator (probably a turbine) as the system, where the
water enters with a velocity V1 and leaves with a negligible velocity, which will be approximated
as 0. Assume that potential energy changes can be neglected. Furthermore, assume that the
temperature of the water does not change in the process, so the change in internal energy is zero.
Also, view the process as adiabatic. The energy balance reduces to
E K W river
where W river is the power of the flowing water. The actual power being provided by the stream
can be calculated using the efficiency information. Let represent the efficiency.
59
W
s
Wriver
W 5 kW
W river s 10 kW
0.5
The value of W river should be negative since the water is supplying work that is stored electrical
energy. Therefore, the energy balance becomes
E K 10000 W
1 2 2
2
m V2 V1 10000 W
kg
m 200
s
m
V2 0
s
Therefore,
m
V1 10
s
There are a number of reasons for the low conversion efficiency. A possible potential energy
loss inherent in the design of the energy conversion apparatus decreases the efficiency. Heat is
lost to the surroundings during conversion. Some of the energy is also lost due to friction (drag)
effects.
60
2.33
Considering the turbine to be the system, rearrangement of the steady-state, open system energy
balance provides
n in n1 n out n 2
Assuming the rate of heat transfer and potential energy effects are negligible and realizing that
there is one inlet and one outlet allows the simplification of the above equation to
W s n 2 h2 h1 e K , 2 e K ,1
h2 h1 c p dT R A BT CT 2 DT 2 ET 3 dT
T2
T1
Since the quantity ek ,2 ek ,1 is multiplied by n , it is rewritten as follows for dimensional
homogeneity
e k ,2
2 2
ek ,1 ( MW ) air V2 V1
1
2
To solve for n , the ideal gas law is used
P2V2 n 2 RT21
P V
n 2 2 2
RT2
To solve for the volumetric flow rate, the fluid velocity must be multiplied by the cross-sectional
area
D2 2V2
V2
4
61
D 2V T2 2 2
Ws
P2
RT2
4
2 2 2
3 1
2
R A BT CT DT ET dT ( MW ) air V2 V1
2
T1
Substituting values from Table A.2.1 and the problem statement results in
62
2.34
First, a sketch of the process is useful:
30 bar 20 bar
100 oC 150 oC
To find the heat in we will apply the 1st law. Assuming steady state, the open system energy
balance with one stream in and one stream out can be written:
0 n h1 h2 Q
Q
h2 h1
n
Thus this problem reduces to finding the change in the thermodynamic property, enthalpy from
the inlet to the outlet. We know 2 intensive properties at both the inlet and outlet so the values
for the other properties (like enthalpy!) are already constrained. From Table A.2.1, we have an
expression for the ideal gas heat capacity:
cp
1.424 14.394 10 3 T 4.392 10 6 T 2
R
with T in (K). Since this expression is limited to ideal gases any change in temperature must be
under ideal conditions. From the definition of heat capacity:
T2
Q
h2 h1 1.424 14.394 10 3 T 4.392 10 6 T 2 dT
n
T1
Q J
5590
n mol
63
2.35
A schematic of the process follows:
Ws
T1 = 350 oC T1 = 308 oC
To solve for W s / n we need a first law balance. With negligible eK and eP, the 1st law for a
steady state process becomes:
0 n h1 h2 Q W s
W s
h
n
We can calculate the change in enthalpy from ideal gas heat capacity data provided in the
Appendix.
T2 T2
W s
h c p dT R 1.213 28.785 10 3 T 8.824 10 6 T 2 dT
n
T1 T1
W s J
5358
n mol
64
2.36
(a)
First start with the energy balance. Nothing is mentioned about shaft work, so the term can be
eliminated from the energy balance. The potential and kinetic energy effects can also be
neglected. Since there is one inlet and one outlet, the energy balance reduces to
Q n 2 h2 n1h1
n 2 n1
Q n1 h2 h1
Using the expressions from Appendix A.1, the energy balance becomes
T2
Q n1R A BT CT 2 DT 2 ET 3 dT
T1
Using
A 3.376
B 0.557 10 3
C 0
D 0.031 105
E 0
J
R 8.314
mol K
mol
n1 20
s
T1 373.15 K
T2 773.15 K
gives
Q 245063 W 245.1 kW
65
(b)
To answer this question, think about the structure of n-hexane and carbon monoxide. N-hexane
is composed of 20 atoms, but carbon monoxide has two. One would expect the heat capacity to
be greater for n-hexane since there are more modes for molecular kinetic energy (translational,
kinetic, and vibrational). Because the heat capacity is greater and the rate of heat transfer is the
same, the final temperature will be less.
66
2.37
First start with the energy balance around the nozzle. Assume that heat transfer and potential
energy effects are negligible. The shaft work term is also zero. Therefore, the energy balance
reduces to
n 2 (h e K ) 2 n1 (h e K )1 0
n1 n 2
1
h2 h1 eK ,1 eK ,2 V12 V2 2
2
Since the steam outlet velocity is much greater than the velocity of the inlet, the above
expression is approximately equal to
1
h2 h1 V2 2
2
J
h1 2827.9 103 (10 bar, 200 C)
kg
J
h2 2675.5 103 (sat. H2O(v) at 100 kPa)
kg
Therefore,
m
V2 552
s
67
m3
v2 1.6940
kg
A 3.07 10 3 m 2
68
2.38
First start with the energy balance around the nozzle. Assume that heat transfer and potential
energy effects are negligible. The shaft work term is also zero. Therefore, the energy balance
reduces to
n 2 (h ek ) 2 n1 (h ek )1 0
The molar flow rates can be eliminated from the expression since they are equal. Realizing that
e K ,2 e K ,1 since the velocity of the exit stream is much larger than the velocity of the inlet
stream simplifies the energy balance to
h2 h1 ek ,2
T2
1
h2 h1 R A BT CT 2 DT 2 ET 3 dT ( MW ) C 3 H 8 V22
2
T1
A 1.213
B 28.785 10 3
C 8.824 10 6
D0
E 0
It is also important that the units for the molecular weight and universal gas constant are
consistent. The following values were used
J
R 8.314
mol K
kg
( MW ) C3 H 8 0.0441
mol
T2 419.2 K
69
2.39
First start an energy balance around the diffuser. Assume that heat transfer and potential energy
effects are negligible. The shaft work term is also zero. The energy balance reduces to
n 2 (h ek ) 2 n1 (h ek )1 0
n1 n 2
The molar flow rates can be eliminated from the expression. Using the definitions of enthalpy
and the kinetic energy, the equation can be rewritten as
T2
1 2 2
c P dT
2
( MW ) air V2 V1
T1
The temperature and velocity of the outlet stream are unknown, so another equation is needed to
solve this problem. From the conservation of mass,
P1V1 P1 A1V1 P2 A2V2
T1 T1 T2
where A2, the cross-sectional area of the diffuser outlet, is twice the area of the inlet. Therefore,
1 P T
V2 1 2 V1
2 P2 T1
Using Appendix A.2 and the above expression, the energy balance becomes
T2 2
2 1 1 P1 T2
R A BT CT DT 2
ET dT ( MW ) air V1 V12
3
2 2 P2 T1
T1
Substituting values from the problem statement provides an equation with one unknown:
T2 381 K
Therefore,
1 1 bar 381 K m
V2 300 m/s 111
2 1.5 bar 343.15 K s
70
2.40
To find the minimum power required for the compressor, one must look at a situation where all
of the power is used to raise the internal energy of the air. None of the power is lost to the
surroundings and the potential and kinetic energy effects must be neglected. Therefore, the
energy balance becomes
0 n1h1 n 2 h2 W s
n1 n 2
W s n1 h2 h1
T2
W s n1R A BT CT 2 DT 2 ET 3 dT
T1
Table A.2.1 and the problem statement provide the following values
A 3.355
B 0.575 10 3
C 0
D 1600
E 0
mol
n1 50
s
T1 300 K
To find the work, we still need T2. We need to pick a reasonable process to estimate T2. Since
the heat flow is zero for this open system problem, we choose an adiabatic, reversible piston
situation. For this situation,
PV k const.
Since we are assuming the air behaves ideally, we can rewrite the equation as
71
k k
n RT n RT
P1 1 1 P2 2 2
P1 P2
P11 k T1k P21 k T2k
1 7 / 5 5/7
7 / 5 1 bar
T2 300 K 579 K
10 bar 1 7 / 5
Substitute this value into the expression for the work and evaluate:
W s 417.4 kW
72
2.41
(a)
Perform a mass balance:
n1 n 2 n out
P1V1 P2V2
n out
RT1 RT2
(b)
No work is done on the system, and we can neglect potential and kinetic energy effects. We will
assume the process is also adiabatic. The energy balance reduces to
We can calculate the enthalpy difference from the given ideal heat capacity:
Again, we must calculate the molar flow rates from the ideal gas law. Upon substitution and
evaluation, we obtain
Tout 329 K
73
2.42
(a)
Since the temperature, pressure, and volumetric flow rate are given, the molar flow rate is
constrained by the ideal gas law.
Note: 1.67 10 8 m 3 /s 1 cm 3 / min
n
PV 1.0135 10 5 Pa 1.67 10 8 m 3 /s 7.45 10 7 mol/s
RT 8.314 J/mol K 273.15 K
To recap, we have shown
(b)
Assumptions: N2 is an ideal gas
All power supplied by the power supply is transferred to the N2
Uniform temperature radially throughout sensor tube
Kinetic and potential energy effects negligible in energy balance
Let x represent the fraction of N2 diverted to the sensor tube, and n s represent the molar flow
rate through the sensor tube. Therefore, the total molar flow rate, ntotal , is
n
ntotal s
x
We can use temperature and heat load information from the sensor tube to find the molar flow
rate through the sensor tube. First, perform an energy balance for the sensor tube:
H n s hout hin Q
Q
n s
T2
cP dT
T1
74
Q
ntotal
T2
x c P dT
T1
To find the flow rate in standard cubic centimeters per minute, apply the conversion factor found
in Part (a)
Q
1 SCCM
vtotal (SCCM )
T2 7.45 10 7 mol/s
x c P dT
T1
(c)
To find the correction factor for SiH4, re-derive the expression for flow rate for SiH4 and then
divide it by the expression for N2 for the same power input, temperatures, and fraction of gas
diverted to the sensor tube.
Q
1 SCCM
T2 7.45 10 7 mol/s T2
x c P, SiH 4 dT c P, N 2
dT
vtotal , SiH 4 T1 T1
Factor
vtotal , N 2 Q
1 SCCM
T2
T2 7.45 10 7 mol/s
c P, SiH 4
dT
x c P, N 2 dT T1
T1
If we assume that heat capacities are constant, the conversion factor simplifies:
c P, N 2
Factor
c P, SiH 4
Factor 0.67
75
2.43
(a)
It takes more energy to raise the temperature of a gas in a constant pressure cylinder. In both
cases the internal energy of the gas must be increased. In the constant pressure cylinder work,
Pv work must also be supplied to expand the volume against the surroundings pressure. This is
not required with a constant volume.
(b)
As you perspire, sweat evaporates from your body. This process requires latent heat which cools
you. When the water content of the environment is greater, there is less evaporation; therefore,
this effect is diminished and you do not feel as comfortable.
76
2.44
From the steam tables at 10 kPa:
T(K) h
323.15 2592.6
h vs. T
373.15 2687.5
423.15 2783 6000
473.15 2879.5 y = 0.0003x
2
+ 1.6241x + 2035.7
5000 2
523.15 2977.3 R =1
4000
573.15 3076.5
673.15 3279.5 3000
dh kJ J
cP 1.6241 0.0006624T kg K 3.516 0.001434T R mol K
dT P
A B
Steam Tables 3.516 0.001434
Appendix A 3.470 0.001450
% difference 1.3 1.4
77
2.45
For throttling devices, potential and kinetic energy effects can be neglected. Furthermore, the
process is adiabatic and no shaft work is performed. Therefore, the energy balance for one inlet
and one outlet is simplified to
n1 h1 n 2 h2
which is equivalent to
m 1 h1 m 2 h2
h1 m 2 h2
kJ
h1 3398.3 (8 MPa, 500 C)
kg
kJ
h2 3398.3
kg
Now that we know u 2 and P2 , T2 is constrained. Linear interpolation of steam table data gives
T2 457 C
78
2.46
(a)
An expression for work in a reversible, isothermal process was developed in Section 2.7.
Equation 2.77 is
P
W nRT ln 2
P1
Therefore,
P
w RT ln 2
P1
J
R 8.314
mol K
T 300 K
P2 100 kPa
P1 500 kPa
gives
J
w 4014
kg
(b)
Equation 2.90 states
w
R
T2 T1
k 1
5
cP R
2
3
cv R
2
and
79
5
k
3
T2 can be calculated by applying the polytropic relation derived for adiabatic expansions. From
Equation 2.89
PV k const
Pv k const
RT
v
P
Since R is a constant, substitution of the expression for P into the polytropic relation results in
P 1 k T k const
P11 k T 1 k P2 1 k T2 k
T2 157.6 K
J
w 1775.9
kg
80
2.47
(a)
The change in internal energy and enthalpy can be calculated using
We would like to calculate these values using the steam tables; however, the appendices dont
contain steam table data for liquid water at 0.0 C and 1 atm. However, information is provided
for water at 0.01 C and 0.6113 kPa. Since the enthalpy and internal energy of liquid water is
essentially independent of pressure in this pressure and temperature range, we use the steam
table in the following way
kJ
h l 1 atm, 100 C 419.02
kg
kJ
h l 1 atm, 0 C h l 0.6113 kPa , 0.01 C 0
kg
kJ
u l 1 atm, 100 C 418.91
kg
kJ
u l 1 atm, 0 C u l 0.6113 kPa , 0.01 C 0
kg
Therefore,
kJ
h 419.02
kg
kJ
u 418.91
kg
(b)
The change in internal energy and enthalpy can be calculated using
81
kJ
h v 2676.0
kg
kJ
h l 419.02
kg
kJ
u v 2506.5
kg
kJ
u l 418.91
kg
Therefore,
kJ
h 2256.99
kg
kJ
u 2087.59
kg
The change in internal energy for the process in Part (b) is 5.11 times greater than the change in
internal energy calculated in Part (a). The change in enthalpy in Part (b) is 5.39 times greater
than the change in enthalpy calculated in Part (a).
82
2.48
To calculate the heat capacity of Ar, O2, and NH3 the following expression, with tabulated values
in Table A.2.1, will be used,
cP
A BT CT 2 DT 2 ET 3
R
T 300 K
J
R 8.314
mol K
The values are not listed for Ar since argon can be treated as a monatomic ideal gas with a heat
capacity independent of temperature. The expression for the heat capacity is
5
c P, Ar R
2
Now that expressions exist for each heat capacity, evaluate the expressions for T 300 K .
J
c P, Ar 20.785
mol K
J
c P, O2 29.420
mol K
J
c P, NH 3 35.560
mol K
By examining the heat capacity for each molecule, it should be clear that the magnitude of the
heat capacity is directly related to the structure of the molecule.
83
Ar
Since argon is monatomic, translation is the only mode through which the atoms can
exhibit kinetic energy.
O2
Translation, rotation, and vibration modes are present. Since oxygen molecules are
linear, the rotational mode of kinetic energy contributes RT per mol to the heat capacity.
NH3
Translation, rotation, and vibration modes are present. Ammonia molecules are non-
linear, so the rotation mode contributes 3RT/2 per mole to the heat capacity.
The vibration contributions can also be analyzed for oxygen and ammonia, which reveals that the
vibration contribution is greatest for ammonia. This is due to ammonias non-linearity.
84
2.49
For a constant pressure process where potential and kinetic energy effects are neglected, the
energy balance is given by Equation 2.57:
Q H
H mh2 h1
kJ
h2 hsat. vapor at 10 kPa 2584.6
kg
kJ
h1 hsat. liquid at 10 kPa 191.81
kg
Therefore,
kJ kJ
Q 2 kg 2584.6 191.81
kg kg
Q 4785.6 kJ
Vf
W PE dV
Vi
The pressure is constant, and the above equation can be rewritten as follows
W PE mv2 v1
m3
v2 vsat. vapor at 10 kPa 14.674
kg
m3
v1 vsat. liquid at 10 kPa 0.00101
kg
Therefore,
85
m3
W 10000 Pa 2 kg 14.674 0.00101 293.5 kJ
kg kg
86
2.50
First, perform an energy balance on the system. Potential and kinetic energy effects can be
neglected. Since nothing is mentioned about work in the problem statement, W can be set to
zero. Therefore, the energy balance is
Q U
m2 m1
where
m1 m1l m1v
Q m1 (u 2 u1 ) m1u 2 m1l u1l m1v u1v
Since two phases coexist initially (water is saturated) and P1 is known, state 1 is constrained.
From the saturated steam tables
kJ
u1l 191.79
kg
(sat. H2O at 10 kPa)
kJ
u1v 2437.9
kg
As heat is added to the system, the pressure does not remain constant, but saturation still exists.
One thermodynamic property is required to constrain the system. Enough information is known
about the initial state to find the volume of the container, which remains constant during heating,
and this can be used to calculate the specific volume of state 2.
l l v v
V2 V1 m1 v1 m1 v1
v2
m2 m1 m1l m1v
87
m3
v1l 0.001010
kg
(sat. H2O at 10 kPa)
m3
v1l 14.674
kg
Therefore,
m3
v2 1.335
kg
The water vapor is now constrained. Interpolation of steam table data reveals
kJ
u 2 2514.6
kg
Now that all of the required variables are known, evaluation of the expression for Q is possible.
Q 11652 kJ
88
2.51
Let the mixture of ice and water immediately after the ice has been added represent the system.
Since the glass is adiabatic, no work is performed, and the potential and kinetic energies are
neglected, the energy balance reduces
H 0
We can split the system into two subsystems: the ice (subscript i) and the water (subscript w).
Therefore,
H mi hi mw hw 0
and
mi hi mw hw
mw
Vw
0.0004 m 3
0.399 kg
vw m3
0.001003
kg
nw
mw
0.399 kg 22.15 mol
MW H 2O 0.0180148 kg
mol
5.55 mol
mi
ni
MW H 2O
Now, lets assume that all of the ice melts in the process. (If the final answer is greater than 0
C, the assumption is correct.) The following expression mathematically represents the change
in enthalpy.
ni c P, i 0 10 C h fus c P, w T f 0 C nwc P, w T f 25 C
Note: Assumed the heat capacities are independent of temperature to obtain this
expression.
c P ,i 4.196 R
c P ,w 9.069 R
and
kJ
h fus 6.0
mol
89
Substitution of values into the above energy balance allows calculation of Tf.
T f 3.12 C
(b)
To obtain the percentage of cooling achieved by latent heat, perform the following calculation
Fractionlatent
ni h fus
nw c P, w T f Tw, i
5.55 mol 6000
J
mol
Fractionlatent 0.911
J
22.15 mol9.069 8.314 3.12 C 25 C
mol K
Percentlatent 91.1%
90
2.52
A mass balance shows
n2 n1l n1v
To develop the energy balance, neglect kinetic and potential energy. Also, no shaft work is
performed, so the energy balance becomes
U Q
Q n1l n1v u 2 n1l u1l n1v u1v
If the reference state is set to be liquid propane at 0 C and 4.68 bar, the internal energies become
u1l 0
u1v u vap 0 C
T2
u 2 u vap 0 C c P R dT
273 K
Once the change in internal energy for vaporization and temperature of state 2 is determined, Q
can be solved. As the liquid evaporates, the pressure increases. At state 2, where saturated
propane vapor is present, the ideal gas law states
RT2
P2
v2
To find v2 , assume that v1v v1l . The volume of the rigid container is
RT1 m3
V n1v v1v n1v 0.00485
P1 mol
Therefore,
V m3
v2
n1l n1v
0.00243
mol
91
Also, since the propane is saturated, P2 and T2 are not independent of each other. They are
related through the Antoine Equation,
ln P sat A
T
B
sat
C
where
Substitution provides,
RT B
ln 2 A
v2 T2 C
bar m 3
Using values from Table A.1.1 and R 8.314 10 5
mol K
T2 301.7 K
hvap h v h l (u Pv) v u Pv l
Since v v v l , the above the change in internal energy of vaporization can be written as
Therefore,
kJ
u vap 0 C 14.39
mol
Evaluation of the following equation after the proper values have been substituted from Table
A.2.1
301.7 K
Q n1l n1v u vap 0 C c P R dT n1l 0 n1v uvap 0 C
273 K
92
gives
Q 18.1 kJ
93
2.53
The equation used for calculating the heat of reaction is given in Equation 2.72. It states
hrxn vi h f i
This equation will be used for parts (a)-(e). Since the heat of reaction at 298 K is desired, values
from Appendix A.3 can be used.
(a)
First the stoichiometric coefficient must be determined for each species in the reaction.
vCH 4 ( g ) 1
vO2 ( g ) 1
vCO2 ( g ) 1
v H 2 O( g ) 1
hf ,298 O ( g ) 0 mol
2
kJ
hrxn
,298 vCH 4 ( g ) h f ,298 CH ( g ) vO2 ( g ) h f ,298 O ( g ) vCO2 ( g ) h f ,298 CO ( g ) v H 2 O ( g ) h f ,298 H O ( g )
4 2 2 2
kJ kJ kJ kJ
hrxn , 298 74.81 0 393.51 241.82
mol mol mol mol
kJ
hrxn , 298 560.52
mol
Now that a sample calculation has been performed, only the answers will be given for the remaining
parts since the calculation process is the same.
94
(b)
kJ
hrxn , 298 604.53
mol
(c)
kJ
hrxn , 298 206.12
mol
(d)
kJ
hrxn , 298 41.15
mol
(e)
kJ
hrxn , 298 905.38
mol
95
2.54
The acetylene reacts according to the following equation
(a)
First, choose a basis for the calculations.
nC 2 H 2 1 mol
Calculate the heat of reaction at 298 K using Equation 2.72 and Appendix A.3
hrxn i
vi h f
hrxn h f O 2hf CO hf H O
2.5 h f
C H2 2 2 2 2
J
hrxn 1.255 10 6
mol
H rxn,298 nC 2 H 2 hrxn
1.255 10 6 J
nO 1 2.5(nC H
2 2 2
)1 2.5 mol
96
T2
H rxn,298 ni 2 c P i dT 0
298 i
T2 6169 K
(b)
The calculations follow the procedure used in Part (a), but now nitrogen is present. The basis is
nC 2 H 2 1 mol
The heat of reaction is the same as in Part (a), but the gas composition is different. Since
stoichiometric amount of air is used,
nO 1 2.5 mol
2
yN
n N 2 nO2 2 9.40 mol
1 1 y
O2 air
Therefore,
T2 2792 K
97
(c)
Now excess air is present, so not all of the oxygen reacts. The heat of reaction remains the same
because only 1 mole of acetylene reacts. Since the amount of air is twice the stoichiometric
amount
nO 1 5 mol
2
yN
n N 2 nO2 2 18.80 mol
1 1 y
O2 air
The table of heat capacity data in Part (b) will be used for this calculation. Using the expression
shown in Part (a)
T2 1787 K
98
2.55
(a)
The combustion reaction for propane is
For all subsequent calculations, the basis is one mole of propane. The heat of reaction is
calculated as follows
hrxn C H 5hf O 3hf CO 4hf H O
h f
3 8 2 2 2
J
hrxn 2.044 10 6
mol
H rxn,298 nC 3 H 8 hrxn
2.044 10 6 J
nO 1 5(nC H
2 3 8
)1 5 mol
yN
n N 2 nO2 2 18.81 mol
1 1 y
O2 air
The stream compositions are listed below
Now all of the necessary variables for the following equation are known, except T2 .
T2
H rxn,298 ni 2 c P i dT 0
298 i
T2 2374 K
99
(b)
The combustion reaction for butane is
For all subsequent calculations, the basis is one mole of butane. The heat of reaction is
calculated as shown in Part (a)
H rxn,298 2.657 10 6 J
The moles of nitrogen and oxygen in the feed stream are calculated according to the method in
Part (a). The compositions are
The c P data listed in Part (a) can also be used for this reaction since there is no remaining
butane.
T2 2376 K
(c)
The combustion reaction for pentane is
The basis is one mole of pentane. The heat of reaction is calculated as shown in Part (a).
H rxn,298 3.272 10 6 J
The moles of nitrogen and oxygen in the feed stream are calculated according to the method in
Part (a). The compositions are listed below
Substitution of the values into the expression used to find T2 and subsequent evaluation results
in
T2 2382 K
The adiabatic flame temperatures are nearly identical in all three cases.
100
2.56
The equation for the combustion of methane is
J
hrxn 8.02 105
mol
nCH 1 1 mol
4
Also, let represent the fractional conversion of methane. Therefore, the composition of the
product gas leaving the reactor is
nCH 2 11 mol
nO 2 21 mol
4
nN 2 7.52 mol
2
nCO 2 mol
2
nH O 2 2 mol
2
H rxn,298 8.02 105 J
After substitution of the outlet composition values, heat capacity data, and the heat of reaction
into the following equation
101
1273
H rxn, 298 n c dT 0
298 i
i 2 P i
integration provides an equation with one unknown: . Solving the equation gives
=0.42
Since the fractional conversion is 0.42, 58% of the methane passed through the reactor unburned.
102
2.57
For the entire cycle,
U11 U12 U 23 U 31 0
U12 50 kJ
U 23 Q23 W23
Q23 0 kJ
U 31 Q31 W31
W31 250 kJ
Process U kJ W kJ Q kJ
State 1 to 2 -50 -400 350
State 2 to 3 800 800 0
State 3 to 1 -750 -250 -500
To determine if this is a power cycle or refrigeration cycle, look at the overall heat and work,
W11 and Q11 .
Since work is done on the system to obtain a negative value of heat, which means that heat is
leaving the system, this is a refrigeration cycle.
103
2.58
Refer to the graph of the Carnot cycle in Figure E2.20. From this graph and the description of
Carnot cycles in Section 2.9, it should be clear that state 3 has the lowest pressure of all 4 states,
and state 1 has the highest pressure. States 1 and 2 are at the higher temperature. States 3 and 4
have the lower temperature. Since both the temperature and pressure are known for states 1 and
3, the molar volume can be calculated using
RT
v
P
State T K P bar
v m 3 /mol
1 1073 60 0.00149
2 1073
3 298 0.2 0.124
4 298
For each step of the process, potential and kinetic energy effects can be neglected. The step from
state 1 to state 2 is a reversible, isothermal expansion. Since it is isothermal, the change in
internal energy is 0, and the energy balance becomes
Q12 W12
P
W12 nRT1 ln 2
P1
where W12 is the work done from state 1 to state 2. The value of P2 is not known, but
recognizing that the process from state 2 to 3 is an adiabatic expansion provides an additional
equation. The polytropic relationship can be employed to find P2 . A slight modification of
Equation 2.89 provides
Pv k const
RT
P
v
Combining this result with the polytropic expression and noting that R is constant, allows the
expression to be written as
104
Tv k 1 const
Therefore,
1
T k 1
v2 3 v3k 1
T2
m3
v2 0.00504
mol
17.7 bar
RT2
P2
v2
W12 10.89 kJ
Calculation of W34 follows a completely analogous routine as calculation for W12 . The
following equations were used to find the necessary properties
1
T k 1 m3
v4 1 v1k 1 0.0367
T4 mol
0.675 bar
RT4
P4
v4
P
W34 nRT3 ln 4
P3
which gives
W34 3.01 kJ
105
For an adiabatic, reversible process, Equation 2.90 states
W
nR
T2 T1
k 1
This equation will be used to calculate the work for the remaining processes.
W23
nR
T3 T2 16.11 kJ
k 1
W41
nR
T1 T4 16.11 kJ
k 1
To find the work produced for the overall process, the following equation is used
Wnet 7.88 kJ
Wnet 7.88
0.72
QH 10.89
TC
1
TH
0.72
106
2.59
Since this is a refrigeration cycle, the direction of the cycle described in Figure 2.17 reverses.
Such a process is illustrated below:
Constant T H
QH
Isothermal
compression
Adiabatic
compression Adiabatic
State 2 State 1 expansion
Well - insulated
T2=T1 T1
P2 P1
State 3 State 4
T3 T4=T3 State 2 State 3
State 1 State 4
T1 T4 P3 P4 T2 T3
P1 P4 P2 P3
Isothermal
expansion
QC
Constant T C
States 1 and 2 are at the higher temperature. States 3 and 4 have the lower temperature. Since
the both the temperature and pressure are known for states 2 and 4, the molar volume can be
calculated using
RT
v
P
State T K P bar
v m 3 /mol
1 1073
2 1073 60 0.00149
3 298
4 298 0.2 0.124
For each step of the process, potential and kinetic energy effects can be neglected. The process
from state 1 to state 2 is a reversible, isothermal expansion. Since it is isothermal, the change in
internal energy is 0, and the energy balance becomes
Q12 W12
107
From Equation 2.77,
P
W12 nRTH ln 2
P1
where W12 is the work done from state 1 to state 2. The value of P1 is not known, but
recognizing that the process from state 4 to 1 is an adiabatic compression provides an additional
relation. The polytropic relationship can be employed to find P1 . A slight modification of
Equation 2.89 provides
Pv k const
RT
P
v
Combining this result with the polytropic expression and noting that R is constant allows the
expression to be written as
Tv k 1 const
Therefore,
1
T k 1
v1 4 v 4k 1
T1
1
T k 1 m3
v1 4 v 4k 1 0.00504
T1 mol
P1 1 17.7 bar
RT
v1
W12 10.9 kJ
Calculation of W34 follows a completely analogous routine as the calculation for W12 . The
following equations were used to find the necessary properties:
108
m3
v3 0.0367
mol
0.675 bar
RT3
P3
v3
P
W34 nRTC ln 4
P3
which gives
W41 3.0 kJ
Equation 2.90 can be used to determine the work for adiabatic, reversible processes. This
equation will be used to calculate the work for the remaining processes.
W23
nR
TC TH 16.11 kJ
k 1
W41
nR
TH TC 16.11 kJ
k 1
To find the work produced for the overall process, the following equation is used
Wnet 7.88 kJ
Therefore, 7.88 kJ of work is obtained from the cycle. The coefficient of performance is defined
in Equation 2.99 as follows
QC
COP
Wnet
where QC is the equal to Q34 . From the energy balance developed for the process from state 3
to state 4
109
Q34 W34 3.01 kJ
Therefore,
3.01 kJ
COP 0.382
7.88 kJ
110
2.60
(a)
The Pv path is plotted on log scale so that the wide range of values fits (see Problem 1.13)
logP
4 1
100
3
0.075
2
log v
v
(b)
The work required to compress the liquid is the area under the Pv curve from state 3 to state 4.
Its sign is positive. The power obtained from the turbine is the area under the curve from state 1
to 2. Its sign is negative. The area under the latter curve is much larger (remember the log
scale); thus the net power is negative.
(c)
First, perform a mass balance for the entire system:
m 1 m 2 m 3 m 4 m
Since no work is done by or on the boiler, the energy balance for the boiler is
m h1 m h4 Q H
m h3 m h2 Q C
To find the necessary enthalpies for the above energy balances, we can use the steam tables:
kJ
h1 3424.5 (520 C, 100 bar)
kg
111
kJ kJ kJ
h2 0.10168.77 0.90 2574.8 2334.2
kg kg kg
(sat. liq at 7.5 kPa) (sat. vap. at 7.5 kPa)
kJ
h3 168.77 (sat. liquid at 0.075 bar)
kg
kJ
h4 342.81 (subcooled liquid at 80 C, 100 bar)
kg
kJ kJ
Q H 100 kg/s 3424.5 342.81 308169 kW
kg kg
kJ kJ
Q C 100 kg/s 168.77 2334.2 216543 kW
kg kg
(d)
Use Equation 2.96:
W net Q net 0
Therefore,
W net 91626 kW
(e)
Using the results from Parts (c) and (d):
91626 kW
0.297
308169 kW
112
Chapter 3 Solutions
Engineering and Chemical Thermodynamics
Wyatt Tenhaeff
Milo Koretsky
[email protected]
3.1
Since entropy is a state function
Step 1 is a constant volume process. Therefore, no work is done. After neglecting potential and
kinetic energy effects, the energy balance for a reversible process becomes
u q
du q
cv dT q
q rev
final T2
cv dT
s sys ,step1
initial
T
T
T1
Step 2 is an isothermal process. For and ideal gas u is zero and the energy balance is (PE and
KE neglected)
q w
dq w
RT
wrev Pdv dv
v
Therefore,
qrev v2
final v2
R
s sys ,step 2
initial
T
v dv R ln v
v1
1
T2
cv dT v
s sys T
R ln 2
v1
T1
2
3.2
T2
cP P
s sys T dT R ln 2
P1
T1
T2
A BT CT 2 P
s sys T
dT R ln 2
P1
T1
T
P
s sys A ln 2 BT2 T1 T22 T12 R ln 2
C
T1 2 P1
3
3.3
Equation 3.3 states
where
S sys m( s2 s1 )
Qsurr
S surr (the temperature of the surroundings is constant)
Tsurr
kJ
s1 7.1228 (300 C, 10 bar)
kg K
v2 v1
m3
v1 0.25794 (300 C, 10 bar)
kg
m3
v2 0.25794
kg
To find the number of phases present in state 2, compare the specific volume of state 2 to the
specific volume for saturated water and saturated water vapor at 1 bar. From the steam tables,
m3
v2l , sat 0.001043 (sat. H2O(l) at 1 bar)
kg
m3
v2v, sat 1.6940 (sat. H2O(v) at 1 bar)
kg
Since v2l , sat v2 v2v, sat , two phases are present. The quality of the water can be calculated as
follows
Therefore,
4
x 0.152
kJ
s2l , sat 1.3025 (sat. H2O(l) at 1 bar)
kg K
kJ
s2v, sat 7.3593 (sat. H2O(v) at 1 bar)
kg K
kJ kJ
S sys 10 kg 2.2231 49.0
kJ
7.1228
kg K kg K K
To find the change in entropy of the surroundings, an energy balance will be useful. Since no
work is done, the energy balance is
U Q
We also know
Qsurr Q
U m 1 x u 2l , sat xu 2v, sat u1
kJ
u1 2793.2 (300 C, 10 bar)
kg
kJ
u 2l , sat 417.33 (sat. H2O(l) at 1 bar)
kg
5
kJ
u 2v, sat 2506.1 (sat. H2O(v) at 1 bar)
kg
U 20583.77 kJ
Therefore,
Q 20583.77 kJ
and
Qsurr 20583.77 kJ kJ
S surr 70.25
Tsurr 293 K K
kJ kJ kJ
S univ S sys S surr 49.0 70.25 21.25
K K K
6
3.4
Entropy Balance:
Considering the copper block to be the system, no work is done on the system; thus, the energy
balance is
T2
q rev
s sys T
T1
T2
cv
s sys T dT
T1
c P 2.723R
cv c P 2.723R (for liquids and solids)
Therefore,
J 280 K
T2
2.723R T
s sys T
dT 2.723R ln 2 2.723 8.314
T1
ln
mol K 373.15 K
T1
J
s sys 6.50
mol K
10 kg 6.50 J 1024 J
S sys
kg mol K
K
0.063465 mol
Since the temperature of the lake remains constant, the change in entropy of the surroundings can
be calculated as follows
7
J
2.723 8.314 280 K 373.15 K
10 kg mol K
S surr
kg 280 K
0.063546 mol
J
S surr 1184
K
J J J
S univ S sys S surr 1024 1184 160
K K K
8
3.5
(a)
Since the process is adiabatic and reversible,
S sys 0
(b)
The change in entropy is calculated by
kJ
1 mol 8.2
qrev nh mol
liquid
J
vap
S sys n
vapor
T
Tb
111 K
73.9
K
The sign is negative because there is less randomness in the liquid phase.
(c)
For this situation
q c P dT
Therefore,
q rev J 273.15 K
Tf 273.15 K
S sys m
cP
dT 18.0148 g
4.2 g K ln 373.15 K
Ti
T 373.15 K
T
J
S sys 23.6
K
The sign is negative because as the water cools, less translational energy states are occupied by
the molecules. Therefore, the randomness decreases.
(d)
First, calculate the temperature at which the blocks (block A and block B) equilibrate. The
energy balance for the process is
n Ac P T2 373.15 K n B c P T2 473.15 K 0
Since the heat capacities number of moles of A and B are equal, we find that
T2 423.15 K
9
S sys S sys, A S sys, B
T T
S sys n Ac P ln 2 n B c P ln 2
T1, A T1, B
J
S sys 0.337
K
The sign is positive because two objects at different temperatures will spontaneously equilibrate
to the same temperature when placed together.
10
3.6
The solution below compares problems 2.14 and 2.15, the calculation of 2.13 was erroneously
included in the problem statement of the first printing and is shown at the end of this problem.
Problem 2.14
Since the system is well-insulated no heat is transferred with the surroundings. Therefore, the
entropy change of the surroundings is zero and
S univ S sys
T P
S sys n c P ln 2 R ln 2
T1 P1
P1 2 bar
P2 1 bar
PV
T1 1 1
2 bar 10 L 240.6 K
nR L bar
0.08314 1 mol
mol K
Solving for T2 is slightly more involved. The energy balance for this system where potential and
kinetic energy effects are neglected is
U W
n1 n2
Let n n1 n2
T2 V2
n cv dT PE dV
T1 V1
11
Since cv and PE are constant
ncv T2 T1 PE V2 V1
nRT2
V2
P2
PV
T1 1 1
nR
Substituting these expressions into the energy balance, realizing that PE P2 , and simplifying
the equation gives
5
P2 P1 V1
T2
2
7
nR
2
P1 2 bar
P2 1 bar
V1 10 L
n 1.0 mol
L bar
R 0.08314
mol K
results in
T2 206 K
Since both states are constrained, the entropy can be calculated from Equation 3.63:
J 206 K J 1 bar
S univ S sys 1.0 mol 3.5 8.314 ln 8.314 ln
mol K 240.6 K mol K 2 bar
J
Suniv S sys 1.24
K
12
Problem 2.15
Since the initial conditions in the piston-cylinder assembly are equal to the initial conditions of
Problem 2.14, T1 and P1 are known. Moreover, P2 is known, so we only need to find T2 in
order to calculate the entropy change. For adiabatic, reversible processes, the following
relationship (Equation 2.89) holds:
PV k const
1
P k
V2 1 V1k
P2
c 7
Noting that k P and substituting the proper values results in
cv 5
V2 16.4 L
TV k 1 const.
Therefore,
T1V1k 1
T2
V2k 1
T2 197.4 K
J 197.4 K J 1 bar
S univ S sys 1.0 mol 3.5 8.314 ln 8.314 ln
mol K 240.6 K mol K 2 bar
J
S univ S sys 0.004
K
The value of the entropy shown represents round-off error. Since the process is reversible and
adiabatic, we know from Table 3.1 and the related discussion in Section 3.3 that the entropy
changes of the system, surroundings, and universe will be zero.
13
Problem 2.13
Qsurr
S univ S sys S surr m( s2 s1 )
Tsurr
where the temperature of the surroundings is constant. First, we will determine s 2 . The first law
can be applied to constrain state 2. With potential and kinetic energy effects neglected, the
energy balance becomes
U Q W
The value of the work will be used to obtain the final temperature. The definition of work
(Equation 2.7) is
V2
W PE dV
V1
Since the piston expands at constant pressure, the above relationship becomes
W PE V2 V1
kJ
s1 6.2119 (10 MPa, 400 C)
kg K
m3
v1 0.02641 (10 MPa, 400 C)
kg
m3
V1 m1v1 (3 kg) 0.02641 0.07923 m
3
kg
W
V2 V1
PE
0.07923 m3
748740 J
6
2.0 10 Pa
0.4536 m3
V
v2 2
0.4536 m 3
0.1512
m3
m2 3 kg kg
14
Since v2 and P2 are known, state 2 is constrained. From the steam tables:
kJ 3
s2 7.1270 20 bar, 0.1512 m
kg K kg
Now U will be evaluated, which is necessary for calculating Qsurr . From the steam tables:
kJ 3
u 2 2945.2 20 bar, 0.1512 m
kg kg
kJ
u1 2832.4 100 bar, 400 C
kg
kJ kJ
U m1 u 2 u1 3 kg 2945.2 2832.4 338.4 kJ
kg kg
Substituting the values of U and W into the energy equation allows calculation of Q
Q U W
Q 338400 [J] 748740 J 1.09 10 6 J Qsurr
so
Qsurr
S univ m( s2 s1 )
Tsurr
kJ kJ 1.09 103 kJ kJ
3 kg 7.1270 6 . 2119 1.126
kg K
kg K 673.15 K K
15
3.7
(a)
T P
s sys c P ln 2 R ln 2
T1 P1
J 7 500 K 0.5 bar J
s sys 8.314 ln ln 20.63
mol K 2 300 K
1 bar mol K
(b)
From Equation 3.65:
T v
s sys cv ln 2 R ln 2
T1 v1
3
0.025 m
J 5 500 K mol 4.85 J
s sys 8.314
ln
ln
mol K 2 300 K m 3 mol K
0.05
mol
(c)
First, we can find the molar volume of state 1 using the ideal gas law.
RT1 m3
v1 0.025
P1 mol
T v
s sys cv ln 2 R ln 2
T1 v1
3
0.025 m
J 5 500 K mol
10.62
J
s sys 8.314
ln
ln
mol K 2 300 K m 3 mol K
0.025
mol
16
3.8
(i). We wish to use the steam tables to calculate the entropy change of liquid water as it goes
from its freezing point to its boiling point. The steam tables in Appendices B.1 B.5 do not
have data for subcooled water at 1 atm. However, there is data for saturated water at 0.01 C and
a pressure of 0.6113 kPa. If we believe that the entropy of water is weakly affected by pressure,
then we can say that the entropy of water at 0.01 C and 0.6113 kPa is approximately equal to the
entropy at 0 C and 1 atm. The molar volumes of most liquids do not change much with pressure
at constant temperature. Thus, the molecular configurations over space available to the water
molecules do not change, and the entropy essentially remains constant. We do not need to
consider the molecular configurations over energy since the temperature difference is so slight.
So, from the steam tables:
kJ
s0 C , 1 atm s0.01 C , 0.6113 kPa 0
kg K
kJ
s100 C , 1 atm 1.3068
kg K
Therefore,
kJ kJ
s s100 C , 1 atm s0 C , 1 atm 1.3068 0
kg K kg K
kJ
s 1.3068
kg K
kJ
sl 100 C , 1 atm 1.3068
kg K
kJ
s v 100 C , 1 atm 7.3548
kg K
Therefore,
kJ kJ
s s v 100 C , 1 atm s l 100 C , 1 atm 7.3548 1.3068
kg K kg K
kJ
s 6.048
kg K
The change in entropy for process (ii) is 4.63 times the change in entropy for process (i). There
are many ways to reconcile this difference, but think about it from a molecular point of view. In
process (i), the available molecular configurations over energy are increased as the temperature
increases. As the temperature increases, the molar volume also increases slightly, so the
17
available molecular configurations over space also increase. Now consider process (ii), where
the molecules are being vaporized and entering the vapor phase. The molecular configuration
over space contribution to entropy is drastically increased in this process. In the liquid state, the
molecules are linked to each other through intermolecular interactions and their motion is
limited. In the vapor state, the molecules can move freely. Refer to Section 3.10 for a discussion
of entropy from a molecular view.
18
3.9
Before calculating the change in entropy, we need to determine the final state of the system. Let
the mixture of ice and water immediately after the ice has been added represent the system.
Since the glass is adiabatic, no work is performed, and the potential and kinetic energies are
neglected, the energy balance reduces
H 0
We can split the system into two subsystems: the ice (subscript i) and the water (subscript w).
Therefore,
H mi hi mw hw 0
and
mi hi mw hw
mw
Vw
0.0004 m 3 0.399 kg
vw m3
0.001003
kg
nw
mw
0.399 kg 22.15 mol
MW H 2O 0.0180148 kg
mol
5.55 mol
mi
ni
MW H 2O
Now, lets assume that all of the ice melts in the process. (If the final answer is greater than 0
C, the assumption is correct.) The following expression mathematically represents the change
in internal energies (cP=cv).
ni c P ,i 0 10 C h fus c P ,w T f 0 C nw c P ,w T f 25 C
Note: Assumed the heat capacities are independent of temperature to obtain this
expression.
c P ,i 4.196 R
c P ,w 9.069 R
and
kJ
h fus 6.0
mol
19
Substitution of values into the above energy balance allows calculation of Tf.
T f 3.12 C
s surr 0
suniv s sys
We will again break the system into two subsystems: the ice and the water. The change in
entropy of the universe can be calculated as follows
S univ ni si nw s w
final
qrev
s T
initial
Assuming the heat capacities of ice and water are independent of temperature, the expressions
for the change in entropy of the subsystems are
Therefore
20
J
Suniv 6.59
mol
21
3.10
(a)
The maximum amount of work is obtained in a reversible process. We also know the entropy
change for the universe is zero for reversible processes. From the steam tables
kJ
s1 6.2119 (400 C, 100 bar)
kg K
kJ
s2 8.5434 (400 C, 1 bar)
kg K
Since the process is isothermal, we know that the temperature of the surroundings is constant at
400 C. Therefore,
Qsurr kJ
S surr 1.1658
Tsurr K
Qsurr 784.76 kJ
To find the work obtained, perform an energy balance. An energy balance where potential and
kinetic energy effects are neglected is
U Q W
We can calculate the change in internal energy from the steam tables, and we also know that
Q Qsurr . From the steam tables:
kJ
u1 2832.4 (400 C, 100 bar)
kg
kJ
u 2 2967.8 (400 C, 1 bar)
kg
Therefore,
22
W mu 2 u1 Q
kJ kJ
W 0 . 5 kg 2967 . 8 2832 . 4 784 . 76 kJ
kg kg
W 717.06 kJ
(b)
If steam is modeled as an ideal gas, the change in internal energy is zero. Therefore, the energy
balance is
Q W
The work can be found using Equation 2.77 which is developed for reversible, isothermal
processes.
P
W nRT ln 2
P1
W
0.5 kg J
8.3145
1 bar
673.15 K ln
kg mol K 100 bar
0.0180148
mol
W 715332 J 715.3 kJ
Therefore,
Q 715.3 kJ
and
Qsurr 715.3 kJ
kJ
S surr 1.06
K
23
3.11
(a)
The initial pressure is calculated as follows
P1 Psurr
25000 kg 9.81 m/s 2 24.6 105 Pa
2
0.05 m
P2 Psurr
35000 kg 9.81 m/s 2 34.4 105 Pa
0.05 m 2
(b)
The temperature should rise, which can be understood by considering an energy balance.
Because the system is insulated, the work done by the mass being added to the piston is
transformed into molecular kinetic energy.
(c)
The final temperature can be calculated with an energy balance:
U W
ncv T2 T1 W
Since the pressure of the surroundings is constant after the block is added to the piston, the work
is calculated as follows:
W P2 V2 V1
nRT2 nRT1
V2 V1 0.00169 m 3
P2 P1
nRT2
n5 / 2 R T2 T1 P2 V1
P2
T2 557 K
24
(d)
Since the system is well-insulated
s surr 0
T V T T P
suniv s sys cv ln 2 R ln 2 5 / 2 R ln 2 R ln 2 1
T1 V1 T1 T1P2
J
s sys 0.354
mol K
(e)
Because the change in entropy of the universe is equal to change in entropy of the system, which
is positive in this situation, the second law is not violated.
25
3.12
(a)
To calculate ssys, we need to take the gas from state 1 to state 2 using a reversible process. The
process in part a can be drawn as follows:
rev.
State 1 adiabatic State 2
rev.
isothermal
PI
T I = 557 K
State I
We need to find the intermediate pressure, PI where we end up at the temperature in state 2, T2.
To find it, we can use the results from pages 78-79 of the text:
Therefore,
k
T 1k
PI 1 P1 35.9 [bar] (I)
T2
20
490 560
T
26
The change in entropy can be represented as follows:
For the reversible adiabatic process, the change in entropy is zero. Therefore,
s sys sisothermal
RT
q rev wrev dP
P
Therefore,
q P
final P2
R J
s sys
initial
T
PI
P
dP R ln 2 0.354
PI mol K
If we substitute relation (I) in the expression above, we get the same expression in the book:
P k T T P
s sys R ln 2 R ln 2 c P ln 2 R ln 2
P1 1 k T1 T1 P1
(b)
For this construction, we use the path diagrammed on the following PT diagram:
40
P
Isothermal
compression
Isobaric heating
20
490 560
T
We write
27
s sys sisobar sisothermal
q rev dh c P dT
Therefore,
T2
c T
sisobar TP dT c P ln T12
T1
Therefore,
v 2 RT2 / P2 v 2 RT2 / P2
P R P P
sisothermal T
dv v
dv R ln 1 R ln 2
P2 P1
v I RT2 / P1 v1 RT2 / P1
T P J
s sys c P ln 2 R ln 2 0.354
T1 P1 mol K
(c)
For this construction, we use the path diagrammed on the following PT diagram:
28
40
P
Isothermal
compression
Isochoric heating I
20
490 560
T
For isochoric heating followed by an isothermal expansion, the entropy can be expressed as
follows:
For the reversible, isothermal expansion, we obtain (refer to Parts (a) and (b) to see how this is
derived)
P
sisothermal R ln 2
PI
However, we can relate the intermediate pressure to the initial pressure through the ideal gas law:
PI P1
T2 T1
or
PI 27.4 [bar]
P2 P2 T P
sisothermal R ln R ln R ln 2 R ln 2
PI T2 T1 P1
P1
T1
q rev cv dT
Therefore,
29
T2
cv T
sisochoric T dT cv ln T12
T1
T P T P J
s sys cv R ln 2 R ln 2 c P ln 2 R ln 2 0.354
T1 P1 T1 P1 mol K
Using any of these three reversible paths, we get the same answer!
30
3.13
(a)
Since the process is reversible and adiabatic, the entropy change for the process is zero.
s sys 0
Furthermore, since the process is adiabatic, the changes in entropy of the surroundings and the
universe are also zero.
(b)
Since this process is isentropic (s=0), we can apply an expression for the entropy change of an
ideal gas.
T2
c P
s 0 P dT Rln 2
T T 1
P1
s 0 R dT R ln 2
T1
T P1
T1 250 K
P1 1 bar
P2 12.06 bar
yields
T2 482 K
(c)
An energy balance gives
T2 T2
w u cv dT
T1
c
T1
P 1dT
31
w R 2.639 0.506 10 3 T 0.227 10 6 T 2 dT 5390
T2
J
T1 mol
(d)
Reversible processes represent the situation where the minimum amount of energy is required for
compression. If the process were irreversible, more work is required, and since the process is
adiabatic, the change in internal energy is greater. Since the change in internal energy is greater,
so too is the change in temperature. Therefore, the final temperature would be higher than the
temperature calculated in Part (b).
32
3.14
The problem statement states that the vessel is insulated, so we can assume that heat transfer to
the surroundings is negligible. Therefore, the expression for the entropy change of the universe
is
suniv s sys
An energy balance will help us solve for the entropy change. Neglecting potential and kinetic
energy effects, the energy balance is
U Q W
Since the vessel is insulated, the heat term is zero. Furthermore, no work is done, so the energy
balance is
U mu 2 u1 0
u1 u 2
The values for the initial pressure and temperature constrain the value of specific energy at state
1. From the steam tables,
kJ
u1 2619.2 (400 C, 200 bar)
kg
kJ
u 2 2619.2
kg
The problem statement provides the pressure of state 2 (100 bar). Since we know the pressure
and internal energy at state 2, the entropy is constrained. Information given in the problem
statement also constrains state 1. From the steam tables,
kJ
s1 5.5539 (400 C, 200 bar)
kg K
kJ kJ
s2 5.7754 100 bar, u 2 2619.2
kg K kg
Therefore,
kJ
suniv ssys s2 s1 0.2215
kg K
kJ
S univ S sys 0.2215
K
33
3.15
The subscript 2 refers to the final state of the system. 1 refers to the gas initially on one side
of the partition. If you take the system to be the entire tank (both sides of the partition), then no
net work is performed as the gas leaks through the hole. Furthermore, the tank is well insulated,
and kinetic and potential energy effects can be neglected. Thus, the energy balance is
u 0
and
T2 T1 (ideal gas)
Initially, the partition is divided into two equal parts. The gas fills the entire volume for the final
state ( V2 2V1 ). The ideal gas law can be used to calculate the moles of gas present.
n
P1V1
10 105 Pa 0.5 m 3 200.5 mol
RT1 J
8.314 300 K
mol K
Now, the change in entropy can be solved using a slight modification of Equation 3.65:
T v 3 T nv
S sys n cv ln 2 R ln 2 n R ln 2 R ln 2
T1 v1 2 T1 nv1
J 3
S sys 200.5 mol 8.314 ln1 ln2
mol K 2
J
S sys 1155
K
34
3.16
N2 N2 O2 N2 N2
N2
Mixing
Process
N2 N2 O2
N2
N2 O2 N2
N2 N2
N2 N2 O2 N2
N2
The tank is insulated. The change in entropy of the universe can be rewritten as
since the tank is well-insulated. After the partition ruptures, the pressure and temperature will
remain constant at 1 bar and 298 K, respectively. This can be shown by employing mass and
energy balances. Therefore, the partial pressures in the system are
298 K
cP P2 J 0.79 bar
s sys, N 2 T dT R ln 8.314
P1 N
ln
mol K 1 bar
298 K 2
J
s sys, N 2 1.96
mol K
298 K
cP P J 0.21 bar
s sys, O2 T
dT R ln 2 8.314
P1 O
ln
mol K 1 bar
298 K 2
J
s sys ,O2 12.98
mol K
Therefore,
J J
S univ 0.79 mol1.96 0.21 mol12.98
mol K mol K
35
J
S univ 4.27
K
N2 N2 N2
N2 N2 N2
N2 N2 O2 N2 N2
N2 Mixing
N2 N2 Process N2
N2 N2
N2 N2
N2 N2 O2
N2 O2 O2 N2
N2 N2
N2 N2 N2 N2 N2 N2
N2 N2
Before calculating the change in entropy, we need to find how the temperature and pressure
change during the process. The energy balance simplifies to
U 0
nO 2
n N 2 cvT2 n N 2 cv 298 K nO2 cv 298 K 0
T2 298 K
We can find the pressure after the rupture by recognizing that the tank is rigid. Therefore,
VO2 V N 2 Vtot
By employing the ideal gas law, we get the following equation (it has been simplified):
nN 2
nO2
nO 2
nN 2
P1, N 2 P1, O2 P2
P2 1.65 bar
Now, use Equation 3.62 to calculate the entropies as was done in Part (a):
36
298 K
cP P J 1.30 bar
s sys, N 2
T
dT R ln 2 8.314
P1 N
ln
mol K 2 bar
298 K 2
J
s sys, N 2 3.58
mol K
298 K
cP P J 0.347 bar
ssys ,O2
298 K
T
dT R ln 2 8.314
P1 O2
ln
mol K 1 bar
J
s sys ,O2 8.8
mol K
Therefore,
J J
S univ 0.79 mol 3.58 0.21mol 8.8
mol K mol K
J
Suniv 4.68
K
37
3.17
First start with the energy balance for the throttle. Potential and kinetic energy effects can be
neglected. During the throttling process, no shaft work is performed and the rate of heat transfer
is negligible. Therefore,
m 2 h2 m 1h1
h2 h1
kJ
h1 3398.3 (500 C, 8 MPa)
kg
kJ
h2 3398.3
kg
kJ
ssurr 0
kg K
The steam tables can be used to calculate the change in entropy of the system.
kJ kJ
s2 8.7098 100 kPa, h2 3398.3
kg K kg
kJ
s1 6.7239 (8 MPa, 500 C)
kg K
Thus,
kJ
suniv ssys s2 s1 1.9859
kg K
38
3.18
For this process to work, conservation of mass and the first and second laws of thermodynamics
must hold. The subscript 1 refers to the inlet stream, 2 refers to the cold outlet, and 3
refers to the hot outlet. To test the conservation of mass, perform a mass balance
m 1 m 2 m 3
kg kg kg
2 0.5 1.5
s s s
Clearly, the conservation of mass holds. Now test the first law of thermodynamics by writing an
energy balance. Since there is no heat transfer or work, the energy balance becomes
Using the conservation of mass we can rewrite the mass flow rate of stream 1 in terms of streams
2 and 3:
m 2 h2 h1 m 3 h3 h1 0
If we assume that the heat capacity of the ideal gas is constant, the equation can be written as
follows:
m 2 cP T2 T1 m 3cP T3 T1 0
Therefore,
kg kg
0.5 cP (60 K) 1.5 cP (20 K) 0
s s
This proves that the first law holds for this system. For the second law to be valid, the rate of
change in entropy of the universe must greater than or equal to 0, i.e.,
dS
0
dt univ
Assuming the process is adiabatic, we can write the rate of entropy change of the universe for
this steady-state process using Equations 3.48-3.50:
dS
m 2 s2 m 3 s3 m 1s1 m 2 s2 s1 m 3 s3 s1
dt univ
39
where a mass balance was used. For constant heat capacity, we can calculate the entropy
differences using Equation 3.63:
T P 5 T P
s2 s1 cP ln 2 R ln 2 R ln 2 ln 2
T1 P1 2 T1 P1
T P 5 T P
s3 s1 cP ln 3 R ln 3 R ln 3 ln 3
T1 P1 2 T1 P1
dS
1
2.77 R 0.057cP
dt univ MW
Therefore, the second law holds if cP 48.6 R ; this value is far in excess of heat capacity for
gases, so this process is possible.
40
3.19
(a)
T1 = 640 oC
P1 = 4 MPa P2 = 100 kPa
J
S sys 0
K
(c)
The steam tables can be used to determine the final temperature. From the interpolation of steam
table data
kJ
s2 s1 640 C, 4 MPa 7.4692
kg K
Now two thermodynamic properties are known for state 2. From the steam tables:
kJ
T2 121.4 C P2 0.1 MPa, s2 7.4692
kg K
(d)
An energy balance is required to calculate the exit velocity. The energy balance for the nozzle is
0 m 2 h eK 2 m 1 h e K 1
eK ,2 h e K 1 h2
which is equivalent to
V2 2 h 0.5V 2 1 h2
Substituting the following values
41
m
V1 20
s
J
h2 2719100 (0.1 MPa, 121.4 C)
kg
J
h1 3767000 (4 MPa, 640 C)
kg
yields
m
V2 1448
s
Note: the velocity obtained is supersonic; however, the solution does not account for this type of
flow.
42
3.20
A schematic of the process follows:
Propane in
Since this process is isentropic (s=0), we can apply an expression for the entropy change of an
ideal gas. We must be careful, however, to select an expression which does not assume a
constant heat capacity.
T2
c P
s 0 P dT Rln 2
T T 1
P1
s 0 R dT R ln 2
T1
T P1
We can also relate P1 to T1 and v1 (which are known) through the ideal gas law:
RT1
P1
v1
T Pv
s 0 1.213 ln 2 28.785 10 3 (T2 T1 ) 4.412 10 6 T22 T12 ln 2 1
T1 RT1
Using
T1 = 623 K,
v1 = 600 cm3/mol,
P2 = 1 atm
43
R=82.06 cm3 atm/mol K,
we get:
T2 = 454 K
W S
To solve for we need a first law balance, with negligible ke and pe, the 1st law for a steady
n
state process becomes:
0 n h1 h2 Q W S
W S
T2 T2
n T1 T1
integrating:
W S
n
R 1.213T2 T1 14.393 10 3 T22 T12 2.94110 6 T23 T13
Substituting
J
R 8.314
mol K
T2 454 K
T1 623 K
results in
W S J
19860
n mol
44
3.21
A reversible process will require the minimum amount of work. For a reversible process,
suniv 0
s sys s surr
For this process, the change in entropy of the system can be calculated as follows
s sys sO2 s N 2
where
T2 c
P O2 P
sO2 T
dT R ln 2
P1
T1 O2
T2 c
P N 2 P
s N 2 T
dT R ln 2
P1
T1 N2
Since the temperature of the streams do not change, the expressions reduce to
P
sO2 R ln 2
P1 O2
P
s N 2 R ln 2
P1 N 2
Therefore,
P2 P
s surr R ln ln 2
P1 O2 P1 N 2
q
s surr surr
Tsurr
where
qsurr q
45
Therefore,
P P2
q surr Tsurr s surr RTsurr ln ln
2
P1 O P1 N 2
2
J 1 bar 1 bar
q surr 8.314
293 . 15 K ln ln
mol K 0.21 bar O2 0.79 bar N 2
J
q surr 4378.2
mol
and
J
q 4378.2
mol
Energy Balance:
n
in
h nout hout Q W S 0
in in
out
Because the temperature of the oxygen and nitrogen doesnt change in the process, the energy
balance per mole of feed becomes
wS q
J
wS 4378.2
mol
46
3.22
(a)
Take the entire container to be the system and assume no heat or work crosses the system
boundary.
Energy Balance:
U 0
After the oscillation cease, the temperature and pressure on both sides of the piston will be the
same assuming the metallic piston is very thin and the thermal conductivity coefficient is large.
Now, lets rewrite the energy balance.
or
n1, A
ugas A ugas B 0
n1,B
n1, A PA,1VA,1TB ,1 P T
A,1 B ,1 2.41
n1,B PB ,1VB ,1TA,1 2 PB ,1TA,1
since
2V A,1 VB,1
so
T
3
2.41 RT2 773.15 K R 5/2 1.5 10-3 T dT 0
2
2 373.15 K
T2 585 K
Mass Balance:
47
PA,1VA,1 PB ,1VB ,1 PA, 2VA, 2 PB , 2VB , 2 P2Vtot
RTA,1 RTB ,1 RTA, 2 RTB , 2 RT2
since the pressure and temperature of state 2 are equal on both sides. But
P2 3.56 bar
(b)
The container is well-insulated, so
suniv s sys
The entropy change of the system can be split into two subsystems:
n1, A n1,B
ssys s gas A s gas B
n1, A n1,B n1, A n1,B
2.41 T P J
sgas A 2.5 R ln 2 R ln 2 1.96
3.41 TA,1 PA,1 mol K
1 2 3.5 1.5 10 3 T P
T
J
sgas B R dT R ln 2 1.51
3.41 TB ,1
T PB ,1 mol K
Therefore,
J J J
suniv ssys 1.96 1.51 3.47
mol K mol K mol K
48
3.23
Lets first calculate S sys . Before this problem is solved, a few words must be said about the
notation used. The system was initially broken up into two parts: the constant volume container
and the constant pressure piston-cylinder assembly. The subscript 1 refers to the constant
volume container; 2 refers the piston-cylinder assembly. i" denotes the initial state before the
valve is opened, and f denotes the final state.
First, the mass of water present in each part of the system will be calculated. The mass will be
conserved during the expansion process. Since the water in the rigid tank is saturated and is in
equilibrium with the constant temperature surroundings (200 C), the waters entropy is
constrained. From the steam tables,
kJ
v1l, i 0.001156
kg
kJ
v1v, i 0.12736
kg
kJ
s1l, i 2.3308 (Sat. water at 200 C)
kg
kJ
s1v, i 6.4322
kg
P sat 1553.8 kPa
Knowledge of the quality of the water and the overall volume of the rigid container can be used
to calculate the mass present in the container.
V1 0.05m1 v1l, i 0.95m1 v1v,i
Using the values from the steam table and V1 0.5 m 3 provides
m1 4.13 kg
49
m1v 0.95m1 3.92 kg
m1l 0.05m1 0.207 kg
For the piston-cylinder assembly, both P and T are known. From the steam tables
m3
v2, i 0.35202
kg
(600 kPa, 200 C)
kJ
s2, i 6.9665
kg K
Now enough information is available to calculate the mass of water in the piston assembly.
m2 2 0.284 kg
V
v2
Now the final state of the system must be determined. It helps to consider what physically
happens when the valve is opened. The initial pressure of the rigid tank is 1553.8 kPa. When
the valve is opened, the water will rush out of the rigid tank and into the cylinder until
equilibrium is reached. Since the pressure of the surroundings is constant at 600 kPa and the
surroundings represent a large temperature bath at 200 C, the final temperature and pressure of
the entire system will match the surroundings. In other words,
kJ
s f s2, i 6.9665 (600 kPa, 200 C)
kg
kJ
S sys 3.05
K
Now we calculate the change in entropy of the surroundings. Since the temperature of the
surroundings is constant,
Qsurr Q
S surr
Tsurr Tsurr
After neglecting potential and kinetic energy effects, the energy balance becomes
50
U Q W
The change in internal energy and work will be calculated in order to solve for Q. The following
equation shows how the change in internal energy can be calculated.
kJ
u1l, i 850.64
kg
(sat. H2O at 200 C)
kJ
u1v, i 2595.3
kg
kJ
u f u 2, i 2638.9 (600 kPa, 200 C)
kg
U 541.0 kJ
Calculating the work is relatively easy since the gas is expanding against a constant pressure of
600 kPa (weight of the piston was assumed negligible). From Equation 2.7,
Vf
W PE dV PE (V f Vi )
Vi
where
PE 600000 Pa
V f (m2 m1v, i m1l,i )v2, i 1.55 m 3
Vi 0.1 m 3 0.5 m 3 0.6 m 3
Note: v2, i was used to calculate V f because the temperature and pressure are the same
for the final state of the entire system and the initial state of the piston-cylinder assembly.
W 570 kJ
51
The energy balance can be used to obtain Q.
Therefore,
1111 kJ kJ
S surr 2.35
473.15 K K
and
kJ kJ kJ
Suniv S sys S surr 3.05 2.35 0.70
K K K
52
3.24
Let the subscript 1 represent the state of the system before the partition is removed. It has two
components: component a and component b. Subscript 2 represents the system after the
partition is removed.
Mass balance:
n1, a n1, b n2
U sys W
n1, a cv, a T2 T1, a n1, b cv, b T2 T1, b P V2 V1, a V1, b
The external pressure, P, for the above energy balance is equal to P1,a. To find the final
temperature, first find the volumes using the ideal gas law.
2 mol 8.314 J
300 K
mol K
V1, a 0.05 m 3
1000 kg 9.8 m s - 2
0.098 m 2
4 mol 8.314 J
T2
mol K
3.3256 10 4 T2 m 3
V2
1000 kg 9.8 m s
- 2
0.098 m 2
Substitute the volumes into the energy balance and solve for T2 :
T P T P
S sys na c P ln 2 R ln 2 nb c P ln 2 R ln 2
T1 P1 a T1 P1 b
The only unknown in the above equation is P1,b, so we can calculate it with the ideal gas law:
53
2 mol 8.314 J
300 K
mol K
P1, b 0.5 10 5 Pa
3
0 .1 m
54
3.25
First perform an energy balance on the process:
U sys W
U n A, f u A, f nB , f u B , f n A,i u
n A,i n A , f nB , f
U n A, f cv TA, f Ti nB , f cv TB , f Ti W
5
Using the ideal gas law and cv R , we obtain
2
5 5 5
PB, f VB, f PA, f V A, f PA, iV A, i W
2 2 2
Vf
W PE dV
Vi
kx mg k VB VB ,i mg
PE Patm Patm
A A A2 A
PE 1.06 10 8.284 10 VB Pa
5 5
Integrating
W 1.06 105VB, f 4.142 105VB2, f
5 5 5
PB, f VB, f PA, f V A, f PA,iV A, i 1.06 105VB, f 4.142 105VB2, f
2 2 2
55
This expression can be simplified by recognizing that PB, f PA, f and V A, f V A, i :
5
2
5
PB, f VB, f V A, i PA, iV A, i 1.06 105VB, f 4.142 105VB2, f
2
We still have two unknowns: PB, f and VB, f . We can eliminate PB, f by writing a force
balance for the final state:
5
2
5
1.06 105 8.284 105VB, f VB, f V A, i PA, iV A, i 1.06 105VB, f 4.142 105VB2, f
2
VB, f 0.268 m 3
S sys, A 0
Therefore,
T P
0 c P ln 2 R ln 2
T1 P1
J 7 T2 3.28 10 5 Pa
0 8.314 ln ln
mol K 2 313.15 K 7.0 10 5 Pa
T2 252 K
56
3.26
A schematic of the process is drawn:
Q surr=69.86 kW
Steam
in W
P1 = 60 bar Turbine
T1 = 500 oC Steam
out P2 = 10 bar
Tsurr=300 K
dS dS dS
0
dt univ dt sys dt surr
or
Q surr
m s2 s1 0
Tsurr
kJ
We can look up property values of state 1 from the steam tables s1 6.8802 and
kgK
kJ
h1 3,422.13 . Converting the units of mass flow rate gives:
kg
kg 1 hr kg
m 4,500 1.25
hr 3,600 s s
so
Q surr kJ
s2 s1 6.6939
m Tsurr kgK
We now know two properties of steam ( s2 and P2) From the steam tables:
T2 200 oC
Using the steam tables for superheated steam at 1 MPa, we find that when water has this value of
entropy
57
kJ
h2 2827.9
kg
Energy balance:
m h2 h1 Q W s
Calculate Ws:
kJ kJ
W s m h2 h1 Q 1.25 kg/s 2827.9 3422.1 69.86 kW
kg kg
W s 673 kW
(c)
The isentropic efficiency is defined as follows:
W S actual
W S reversible
For this situation,
(e)
Use the energy balance:
m h2 h1 Q W s
Q W s 69.86 kW -447.5 kW kJ kJ
h2 h1 3422.1 3008.2
m 1.25 kg/s kg kg
T2 actual 280.2 C
58
3.27
Since the compressor is adiabatic, the energy balance after neglecting potential and kinetic
energy becomes
n h2 h1 W S
Using the ideal gas law and Appendix A.2, the above equation becomes
T
P1V1 2
A BT CT 2 DT 2 ET 3 dT
T1
WS
T1
P1 1 bar 1 105 Pa
m3
V1 1
s
T1 293.15 K
T2 473.15 K
A 3.355; B 0.575 10 3 ; C 0; D 0.016 105 ; E 0 (Table A.2.2)
gives
W S 218.8 kW
This is the work of our real turbine (with 80% isentropic efficiency). We can use the isentropic
efficiency to calculate the work of an equivalent 100% efficient process.
W
compressor S reversible
WS compressor
W S reversible W S compressor compressor 0.8218.8 kW
W S reversible 175.4 kW
Since the process is adiabatic, we can use the following equation again to calculate what the final
temperature would be in a reversible process.
T
P1V1 2
WS A BT CT 2 DT 2 ET 3 dT
T1
T1
59
A BT CT
T2
P V
175400 W 1 1 2
DT 2 ET 3 dT
T1
T1
To obtain the pressure the final pressure for the real compressor, we can calculate the final
pressure for the reversible process because the final pressure is the same in both cases. For the
isentropic expansion
T2
cp P
s 0 T
dT R ln 2
P1
T1
A BT CT 2 DT 2 ET 3
T2
P
s 0 R dT R ln 2
T1
T P1
T2 438 K
s 0 R
3.355 0.575 10 3 T 1600T 2
dT ln
P2
1 bar
T
1
T 293 .15 K
P2 4.16 bar
60
3.28
Isentropic efficiency for a turbine is defined as
ws actual
ws rev
If the rate of heat transfer is assumed negligible, the energy balance for this process is
0 m h1 h2 W s
s sys 0
s2,rev s1 500 C, 10 MPa
kJ
s2,rev s1 6.5965
kg K
Since s2l , sat s2 s2v, sat , a mixture of liquid and vapor exists. The quality of the water can be
calculated as follows
Therefore,
kJ
h1 3373.6 (500 C, 10 MPa)
kg
61
The reversible work is calculated as:
W S rev h kJ kJ kJ
h1 2391.0 3373.6 982.6
2
m 1 kg kg kg
W S actual
W S rev 835 kJ
m 1 m 1 kg
The exit temperature is calculated by determining the enthalpy of the actual exit state:
W S actual kJ
2, actual h1 2538.4
h
m 1 kg
This state is still saturated (although the quality is higher), so the temperature is
T2 = 99.6 oC
62
3.29
The subscripts 2 and 3 represent the two outline streams, and 1 represents the inlet stream.
First, perform a mass balance:
2 1
n1 n 2 n3 n1 n1
3 3
0 n1h1 n2 h2 n3 h3 Q W S
Since the system is insulated and there is no shaft work, the energy balance can be rewritten as:
2 1
n1h2 n1h3 n1h1 0
3 3
2
h2 h1 1 h3 h1 0
3 3
Substituting expressions for heat capacity A.2.2, we obtain the following expression:
400 K T
2 3 4000 1 3
3 4000
3
R 3.280 0.593 10 T T 2 dT 3 R 3.280 0.593 10 T T 2 dT 0
300 K 300 K
Solving, we find
T3 100.6 K
Now set up an entropy balance for the process. The minimum pressure is required for a
reversible process.
2 1
n1s2 n1s3 n1s1 0
3 3
2
s2 s1 1 s3 s1 0
3 3
Assuming ideal gas behavior, we can express the changes in entropies using Equation 3.62:
63
Substitute these expressions into the entropy balance and solve for P1:
P1 1.85 bar
64
3.30
A schematic of the process is illustrated below:
reversible Tsurr = 25 oC
Psurr = 1 bar
Pure H 2O
P2 = 20 bar
Pure H2O
T 2 = ? oC
P1 = 60 bar
m = 5 kg
T1 = 540 oC
m = 5 kg
well-insulated
(a)
Since the process is adiabatic and reversible,
suniv 0
s sys 0
s surr 0
(b)
We can obtain the final temperature using the steam tables.
kJ
s1 6.999 (540 C, 60 bar)
kg K
kJ
s2 6.999 (20 bar)
kg K
The pressure and entropy of state 2 can be used to back out T2 from the steam tables.
T2 362.5 C
(c)
To obtain the value of work, perform an energy balance. The process is adiabatic, and potential
and kinetic energy effects can be neglected. Therefore, the energy balance is
U mu 2 u1 W
kJ
u2 2881.15 (362.5 C, 20 bar)
kg
65
kJ
u1 3156.12 (540 C, 60 bar)
kg
Hence,
kJ kJ
W 5 kg 2881.15 3156.12 1373 kJ
kg kg
(d)
The specific volume can be found in the steam tables
m3
v2 0.14173
kg
Therefore,
m3
V2 5 kg 0.14173
0.709 m
3
kg
66
3.31
A schematic of the process is shown below:
very tiny Tsurr = 25 oC
leak hole
Psurr = 1 bar
P1 = 60 bar P2 = 20 bar
T1 = 540 oC T ' 2 = ? oC
m = 5 kg
well-insulated
In order to leave the system, the gas must do flow work on the surroundings. The initial state is
the same as for Problem 3.30 and the final pressures are the same. Since the water only expands
against 1 bar, the work is lower than that for the differential process described in Problem 3.30.
Thus, this adiabatic process looses less energy, leading to a higher final temperature.
Another way to view this argument is to look at this process as a closed system. This depiction
is the expansion analog of the compression process depicted for Example 2.5 in Figure E2.5B
(page 57). We can represent this process in terms of two latches, one that keeps the process in its
initial state at 60 bar and one that stops the expansion after the pressure has reduced to 20 bar.
The process is initiated by removal of the first latch and ends when the piston comes to rest
against the second latch. Such a process is depicted as Problem 3.31 below. The
corresponding reversible process of Problem 3.30 is shown next to it for comparison. Clearly the
process on the left does less work, resulting in a greater final temperature.
T surr = 25 oC
Psurr = 1 bar reversible Tsurr = 25 oC
Psurr = 1 bar
latches
Pure H2O
Pure H 2O
Pure H2O irreversible P2 = 20 bar P2 = 20 bar
Pure H2O
P1 = 60 bar T ' 2 = ? oC T 2 = ? oC
P1 = 60 bar
T1 = 540 oC m = 5 kg m = 5 kg
T1 = 540 oC
m = 5 kg
m = 5 kg
well-insulated
67
3.32
(a)
Consider the tank as the system.
Mass balance
dm
m in m out m in
dt
dm m in dt
m1 0
or
t
m2 m1 m in dt
0
Energy balance
Since the potential and kinetic energy effects can be neglected, the open system, unsteady state
energy balance is
dU
m out hout m in hin Q W s
dt sys out in
The process is adiabatic and no shaft work is done. Furthermore, there is only one inlet stream
and not outlet stream. Therefore, the energy balance simplifies to
dU
m in hin
dt sys
U2 t t
dU m
U1 0 0
h dt hin m in dt
in in
0
U 2 m2u 2 m2 hin
u 2 hin
68
From the steam tables,
kJ
hin 3456.5 (3 MPa, 773 K)
kg
Now the water in the tank is constrained. From the steam tables:
kJ
s2 7.743
kg K kJ
3 MPa, u 2 3456.5
m 3
kg
v2 0.14749
kg
dS dS
m in sin
dt univ dt sys
t
S univ m2 s2 sin m in dt m2 s2 sin
0
kJ
sin 7.2337 (3 MPa, 773 K)
kg K
Therefore,
0.05 m 3 kJ kJ
0.173
kJ
S univ 7.743
7.2337
m 3 kg K kg K K
0.14749
kg
(b)
If it tank sits in storage for a long time and equilibrates to a final temperature of 20 C, some or
all of the vapor will condense and exchange heat with the surroundings. Let the subscript 3
designate the final state of the water when it has reached a temperature of 20 C. The change in
entropy between state 2 and state 3 is given by the following equation
69
S univ m2 s3 s2
Qsurr
T
Let the subscript 3 designate the final state of the water when it has reached a temperature of
20 C. First, find quality of the water. From the steam tables:
m3
v3sat
, l 0.001002 kg
(sat. water at 20 C)
m3
v3sat
,v 57 . 79
kg
v3 v2 1 x v3sat sat
, l xv3, v
m3 m3 m3
0.14749 1 x 0.001002 x 57.79
kg kg kg
x 0.0025
s3 1 x s3sat sat
, l xs3, v
kJ kJ
s3 1 0.0025 0.2966
kg K 0.0025
8.6671 kg K
kJ
s3 0.3175
kg K
Q surr Q m2 u 3 u 2
u 3 1 x u 3sat sat
,l xu 3,v
70
kJ
u 3 89.74
kg
Therefore,
0.05 m 3 kJ kJ
Q surr 89.74 3456.5 1141 kJ
m 3 kg
kg
0 . 14749
kg
0.05 m 3
kJ kJ 1141 kJ
S univ kg K kg K 293 K
3
0 .3175 7 . 743
0.14749
m
kg
kJ
S univ 1.38
K
The entropy change for both processes can be fount by adding together the entropy change from
Part (a) and Part (b):
71
3.33
A schematic is given below
valve maintains
pressure in system
constant
v
T1 = 200 oC
x1 = 0.4
V = 0.01 m3
Mass balance
dm
m in m out m out
dt
dm m
m1 0
out dt
or
t
m2 m1 m out dt
0
Entropy Balance:
dS Q dS Q
Suniv m out sout surr m out sout
dt sys Tsurr dt sys Tsurr
72
kJ
s1 (1 x)u f xu g 0.6 2.3309 0.4 6.4323 3.9715
kg K
kJ kJ
s 2 6.4323 is the same as sout 6.4323 s2
kg K kg K
Q
Suniv m1 ( s 2 s1 ) (2)
Tsurr
dU
m out hout Q
dt sys
Integrating
dU
m 2 u 2 t t t
m out hout Q dt hout m out dt Q dt
m1u1 0 0 0
Q m 2 u 2 m 1 u 1 m 2 m 1 h out
m3
v1 (1 x)v f xv g 0.6 .001 0.4 0.1274 0.051
kg
m3
v2 0.1274
kg
V
m1 1
0.01 m 3
0.196 kg
V
and m2 2
0.01 m 3
0.0785 kg
v1 m3 v2 m3
0.051 0.1274
kg kg
73
kJ
u1 (1 x)u f xu g 0.6 850.64 0.4 2597.5 1549
kg
kJ kJ
u 2 2595.3 and hout 2793.2
kg kg
Plugging in:
74
3.34
The process can be represented as:
V1 = 1 m3 Process V2 = 10 m3
P1 = 10 bar P2 = ?
T1 = 1000 K Work out T2 = ?
Q=0
State 1 State 2
Solving for the number of moles:
PV
n 1 1 120.3 mol
RT1
The maximum work is given by a reversible process. Since it is also adiabatic, the entropy
change of the system is zero:
T P
s 0 cP ln 2 R ln 2
T1 P1
5
Since cP R
2
5/ 2
T2 P2 nRT2
T1 P1 V2 P1
2/3
nRT 5 / 2
T2 1
215 K
V2 P1
T2 T2
U = Q W W n c v dT n cP R dT
T1 T1
Solving for work, we get
3
W n RT2 T1 1.18 10 J
6
2
75
3.35
The maximum efficiency is obtained from a Carnot cycle. From Equation 3.32
T
n 1 C
TH
298.15 K
n 1 0.614
773.15 K
76
3.36
The labeling described in Figure 3.8 will be used for this solution. First, a summary of known
variables is provided.
From our knowledge of ideal Rankine cycles, the table can be expanded as follows
The saturation condition constrains state 3. First, start with the turbine. Since the rate of heat
transfer is negligible and the expansion occurs reversibly in an ideal Rankine cycle, the
following is known
s1 s2
and
W s m h2 h1
From the steam tables:
kJ
s1 7.2337 (500 C, 30 bar)
kg K
kJ
s2l 0.6492
kg K
(sat. H2O at 0.1 bar)
kJ
s2v 8.1501
kg K
kJ
h1 3456.5 (500 C, 30 bar)
kg
kJ
h2l 191.81
kg
(sat. H2O at 0.1 bar)
kJ
h2v 2584.6
kg
77
The quality of the water can be calculated as follows
s1 (1 x) s2l xs2v
x 0.878
h2 (1 x)h2l xh2v
Substituting the steam table data and the quality calculated above, we get
kJ
h2 2292.7
kg
Therefore,
kg kJ kJ
W s 100 2292.7 3456.5 116.38 MW
s kg kg
At state 3,
kJ
h3 191.81
kg
(sat. H2O(l) at 0.1 bar)
m3
v3 0.001010
kg
Therefore,
Q C m h3 h2 210.09 MW
Since the molar volume of water does not change noticeably with pressure and the compressor is
adiabatic, the work required for compressing the fluid can be calculated as follows
W c m v3 P4 P3
kg m3
W c 100 0.00101 30 105 0.1 105 Pa 3.02 10 5 W 0.302 MW
s kg
78
Now we find the rate of heat transfer for the evaporator. For the entire Rankine cycle, Equation
3.84 gives
Q net W net Q C Q H W s WC 0
We have
Q H W s WC Q C
Q H 326.2 MW
Equation 3.84:
W net W S WC
W net 116.08 MW
W 116.08 MW
net 0.356
QH 326.02 MW
79
3.37
To make the Rankine cycle more efficient, we need to increase the area that represents the net
work in Figure 3.8. This can be done in a variety of ways:
1. Increase the degree of superheating of steam in the boiler. This process is sketched in the
upper left hand Ts diagram below. This change reduces the moisture content of steam
leaving the turbine. This effect is desirable since it will prolong the life of the turbine;
however, if the steam is heated too high, materials limitations of the turbine may need to be
considered.
2. Lower the condenser pressure. Lowering the pressure in the condenser will lower the
corresponding saturation temperature. This change will enlarge the area on the Ts diagram,
as shown on the upper right below. Thus, we may want to consider lowering the pressure
below atmospheric pressure. We can achieve this since the fluid operates in a closed loop.
However, we are limited by the low temperature off the heat sink that is available.
increase the superheating of steam in the boiler lower the condenser pressure
T T
1new
1
wnet wnet
4 4new
2new
3
3new 2new
s s
1 Turbine 2
T T
1 3
6
2
wnet wnet
4
4
3 2 5
s s
80
3. Increase the boiler pressure. This will increase the boiler temperature which will increase the
area as shown on the bottom left Ts diagram.
4. We can be more creative about how we use the energy available in the boiler. One way is to
divide the turbine into two stages, Turbine 1 and Turbine 2, and reheat the water in the boiler
between the two stages. This process is illustrated below.
Rankine cycle with reheat
Ws Ws
turbine Tur-
1 bine
1 2
2
Fuel
Air
4
QH 3
Boiler QC cooling
water
6
condenser
5
Wc
compressor
The pressure in Turbine 1 will be higher than the pressure in turbine 2. This process is
schematically shown on the bottom right Ts diagram. This process leads to less moisture
content at the turbine exit (desirable) and limits the temperature of the superheat (desirable).
81
3.38
The labeling described in Figure 3.8 will be used for this solution. First, a summary of known
variables is provided.
From our knowledge of ideal Rankine cycles, the table can be expanded as follows
The saturation condition constrains state 3. First start with the turbine (state 2). Since the rate of
heat transfer is negligible and the expansion occurs reversibly in an ideal Rankine cycle, the
following is known
s1 s2
kJ
s1 6.5965 (500 C, 100 bar)
kg K
kJ
s2l 1.3025
kg K
(sat. H2O at 1 bar)
kJ
s2v 7.3593
kg K
s1 (1 x) s2l xs2v
x 0.874
82
States 1, 2, and 3 are constrained; therefore, the enthalpies can be determined using the steam
tables. The enthalpy of state 4 can be calculated from Equation 3.80:
h4 h3 vl P4 P3
m3
where vl v3 0.001043
since specific volume of liquids are insensitive to pressure
kg
changes. From the steam tables:
State 1 2 3 4
h kJ kJ kJ kJ
3240.8 2391.0 417.44 427.8
kg kg kg kg
Note: h2 (1 x)h2l xh2v
Equation 3.84:
W net Q H Q C
Therefore,
W net m h1 h4 m h3 h2 m h1 h4 h3 h2
and
Now consider the non-ideal turbine and compressor. Use the definition of isentropic
efficiencies:
turbine
W s actual
Ws reversible
W S actual turbine W S reversible turbine m h2 h1 reversible
83
compressor
WC reversible
WC actual
WC reversible m h4 h3 reversible
WC actual
compressor compressor
W net W s WC
W net m turbine h2 h1
h4 h3 reversible
reversible compressor
Substituting the enthalpy values from the table shown above and the efficiency values gives
kg
m 147
s
84
3.39
The labeling shown in Figure 3.8 will be used for this solution. First, a summary of known
variables is provided.
From our knowledge of ideal Rankine cycles, the table can be expanded as follows
The enthalpy of states 1 and 3 can be found using the NIST website. For states 2 and 4, use the
following expressions to find the enthalpies
h2 (1 x)h2l xh2v
h4 h3 vl P4 P3 h3 v3 P4 P3
s1 s2
s1 (1 x) s2l xs2v
J
s1 146.38 (sat. vapor at 1.7 MPa)
mol K
J
s2l 90.45
mol K
(sat. mixture at 0.7 MPa)
J
s2v 150.07
mol K
85
x 0.938
J
h1 36046 (sat. vapor at 1.7 MPa)
mol
J
h2l 18416
mol
(sat. mixture at 0.7 MPa)
J
h2v 35353
mol
J
h3 18416
mol
(sat. liquid at 0.7 MPa)
5 m
3
v3 6.954 10
mol
J
h2 34302.9
mol
J
h4 18485.5
mol
h2 h1 h4 h3
rankine
h1 h4
J J J J
34302.9 36046 18485.5 18416
mol mol mol mol
rankine
J J
36046 18485.5
mol mol
rankine 0.095
86
3.40
(a)
Using the information in the problem statement and our knowledge of ideal vapor compression
cycles, the following table can be created
Q C n h2 h1
Q C n h2 h4
kJ
h2 39.295 (sat. vapor at 0.12 MPa)
mol
kJ
h4 24.181 (sat. liquid at 0.7 MPa)
mol
Therefore,
mol kJ kJ
Q C 0.5
s
39.295 mol 24.181 mol 7.557 kW
(b)
The power input to the compressor can be calculated with Equation 3.86:
WC n h3 h2
s 2 s3
87
From the NIST website:
J
s3 s2 sat. vapor at 0.12 MPa 177.89
mol K
kJ J
h3 43.031 0.7 MPa, s3 177.89
mol mol K
Therefore,
mol kJ kJ
WC 0.5 43.031 39.295 1.87 kW
s mol mol
(c)
Equation 3.90 states:
kJ kJ
39.295 24.181
h h h h
COP 2 1 2 4 mol mol
4.05
h3 h2 h3 h2 kJ kJ
43.031 39.295
mol mol
88
3.41
Using the information in the problem statement and our knowledge of ideal vapor compression
cycles, the following table can be created
Q C h2 h1
COP
W C h3 h2
The enthalpy of state 2 can be found directly from the NIST website, but the enthalpies of states
1 and 3 require the use of additional information.
kJ
h2 39.295 (sat. vapor at 0.12 MPa)
mol
s 2 s3
J
s3 s 2 177.89 (sat. vapor at 0.12 MPa)
mol K
kJ J
h3 43.031 0.7 MPa, s3 177.89
mol mol K
J
s1 s 4 115.09 (sat. liquid at 0.7 MPa)
mol K
89
s1 1 x s1l xs1v
where
J
s1l 90.649 (sat. R134a(l) at 0.12 MPa)
mol K
J
s1v 177.89 (sat. R134a(v) at 0.12 MPa)
mol K
Therefore,
J J J
115.09 1 x 90.649 x177.89
mol K mol K mol K
x 0.280
Using the quality of R134a, the enthalpy of state 1 can be calculated as follows
h1 1 x h1l xh1v
kJ
h1l 17.412 (sat. R134a(l) at 0.12 MPa)
mol
kJ
h1v 39.295 (sat. R134a(v) at 0.12 MPa)
mol
Therefore,
kJ kJ
h1 1 0.2817.412 0.28 39.295
mol mol
kJ
h1 20.74
mol
kJ kJ
39.295 20.74
h h
COP 2 1 mol mol 4.97
h3 h2 kJ kJ
43.031 39.295
mol mol
90
Is this modification practical?
No. An isentropic turbine adds significant level of complexity to the cycle. Turbines are
expensive and wear over time. Furthermore, the real turbine added to the cycle will not be 100%
efficient, so the COP will not increase as much. The cost of the turbine is not justified by the
increase in COP.
91
3.42
The subscript h refers to the hotter cycle, while c refers to the cooler cycle. The
(a)
The flow rate of the cooler cycle can be found by performing an energy balance on the
condenser/evaporator shared between the two cycles. An energy balance shows:
Q H , c Q C , h Q
n c h8 h7 n h h2 h1
We can find the enthalpies for positions 2 and 8 directly from the thermodynamic tables because
the fluid is saturated in at these positions. To find the other enthalpies we must use the following
relationships:
s 6 s7
h1 h4
J
h1 24181 (sat. R134 liquid at 0.7 MPa)
mol
J
h2 40967 (sat. R134 vapor at 0.35 MPa)
mol
J J
h7 41510 (R134 vapor at 0.35 MPa with s 177.89 )
mol mol K
J
h8 21099 (sat. R134 liquid at 0.35 MPa)
mol
J J
0.5 mol/s 40967 24181
mol mol
n c
J J
21099 41510
mol mol
mol
nc 0.411
s
92
(b)
Performing an energy balance we find
Q C n c h6 h5
J
h5 h8 21099
mol
J
h6 39295 (sat. R134 vapor at 0.12 MPa)
mol
Therefore,
mol J J
Q C 0.411 39295 21099 7.48 kW
s mol mol
(c)
The power input is calculated as follows:
WC , total WC , c WC , h
where
WC , c nc h7 h6
WC , h n h h3 h2
We have all of the required enthalpies except h3. State 3 is constrained because
J
s3 s 2 175.95 (sat. R134a vapor at 0.35 MPa)
mol K
P3 0.7 MPa
J J
h3 42428 (R134a vapor at 0.7 MPa with s 175.95 )
mol mol K
93
J J J J
WC , total 0.411 mol/s 41510 39295 0.5 mol/s 42428 40967
mol mol mol mol
W C , total 1.64 kW
(d)
The coefficient of performance is calculated using the following equation:
Q C 7.48 kW
COP 4.56
WC , total 1.64 kW
(e)
The COP for the cascade is 4.56, while the COP is 4.05 in Problem 3.40. The cascade systems
COP is 12.6% greater.
94
3.43
One design follows; your design may differ:
In order to cool a system to -5 C, the temperature of the fluid must be colder than -5 C so that
heat transfer will occur. We arbitrarily specify that the working fluid evaporates at -10 C.
Similarly, in order to eject heat to the 20 C reservoir, the fluid must condense at a temperature
greater than 20 C. Arbitrarily, we choose 25 C. One possible refrigeration cycle is presented
below.
valve
T 3
1
evaporator QH
4
refrigerated 4
unit at -5 C QC
QH 20 C
reservoir
2
1 2
condenser QC
3
Wc s
compressor
A number of fluids will work sufficiently for this system, but the design process will be
illustrated using R134a. For states 3 and 4, the pressure is constant, and for states 1 and 2, the
pressure is constant at a different value. From the NIST website, we find
Now that the pressures are known, we can compute the required flow rate required in order to
provide 20 kW of cooling.
Q C n h2 h1
We can get h2 from the thermodynamic tables for saturated R134a. In order to find h1, we can
use the following relationship:
h1 h4
J
h2 40064 (sat. R134 vapor at 0.201 MPa)
mol
95
J
h4 23931 (sat. R134 liquid at 0.665 MPa)
mol
Q C 20,000 J/s
n
h2 h1 J J
40064 23931
mol mol
mol
n 1.24
s
96
3.44
The schematic of the process and the corresponding Ts is presented below
T = 5 oC
valve
Fridge 6
T
2 QC,R 1
valve 1 QH
evaporator
Wc
3 6
QC,R
5 oC
T = -15 oC Freezer evaporator QH high T 2 3
QC,F reservoir QC,F
4 -15 oC 4 5
5 condenser
Wc s
compressor
A number of refrigerants will work for this system, but the design process will be illustrated for
R134a only. To define each state, we need to thermodynamically constrain each state. You
should note that the problem doesnt state what temperature the fluid condenses at. Therefore,
we can assume the temperature is 25 C. By using information in the Ts diagram and from the
NIST website, we find
Before solving for additional pressures and temperatures, we will list the known temperatures
and pressures.
States 4, 5, and 6 are completely constrained as confirmed by Gibbs phase rule. Now, we need
to find the liquid and vapor compositions of states 1, 2, and 3 to completely constrain the states.
Since the heat duties are equal for the refrigerator and the freezer, we have the following
relationship:
Q C , F Q C , R
97
which upon performing an energy balance becomes
h4 h3 h2 h1
h1 h6
h2 h3
Substituting these relationships into the expression that equates the heat loads results in
h h6
h2 h3 4
2
J
h4 39755
mol
J
h6 23931
mol
J
h1 h6 23931
mol
J J
39755 23931
h2 h3 mol mol 31843 J
2 mol
Technically, the states are now all constrained, but we would like also like to know the vapor and
liquid compositions of each state. The compositions of states 4, 5, and 6 are already known; we
can calculate the compositions of states 1, 2, and 3 as follows:
h1 1 x1 h1sat sat
, l x1h1, v
h2 1 x2 h2sat sat
, l x2 h2, v
h3 1 x3 h3sat sat
, l x3 h3, v
98
J J
h1sat
, l 21095 mol h1sat
, v 40965 mol
J J
h2sat
, l 21095 mol h2sat
, v 40965 mol
J J
h3sat
, l 18380 mol h3sat
, v 39755
mol
99
3.45
In this case, the working material is a solid. The four states of the magnetic material are shown
of the sT diagram below. Note the axis are shifted from the usual manner.
3
2
QH
QC
1 4
(a) The heat expelled by the cold reservoir can be approximated by:
qC TC s 2 s1
where T is the average temperature between states 1 and 2, which is approximately 1 K. Taking
values of S/R from the sT diagram, we get:
J
qC R1.9 1.3 5.0
mol K
(b) Similarly, the heat absorbed by the hot reservoir can be approximated by:
J
q H TH s 4 s 4 8.8R1.2 1.9 51
mol K
qC qC
COP 0.11
wC ,total q H qC
(d) The value of COP is much lower than a typical refrigeration process (COP= 4-6); in general,
refrigeration processes at these low temperatures are much less efficient.
100
3.46
When the polymer is unstretched it is in a more entangled state. When stretched the polymer
chains tend to align. The alignment decreases the spatial configurations the polymer can have,
and therefore, reduces that component of entropy. If the process is adiabatic, the entropy of the
system cannot decrease. Consequently, its thermal entropy must increase. The only way this can
be accomplished is by increased temperature.
101
3.47
Assume the temperature is 298 K. The following data was taken from Table A.3.2:
The data listed above were used to create the next table using
h f g f
s f
T
Species s f (kJ/mol K)
Cu2O (s) -0.0766
O2 (g) 0
CuO (s) -0.0932
Now the change in entropy of the reaction can be calculated in a method analogous to Equation
2.72
s rxn
i
kJ
vi s f 2 - 0.0766
kJ
4 - 0.0932
mol K
mol K
kJ
s rxn 0.22
mol
Does this violate the second law of thermodynamics? This problem shows that the entropy
change of the system is negative, but nothing has been said about the entropy change of the
universe. We must look at the change in entropy of the surroundings to determine if the second
law is violated. By looking at the enthalpies, we see that the reaction is exothermic, which
means that heat is transferred from the system to the surroundings. Therefore, the entropy of the
surroundings will increase during this reaction.
102
3.48
The temperature and pressure terms in the equation for entropy do not contribute anything to the
entropy change because they are constant. Therefore, the only remaining entropy contribution,
the randomness of the atomic arrangements, must be considered. The randomness does not
increase when CdTe forms. In pure crystals of Cd and Te, the location of each atom is known
because the crystal lattice constrains the atomic locations. In CdTe, the crystal lattice still
defines the location of each atom, so the randomness has not increased. Therefore, the change in
entropy is zero.
3.49
This argument is not scientifically sound. Morris is arguing that since evolution results in more
order, the second law of thermodynamics is violated, so evolution must be impossible. However,
the flaw in this argument is caused by ignoring the entropy change of the entire universe. The
second law states that the entropy of the universe will remain constant or increase for any
process. Morris argument was based on the entropy of the system undergoing evolution not
the entropy of the entire universe. A system can decrease in entropy if the entropy of the
surroundings increases by at least that much.
103
3.50
Entropy is related to the number of configurations that a state can have. The greater the number
of configurations, the more probable the state is and the greater the entropy. We can
qualitatively relate this concept to the possible hands in a game of poker, but it is more
interesting to quantify the results using some basic concepts of probability.
We consider a hand of poker containing 5 cards randomly draw from a deck of 52 cards.. There
are a finite number of permutations in which we can arrange a 52 card deck in 5 cards. For the
first card in the hand, we select from 52 cards, the second card can be any other card so we select
from 51 cards, the third card has 50 cards, and so on. Thus the number of permutations of 5
cards is:
P 52 51 50 49 48 311,875,200
However, we do not care the order in which the cards are dealt, so we must divide this number
by the number of ways we can come up with the same hand. We do this math in a similar way.
For a given hand there are five cards we can pick first, four we can pick second, . So the
number of ways we can make the same hand is:
N 5 4 3 2 1 120
The number of unique configurations can be found by dividing P by N . Thus, the cards can
display
P
C 2,598,960
N
or 2,598,960 unique configurations. To find the entropy of a given hand, we need to find out
how many of these unique configurations belong to the hand
Consider a four of a kind. There are 13 different possible ranks of four of a kind, one for each
number A, 2, 3, 4, 5, 6, 7,8, 9, 10, J, Q, K. The fifth card in the hand could be any of the other
48 cards. Therefore, the number of combinations of four of a kind is: 13 x 48 = 624
The probability of a four of a kind is: 624 / 2598960 = 1 / 4165 which is very unlikely. Thus the
entropy of this hand is low. In contrast, there are 1,098,240 to have a hand that has one pair;
therefore the probability of getting this hand is much greater, 1/2.4, and its entropy is high.
The probability of having nothing is 1/2 which represents the most likely hand in poker, or the
hand with the highest entropy. In fact, the rules of poker are defined so the hand of lower
entropy always beats the hand of higher entropy.
104
Chapter 4 Solutions
Engineering and Chemical Thermodynamics
Wyatt Tenhaeff
Milo Koretsky
1
4.1.
(a)
Yes, the form of the equation is reasonable. This can be rewritten:
RT a
PA n
v A v A
It is equivalent to Pv RT , except the pressure term has been corrected to account for the ions
intermolecular forces. The coulombic forces between the gas molecules affect the system
pressure. This modification is similar to the van der Waals equation. Since we are limited to 1
parameter, we need to choose the most important interaction. Since net electric point charges
exert very strong forces, this effect will be more important than size.
(b)
The sign should be negative for a because the positively charged gas molecules repel each other
due to coulombic forces. Therefore, the system pressure increases, i.e.,:
PA Pideal
Coulombic forces are much stronger than van der Waals interactions, so a will be large much
larger than a values for the van der Waals EOS.
(c)
Coulombic repulsion is the primary intermolecular force present in the gas. Coulombic potential
energy is proportional to r-1. v is proportional to r3, so the coulombic potential goes as v-1/3.
Therefore,
1
n
3
a must have must have the following units to maintain dimensional homogeneity with the
pressure term:
a N J
v1/ 3 m 2 m 3
so
a J
1/3
m mol
2
2
4.2
The attractive interactions are described by van der Waals forces. For both O2 and propane, the
dipole moments are zero. Therefore, the expression for the interactions reduce to
3 i j Ii I j
2 r6 Ii I j
aa
3 16 10 25 cm 3 1.933 10
2 11
erg 1.933 10 11 erg
11
2 r6 1.933 10 erg 1.933 10 11 erg
aa
3.71 10 59
r 6
cm 6
erg
ab
3 16 10 25 cm 3 62.9 10 25 cm 3 1.933 10 11
erg 1.753 10 11 erg
2 r6 1.933 10 11 erg 1.753 10 11 erg
ab
1.39 10 58
r 6
cm 6
erg
bb
3 62.9 10 25 cm 3 1.753 10
2 11
erg 1.753 10 11 erg
2 r6 1.753 10 11 erg 1.753 10 11 erg
bb
5.20 10 58
r 6
cm 6
erg
(b)
Calculation:
aa bb
3.71 10 59
r6
cm erg
6
5.20 10 58
r6
cm 6 erg
3
aa bb
r
1.39 10 58
6
cm
erg 6
Note: Disregarded the positive value.
The value of aa bb is equal to ab . The values are equal because the ionization energies are
similar.
(c)
An expression for the average intermolecular attraction in the mixture can be found using the
mixing rules
mix y a2 aa 2 y a yb ab yb2 bb
mix
3.71 10 59 y a2 1.39 10 58 y a yb 5.20 10 58 yb2
cm 6 erg
6
r
4
4.3
(a)
300 K, 10 atm. The intermolecular distance of molecules is greater at lower pressures.
Therefore, fewer intermolecular interactions exist, which cause less deviation from ideality.
(b)
1000 K, 20 atm. At higher temperatures, the kinetic energy of the molecules (speed) is greater.
The molecules interact less; thus, the compressibility factor is closer to unity.
(c)
Let subscript 1 denote BClH2 and 2 denote H2. For the mixture, we calculate the
compressibility factor as follows
B
z 1 mix
v
where
Bmix y12 B1 2 y1 y2 B12 y 22 B2
B1 Bmix B2
0.9
zBClH2 lowest
0.8
z
0.7
0.6
0 0.2 0.4 0.6 0.8 1
yBClH2
5
4.4
(a)
The intermolecular attractions and volume occupied by the styrene monomers will contribute to
the deviations from ideality. Since styrene monomers are essentially non-polar, the order of
importance is as follows
The van der Waals EOS is appropriate since it accounts for the occupied molar volume and
intermolecular forces, but it should be noted that more modern EOSs will give more accurate
results.
(b)
Use critical data to calculate the a and b parameters:
2
J
8.314 647.15 K
27 RTc mol K 3.13 J m
2 3
27
a
64 Pc 64 39 105 Pa 2
mol
J
8.314 647.15 K
RTc mol K m3
1.72 10 4
b
8 Pc
8 39 105 Pa
mol
Now we can use the van der Waals EOS to solve for temperature.
RT a
P
v b v2
J J m3
8.314 T 3.13 2
mol K mol
10 10 Pa
5
4 m m3 2
m3
3
3.0 10 1.72 10 4 3.0 10 4
mol
mol
mol
T 550.8 K 277.65 C
(c)
The a parameter is related to attractive intermolecular forces. Dispersion is the controlling
intermolecular force in this system, and its magnitude is directly related to the size of the
molecules (polarizability component of dispersion). For the 5-monomer long polymer chain, a is
5 times the a value in Part (b). The b parameter is also related to the size of the molecule since it
accounts for the volume occupied by the molecules. Again, the b parameter for the reduced
polymer chain is 5 times the b value from Part (b).
6
J m3
a 15.65 2
mol
m3
b 8.6 10 4
mol
(d)
We must realize that if we initially believed there were 100 moles of styrene in the reactor, then
there can only be 20 moles of the 5-monomer long polymer chain. Therefore,
m3
v 1.5 10 3
mol
and
J J m3
8.314 T 15 . 65 mol 2
mol K
10 10 Pa
5
3
3
2
3 m
m m 3
1.5 10 3 4
mol 1.5 10
8 . 6 10
mol mol
T 612.41K 339.26 C
7
4.5
There are many ways to solve this problem, and the level of complexity varies for each method.
To illustrate the principles in Chapter 4, two of the simplest solution methods will be illustrated.
The polarizability of a molecule scales with the number of atoms; the polarizabilities of
individual atoms are additive. Using the first two molecules, solution of the following system of
equations
1 C 4 H 1 CH 4
2 C 6 H 1 C2 H 6
gives
C 11.4 10 25 cm 3
H 3.65 10 25 cm
3
Now, the polarizability of the chlorine atom can be found. For chloroform,
1 C 4 Cl 1 CCl 4
Therefore,
Cl 23.4 10 25 cm 3
The values of C3 H 8 , CH 3Cl , CH 2Cl2 , and CHCl3 are calculated with the polarizabilities found
above as follows, and compared to the values given.
All agree reasonably well. Now, the polarizabilities of C 4 H10 and C 2 H 5Cl will be calculated.
C 4 H 10
4 11.4 10 25 cm 3 10 3.65 10 25 cm 3 82.1 10 25 cm 3
8
C 2 H 5 Cl
2 11.4 10 25 cm 3 5 3.65 10 25 cm 3 1 23.4 10 25 cm 3 64.45 10 25 cm 3
For this method, we calculate the molecules polarizabilities by adding the polarizability of each
bond, instead of the atoms. For the methane molecule
4 C H 1 CH 4
C H
26 10 25 cm 3
6.5 10 25 cm 3
4
6 C H C C 1 C 2 H 6
C C 44.7 10 25 cm 3 6 6.5 10 25 cm 3 5.7 10 25 cm 3
The C-C and C-H polarizability calculated above predict the given polarizability of propane well.
The polarizability of C-Cl bonds is calculable with the polarizability of chloroform.
4 C Cl 1 CCl 4
C Cl
105 10 25 cm 3
26.25 10 25 cm 3
4
The values of C3 H 8 , CH 3Cl , CH 2Cl2 , and CHCl3 are calculated with the polarizabilities found
above as follows, and compared to the values given.
All agree reasonably well. This value predicts the polarizabilities of the other species in the table
reasonably well. Now, the polarizabilities of C 4 H10 and C 2 H 5Cl will be calculated.
1 C 4 H 10 3 C C 10 C H
C 4 H 10
3 5.7 10 25 cm 3 10 6.5 10 25 cm 3 82.1 10 25 cm 3
9
1 C 2 H 5 Cl 1 C C 5 C H 1 C Cl
C 2 H 5 Cl
1 5.7 10 25 cm 3 5 6.5 10 25 cm 3 1 26.25 10 25 cm 3 64.45 10 25 cm 3
Note both the atom method and the bond method yield identical results. More accurate values
for the polarizabilities can be calculated using more of the data given in the problem.
10
4.6
(a)
increases with molecular size. Therefore,
I S O
2 2 2
scales with magnitude of van der Waals interactions. Because these are non-polar, diatomic
molecules, only dispersion forces are present. Dispersion forces depend on the first ionization
potential and polarizability. Ionization energy is approximately equal for each molecule. The
polarizability scales with molecular size, so
I S O
2 2 2
(b)
increases with molecular size. Diethylether and n-butanol have the same atomic formula and
similar spatial conformations. Therefore, they should be about equal in size. Methyl ethyl
ketone has fewer atoms, but has two exposed electron pairs on the double-bonded oxygen. The
size of the molecular electron orbital of methyl ethyl ketone is approximately equal to the sizes
of diethyl ether and n-butanol, so
Methyl ethyl ketone and n-butanol are much more polar than diethyl ether due to their greater
asymmetry, so their values are greater than diethyl ethers. Now we must determine if there is
greater attraction in n-butanol or methyl ethyl ketone. scales with the magnitude of van der
Waals interactions. Since induction and dispersion forces are similar in these molecules, we
must consider the strength of dipole-dipole forces. There is greater charge separation in the
double bond of ketone, so
11
4.7
(a)
At 30 bar, the water molecules are in closer proximity than they are at 20 bar. Intermolecular
attractions are greater, so the magnitude of molecular potential energy is greater. The potential
energy has a negative value for attractive interactions. The molecular kinetic energy is identical
since the temperature is the same. Hence, the internal energy, the sum of kinetic and potential
energies, is less at 30 bar.
(b)
The key to this phenomenon is hydrogen bonding. At 300 K and 30 bar, isopropanol and n-
pentane are both liquids. The hydrogen bonding and dipole-dipole interactions are present in
isopropanol, and dispersion is present in n-pentane. The intermolecular forces are greater in the
isopropanol, so the compressibility factor is smaller for isopropanol.
At 500 K and 30 bar, both species are gases. In the gas phase, hydrogen bonding does not play a
significant role. The dispersion forces in n-pentane are stronger than the dipole-dipole forces of
isopropanol. Therefore, the compressibility factor is smaller for n-pentane.
12
4.8
(a)
Ideal: For ideal NH3, the compressibility factor is equal to one. For real NH3, the strong
intermolecular forces (dipole-dipole and dispersion) cause the molar volume to decrease. They
outweigh the volume displaced by the physical size of NH3; thus, z will be less than one.
(b)
Internal energy will be greater for the ideal gas. In the real gas, intermolecular attractions are
present. Internal energy value is the sum of potential and kinetic energies of the molecules. The
absolute values of the kinetic energies are identical at identical temperature: however, the
potential energy decreases for real NH3 due to attractive interactions - so the internal energy is
less for the real NH3.
(c)
The entropy will be greater for the ideal gas. Entropy is a measure of possible molecular
configurations or randomness. Ammonia has an electric dipole in which positive and negative
charge are separated. The intermolecular forces in the real gas cause the molecules to align so
that the positive charge in one molecule is adjacent to a negative charge in a neighboring
molecule to reduce potential energy. Therefore, fewer possible configurations exist, which
creates less randomness and lower entropy.
13
4.9
(a)
We can determine which case has the higher compressibility factor by comparing the molar
volumes at constant T and P. With Ne, very weak intermolecular attractions are present, so
volume displacement becomes important. The compressibility factor will be slightly greater than
one. In NH3, the strong intermolecular attractions decrease the molar volume, so z is less than
one. The compressibility factor is greater for Ne.
(b)
Entropy is a measure of randomness. Both species are gases at these conditions. The
intermolecular attractions present in NH3 reduce the number of possible configurations. The
weak forces present in Ne have a much smaller effect. However, NH3 is asymmetrical, while Ne
is symmetrical. The asymmetry of NH3 results in more possible configurations that NH3 can
have. Therefore, it is difficult to qualitatively show for which case the entropy is greater. Since
both species are gases, intermolecular interactions are relatively weak, and we can guess that
entropy is greater in NH3.
14
4.10
(a)
To find the average distance between the two atoms of Ar, we can find the volume that each
atom occupies. The molar volume can be found from the compressibility factor. At 300 K and
25 bar,
Tr 1.99
Pr 0.513
z 0.9859
Therefore,
J
8.314 300 K
mol K 4
m3
v 0.9859 9.84 10
25 105 Pa mol
m3
v 1.63 10 27
atom
The second number was found by dividing the molar volume by Avogadro's number. To
estimate the distance between each atom, we note that the distance between molecules can be
related to the volume by:
r3 v
3
v
r
Ar Ar
A rough estimate of r is
r 1.63 10 27 m 3 1 / 3 1.18 109 m
15
r 11.8 A
(b)
The potential energy due to gravity can be calculated as follows
G m 2Ar
G
r
m2
where G is the gravitation constant G 6.67 10 11 and mAr is the mass of an argon
2
kg s
atom. The mass of an argon atom is calculated as follows
39.948 g Ar 1 mol 1 kg
kg
m Ar 6.633 10 - 26
1 mol Ar 6.023 10 atoms 1000 g
23 atom
G 2.49 10 -52 J
(c)
Equation 4.13 quantifies the potential energy due to London interactions
3 Ar Ar I Ar I Ar
Ar Ar
2 r6 I Ar I Ar
Using the distance calculated in Part (a) r 1.18 10 7 cm , we get
(d)
The potential energy due to London interactions is around 10 35 times greater than the potential
energy due to gravity. Clearly, London interactions are much more important, and gravitational
effects can be neglected.
16
4.11
We want to compare the Lennard-Jones potential to one with an exponential repulsion term. As
provided in the text, Equation 4.19, the Lennard-Jones potential is
12 6
4
r r
1 1
*
4 r *12 r *6
where
r
r*
c2 1
*
exp c1 exp
r * r * 12
We have two adjustable parameters, c1 and c2, to match the first (LJ) potential to the second
(exp) potential. We need to choose reasonable criteria to specify. For this solution we choose
equal well depths and equal values at =1. Other choices may be just as valid; you should
realize that you have two parameters to fit and so must specify two features.
This solution is shown at two magnifications in the plots on the following page:
17
Potential functions
0.3
0.2
Lennard-Jones
exponential
0.1
0
0.9 1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2 2.1
-0.1
-0.2
-0.3
Potential functions
20
15 Lennard-Jones
exponential
10
0
0.7 1 1.3 1.6 1.9
-5
1. The most stable configuration (the bottom of the well) occurs at a greater separation for
the exp model.
2. The Lennard-Jones potential increases more steeply at small radii, i.e., it behaves more
like the hard-sphere potential.
3. The two models are in reasonable qualitative agreement
18
4.12
(a)
The bond strength of a sodium ion can be viewed as the amount of energy it would take to
remove the sodium ion from the crystal lattice. The interaction between the chlorine and sodium
ions is Coulombic attraction.
6
QNa QCl
6
4.803 10 10 esu 4.803 10 10 esu
5.01 10 11 erg
r
2.76 10 -8 cm
31.3 eV
(b)
QNa QCl Q Q
6 12 Na Na
r r
5.01 10 11 erg 12
4.803 10 10
esu 4.803 10 10 esu
3.90 10 -8 cm
2.089 10 11 erg 13.0 eV
(c)
QNa QCl Q Q Q Q
6 12 Na Na 8 Na Cl
r r r
2.089 10 11 erg 8
4.803 10 10 esu 4.803 10 10 esu
4.78 10 -8 cm
1.77 10 11 erg 11.1 eV
(d)
QNa QCl Q Q Q Q Q Q
6 12 Na Na 8 Na Cl 6 Na Na
r r r r
1.77 10 erg 6
11
4.803 10 esu 4.803 10 10 esu
10
5.52 10 -8 cm
7.37 10 12 erg 4.60 eV
19
4.13
The boiling points of the halides depend upon the strength of intermolecular attractions. The
stronger the intermolecular attraction, the higher the boiling point. Dispersion and dipole-dipole
interactions are present in all five species listed. The magnitude of the dipole-dipole interactions
is similar so the pertinent intermolecular force in these molecules is dispersion. The molecular
size increases from left to right. Polarizabilities are greater in larger molecules, which manifests
in larger dispersion forces. Therefore, the boiling point increases from left to right.
20
4.14
van der Waals forces hold the Xe atoms together in a molecule of Xe2. The potential energy can
be quantified with the Lennard-Jones potential function. The bond length is the r value where
the potential is a minimum. From Table 4.2:
229 K
k
and
4 .1 A
12 6
4
r r
d 12 12 6 6
4 0
dr r r r r
where we set this derivative equal to zero to find the minimum. Solving gives
6
1
r 2
or
r 6 2 1.12 4.60 A
21
4.15
The data in the following table were taken from Table A.1.1
Species Tc K Pc 10 -5 Pa
He 5.19 2.27
CH4 190.6 46.00
NH3 405.6 112.77
H2O 647.3 220.48
The van der Waals a parameter can be calculated using Equation 4.39.
27 RTc 2
a
64 Pc
The van der Waals b parameter can be calculated using Equation 4.40.
RTc
b
8 Pc
Using the data table and equations listed above, the following table was created; the values of
dipole moment and polarizability reported in Table 4.1 are also included.
J m3 m3
Species a b 105
2
mol mol [D] [cm x 1025]
3
The values of a for helium is two orders of magnitude less than the other species since it only has
weak dispersion forces (small atom, small The values of a for methane, ammonia, and water
are of the same magnitude because the sums of the intermolecular attractions for each molecule
are similar. All three molecules have comparable dispersion forces; although slightly weaker in
the ammonia and water. However, unlike methane, these two molecules also have dipole-dipole
and induction forces. In fact, the strong dipole in water gives it the largest value
The values of b are all of the same magnitude, as expected since b scales with size according to
the number of atoms in the molecule.
Size of Molecules: CH 4 NH 3 H 2O He
Value of b: bCH 4 bNH 3 bH 2O bHe
22
4.16
The van der Waals b parameter can be calculated using Equation 4.40.
RTc
b
8 Pc
Critical point data can be found in Table A.1.1. The following table was made:
Species Tc K Pc 10 -5 Pa m3
b 10 5
mol
CH4 190.6 46.00 4.306
C6H6 562.1 48.94 11.94
CH3OH 512.6 80.96 6.58
J
Note: Used R 8.314 for calculation of a.
mol K
3b
3
2N A
The table listed below was created using this equation, and data from Table 4.2 are included
along with the percent difference.
23
4.17
The van der Waals a parameter can be calculated using Equation 4.39. The values for the above
equation were taken from Table A.1.1, and the following table was made:
J m3
Species Tc K Pc 10 -5 Pa a
2
mol
CH4 190.6 46.00 0.2303
C6H6 562.1 48.94 1.883
CH3OH 512.6 80.96 0.9464
J
Note: Used R 8.314 for calculation of a.
mol K
2N a2C6
a
3 3
3a 3
C6
2N a2
The values can be found in Table 4.2. The following table can now be created
Species 1010 m
C 6 10 77 J m 6
CH4 3.8 1.66
C6H6 5.27 36.3
CH3OH 3.6 5.82
To compare the values obtained from Equation 4.13, first calculate the potential energy with
Equation 4.13:
3 i j Ii I j
2 r6 Ii I j
Therefore
3 II
C6 i j i j
2 Ii I j
We obtain the following values using this equation, and the corresponding percent differences
were calculated.
24
Species
C 6 10 77 J m 6
Percent Difference
Page 187 Equation 4.13
CH4 1.66 1.021 63
C6H6 36.3 12.0 202
CH3OH 5.82 1.36 328
The values for the van der Waals a constant have the correct qualitative trends and order of
magnitude; however, those values predicted from basic potential theory vary significantly from
corresponding states. The basic potential result presented in the text assumes that the species are
evenly distributed throughout the volume. It does not take into account the structure given to the
fluid through intermolecular forces. In fact, a more careful development includes a radial
distribution function, which describes how the molecular density of the fluid varies with r. The
radial distribution function depends on pressure and temperature of the fluid.
25
4.18
RT a
van der Waals: P
v b v2
RT a
Redlich-Kwong P 1/2
v b T v(v b)
RT a(T )
Peng-Robinson: P
v b v(v b) b(v b)
As you may have discovered, these equations are largely empirical with no theoretical
justification. They simply represent experimental data better. We can use our knowledge of
intermolecular forces, however, to explain why these may work better.
If we compare these equations to the van der Waals equation, we note that the first term
on the right hand side of all three equation is identical; we accounted for this term as a correction
for finite molecular size (or alternatively repulsive interaction due to the Pauli Exclusion
Principle). This form represents a hard sphere model.
The second term, that which deals with intermolecular attraction, is different in all three
models. Both of the later equations include a temperature dependence in this term. We have
seen that if attractive forces depend on orientation (dipole-dipole), they fall off with T as a result
of the averaging process (recall discussion of Equation 4.11).
Another explanation goes as follows: as T increases, the molecules move faster, reducing
the effect of intermolecular forces. If we say that the potential energy between two molecules
depends on the amount of time that they spend close to each other, then it would be inversely
related to velocity (The faster molecules are moving, the less time they spend in the vicinity of
other species). In this case, the correction term would go as V-1, where V is the molecular
velocity. If we relate molecular velocity to temperature
2
1
2 mV 32 kT
Therefore the correction term goes as T-1/2, as shown in the Redlich-Kwong equation.
The inclusion of a "b" term in the second term may help relax van der Waal's "hard
sphere" model with a more realistic potential function, i.e., something closer to a Lennard Jones
potential (Figure 4.8) than the Sutherland potential (Figure 4.7) upon which the van der Waals
equation is based. It makes sense that this should be included in the force correction since this is
taking into account repulsive forces. One example of a more detailed explanation follows:
If we look at the Redlich-Kwong equation, it says that if we have 2 species with the same
attractive strength (same a -> same magnitude of van der Waals forces), the larger species will
have less of an effect on P. The following sketch illustrates how 2 species could have the same
van der Waals attractive forces:
26
London
- + + = London
Dipole
Species 1 Species 2
smaller and polar larger and non-polar
(smaller b, same a) (larger b, same a)
In the case above, when the two species are the same distance apart, they have the same
attractive force; however, the smaller species (1) can get closer before its electron cloud
overlaps. Thus it has more opportunity for attractive interactions than the larger species. The
Peng-Robinson equation exhibits the most complicated form in an attempt to better fit
experimental data.
27
4.19
(a)
Work is defined as follows
W PdV
nRT
W dV
V
J 1 L
W 2.0 mol 8.314 1000 K ln
mol K 10 L
W 38.29 kJ
(b)
Instead of substituting the ideal gas law into the definition of work, the Redlich-Kwong equation
is used:
v2
RT a
W n dv
v b T v (v b)
1 / 2
v1
Substituting
J
R 8.314
mol K
T 1000 K
n 2 mol
J K1/2 m 3
a 14.24
2
mol
3
5 m
b 2.1110
mol
m3
v1 0.0005
mol
m3
v2 0.005
mol
28
and evaluating the resulting formula gives
W 37.35 kJ
(c)
Energy balance:
u q w
q Ts
and
w u Ts
To use the steam tables conveniently, we need the initial and final pressures. We can calculate
these with the Redlich-Kwong EOS:
P1 1.65 MPa
P2 15.6 MPa
kJ kJ
u1 3522.6 s1 8.101
kg kg K
kJ kJ
u 2 3462.2 s2 7.00
kg kg K
Therefore,
kJ kJ kJ kJ
w 3462.2 3522.6 1000 K 7.00 8.101
kg kg kg K kg K
kJ
w 1040
kg
W 37.5 kJ
The answers from the three parts agree very well. Part (a) is not as accurate as Part (b) and Part
(c) because water is not an ideal vapor. The value from Part (b) is 1.1 % smaller than the value
from Part (c). Clearly, the Redlich-Kwong EOS or steam tables are appropriate for this
calculation.
29
4.20
First, we can obtain an expression for the compressibility factor with the van der Waals equation.
Pv v a 1 a
RT v b RTv b RTv
1
v
1
1 x x 2 x 3 ...
1 x
Therefore,
2
1 b b
1 ...
b v v
1
v
and
a
b 2
Pv RT b
1 ...
RT v v
This expression can also be expanded in pressure. From Equation 4.58, we know that
B
B'
RT
C B2
C'
RT 2
where
Pv
1 B' P C ' P 2 ....
RT
bRT a 2abRT a 2
B' and C '
RT 2 RT 4
Therefore,
30
bRT a
P 2abRT a P 2 ....
Pv 2
1
RT RT 2 RT 4
31
4.21
Using the virial expansion, the pressure can be written as follows
1 B C
P RT ...
v v 2 v3
P
RT
v
1 B' P C ' P 2 ...
1 B C
2
1 B C RT 1 B C
RT ... 1 B ' RT C ' RT
v v 2 v3 v v v 2 v 3 v v 2
v3
If we combine like terms on the right side and set them equal to terms on the left, we find
B B' RT
C BB' RT C ' RT 2
Substitute the expression for B into the expression for C and solve for C:
CB
C'
RT 2
and
B
B'
RT
32
4.22
At the critical point we have:
RTc a
Pc (1)
vc b Tc vc2
P RTc 2a
0 (2)
v Tc vc b Tc vc3
2
2
P 2 RTc 6a
2 0
v
Tc
vc b 3
Tc vc4 (3)
4a 6avc b
0
vc3 vc4
vc 3b
9
a vc RTc2
8
Finally, if we plug this back into Equation 1 we can solve for the Berthelot constants in terms of
the critical temperature and the critical pressure:
27 R 2Tc2 RTc
a and b
64 Pc 8Pc
33
4.23
Before we can find the reduced form, we need to find expressions for a and b in terms of critical
point data (Problem 4.22). We find
vc 3b (1)
27 R 2Tc3
a (2)
64 Pc
RTc
b (3)
8 Pc
RT a
P
v b Tv 2
First, substitute Equation 1 into the first term on the right-hand side and rearrange to get
RT
a
P b
3vr 1 Tv 2
8 PcTr a
P
3v r 1 Tv 2
Substitute Equation 2:
27 R 2Tc3
8 PcTr 64 Pc
P
3v r 1 Tv 2
From Equation 3:
R 2Tc2 64 Pc2 b 2
8PcTr 27 PcTc b 2
P
3v r 1 Tv 2
34
Finally, substitute Equation 1:
8Tr 3
Pr
3vr 1 Tr vr2
35
4.24
(a)
The virial equation is:
Pv B C D
z 1 2 3 ... (1)
RT v v v
Pv C
z 1v 1v B ... (2)
RT v
PvT data from the steam tables are given in the steam tables and corresponding values of (z-1)v
and 1/v can be calculated. An illustrative plot of (z-1)v vs. 1/v for T = 300 oC is shown below.
-120
(z -1)v (cm3 / mol)
-125
(z -1)v
-130
-135
0
2
0.5
1.5
1/v (mol/l )
We see that at pressures above 1.5 MPa (15 bar) we see that the plot reaches a constant value of
around -118 cm3/mol. We may choose to report this value as B. A more careful examination of
Equation 2 suggests another possibility.
If we plot (z-1)v vs. 1/v, we should get a straight line at low to moderate pressures with an
intercept equal to the second virial coefficient, B, and a slope equal to the third virial coefficient,
C. At very low pressures, the curve is indeed linear and gives an intercept of -140 cm3/mol, as
shown on the next page. The slope of this region would yield the third virial coefficient, C.
36
The first value uses much more data while the second method uses limited data in a better range.
What value would you be more apt to use? Water exhibits this odd general behavior at all
temperatures.
-120
(z -1)v (cm3 / mol)
y = 177.058x - 140.373
-125
(z -1)v
-130 lin
-135
0
2
0.5
1.5
1/v (mol/l )
If we do this at other values of T, we can compile a set of second virial coefficients vs.
temperature and get an idea of the differences in the two approaches. Temperatures of 200, 250,
300, 350, 400, 500, 600, 700, 800, 900, 1000, and 1200 oC analyzed in this manner. The results
are reported in the table and figure below. Also reported were the average value and the value at
the lowest pressure used. Note that the average value is indicative of the 1st method above while
the value at the lowest pressure used is indicative of the second value used.
37
Second Virial Coefficient of Water by 4 Methods
-50
B (water) cm 3 /mol
-100
B (average)
B (lowest P)
-150
B (Level)
-200 B (Linear)
CRC
-250
1000
1250
1500
250
500
750
T (K)
For comparison, values of -142.2 cm3/mol and -7160 cm6/mol2 are reported for B and C,
respectively, at 260 oC1. Values of B reported in the CRC are also shown on the summary plot.
They agree most closely with the first (level) method.
(b)
There are several alternatives how to solve this problem, each of which comes up with a slightly
different result:
Alternative 1:
Rewrite the virial equation:
1 BH 2 O
P RT ...
v v2
P 1 2 BH 2 O
0 RTc 2
v Tc v v 3c
c
so
v cm3
BH 2O c 28
2 mol
1
J.H. Dymond and E.B. Smith, The Virial Coefficients of Pure Gases and Mixtures, A Critical Compilation, Oxford
University Press, Oxford, 1980.
38
Alternative 2:
From the virial equation:
Pv cm3
BH 2O c c 1vc 43.2
RTc mol
Alternative 3:
1 BH 2 O C H 2 O
P RT ...
v v2 v3
so
P 1 2 BH 2O 3C H 2O
0 RTc 2 (1)
v Tc v v 3
v 4c
c c
and
2P 6B 12C H 2O
0 RTc 2 H 2O (2)
v 2 v3 4
vc
5
Tc c
vc
2 2 BH 2 O
0
v 3c v 4c
so
cm3
BH 2O vc 56
mol
39
4.25
A trial-and-error is the easiest method for solving this problem. The general method is as
follows
1. Guess P sat
2. Calculate values of molar volumes that result at the chosen P sat : vlow , vmid , vhigh
3. Find areas under the curves with the following expressions
P
v mid
sat
f v dv
v low
f v P dv
v high
sat
v mid
[ f v represents the particular EOS implicit in molar volume being used.]
4. If the values of the expressions are not equal, repeat the process until they are.
(a)
1. Guess P sat :
P sat 6 bar
J K1/2 m 3
J
363.15 K
8.314 41 .851 2
6 105 Pa
mol K mol
m3 3
v 0.0001 363.15 K 1 / 2 v v 0.0001 m
mol mol
m3 m3 m3
vlow 0.000152 mid
v 0.000551 high
v 0.00459
mol mol mol
3.
0.000551
sat RT a dv 1157.19
P
v b T 1 / 2 v (v b )
0.000152
0.00459
RT a P sat dv 1325.78
v b T 1 / 2 v (v b )
0.000551
40
4. Repeat this process until the areas are equal. This occurs approximately at
P sat 6.38 bar
(b)
1. Guess Psat :
Psat 6 bar
RT a T
P sat
v b v(v b) bv b
using
J m3 m3
a 2.066
2
, b 9 10 6 , T 1.188
mol mol
m3 m3 m3
vlow 0.000152 mid
v 0.000551 high
v 0.00459
mol mol mol
3.
a T
0.000551
sat RT
P v b v(v b) bv b dv 1568.48
0.000152
a T
0.00459
RT sat
v b v(v b) bv b P dv 1100.19
0.000551
4. Repeat this process until the areas are equal. This occurs approximately at
P sat 4.945 bar
The value calculated with the Redlich-Kwong EOS is 11.9% greater than the measured value of
5.7 bar. The Peng-Robinson EOS results in a value that is 13.25% less than the measured value.
41
4.26
First, calculate a , b , and T . From Table A.1.1:
2
J
0.45724 8.314 305.4 K
mol K J m3
a 0.6036
48.84 105 Pa 2
mol
J
0.07780 8.314 305.4 K
mol K m3
b 4.045 10 5
48.84 10 Pa
5
mol
2
243.15 K
T 1 0.37464 1.542260.099 0.269920.0992 1
305.4 K
T 1.12
RT a T
P sat
v b v(v b) bv b
we get
m3 m3 3
3 m
v 7.199 10 5 , v 2.167 10
4
, v 1.578 10
mol mol mol
The molar volume of saturated ethane liquid is the smallest value from the list above, while the
molar volume of saturated ethane vapor is the largest value. Therefore,
42
g
MW ethane 36.0694
mol 1 m
3
0.501 g
liq 3
v liq
5 m
3 100 cm
3 3
cm
7.199 10
mol
g
MW ethane 36.0694
mol 1 m
3
0.0229 g
vap 3
v vap 3 m
3 100 3 cm 3 cm
1.578 10
mol
Both of the densities calculated with the Peng-Robinson EOS are larger than the reported values.
The liquid density is 7.05% larger, and the vapor density is 18.7% greater.
43
4.27
(a)
We will use the Redlich-Kwong EOS in order to obtain an accurate estimate. First, calculate the
a and b parameters:
0.42748R 2Tc2.5 J m 3 K1 / 2
a 1.56
Pc mol 2
0.08664 RTc m3
b 2.69 10 5
Pc mol
(Note: Critical data for nitrogen obtained in Table A.1.1.)
RT a
P
v b T vv b
1 / 2
Assume the temperature of the gas is 22 C, substitute values, and find the molar volume with a
solver function:
m3
v 1.97 10 4
mol
1290 m 3
n 6.55 10 6 mol
m 3
1.97 10 4
mol
m 6.55 10 6 mol 0.02801 kg/mol 1.83 105 kg
Value $6.1 / kg 1.83 105 kg $1,116,300
44
If we use the ideal gas law, we find:
J
8.314 295 K
RT mol K m3
v 1.98 10 4
P 12400000 Pa mol
Value $1,113,200
The value calculated using the ideal gas law is $3100 less than the value calculated using the
Redlich-Kwong EOS.
(b)
First, calculate the a and b parameters:
0.42748R 2Tc2.5 J m 3 K1 / 2
a 1.74
Pc mol 2
0.08664 RTc m3
b 2.21 10 5
Pc mol
(Note: Critical data for oxygen obtained in Table A.1.1.)
Assume the temperature of the gas is 22 C, substitute values, and find the molar volume with a
solver function:
m3
v 1.53 10 4
mol
Value $2,428,000
m3
v 1.64 10 4
mol
which provides
Value $2,265,000
45
The value found with the ideal gas law is $163,000 less than the value found using the Redlich-
Kwong EOS.
46
4.28
First, rewrite the Redlich-Kwong equation in cubic form.
RT 2 a RT ab
v3 v b b 2 v 0
P PT 1 / 2 P PT 0.5
RTc 2 a RT ab
v3 v c b b 2 v 0
Pc 0 .5 P 0.5
c c
P T c P T
c c
Now expand v vc 3 0
v 3 3v 2 vc 3vvc2 vc3 0
RTc
3vc (1)
Pc
a RTc
3vc2 b b2 (2)
PcTc0.5 Pc
ab
vc3 (3)
Pc Tc0.5
RT
a 3vc2 C b 2 Pc Tc0.5
PC
b 21 / 3 1 vc (4)
47
0.08664 RTc
b (Equation 4.48)
Pc
Now, substitute Equation 1 and Equation 4.48 into the expression for a to get
0.42748RTc2.5
a (Equation 4.47)
Pc
Pv
To verify Equation 4.50, substitute Equation 1 into z c c c
RTc
RT
Pc c
zc 1
3Pc
(Equation 4.50)
RTc 3
RT a
P
v b T 0.5vv b
0.2599 RT
b a
P
vr 0.2599 T 0.5vvc vr 0.2599
3Tr Pc a
P
v r 0.2599 T vvc v r 0.2599
0. 5
Pc2 b 2
R 2Tc2
0.086642
48
which upon substitution yields
3Tr Pc 0.42748 Pc b 2
P
v r 0.2599 0.08664 2 Tr0.5 vvc v r 0.2599
Therefore,
3Tr 1
Pr (Equation 4.49)
vr 0.2599 0.2599Tr vr vr 0.2599
0.5
49
4.29
Since we are concerned with liquid water, we can base all our calculations on saturated liquid
water since the molar volume of liquids are weakly dependent on pressure. Therefore, our
results are also applicable to sub-cooled water (the pressure is greater than the saturation
pressure). To find the thermal expansion coefficient, we will use the following approximation
m3 m3
v sat 15 C 0.001001 v sat 95 C 0.001040
kg kg
m3 m3
v sat 20 C 0.001002 v sat 100 C 0.001044
kg kg
m3 m3
v sat 25 C 0.001003 v sat 105 C 0.001047
kg kg
m3 m3
1 0 . 001003 0 . 001001
m3
kg
kg
20 C 0.001002
kg 10 K
m3 m3
1 0.001047 0.001040
m3
kg
kg
100 C 0.001044
kg 10 K
20 C 2.0 10 4 K -1
100 C 6.71 10 K
4 -1
The accuracy of these values may be limited since they are based on the small differences
between liquid volumes.
50
dv sat
1
v sat T dP
T
0.001050
v = -5E-07*P + 0.001044
R2 = 0.9949
Specific Volume [m /kg]
0.001040
3
0.001030
T = 20C
0.001020 T = 100 C
Linear (T = 100 C)
Linear (T = 20C)
0.001010
v = -5E-07*P + 0.001002
R2 = 0.9979
0.001000
0.000990
0 5 10 15 20 25
Pressure [MPa]
It is clear that
dv sat dv sat m3
5 10 7
dP dP
T 20 0 C T 100 0 C kg MPa
Therefore,
1
m3 m3
20 C 0.001002 5 10 7
kg kg MPa
1 1
20 C 4.99 10 4 4.99 10 10
MPa Pa
1
m3 m3
100 C 0.001044 5 10 7
kg kg MPa
1 1
100 C 4.79 10 4 4.79 10 10
MPa Pa
51
4.30
(a)
Substitute critical data into the Rackett equation:
vcalc
8.314190.6 0.29056 0.087750.0081 1111 / 190
2/7
5
46 10
cm 3
vcalc 38
mol
v
calc vexp
Error 100 %
vexp
Error 0.8%
(b)
cm 3
vcalc 54.2 Error 1.1%
mol
(c)
cm 3
vcalc 161.2 Error 10.8%
mol
(d)
cm 3
vcalc 20.5 Error 13.5%
mol
(e)
cm 3
vcalc 68.5 Error 19.7%
mol
1-hexanol had the largest percent error, while ethane and methane had the smallest percent
errors. The alkenes have low percentile errors due to their nonpolar nature, which the alcohols
and acids had large percentile errors due to being polar substances and exhibiting hydrogen
bonds.
52
4.31
(a)
Since we are given the temperature and pressure, we can make use of the generalized
compressibility charts. Using data from Table A.1.1, the required quantities can be found.
T 343.15 K
Tr 1.12
Tc 305.4 K
P 30 bar
Pr 0.616
Pc 48.74 bar
w 0.099
z 0 0.8376
z 1 0.0168
and
z z 0 wz 1 0.8376 0.099 0.0168 0.8393
Therefore,
30 kg
0.8393 J
8.314 343.15 K
znRT 0.03007 kg/mol mol K
V
P
30 10 5
Pa
V 0.796 m 3
(b)
The Redlich-Kwong EOS should give reasonably accurate results. Room temperature was
assumed to be 25 C. The molar volume is required for the calculations, so
v V
V
0.1 m 3
7.518 10 5
m3
n 40 kg mol
0.03007 kg/mol
Using data from Table A.1.1 and Equations 4.47 and 4.48,
J m 3 K1 / 2
a 9.88
2
mol
53
3
5 m
b 4.51 10
mol
Substituting these values into the Redlich-Kwong EOS and evaluating gives
P 191.3 bar
54
4.32
(a)
For propane
MW propane 0.0441 kg
mol
50 kg
n 1130 mol
kg
0.0441
mol
V 0.870 m 3
(b)
The Redlich-Kwong EOS is
RT a
P
v b T 1 / 2 vv b
where
0.42748R 2Tc2.5 J m 3 K1 / 2
a 18.33
Pc mol 2
0.08664 RTc m3
b 6.28 10 5
Pc mol
(Note: Critical data for propane obtained in Table A.1.1.)
Now that these values are known, there is only one unknown in the Redlich-Kwong EOS: v .
Using a numerical technique, e.g., the solver function on a graphing calculator
m3
v 0.000109
mol
V 0.124 m 3
55
(c)
The Peng-Robinson EOS is
RT a T
P
v b vv b bv b
where
1 1 Tr 2
Tr
T
50 273.15 K 0.873
Tc 370 K
0.37464 1.542260.152 0.269920.11522 0.605
1.081
0.45724RTc 2 J m3
a 1.02
Pc mol
0.07780 RTc m3
b 5.64 10 5
Pc mol
m3
v 0.0000942
mol
V 0.107 m 3
(d)
We must calculate the reduce temperature and pressure to use the compressibility charts:
T 323.15 K
Tr 0.873
Tc 370 K
P 35 bar
Pr 0.825
Pc 42.44 bar
z 0 0.1349
z 1 0.052
Therefore,
56
z z 0 wz 0 0.1349 0.152 0.052 0.127
Pv
RT
0.127 RT m3
v 0.00009749
P mol
V 0.111 m 3
(e)
From ThermoSolver
Using the Peng-Robinson equation,
m3
v 0.00009429
mol
V 0.107 m 3
m3
v 0.0000971
mol
Therefore,
V 0.110 m 3
57
4.33
Note: Multiple possibilities exist for which substance to use in the vial. The solution using one
possibility is illustrated below.
(a)
The substance must have a critical temperature above room temperature but below the
temperature of ones hand. From the Appendix A.1.2, we that for carbon dioxide
Tc 304.2 K 31.1 C
(b)
The vial must be able to withstand the pressure of the substance at its critical point. Therefore,
the vial must withstand 73.76 bar (the critical pressure of carbon dioxide).
(c)
Since the substance passes through its critical point, the molar volume at that state is constrained
by the critical temperature and pressure. To estimate the molar volume, we can use the Peng-
Robinson EOS. Calculate the necessary parameters for the EOS:
2
J
8.314 304.2 K 2
mol K J m3
a 0.45724 0.40
73.76 10 5 Pa 2
mol
J
8.314 304.2 K
mol K m3
b 0.07780 3.4 10 4
73.76 10 Pa
5
mol
1
Substitute the above parameters and solve for the molar volume:
m3 cm 3
vc 1.18 10 4 118
mol mol
n
V
100 cm 3
0.847 mol
vc cm 3
118
mol
58
(d)
If the vial contains less substance than needed, the molar volume will be greater. Therefore, the
substance will not pass through the critical point. As the substance is heated, it will go through a
transition from a saturated liquid and vapor mixture to superheated vapor. Then, it will become a
supercritical fluid once the critical temperature is exceeded.
critical
point
liquid
part d
not enough
part c CO2
vapor
liquid - vapor
59
4.34
Methane:
The following quantities are required to calculate the molar volume with the Peng-Robinson
equation.
T 1 1 Tr 2
T 0.9626
0.45724RTc 2 J m3
a 0.25
Pc mol
0.07780 RTc m3
b 2.68 10 5
Pc mol
T Tr Tc 209.66 K
P Pr Pc 55.2 10 5 Pa
RT a T
P
v b vv b bv b
4
m3
v 1.83 10
mol
Therefore,
5
55.2 10 Pa 1.83 10
3
4 m
mol
z 0.58
J
8.314 209.66 K
mol K
z 0 0.5984
z 1 0.0897
Therefore,
z z 0 wz 1 0.5984 0.0080.0897
z 0.599
60
Methanol:
m3
v 3.14 10 4
mol
z 0.651
z 0 0.5984
z 1 0.0897
w 0.559
Therefore,
Summary: Methane
z 0.58 (Peng-Robinson)
z 0.599 (Compressibility charts)
The value from the Peng-Robinson EOS is 3.2% smaller than the value
from the charts.
Methanol
z 0.651 (Peng-Robinson)
z 0.649 (Compressibility charts)
The value from the Peng-Robinson EOS is 0.31% smaller than the value
from the charts.
61
4.35
From Table A.1.1:
Tc 405.6 K
Pc 112.77 bar
w 0.25
Tr
T
92 273.15 K 0.9
Tc 405.6 K
P 306.5 bar
Pr 2.718
Pc 112.77 bar
z 0 0.4133
z 1 0.1351
Therefore,
z z 0 wz 1 0.38
Pv
RT
and
0.38RT m3
v 3.76 10 5
P mol
Since the temperature of the ammonia in this system is below the critical temperature, we know
the ammonia is not a supercritical fluid. The pressure is greater than the critical pressure, so the
ammonia is a liquid.
62
4.36
Let the subscript ace represent acetylene and but represent n-butane. First, convert the given
quantities of acetylene and n-butane into moles.
30 kg
nace 1152.2 mol
-3 kg
26.038 10 mol
50 kg
nbut 860.2 mol
-3 kg
58.123 10 mol
Therefore,
y ace 0.573
ybut 0.427
We can also calculate a and b parameters for pure species using the following equations
0.42748 R 2Tc2.5
a
Pc
0.08664 RTc
b
Pc
J K1/2 m3 m3
aace 8.03 bace 3.62 10 5
2
mol mol
J K1/2 m 3 m3
abut 29.07 bbut 8.08110 5
2
mol mol
Now, we can use mixing rules to calculate the parameters for the mixture.
63
5
m3
bmix y1b1 y 2b2 5.53 10
mol
Substitution of the mixture parameters into the Redlich-Kwong EOS results in an equation with
one unknown.
m3
v 0.00111
mol
m3
V nace nbut v 1152.2 mol 860.2 mol 0.00111
mol
V 2.23 m 3
64
4.37
First, we need to calculate the a and b parameters for species (1) and (2). From Table A.1.1 and
A.1.2:
Tc 304.2 K
CO2 (1):
Pc 73.76 bar
0.42748R 2Tc2.5 J m 3 K1/2
a1 6.466
Pc mol
0.08664 RTc m3
b1 2.97 10 5
Pc mol
Tc 591.7 K
Toluene (2):
Pc 41.14 bar
0.42748R 2Tc2.5 J m 3 K1/2
a2 61.17
Pc mol
0.08664 RTc m3
b1 1.036 10 4
Pc mol
2 mol m 3 3 mol m3
1.036 10
5 4
bmix 2.97 10
5 mol mol 5 mol mol
3
5 m
bmix 7.40 10
mol
Now we can substitute the known values into the Redlich-Kwong EOS. Note that the molar
volume can be calculated as follows
v
V
0.01 m 3 m3
0.002
ntot 5 mol mol
J m 3 K1/2
a mix 31.0
mol
65
From the mixing rules, we know
J m3 K1/2 2 2
amix 31.0 y1 a1 2 y1 y2 a12 y2 a2
mol
Therefore,
J m 3 K1/2
a12 16.55
mol
k12 0.17
66
4.38
(a)
We can match the species by looking at the b parameter alone. The b parameter is directly
related to the size of the molecule. Because sizeC 2 H 6 size H 2 O size H 2 ,
bC 2 H 6 bH 2 O bH 2 . Therefore,
J m3 m3
a b Species
mol mol
0.564 6.38 10 5 C2H6 (3)
0.025 2.66 10 5 H2 (1)
0.561 3.05 10 5 H2O (2)
These values are also consistent with what we would expect for the magnitude of van der Waals
interactions given by a.
(b)
The van der Waals parameters for the mixture can be calculated according to the mixing rules.
From Equations 4.81 and 4.82
n1 5 mol
y1 0.5
n1 n2 n3 10 mol
y 2 0 .4
y3 0.1
Since binary interaction parameters are not available, we must use Equation 4.78 to calculate
a12 , a13 , a23 .
J m3
a12 a1a2 0.118
mol
J m3
a13 0.119
mol
67
J m3
a23 0.562
mol
Now we can find the numerical values for the van der Waals parameters.
J m3
amix 0.206
mol
m3
bmix 3.19 105
mol
(c)
For the mixture, the van der Waals equation is
RT a
P mix
v bmix v2
v
V
0.0125 m 3 m3
0.00125
n1 n2 n3 10 mol mol
Therefore,
J m3
J
300 K
8.314 0.206
mol K
mol
P
3 3 2
0.00125 m 3.19 10 5 m m3
0.00125
mol mol
mol
P 1.92 MPa
68
4.39
(a)
To calculate the molar volume of the mixture, we will use the virial expansion in pressure since
it is more accurate at moderate pressures. But first, we must calculate the second virial
coefficient for the mixture.
Therefore,
RT Bmix RT
v 1 P Bmix
P RT P
cm 3 m3
v 2445.6 0.002446
mol mol
(b)
We can estimate k12 using an EOS where the a parameter is the only unknown. b can be
calculated as follows
m3 m3 3
bmix 0.25 5.83 10 5 0.75 2.97 10
5
3.69 10
5 m
mol mol mol
Substitute the known properties into the Redlich-Kwong EOS and solve for a.
J K1/2 m 3
amix 8.693
2
mol
69
J K1/2 m 3
a1 28.03
2
mol
(Calculated using Equation 4.47)
J K1/2 m 3
a2 6.466
2
mol
J K1/2 m 3
a12 8.81
2
mol
70
4.40
During this process the gas in Tank A undergoes an expansion from 7 bar to 3.28 bar. During
this isentropic process, the gas cools. If attractive forces are present, they will manifest
themselves more in state 1 at the higher pressure than in state 2. (While T has some effect to
counter this trend, we expect it will be secondary to the effect of pressure). Thus, in the non-
ideal case, additional energy must be supplied to pull the molecules apart. Since the tank is
insulated, the only place that this can come from is the kinetic energy of the molecules. Thus,
they slow down even more than in the ideal gas case, and the final temperature is lower. If we
had an equation of state that appropriately described the non-ideal behavior, we could apply the
concepts of the thermodynamic web that we have learned in this chapter to solve for the final
temperature.
71
4.41
Select the Equation of State Solver from the main menu. Enter the pressure and temperature
and solve for the molar volume using the Lee-Kesler EOS. The program provides
m3
v 1.07 10 4
mol
To calculate the volume occupied by 10 kilograms of butane, we need to know the number of
moles present in 10 kg.
10 kg
n 172.1 mol
0.05812 kg/mol
Therefore,
V 172.1 mol1.07 10 4
m3
0.018 m
3
mol
The answers provided in Example 4.9 agree well with this answer. The answer found using the
Redlich-Kwong EOS is 16.7% larger than the answer from ThermoSolver, but the answer from
the compressibility charts is only 5.56% larger.
72
4.42
Using ThermoSolver should be straightforward; thus, only the answers are provided.
(a)
Tc 305.4 K
Pc 48.84 bar
kJ
h f ,298 84.68
mol
(b)
Tsat 299.42 K
(c)
The percent difference is calculated as follows
v LK v PR z LK z PR
% Difference or 100 %
v LK z LK
(i).
(i).
Quantity Lee Kessler EOS Peng Robinson EOS % Difference
m3 m3
v 3.68 10 4 3.58 10 4 2.7
mol mol
z 0.5859 0.5711 2.5
73
Chapter 5 Solutions
Engineering and Chemical Thermodynamics
Wyatt Tenhaeff
Milo Koretsky
(a)
Following the example given by Equation 5.5a in the text
u u
du dT dP
T P P T
(b)
u u
du dT ds
T s s T
(c)
u u
du dh ds
h s s h
2
5.2.
The internal energy can be written as follows
u u
du dT dv
T v v T
u u P
cv and T P
T v v T T v
P
du cv dT T P dv
T v
P R
T v v
Therefore,
RT
du cv dT P dv
v
RT
which upon noting that P for an ideal gas, becomes
v
du cv dT
3
5.3
The heat capacity at constant pressure can be defined mathematically as follows
h u Pv u v
cP P
T P T P T v T P
v R
T P P
Therefore,
u
cP R
T v
u u
du dT dP
T P P T
u
We can now rewrite :
T v
u u T u P
cv
T v T P T v P T T v
u
0
P T
so
u
cv
T P
c P cv R
4
5.4
In terms of P, v, and T, the cyclic equation is
P T v
1
T v v P P T
Pv RT
P R
T v v
T P
v P R
v RT v
2
P T P P
Therefore,
P T v R P v
1
T v v P P T v R P
5
5.5
For a pure species two independent, intensive properties constrains the state of the system. If we
specify these variables, all other properties are fixed. Thus, if we hold T and P constant h cannot
change, i.e.,
h
0
v T , P
6
5.6
Expansion of the enthalpy term in the numerator results in
h Ts vP
T s T s
h P
v
T s T s
h s
v
T s v P
h s T
v
T s T P v P
s c
P (use thermodynamic web)
T P T
T 1 RT 2a 2ab
(differentiate van der Waals EOS)
v P R v b v 2 v 3
Therefore,
h vc RT 2a 2ab v 2a b
P c P 1
T s RT v b v 2 3
v v b vRT v
7
5.7
h
:
P T
h Ts vP s
T v
P T P T P T
s
P T
v R
1 B' P C ' P
T P P
2
h
T
P
RT
P
1 B ' P C ' P 2 v v v
h
0
P T
h
:
P s
h Ts vP
v
P s P s
h
P s
RT
P
1 B' P C ' P 2
h
:
T P
h
cP (Definition of cP)
T P
h
:
T s
h Ts vP P
v
T s T s T s
P s P c T c P
T s
P
T P s T
P
T v P T R 1 B' P C ' P 2
8
h Pv 1
1 B' P C ' P 2
T s
c P
RT 1 B' P C ' P 2 c P
1 B' P C ' P 2
h
cP
T s
9
5.8
(a)
A sketch of the process is provided below
well m
insulated
s = 0
T1
P1 T2
P2
The diagram shows an infinitesimal amount of mass being placed on top of the piston of a
piston-cylinder assembly. The increase in mass causes the gas in the piston to be compressed.
Because the mass increases infinitesimally and the piston is well insulated, the compression is
reversible and adiabatic. For a reversible, adiabatic process the change in entropy is zero.
Therefore, the compression changes the internal energy of the gas at constant entropy as the
pressure increases.
(b)
To determine the sign of the relation, consider an energy balance on the piston. Neglecting
potential and kinetic energy changes, we obtain
U Q W
U W
As the pressure increases on the piston, the piston compresses. Positive work is done on the
system; hence, the change in internal energy is positive. We have justified the statement
u
0
P s
10
5.9
(a)
By definition:
1 v
v T P
and
1 v
v P T
Dividing, we get:
v
T P v P
v T P v T
P T
v P T
1
T P v T P v
Hence,
P
T v
(b)
If we write T = T(v,P), we get:
T T
dT dv dP (1)
v P P v
cv P c v
ds dT dv P dT dP
T T v T T P
11
T P T v
dT dv dP (2)
c P cv T v c P cv T P
For Equations 1 and 2 to be equal, each term on the left hand side must be equal. Hence,
T T P
v P c P cv T v
or
P v v
c P cv T T
T v T P T P
where the result from part a was used. Applying the definition of the thermal expansion
coefficient:
P v Tv 2
c P cv T
T v T P
12
5.10
We need data for acetone, benzene, and copper. A table of values for the molar volume, thermal
expansion coefficient and isothermal compressibility are taken from Table 4.4:
m3
Species v 10 6 103 K -1 1010 Pa -1
mol
Acetone 73.33 1.49 12.7
Benzene 86.89 1.24 9.4
Copper 7.11 0.0486 0.091
We can calculate the difference in heat capacity use the result from Problem 5.9b:
vT 2
c P cv
or
73.33 10 6
m3
293 K 1.49 10 K
3 -1 2
mol J
c p cv
12.7 10 10
Pa -1
37.6
mol K
J J
c p cv cp
mol K mol K
Species % difference
We can compare values to that of the heat capacity given in Appendix A2.2. While we often
assume that cP and cv are equal for condensed phases, this may not be the case.
13
5.11
We know from Equations 4.71 and 4.72
1 v 1 v
and
v T P v P T
Maxwell relation:
s P
v T T v
P P v
T v v T T P
Therefore,
s
v T
Maxwell Relation:
s v
P T T P
v
v
T P
Therefore,
s
v
P T
14
5.12
(a)
An isochor on a Mollier diagram can be represented mathematically as
h
s v
h Ts vP P
T
s v s v s v
h T s
T v
s v s v v T
We know
T T s P
and
s v c v v T T v
s P R
v T T v v
Therefore,
h T R R
T v T 1
s v cv v cv
(b)
In Part (a), we found
h T P
T v
s v cv T v
15
P R
T v v b
Therefore,
h RT v
T
s v cv v b
16
5.13
(a)
The cyclic rule can be employed to give
T T s
P s s P P T
T T v
P s c P T P
v R
T P P
Therefore,
T RT 1 v
P s P cP cP
(b)
Separation of variables provides
T R P
T cP P
Integration provides
R
T P cP
ln 2 ln 2
T1 P1
R
T2 P2 cP
T1 P1
17
R
P2 v 2 P2 cP
P1v1 P1
R R
1 1
P2 P v P1 P v
c c
2 1
where
R c P R cv 1
1
cP cP cP k
P2 v 2k P1v1k
Pv k const.
(c)
In Part (a), we found
T T v
P s c P T P
Using the derivative inversion rule, we find for the van der Waals equation
v Rv 3 v b
T P RTv 3 2av b 2
Therefore,
T 1 RTv 3 v b
P s c P RTv 3 2av b 2
18
5.14
The development of Equation 5.48 is analogous to the development of Equation E5.3D. We
want to know how the heat capacity changes with pressure, so consider
c P
P T
c P h h
P T P T P T T P T P
h
Consider the term:
P T
h Ts vP s v
T v T v
P T P T P T T P
c
Substitution of this expression back into the equation for P results in
P T
c P v
T v
P T T T P P
c P T v 2v
T v
P T T T P
T P T P
2
c P 2v
T 2
P T T
P
Therefore,
c Preal Preal 2v
dc P T dP
T
2
c Pideal Pideal P
and
Preal
2v
c Preal c P
ideal
T dP
2
Pideal T
P
19
5.15
In order to solve this problem we need to relate the change in entropy from 10 to 12 bar to the
change in molar volume (for which we have complete data). First, we can rewrite the change in
entropy as
12 bar
s
s s 2 s1 dP
P T
10 bar
Applying a Maxwell relation, we can relate the above equation to the change in molar volume:
12 bar 12 bar
s v
s 2 s1 P T dP s1 T P dP
10 bar 10 bar
As 10 bar:
v v 4
m3
5.60 10
T P T P kg K
At 12 bar:
v v 4
m3
4.80 10
T P T P kg K
v
To integrate the above entropy equation, we need an expression that relates to pressure.
T P
Thus, we will fit a line to the data. We obtain
v m3 3
4.0 10 P 9.6 10
10 4 m
T P
kg K Pa kg K
20
5.16
A schematic of the process follows:
Propane in
We also know the ideal gas heat capacity from Table A.2.1:
cP
1.213 28.785 10 3 T 8.824 10 6 T 2
R
Since this process is isentropic (s=0), we can construct a path such that the sum of s is zero.
Ideal
step 2 Gas
v2 ,T2
s2
volume
s1
step 1
s=0
v1 ,T 1
Temperature
or in mathematical terms:
s s
ds dT dv 0
T v v T
21
c P
ds v dT dv
T T v
To get s1
RT a P R
P 2 so
vb v T v v b
and
dv 2 R dv R ln v2 b
v 2 P v
s1 ds
v1 T v v1 v b v1 b
or, using the ideal gas law, we can put s1 in terms of T2:
RT2 b
P
s1 R ln 2
v1 b
For step 2
s s1 s 2 0
RT2
P b
ln 2 0.213 ln
v1 b
T2
623.15
28.785 10 3
T2 623.15 K
8.824 10 6 2
2
T2 623.15 K 2
Substitute
T1 623.15 K
v1 600 cm 3 /mol
P2 1 atm
cm 3 atm
R 82.06
mol K
22
and solve for T2:
T2 448.3 K
Choosing T and P as the independent variables, (and changing T under ideal gas conditions), we
get:
P1,T1
s=0
Pressure
s1
step 1
Ideal
step 2
Gas
P2,T2 s2
Temperature
s s
ds dT dP 0
T P P T
c v
ds P dT dP 0
T T P
R
v b
v
T P RT 2a
v b
2
v 3
Therefore,
23
R
T2 P2 v b
cP
s T
dT RT 2a
dP 0
P1
T1
v b 2 v 3
We cant integrate the second term of the expression as it is, so we need to rewrite dP in terms
of the other variables. For the van der Waals equation at constant temperature:
2a RT
dP dv
v
3
v b 2
Substituting this into the entropy expression, we get
Upon substituting
T1 623.15 K
v1 600 cm 3
RT2
v2 (gas acts ideally at 1 atm)
P2
cm 3
b 91
mol
cm 3 atm
R 82.06
mol K
T2 448.3 K
24
5.17
(a)
Attractive forces dominate. If we examine the expression for z, we see that at any absolute
temperature and pressure, z 1. The intermolecular attractions cause the molar volume to
deviate negatively from ideality and are stronger than the repulsive interactions.
(b)
Energy balance:
h2 h1 q
P [bar]
P,T1 P,T2
q = h
50
step 1
step 3
h1 h
3
h2 step 2
0
ideal gas
h h
dh dT dP
T P P T
v
dh c P dT T v dP
T P
1
v R aT 1 / 2
P
25
Taking the derivative gives:
v R 0.5
0.5aRT
T P P
so
dh cP dT 0.5aRT 0.5 dP
For step 1
3.58 3.02 10
500 K
3 J
h2 R T 0.875T 0.5 dT 7961
300 K
mol
For step 3:
J
q h1 h2 h3 7888
mol
h c Preal
P real 2v
c Preal c Pideal T 2
T dP
P ideak
P
For the given EOS
26
1
v R aT 1 / 2
P
Therefore,
2v
0.25aRT 1.5
T 2
P
and
P real 2v P 50 bar
real
T 2
T dP
0.25aRT 0.5 dP 0.875 K 1/2 RT 0.5
P ideak
P P ideak 0 bar
We can combine this result with the expression for c Preal and find the enthalpy change.
3.58 3.02 10
500 K
3
h R T 0.875T 0.5 dT
300 K
J
q h 7888
mol
27
5.18
(a)
Calculate the temperature of the gas using the van der Waals equation. The van der Waals
equation is given by:
RT a
P 2
v b v
V
v1 1
Al 0.1 m 2 0.4 m m3
0.00016
n n 250 mol mol
10000 kg 9.81 m2
s
1.01325 10 5 Pa 1.08 10 6 Pa
mg
P1 Patm
A
0.1 m 2
Substituting these equations into the van der Waals equation above gives
J m3
J
8.314 T1 0.5
mol K mol
1.08 10 Pa
6
m3 5 m
3
3
2
0.00016 4 10 0.00016 m
mol mol mol
T1 297.5 K
Since the process is isothermal, the following path can be used to calculate internal energy:
v2 ,T2 = T 1
v
u
s
v1 ,T1
28
u u u
du dT dv dv
T v v T v T
v2
P
u T T v P dv
v1
RT a
P 2
v b v
so
P R
T v v b
Therefore,
v2
a
u v 2 dv
v1
We can assume the gas in state 2 is an ideal gas since the final pressure is atmospheric.
Therefore, we calculate v 2 ,
RT2 m3
v2 0.0244
P2 mol
and
J m3
0.0244
0.5
mol J
u dv 3104.5
0.00016 v2 mol
or
U 776.1 kJ
29
(b)
From the definition of entropy:
s s
ds sys dT dv
T v v T
s
ds sys dv
v T
v2
P
s sys dv
T v
v 1
P R
T v v b
v2
R
s sys v b dv
v1
J
8.314
mol K dv 44.17 J
0.0244
s sys mol K
3
0.00016 v 4 10 5 m
mol
J
S sys 11042.5
K
Qsurr
s surr
Tsurr
where
30
Qsurr Q (Q is the heat transfer for the system)
Q U W
We know the change in internal energy from part a, so lets calculate W using
v2
W n Pdv
v1
Therefore,
1.39 10 6 J J
S surr 4672
297.5 K K
J J J
S univ 11042.5 4672 6370.5
K K K
31
5.19
First, calculate the initial and final pressure of the system.
To find the final temperature, we can perform an energy balance. Since the system is well-
insulated, all of the work done by adding the third block is converted into internal energy. The
energy balance is
u w
To find the work, we need the initial and final molar volumes, which we can obtain from the
given EOS:
vi 8.37 10 4 m 3 /mol
vf
8.314T f
3.2 10 -5 m 3 /mol
6.89 10 6 1
25
Tf
step 3
vi ,Ti
u = -P (v -v )
f f i
vf ,Tf
500 Tf T
32
For step 1, we calculate the change in internal energy as follows
v RT / Plow v RT / Plow
u P
u1 dv
v T
T P dv
T v
vi vi
v RT / Plow 2
aRTi dv aRTi ln RTi / Plow b
2
u1 T a 2 v b T a 2 vi b
vi i i
vf v
u
f
P
u3 v T dv T v dv
T P
v RT / Plow v RT / Plow
vf aRT 2 dv aRT f 2 vf b
T a 2 v b T a 2 RT f / Plow b
f
u3
ln
v RT / Plow f f
Insert the expression for the final molar volume into the equation for u3 :
aRT f 2 8.314T f
ln
u3
T f a 2 6.89 10 6 1 25 / T f RT f / Plow b
Since the pressure is low (molar volume is big) during the second step, we can use the ideal heat
capacity to calculate the change in internal energy.
Tf Tf
u 2 cv dT 20 R 0.05T dT
Ti 500 K Ti 500 K
u 2 11.686 T f 500 0.025 T f2 500 2
If we set the sum of the three steps in the internal energy calculation equal to the work and
choose an arbitrary value for Plow, 100 Pa for example, we obtain one equation with one
unknown:
33
aRTi 2 RT / P b
ln i low
Ti a vi b
2
11.686 T f 500 0.025 T f2 500 2
aRT f 2 8.314T f
ln
Tf a 2
6
6.89 10 1 25 / T f RT f / Plow b
6.89 10 6 Pa 8.314T f
3.2 10 -5 8.37 10 4
6.89 10 1
6
25
Tf
T f 536.2 K
suniv s sys
Since the gas in the cylinder is not ideal, we must construct a hypothetical path, such as one
shown below, to calculate the change in entropy during this process.
P
Pf ,Tf
6.9
4.9 Pi ,Ti
step 3
step 1
ideal gas
Plow
step 2
500 536 T
Plow Plow
s v
s1 dP dP
Pi
P T Pi T P
34
Pf P
s
f
v
s3 dP dP
P T P
Plow
T
Plow
For step 2
Tf Tf Tf
s c 20
s 2 dT P dT 0.05 dT
Ti
T P Ti
T
Ti
T
Tf
s2 20 ln
0.05 T f Ti
i
T
Sum all of the steps to obtain the change in entropy for the entire process
suniv s sys s1 s 2 s3
suniv
RTi 2a Ti Plow
ln
Tf
20 ln
0.05 T f Ti
RT f 2a T f
Pf
T ln P
a Ti 2 Pi i a T f
2
low
Arbitrarily choose Plow (try 100 Pa), substitute numerical values, and evaluate:
J
suniv 0.388
mol K
J J
S univ 2 mol 0.388 0.766
mol K K
35
5.20
A schematic of the process is given by:
well
insulated
V=1L V=2L
V=1L
T = 500 K T=?
Vacuum
nCO =1 mole
nCO=1 mole
State i State f
(a)
The following equation was developed in Chapter 5:
v 2P
cvreal cvideal T T 2 dv
v ideal v
2P
0
T 2
Therefore,
cvreal cvideal
3100 J
cvreal R 3.376 5.57 10 4 500 R 22.0
500 2 mol K
(b)
As the diaphragm ruptures, the total internal energy of the system remains constant. Because the
volume available to the molecules increases, the average distance between molecules also
increases. Due to the increase in intermolecular distances, the potential energies increase. Since
the total internal energy does not change, the kinetic energy must compensate by decreasing.
Therefore, the temperature, which is a manifestation of molecular kinetic energy, decreases.
36
(c)
Because the heat capacity is ideal under these circumstances we can create a two-step
hypothetical path to connect the initial and final states. One hypothetical path is shown below:
v
vf ,Tf
step 2
2
s2
step 1
vi ,Ti
u1, s1
Tf 500 T
Tf Tf
u1 cvreal dT cvideal dT
Ti Ti
Tf
4 3100
u1 R
2.376 5.57 10 T 2 dT
500 K i T
u1 2.32 10 3 T f2 19.75T f
25773.4
Tf
10507.4
vf
u
u 2 v T dv
vi
vf
P
u 2 T v dv
T P
vi
RT a
For the van der Waals EOS, P 2 ,
v b v
37
P R
T v v b
Therefore,
v f 0.002 v f 0.002
RT a J
u1 v b P dv v 2 dv 73.7 mol
v i 0.001 v i 0.001
Now set the sum of the two internal energies equal to zero and solve for Tf:
u1 u 2 2.32 10 3 T f2 19.75T f
25773.4
Tf
10507.4 73.7 0
T f 497 K
(d)
Since the system is well-insulated
suniv s sys
s s
ds sys dT dv
T v v T
s cv
T v T
s P
v T T v
P R
T v v b
38
497 K 0.002
cv R
s sys T
dT v b dv
500 K 0.001
497 K 2.376
0.002
3100 dv
s sys R 5.57 10 4
dT
500 K T 5
T3 0.001 v 3.95 10
J
suniv s sys 5.80
mol K
39
5.21
A schematic of the process is given by:
well
insulated
State i State f
Energy balance:
u 0
Because the gas is not ideal under these conditions, we have to create a hypothetical path that
connects the initial and final states through three steps. One hypothetical path is shown below:
step 2
v [m 3/mol] ideal gas
step 3
step 1
0.2/400 vf ,Tf
0.1/400 vi ,Ti
Tf 300 T [K]
v
u
u1 dv
v T
vi
v
P
u1 T P dv
T v
vi
40
For the van der Waals EOS
P R a
T v v b T v 2
2
Therefore,
v v
RT a 2a J
u1 4 v b Tv 2 dv 4 T v 2 dv 1120 mol
P
vi 2.510 vi 2.510 i
v f 510 4 v f 510 4
RT a 2a 168000 J
u3 v b 2 P dv
Tv
Tf v 2
dv
Tf mol K
v v
For step 2, the molar volume is infinite, so we can use the ideal heat capacity given in the
problem statement to calculate the change in internal energy:
u 2
3
2
R T f 300 K
If we set sum of the changes in internal energy for each step, we obtain one equation for one
unknown:
J 3
u1 u 2 u3 1120
J
8.314 T f 300 K
168000
0
mol 2 mol K Tf
T f 261.6 K
41
5.22
A schematic of the process is shown below:
T surr = 293 K
Initially:
vacuum
(a)
Consider the tank as the system. Since kinetic and potential energy effects are negligible, the
open system, unsteady-state energy balance (Equation 2.47) is
dU
nin hin n out hout Q W s
dt sys in out
The process is adiabatic and no shaft work is done. Furthermore, there is one inlet stream and no
outlet stream. The energy balance reduces to
dU
nin hin
dt sys
U2 t
dU n
U1 0
h dt
in in
t
n2 u 2 n1u1 hin n in dt nin hin (n2 n1 )hin
0
Since the tank is initially a vacuum, n1=0, and the relation reduces to:
u 2 hin
42
As is typical for problems involving the thermodynamic web, this problem can be solved in
several possible ways. To illustrate we present two alternatives below:
or
J
Pin vin RT 1 B' P 8.314
552 K 1 2.8 10 8 3 10 6 Pa 3,800 J
mol (2)
mol K
The change in internal energy can be found from the following path:
3 MPa
step 1
step 3
u1 u3
u2
Plow
step 2 ideal gas
500 K T2 T
For steps 1 and 3, we need to determine how the internal energy changes with pressure at
constant temperature: From the fundamental property relation and the appropriate Maxwell
relation:
u s v v v
T P T P
P T P T P T T P P T
u
RT
1 B'P P RT2 B' RT
P T P P
43
So for step 1:
u
0 0
u1 dP B ' RTdP B ' RT Pin 349 [J/mol] (3)
Pin
P T Pin
u
P2 P2
u 3 dP B ' RTdP B ' RT P2 0.7T (4)
0
P T 0
For step 2
u
T T
h Pv T
u 2 dT dT c P R dT
500
T P 500
T P T P 500
or
0.131 19.225 10
T2
3
u 2 R T 5.56110 6 T 2 dT (5)
T1 500 K
T2 552 K
u 2 hin
h2 P2 v2 hin
h2 hin P2 v2
To calculate the enthalpy difference, we can use the real heat capacity
P 2
v
c Preal c P
ideal
T dP
2
Pideal T P
44
2v
0
T 2
P
Therefore,
c Preal c P
ideal
Now, we can calculate the change in enthalpy and equate it to the flow work term.
T2
cP
ideal
dT P2 v2
T1 500 K
1.131 19.225 10
T2
R 3
T 5.561 10 6 T 2 dT P2 v2 RT2 1 B' P2
T1 500 K
T2 552 K
(b)
In order to solve the problem, we will need to find the final pressure. To do so, first we need to
calculate the molar volume. Using the information from Part (a) and the truncated virial
equation to do this
J
8.314 552 K
v
RT
P
1 B' P mol K
3 10 Pa6
1 2.8 10 3 10 Pa
8 6
m3
v 0.0014
mol
This quantity will not change as the tank cools, so now we can calculate the final pressure.
m3
P2 0.0014
J
mol
1 2.8 10 8 P2
8.314 293 K
mol K
Solve for P2 :
45
P2 1.66 10 6 Pa
To solve for the change in entropy of the system start with the following relationship:
s s
ds sys dT dP
T P P T
cP v
ds sys dT dP
T T P
3 MPa
s3
step 3
1.66 MPa
step 1
s 1
s 2
Plow
step 2 ideal gas
500 K T2 T
v
dP 1 B P dP
1 Pa 1 Pa
R
s1 '
1.66 MPa
T P 1.66 MPa
P
For step 3,
v
3 MPa 3 MPa
R
s1
1 Pa
dP
T P
1 Pa
P
1 B ' P dP
46
For step 2:
1.131
293 K 293 K
cP
s 2
552 K
T
dT R
552 K
T
19.225 10 3 5.561 10 6 T dT
J
s sys 46.9
mol K
______________________________________________________________________________
c real v
ds sys P dT dP
T T P
v 1
R B '
T P P
Now, substitute the proper values into the expression for entropy and integrate:
293 K 1.6610 6 Pa
1.131 3 6 1
s sys R T 19.225 10 5.561 10 T dT R B' dP
P
552 K 310 ^ Pa
J
s sys 46.9
mol K
______________________________________________________________________________
In order to calculate the change in entropy of the surroundings, first perform an energy balance.
u q
h Pv q
47
Since the real heat capacity is equal to ideal heat capacity and the molar volume does not change,
we obtain the following equation
Tf
cP
ideal
dT v P f Pi q
Ti
T f 293K
R 1.131 19.225 10 3
m3
T 5.561 10 6 T 2 dT 0.0014
1.66 10 Pa - 3 10 Pa q
6 6
Ti 552 K mol
J
q 15845
mol
Therefore,
J
q surr 15845
mol
and
J
s surr 54.08
mol K
Before combining the two entropies to obtain the entropy change of the universe, find the
number of moles in the tank.
n
0.05 m 3 75.7 mol
m3
0.0014
mol
J J
S univ 75.7 mol 54.08 46.9
mol K mol K
J
S univ 544
K
48
5.23
First, focus on the numerator of the second term of the expression given in the problem
statement. We can rewrite the numerator as follows:
uTr , v r uTideal
r ,vr
uTr , v r uTideal
r ,vr
uTideal
r ,vr
uTideal
r ,vr
For an ideal gas, we know
uTideal uTideal 0
r ,vr r ,vr
Therefore,
uTr , v r uTideal
,vr
uTr , v r uTideal
,vr
r r
Substitute this relationship into the expression given in the problem statement:
uTdep
,v
uTr , v r uTideal
,vr
uTr , v r uTideal
,vr
r r
r
r
Now, we need to find an expression for uTr , vr uTideal . Note that the temperature is constant.
r ,vr
P
duT T P dv
T v
zRT
P
v
and substituted into the expression for the differential internal energy
RT z RT zRT RT 2 z
duT T dv dv
v T v v v v T v
v
49
duTr T 2 z
r dvr
RTc vr Tr v r
v uTr , v r uTideal
duTr vr 2
T z
r, vv
RTc RTc
r
v r Tr v
dv r
v v r
Therefore,
uTdep
,v
uTr , v r uTideal
, v
uTr , v r uTideal
,vr
vr 2
T z
r r
r r
r
r dv r
RTc RTc RTc v r Tr v
v r
50
5.24
We write enthalpy in terms of the independent variables T and v:
h h
dh dT dv
T v v T
dh Tds vdP
P P
dhT T v dv
T v v T
P R 1 a
T v v b 2 T 3 / 2
vv b
P RT a a
v T v b 2 T 1 / 2 v 2 v b T 1 / 2 vv b 2
Therefore,
RT RTv 3 a a
dhT dv
v b v b 2 2 T 1 / 2 vv b T 1 / 2 v b 2
v v RT RTv 3 a a
dhT v b v b 2 2 T 1 / 2 vv b T 1 / 2 v b 2 dv
dep
h
v v
RTb 3a v a
h dep ln
v b 2bT 1 / 2 v b T 1 / 2 v b
To calculate the entropy departure we need to be careful. From Equation 5.64, we have:
sT , P sTideal
,P
gas
sT , P sTideal gas
ideal gas
, P 0 sT , P sTideal gas
, P 0
51
However, since we have a P explicit equation of state, we want to put this equation in terms of v.
Lets look at converting each state. The first two states are straight -forward
sT , P sT , v
and
For the third state, however, we must realize that the ideal gas volume v at the T and P of the
system is different from the volume of the system, v. In order to see this we can compare the
equation of state for an ideal gas at T and P
RT
P
v'
RT a
P
vb T vv b
The volume calculated by the ideal gas equation, v, is clearly different from the volume, v,
calculated by the Redlich-Kwong equation. Hence:
Thus,
sT , P sTideal
,P
gas
sT , v sTideal
gas ideal gas
, v sT , v sTideal
gas ideal gas
, v sT , v ' sTideal
,v
gas
ds P
dv T T v
Therefore,
P
dsT dv
T v
52
For the Redlich-Kwong EOS
P R 1 a
T v v b 2 T 3 / 2
vv b
so
s
v
R 1 a
ideal gas
T , v sT , v dv
v
v b 2 T 3/ 2
v v b
P R
T v v
so
v
R
sTideal
gas
,v sTideal gas
,v v dv
v
Finally:
v'
dv v' RT
ideal gas ideal gas
s ' sT , v R R ln R ln
T ,v v v Pv
v
s dep R ln
v b a v
ln R ln
RT
v 2bT 3 / 2 v b Pv
53
5.25
Calculate the reduced temperature and pressure:
Tc 647.3 K
Pc 220.48 bar (Table A.1.2)
w 0.344
300 bar
Pr 1.36
220.48 bar
673.15 K
Tr 1.04
647.3 K
( 0) (1)
h dep h dep
Tr , Pr 2.921 Tr , Pr 1.459
RTc RTc
( 0) (1)
s dep s dep
Tr , Pr 2.292 Tr , Pr 1.405
R R
dep
( 0)
h dep
(1)
hTr , Pr
h dep RTc w
Tr , Pr
RTc
RTc
J J
h dep 8.314 647.3 2.921 0.344 1.459 18421
mol K mol
dep
( 0)
s dep
(1)
sTr , Pr
s dep R w
Tr , Pr
R
R
J J
s dep 8.314 2.292 0.344 1.405 23.07
mol K mol K
54
To use the steam tables for calculating the departure functions, we can use the following
relationships.
h dep hT , P hTideal
,P
s dep sT , P sTideal
,P
kJ kJ
hT , P 2151.0 and sT , P 4.4728
kg kg K
We need to calculate the ideal enthalpies and entropies using the steam tables reference state.
673.15 K
hTideal
, P h
vap
0 .01 C
c ideal
p dT
273.16 K
kJ
We can get h vap 45.1 from the steam tables and heat capacity data from Table A.2.2.
mol
Using this information, we obtain
673.15 K
kJ kJ 3 0.121 10 5
hTideal
, P 45.1
0.008314
mol K
3.47 1.45 10
273.16 K
T
T2
dT
mol
kJ
, P 59.14
hTideal
mol
673.15 K c ideal
P
sTideal
, P s
vap
0 .01 C
p
dT R ln 2
273.16 K
T P1
kJ
s vap 0.01 C 0.165
mol K
55
kJ
673.15 K 5
3.47 3 0.121 10 30
sTideal
, P 0.165 0.008314
mol K 273.16 K T
1.45 10
T3
dT ln
0.000613
J
sTideal
, P 107 mol K
kJ kJ kJ
h dep 2151.0 0.0180148 kg/mol 59.14 20.4
kg mol mol
kJ
s dep 4.4728 0.0180148 kg/mol 0.107 kJ 0.0264 kJ
kg K mol K mol K
Table of Results
Generalized Percent Difference
Steam Tables
Tables (Based on steam tables)
kJ
h dep
mol
-18.62 -20.4 9.9
kJ
s dep
mol K
-0.0231 -0.0264 12.5
56
5.26
State 1 is at 300 K and 30 bar. State 2 is at 400 K and 50 bar. The reduced temperature and
pressures are
30 bar 50 bar
P1, r 0.616 P2,r 1.026
48.74 bar 48.74 bar
300 K 400 K
T1, r 0.982 T2,r 1.31
305.4 K 305.4 K
and
0.099
( 0) (1)
h dep h dep
T1, r , P1, r T1, r , P1, r
RT 0.825 RT 0.799
c c
( 0) (1)
h dep h dep
T2 , r , P2 , r T2 , r , P2 , r
RT 0.711 RT 0.196
c c
Therefore,
h dep
T1, r , P1, r
RT 0.825 0.099 0.799 0.904
c
h dep
T2 ,r , P2 ,r 0.711 0.099 0.196 0.730
RTc
The ideal enthalpy change from 300 K to 400 K can be calculated using ideal cP data from Table
A.2.1.
400 K
3
hTideal
T
R 1.131 19.225 10 T 5.56110 6T 2 dT 717.39 R
1 2
300 K
57
dep dep
h hT , P hTideal
T2
hT , P
1, r 1, r 1 2,r 2,r
s dep
T1, r , P1, r 0.601 0.099 0.756 0.676
R
s dep
T2 , r , P2 , r
R 0.394 0.099 0.224 0.416
s ideal 1.542 R
Therefore,
J
s 14.98
mol K
58
5.27
The turbine is isentropic. Therefore, we know the following
dep dep
s sT , P s ideal sT , P 0
1, r 1, r 2,r 2,r
Using the van der Waals EOS, we can find P1,r, which leaves one unknown in the above
equation: T2.
3 atm cm 3
82.06 cm atm 623.15 K 91 10 5
mol K
2
mol
P1
3 3 3
2
600 cm 91 cm 600 cm
mol mol mol
P1 75.19 atm 76.19 bar
Calculate reduced temperature and pressures using data from Table A.1.1
Also,
0.152
s dep
T1, r , P1, r 0.327 0.152 0.102 0.343
R
59
s dep
1.213 28.785 10 3 T 8.824 10 6 T 2
T2
T2 , r , P2 , r
R 0.343 T
dT 4.32 R s
623.15 K
We can solve this using a guess-and-check method
Therefore,
T2 446.6 K
60
5.28
A reversible process requires the minimum amount of work. Since the process is reversible and
adiabatic
s 0
dep dep
s sT , P s ideal sT , P 0
1, r 1, r 2,r 2,r
Calculate reduced temperature and pressures using data from Table A.1.1
1 bar 10 bar
P1, r 0.0217 P2,r 0.217
46.0 bar 46.0 bar
300 K
T1, r 1.57
190.6 K
s dep
T1, r , P1, r
R 0.00457 0.008 0.0028 0.0046
Therefore,
T2 3 6 2 s dep
1.213 28.785 10 T 8.824 10 T T2 , r , P2 , r
s R 0.0046 dT 2.303
T R
300 K
We can solve using a guess-and-check method
61
J
s 1.42
mol K
Therefore,
T2 379 K
h2 h1 h ws
We can calculate the enthalpy using departure functions. From Tables C.3 and C.4:
h dep
T1, r , P1, r
RT 0.0965 0.008 0.011 0.0966
c
h dep
T2 , r , P2 , r
RT 0.0614 0.00890.015 0.0613
c
Ideal heat capacity data can be used to determine the ideal change in enthalpy
379 K
h ideal R 1.702 9.081 10 3 T 2.164 10 6 T 2 dT
300 K
Therefore,
J
379 K
h 8.314 190.6 K 0.0966 0.0613 1.702 9.081 10 3
T 2.164 10 6 2
T dT
mol K 300 K
J
w s h 3034.2
mol
and
mol J
W S 1 / 30
3034.2 mol 101.1 W
s
62
5.29
Equation 4.71 states
1 v
v T P
v
v
T P
vT 1
JT
cP
63
5.30
For an ideal gas
v R
T P P
Therefore,
RT
v
JT
P v v 0
cP cP
64
5.31
The van der Waals equation is given by:
RT a
P (1)
v b v2
1 v 1 T
(2)
v T P v v P
a v b
T P 2
v R
T a 1 v b 2a Pv 3 av 2ab
P 2 (3)
v P v R R v3 Rv 3
Rv 2
Pv 3 av 2ab
Rv 2 v b
RTv 3 2av b
2
1 v 1 P
v P T v v T
P 2a RTv 3 2av b
2
RT
v T v b 2 v 3 v 3 v b
2
so
65
v 2 v b
2
RTv 3 2av b
2
v
T v
JT T P
Preal
v
2
c Pideal T 2 dP
T P
Pideal
2T T 1 T 2a 6av b
2 3
v P v b v b v P Rv
2
Rv 4
2T 2a v 3b
2 (5)
v P Rv 4
bRTv 3 2avv b
2
RTv 3 2av b
2
JT
RTv 4
Preal
c Pideal
2av 3b
dP
Pideal
At a given temperature the integral in pressure can be rewritten in terms of volume using the van
der Waals equation to give:
66
bRTv 3 2avv b
2
RTv 3 2av b
2
JT
RTv RTv 3 2av b 2
vreal
c Pideal
2av 3b
dv
videal v b 2
67
5.32
We can solve this problem by using the form of the Joule-Thomson coefficient given in Equation
5.75. The following approximation can be made
v v
T P T P
At 300 C,
h
c P h
dT P T P
At 300 C,
68
3 3
v 573.15 K 0.0005 m 0.25794 m
T
T v
kg K
kg
JT P
c P kJ
2.15
kg K
m3 K
JT 0.0133
kJ
69
5.33
At the inversion line, the Joule-Thomson coefficient is zero. From Equation 5.75:
v
T v
JT T P 0
Preal
v
2
c Pideal T 2 dP
T P
Pideal
v
T v 0
T P
RT a
P 2
v b v
Solving for T:
a v b
T P 2
v R
so
T a 1 v b 2a Pv 3 av 2ab
P 2
v P v R R v3 Rv 3
Substituting for P:
T RTv 3 2av b
2
v P v b Rv 3
Hence,
v bRTv 3 2avv b
2
T v 0
T P RTv 3 2av b
2
Solving for T:
2avv b
2
T (1)
bRv 3
70
Substituting this value of T back into the van der Waals equation gives
2avv b a a2v 3b
P 2 (2)
bv 3 v bv 2
We can solve Equations 1 and 2 by picking a value of v and solving for T and P. For N2, the
critical temperature and pressure are given by Tc = 126.2 [K] and Pc = 33.84 [bar], respectively.
Thus, we can find the van der Waals constants a and b:
2
RT Jm 3
a 27 c = 0.137 2
64 Pc mol
RTc m3
b = 3.88 10 -5
8 Pc mol
Using these values in Equations (1) and (2), we get the following plot:
1000
800
600
400
200
0
0 100 200 300 400
T [K]
71
5.34
We can solve this problem using departure functions, so first find the reduced temperatures and
pressures.
50 bar 10 bar
P1, r 0.99 P2,r 0.2
50.36 bar 50.36 bar
273.15 K
T1, r 0.967
282.4 K
Since the ethylene is in two-phase equilibrium when it leaves the throttling device, the
temperature is constrained. From the vapor-liquid dome in Figure 5.5:
T2,r 0.76
T2 214.6
dep dep
h hT , P hTideal
T2
hT , P 0
1, r 1, r 1 2,r 2,r
Therefore,
dep dep
hT , P hT , P hTideal
T2
2,r 2,r 1, r 1, r 1
1.424 14.394 10
214.6 K
3
hTideal R T 4.392 10 6 T 2 dT
1 T2
273.15 K
h dep
T1,r , P1,r 3.678 0.085 3.51 3.976
RT
c
72
hTdep, P 214.6 K
1 3 6 2
2,r 2,r
3.976 1.424 14.394 10 T 4.392 10 T dT
RTc 282.4 K
273.15 K
hTdep, P
2,r 2,r
3.01
RTc
We can calculate the quality of the water using the following relation
hTdep,, liq
4.6 0.085 5.5 5.068
2,r P 2,r
RTc
hTdep,, vap
0.4 0.0859 0.75 0.464
2,r P 2,r
RTc
Thus,
x 0.447
73
5.35
Density is calculated from molar volume as follows:
MW
v
2
Substitute the above into the expression for Vsound :
Vsound
2 P 1 P
MW MW 1 / v s
v s
1 v
v v2
Therefore,
2 P v 2 P
Vsound
s MW v s
74
5.36
From Problem 5.35:
2 P v 2 P
Vsound
s MW v s
P P s T P T s P s
v s s v v P v s T s s v v T s T T P
P T s P c P c P T
P
v s cv v T s T T cv T v v P
cP 7 P R T P
cv 5 T v v v P R
Therefore,
P 7 P
v s 5 v
and
v 2 7 P 7 RT
Vsound
MW 5 v 5 MW
75
5.37
From Problem 5.35:
2 P v 2 P
Vsound
s MW v s
P P s T P T s P s
v s s v v P v s T s s v v T s T T P
so
P T P T c P
v s cv T v v P T
P P T
v s T v v P
P T v
1
T v v P P T
So
P P
v s v T
m3
P = 2.34 kPa and v .001002
kg
m3
P = 5 MPa and v . 0009995
kg
76
So
and
P
Vsound v 2 1414 m/s
v s
77
5.38
(a)
The fundamental property relation for internal energy is
dU Qrev Wrev
dU TdS Fdz
dA dU d TS dU TdS SdT
dA Fdz SdT
(b)
First, relate the entropy differential to temperature and length.
S S
dS dT dz
T z Z T
u TS Fz S
nc z n T
T z T z T z
Therefore,
S nc a
z n b
T z T T
The following statement is true mathematically (order of differentiation does not matter):
A A
Z T z T T Z T z
Furthermore,
78
A S
Z T z T Z T
A F
T Z T z T z
S F
k z z0
Z T T z
Substituting the expressions for the partial derivatives into the expression for the entropy
differential, we obtain
a
dS n b dT k z z 0 dz
T
(c)
First, start with an expression for the internal energy differential:
U U
dU dT dz
T z Z T
U
na bT
T z
Using the expression for internal energy developed in Part (a) and information from Part (b)
U S
T F T k z z 0 kT z z 0 0
Z T z T
Therefore,
dU na bT dT 0dz na bT dT
(d)
We showed in Part (c) that
U
FU 0
z T
79
S
Using the expression for developed in Part B, we obtain
z T
S
FS T kT z z0
z T
(e)
First, perform an energy balance for the adiabatic process.
dU W
na bT dT Fdz kT z z 0 dz
Rearrangement gives
dT kT z z 0
dz na bT
The right-hand side of the above equation is always positive, so the temperature increases as the
rubber is stretched.
80
5.39
The second law states that for a process to be possible,
suniv 0
To see if this condition is satisfied, we must add the entropy change of the system to the entropy
change of the surroundings. For this isothermal process, the entropy change can be written
cv P P
ds dT dv dv
T T v T v
R
ds dv
v b
Integrating
v2 J
s sys R ln 11.5
v1 mol K
For the entropy change of the surroundings, we use the value of heat given in Example 5.2:
J
q q surr 600
mol
q surr 600 J
s surr 1.6
Tsurr 373 mol K
and
J
suniv s sys ssurr 9.9
mol K
Since the entropy change of the universe is positive we say this process is possible and that it is
irreversible.
Under these conditions propane exhibits attractive intermolecular forces (dispersion). The closer
they are together, on average, the lower the energy. That we need to put work into this system
says that the work needed to separate the propane molecules is greater than the work we get out
during the irreversible expansion.
81
5.40
A schematic of the process is given by:
Gas A in
h wS
Because the gas is not ideal under these conditions, we have to create a hypothetical path that
connects the initial and final states through three steps. One hypothetical path is shown below:
Pi ,Ti
100
ideal gas
P [bar]
step 1
Pf ,Tf
20
step 3
step 2
Plow
P 0
h
h1 dP
P T
P1
P 0
v
h1 Ti v dP
T P
Pi
82
v R aP
2
T P P Ti
Therefore,
P 0 v 0
2aP
dP 2467
aP J
h1 5 P v
dP
b
mol
Pi 10010 5 i
Ti T
Pi 10010
Pf 2010 5
2aP J
b dP 250
h3 Tf
mol
P 0
For step 2, the pressure is zero, so we can use the ideal heat capacity given in the problem
statement to calculate the enthalpy change.
Tf 445 K
J
h2 c P dT 30 0.02T dT 6270 mol
Ti 600 K
J
h h1 h2 h3 8487
mol
J
ws 8487
mol
In other words, for every mole of gas that flows through the turbine, 8487 joules of work are
produced.
83
Chapter 6 Solutions
Engineering and Chemical Thermodynamics
Wyatt Tenhaeff
Milo Koretsky
[email protected]
6.1
(a)
The Clausius-Clapeyron equation:
so
h vap 1 1
Pi sat
101 kPaexp i
R T 373 [K]
(b) and (c)
Using
kJ
h vap 40.626
mol
2
The logarithmic trend is well-represented. However, at lower temperatures the Clausius-
kJ
Clapeyron equation is up to 37% off. The actual heat of vaporization changes from 2501.3
kg
kJ
at 0.01 oC to 2257.0 at 100 oC, a difference of around 10%.
kg
100.00
10.00
Eqn 6.23
Tables
1.00
0.10
273 283 293 303 313 323 333 343 353 363 373
Temperature [K]
(d)
For 100 C to 200 C, we obtain the following table:
3
468.15 1445.10 1397.8 3.38%
473.15 1613.62 1553.8 3.85%
10000.00
Eqn 6.23
1000.00
Tables
100.00
373 383 393 403 413 423 433 443 453 463 473
Temperature [K]
Over this range the Clausius-Clapeyron equation represents the data well and is no more than 4
% off. The actual heat of vaporization changes from 2257.0 kJ/kg at 100 oC to 1940.7 kJ/kg
at 200 oC, a difference of around 15%.
(e)
The heat of vaporization can be corrected for temperature as follows
Tb T
hvap T
l
cP dT hvap Tb c Pv dT
T Tb
We can acquire heat capacity data from Appendix A.2, but to simplify the analysis, we will use
an average heat capacity for the vapor.
dPisat
56026 41.27T dT
Pisat RT 2
Integrate:
4
1 1 1 T
Pisat 101.35 kPa exp 56026 41.27 ln
R T 373.15 373.15
5
468.15 1283.56 1397.8 8.17
473.15 1417.67 1553.8 8.76
10000.00
Equation Steam Tables
1000.00
100.00
Psat (kPa)
10.00
1.00
250 300 350 400 450 500
0.10
T (K)
The agreement between the two values at lower temperatures improves significantly at lower
temperatures, but actually worsens at higher temperatures. The agreement could potentially be
improved by not averaging the heat capacity.
6
6.2
We can find the required pressure by applying the Clapeyron equation:
dP hl h s
dT
vl v s T
We can find the molar volume of water ice from any number of reference books. At 0 C and 1
bar:
kg 5
m3
1000
l l
v 1.80 10
m3 mol
kg m3
s 917 v s 1.97 10 5
m3 mol
J
h l h s 6010
mol
If we assume that h l h s and v l v s are independent of temperature and pressure, we can
separate variables in the Clapeyron equation and integrate.
hl h sT
ln 2
P2 P1
v vl s
T1
6010 J/mol 268.15 K
P2 1 10 5 Pa
1.80 10 5
1.97 10 5
3
ln
m / mol 273.15 K
so
P2 66.1 bar
7
6.3
(a)
At 1 bar, the gas will act as an ideal gas.
J
8.314 300 K
RT mol K m3
v 0.0249
P 1 10 5 Pa mol
The number of moles of vapor are found as follows (neglect molar volume of liquid)
vV
n
0.001 m 3
v m3
0.0249
mol
n v 0.0402 mol
(b)
At 21 bar, the gas will not behave ideally. Since we are assuming that the molar volume of
liquid is negligible and the heat of vaporization is independent of temperature, the Clapeyron
equation becomes
dP h vap
v
dT vT
vv
RT 1 B ' P
RT
1
B'
P P
dP h vap
dT 1
R B ' T 2
P
P2 21 105 Par
h vap
T2
1 ' dT
B dP
P R
T1 300 K T
2
P1 1 105 Pa
8
P h vap 1 1
ln 2 B ' P2 P1
P1 R T2 T1
We can substitute values for given quantities and constants to solve for T2.
T2 523.3 K
(c)
Using the virial equation,
1 J 7 m
3
v v RT B ' 8.314
523.3 K
1
1 10
P mol K 21 10 Pa
5
J
m3
v v 0.00164
mol
nv
V
0.001 m 3
vv m3
0.00164
mol
n v 0.61 mol
9
6.4
We can use the following computational path to solve for pressure at which graphite and
diamond are in equilibrium at 25 oC.
graphite diamond
P = 1 [atm]
J
g 2866
mol
P 1
g1 v graphdP g 3 v diam dP
1 P
graphite g 2 0 diamond
J
g 1[atm] 2866 g1 g 2 g 3
mol
To find the change in Gibbs energy with pressure, we apply the fundamental property relation,
Equation 5.9. At constant temperature:
0
dgi v idP sidT
g i vi dP vi P
Solving
J 1
2866
1 cm 3 g 1 m 3
12
mol 2.26 3.51 g mol 10 cm
6
5
3 P 1.01 10 Pa
or
10
6.5
From the Clausius-Clapeyron equation:
dP h lAl h Al
s hTfus
(I)
dT v lAl v Al
s
T v lAl v Al
s
T
where h Tfus is the enthalpy of fusion at temperature T. We can get the molar volumes from the
densities:
kg
0.027
v lAl
MW
mol 1.17 10 5 m3
l kg mol
2,300
m3
and
kg
0.027
s
MW
mol 1.00 10 5 m3
v Al
s kg mol
2,700
m3
so
m3
v lAl v Al
s
1.7 10 6
mol
solid liquid
T
h Tfus
933 T
cP dT
s
cPl dT
T 933
T = 933 [K]
solid J liquid
h fus 10,711
mol
11
933.45 T
J
hTfus
s l
c P dT 10,711 cP dT
T 933.45
mol
J J
Using c Ps 20.608 0.0138T l
and c P 31.748 , we get:
mol K mol K
Integrating:
or
T = 934.91 [K]
12
6.6
We can assume that silver acts as an ideal gas at 1500 K. We can also assume the molar volume
of the vapor is much greater than the molar volume of liquid. Therefore, we can use Equation
6.22
dP sat h vap dT
P sat RT 2
We can differentiate the expression for pressure in the problem statement to obtain
Therefore,
14260 0.458 h
vap
T T RT 2
2
Therefore,
At 1500 K,
J
h vap 8.314 14260 0.4581500 K
mol K
kJ
h vap 112.8
mol
13
6.7
For a single component system:
Gi g i g
dg sdT vdP
We can identify a phase transition from the vertical line of the g vs. T plot, as indicated below.
Since this transition is vertical, i.e., the temperature is constant, the pressure must also be
constant. Thus, we can differentiate the Gibbs energy with respect to temperature at constant
pressure to get:
g
s
T P
Hence the slope of a plot of g (or ) vs. T at any temperature must be the negative of the value of
entropy on the plot for s vs. T. The resulting curve is sketched below.
phase transition;
vertical line indicates
P = const
14
6.8
The ferrite phase has stronger bonds. At room temperature, iron is in the ferrite phase. The
heating to 912 C has the effect of increasing the entropy contribution to the Gibbs energy. At a
high enough temperature, the austenite phase becomes stable, so that its entropy must be greater
than the ferrite phase. If the entropy of the austenite phase is greater, the enthalpy of the ferrite
phase must be greater or else the austenite phase would be stable over the entire temperature
range. Hence, the ferrite phase has stronger bonds.
15
6.9
Since the pressures are low, we can assume ideal gas behavior. We can also assume that the
molar volume of the vapor is much greater than the molar volume of liquid and the heat of
vaporization is independent of temperature. Therefore, we can rearrange Equation 6.24 to obtain
J 760 torr
8.314 ln
mol K 400 torr
h vap
1 1
353.25 K 333.75 K
kJ
h vap 32.3
mol
16
6.10
(a)
The freezing point occurs where there is a discontinuity in the g vs. T plot, as indicated below.
The liquid is at a temperature higher than the freezing point and the solid at lower temperature.
These are demarked below. The melting temperature is 250 K, which occurs at a value g = 3,000
[J/mol]
New Freezing
point
4,000
solid
Freezing
point
2,000 1.2
units
liquid up
Temperature, T (K)
(b)
At constant pressure, the entropy can be found from Equation 5.14. For the solid we have:
g g 1,000 J
s 10
T P T 10 mol K
g g 2,000 J
s 40
T P T 50 mol K
(c)
As we change pressure, we can see how the Gibbs energy changes at any given temperature by
Equation 5.14:
17
g
v
P T
Assuming the molar volumes of the liquid and vapor stay constant over the temperature range
around the melting point, we see that the Gibbs energy of the liquid increases by 1.2 times the
Gibbs energy of the solid, since the molar volume of the liquid is 20% larger. The Gibbs energy
of the new freezing point at higher pressure is schematically drawn on the plot above. For
convenience, we choose the solid to increase by 1 unit on the plot. Thus, the liquid increases by
1.2 units. As the sketch shows, the freezing point, where the two lines intersect, will shift to
higher temperature.
18
6.11
For a single component system the fundamental property relation, Equation 5.9, gives:
dg sdT vdP
We can identify a phase transition from the vertical line of the g vs. P plot, as indicated below.
Since this transition is vertical, i.e., the pressure is constant, the temperature must also be
constant. Thus, we can differentiate the Gibbs energy with respect to pressure at constant
temperature to get:
g
v
P T
Hence the slope of a plot of g vs. P must have a slope that matches the plot for v vs. T. Since the
molar volume of phase is about twice the value of phase , its slope should be twice as big.
The resulting curve is sketched below.
phase transition;
v vertical line indicates
v T = const
P* P
19
6.12
The saturation pressure can be found using the Clausius-Clapeyron equation with the assumption
that the heat of vaporization is independent of temperature. First, we need to use the given data
for the 63.5 C and 78.4 C to find the heat of vaporization.
P sat
R ln 2 J 760 torr
P sat 8.314 ln
vap 1 mol K 400 torr
h
1 1 1 1
351.55 K 336.65 K
T2 T1
kJ
h vap 42.39
mol
kJ
42.39
h vap 1 mol
760 torr exp
1 1 1
P3 P2 exp
sat sat
R T3 T2 kJ 373.15 K 351.55 K
0.008314
mol K
P3sat 1760 torr 2.32 atm
In comparison, ThermoSolver gives a value of 2.23 atm, using the Antoine equation.
20
6.13
We can show using the Chain Rule that
gi Tg i g i T
T T2 1 g i g i
T T T T P T 2
P P
g i
si
T P
Therefore,
gi
T Tsi g i Tsi hi Tsi hi
T T2 T2 T2
P
21
6.14
Let T1 922 K, T2 1,300 K
g 2 h2 T2 s 2
dh c P dT
1,300 K
h2 h1 c P dT
922 K
J
h2 h1 c P T2 T1 39,116
mol
c
ds P dT
T
T J
s2 s1 c P ln 2 85.10
T1 mol K
J
g 2 h2 T2 s2 71,500
mol
d g
T hT h1 c P T T1 We must leave h as a function of T
dT T2 T2
g
T 2
h c T T
1300 K
g
d T 1 PT 2 1 dT
g 922 K
T 1
g g 1 1 T J
h1 c PT1 c P ln 2 55.01
T 2 T 1 T2 T1 T1 mol K
J
g 2 71,500
mol
22
6.15
A possible hypothetical solution path is presented below:
monoclinic orthorhombic
T = 298 [K]
g g
1 3
T = 368.3 [K]
monoclinic g = 0 orthorhombic
2
From the diagram, we see that the Gibbs energy for steps one and three can be calculated as
follows:
g m
368 K
g1 dT
298 K
T P
and
g o
298 K
g 3 T dT
368 K P
respectively. We can apply Equation 5.14 from the thermodynamic web
g
s
T P
g 3m68o
K 0
mo
Using these above relationships, the expression for g 298 K becomes
Therefore, the transition from the monoclinic to orthorhombic state occurs spontaneously. The
orthorhombic state is more stable.
23
6.16
At the phase transition, the following is true
g g
T Sr ( s ) T Sr (l )
Using the thermodynamic web, the following can be shown (see Problem 6.13)
g / T h
2
T P T
T
h l T 49179 35.146 dT 35.146T 3540
1500 K
T
h s T 20285 37.656 dT 37.656T 16305.4
900 K
g
can also be calculated at 900 K for solid Sr and 1500 K for liquid Sr.
T
gl J
83.85
T
ref mol K
gs J
68.68
T
ref mol K
g
We can find at any temperature using the differential equation as follows
T
g g / T h
d T T P
dT 2 dT
T
g l /T T
g 35.146T 3540
d T T 2 dT
83.85 1500 K
24
gl
35.146 ln(T ) 3540 176.04
T T
g s /T T
g 37.656T 16305.4
d T T 2 dT
68.68 900 K
gs
37.646 ln(T ) 16305.4 205.5
T T
gl gl
Set and solve for T:
T T
T melt 1059.8 K
h fus h s T melt h l T melt
Using the expressions developed above
kJ kJ kJ
h fus 26.30 33.71 7.41
mol mol mol
25
6.17
At the phase transition, the temperature and Gibbs energy of both phases must be equal.
Mathematically, this is equivalent to
g g
T SiO2 ( s ) T SiO 2 (l )
Using the thermodynamic web, the following can be shown (see Problem 6.13)
g / T h
2
T P T
hT T
dh c P dT
href Tref
T
hT href c P dT
Tref
T
h T 738440
l
85.772 dT
2500 K
g
can also be calculated at 1100 K for solid SiO2 and 2500 K for liquid SiO2.
T
gl J
487.3
T
ref mol K
gs J
903.5
T
ref mol K
26
and hT into the above differential equation and
g
We can substitute our expressions for
T
separate variables to obtain
T
l
g /T T
738440 85.772 dT
g
T 2500 K
d 2
dT
487.3 2500 K T
T
5 2 9 3
s
g /T T
856840 53.466 0.02706T 1.27 10 T 2.19 10 T dT
g
T 1100 K
d 2
dT
903.5 1100 K T
Integration provides
gl gs
If we plot vs. T, we obtain the following:
T T
27
There are three solutions, but only the solution between 1100 K and 2500 K is physically
meaningful. If we magnify the plot near the middle solution, we find
T 1983 K
The enthalpy of fusion is defined as
h fus h s T melt h l T melt
Using the expressions developed above
kJ kJ kJ
h fus 792.5 782.78 9.72
mol mol mol
28
6.18
From the Clausius-Clapeyron equation
sat
dPCS h vap
2
dT T vv vl
Assuming:
v v v l
we get
sat
dPCS h vap
2
(I)
dT Tv v
4.7063 10 3
sat
ln PCS 62.7839 6.7794 ln T 8.0194 10 3 T (II)
2 T
sat sat
d ln PCS 1 dPCS 4.7063 10 3 6.7794
2
2
8.0194 10 3 (III)
sat
dT PCS dT T2 T
2
vap kJ
Solving for vv using hCS 24.050 gives:
2
mol
1
h vap 4.7063 10 3 6.7794 m3
v
v 8.0194 10 3 6.08 10 3
sat
TPCS T2 T mol
2
Pv B
z 0.878 1
RT v
29
or
m3 cm3
B 7.4 10 4 740
mol mol
Alternative solution:
sat
dPCS h vap
2
dT
1
RT 2 B '
sat
PCS 2
1 h vap
sat B' dPCS
sat
2 dT
P 2 RT
CS2
We must be careful about the limits of integration. We need to pick a value of T close so
enthalpy of vaporization is not too different, but far enough away to avoid round off error. If we
5.04 10 Pa. Integrating:
sat 5
choose T = 378 K, Equation I gives PCS
2
4.4810 5 Pa 373
h vap
sat B dPCS
1 '
sat
2 2 dT
P
5.04 10 Pa CS2
5 378 RT
or
m 3
B B ' RT 7.9 10 4
mol
30
6.19
Calculate vA, vB, v, VA, VB, and V from the ideal gas law:
RT
vA 0.05 m 3 / mol V A n A v A 0.1 m 3
P
RT
vB 0.05 m 3 / mol V B n B v B 0.15 m 3
P
RT
v 0.05 m 3 / mol V ntot v 0.25 m 3
P
V
V A n A n B RT RT 0.05 m 3 / mol
n A T , P, n B n A P P
V
V B n A n B RT RT 0.05 m 3 / mol
n B T , P, n A n B P P
To find the remaining quantities, we can apply Equations 6.44 and 6.46
Vmix n A V A v A n B V B v B
Vmix 20.05 0.05 30.05 0.05 0
vmix x A V A v A x B V B v B
vmix 0
31
6.20
(a)
For a pure species property
va v y a 1
Substitution yields
cm 3
va 1001 800 2.510 100
mol
(b)
From Equation 6.29
V
Va
na n b ,T , P
We can find V by multiplying the given expression for molar volume by the total number of
moles.
y a yb
V na nb 100 ya 80 yb 2.5
n n
100na 80nb 2.5 a b
y a yb na nb
n n nb na nb
Va 100na 80nb 2.5 a b 100 2.5 2.5
na na nb n
b
na nb na nb 2
so
To find the molar volume at infinite dilution, we can use the following relation
Va lim Va
ya 0
cm 3
Va 102.5
mol
32
(c)
Since species A contributes more to a mixture than to a pure species,
v mix 0
Note: The Gibbs-Duhem equation says that species B also contributes more.
33
6.21
Calculate mole fractions:
n1 1 mol
y1 0.2 y 2 0.4 y3 0.4
ntot 5 mol
Calculate v.
Obtain an expression for v:
2 A B
1 P RT y1 y 2 RT
RT
v
P
Substitute values:
3
82.06 cm atm 500 K
mol K
v
50 atm
1 50 9.0 10
2 5
0.2 0.4 3.0 10 5
cm 3
v 919.0
mol
Calculate V.
V ntot v
cm 3
V 5 mol 919.0
mol
V 4595 cm 3
Calculate v1.
The value of v1 can be found by substituting y1=1 into the expression for v1.
3
82.06 cm atm 500 K
mol K
v1
50 atm
1 50 9.0 10
2 5
1 0 3.0 10 5
cm 3
v1 698
mol
34
Calculate v2.
3
82.06 cm atm 500 K
v2 mol K
50 atm
1 50 9.0 10
2 5
0 1 3.0 10 5
cm 3
v2 1067
mol
Calculate v3.
3
82.06 cm atm 500 K
mol K
v3
50 atm
1 50 3.0 10
2 5
cm 3
v3 882
mol
Calculate V1 .
From Equation 6.29:
V nv
V1
n1 n 2 , n 3 ,T , P n1 n 2 , n 3 ,T , P
We can substitute the expression for V into this derivative and use the fact that
ntot n1 n2 n3 to obtain
RT 2 A
V1 n n n P RT n n
B
n n n
n1 P
1 2 3 1 2 1 2 3
RT n 2 , n 3 ,T , P
Differentiating we get
RT 2 A B
V1 1 P RT RT
P
Substitute values:
35
3
82.06 cm atm 500 K
mol K
V1
50 atm
1 50 9 10
2 5
3 10 5
cm 3
V1 697.5
mol
36
6.22
(a)
By definition:
H
H a
na T , P , n b , n c
n na nb nc
na nbnc
H nh 5,000na 3,000nb 2,200nc 500
na nb nc 2
H
Ha 5,000 500xb x c 1 2xa J/mol
na T ,P,n ,n
b c
(b)
1
xa xb xc
3
Ha 5,018.5 J/mol
(c)
xa 1 , xb xc 0
Ha 5,000 J/mol
(d)
xb 1 , xa xc 0
Hb hb 3,000 J/mol
37
6.23
Let the subscript 1 designate CO2, and 2 designate propane. To calculate the partial molar
volumes, the following formulas will be used:
dv
V1 v y 2
dy 2
v y1V1 y 2V2
Expressions cant be obtained for the molar volume with the van der Waals EOS; therefore, the
problem will be solved graphically. First, obtain an expression for the pressure that contains the
mole fractions of CO2 and propane:
RT y12 a1 2 y1 y 2 a1a2 y 22 a2
P
v y1b1 y 2 b2 v2
J m3 mol
a1 0.366 b1 4.29 10 5
mol 2 m3
J m3 mol
a 2 0.941 b2 9.06 10 5
2
mol m3
y1 y2 amix bmix v
The last column contains the molar volumes obtained by solving the van der Waals equation
with the spreadsheets solver function. After the table is completed, we create the following
graph.
38
v vs. y 2
1.48E-03
1.46E-03
1.44E-03
1.42E-03
v (m 3 mol -1)
1.40E-03
1.38E-03
1.36E-03
y = -0.00008x 2 - 0.00008x + 0.00147
1.34E-03
R2 = 0.99946
1.32E-03
1.30E-03
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
y2 (Mole Fraction)
v 8 10 5 y 22 8 10 5 y 2 0.00147
Therefore,
dv
1.6 10 4 y 2 8 10 5
dy 2
and
V1 v y 2 1.6 10 4 y 2 8 10 5
We can find the partial molar volume of propane from the following relationship
v y1V1 y 2V2
v y1V1
V2
y2
Tabulate the values of the partial molar volumes in the spreadsheet and create the following
graph
39
Partial Molar Volumes as a Function of Carbon Dioxide Mole
Fraction
0.0016
V1
V2
0.00155
mol -1)
0.0015
3
Partial Molar Volume (m
0.00145
0.0014
0.00135
0.0013
0.00125
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
y1 (CO 2 mole fraction)
40
6.24
(a)
kJ
ga 40
mol
dg
Ga g x b
dx b
dg
40 60 RT ln xa 1 ln x b 1 5x a 5xb
dx b
Ga 40xa 60x b RT xa ln x a x b ln xb 5xa x b
x b40 60 RT ln xa 1 ln x b 1 5x a 5xb
Ga 40xa xb RT xa xb ln x a 5xb2
Ga 40 RT ln x a 5 xb2
8.314(300) kJ kJ
Ga 40 ln 0.2 5(0.64) 40.8
1000 mol mol
Ga
Gmix ng x a ga xb gb
g xa ga x bgb RT x a ln xa xb ln x b 5 x a xb
(b)
hmix 5x a xb 0
so
h hmix hsensibleheat 0
so
hsensibleheat 0 and T goes down
41
6.25
To find V1 and V2 , we can read values directly from the graphs. Calculate mole fractions
1
x1 0.2
5
At x1 0.2 ,
cm 3
V1 46.5
mol
cm 3
V2 69.8
mol
The following relationships are employed to calculate the molar volumes of pure species
v1 lim V1
x1 1
v 2 lim V2 lim V2
x 2 1 x1 0
cm 3
v1 50
mol
cm 3
v2 70
mol
Therefore,
V1 50 cm 3
V2 280 cm 3
V n1V1 n2V2
V 146.5 469.8 325.7 cm3
42
Therefore
V 325.7
v
ntot 1 4
cm 3
v 65.14
mol
Vmix 4.3 cm 3
43
6.26
(a)
Expression for hmix :
Therefore,
H mix
H mix Cd H Cd hCd
nCd n Sn ,T , P
(b)
We can show by repeating Part (a) for Sn that
H mix Sn H Sn hSn
Equation 6.65:
dhmix
H mix Cd hmix X Sn
dX Sn
Since,
hmix 13000 X Cd X Sn
We get,
dhmix
d
13000 X Cd X Sn 13000 X Cd X Sn
dX Sn dX Sn
dhmix
d
13000 X Cd X Sn 13000 X Sn X Cd
dX Cd dX Cd
44
H Cd hCd H mix Cd 13000 X Cd X Sn X Sn13000 X Cd X Sn 13000 X Sn
2
J
H Cd hCd H mix Cd 2080
mol
and
J
H Sn hSn H mix Sn 13000 X Cd
2
4680
mol
(c)
Equation 6.37:
d H mix Cd d H mix Sn
nCd n Sn 0
dxCd dxCd
where
d H mix Cd
26000 X Sn 26000 26000 X Cd
dX Cd
d H mix Sn
26000 X Cd
dX Sn
Therefore,
d H mix Cd d H mix Sn
nCd nSn xCd ntot 26000 26000 xCd 1 xCd ntot 26000 xCd
dxCd dxCd
d H mix Cd d H mix Sn
nCd n Sn 0
dxCd dxCd
(d)
A graphical solution can be found using the tangent-slope method discussed on pages 285-287:
A plot of a line tangent to the enthalpy of mixing curve at XCd = 0.6, is given below:
45
Heat of mixin g in cadmium (Cd )-Tin (Sn) s ystem
6000
5400
4800
4200
mol
3600
J
data
3000
hmix
2400
1800
1200
600
0
0 0 .2 0 .4 0 .6 0 .8 1
fit to: xCd
J
hmix 13, 000 xCd x Sn
mol
J
H Cd hCd 2050
mol
J
H Sn hSn 4800
mol
The values using the graphical method are reasonably close to the analytical method.
46
6.27
The following can be shown with the Gibbs-Duhem equation
0 x1V1 x 2V2
dV1 dV
0 x1 x2 2
dx1 dx1
If the partial molar volume of species 1 is constant, the Gibbs-Duhem equation simplifies to
dV2
0
dx1
Note that in this case, since the partial molar volume of species 1 is constant:
V1 v1
` V2 v2
47
6.28
(a)
Let species 1 represent HCl and species 2 represent H2O. An expression for the enthalpy of the
solution is
h x1h1 x 2 h2 hmix
Therefore,
~
H n1h1 n2 h2 n1hs
To use the heat of solution data in Table 6.1, we need to determine the values of n1 and n2
consistent with the convention used in the table. As seen in Example 6.6,
n1 1
x1
n1 n2 1 n
For this problem
x1 0.2
Therefore,
n1 1
n2 n 4
H
H H 2 O H 2
2 n1 , T , P
n
~
dhs
H H 2 O H 2 h2 n1
dn2
~
dhs
H H 2 O hH 2 O H 2 h2 n1
dn2
48
~ J
hs 61,204
mol solute
56852
J J
64049
dhs hs n 5 hs n 3
~ ~ ~
mol solute mol solute
dn2 5 mol 3 mol 5 mol 3 mol
~
dhs J
3598.5
dn2 mol mol solute
Therefore,
J
H H 2 O hH 2 O 1 mol solute 3598.5
J
3599
mol mol solute mol
H HCl hHCl
hmix x H 2 O H H 2 O hH 2 O xHCl h~s xH O H H O hH O
2 2 2
x HCl x HCl
J J
0.2 61204 0.8 3599
mol mol J
H HCl hHCl 46808
0.2 mol
(b)
For n1 2 and n2 80 ,
n n2 82
n 1 1 1 40
n1 2
The new values for the number of moles consistent with Table 6.1
n1 1 mol
n2 40 mol
dhs hs n 50 hs n 30
~ ~ ~
dn2 50 mol 30 mol
Interpolating the data in Table 6.1
49
hs n 30 72428
~ J
mol solute
Therefore,
72428
J J
~ 73729
dhs
mol solute mol solute 65.05 J
and
dn2 50 mol H 2O 30 mol H 2O
mol mol solute
H H 2 O hH 2 O 1 mol solute 65.05
J J
65.05
mol mol solute mol H 2 O
50
6.29
First perform an energy balance on the mixing process.
hmix q
We can calculate hmix using data from Table 6.1. Referring to Equation E6.7A, we find
~
hmix x HCl hs
Calculate x HCl :
wHCl 0.30
MW HCl 36.46
x HCl 0.175
wHCl w 0 . 30 0.70
H 2O
MW HCl MW H 2O 36.46 18.0148
Heats of data are tabulated for a solution containing one mole of the solute for various amounts
of water. Thus, we need to calculate how many moles of water must be added to HCl to obtain
the above mole fraction.
1 mol HCl
x HCl ; where n is the number of moles of H2O
1 mol HCl n
n 4.71 mol H 2 O
~ J
hs 63224 (for n 4.71 )
mol
Therefore,
J J
hmix 0.175 63224 11064
mol mol
and
J
q hmix 11064
mol
51
6.30
To calculate the enthalpy of mixing from Table 6.1, we must use the following expression
~
hmix x H 2 SO 4 hs
1
x H 2 SO4
1 n
hmix 74.4 x H 2 SO 4 x H 2 O 1 0.561x H 2 SO
4
For n 1 , x H 2 SO 4 0.5 and x H 2 O 0.5
~ J
Table 6.1: hs 31087
mol
J J
hmix 0.5 31087 15543.5
mol mol
52
As you can see, the percent difference between the two methods decreases as the mole fraction
of sulfuric acid decreases. Although Equation 6.47 fit data at 21 C, while the Table 6.1
tabulates data taken at 25 C, we do not expect the temperature dependence to account for all the
observed difference. The table and equation come from different experimental data sets, and
also represent measurement uncertainty. Nevertheless, the agreement is reasonable.
53
6.31
To calculate the enthalpy of mixing from Table 6.1, we must use the following expression
~
hmix x HCl hs
1
x HCl
1 n
where n is the number of moles of water. The following table was made using these two
equations.
~ J J
n [mol H2O] x1 hs hmix
mol HCl mol
1 0.5 -26225 -13112.5
2 0.333 -48819 -16273
3 0.25 -56852 -14213
4 0.2 -61204 -12240.8
5 0.167 -64049 -10674.8
10 0.091 -69488 -6317.09
20 0.048 -71777 -3417.95
50 0.020 -73729 -1445.67
100 0.0099 -73848 -731.168
54
6.32
A schematic for the process is given below. The inlet streams are labeled 1 and 2 and the
exit stream 3.
q
Stream 1 Stream 3
Stream 2 H 2O
Q H 3 H 2 H1
wNaOH 0.50
MW NaOH 40
x NaOH ,1 0.311
wNaOH w H 2O 0.50 0.50
MW NaOH MW H 2O 40 18.0148
wNaOH 0.10
MW NaOH 40
x NaOH ,3 0.048
wNaOH wH 2 O 0.10 0.90
MW NaOH MW H 2O 40 18.0148
We now calculate the moles of water per mole of NaOH so that we can use Table 6.1:
1
x NaOH
1 nH 2 O
n H 2 O,1 2.21
55
n H 2 O,3 19.8
Since enthalpy is a state function, we can choose any hypothetical path to calculate the change in
enthalpy. One such path is shown below. The box in our original schematic is depicted with
dashed lines below. We pick a basis of 1 mole NaOH. In step A, the inlet stream is separated
into its pure components. In step B, 17.6 additional moles of water are added to the pure water
stream. Finally the H2O and NaOH streams are remixed
Step A Step C
1 mol NaOH 1 mol NaOH 1 mol NaOH
2.2 mol H2O 19.8 mol H2O
Step B
hmix=0
2.2 mol H2O 19.8 mol H2O
H 3 H 2 H 1 H A H B H C
The enthalpies of mixing for steps A and C can be related to enthalpy of solution data from Table
6.1:
~ J
hs ,1 23906
mol NaOH
~ J
hs ,3 42858
mol NaOH
For step A, we need the negative value of the heat of solution of stream 1. Thus for a basis of 1
mole NaOH:
H A 23906 J
56
while for step C:
H C 42858 J
To get the total heat that must be removed per mole of product solution, we divide by the number
of moles of product per mol of NaOH:
H J
q 910
n NaOH mol
57
6.33
The partial molar property can be written as follows:
n K
K1 T
n1 T , P ,n2 ,n3
n k
K1 k T nT
n1 T , P ,n2 ,n3 n1 T , P ,n2 ,n3
k
K1 k nT (1)
n1 T , P ,n2 ,n3
k
Now focus on nT . At constant T and P, we can write,
1 T , P ,n2 ,n3
n
k k k
dk dx1 dx2 dx3
x1 T , P, x2 , x3 x2 T , P, x1 , x3 x3 T , P, x1 , x2
Therefore,
k k x1 k x 2
n1 T , P ,n2 ,n3 x1 T , P , x2 , x3 n1 T , P ,n2 ,n3 x2 T , P , x1 , x3 n1 T , P ,n2 ,n3
(2)
k x3
x3 T , P , x1 , x2 n1 T , P ,n2 ,n3
but,
n1
x1
n1 n2 n3
so
x1
1 x1
1 n1 1
(3)
n1 T , P ,n2 ,n3 n1 n2 n3 n1 n2 n3
2
nT
Similarly,
58
x 2 n2 x
2 (4)
n1 T , P ,n2 ,n3 n1 n2 n3 2
nT
and
x3 n3 x
3 (5)
n1 T , P ,n2 ,n3 n1 n2 n3 2
nT
k
K1 k 1 x1 k x2 k x3
x1 T , P, x 2 , x 3 x 2 T , P, x1 , x 3 x3 T , P, x1 , x 2
k k
dk dx1 dx 2
x1 T , P, x2 , x3 x 2 T , P, x1 , x3
k k dx1 k dx2
x 2 T , P , x3 x1 T , P , x2 , x3 dx2 x2 T , P , x1 , x3 dx2
Thus,
k k k
(7)
x2 T ,P , x3 x1 T ,P , x2 , x3 x2 T ,P , x1 , x3
k k k
(8)
x3 T , P , x2 x1 T , P , x2 , x3 x3 T , P , x1 , x2
59
k k
K1 k x 2 x3
x2 T , P , x3 x3 T , P , x2
k
K i k xm
m i xm T , P , x ji ,m
60
6.34
The expression can be found by employing the Gibbs-Duhem equation:
0 n1V1 n2V2
Differentiate with respect to x1 and then divide by the total number of moles:
dV1 dV dV dV
0 x1 x 2 2 x1 1 1 x1 2
dx1 dx1 dx1 dx1
dV1
5.28 5.28 x1
dx1
Substitution of this result into the Gibbs-Duhem equation and rearrangement yields
Integrate:
cm 3
V2 2.64 x12 C
mol
To determine C, we can use the density information given in the problem statement.
V2 x 2 1 V2 x1 0 C v 2
where
1 1 cm 3
C v2 109.58
2 0.768 g/cm 3
mol
MW2 84.16 g/mol
Therefore,
cm 3
V2 2.64 x12 109.58
mol
61
6.35
An expression for the enthalpy of the solution is
h x1h1 x 2 h2 hmix
which is equivalent to
~
h x1h1 x 2 h2 x1hs
To use the heat of solution data in Table 6.1, we need to determine the values of n1 and n2
consistent with the convention used in the table. As seen in Example 6.6,
n1 1
x1
n1 n2 1 n
For this problem
x1 0.33
Therefore,
n1 1
n2 n 2
H
H H 2 O H 2
n2 n1 ,T , P
~
dhs
H H 2 O H 2 h2 n1
dn2
~
dhs
H H 2 O hH 2 O H 2 h2 n1
dn2
62
812
J J
2,787
dhs hs n 3 hs n 1
~ ~ ~
mol solute mol solute
dn2 3 mol 1 mol 3 mol H 2O 1 mol H 2O
~
dhs J
987.5
dn2 mol mol solute
Therefore,
J
H H 2 O hH 2 O 1 mol solute 987.5
J
987.5
mol mol solute mol
Calculate the partial molar enthalpy:
J
H H 2O hH 2 O 987.5
mol
kJ
hH 2 O 104.87
kg
kJ
hH 2 O 1.89
mol
kJ J J
H H 2 O 1.89 0.988 0.90
mol H 2 O mol H 2 O mol H 2 O
63
6.36
(a)
Calculate the mole fraction of sulfuric acid
w1 0.20
x1
MW 1 98.078 0.044
w1 w2 0.20 0.80
MW 1 MW 2 98.078 18.0148
1
x1 0.044
1 n
n 21.7 mol H 2 O
~ J
hs 74621
mol
J
q hmix x1hs 0.044 74621
~
mol
J
q 3283
mol
(b)
Calculate the mole fraction of pure sulfuric acid. Consider a mixture of 20 kg of 18 M sulfuric
acid and 80 kg of water. Find the mass of sulfuric acid present.
mH 2 SO4 20 kg
VH 2 SO4 10.9 L
H 2 SO4 1.84 kg/L
nH 2 SO4 VM 10.9 L 18 mol/L 196.2 mol
m H 2 SO 4 n H 2 SO4 MW H 2 SO 4 196.2 mol0.098 kg/mol 19.2 kg
Since both the initial (i) and final (f) states contain mixtures, to get the enthalpy of mixing, we
need to calculate the relative differences follows:
64
~ ~
q hmix x1, f hs, f x1,i hs ,i
w1 0.192
MW 1 98.078
x1, f 0.042
w1 w2 0.192 0.808
MW 1 MW 2 98.078 18.0148
Calculate n to use in Table 6.1:
1
x1, f 0.042
1 n
n 22.8 mol H 2 O
~ J
hs , f 74689
mol
w1 0.192
MW 1 98.078
x1,i 0.81
w1 w2 0.192 0.008
MW 1 MW 2 98.078 18.0148
Calculate n to use in Table 6.1:
1
x1,i 0.81
1 n
n 0.23 mol H 2 O
We must extrapolate from Table 6.1. To do this we wish to extend the trend at low water
concentration. A plot of the data in Table 6.1 is useful. A semi-log plot follows:
65
0
3200 [J/mol]
-10000
-30000
-40000
-50000
-60000
-70000
-80000
-90000
0.1 1 10 100
n
~ J
hs ,i 3200
mol
J J
q hmix x1, f hs , f x1,i hs,i 0.042 74689 0.81 3200
~ ~
mol mol
J
q 710
mol
(c)
Calculate the mole fraction of sodium hydroxide
w1 0.20
x1
MW 1 40 0.101
w1 w2 0.20 0.80
MW 1 MW 2 40 18.0148
1
x1 0.101
1 n
n 8.9 mol H 2 O
~ J
hs 41458
mol
66
Now calculate the heat transfer
J
q hmix x1hs 0.101 41458
~
mol
J
q 4187
mol
(d)
Calculate the mole fraction of ammonia
w1 0.20
x1
MW 1 17.03 0.209
w1 w2 0.20 0.80
MW 1 MW 2 17.03 18.0148
1
x1 0.209
1 n
n 3.78 mol H 2 O
~ J
hs 33153
mol
J
q hmix x1hs 0.209 33153
~
mol
J
q 6929
mol
67
6.37
Let species 1 designate ethanol and species 2 designate water. We need to obtain an expression
for the molar volume, so first, convert the given the mass fractions and densities to mole
fractions and molar volumes.
w1
Mole fractions: x1
MW 1
w1 w2
MW 1 MW 2
The following graph plots the data. The trendline relates v to x1.
60.00
50.00
40.00
v (ml/mol)
2
v = 4.5491x1 + 35.918x1 + 17.957
30.00
R2 = 1
20.00
10.00
0.00
0.000 0.200 0.400 0.600 0.800 1.000
Mole Fraction EtOH (x1)
68
Now, we can calculate V1 .
V nv
V1
n1 n 2 , n 3 ,T , P n1 n 2 , n 3 ,T , P
We can substitute the trendline for V into this derivative and use the fact that ntot n1 n2 to
obtain
n12
V1 4.5491 35.918n1 17.957n1 n2
n1 n1 n2
n 2 ,T , P
Differentiating we get
2n n n n 2
V1 4.5491 1 1 2 1 53.875 ml
mol
n1 n2 2
ml
V1 4.5491 2 x1 x12 53.875
mol
Calculate V2 :
2
V2 4.5491 n1 35.918n1 17.957n1 n2
n2
n1 n2
n1 ,T , P
Differentiating we get
n2 ml
V2 4.5491 1
17.957
n1 n2 2 mol
ml
V2 4.5491x12 17.957
mol
69
Partial Molar Volumes vs. EtOH Mole Fractions
60 20
59 19
EtOH Partial Molar
Volume (ml/mol)
58 18
57 17
EtOH
56 16
H2O
55 15
54 14
53 13
52 12
0 0.2 0.4 0.6 0.8 1
Mole Fraction EtOH (x1)
(b)
v mix can be calculated using Equation 6.46
v mix x H 2 O V H 2 O v H 2 O x EtOH V EtOH v EtOH
ml
v EtOH v x1 1 58.36
mol
ml
v H 2 O v x1 0 18.05
mol
ml
V1 4.5491 20.5 0.52 53.875 57.29
mol
ml
V2 4.54910.52 17.957 16.82
mol
Therefore,
ml ml
v mix 0.516.82 18.05 0.5 57.29 58.36
mol mol
ml
vmix 1.15
mol
70
6.38
We can use the density data given in the problem statement to determine the pure species
properties. For pure ethanol x1 1 :
MW1 46 g/mol cm 3
v1
1
0.7857 g/cm 3
58.55
mol
cm 3
V1 n1v1 3 mol 58.55 175.7 cm
mol
3
MW2 45 g/mol cm 3
v2
2
1.1314 g/cm 3
39.77
mol
cm 3
V2 n2 v 2 1 mol 39.77 39.77 cm
3
mol
To calculate V and v, interpolate in the data table to obtain the density of the mixture
when x1 0.75 :
0.8550 g/cm 3
Therefore,
MW
v
where
MW 0.7546 g/mol 0.2545 g/mol 45.75 g/mol
v
45.75 g/mol
0.8550 g/cm 3
53.51 cm 3 /mol
cm 3
V nv 4 mol 53.51 214.0 cm
3
mol
Now, we can calculate the volume change of mixing by using Equation 6.43:
71
Vmix V V1 V2 214 cm 3 175.7 cm 3 39.77 cm 3
Vmix 1.47 cm 3
The intensive volume change of mixing:
v mix
Vmix 1.47 cm 3
0.368 cm 3 / mol
n 4
We can use Equation 6.65 to determine the partial molar volume of formamide:
dv v
V2 v x1 v x1
dx1 x1
45.8 45.7
v
0.8401 0.8701 19.50 cm 3 / mol
x1 0.8009 0.6986
Therefore,
V2 53.51 cm 3 / mol 0.75 19.5 cm 3 / mol 38.89 cm 3 / mol
Now calculate the partial molar volume of ethanol:
V n1V1 n2V2
V1
214.0 cm 3 1 mol 38.89 cm 3 /mol
58.37 cm 3 / mol
3 mol
72
6.39
Using the definition of G, the Gibbs energy of mixing of an ideal gas can be rewritten in terms of
the enthalpy of mixing and the entropy of mixing:
g mix
ideal gas
hmix
ideal gas
Tsmix
ideal gas
ideal gas
hmix 0
ideal gas
smix R xa ln xa xb ln xb
so
ideal gas
g mix RT xa ln xa xb ln xb
To find the partial molar Gibbs energy of mixing of species a, we apply Equation 6.29:
nb
RT na ln na nb ln nb na nb lnna nb
na
ng mix RT na ln nb ln
na nb na nb
where the mathematical relation of logarithms was used. Thus,
Gmix a a g a
As chemical engineers, we are often interested in the limiting case of infinite dilution. We see
that even for ideal gas mixtures the chemical potential in this limit is not mathematically well-
.behaved, In Chapter 7, we will develop a different function, the fugacity, which behaves better.
73
6.40
At equilibrium
H
l
O
H
v
2 2O
or
G Hl O G Hv O
2 2
GHl 2O g Hl 2O hHl 2O Ts Hl 2O
The enthalpy and entropy of liquids are not sensitive to pressure changes. We can use data from
the saturated steam tables at 25 C to determine the Gibbs energy.
l
g H 104.87 kJ/kg 298.15 K 0.3673 kJ/kg K
2O
l
g H 4.640 kJ/kg
2O
gHl
2O
l
MWH 2 O g H
2O
0.0180148 kg/mol 4.640 kJ/kg 0.0836 kJ/mol
Therefore,
H
v l
gH 83.6 J/mol
2O 2O
74