第二章 流体的P-V-T关系和状态方程

为什么要研究流体的P-V-T 关系？ 流体的P-V-T 关系可直接用于设计 如：1）一定T、P 下， ρ ？ Vm ? 2）管道直径的选取： 流量 3）储罐的承受压力：P 利用可测的热力学性质（T，P，V，CP）计算不可测的热力学性质（H，S，G，f，φ ，α ，γ ）

存在的问题 热力学最基本性质有两大类 P, V, T, Cp, x U, H, S, G 易测 难测 存在问题： 从易测的物性数据推算难测量的数据 易测 难测 存在问题： 测定所有流体的P-V-T数据费时耗资，且不现实； 有限的P-V-T数据，无法全面了解流体的P-V-T行为； 离散的P-V-T数据，不便于求导和积分，无法获得数据点以外的P-V-T和U, H, S, G数据。

本章内容 定性认识流体P-V-T 行为； 掌握描述流体P-V-T 关系的模型化方法，了解几种常见的状态方程； 掌握对比态原理和普遍化状态方程

2.1 纯流体的P-V-T关系 纯物质P-V-T 图 凝固时收缩 凝固时膨胀 图2-1(a) 纯物质的p-V-T相图 P P V V T T 气 临界点 液-汽 液 固 固-汽 汽 三相线 T V P 固 液 液-汽 汽 气 临界点 三相线 固-汽 P V T 凝固时收缩 凝固时膨胀 图2-1(a) 纯物质的p-V-T相图

图2-1(b) 纯物质的p-V-T相图

p-V 图

P-V图的特征 单相区（V, G, L） 两相共存区（V/L） 饱和线（饱和液体线、饱和气体线） 过热蒸汽 过冷液体 等温线（T=Tc、T>Tc、T<Tc） 临界等温的数学特征 超临界流体（T>Tc和P>Pc）

Tc是加压使气体液化所允许的最高温度 1）当T>Tc时，无论施加多大的压力都不能使该气体液化； 所谓“永久气体”也是由于室温比其Tc高所致。1908年7月10日，荷兰物理学家昂纳斯(H. K. Onnes)实现了1.15K的低温，将氦液化了，从而消除了最后一种“永久气体”。 3）当T=Tc时，只需施加比Pc稍大的压力就可使该气体液化。

p-T 图

P-T图的特征 单相区 两相平衡线（饱和曲线） 汽化曲线、熔化曲线、升华曲线 三相点和临界点（Tc, Pc, Vc） 等容线 临界等容线V=Vc、V>Vc、V<Vc

例2-1： 1）过热蒸汽等温冷凝为过冷液体； 2）过冷液体等压加热成过热蒸汽； 3）饱和蒸汽可逆绝热膨胀； 4）饱和液体恒容加热； 5）在临界点进行的恒温膨胀。

纯物质PVT关系的应用 超临界技术 液化气体成分的选择 制冷剂的选择 气体液化和低温技术

临界点及超临界流体 临界点：气液两相共存的最高温度或最高压力； 超临界流体（Supercritical Fluid, SCF) 在T>Tc和P>Pc区域内，气体、液体变得不可区分，形成一种特殊状态的流体，称为超临界流体。 超临界流体的特征 1）具有液体一样的密度、溶解能力和传热系数，具有气体一样的低粘度和高扩散系数； 2）物质的溶解度对T、P的变化很敏感，特别是在临界状态附近，T、P的微小变化会导致溶质的溶解度发生几个数量级的突变。超临界流体正式利用这一特征，通过对T、P的调控来进行物质的分离。

0.7

超临界萃取技术的工业应用 现研究较多的超临界流体包括：CO2、H2O、NH3、甲醇、乙醇、戊烷、乙烷和乙烯等； 1）临界条件温和：Tc=31℃，Pc=7.4MPa。萃取温度低（30~50 ℃ ），能保留产品的天然有效活性。溶解能力强； 2）惰性（不污染产品）、廉价易得、选择性良好以及产物易从反应混合物中分离出来。

