化工热力学逸度计算_Matlab程序

例：求410C H 在460K ，1.520MPa 下的逸度

方法 1. I.g 2. R-K 3. 普遍化关联

解： R-K 法：

22.522.5666

425 3.800.193

0.42748/0.427488.314425/3.810 2.8955100.08664/0.086648.314425/3.81080.5610c c c c c c T K P MPa w a R T P b RT P -=====⨯⨯⨯=⨯==⨯⨯⨯=⨯

v=Z*R*T/P

输入命令开始迭代

Z=0.5;ff=myfun1(Z);

while ff>0.001

Z=Z+0.0001;

if Z<1

ff=myfun1(Z)

myfun2(Z)

end

end

结果为：压缩因子：0.8835

正丁烷的摩尔体积：0.002223

正丁烷的逸度系数：0.89431

正丁烷的逸度：1359345.5333

2．普遍化关联法

T=460;P=1.52e+6; %普遍化关联法

Tc=425;Pc=3.8e+6;w=0.193;R=8.314;

Tr=T/Tc;Pr=P/Pc;

BO=0.083-0.422/(Tr^1.6);

B1=0.139-0.172/(Tr^4.2);

B=BO+w*B1;

xishu=exp(B*Pr/Tr);

f=P*xishu;

disp(['正丁烷的逸度系数：',num2str(xishu)]);

disp(['正丁烷的逸度：',num2str(f)]);

正丁烷的逸度系数：0.89977

正丁烷的逸度：1367646.3748