《化工热力学》第二章流体的pVT关系课件.ppt

《《化工热力学》第二章流体的pVT关系课件.ppt》由会员分享，可在线阅读，更多相关《《化工热力学》第二章流体的pVT关系课件.ppt（140页珍藏版）》请在三一办公上搜索。

1、第二章 流体的p-V-T关系,2,2.1 纯物质的 p-V-T 关系 2.2 气体的状态方程 2.3 对比态原理及其应用 2.4 真实气体混合物的 p-V-T 关系 2.5 液体的 p-V-T 关系,1.教学基本要求（1）掌握纯物质的P-V-T性质；（2）掌握真实气体的状态方程式：Virial方程式、两常数状态方程式及其修正式、多常数状态方程式；（3）了解真实气体特性；2.教学重点和难点教学重点、难点是纯物质的P-V-T性质；真实气体的状态方程式,第二章 流体的p-V-T 关系,2,在众多的热力学性质中，流体的压力p，摩尔体积V 和温度T 是可以通过实验测量的；利用流体p-V-T 数据和热力学

2、基本关系式可计算不能直接由实验测得的其他性质，如焓H，内能U，熵S，自由焓G 等。因此，流体p-V-T 关系的研究是一项重要的基础工作。,2.1 纯物质的P-V-T 关系,2.1,图2-1,“固-汽”，“固-液”和“液-汽”分别表示固汽，固液和液汽平衡共存的两相区。,2.1,图2-1的三维曲面显示了处于平衡态下的纯物质 p-T-V关系，曲面以上或以下的空间是不平衡区。,三维曲面上“固”，“液”和“气（汽）”分别代表固体，液体和气体的单相区；,图2-1,图2-1的投影图,2.1,2.1,两相区在p-T图（图2-2）上的投影是三条相平衡曲线，升华线、熔化线和汽化线，三线的交点是三相点。,汽化线的另

3、一个端点是临界点C，它表示汽液两相能共存的最高压力和温度，即临界压力pc和临界温度Tc。,图2-2,气相,液相,固相,流体p-V-T关系，还可以用以T为参变量的p-V图表示，见图2-3。图中高于临界温度的等温线T1、T2，曲线平滑且不与相界线相交。,2.1,近于双曲线，即pV=RT=常数。（相当于理想气体的情形）。 小于临界温度的等温线T3 、T4有三个不同部分组成。,图2-3,沿着等温线，随着压力的提高，气体的体积逐渐缩小，变为这个温度下的饱和蒸汽。中间水平线段表示汽液平衡共存，在给定温度对应一个确定不变的压力，即该纯物质的饱和蒸汽压。汽液平衡混合物的组成从右端100%蒸汽变化到左端100%

4、液体。再压缩，液体的体积变化很小，显示了液体的不可压缩性。,2.1,曲线AC为饱和液体线，曲线BC为饱和蒸汽线。曲线ACB下面是两相区，其左、右面分别为液相区和气相区。,C,图2-3,气体,液体,等温线在两相区中的水平线段随着温度升高而缩短，最后在临界温度时缩成一点C。从图2-3上看出，临界等温线在临界点上的斜率和曲率都等于零。,2.1,数学上表示为:,(2-1),(2-2),图2-3,式(2-1)和式（2-2）提供了经典的临界点定义。Martin和侯虞钧在研究气体状态方程时发现，在临界点P对V的三阶和四阶导数也是零或是很小的数值。,2.1,图2-3,随着温度变化，饱和液体和饱和气体的密度迅速

5、改变，但两者改变的总和变化甚微。,2.1,Cailleter和Mathias注意到，当以饱和液体和饱和蒸汽密度的算术平均值对温度作图时，得一近似的直线，如图2-4所示。这结果称为直线直径定律，常用于临界密度的实验测定。,图2-4,2.2 气体的状态方程,2.2,用数学形式表示的体系 P-V-T 之间的关系式，称为状态方程（EOS） 。,(2-3),据相律可知，单相纯流体的 p、V、T 性质中任意两个确定后，体系的状态也就确定了。,2.2,状态方程的重要价值表现为：,用状态方程可精确地代表相当广泛范围内的 p-V-T 数据，从而大大减少实验测定的工作量；即 关联和推算的方法。根据有限的实验数据关

