Volumetric Properties of Pure Fluids

A pressure-temperature diagram
• the sublimation curve
• the fusion curve
• the vaporization curve
• the triple point
• the critical point

Fig 3.1
 A pressure-volume diagram
• The isotherms
– the subcooled-liquid and the
superheated-vapor regions
– isotherms in the subcooled-
liquid regions are steep
because liquid volumes change
little with large change in
pressure
• The two-phase coexist region
• The triple point is the horizontal
line
• The critical point
Fig 3.2
 • An equation of state exists relating pressure, molar or specific
volume, and temperature for any pure homogeneous fluid in
equilibrium states.
• An equation of state may be solved for any one of the three
quantities P, V, or T as a function of the other two.
• Example:  V   V 
dV =   dT +   dP
 T P  P T
1  V  1  V 
Volume expansivity:     Isothermal compressibility:   −  
V  T  P V  P T
dV
= dT − dP
V
– For incompressible fluid, both β andκ are zero.
– For liquids β is almost positive (liquid water between 0°C and 4°C is
an exception), and κ is necessarily positive.
– At conditions not close to the critical point, β andκ can be assumed
constant: V2
ln =  (T2 − T1 ) −  ( P2 − P1 )
V1
 Virial equations of state
• PV along an isotherm:
• PV = a + bP + cP 2 = a(1 + BP + C P 2 + DP 3 + ...)
– The limiting value of PV as P →0 for all the gases:
– (PV )* = a = f (T )
– (PV )* = a = RT , with R as the proportionally constant.
– Assign the value of 273.16 K to the temperature of the triple point
of water: (PV )* = R  273.16
t
• Ideal gas:
– the pressure ~ 0; the molecules are separated by infinite distance;
the intermolecular forces approaches zero.

– R=
(PV )*t
= 83.1447
cm3 bar
273.16 mol K
 The virial equations
PV
• The compressibility factor: Z
RT

• the virial expansion: Z = (1 + BP + C P 2 + DP 3 + ...)

– the parameters B’, C’, D’, etc. are virial coefficients,

accounting of interactions between molecules.
– the only equation of state proposed for gases having a
firm basis in theory.
– The methods of statistical mechanics allow derivation
of the virial equations and provide physical significance
to the virial coefficients.
 Ideal gas
• No interactions between molecules.
• gases at pressure up to a few bars may often be considered
ideal and simple equations then apply
• the internal energy of gas depends on temperature only.
– Z = 1; PV = RT
– U = U (T)
– C   U  = dU = C (T ) C   H  = dH = dU + R = C + R = C (T )
 T V  T  P
V V P V P
dT dT dT

U =  CV dT H =  CP dT
– Mechanically reversible closed-system process, for a unit mass or a
mole:
CV C
dQ + dW = CV dT dQ = VdP + P PdV
dW = − PdV R R
• For ideal gas with constant heat capacities undergoing a
mechanically reversible adiabatic process:
R
T2  P2  CP
dQ + dW = CV dT =  
T1  P1 

dW = − PdV PV  = const.
C
dT R dV T2  V1 
R
CV  P
=− =   CV
T CV V T1  V2 
– for monatomic gases,  = 1.67
– for diatomic gases,  = 1.4
– for simple polyatomic gases, such as CO2, SO2, NH3, and CH4,  = 1.3
• The work of an irreversible process is calculated:
– First, the work is determined for a mechanically reversible process.
– Second, the result is multiple or divided by an efficiency to give
the actual work.
Air is compressed from an initial condition of 1 bar and 25°C to a final state of 5
bar and 25 °C by three different mechanically reversible processes in a closed
system. (1) heating at constant volume followed by cooling at constant pressure; (2)
isothermal compression; (3) adiabatic compression followed by cooling at constant
volume. Assume air to be an ideal gas with the constant heat capacities, CV =
(5/2)R and CP = (7/2)R. Calculate the work required, heat transferred, and the
changes in internal energy and enthalpy of the air in each process.

Fig 3.7
Choose the system as 1 mol of air, contained in an imaginary frictionless piston
/cylinder arrangement.

For R = 8.314 J/mol.K, CV = 20.785, CP = 29.099 J/mol.K

The initial and final molar volumes are: V1 = 0.02479 m3 and V2 = 0.004958 m3
The initial and final temperatures are identical: ΔU = ΔH = 0