大规模超临界流体萃取的兴起于： 用超临界CO2成功地从咖啡中提取咖啡因； 用超临界戊烷流体从石油中分离重油组分。 现在用于提取油脂、香精、色素、药物、酶的有效成分。特别是中药的提取方兴未艾： 从红豆杉树皮叶中获得的紫杉醇是抗癌药物； 从鱼内脏和骨头提取的鱼油； 从银杏叶中提取的银杏黄酮； 从蛋黄中提取的卵磷脂等治疗心脑血管疾病。 可见，超临界流体在化学工业的应用将会越来越广泛。

临界点数据的应用——液化气成分的选择 例2-2：液化石油气的成分为何是丙烷、丁烷和少量的戊烷，而不是甲烷、乙烷和正己烷？

2.2 流体的状态方程 纯物质的状态方程(Equation of State, EOS) 是描述流体p-V-T性质的关系式，即： f( p, T, V ) = 0 状态方程类型：立方型、多常数型、理论型； 混合物的状态方程从纯物质出发，通过引入混合规则，来计算混合物的热力学性质。

状态方程的应用主要表现在： 1.用一个状态方程即可精确地代表相当广泛范围内的p、V、T实验数据，借此可精确地计算所需的p、V、T数据。 2.用状态方程可计算不能直接从实验测定的其它热力学性质。 3.用状态方程可进行相平衡、化学反应平衡计算和进行化工过程能量分析。

2.2.1 理想气体状态方程 理想气体状态方程的形式为： 理想气体两个假设： 1）气体分子之间无相互作用力 2）气体分子本身不占有体积

掌握理想气体状态方程需明确的三个问题： 理想气体本身是假设的，实际上是不存在的。但它是一切真实气体当P →0 时可以接近的极限，因而该方程可以用来判断真实气体状态方程的正确程度，即：真实气体状态方程在P →0 时，应变为：PV = RT 2. 低压下的气体（特别是难液化的N2，H2，CO，CH4，…），在工程设计中，在几十个大气压（几个MPa）下，仍可按理想气体状态方程计算P、V、T： 而对较易液化的气体，如NH3，CO2，C2H4（乙炔）等，在较低压力下，也不能用理想气体状态方程计算。 3. 理想气体状态常被作为真实流体的参考态或初值。

理想气体EOS只适合压力非常低的气体，不适合真实气体。 例2-3：将1kmol甲烷压缩储存于容积为0.125m3，温度为323.16K的钢瓶内，问此时甲烷产生的压力多大？（实验值1.875×107Pa) 解： 理想气体EOS只适合压力非常低的气体，不适合真实气体。

2.2.2 气体的非理想性 真实气体分子有大小、分子间有相互作用力是造成气体非理想性的原因。 真实气体对理想气体的偏离程度可以用压缩因子Z来表达： 分子间吸引力促使Z<1; 分子间排斥力使Z>1； 吸引力和排斥力的平衡暗指Z=1。

2.2.3 立方型状态方程 2.2.3.1 Van der Waals (vdW)范德华方程 vdW方程的优点： 压力修正项， 体积修正项， b为有效分子体积，斥力参数 vdW方程的优点： 1）1873年范德华在其著名的论文“关于气态和液态的连续性”中提出，是第一个有实用意义的状态方程。1910年曾获诺贝尔奖。 2）是第一个能同时计算汽、液两相和临界点的方程。

范德华方程常数a、b的确定 应用临界等温的数学特征得： 联立方程组解得a、b： 代入vdW方程，并应用于临界点，得：

vdW EOS的缺点： vdW常数常用表达形式： 两项修正项过于简单，准确度低，不能在任何情况下都能精确描述真实气体的PVT关系。 给出了一个固定的临界压缩因子Zc=0.375，与实际流体Zc在0.23~0.29之间有明显差别，故由vdW计算的Vc大于实际流体的Vc。