6、联出状态方程，有了状态方程就可推算更大温度、压力范围的数据；用状态方程可计算不能直接从实验测定的其他热力学性质； 用状态方程可进行相平衡计算，如计算饱和蒸气压、混合物汽液平衡、液液平衡等。尤其在计算高压汽液平衡时的简捷、准确、方便，是其他方法不能与之相比的。,2.2,总之，离散的p-V-T 实验数据点，经状态方程函数化后，在化工过程开发和设计中，不但可避免传统查图、查表的麻烦，而且借助电子计算机可实现准确快速的计算，极大提高工作效率。,2.2,一个优秀的状态方程应是形式简单，计算方便，适用范围广，计算不同热力学性质均有较高的准确度。但已发表的数百个状态方程中，能符合这些要求的为数不多。因此有关

7、状态方程的深入研究尚在继续进行。,2.2,状态方程的分类：,分为三类：非解析型状态方程 密度为三次方的立方型方程 非立方型方程（多常数Virial型方程）另一种分为两类：非解析型状态方程 解析型状态方程 本教材介绍重要而常用的解析型状态方程。,2.2.1,2.2.1 理想气体方程,理想气体的定义： 气体分子之间无相互作用力 。 气体分子本身不占有体积。,什么是理想气体？,2.2.1,2.2.1 理想气体方程,R的值：,PV=nRT 若n=1 则 PV=RT,(2-4),2.2.1,理想气体是极低压力和较高温度下各种真实气体的极限情况，实际上并不存在。理想气体方程除了在工程设计中可用作近似估算外

8、，更重要的是为判断真实气体状态方程的正确程度提供了一个标准。,任何真实气体状态方程都应还原为理想气体方程。使用状态方程时，应注意通用气体常数R的单位必须和 p, V, T的单位相适应。,2.2.1,随着化工生产技术的进展，系统压力的提高和温度的降低，真实气体与理想气体的偏离越来越显著。而且随着化工生产的大型化，对计算精度的要求越来越高，因此有必要来研究真实气体的状态方程。,要提出真实气体的状态方程，必须先了解真实气体的特性。,2.2.1,实验技术的不断改进，使人们掌握了真实气体偏离理想气体的许多事实与实验数据。而真实气体的特性正是通过实验，积累数据，并对这些数据进行分析和研究，才为人们所逐步认

9、识的。在此基础上，通过对实验数据的归纳、拟合，提出真实气体的状态方程。评价提出的状态方程的计算精度以及应用范围，也要靠实验数据来检验。因此实验在状态方程的开发与应用方面具有非常重要的作用。 下面介绍几种常用的真实气体状态方程。,2.2.2,2.2.2 立方型状态方程,所谓立方型状态方程是因为方程可展开为体积（或密度）的三次多项式。van der Waals 方程（1873年）是第一个适用真实气体的立方型方程，是对理想气体方程（2-4）的校正。,(2-5),方程中常数a -分子间存在相互作用的校正，常数b-分子本身占有体积的校正。,2.2.2,利用临界点 ， 的条件可以确定,虽然van der

10、Waals 方程准确度不高，无很大实用价值，但建立方程的理论和方法对以后立方型方程的发展产生了重大影响。目前工程上广泛采用的立方型方程基本上都是从van der Waals方程衍生出来的。,2.2.2,其中有代表性的有：,(1) Ridlich-Kwang方程（1949年）,Ridlich-Kwang方程简称RK方程，其形式为,(2-6),式中a ,b是方程常数，与流体的特性有关。,2.2.2,由纯物质临界性质计算,RK方程(Zc=1/3=0.333)适用非极性和弱极性化合物，计算准确度比van der Waals方程(Zc=3/8=0.375)有很大提高，但对多数强极性化合物有较大偏差。,(

11、2-7a),(2-7b),2.2.2,（2）Soave-Ridlich-Kwang方程（1972年）,(2-8),该方程简称 SRK方程。,2.2.2,Soave对RK方程的改进是将原方程中的常数a作为温度函数，即,(2-9a),(2-9d),(2-9c),(2-9b),式中为偏心因子。,2.2.2,SRK方程提高了对极性物质和量子化流体 p-V-T 计算的准确度。更主要的是Soave对方程的改进使方程可用于饱和液体密度的计算。在此基础上，用单一的SRK方程便可较精确地计算汽液平衡，拓宽了方程应用的领域。,2.2.2,（3）Peng-Robinson方程（1976年）,(2-10),其中,(2

12、-11a),(2-11b),(2-11c),(2-11d),2.2.2,Peng-Robinson(PR)方程中常数a仍是温度的函数，对体积表达的更精细的修正目的是为了提高方程计算Zc和液体密度的准确性。因此PR方程在计算饱和蒸汽压，饱和液相密度方面有更好的准确度。PR方程和SRK方程一样，是工程相平衡计算中最常用的方程之一。,2.2.2,(4) Patel-Teja方程（1982年）,(2-12),式中,(2-13a),(2-13b),(2-13c),(2-13d),2.2.2,a、b和c的计算方法如下：,而b是下式中最小的根,上述诸式中的 及F 是两个经验参数，由纯物质的饱和性质求得。,2

