第二章 流体的p-V-T关系
0.42748??8.314???304.19?a?T??a???T??0.42748RT/pc???T?＝?1.08967.382?10＝0.40335Pa?m3?mol?12222cb?0.08664RTc/pc＝0.08664?8.314?304.19＝2.968?10?5m3?mol?1 67.382?10p?RTa?T?8.314?273.150.40335?＝－V?bV?V?b?550.1?10?6－29.68?10?6550.1?10?6??550.1＋29.68??10?63.090?3.099?100％＝?0.291％
3.090＝3.099?106Pa＝3.099MPa误差％＝
比较几种方程的计算结果，可见，van der Waals方程的计算误差最大，RKS方程的计算精度最好。RK方程的计算精度还可以。
2－18．一个体积为0.3m3的封闭储槽内贮乙烷，温度为290K、压力为25×105Pa，若将乙烷加热到479K，试估算压力将变为多少？
解：乙烷的临界参数和偏心因子为：Tc＝305.32K，pc＝4.872MPa，?＝0.099 因此：Tr1?T1/Tc?290/305.32?0.95 pr1?p1/pc?2.5/48.72?0.513 故使用图2-11，应该使用普遍化第二维里系数计算
B(0)?0.083?0.422/Tr1.6?0.083?0.422?0.95?1.6??0.375
B(1)?0.139?0.172/Tr4.2?0.139?0.172?0.95?4.2??0.074
Z?1?Bp?1?B?0???B?1?RT?p?????T?rr?0.513????1??0.375?0.099?0.074??0.7935 ?0.95?V?ZRT0.7935?8.314?290?53?1??76.5?10m?mol 5p25?10??n?V总0.3??392.2?mol? ?5V76.5?10加热后，采用RK方程进行计算。
其中：T=479K，摩尔体积仍然为V?76.5?10m?mol，首先计算：
?53?1a?0.42748RT22.5c0.42748??8.314???305.32?60.5-2/pc＝＝9.879Pa?m?K?mol 64.872?1022.5?? 11
第二章 流体的p-V-T关系
b?0.08664RTc/pc＝代入RK方程：
0.08664?8.314?305.32?53?1＝4.514?10m?mol 64.872?10??RTa?0.5V?bTV(V?b)8.314?4799.879 ＝－?6?60.5?6?6765.0?10－45.14?10?479??765.0?10??765.0＋45.14??10p?＝4.804?106Pa＝4.804MPa2－19．如果希望将22.7kg的乙烯在294K时装入0.085m3的钢瓶中，问压力应为多少？ 解：从附录三查得乙烯的临界参数为：Tc＝282.34K，pc＝5.041MPa，?＝0.085
m22.7?103n???810.7?mol?
M28V总0.085??104.8?10?6?m3?mol?1? 摩尔体积V?n810.7采用RK方程进行计算。 首先计算：
a?0.42748RT22.5c0.42748??8.314???282.34?/pc＝＝7.851Pa?m6?K0.5?mol-2 65.041?1022.5??b?0.08664RTc/pc＝代入RK方程：
0.08664?8.314?282.34＝4.034?10?5m3?mol?1 65.041?10??RTa?0.5V?bTV(V?b)8.314?2947.851 ＝－104.8?10?6－40.34?10?6?294?0.5?104.8?10?6??104.8＋40.34??10?6p?＝7.817?106Pa＝7.817MPa2-20（由于较简单省略了，忽略不计了）
2－21．用Pitzer的普遍化关系式计算甲烷在323.16K时产生的压力。已知甲烷的摩尔体积
3?1m?mol为1.25×10－4，压力的实验值为1.875×107Pa。
解：从附录三查得甲烷的临界参数为：Tc＝190.56K，pc＝4.599MPa，?＝0.011
Tr?T/Tc?323.16/190.56?1.696；但是pr不能直接计算，需要试差计算 p?ZRT/V?Z?8.314?323.16/1.25?10?4?2.149?107Z
?? 12
第二章 流体的p-V-T关系
并且 p?pc?pr?4.599?106pr
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（3） 将查图得到的Z?0??0?和Z?1?；
和Z?1?值代入（B）式求得Z值；
（4） 比较Z的计算值与实验值，如果相差较大，则代入（A）式重新计算，直到迭代收
敛。
依据上述迭代结果为：pr＝4.06时，Z＝0.877
则：p?pc?pr?4.599?106pr＝4.599?106?4.06＝1.867?107Pa 误差： ?1.875?1.867??107/1.875?107?0.43%
2－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? 13
第二章 流体的p-V-T关系
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 14
第二章 流体的p-V-T关系
混合物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
（3）使用三项维里方程，维里系数实验值为B11?14?10?6，B22??265?10?6，
3?1（B的单位为m?mol）。CB12??9.5?10?6，?1.3?10?9，C222?3.025?10?9，111。已知氮及正丁烷的临界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
15