R-K 方程是1949 年由Redlich 和Kwong共同研究提出的，是目前公认的最准确的双参数EOS。 2.2.3.2 Redlich-Kwong (RK) 方程 vdW方程： R-K 方程是1949 年由Redlich 和Kwong共同研究提出的，是目前公认的最准确的双参数EOS。 对方程的引力项进行修正，以使计算的V减小，提高计算的准确性，是真正实用的EOS。

R-K 方程中的常数a,b 的求取 用同于vdW方程的方法得到常数a,b值， 即临界等温线在临界点的条件得到：

R-K 方程的特点 1.RK方程计算气相体积准确性有了很大提高。 2.RK方程能较成功地用于气相P-V-T的计算，但计算 液相体积的准确性不够，不能同时用于汽、液两相。 3.RK方程用于烃类、氮、氢等非极性气体时，即使在几百大气压精度都较高，误差仅在2%左右；但对于氨、水蒸气等极性较强的气体则精度较差，误差在10-20%。 4.对RK方程进行修正，但同时降低了RK的简便性和易算性。成功的有Soave的修正式（SRK）

2.2.3.3 Soave-Redlich-Kwong (SRK) 方程 与RK方程相比，SRK方程可计算极性物质，更主要的是可计算饱和液体密度，使之能用于混合物的汽液平衡计算，故在工业上获得了广泛应用。

2.2.3.4 Peng-Robinson (PR) 方程 PR方程预测液体摩尔体积的准确度较SRK方程有明显改善，而且也可用于极性物质。 能同时适用于汽、液两相，在工程相平衡计算中广泛应用。

2.2.3.5 立方型状态方程的通用形式 前面介绍的都是立方型状态方程，其中 里程碑式的状态方程：vdW方程，RK方程； 在工程上有广泛应用的状态方程：SRK方程，PR方程； 其通用形式为：

2.2.3.6 立方型状态方程根的求取 以vdW方程为例： 当T > Tc 时，立方型状态方程有一个实根，它是气体容积。 当T<Tc时，高压下立方型状态方程有一个实根，它是液体容积。低压存在三个不同实根，最大的V值是蒸气容积，最小的V值是液体容积，中间的根无物理意义

2.2.3.7 立方型状态方程的解题方法 RK方程： 1）已知T，V，求P； 显压型，直接计算 2）已知P，T，求V； 工程中常见情况，迭代求解。 3）已知P，V，求T。 迭代求解 直接迭代法 牛顿迭代法

已知P，T，求V 蒸汽的摩尔体积 得: 方程两边乘以 初值取

液体的摩尔体积 将方程写成三次展开式 初值取

牛顿迭代法

已知P，V，求T—牛顿迭代法

例2-3：将1kmol甲烷压缩储存于容积为0. 125m3，温度为323. 16K的钢瓶内，问此时甲烷产生的压力多大？（实验值1 例2-3：将1kmol甲烷压缩储存于容积为0.125m3，温度为323.16K的钢瓶内，问此时甲烷产生的压力多大？（实验值1.875×107Pa) 解：1) 理想气体 误差高达14.67%!!! 2) RK方程 查附录三得：Tc=190.6K，Pc=4.600MPa=4.6×106Pa R=8.314×103 m3.Pa/kmol.K 误差仅为1.216%!!!

应用Excel和现成软件计算 迭代法用手工计算来完成是相当繁琐的，可以自己编程序，或用诸如Mathcad或Maple的软件包求解，当然必须给定初值或求解范围。 推荐一个比较简单的方法—应用Excel的“单变量求解”工具，它将牛顿迭代法固化Excel中，能大大简化使用者的手工计算量； 当然最简单的是利用网上免费状态方程计算软件（http://www.cheng.cam.ac.uk/~pjb10/thermo/pure.html)来求解。 Thermalcal热力学性质计算程序

应用Excel求解 以RK方程为例： 将其变为目标函数 工具 单变量求解： 目标单元格：目标函数 目标值： 0 可变单元格：需要求解的参数

应用状态方程时要注意： （1）式中的V是摩尔体积！SI制单位：[m3/mol] （2）式中的T是绝对温度，[K] （3）式中P是绝对压力，SI制单位： [pa] （4）式中的n为1mol （5）通用气体常数R的单位必须和P,V,T的单位相适应。 建议各物理量均换算为SI制单位，再代入方程式中进行计算。