13、.2.2,Patel-Teja（PT）方程中引进了新的常数C，常数个数达到三个。常数多有利于提高方程的准确度，但也给方程的简明性和易算性带来损失。 用PT方程计算了一些极性和非极性纯物质的饱和气体和液体密度，其平均偏差分别为1.44%和2.94%（1070个数据点）。用PT方程计算轻烃、醇水等体系的汽液平衡也取得较好效果。,2.2.2,立方型方程形式简单，方程中一般只有两个常数，且常数可用纯物质临界性质和偏心因子计算。,表2-2 列出了立方型方程的三次展开式。摩尔体积V 和压缩因子Z 的三次展开式 (2-14)(2-21),2.2.2,2.2.2,在方程的应用中，准确地求取方程的体积根是一个重

14、要的环节。三次方程的求根方法有三次求根公式（表22）和数值计算法两大类。,该方法的思路是：,该方法计算过程稳定，缺点是耗时太多。,2.2.2,图2-5给出了数值算法之一的对分法求根的计算框图。,2.2.2,若将立方型方程改写成表2-3中的形式，则有可能仅仅通过手工计算就可求取方程的根，而不必借助计算机，迭代（试算）步骤是：,设初值Z（可用理想气体为初值，即取Z=1）；将Z值代入式（2-25）计算h，h;将h，h值代入表2-3中的状态方程计算Z 值；比较前后两次计算的Z值，若误差已达到允许范围，迭代结束；否则返回步骤2，再行计算。,请注意，该方法不能用于液相体积根的计算。,2.2.2,表2-3

15、立方型方程的另一形式,(2-22),(2-23),(2-24),(2-25),例题2-1 2-2,2.2.3,2.2.3 多常数状态方程,与简单的状态方程相比，多常数方程的优点是应用范围广，准确度高；缺点是形式复杂，计算难度和工作量都较大。由于电子计算机的日益普及，克服这些缺点已不成问题，因此多常数方程正越来越多地在工程计算中得到应用。,2.2.3,(1) Virial方程（1901年）,式中B(B)、C(C)、D(D)、.分别称为第二、第三、第四、维里系数。对一定的物质来说，这些系数仅仅是温度的函数。,(2-26),或者,(2-27),2.2.3,两组维里系数间的关系：（当维里方程取有限项时

16、的近似关系）,Virial方程是理论方程，具有坚实的理论基础，其系数有着确切的物理意义。,2.2.3,如第二Virial系数是考虑到二个分子碰撞或相互作用导致的与理想行为的偏差，第三Virial系数则是反映三个分子碰撞所导致的非理想行为。因为二个分子间的相互作用最普遍，而三分子相互作用、四分子相互作用等等的概率依次递减。,因此方程中第二Virial系数最重要， 在热力学性质计算和相平衡中都有应用。高次项对Z的贡献逐项迅速减小，只有当压力较高时，更高的Virial系数才变得重要。,2.2.3,方程式（2-26）和式（2-27）为无穷级数，如果以舍项形式出现时，方程就成为近似式。从工程实用上来讲，

17、在中、低压时，取方程式（2-26）和式（2-27）的二项或三项即可得合理的近似值。Virial方程的二项截断式如下：,式（2-28）可精确地表示低于临界温度、压力为1.5MPa左右的蒸汽的p V T 性质。,或,（2-28a）,（2-28b）,2.2.3,当压力超过适用范围而在5MPa以上时，需把Virial方程舍项成三项式，方能得到满意结果：,由于对第三Virial系数以后的Virial系数知道很少，且高于三项的Virial式使用起来不方便，所以对于更高的压力，通常都采用其他状态方程。,（2-29）,2.2.3,由于分子间相互作用十分复杂，至今建立的分子间位能函数仅对简单分子有较好精度，大多