(1) Q = CVΔT + CPΔT = -9915 J; W = ΔU - Q = 9915 J

 P1 
(2) Q = −W = RT ln   = −3990 J
 P2 
 −1 
 V1  V 
(3) adiabatic compression: 2 1  
T = T = 567.57 K P2 = P1  1  = 9.52 bar
 V2   V2 
W = CV T = 5600 J

cooling at constant V, W = 0.
overall, W = 5600 J, Q = ΔU - W = -5600 J.
An ideal gas undergoes the following sequence of mechanically reversible
processes in a closed system:
(1) From an initial state of 70°C and 1 bar, it is compressed adiabatically to 150 °C.
(2) It is then cooled from 150 to 70 °C at constant pressure.
(3) Finally, it is expanded isothermally to its original state.
Calculate W, Q, ΔU, and ΔH for each of the three processes and for the entire cycle.
Take CV = (3/2)R and CP = (5/2)R. If these processes are carried out irreversibly but
so as to accomplish exactly the same changes of state (i.e. the same changes in P, T,
U, and H), then different values of Q and W result. Calculate Q and W if each step
is carried out with an efficiency of 80%.

Fig 3.8
Choose the system as 1 mol of air, contained in an imaginary frictionless piston
/cylinder arrangement. For R = 8.314 J/mol.K, CV = 12.471, CP = 20.785 J/mol.K
(1) For an ideal gas undergoing adiabatic compression, Q = 0
ΔU = W = CVΔT = 12.471(150 – 70) = 998 J
ΔH = CPΔT = 20.785(150 – 70) = 1663 J

 T2  ( −1)
P2 = P1   = 1.689 bar
 T1 
(2) For the constant-pressure process:
Q = ΔH = CPΔT = 20.785(70 – 150) = -1663 J
ΔU = CVΔT = 12.471(70 – 150) = -998 J
W = ΔU – Q = 665 J
(3) Isotherm process, ΔU and ΔH are zero:
 P3 
Q = −W = RT ln   = 1495 J
 P1 
(4) Overall: Q = 0 – 1663 + 1495 = -168 J
W = 998 + 665 – 1495 = 168 J
ΔU = 0
ΔH = 0
Irreversible processes:
(1) For 80% efficiency:
W(irreversible) = W(reversible) / 0.8 = 1248 J
ΔU(irreversible) = ΔU(reversible) = 998 J
Q(irreversible) = ΔU – W = -250 J
(2) For 80% efficiency:
W(irreversible) = W(reversible) / 0.8 = 831 J
ΔU = CVΔT = 12.471(70 – 150) = -998 J
Q = ΔU – W = -998 – 831 = -1829 J
(3) Isotherm process, ΔU and ΔH are zero:
W(irreversible) = W(reversible) x 0.8 = -1196 J
Q = ΔU – W = 1196 J
(4) Overall: Q = -250 – 1829 + 1196 = -883 J
W = 1248 + 831 – 1196 = 883 J
ΔU = 0
ΔH = 0
The total work required when the cycle consists of three irreversible steps is
more than 5 times the total work required when the steps are mechanically
reversible, even though each irreversible step is assumed 80% efficient.
A 400g mass of nitrogen at 27 °C is held in a vertical cylinder by a frictionless
piston. The weight of the piston makes the pressure of the nitrogen 0.35 bar higher
than that of the surrounding atmosphere, which is at 1 bar and 27°C. Take CV =
(5/2)R and CP = (7/2)R. Consider the following sequence of processes:
(1) Immersed in an ice/water bath and comes to equilibrium
(2) Compressed reversibly at the constant temperature of 0°C until the gas volume
reaches one-half the value at the end of step (1) and fixed the piston by latches
(3) Removed from the ice/water bath and comes to equilibrium to thermal
equilibrium with the surrounding atmosphere
(4) Remove the latches and the apparatus return to complete equilibrium with its
surroundings.
Nitrogen may be considered an ideal gas. Calculate W, Q, ΔU, and ΔH for each step
of the cycle.
The steps: m
(1) 27 C , 1.35bar ⎯const n= = 14.286 mol
⎯⎯ P
→ 0 C , 1.35bar M
1
(2) 0 C , V2 ⎯const
⎯⎯ T
→ 0 C , V3 = V2
2
(3) 0 C , V3 ⎯const
⎯⎯ V
→ 27 C , V4 = V3
4 =T1
(4) 27 C , P4 ⎯T⎯⎯→ 27 C , 1.35bar
 Fig 3.9