例2-4：某企业要求以气体的形式存储10℃、1atm的丙烷35000kg。有两种方案争论不休： 你作为企业的工程师将采用何种方案，请比较这两种方案各自的优缺点。必要时可采用定量分析的方法。

2.2.4 Virial（维里）方程 Virial Equation，1907，荷兰Leiden 大学，Onness 令： 当P→0时，真实气体趋近于理想气体，PV=a,a=RT

Virial 系数 物理意义： B, B'：第二维里系数，它表示对一定量的真实气体，两个分子间的作用所引起的真实气体与理想气体的偏差。 C,C'：第三维里系数，它表示对一定量的真实气体，三个分子间的作用所引起的真实气体与理想气体的偏差。 维里系数= f （物质，温度） 二种形式的Virial方程是等价的，其系数之间也有 相互关系。

Virial 截断式 维里系数的获得： 1、由统计热力学计算，应用很少； 2、实验测定；精度较高 3、用普遍化关联式计算，简便； 1.两项维里截断式 2.三项维里截断式

⑴用于气相 PVT 性质的计算，对液相不适用； ⑵不适合高压； 最常用的维里截断式： ⑴用于气相 PVT 性质的计算，对液相不适用； ⑵不适合高压； ⑶理论意义大于实际应用价值，已超出PVT的应用，能描述气体的粘度、声速和热容。 ⑷是其他多参数方程BWR方程、MH方程的基础。

2.2.5 多参数状态方程 BWR方程（1940年提出） 应用条件： 式中，ρ为密度； 是8个常数，由纯物质的p-V-T 数据和蒸汽压数据确定。

Martin-Hou方程 （1955年提出） “新的 MH 方程（33参数）…，其计算的准确度已达到或超过了33参数 BWR 方程。” 方程参数：式中， 皆为方程的参数，这些参数可从纯 物质临界参数及饱和蒸汽压曲线上的 一点数据求得。 MH-55方程中， MH-81型方程中， “新的 MH 方程（33参数）…，其计算的准确度已达到或超过了33参数 BWR 方程。” ——高校化学工程学报 第15卷 第2期 2001第107页

应用范围： MH-55型计算精度高，气相1%，液相<5%;并且常数易确定，仅需要临界点和常压下的数据；但对液相极性物质计算的误差较大，最大可达16%。 MH-81型状态方程能同时用于汽、液两相，方程准确度高，适用范围广，能用于包括非极性至强极性的物质（如NH3、H2O)，对量子气体H2、He等也可应用，在合成氨等工程设计中得到广泛使用。

状态方程的选用 真实流体Zc=0.23~0.29

2.3 对应态原理和普遍化关联式 状态方程存在的问题 真实气体状态方程都涉及到物性常数，方程的通用性受到了限制。在实际工作中，当研究的物质其热力学性质既没有足够的实验数据，又没有状态方程中固有的参数时，计算便会产生困难，因此十分需要研究能用于真实气体的普遍化方法。

2.3.1 对应态原理 气体的对应态原理 对比参数定义为： 真实气体 实验发现所有气体的临界压缩因子ZC 值相近，表明所有气体在临界状态具有与理想气体大致相同的偏差： 对多数非极性物质Zc≈0.27

以临界状态为起点，将温度、压力、体积表示为对比参数： 对应态原理（The principle of corresponding states)：所有的物质在相同的对比态（Tr,Pr)下，表现出相同的性质（Vr或Z近似相等），即：组成、结构、分子大小相近的物质有相近的性质。

2.3.2 两参数对应态原理 范德华第一个提出两参数对应态原理： 在实际中，人们又发现对于大多数物质（约60%）的临界压缩因子Zc 在0.26-0.29 之间，因而一般取Zc=0.27。把临界压缩因子看作常数： 意义：对应态原理是一种特别的状态方程，也是预测流体性质最有效的方法之一 。

两参数普遍化压缩因子图（低压段）