18、数物质的Virial系数还需要通过实验测定。许多气体的第二Virial系数可从文献或有关手册中查到。但已测得的第三、第四Virial系数比较少。因而影响了Virial方程的应用。,随着分子相互作用的理论的进展，将可能从有关物质分子的基本性质精确计算维里系数，维里方程的应用还是有希望的。,2.2.3,(2) Martin - Hou方程（1955年）,该方程1955年由Martin和侯虞钧提出，简称MH方程。1959年对该方程作了进一步的改进，提高了其在较高密度区的精确度。 1981年侯虞钧等又将方程的适用范围扩展到液相区。,Martin - Hou方程的通式为,(2-32),式中,2.2.3,

19、81型MH方程的展开式为,式中A2，B2，C2，A3，B3，C3，A4，B4，B5及b（10个常数）皆为方程的常数，可从纯物质临界参数Tc、Pc、Vc及饱和蒸气压曲线上的一点数据（TS、pS）求得。,（2-33）,2.2.3,81型MH方程用于烃类和非烃类气体均令人十分满意，一般误差小于1%。对许多极性物质如NH3、H2O在较宽的温度范围和压力范围内，都可以得到精确的结果，对量子气体H2、He等也可应用，目前它已成功地用于合成氨的工艺计算。 1974年曾对合成氨系统的热量衡算和化学反应平衡中所需的气体混合物组成的热力学性质作了计算。并对年产三十万吨合成氨装置进行了热力学核算。认为该方程是满意的

20、。,2.2.3,（3）Benedict-Webb-Rubin方程（1940年）,该方程属于Virial型方程，简称BWR方程，在计算和关联轻烃及其混合物的液体和气体热力学性质时极有价值。与三参数对应态原理相结合，用来计算气体的热力学性质是BWR方程的一个重要应用。,其表达式为,(2-34),2.2.3,式中为密度；A0、B0、C0、a、b、c、和等8个常数由纯组分pVT 数据和蒸气压数据确定。作者在提出方程时，给出了12个轻组分的常数值，1967年Cooper和Goldfrank推荐了33种物质的常数值，1976年Holub又补充了8个组分的数据。不同来源的常数不能凑成一套使用，8个常数均有物

21、理量纲，使用时须采用一致的单位。 1972年Starling在BWR方程基础上提出11个常数的SHBWR方程：,（2-35）,2.2.3,修正式增加了D0、E0、d 三个常数，应用范围扩大，对比温度可以低到Tr=0.3，在比临界密度高达3倍的条件下也能用来计算气体的p-V-T关系。对轻烃气体、CO2、H2S和N2的广度性质作计算，误差范围在0.5%2.0%之间，对于液化天然气和液化石油气等类型的轻烃混合物尤其成功。,用多常数状态方程计算流体的热力学性质，必须借助计算机进行。,2.3.对比态原理及其应用,对比态原理认为，在相同的对比态下（Tr、Pr、Vr 三个对比参数中有两个相同），所有的物质表

22、现出相同的性质（真实气体的共性）。,2.3.1 对比态原理,2.3.1,这个原理虽然不十分严格，却是一个很有用的近似。,2.3.1,根据对比态原理，对不同的气体，若两个对比参数相同，则第三个对比参数也应相同。而Zc 基本是常数，则Z 应相同。上式只有在各种气体的临界压缩因子Zc 相等的条件下，才能严格成立。实际上，物质的Zc 在0.2 到0.3 范围内变动，不是常数。 因此简单的二参数（Pr ，Tr）对比态原理仅能应用于球形非极性的简单分子和组成、结构、分子大小近似的物质。,2.3.1,拓宽应用范围和提高计算准确性的有效方法是在简单对比态关系式中引入第三参数。 第三个参数可以是Zc ，也可采用

23、物质其他具普遍性的性质，如偏心因子，Riedel因子ac 等。在工程中使用较多的是以偏心因子为第三参数的对比态关联式。,偏心因子表达了气体分子间相互作用的位能与简单流体（Ar、Kr、Xe）的偏差。 偏心因子是根据某一规定对比温度下的对比蒸汽压确定的。 物质的对比蒸汽压的对数与对比温度有近似线性关系，即：,2.3.2 以偏心因子为第三参数的对比态原理,2.3.2,式中， 为对比饱和蒸汽压。,2.3.2,在临界点处，,此时，对比蒸汽压方程变成 0 = a - b，或 a = b。于是，对比蒸汽压方程可以表示为,因此，当以 对 作图时，得到一直线，a 是对比蒸汽压线的负斜率。,2.3.2,根据对比态