(1) W = −n PdV = −nPV = −nRT = 3207 J Q = nH = nC T = 11224 J

1  1 1 P

nU1 = Q1 + W1 = −11224 + 3207 = −8017 J

(2) U = H = 0 V3
2 2
Q2 = −W2 = nRT ln = −22487 J
V2

(3) W = 0 Q = nU = nC T = 8017 J nH = nC T = 11224 J

3 3 3 V 3 P

(4) the oscillation of the piston

U 4 = H 4 = 0 Q4 = −W4
Air flows at a steady rate through a horizontal insulated pipe which
contains a partly closed valve. The conditions of the air upstream
from the valve are 20°C and 6 bar, and the downstream pressure is 3
bar. The line leaving the valve is enough larger than the entrance line
so that the kinetic-energy change as it flows through the valve is
negligible. If air is regarded as an ideal gas, what is the temperature
of the air some distance downstream from the valve?
Flow through a partly closed valve is known as a throttling process.
For steady flow system:
d (mU ) cv  1  
+   H + u 2 + zg m  = Q + W H = 0
dt  2   fs

Ideal gas: H =  CP dT T2 = T1

The result that ΔH = 0 is general for a throttling process.

If the flow rate of the air is 1 mol/s and if the pipe has an inner diameter of 5 cm,
both upstream and downstream from the valve, what is the kinetic-energy change
of the air and what is its temperature change? For air, CP = (7/2)R and the molar
mass is M = 29 g/mol.

Upstream molar volume:

RT1 83.14  293.15 −6 −3 m 2 1 V
V1 = = 10 = 4.062 10 u = n = n = 2.069 m
mol A s
1
P1 6 A
Downstream molar volume:
V2 = 2V1 u = 2u = 4.138 m
2 1 s
The rate of the change in kinetic energy:
1 2 1 2 −3 ( 4.138 − 2.069 )
2 2
m  u  = n M  u  = (1 29  10 ) = 0.186 J
2  2  2 s

d (mU ) cv  1   C  1 
+   H + u 2 + zg m  = Q + W m  P T +  u 2   = 0
dt  2   fs M  2 

T = −0.0064 K
 Application of the virial equations
• Differentiation:
 Z   Z 
  = B + 2C P + 3DP + ...
2
  = B
 P T  P T ; P =0

• the virial equation truncated to two terms satisfactorily

represent the PVT behavior up to about 5 bar
BP B
Z = 1 + BP Z = 1+
RT
= 1+
V
• the virial equation truncated to three terms provides good
results for pressure range above 5 bar but below the critical
pressure
B C
Z = 1 + BP + C P 2 Z = 1+ + 2
V V
Reported values for the virial coefficients of isopropanol vapor at 200°C are:
B = -388 cm3/mol and C = -26000 cm6/mol2. Calculate V and Z for isopropanol
vapor at 200 °C and 10 bar by (1) the ideal gas equation; (2) two-term virial
equation; (3) three-term virial equation.

(1) For an ideal gas, Z = 1:

RT 83.14  473.15 3
V= = = 3934 cm
P 10 mol
(2) two-term virial equation:
RT 3 PV
V= + B = 3934 − 388 = 3546 cm Z= = 0.9014
P mol RT
(3) three-term virial equation:
RT  B C   388 − 26000  cm3
Vi +1 = 1 + + 2  = 39341 + + 
2 
= 3539 1st iteration
 3934 (3934)  mol
P  Vi Vi 

Ideal gas value

...
3 PV
After 5 iterations V4 ~ V5 = 3488 cm Z= = 0.8866
mol RT
 Cubic equations of state
• Simple equation capable of representing both
liquid and vapor behavior.
RT a
• The van del Waals equation of state: P = −
V −b V 2
– a and b are positive constants
– unrealistic behavior in the two-phase region. In reality,
two, within the two-phase region, saturated liquid and
saturated vapor coexist in varying proportions at the
saturation or vapor pressure.
– Three volume roots, of which two may be complex.
– Physically meaningful values of V are always real,
positive, and greater than constant b.
Fig 3.12
A generic cubic equation of state
• General form: RT  (V −  )
P= −
V − b (V − b)(V 2 + V +  )

– where b, θ, κ,λ and η are parameters depend on temperature and