两参数普遍化压缩因子图（中压段）

两参数普遍化压缩因子图（高压段）

表2-3-1 不同气体的Zc Ar 0.291 N2 0.290 NH3 0.242 CH4 0.288 乙炔 0.271 He 0.301 乙炔 0.271 He 0.301 H2 0.305 H2S 0.284 乙烷 0.285 甲醇 0.224 Ne 0.311 O2 0.288 HCl 0.249 丙烷 0.281 乙醇 0.248 Ke 0.288 Cl2 0.275 H2O 0.229 正丁烷0.274 乙醛 0.22 Xe 0.286 Br2 0.270 SO2 0.268 异丁烷0.283 丙酮 0.232 CO 0.295 NO 0.25 乙烯 0.276 醋酸 0.2 CO2 0.274 NO2 0.48 丙烯 0.275 苯 0.271 丁烯 0.277 甲苯 0.264

2.3.3 三参数对应态原理 z = f (Tr,Pr,第三参数) 两参数对应态原理只能适合于简单的球形流体； 为提高对比态原理的精度，需引入第三参数； 第三参数的特性：最灵敏反映物质分子间相互作用力的物性参数，当分子间的作用力稍有不同，就有明显的变化。 z = f (Tr,Pr,第三参数) 1955 年，K.S.Pitzer 提出了以偏心因子作为第三因子的关系式： z = f (Tr,Pr,ω )。把压缩因子看作是对比温度、对比压力和偏心因子的函数。

偏心因子 在低压下，克-克方程表示为： 实验发现，纯态流体对比饱和蒸气压的对数 与对比温度Tr的倒数近似于直线关系。 简单流体(氩、氪、氙)作 lgprS～1/Tr图，其斜率相同，且通过点（Tr=0.7，lgprS=-1）;

表征了一般流体与简单流体分子间相互作用的差异。 非球形分子的直线都位于球形分子的下面，物质的极性越大，其偏离程度也越大。在Tr=0.7时， lgprs<-1。 表征了一般流体与简单流体分子间相互作用的差异。

三参数对应态原理 在相同的Tr和Pr下，具有相同ω值的所有流体都具有相同的压缩因子Z，而且它们偏离理想气体的程度都相同； 这比原始的两参数对应态原理有了很大的进步； 由此可知：气体偏离理想气体的行为不是单由T、P决定的，而是由对比温度、对比压力以及偏心因子共同决定的。 根据以上结论，Pitzer 提出了两个非常有用的普遍化关系式： 1、压缩因子的多项式表示的普遍化关系式（简称普压法） 2、两项维里系数表示的普遍化关系式（简称为普维法）

2.3.4 普遍化压缩因子图法 对于所有相同的流体，若处在相同的Tr和pr下，其压缩因子必定相等。 Z(0) ——简单流体的压缩因子 Z(0)和Z(1)的都是Tr和pr的复杂函数，计算时可直接查图 Pitzer普遍化关系式对于非极性或弱极性气体能够提供可靠结果，误差<3%，对强极性气体则误差达5~10%；而对于缔合气体和量子气体，误差较大，使用时应当特别注意。

Z（0）的普遍化关系

Z（1）的普遍化关系

例2-5：计算 1Kmol乙烷在382K、21.5MPa时的体积。 解： 若已知T,V, 如何求Z?

例2-3：将1kmol甲烷压缩储存于容积为0. 125m3，温度为323. 16K的钢瓶内，问此时甲烷产生的压力多大？（实验值1 例2-3：将1kmol甲烷压缩储存于容积为0.125m3，温度为323.16K的钢瓶内，问此时甲烷产生的压力多大？（实验值1.875×107Pa) 解：1、查附录三得 Tc=191.6 K , Pc=4.60MPa, 所以采用普遍化压缩因子方法。

误差仅为0.535!!!

2.3.5 普遍化第二virial系数法 普遍化第二维里系数 Pitzer等提出了如下关联式 Tsonopoulos提出的关联式 ：