24、原理，如果这一原理准确，则所有物质应该具有相同的对比蒸汽压曲线，斜率a对所有物质都应该相同。但实际情况并非如此,每种物质都有一定的斜率。,图2-6,2.3.2,Pitzer注意到氩，氪，氙（球形对称的单原子分子）的数据全都位于同一根对比蒸气压曲线上，并且这条线通过 lgprs= -1和对比温度Tr=0.7 (1/Tr=1/0.7=1.43) 这一点。很明显，其他流体在Tr=0.7 处的纵坐标 值与氩、氪和氙在同一条件下的 lgprs 值的差能够表征该物质的某种特性，Pitzer就把这个差值定义为偏心因子，即,（2-37）,常见物质的Pc、Vc、Tc、Zc 和值见附录二。,2.3.2,由的定义,

25、 简单流体(氩,氪,氙)的值等于零，这些气体的压缩因子仅是Tr 和Pr 的函数。而对所有值相同的流体来说，若处于相同Tr、Pr下，其压缩因子必定相等。这就是Pitzer提出的普遍化三参数压缩因子关系式，表示为:,（2-38）,2.3.2,式中 Z0 和Z1 是Tr 和Pr 二者的复杂函数。根据实验数据得到的Z0 和Z1 与Tr 和Pr 的函数关系分别示于图2-7(a)、(b)和图2-8(a)、(b)。附录三给出了不同Pr（从0.2到9.0）和Tr（从0.35到4.0）下的Z0 和Z1 值，可供工程计算使用。可根据Pr，Tr 查Z0，Z1表。 Pitzer关系式对于非极性或弱极性的气体能够提供可

26、靠的结果，其误差在3以内；应用于极性气体时，误差达(510)；对于缔合气体，其误差要大得多；对量子气体，如氢，氦等，普遍化关系得不到好的结果。应当指出，普遍化关系并不能用来代替p -V-T 的可靠实验数据。,2.3.2,图2-7(a),Pr,Z0,2.3.2,图2-7(b),Z0,Pr,2.3.2,图2-8(a),Z1,Pr,2.3.2,例题2-3,图2-8(b),Z1,Pr,以下介绍Pitzer提出的普遍化第二Virial系数关系式。该式是一个解析计算式，计算时不需要查图，在工程应用中受到欢迎。,2.3.3 普遍化状态方程,2.3.3,将 ， ，代入舍项Virial方程（2-28a）中得到,

27、（2-42）,变量( BPc/RTc )是无因次的，可以看成对比第二维里系数。对于指定的气体来说，B仅仅是温度的函数，B的普遍化关系只与对比温度有关，而与对比压力无关。因此，Pitzer提出了如下的关联式,2.3.3,（2-43）,式中B0 和B1 只是对比温度的函数，用下述关系式表示：,（2-44b）,（2-44a）,上述关系式适用范围位于图2-9所示曲线上面的区域。该线是根据对比体积Vr2绘制的。当对比温度高于Tr = 4时，则对压力没有什么限制，但是须Vr2。对于较低的对比温度，允许的压力范围随着温度的降低而降低。但在对比温度约为0.9这一点之前，压力范围则受饱和（冷凝）蒸气压的限制。图

28、2-9的虚线表示饱和线。图2-9曲线以下的条件范围内，适用Pitzer提出的普遍化三参数压缩因子关系式（238）。,2.3.3,2.3. 3,图2-9,Tr,Pr,例题2-4,Vr 2,Vr 2,第三节完,上次课内容回顾,1、多常数状态方程Virial方程,2、对比态原理 对比态原理认为，在相同的对比态下（Tr、Pr、Vr 三个对比参数中有两个相同），所有的物质表现出相同的性质（真实气体的共性）。简单的二参数（Pr ，Tr）对比态原理仅能应用于球形非极性的简单分子和组成、结构、分子大小近似的物质。,3、以偏心因子为第三参数的对比态原理（1）偏心因子的概念 偏心因子表达了气体分子间相互作用的位能

