C2.2第二节状态方程

1、 从从2.1可知，流体的可知，流体的p-V-T数据是多么的重要。数据是多么的重要。 毫无疑问，流体的毫无疑问，流体的p-V-T数据是可以直接测定的。现已积累了大量数据是可以直接测定的。现已积累了大量水、空气、氨和氟里昂等的水、空气、氨和氟里昂等的 p-V-T数据。数据。 但测定数据费时耗资，而且测定但测定数据费时耗资，而且测定所有流体所有流体的的所有所有p-V-T数据数据是不现是不现实的，特别是高温高压下的实的，特别是高温高压下的p-V-T数据测定不但技术上相当困难而数据测定不但技术上相当困难而且非常危险。且非常危险。 仅从有限的仅从有限的p-V-T测定数据不可能全面地了解流体的测定数据不可能

2、全面地了解流体的p-V-T行为。行为。 流体的流体的p-V-T数据更重要作用是通过它推算数据更重要作用是通过它推算内能、焓、熵等信息。内能、焓、熵等信息。 需要需要P-V-T关系的解析形式！关系的解析形式！状态方程！状态方程！2.2 气体的状态方程2.2.1 状态方程（状态方程（EOS）的定义）的定义2.2.2 理想气体的状态方程理想气体的状态方程2.2.3 真实气体的状态方程真实气体的状态方程 van der Waals范德华状态方程范德华状态方程 Redlich-Kwong状态方程状态方程 Soave- Redlich-Kwong状态方程状态

3、方程 Peng-Robinson状态方程状态方程 Virial（维里）状态方程（维里）状态方程2.2.4 状态方程的小结状态方程的小结2.2.1 状态方程（状态方程（EOS）的定义）的定义n EOS是计算热力学性质是计算热力学性质最重要的模型之一最重要的模型之一。 . EOS是物质是物质P-V-T关系的解析式即关系的解析式即用一个用一个EOS即可精确地代表相当广泛范即可精确地代表相当广泛范围内的围内的P、V、T实验数据，借此可精确地计算所需的实验数据，借此可精确地计算所需的P、V、T数据。数据。2. 用用EOS可计算不能直接从实验测定的其它热力学性质可计算不能直

4、接从实验测定的其它热力学性质（H，S，G）数据。）数据。3. 用用EOS可进行相平衡和化学反应平衡计算可进行相平衡和化学反应平衡计算n 定义：定义：根据相律，纯物质在单相区的自由度根据相律，纯物质在单相区的自由度F=C-P+2为为2，因此因此P,V,T中任意两中任意两个指定后，状态就确定了。个指定后，状态就确定了。 P=f (V,T)或或 V=f (P,T)或或T=f (V, P) f(P,V,T) =0 状态方程状态方程 Equation of State(EOS) 目前已有目前已有150多种多种EOS。但没有一个。但没有一个EOS能描述在工程应用范围能描述在工程应用范围内内任何气体任何气体

5、的行为。的行为。 状态方程包含的状态方程包含的规律愈多规律愈多，方程就，方程就愈可靠愈可靠；准确性越高，范围准确性越高，范围越广，模型越有价值。越广，模型越有价值。 状态方程的准确度和方程型式的简单性是一对矛盾。状态方程的准确度和方程型式的简单性是一对矛盾。 建立建立EOS的方法：或以理论法为主、或以经验法为主。实际应的方法：或以理论法为主、或以经验法为主。实际应用以用以半经验半理论半经验半理论和和纯经验纯经验的的EOS为主。为主。 我们介绍各种我们介绍各种EOS的特点和应用范围，并不要求建立。的特点和应用范围，并不要求建立。2.2.2 理想气体的状态方程理想气体的状态方程n理想气体：理想气体

6、：假定分子的大小如同几何点一样，气体分子体积可忽假定分子的大小如同几何点一样，气体分子体积可忽略；分子间不存在相互作用力略；分子间不存在相互作用力n理想气体状态方程：理想气体状态方程： PV=RT; Z=PV/RT=1理想气体理想气体EOS是是f(P,V,T) =0 中中最简单的一种形式。最简单的一种形式。n应用状态方程时要注意：应用状态方程时要注意： (1)式中式中p是绝对压力，是绝对压力，V是摩尔体积，方程式中的是摩尔体积，方程式中的T是绝对温度是绝对温度 (2)通用气体常数通用气体常数R的单位必须和的单位必须和P,V,T的单位相适应。的单位相适应。 理想气体理想气体EOS只适合压力非常低