例2-6：试用下列三种方法计算510K、2.5MPa下正丁烷的摩尔体积。已知实验值为1.4807m3.kmol-1。 (1) 用理想气体方程； (2) 用普遍化压缩因子关联； (3) 用普遍化Virial系数关联。 解： 从附录1.1查得正丁烷的临界参数 （1）用理想气体方程

（2）用普遍化压缩因子关联 查压缩因子表，内插查得 （3）用普遍化Virial系数关联

此值和压缩因子关联结果相比，误差小于1%。

注意事项： 选用方程进行计算时，精度的大小对于工程技术人员来说也是一个很重要的指标； 三参数普遍化关系是能够很好的满足工程需要，一般对于非极性和弱极性物质，误差约3%，强极性物质为5-10%； 在这里要提醒大家的是，在工作中要计算PVT性质时，首先必须会查找手册，查出实验数据，只有实验数据才是最为可靠的； 如果确实找不到实验数据，就要进行计算，计算方法就是我们前面介绍的，但并不仅仅是这些，有些我们没有讲到的方法也是很有价值的； 在选取方程式计算时，一定要注意你所选取的方程是否适用于你所研究的范围，切不可没有原则的乱用。

2.4 真实气体混合物的P-V-T关系 2.4.1 混合规则 真实气体混合物的非理想性：一方面是由纯气体的非理想性，另一方面是由于混合作用所引起的非理想性。 真实气体混合物PVT性质的计算方法：一种是状态方程法，另一种是普遍化关系法。 研究思路 1）状态方程是针对纯物质提出的； 2）把混合物看成一个虚拟的纯物质，算出虚拟的特征参数，然后代入纯物质的状态方程计算混合物的性质； 3）混合物的虚拟特征参数是纯物质参数项和组成的函数，即混合规则： 4）混合规则是计算混合物性质中最关键的一步。 真实气体混合物的非理想性：一方面是由纯气体的非理想性，另一方面是由于混合作用所引起的非理想性。 真实气体混合物PVT 性质的计算方法：一种是状态方程法，另一种是普遍化关系法。 研究思路 1）状态方程是针对纯物质提出的； 2）把混合物看成一个虚拟的纯物质，算出虚拟的特征参数，然后代入纯物质的状态方程计算混合物的性质； 3）混合物的虚拟特征参数是纯物质参数项和组成的函数，即混合规则： 4）混合规则是计算混合物性质中最关键的一步。

2.4.2 虚拟临界参数法和Kay规则 该法是由 W.B.Kay 提出的，其主体思想是把混合物人为地看作是一种虚拟的纯物质，从而将纯物质对比态原理的计算方法用混合物上。 Kay 规则是最简单的混合规则，它将混合物的虚拟临界特征参数表示为： 适用范围： Prausnitz和Gunn提出了一个改进式

按纯组分气体PVT性质的计算方法进行计算。具体计算过程如下：普压法、普维法 虚拟对比参数 按纯组分气体PVT性质的计算方法进行计算。具体计算过程如下：普压法、普维法 缺点： 均没有涉及组分间的相互作用项； 对组分结构差异较大尤其是有极性和缔合作用的体系会产生较大的误差。

2.4.3 气体混合物的第二Virial系数 由统计力学可导出气体混合物的第二Virial系数与组成的关系为： 当i≠j时，Bij 为交叉第二Virial系数，且Bij = Bji 。i=j 时为纯组分i 的第二Virial系数。对二元混合物：

2.4.4 气体混合物的立方型状态方程 立方型状态方程（vdW、RK、SRK、PR）用于气体混合物，方程中参数a和b常采用以下混合规则： bi 是纯组分的参数，没有b的交叉项；aij 既包括纯组分参数(i=j)，也包括交叉项i≠j 。交叉项aij 按下式计算： Kij 为经验的二元相互作用参数，一般从混合物的实验数据拟合得到，对组分性质相近的混合物或近似计算可取 Kij = 0 。 Prausnitz 等人建议用下式计算交叉项aij