29、与简单流体（Ar、Kr、Xe）的偏差。,（2）普遍化三参数压缩因子关系式对所有值相同的流体来说，若处于相同Tr、Pr下，其压缩因子必定相等。,（3）普遍化第二Virial系数关系式,2.4 真实气体混合物的P-V-T关系,化工生产中，处理的物系往往是多组分的真实气体混合物。目前虽然已收集，积累了许多纯物质的数据，但混合物的实验数据很少，不能满足工程设计的需要。解决这一问题的有效方法是在纯组分性质和混合物性质之间建立起联系，用纯物质性质来预测或推算混合物性质。描述纯物质性质和混合物性质之间联系的函数式称为混合规则。纯气体的P-V-T关系式借助于混合规则便可推广到气体混合物。,2.4,目前使用的混

30、合规则绝大部分是经验的，是从大量实际应用中总结归纳后建立起来的。 对于真实气体混合物，通常依赖某种混合规则，将混合物视为假想纯物质，计算求出混合物气体的虚拟临界参数，这样就可以把纯物质的对比态方法应用到混合物上。,2.4.1 混合规则与虚拟临界参数,2.4.1,Kay提出的最简单的混合规则，将混合物虚拟临界参数表示为,2.4.1,（2-48）,式中Tcm、Pcm为虚拟临界温度和压力；yi为组分i 的摩尔分数。求得虚拟临界参数后，混合物即可作为假想的纯物质进行计算。 Kay规则不含有组分间的相互作用项，因此该混合规则不能真正反映混合物性质。对于组分差别很大的混合物，尤其是含有极性组分或有缔合成二

31、聚体倾向的体系，这些混合规则不能适用。,混合物虚拟临界参数的应用,2.4.1,查图或表,假想的纯物质,气体混合物,混合规则,由统计力学可以导出气体混合物第二Virial系数为,2.4.2 气体混合物的第二Virial系数,2.4.2,（2-50）,这是少数从理论推导的混合规则之一。式中yi代表混合物中各组分的摩尔分数；Bij 表示组分i和 j 之间的相互作用，i 和 j 相同，表示同类分子作用，i 和 j 不同表示异类分子作用，且Bij=Bji。,对于二元混合物，式（2-50）展开为,2.4.2,（2-51）,式中B11、B22是纯物质1和2的第二Virial系数。Bij代表混合物性质，称为交

32、叉第二Virial系数，它们都只是温度的函数。纯物质第二Virial系数（例如 B11、B22）可按公式计算，交叉第二Virial系数（例如 B12）按以下的经验计算：,2.4.2,（2-52）,式中B0 和B1 是如同式（2-44a）(2-44b)一样的对比温度Tr的函数，Prausnitz对计算Tcij、Pcij和ij 提出如下的混合规则,（2-54）,（2-53）,2.4.2,（2-55）,（2-57）,（2-56）,在近似计算中Kij 可取为零。,2.4.2,由式（2-60）计算Bij，然后代入式（2-58），求得混合物的B。混合物压缩因子用式（2-28a）计算,气体混合物计算虽包含了

33、许多步骤，但每一步都不复杂，整个计算过程可方便地编成程序由计算机完成。,例题2-6,若把混合物看做一个整体，可以通过混合规则从纯物质状态方程得出混合物状态方程。一个纯物质状态方程可使用不同的混合规则，一个混合规则也可用于不同的状态方程，其优劣程度由实际计算结果检验。,2.4.3 混合物的状态方程,2.4.3,（1）RK型状态方程 RK型状态方程（van der Waals, RK, SRK, PR方程）用于混合物时，下面的混合规则被推荐用来计算方程中的常数 a 和 b：,bi 是纯组分常数，没有 b 的交叉项。aij 既包括纯组分常数（下标i =j），也包括交叉项（下标ij）交叉项aij 按下

34、式计算：,（2-58）,（2-59）,（2-60）,对二元混合物：,对三元混合物：,例题2-7,2.4.3,2.4.3,（2）Benedict-Webb-Rubin方程 该方程应用于混合物时，采用的混合规则为,（2-61）,对8个BWR方程常数，,r 的值分别为,2.4.3,（3）Martin-Hou方程,（2-62）,MH方程用于混合物时，有三种混合规则，临界参数混合规则、方程常数混合规则和温度函数混合规则。本文只介绍温度函数混合规则。,温度函数混合规则的通式为,式中n为正整数；yi为组分的摩尔分数。,2.4.3,若L代表方程常数b，则n=1,若L代表（2-32）中的温度函数fi (T)，则