7、的气体只适合压力非常低的气体,不适合真实气体。不适合真实气体。它可以它可以 作作为衡量真实气体状态方程是否正确的标准之一，为衡量真实气体状态方程是否正确的标准之一， 当压力趋近于当压力趋近于 0或者体积趋于无穷大时，任何真实气体状态方程都或者体积趋于无穷大时，任何真实气体状态方程都应还原为理想气体方程。应还原为理想气体方程。理想气体理想气体EOS只适合压力非常低的气体只适合压力非常低的气体,不适合真实气体。不适合真实气体。例例:将将1kmol甲烷压缩储存于容积为甲烷压缩储存于容积为0.125m3,温度为温度为323.16K的钢瓶内的钢瓶内.问此时甲烷产生的压力多大问此时甲烷产生的压力多大?其实

8、验其实验值为值为1.875x107Pa.解解:%67.14%10010875.110875.1150.210150.2125.016.323314.81000777 ）（误误差差PaVnRTP误差高达误差高达14.67%! 2.2.3 真实气体的状态方程 真实气体分子真实气体分子有大小、分子间有相互作用力有大小、分子间有相互作用力是造成气体非理想性的原因是造成气体非理想性的原因 真实气体对理想气体的真实气体对理想气体的偏离程度偏离程度可以用可以用压缩因子压缩因子Z Z来表达来表达RTPVZ )(理理想想气气体体VVZ (a)分子间吸引力促使分子间吸引力促使Z1 。 (c)吸引力和排斥力的平衡暗

9、指吸引力和排斥力的平衡暗指Z=1 。（注意理想气体是。（注意理想气体是一个特例，既没有吸引力，也没有排斥力）。一个特例，既没有吸引力，也没有排斥力）。更直观的表达为：更直观的表达为： 真实气体偏离理想行为，真实气体偏离理想行为，理想气体状态方程理想气体状态方程不能描述真实气体的不能描述真实气体的状态，因此出现了：状态，因此出现了： van der Waals（ vdW范德华）状态方程范德华）状态方程 Redlich-Kwong状态方程状态方程 Soave-Redlich-Kwong状态方程状态方程 Peng-Robinson状态方程状态方程 Virial（维里）状态方程（维里）状态方程立方型立

10、方型多参数高次型多参数高次型2.2.3 真实气体的状态方程 范德华方程范德华方程理想气体理想气体体积修正项体积修正项b为斥力参数为斥力参数分子间力的修正项分子间力的修正项a为引力参数为引力参数。2VabVRTP VRTP van der Waals 方程是最早提出的描述实际流体体系的立方型状态方程。方程是最早提出的描述实际流体体系的立方型状态方程。由于分子相互吸引力存在，分子撞击器壁的力减小由于分子相互吸引力存在，分子撞击器壁的力减小，造成压力减小。压力减小的数值与撞击器壁的分，造成压力减小。压力减小的数值与撞击器壁的分子成正比；与吸引其分子数成正比，即与气体比容子成正比；与吸引

11、其分子数成正比，即与气体比容的平方成反比。的平方成反比。分子本身占有体积，分子自由活动分子本身占有体积，分子自由活动空间减小，由空间减小，由V变成变成V-b。 1873年荷兰莱顿大学年荷兰莱顿大学范德华范德华，在，在其博士论文其博士论文“关于气态和液态关于气态和液态的连续性的连续性”中提出中提出vdW EOS ； 不仅修正了理想气体不仅修正了理想气体EOS，而且，而且论证了气液态混合物是以连续方论证了气液态混合物是以连续方式互相转化的，是第一个式互相转化的，是第一个同时能同时能计算汽液两相和临界点计算汽液两相和临界点的方程。的方程。 1910年年获诺贝尔化学奖。获诺贝尔化学奖。1）范德华方程常

12、数）范德华方程常数a、b的确定的确定对于对于Van der Waals方程方程应用临界等温的数学特征，即应用临界等温的数学特征，即2RTaPVbV （ 1 1）0)Vp(,0)Vp(ccT22T （2）v临界等温线在临界等温线在C点点的斜率等于零的斜率等于零v临界等温线在临界等温线在C点点的曲率等于零的曲率等于零对对(1)式求关于摩尔体积式求关于摩尔体积V的一阶和二阶偏导数的一阶和二阶偏导数 ，即可得，即可得06)(2)(02)()(432232 cccTccTVabVRTVpVabVRTVpcc（3）（4）联立求解方程（联立求解方程（3）和（）和（4）,得得cccVRTaVb)89( 31

13、(5)将方程（将方程（1）用于临界点，即与）用于临界点，即与(5)式联立，即得式联立，即得2ccccVabVRTP ccPTRa226427 ccPRTb81 vdW常数常用形常数常用形式式将（将（6）式与（）式与（5）联立，即得）联立，即得由于由于VC的实验值的实验值误差大误差大，a，b要要用用Pc,Tc来表示来表示（6）（7）根据气体压缩因子的定义：根据气体压缩因子的定义：330.37588CCCCCCCCp VRTVZRTVRT1 van der Waals方程方程给出了一个固定的压给出了一个固定的压缩因子缩因子Zc ，即，即Zc=0.375。而而多数流多数流体的体的Zc实实验值在验值在

14、0.230.29之间之间，明显低于明显低于 van der Waals 方程方程的的Zc值。可见值。可见van der Waals 方程计算准确性不高。方程计算准确性不高。2 VdW EOS的缺点的缺点 两项修正项过于简单，准确度低，不能在任两项修正项过于简单，准确度低，不能在任何情况下都能精确描述真实气体的何情况下都能精确描述真实气体的P-V-T关系。实际应用少。关系。实际应用少。3 VdW EOS的改进的改进改进形式为改进形式为RK ；SRK ；PR 状态方程状态方程但改进形式均以但改进形式均以vdW状态方程为基础状态方程为基础 改变了方程的引力项改变了方程的引力项Patt，以使得计算的，

15、以使得计算的V减小，试图改进方程减小，试图改进方程计算计算P-V-T的准确性；的准确性； RK方程明显优于方程明显优于vdW方程，方程，是真正实用的是真正实用的EOS。2VabVRTP VdW方程方程: bVVTabVRTP 2/ Redlich-Kwong状态方程状态方程RK方程方程RK方程常数的获取ccPTRa5.2242748.0 ccPRTb08664.0 v用同于用同于vdW方程的方法得到常数方程的方法得到常数a,b值，即临界等温线值，即临界等温线在临界点的条件得到：在临界点的条件得到：R-K 方程中常数值不同于范德华方程中的方程中常数值不同于范德华方程中的a、b值，

16、不能将二者混淆。值，不能将二者混淆。同理，我们可以得到同理，我们可以得到RK方程的方程的ZC10.3338CZ 与与van der Waals方程相比，方程相比，RK方程计算的方程计算的ZC与多数气体与多数气体ZC的实验的实验值更为接近，因此值更为接近，因此RK方程计算精度比方程计算精度比van der Waals方程有了明显的方程有了明显的改善。改善。RK方程的特点方程的特点 1、RK方程计算气相体积准确性有了很大提高。方程计算气相体积准确性有了很大提高。 2、RK方程方程计算液相体计算液相体积的准确性不够，积的准确性不够，不能同时用于汽、液两相。不能同时用于汽、液两相。 3、RK方程用于烃

17、类、氮、氢等方程用于烃类、氮、氢等非极性气体时非极性气体时，即使在几百大，即使在几百大气压精度都较高，误差仅在气压精度都较高，误差仅在2%左右左右；但对于氨、水蒸气等；但对于氨、水蒸气等极极性较强的气体则精度较差性较强的气体则精度较差，误差在，误差在10-20%。 4、 对对RK方程进行修正，但同时降低了方程进行修正，但同时降低了 RK的简便性和易算性。的简便性和易算性。成功的有成功的有Soave的修正式（的修正式（SRK）定义参数定义参数A和和B：22 .52 .50 .4 2 7 4 8rrPa pARTT0.08664rrPbpBRTTrcTTTrcPPPRK方程可以表示成压缩因子方程可

18、以表示成压缩因子Z的三次方表达式：的三次方表达式：322()0ZZABBZABPVZRT RTa TPVbVVb ()11 cccccR Ta TaTTPRTbP 220.427480.08664 220.510.481.5740.1761rTT Soave把把RK方方程中的常数程中的常数a看看作是作是T的函数的函数:偏心因偏心因子子 与与RK方程相比，方程相比，SRK方程可计算方程可计算极性物质极性物质，更主要的是可计算，更主要的是可计算饱和液体密度，使之能用于混合物的汽液平衡计算，故在工业饱和液体密度，使之能用于混合物的汽液平衡计算，故在工业上获得了上获得了广泛应用。广泛应用。2.2.3.

19、3 SRK 方程方程以偏心因子作为第三参数1.1.偏心因子的定义偏心因子的定义：纯物质的偏心因子是根据物质的蒸气纯物质的偏心因子是根据物质的蒸气压来定义的。实验发现，纯态流体对比饱和蒸气压的对数与对压来定义的。实验发现，纯态流体对比饱和蒸气压的对数与对比温度的倒数呈近似直线关系，即符合比温度的倒数呈近似直线关系，即符合ssrcppp 对比饱和蒸气压1lo g1srrpT对于不同的流体对于不同的流体, 其其 具有不同的值。但具有不同的值。但Pitzer发现，简单流发现，简单流体（氩、氪、氙）的所有蒸气压数据落在了同一条直线上，而体（氩、氪、氙）的所有蒸气压数据落在了同一条直线上，而且该直线通过，

20、且该直线通过，Tr0.7， 这一点，如下图这一点，如下图log1srp 1.01.8-1-1.2-1.8lgPrs1/Tr 1 2简单流体简单流体(Ar,Kr,Xe）非球形分子非球形分子1(正癸烷正癸烷)非球形分子非球形分子2对比蒸汽压与温度对比蒸汽压与温度的近似关系的近似关系000.1)lg(7.0 rTsrP 1/Tr=1.43即即Tr=0.7对于给定流体对比蒸气压曲线得位置，可以用在对于给定流体对比蒸气压曲线得位置，可以用在Tr=0.7的流体与的流体与氩氩、氪、氙（简单球形分子）氪、氙（简单球形分子）的的 值之差来表征。值之差来表征。lgsrPPitzer把这一差值定义

21、为把这一差值定义为 偏心因子偏心因子1955年，年，Pitzer在系统研究了蒸气的数据后正式提出了偏心因在系统研究了蒸气的数据后正式提出了偏心因子的概念。子的概念。0.70.7lg()lg() 1lg()rrssrrTsrTppp 简单流体某流体某流体显然，简单流体的偏心因子为零，而其它流体的偏心因子均大显然，简单流体的偏心因子为零，而其它流体的偏心因子均大于零于零(除过除过H2和和He)。偏心因子实际上代表了实际流体与简单流。偏心因子实际上代表了实际流体与简单流体的分子间相互作用的差异。体的分子间相互作用的差异。目前，各种纯流体的偏心因子可通过查表得到，见本教材附录。目前，各种纯流体的偏心因

22、子可通过查表得到，见本教材附录。SRK方程的特点：方程的特点： 在临界点同在临界点同RK，Zc=1/3（偏大）；（偏大）； 计算常数需要计算常数需要Tc，Pc和和 （比（比RK多）多），a是温度的函数；是温度的函数； 除了能计算气相体积之外，能用于表达蒸汽压（除了能计算气相体积之外，能用于表达蒸汽压（汽液平衡汽液平衡），），是一个适用于汽、液两相的是一个适用于汽、液两相的EOS，但计算液相体积误差较大；，但计算液相体积误差较大； 与与RK方程相比，方程相比，SRK方程大大提高了表达纯物质汽液平衡方程大大提高了表达纯物质汽液平衡的能力，使之能用于混合物的汽液平衡计算，故在工业上获得的能力，使之能

23、用于混合物的汽液平衡计算，故在工业上获得了广泛的应用。了广泛的应用。SRK方程可以表示成压缩因子方程可以表示成压缩因子Z的三次方表达式：的三次方表达式： 2220.42748rrPapATR TT0.08664rrPbpBRTT3220ZZABBZAB a TRTPVbVVbb Vb (12) 220.457240.07780cccccR Ta TaTTPRTbP 220.510.374641.542260.269921rTT Peng-Robinson（PR）方程）方程由于由于RK方程和方程和SRK方程均无法准确预测液相摩尔体积，且方程均无法准确预测液相摩尔体积，且ZC数值与

24、实验值相比明显偏大，因此数值与实验值相比明显偏大，因此Peng和和Robinson提出提出了他们的状态方程了他们的状态方程PR方程的特点：方程的特点： Zc=0.307，更接近于实际情况，更接近于实际情况，虽较真实情况仍有差别，但虽较真实情况仍有差别，但PR方程计算液相体积的准确度较方程计算液相体积的准确度较SRK确有了明显的改善；确有了明显的改善； 计算常数需要计算常数需要Tc，Pc和和 ，a是温度的函数；是温度的函数； 能同时适用于汽、液两相，工业中得到广泛应用；能同时适用于汽、液两相，工业中得到广泛应用； 在提供的计算软件在提供的计算软件Thermo-Pro中，用中，用PR作为状态方程模

25、型，作为状态方程模型，用于均相性质、纯物质饱和性质、混合物汽液平衡计算等。用于均相性质、纯物质饱和性质、混合物汽液平衡计算等。 PR方程预测方程预测液体摩尔体积液体摩尔体积的准确度较的准确度较SRK方程有明显的改善。方程有明显的改善。PR方程可以表示成压缩因子方程可以表示成压缩因子Z的三次方表达式：的三次方表达式：2220.45724rrPapATR TT0.07780rrPbpBRTT322231230ZB ZABBZABBB0)(23 PabVPaVPRTbV1) TTc2) T=Tc3) TTc仅有一个实根，对应于仅有一个实根，对应于液体的摩液体的摩尔体积尔体积VL或气体的摩尔体积或气体