(mixture) composition.
– Reduce to the van der Waals equation when η= b, θ= a, and κ=λ=
0.
– Set η= b, θ= a (T), κ= (ε+σ) b, λ = εσb2, we have:

RT a(T )
P= −
V − b (V + b)(V + b)
• where ε and σ are pure numbers, the same for all substances, whereas
a(T) and b are substance dependent.
• Determination of the parameters:
– horizontal inflection at the critical point:
 P   2P 
  =  2  = 0
 V T ;cr  V T ;cr
• 5 parameters (Pc, Vc, Tc, a(Tc), b) with 3 equations, one has:
3 RTc 27 R 2Tc2 1 RTc
Vc = a= b=
8 Pc 64 Pc 8 Pc
PcVc 3
• Zc  =
RTc 8
• Unfortunately, it does not agree with the experiment. Each
chemical species has its own value of Zc.
• Similarly, one obtain a and b at different T.

RT a (T )  (Tr ) R 2Tc2 RTc

P= − a(T ) =  b=
V − b V (V + b) Pc Pc
 Two-parameter and three-parameter
theorems of corresponding states
• Two-parameter theorem: all fluids, when compared at the same
reduced temperature and reduced pressure, have approximately the
same compressibility factor, and all deviate from ideal-gas behavior to
about the same degree.
T P
• Define reduced temperature and reduced pressure: Tr  Pr 
Tc Pc

• Not really enough to describe the state, a third corresponding-states

parameter is required.
– The most popular such parameter is the acentric factor (K.S. Pitzer, 1995):

  −1.0 − log (Prsat )T =0.7

r

• Three-parameter theorem: all fluids having the same value of ω, when

compared at the same reduced temperature and reduced pressure, and
all deviate from ideal-gas behavior to about the same degree.
 • Vapor and vapor-like
RT a(T ) V −b V starts with V(ideal-gas)
V= +b−
P P (V + b)(V + b) and then iteration
a(T )  (Tr ) bP P
q =  = r
bRT Tr RT Tr
Z −
Z = 1 +  − q
( Z +  )(Z +  )
• Liquid and liquid-like
 RT + bP − VP  V starts with V = b and
V = b + (V + b)(V + b)  
 a(T )  then iteration
a(T )  (Tr ) bP P
q =  = r
bRT Tr RT Tr
1+  − Z 
Z =  + ( Z +  )( Z +  ) 
 q 

Equations of state which express Z as a function of Tr and Pr are said to be

generalized, because of their general applicability of all gases and liquids.
 2-parameter/3-parameter E.O.S.
• Express Z as functions of Tr and Pr only, yield 2-
parameter corresponding states correlations:
– The van der Waals equation
– The Redlich/Kwong equation
• The acentric factor enters through function α(Tr;ω)
as an additional paramter, yield 3-parameter
corresponding state correlations:
– The Soave/Redlich/Kwong (SRK) equation
– The Peng/Robinson (PR) equation
Table 3.1
Given that the vapor pressure of n-butane at 350K is 9.4573 bar, find the molar
volumes of (1) saturated-vapor and (2) saturated-liquid n-butane at these conditions
as given by the Redlich/Kwong equation.
350 9.4573
Tr = = 0.823 Pr = = 0.2491
425.1 37.96
 (Tr ) Pr
q= = 6.6048  = = 0.026214
Tr Tr

(1) The saturated vapor

Z −
Z = 1 +  − q Z starts at Z = 1 and converges on Z = 0.8305
( Z +  )(Z +  )
ZRT cm3
V= = 2555
P mol
(2) The saturated liquid
1+  − Z 
Z =  + ( Z +  )( Z +  )  Z starts at Z = β and converges on Z = 0.04331
 q 
ZRT cm3
V= = 133.3
P mol
 Generalized correlations for gases
• Pitzer correlations for the compressibility factor:
Z = Z 0 + Z 1
– Z0 = F0 (Tr, Pr)
– Simple linear relation between Z and ω for given values
of Tr and Pr.
– Of the Pitzer-type correlations available, the Lee/Kesler
correlation provides reliable results for gases which are
nonpolar or only slightly polar (App. E).
– Only tabular nature (disadvantage)
 Pitzer correlations for the 2nd virial
