Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Assignment 6: Solutions

Download as pdf or txt
Download as pdf or txt
You are on page 1of 4

Assignment 6: Solutions Page |1

Q1
From (60) of the summary for Week 6

𝐻𝑅 𝑑 ln 𝛼(𝑇𝑟 ) (A)
=𝑍−1+[ − 1] 𝑞𝐼
𝑅𝑇 𝑑 ln 𝑇𝑟
We need to find Z and other terms on the RHS of (A) to evaluate the residual enthalpy.

The compressibility factor Z is calculated from the EoS (Peng-Robinson) in this case.

From (58) of summary for Week 5, Z is a root of the following cubic equation.

𝑍 3 + 𝑍 2 [(𝜖 + 𝜎)𝛽 − 1 − 𝛽] + 𝑍𝛽[(𝑞 − 𝜖 − 𝜎) − (𝜖 + 𝜎 − 𝜖𝜎)𝛽] − 𝛽 2 (𝑞 − 𝜖𝜎𝛽 − 𝜖𝜎) = 0 (B)

For PR EoS, again from (58), the following equations hold.

𝜶(𝑻𝒓 ) 𝝈 𝝐 𝛀 𝚿
PR 𝛼𝑃𝑅 (𝑇𝑟 ; 𝜔) 1 + √2 1 − √2 0.0778 0.45724

1/2 2 (C)
𝛼𝑃𝑅 (𝑇𝑟 ; 𝜔) = [1 + (0.37464 + 1.54226𝜔 − 0.26992𝜔2 ) (1 − 𝑇𝑟 )]

𝑃𝑟 Ψ𝛼(𝑇𝑟 ) (D)
β=Ω 𝑞=
𝑇𝑟 Ω𝑇𝑟

From the above equations the following calculations follow. (E)

𝑃 25
𝑃𝑟 = = = 0.5123
𝑃𝑐 48.8

𝑇 300
𝑇𝑟 = = = 0.982
𝑇𝑐 305.5

1/2 2
𝛼𝑃𝑅 = [1 + (0.37464 + 1.54226𝜔 − 0.26992𝜔2 ) (1 − 𝑇𝑟 )]
2
=[1 + (0.37464 + 1.54226 × 0.099 − 0.26992 × 0.0992 )(1 − √0.982)]
=1.0095
Assignment 6: Solutions Page |2

𝑃𝑟 0.5122
β=Ω = 0.0778 × = 0.040587
𝑇𝑟 0.982
Ψ𝛼(𝑇𝑟 ) 0.45724 × 1.0095
𝑞= = = 6.0417
Ω𝑇𝑟 0.0778 × 0.982

The cubic in Z, i.e. Eq. (B) reduces to


𝑍 3 + 0.959413 𝑍 2 + 0.159103𝑍 − 0.01167 = 0

The only real root for this equation is, Z = 0.7731

The value of I, from (62) in summary for week 6 is given by

1 𝑍 + 𝜎𝛽
If 𝜖 ≠ 𝜎, 𝑡ℎ𝑒𝑛 𝐼 = 𝑙𝑛 ( )
𝜎−𝜖 𝑍 + 𝜖𝛽
1 0.7731 + 0.040587(1 + √2)
= 𝑙𝑛 ( ) = 0.049964
(1 + √2) − (1 − √2) 0.7731 + 0.040587(1 − √2)

From (65) in summary for week 6,


𝑑 ln 𝛼(𝑇𝑟 )
= −𝜁√𝑇𝑟 ⁄𝛼 = (0.37464 + 1.54226𝜔 − 0.26992𝜔2 )√𝑇𝑟 ⁄𝛼
𝑑 ln 𝑇𝑟
= (0.37464 + 1.54226 × 0.099 − 0.26992 × 0.0992 )√0.982⁄1.0095
= −0.5175

𝐻𝑅 𝑑 ln 𝛼(𝑇𝑟 )
=𝑍−1+[ − 1] 𝑞𝐼 = 0.7731 − 1 + [−0.5175 − 1] × 6.0417 × 0.049964
𝑅𝑇 𝑑 ln 𝑇𝑟
= −0.685
𝐻 𝑅 = −0.685𝑅𝑇 = −0.685 × 8.314 × 300 = −1708.5 J mol−1
Assignment 6: Solutions Page |3

Q2
At the given conditions

𝑃 29.28
𝑃𝑟 = = = 0.6
𝑃𝑐 48.8

𝑇 290.23
𝑇𝑟 = = = 0.95
𝑇𝑐 305.5

0 1
𝑆𝑅 𝑆𝑅
From Lee/Kesler tables we obtain ( ) , ( ) at the desired reduced conditions.
𝑅 𝑅

0 1
𝑆𝑅 𝑆𝑅 𝑆𝑅
= ( ) + 𝜔 ( ) = −0.671 + 0.099(−0.961) = −0.766
𝑅 𝑅 𝑅
𝑆 𝑅 = −0.766 × 8.314 = −6.4 J mol−1 K −1
Assignment 6: Solutions Page |4

Q3
From the Rackett equation

(1−𝑇𝑟 )2/7
𝑉 𝑠𝑎𝑡 = 𝑉𝑐 𝑍𝑐
To find 𝑉𝑐 we will use conditions at the critical point

𝑍𝑐 𝑅𝑇𝑐 0.242 × 83.14 bar cm3 mol−1 K −1 × 405.7 K


𝑉𝑐 = = = 72.16 cm3 mol−1
𝑃𝑐 112.8 bar

𝑇 293.15
The reduced temperature 𝑇𝑟 = = = 0.7226
𝑇𝑐 405.7

So,

(1−𝑇𝑟 )2/7 2/7


𝑉 𝑠𝑎𝑡 = 𝑉𝑐 𝑍𝑐 = 72.1 × 0.242(1−0.7226) = 27.1cm3 mol−1

You might also like