26、的摩尔体积VV。当当PPc时，时，仅有仅有一个实根一个实根VL或或VV当当P=Pc时，三个重实根时，三个重实根 V=VcCPVV SL V x V SVP*TTc3)TTc三个不同实根，发生于两相区三个不同实根，发生于两相区V大大对应于饱和蒸汽摩尔体积对应于饱和蒸汽摩尔体积VSVV小小对应于饱和液体摩尔体积对应于饱和液体摩尔体积VSLV中中无物理意义。无物理意义。 bVVTabVRTP 2/11)已知已知T,V，如何求，如何求P？v显压型，直接计算，很方便。显压型，直接计算，很方便。v在计算时，一定要注意单位，在计算时，一定要注意单位，1atm=0.101325106Pa=0.101325MP

27、a2)已知已知P,T，如何求，如何求V？ 工程上最常用的情况，因为工程上最常用的情况，因为P，T易测。用各种易测。用各种迭代法迭代法求解求解。2.立方型状态方程具体解题方法立方型状态方程具体解题方法以以RK方程为例：方程为例： 考虑到迭代法的收敛问题，需要改变方程的形式，以求方程的特解。考虑到迭代法的收敛问题，需要改变方程的形式，以求方程的特解。 11/ 2nnnna VbRTVbppTVVb 已知已知P,T，如何求，如何求V？直接迭代法直接迭代法 bVVTabVRTP 2/11/ 2a VbRTVbpTpV Vb变换成迭代式变换成迭代式(a)求气体的摩尔体积求气体的摩尔体积迭代初值迭代初值0

