例:求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