blob: e4835b4d4c8b8b1eca2553b1319046a20c417a2c (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
|
// Y.V.C.Rao ,1997.Chemical Engineering Thermodynamics.Universities Press,Hyderabad,India.
//Chapter-3,Example 11,Page 72
//Title:Volume using two paramter and three parameter compressibility factor correlation
//================================================================================================================
clear
clc
//INPUT
T=409.41;//temperature of n-octane in degree celsius
P=4.98;//pressure in bar
Tc=569.4;//critical temperature of n-octane in K
Pc=24.97;//critical pressure of n-octane in bar
w=0.398;//acentric factor (no unit)
R=8.314;//universal gas constant in (Pa m^3)/(mol K)
//CALCULATION
Tr=(T+273.15)/Tc;//calculation of reduced temperature (no unit)
Pr=P/Pc;//calculation of reduced pressure (no unit)
z0=0.96;//value of compressibilty factor at tr and pr obtained from Fig.3.11
V1=(z0*R*(T+273.15))/(P*10^5);//calculation of volume in m^3/mol using the two parameter compressibilty factor correlation
z1=0.01;//value of compressibilty factor at tr and pr obtained from Fig.3.12
z=z0+(w*z1);//calculation of compressibility factor using Eq.3.64
V2=(z*R*(T+273.15))/(P*10^5);//calculation of volume in m^3/mol using the three parameter compressibility factor correlation
//OUTPUT
mprintf('\n The volume occupied by n-octane obtained by the two parameter compressibilty factor correlation= %f m^3/mol\n',V1);
mprintf('\n The volume occupied by n-octane obtained by the three parameter compressibility factor correlation= %f m^3/mol\n',V2);
//===============================================END OF PROGRAM===================================================
|
