6并且 p?pc?pr?4.599?10pr

4.599?106pr?0.214pr （A） 因此，结合上两式得：Z?72.149?10Pitzer的普遍化关系式为：Z?Z(0)??Z(1) （B）

根据（A）、（B）两式进行迭代，过程为： （1） 设Z值，然后代入（A）式求出pr；

（2） 根据Tr和pr值查（2－9）和（2－10）得到Z?0?和Z?1?； （3） 将查图得到的Z?0?和Z?1?值代入（B）式求得Z值；

（4） 比较Z的计算值与实验值，如果相差较大，则代入（A）式重新计算，直到迭代收

敛。

依据上述迭代结果为：pr＝4.06时，Z＝0.877

661.867?107Pa 则：p?pc?pr?4.599?10pr＝4.599?10?4.06＝误差： ?1.875?1.867??10/1.875?10?0.43%

772－22．试用RK方程计算二氧化碳和丙烷的等分子混合物在151℃和13.78 MPa下的摩尔体积。

解：计算所需的数据列表如下： 组元 Tc/K pc?10?5/PaVc?106/(m3?mol-1)Zc 94.0 200 140.4 a/(Pa?m6?K0.5?mol-2)b/(m3?mol-1) CO2（1） 304.2 C3H8（2） 369.8 12 335.4 73.82 42.48 54.72 0.274 6.460 0.277 18.29 11.12 2.968×10-5 6.271×10-5 由（2－51a）和（2－51b）得：

2aM?y12a11?2y1y2a12?y2a22＝0.52?6.460?2?0.5?0.5?11.12?0.52?18.29＝11.75Pa?m?K60.5?mol-2

bM?y1b1?y2b2＝0.5?2.968?10?5?0.5?6.271?10?5＝4.415?10?5m3?mol-1

??aMp11.75?13.78?106A?22.5＝＝0.6322 22.5RT?8.314???424.15?bMp4.415?10?5?13.78?106B?＝＝0.1725

RT8.314?424.15按照式（2-16a）Z?1A?h?1?h????3.6084??＝? （A） 1?hB?1?h?1?h1?h??bB0.1725?? （B） VmZZ和式（2－16b） h?联立求解方程（A）、）（B）进行迭代计算得： 迭代次数 0 1 2 3 4 5 因此：Z＝0.5975，h＝0.2887 混合物得摩尔体积为：V?Z 1 0.6776 0.6093 0.5987 0. 5976 0.5975 h 0.1725 0.2546 0.2831 0.2881 0.2887 0.2887 ZRT0.5975?8.314?424.15?43??1.52?10m?mol 6p13.78?102－23．混合工质的性质是人们有兴趣的研究课题。试用RKS状态方程计算由R12（CCl2F2）和R22（CHClF2）组成的等摩尔混合工质气体在400K和1.0MPa，2.0MPa，3.0 MPa，4.0 MPa和5.0 MPa时的摩尔体积。可以认为该二元混合物的相互作用参数k12＝0（建议自编软件计算）。计算中所使用的临界参数如下表 组元（i） R22（1） R12（2） Tc/K 369.2 385 pc/ MPa 4.975 4.224 ? 0.215 0.176 解：计算过程是先计算两个纯组分的RKS常数，再由混合规则获得混合物的RKS常数

aM,bM后，可以进行迭代计算，也可以求解三次方程的体积根。

建议大家自编程序进行计算。所得的结果列于下表： T/K 组成 RKS方程常数 400 y1?y2?0.5 组分（1）：a＝0.7568 b＝5.346×10-5 组分（2）：a＝1.007 b＝6.565×10-5 混合物a＝0.8774 b＝5.956×10-5 p/MPa 1 2 1442.3 3 877.0 4 585.5 5 399.3 Vcal/（cm3·mol-1） 3114.0 2－24．试用下列方法计算由30％（摩尔％）的氮（1）和70％正丁烷（2）所组成的二元混合物，在462K、69×105Pa下的摩尔体积。 （1）使用Pitzer三参数压缩因子关联式 （2）使用RK方程，其中参数项为：

bi?0.086640RTci

pci2.50.427480R2Tcijaij?pcij

?6?6（3）使用三项维里方程，维里系数实验值为B11?14?10，B22??265?10，

?9?93?1B12??9.5?10?6，（B的单位为m?mol）。C111?1.3?10，C222?3.025?10，

C112?4.95?10?9，C122?7.27?10?9（C的单位为m6?mol?2）。已知氮及正丁烷的临界

参数和偏心因子为

N2 Tc＝126.10K，pc＝3.394MPa，?＝0.040 nC4H10 Tc＝425.12K，pc＝3.796MPa，?＝0.199

解：（1）根据Kay规则求出混合物的虚拟临界参数

Tpc??yiTci?0.3?126.10?0.7?425.12?335.41K

ippc??yipci?0.3?3.394?0.7?3.796?3.675MPa

i???yi?i?0.3?0.040?0.7?0.199?0.1513

i虚拟对比条件为：Tpr?T462??1.377 Tpc335.41ppr?p6.9??1.877 ppc3.675?0?查图2－9和2－10得：Z则： Z?Z(0)?0.77,Z?1??0.19

??Z(1)?0.77?0.1513?0.19?0.7987

V?ZRT0.7987?8.314?462??4.446?10?4m3?mol?1 6p6.9?10??（2） RK方程

Z?1a?h??? 1.5?1?hbRT?1?h?h?bbp ?VZRT2a?y12a11?2y1y2a12?y2a22

b?y1b1?y2b2

.50.42748R2Tc212 a12?pc12组元 Tc/K pc?10?5/PaVc?106/(m3?mol-1)Zc 90.1 255 158.5 a/(Pa?m6?K0.5?mol-2)b/(m3?mol-1) 11 22 12 126.10 425.12 231.53 33.94 37.96 34.37 0.292 1.555 0.274 29.01 0.283 7.012 2.676×10-5 8.067×10-5 aM?0.32?1.555?2?0.3?0.7?7.012?0.72?29.01?17.30Pa?m6?K0.5?mol-2 bM?0.3?2.676?10?5?0.7?8.067?10?5?6.450?10?5m3?mol-1

Z??1a?1?hbRT1.5117.30?h?????1.5?1?h?1?h6.450?10?5?8.314??462??h???1?h??1?h??3.25??1?h1?h??

bbp6.450?10?5?69?1050.1159h????

VZRTZ?8.314?462Z进行试差迭代得：Z?0.746, h=0.156

V?ZRT0.746?8.314?462?43?1??4.15?10m?mol 6p6.9?10??（3） 三项的维里方程为：

Z?pVBC?1??2 RTVV