35、,对于第二项的温度函数f2(T)有,2.4.3,对于二元体系展开为,不同的状态方程，有与之配套的不同的混合规则。通过修正混合规则，可改进和提高原方程应用于气体混合物时的计算精度和应用范围。,式中Qij为二元相互作用参数，由实验数据求得。,2.5 液体的P-V-T性质,液体的P-V-T 关系比较复杂，对液体理论研究远远不如对气体研究那样深入。但是与气体相比，液体的摩尔体积容易实验测定；除临界区外，压力和温度对液体容积性质影响不大。,2.5,2.5.1 经验关联式,(1) Tait方程,2.5.1,该方程的表达式,（2-63）,式中D、E 在给定温度下为常数； 和 为指定温度下，该液体在参比态时的

36、体积和压力。当有足够的数据，求出D、E之后，就可计算沿着等温线的PV关系。此方程可以用于很高的压力。,（2）Chueh-Prausnitz方程,2.5.1,该方程的表达式,（2-64）,式中,和 分别为体系温度下的饱和液体密度和饱和蒸汽压。,（3）修正Rackett方程,2.5.1,（2-65）,式中VS 是饱和液体的摩尔体积；ZRA值可查阅文献，若文献中没有，则用下式估算,（2-66）,该式计算仅仅需要临界参数，所得结果误差最大为7%左右，通常为(12)。但它不能准确预测临界体积VC，除非ZRA=ZC。,2.5.2 普遍化关联式,由Lydersen等人提出的计算液体密度的普遍化关联式为,2.

37、5.2,（2-68）,式中 为ZC0.27 时液体的对比密度，可从图2-10中查取； 是ZC0.27时的液体的对比密度；D为校正系数，可从图2-11中读得。,2.5.2,图2-10,2.5.2,图2-11,液体的对比密度定义为,2.5.2,（2-69）,因此，已知临界数据，便可由式（2-68）和式（2-69）计算液体密度或体积。,例题2-8,2.5.2,第二章完,1、纯物质P-V-T的关系（1）临界点,本章小结,(2)直线直径定律随着温度的变化，饱和液体和饱和蒸汽的密度迅速改变，但两者改变的总和变化甚微。当以饱和液体和饱和蒸汽密度的算术平均值对温度作图时，得到一条近似的直线，这结果称为直线直径

38、定律。常用于临界密度的实验测定。,2、气体的状态方程（1）状态方程定义,用数学形式表示的体系 P-V-T 之间的关系式，称为状态方程（EOS） 。,用状态方程可精确地代表相当广泛范围内的 p-V-T 数据，从而大大减少实验测定的工作量；即 关联和推算的方法。根据有限的实验数据关联出状态方程，有了状态方程就可推算更大温度、压力范围的数据；用状态方程可计算不能直接从实验测定的其他热力学性质； 用状态方程可进行相平衡计算，如计算饱和蒸气压、混合物汽液平衡、液液平衡等。尤其在计算高压汽液平衡时的简捷、准确、方便，是其他方法不能与之相比的。,(2)状态方程价值表现,(3)状态方程,（1）van der

39、Waals方程,(2) Ridlich-Kwang方程,(3)Soave-Ridlich-Kwang方程,(4)Peng-Robinson方程,(4) Patel-Teja方程,2.2.2,表2-3 立方型方程的另一形式,(2-22),(2-23),(2-24),(2-25),3、多常数状态方程Virial方程,4、对比态原理 对比态原理认为，在相同的对比态下（Tr、Pr、Vr 三个对比参数中有两个相同），所有的物质表现出相同的性质（真实气体的共性）。简单的二参数（Pr ，Tr）对比态原理仅能应用于球形非极性的简单分子和组成、结构、分子大小近似的物质。,5、以偏心因子为第三参数的对比态原理（1

40、）偏心因子的概念 偏心因子表达了气体分子间相互作用的位能与简单流体（Ar、Kr、Xe）的偏差。,（2）普遍化三参数压缩因子关系式对所有值相同的流体来说，若处于相同Tr、Pr下，其压缩因子必定相等。,（3）普遍化第二Virial系数关系式,例 题,例题 2-1 例题 2-2 例题 2-3 例题 2-4 例题 2-6 例题 2-7 例题 2-8,例 题 2-1,将1kmol氮气压缩贮于容积为0.04636m3，温度为273.15K的钢瓶内。问此时氮气的压力有多大？ 分别用理想气体方程，RK 方程和SRK方程计算。其实验值为101.33MPa 。(约1000 atm),从附录二中查得氮的临界参数为