28、RTVp代入迭代式，得到的代入迭代式，得到的V值后再代到等式的右边，一直迭代到值后再代到等式的右边，一直迭代到V值的变值的变化很小，达到要求为止，譬如化很小，达到要求为止，譬如 。 4110nnVV 例例:将将1kmol甲烷压缩储存于容积为甲烷压缩储存于容积为0.125m3,温度为温度为323.16K的的钢瓶内钢瓶内.问此时甲烷产生的压力多大问此时甲烷产生的压力多大?其实验值为其实验值为1.875x107Pa.解解:1）理想气体）理想气体2） RK方程方程查附录得查附录得TC=190.6K，PC=4.600MPa=4.6x106Pa R=8.314x103 m3.Pa/kmol.K误差高达误差

29、高达14.67%! ccPTRa5.2242748.0 ccPRTb08664.0 误差仅为误差仅为1.216%! PabVVTabVRTP72/1108978.1 例例 异丁烷是取代氟利昂的环保制冷剂，用于冰箱、冷柜、冷饮机。异丁烷是取代氟利昂的环保制冷剂，用于冰箱、冷柜、冷饮机。现需要将现需要将1kmol、300K、0.3704MPa的异丁烷装入容器，的异丁烷装入容器，请问请问需设计多大的容器需设计多大的容器？试用理想气体方程、？试用理想气体方程、RK、SRK和和PR方程方程分别计算，并与实际值进行比较（实际值为分别计算，并与实际值进行比较（实际值为V=6.081m3/kmol ）kmol

