Introductory Chemical Engineering Thermodynamics: Chapter 7 - Departure Functions
Introductory Chemical Engineering Thermodynamics: Chapter 7 - Departure Functions
Engineering Thermodynamics
Chapter 7 - Departure Functions
U(T,V)-U (T,V) = ∫ ∂V
− dV
ig
∞ T ∂V T
FPR tells us (dU)T = T(dS)T - P(dV)T
For the ideal gas (dUig)T = T(dSig)T -P(dVig)T
ig
where (dSig)T = RdlnV and P(dV )T = RT/V dV = RT dlnV
ig
Substituting we find, (dU )T = T*R dlnV - RT dlnV = 0
∂U
ig
Therefore, ∂V =0
T
∂U ∂S ∂V
Returning to the FPR, ∂V = T ∂V − P ∂V
T T T
∂S ∂P
Maxwell’s Relation ∂V =
T ∂T V
∂P
V
(U − U ) = ∫ T
ig
− P dV
Finally,
∞
∂T V
A(T , V ) − A (T , V ) = ∫
ig
− dV
∞ ∂V T ∂V T
∂A ∂lnV − RT
ig
∂V T ∂V T
T
V V
∂A ∂A ∂A
ig
RT
= −P ⇒ −
= − P +
Also ∂V
T ∂ V T ∂ V T V
Transform to ρ⇒dV = -V dρ/ρ
ρ
A(T ,V ) − A ig (T ,V ) Z −1
⇒ = ∫ dρ
RT 0 ρ
(H − H ) = ig ρ
∂Z dρ (G (T , P ) − G (T , P )) = (Z − 1) dρ + Z − 1 − lnZ
ig
∫
ρ
RT ∫o ∂T ρ + Z − 1
− T
RT o
ρ
(U − U ) ig ρ
∂ Z dρ (A(T , V ) − A (T , V )) = (Z − 1) dρ
ig ρ
RT ∫o ∂T ρ
− T
RT ∫
o
ρ
(S (T , V ) − S (T , V )) =
ig ρ
∂Z
( )
dρ (S (T , P ) − S (T , P )) =
ig ρ
∂Z
( )
dρ
∫o ∂T p
− T − Z − 1 ∫o ∂T p
− T − Z − 1
ρ
+ lnZ
R ρ R
Solution:
Path, for H(190,25) - H(105,5)
[H(190,25) - Hig(190,25)]+[Hig(190,25) - Hig(105,5)]+[Hig(105,5) - H(105,5)]
Similarly for S(190,25) - S(105,5)
[S(190,25 - Sig(190,25)]+[Sig(190,25) - Sig(105,5)]+[Sig(105,5) - S(105.5)]
I. Departure Function + II. Ideal gas + III. Departure function
I. (190,25) → (190,25)ig
190 + 273.15 = 463.15K & 25 bar ⇒ Tr = 1.25135; Pr = 0.58837
Solution: In this example, we are directed to use a reference state such that, for enthalpy:
H2 - H1 = (H2 - Href) - (H1 - Href), and for entropy: S2 - S1 = (S2 - Sref) - (S1 - Sref). Note the
equivalence of this procedure to the way steam tables are computed. Furthermore, the
computation of H2 - Href or S2 - Sref is entirely equivalent to the procedure given in
Example 7.1.
1. REF: Enter the values of Tc , Pc , ω, A, B, C, D and define the T, P, and root of
interest.
2. Press PVTF to enter the pressure and temperature and choose the root of interest.
E.g. at 463.15 K and 2.5 MPa, V = 1369 cm3/mole
3. Press UHSG to compute internal energy, enthalpy, entropy, and Gibbs free energy.
E.g. at 463.15 K and 2.5 MPa, H2 - Href = 36901 J/mole; S2 - Sref = 109.15 J/mole-K
4. Repeat at 378.15 K and 0.5 MPa, H1 - Href = 29586 J/mole; S1 - Sref = 104.13 J/mole-K
5. Subtract ⇒ ∆H = 36901-26756 = 7315 and ∆S = 109.15-104.13 = 5.02 J/mole-K
[ ( )]
2
a = ac 1 + κ 1 − T / Tc ; ac ≡ 0.45724
R 2Tc 2
Pc
where κ = 0.37464 + 1.54226 ω - 0. 26993 ω 2
T − κ 1 −1 / 2
da
dT
= ac 2 1 + κ 1 −
Tc
2
T
T
da
[ (
⇒ dT = − a c 1 + κ 1 − T / Tc κ T / Tc)]( )
Tc
∂Z
−T =
bρ −a
−
ac
− 1[ (
+ κ 1 − T / Tc )](
κ T / Tc
≡) bρ
F (Tr )
(1 + 2bρ − b 2 ρ 2 )
2 2
∂T ρ 1 + 2bρ − b ρ bRT bRT
F(Tr) is shorthand. Also B ≡ bP/RT ⇒ bρ = B/Z and A ≡ aP/R2T2 ⇒ a/bRT = A/B
= RT = bρ ⇒ ∫ − T
∂Z d (bρ ) F (Tr ) Z + 1 + 2 B
= ln
( )
Note Z p 0 ∂ T T b ρ 8 Z + 1 −( 2 B
)
ρRT
(H − H ) = Z − 1 +
ig
1
ln
(1 + 2 )B − a − a [1 + κ (1 −
Z +
)]
Tr κ Tr
(1 − 2 )B bRT bRT
c
nRT 8 Z +
l
( 2 )B A 1 + κ Tr
Z + 1+
= Z - 1 - n
( 2 )B B 8 α
Z + 1-
nRT ∫
o
ρ
(
1 + 1 + 2 b ρ
a )
=ln(1-bρ) – lnZ + Z – 1 + bRT 8 ln 1 + 1 - 2 bρ
( )
(
Z + 1+ 2 B A)
=ln(Z-B)- ln Z + 1 - 2 B B 8
( )