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

例：求C4H10在460K，1.520MPa下的逸度

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

解： R-K法：

Tc 425KPc 3.80MPaw 0.193

.56/23.8 10 2.8955 10a 0.42748R2Tc2./5Pc 0.42748 8.314 2425

b 0.08664RTc/Pc 0.08664 8.314 425/3.8 106 80.56 10 6

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