例2-7：试求二氧化碳和丙烷在311K、1.378MPa下等分子混合物的摩尔体积。 计算所需的临界参数及偏心因子数据列表如下 ij 11 304.2 7.376 0.0942 0.274 0.225 22 369.8 4.246 0.2030 0.281 0.145 12 335.4 5.470 0.1416 0.278 0.185 求出下列各第二Virial系数B值 ij 11 -0.324 -0.018 -0.1125 22 -0.474 -0.218 -0.3660 12 -0.394 -0.098 -0.2100

为了校验所用方法的适用性，计算虚拟临界温度和 虚拟临界压力。

该状态点落在图2-6的曲线之上，所以采用普遍化第二Virial系数法合适。

例2-8：试分别用下述计算二氧化碳和丙烷等分子混合物在444K、13.78MPa下的摩尔体积。 (1) 用RK方程； (2) 用普遍化压缩因子关系式。 已知实验值为

解： (1) 用RK方程 临界性质由例2-7给出，将上表中的有关数据代入公式，得如下结果 ij 11 6.440 22 18.214 12 11.072 － 计算混合常数

用RK方程的普遍化形式

混合物的摩尔体积 (2)用普遍化压缩因子关系式 虚拟临界参数同例2-7 查图

状态方程混合规则是气液两相均适用，但用于液相可靠性较差； 当计算混合物性质时，使用不同的状态方程，应采用不同的混合规则，计算不同的虚拟特征参数。 使用的状态方程类型 计算的虚拟特征参数 普遍化关系式 普遍化压缩因子图法 式（2-47）、（2-48）计算虚拟临界参数Tpc、Ppc 普遍化virial系数法 Virial方程 式（2-50）至（2-53）计算气体混合物的第二virial系数BM 立方型状态方程 式（2-55）至（2-58）计算气体混合物立方型状态方程参数am,bm

2.5 液体的P-V-T关系 2.5.1 饱和液体摩尔体积 Rackett方程 修正的Rackett方程 Vsl是饱和液体的摩尔容积；ZRA值可查阅文献，或用下式估算

例2-9： 计算异丁烷在273. 15K时饱和蒸汽压和饱和液体摩尔体积(实验值分别为152561Pa和100 例2-9： 计算异丁烷在273.15K时饱和蒸汽压和饱和液体摩尔体积(实验值分别为152561Pa和100.1cm3·mol-1)，并估计饱和汽相摩尔体积。 解：(a) 饱和蒸汽压 由Antoine方程计算。由附录查得Antoine方程常数 A＝6.5253，B＝1989.35，C＝-36.31 Ps ＝ 0.15347MPa＝153470Pa 与实验值的相对偏差为0.60％。

(b)饱和液相摩尔体积 用修正的Racket方程计算。 查得Tc＝408.10K，Pc＝3.646MPa, ω ＝0.176 与实验值的相对偏差为-0.17％。

( c )饱和汽相摩尔体积可以用Virial方程计算。

2.5.2 液体摩尔体积 Lyderson等提出了一个估算液体体积的普遍化方法。 液体的普遍化密度关系 液体对比密度：

2.6 状态方程的比较、选用和应用 每一个方程的特点和适用范围，见p.24表2-1. 作为工程师和设计人员得主要任务就是根据对精度的要求来选择状态方程。注意每一个方程的特点和适用范围。 状态方程精度粗略评价：多参数状态方程>立方型状态方程>两项截断维里方程>理想气体状态方程。 立方型状态方程中：PR>SRK>RK>vdW 1）若计算液体体积，则直接用修正的Rackett方程 既简单精度又高，不需要用立方型状态方程来计算；

2)若计算气体体积。 .SRK,PR是大多数流体的首选，无论压力、温度、极性如何，它们能基本满足计算简单。精度较高的要求，因此在工业上已广泛使用。 .对于个别流体或精度要求特别高的，则需要使用对应的专用状态方程或多参数状态方程，如对于CO2、H2S和N2首选BWRS方程； .在没有计算机软件又需要快速估算的情况下，精度要求非常低的可用理想气体状态方程，精度要求稍高可以使用普遍化方法。