30、mbkmolKmkPaa/08058.010648.31 .408314.808664.0/10725.210648.31 .408314.842768.03325 .06435 .22 ccPTRa5.2242748.0 ccPRTb08664.0 bVVTabVRTP 2/1 11/ 2nnnna VbRTVbPPTVVb 141/ 28.3143000.08058370.42.725100.08058370.43000.08058nnnnVVVV 14.2480.080586.8140.08058nnnnVVVV kmolmPRTV/734.64 .370300314.830 198.6

31、08058.0734.6734.608058.0734.6248.4814.61 V 146.608058.0198.6198.608058.0198.6248.4814.62 VkmolmVVV/140.6140.6141.6343 实验值为实验值为V=6.081m3/kmol,误差为：误差为：0.97%v ( 2 ) SRK方程方程7351.01.408300 rTkmolmbkmolmkPaa/08058.036481 .408314.808664.0/7 .16532259.136481 .408314.842768.032622 2259.17351.01176.1760.0176.

32、0574.148.0125.02 T bVVabVRTP bVPVbVabPRTV 1nnnna VbRTVbPPVVb 18.3143000.08058370.41653.70.08058370.40.08058nnnnVVVV 14.4650.080586.8140.08058nnnnVVVV kmolmPRTV/734.64.370300314.830 167.608058.0734.6734.608058.0734.6465.4814.61 V 109.608058.0167.6167.608058.0167.6465.4814.62 VkmolmVVV/101.6101.6102.6