coefficient
BP 0 Pr 1 Pr
• Correlation: Z = 1+ = 1+ B + B Z = Z 0 + Z 1
RT Tr Tr
Pr P
Z 0 = 1+ B 0 Z 1 = B1 r
Tr Tr
– Validity at low to moderate pressures
– For reduced temperatures greater than Tr ~ 3, there
appears to be no limitation on the pressure. 0 0.422
B = 0.083 − 1.6
– Simple and recommended. Tr
– Most accurate for nonpolar species. 0.172
1
B = 0.139 −
Tr4.2
 Determine the molar volume of n-butane at 510K and 25 bar by, (1) the ideal-gas
equation; (2) the generalized compressibility-factor correlation; (3) the generalized
virial-coefficient correlation.
(1) The ideal-gas equation
RT cm3
V= = 1696.1
P mol
(2) The generalized compressibility-factor correlation
510 25 the acentric factor
Tr = = 1.200 Pr = = 0.659  = 0.200
425.1 37.96

the Lee/Kesler correlation

ZRT cm3
Z = 0.865
0
Z = 0.038
1
Z = Z + Z = 0.873 V =
0 1
= 1480.7
P mol
(3) The generalized virial-coefficient correlation
510 0.422 1 0.172
Tr = = 1.200 B = 0.083 − 1.6
0
B = 0.139 −
425.1 Tr Tr4.2

Pr P 3
Z = 1 + B0 + B1 r = 0.879 V = ZRT = 1489.1 cm
Tr Tr P mol
What pressure is generated when 1 (lb mol) of methane is stored in a volume of 2
(ft)3 at 122°F using (1) the ideal-gas equation; (2) the Redlish/Kwong equation; (3)
a generalized correlation .
(1) The ideal-gas equation
RT 0.7302(122 + 459.67)
P= = = 212.4 atm
V 2
(2) The RK equation
581.67  (Tr ) R 2Tc2 atm RTc
Tr = = 1.695 a (T ) =  = 453.94 b =  = 0.4781 ft 3

343.1 Pc ft 6 Pc
RT a(T )
P= − = 187.49 atm
V − b V (V + b)
(3) The generalized compressibility-factor correlation is chosen (high pressure)
ZRT Z (0.7302)(122 + 459.67)
P= = = 212.4Z atm
V 2
P Z 581.67 Z starts at Z = 1 and
Pr = = Tr = = 1.695
45.4 0.2138 343.1 converges on Z = 0.890
P = 189.0 atm
A mass of 500 g of gases ammonia is contained in a 30000 cm3 vessel immersed in
a constant-temperature bath at 65°C. Calculate the pressure of the gas by (1) the
ideal-gas equation; (2) a generalized correlation .
Vt cm3
V= = 1021.2
n mol
(1) The ideal-gas equation
RT
P= = 27.53 bar
V

(2) The generalized virial-coefficient correlation is chosen (low pressure, Pr ~ 3 )

338.15 0.422 0.172
Tr = = 0.834 B 0
= 0.083 − B1 = 0.139 −
405.7 Tr1.6 Tr4.2
27.53
Pr ~ = 0.244 the acentric factor
112.8  = 0.253

(
Z = 1 + B 0 + B1 ) TP
r
= 1 − 0.541
Pr
Tr
r

ZRT
P= = 23.76 bar
V
 Generalized correlations for liquids
• The generalized cubic equation of state (low accuracy)
• The Lee/Kesler correlation includes data for subcooled
liquids
– Suitable for nonpolar and slightly polar fluids
• Estimation of molar volumes of saturated liquids
– Rackett, 1970: V sat = V Z (1−Tr )0.2857
c c

• Generalized density correlation for liquid (Lydersen,

Greenkorn, and Hougen, 1955):  Vc
r  =
c V
 Fig 3.17

V2 = V1 r1
r 2
For ammonia at 310 K, estimate the density of (1) the saturated liquid; (2) the
liquid at 100 bar
(1) Apply the Rackett equation at the reduced temperature

Tr =
310 V = 72.47 Z c = 0.242
= 0.7641 c
405.7

cm3
V sat = Vc Z c(1−Tr ) = 28.33
0.2857

mol
(2) At 100 bar
100
Pr = = 0.887
112.8 Fig 3.17
 r = 2.38 Vc cm3
V= = 30.45
Tr =
310
= 0.7641 r mol
405.7

 r1  r1,310K , saturated liquid 2.34 cm3

V2 = V1 =V 310 K
= 29.14 = 28.65
r 2  r 2,100bar 2.38 mol