41、Tc=126.2K , Pc= 3.394MPa , =0.040,解：,氮气的摩尔体积为,（1） 理想气体方程,误差,例 题 2-1,(2) RK方程,将TC , PC值代入式（2-7a）和式(2-7b)得,代入式（2-6）得,误差,例 题 2-1,(3) SRK方程,将代入式（2-9d）得,从式(2-9a), 式(2-9b)得,=273.15/126.2=2.1644，代入式（2-9c）得,例 题 2-1,例 题 2-1,将上述值代入式(2-8),误差,上述计算表明，在高压低温下(即使是非极性难液化的气体氮)理想气体方程根本不能适用，RK方程亦有较大误差；SRK方程的计算准确度较好。,EN

42、D,例 题 2-2,试用SRK和PR方程分别计算异丁烷在300K，3.704105Pa时饱和蒸气的摩尔体积。已知实验值为V=6.03110-3m3/mol。,由附录二查得异丁烷临界参数为,解：, SRK方程,例 题 2-2,按式（2-22）,和式（2-25）,，,迭代计算。,例 题 2-2, PR方程,例 题 2-2,和式（2-25）,迭代计算。,例 题 2-2,按式（2-23）, 取初值Z=1,迭代计算结果见下表。,SRK方程,误差,例 题 2-2,PR方程,从以上计算看出，迭代收敛速度很快，仅迭代四次即得终点。SRK和PR方程都有较好的准确度，其结果均可为工程界接受。,误差,END,例 题

43、 2-2,例 题 2-3,用pitzer的普遍化关系式计算甲烷在323.16k时产生的压力。已知甲烷的摩尔体积为1.2510m3/mol，压力的实验数值为1.875107Pa。,从附录二查得甲烷的临界参数为,解：,因此,Pr 不能直接计算，需迭代求解。,而,因此，,据式(2-46),（A）,（B）,例 题 2-3,设Z 值代入（A）式求出 pr，根据Tr , pr 值查图2-7 和图 2-8 得Z(0) 和Z(1) ;再将Z(0)、Z(1) 代入（B）式求得Z 值。比较Z 的计算值与假设值，如相差较大则代入（A）式重新计算，直至迭代收敛。,例 题 2-3,假设Z 值,（A）,Pr,Tr,查图得

44、Z(0) 和Z(1),（B）,Z （计算值）,比较,迭代的结果为 Pr = 4.06时，Z0.877,误差,END,例 题 2-3,例 题 2-4,质量为500g的氨气储于容积为0.03m3的钢弹内，钢弹浸于温度为65的恒温浴中试用普遍化第二维里系数计算氨气压力，并与文献值p=2.382MPa比较。,从附录二查得氨的临界参数,解：,求氨的摩尔体积,由于Vr2故适用普遍化第二维里系数关系式。,由式(2-43),所以,例 题 2-4,由式(2-42),则,计算值与实验值2.382MPa是非常接近的，说明在压力不太高时，普遍化第二virial系数用于极性物质（NH3）广度性质计算，结果也较满意。,E

45、ND,例 题 2-4,即,例 题 2-6,试求CO2(1)-C3H8(2)体系在311K和1.50MPa的条件下的混合物摩尔体积，两组分的摩尔比为3:7。,计算所需的有关数据列表如下：,解：,，,由式（2-51）得,END,例 题 2-6,例 题 2-7,试用RK方程计算二氧化碳和丙烷的等分子混合物在151和13.78MPa下的摩尔体积。,临界性质由例 2-6 给出，把这些值代入式（2-7a、b）和式（2-68）得如下结果（取kij=0）：,解：,混合物常数由式（2-58）和式（2-59）求出：,例 题 2-7,表2-2,将上面两个量代入方程式（2-22）和式（2-25）得,（B）,（A）,联立求解方程（A）、（B）得,所以混合物摩尔体积为,END,例 题 2-7,例 题 2-8,估算37的饱和液氨的摩尔体积； 估算37和10.13MPa下液氨的摩尔体积。,由附录二查得氨的临界参数：,解：,，, 采用修正的Rackett方程,查文献得氨的ZRA=0.2465，按式（2-65）,将此值与实验值2.914105比较，误差为0.2。, 采用普遍化密度关系式,例 题 2-8,对比参数,根据Tr，pr 值查图2-10得r = 2.38，查图2-11（Zc 0.27）得D = 7.2，代入式（2-68）得,和实验值2.68105 m3/mol相比，误差为1.7。,END,例 题 2-8,