33、343 实验值为实验值为V=6.081m3/kmol,误差为：误差为：0.38% （3）PR方程：方程： 相同方法可求得：相同方法可求得：V=6.0685m3/kmol 实验值为实验值为V=6.081m3/kmol,误差为：误差为：-0.16% （4）各种状态方程）各种状态方程误差比较误差比较：EOS方程方程误差误差理想气体方程理想气体方程10.74%RK方程方程0.99%SRK方程方程0.38%PR方程方程-0.16%注：实验测定的注：实验测定的误差为误差为0.5%左右，左右，因此因此SRK、 PR方程的方程的误差已小于实验误差，非常了不起！误差已小于实验误差，非常了不起！应用应用Excel

34、Excel和现成软件计算和现成软件计算EOSEOS 迭代法用手工计算来完成是相当繁琐的，可以自己编程序，或迭代法用手工计算来完成是相当繁琐的，可以自己编程序，或用诸如用诸如Mathcad或或Maple的软件包求解，当然必须给定初值或的软件包求解，当然必须给定初值或求解范围。求解范围。 当然最简单的是利用网上免费状态方程计算软件当然最简单的是利用网上免费状态方程计算软件（http:/www.cheng.cam.ac.uk/pjb10/thermo/pure.html）来求）来求解。解。 应用状态方程时要注意： (1)式中的式中的V是是摩尔体积摩尔体积! SI制单位：制单位：m3/mol (2)

35、式中的式中的T是是绝对温度绝对温度, K (3)式中式中P是是绝对压力，绝对压力，SI制单位：制单位：pa (4)式中的式中的n为为1 mol (5)通用气体常数通用气体常数R的的单位必须单位必须 和和P,V,T的单位相适应的单位相适应。 表表2-1 通用气体常数通用气体常数R值值(P.6) 注：注：R=8.314m3.pa/mol.K=8.314 J/mol.K建议各物理量均换算为建议各物理量均换算为SI制单位制单位，再代入方程式中进行计算。，再代入方程式中进行计算。 Virial (维里)方程 21VCVBRTPVZ 21PCPB Virial方程的形式方程的形式Onnes（

36、昂尼斯）提出（昂尼斯）提出 立方型状态方程属于半经验半理论的状态方程，而在状态方程立方型状态方程属于半经验半理论的状态方程，而在状态方程中有严格理论基础的当属中有严格理论基础的当属Virial方程。方程。0PV ，Z为压缩因子为压缩因子1PVZRT 当当 B、C（或（或B、C）称作）称作Virial系数，是有物理意义的。系数，是有物理意义的。 微观上，微观上，Virial系数系数反映了反映了分子间的相互作用，分子间的相互作用，第二维里系数第二维里系数B反映了两个反映了两个分子之间的相互作用；第三维里系数分子之间的相互作用；第三维里系数C反映了反映了三重三重分子的相互作用。分子的相互作用。 宏观上，宏观上，Virial系数系数仅是温度的函数。仅是温度的函数。 Virial 系数的获取系数的获取 ( 1 ) 由统计力学进