第二章放映稿

第二章放映稿第一页，共一百四十二页，2022年，8月28日2推算流体p-V-T行为的途径1）状态方程(EOS)

p-V-T关系的解析式。2）对应态原理(CSP)一种特别的状态方程，以对比参数来表达方程，使流体性质在对比状态下便于比较，并统一到较好的程度。第二页，共一百四十二页，2022年，8月28日3p-V-T关系和状态方程的重要性在计算热力学性质时需要输入流体最基本的性质以及表达系统特征的模型。

状态方程不仅本身是重要的p-V-T关系式，而且从p-V-T的角度反映了系统的特征，是经典热力学中推算其它性质不可缺少的模型之一。

第三页，共一百四十二页，2022年，8月28日4本章主要内容1）纯物质的p-V-T行为2）常见的状态方程3）常用的对应态原理4）混合法则第四页，共一百四十二页，2022年，8月28日§2-2

p-V-T相图SLGCVV/SV/LS/LAB该图是表示纯物质在平衡状态下，压力、摩尔体积与温度关系的p-V-T曲面。第五页，共一百四十二页，2022年，8月28日相图包括：1单相区：S、L和V（G）分别表示固相、液相和蒸汽（气相）；2两相共存区：S/L、V/S和V/L分别代表固/液、汽/固、汽/液两相平衡区第六页，共一百四十二页，2022年，8月28日3临界点C：汽/液共存的最高温度或压力点，该点的温度、压力和摩尔体积分别称为临界温度Tc、临界压力Pc和临界体积Vc。数学上的关系表示为：（在C点）（在C点）第七页，共一百四十二页，2022年，8月28日流体在临界的特性和临界参数在状态方程研究中有重要作用。在T>Tc和p>pc的区域内，气体和液体变得不可区分，称为超临界流体。临界点附近，流体的许多性质有突变的趋势，如密度、溶解其它物质的能力等，已开发的工业过程有超临界分离技术、超临界化学反应等。第八页，共一百四十二页，2022年，8月28日4饱和线：ACB是汽/液两相共存区的边界线。AC为饱和液体线也称为泡点线，BC为饱和蒸汽线也称为露点线。第九页，共一百四十二页，2022年，8月28日5三相线：通过A、B的直线，是三个两相平衡区的交界线。在三相线上有固定的温度、压力，此状态下的纯物质处于气-液-固三相平衡。第十页，共一百四十二页，2022年，8月28日AB三相点

纯物质的p–T图

纯物质的p–V图

第十一页，共一百四十二页，2022年，8月28日

特定条件时，存在过热液体—在一定温度下，当压力低于饱和蒸汽压（或一定压力下，温度高于其沸点），仍能以液体形式存在过冷蒸汽—压力高于同温度下的饱和蒸汽压（或温度低于同压力的沸点），仍能以蒸汽形式存在。过冷蒸汽和过热液体都是亚稳定状态。第十二页，共一百四十二页，2022年，8月28日§2-3状态方程（EOS）状态方程是流体p-V-T的解析表达式。从研究方法上看，状态方程可以分为理论型、经验型和半理论型；从形式上看，又可以分为立方型（可化为V的三次多项式）和高次型。第十三页，共一百四十二页，2022年，8月28日一般采用如下分类：1立方型状态方程，如vanderWaals、RK、SRK、PR等2多常数状态方程，如virial、BWR、MH等3理论型状态方程第十四页，共一百四十二页，2022年，8月28日

第一、第二类直接以工业应用为目标，在分析、探找流体性质规律的基础上，结合一定的理论，由半经验方法建立模型。有若干模型参数需从实验数据确定。本章主要介绍一、二类方程

第三类从微观出发，是分子间相互作用与统计力学结合的结果，离实际使用有差距。第十五页，共一百四十二页，2022年，8月28日

状态方程既有将p作为函数（T，V作自变量）的形式，如p=p（T，V），也有以V为函数（T，p作自变量）的形式，如V=V（T，p）这两种形式所适用的范围有所不同。目前以前者为普遍，也是介绍和应用的重点。第十六页，共一百四十二页，2022年，8月28日如果将以T，V为自变量的状态方程，用于以T，p为独立变量的系统的性质计算，要先计算V（类似于数学上的求反函数）。对于T，p为自变量的情况也是相似的。第十七页，共一百四十二页，2022年，8月28日小结1化工热力学的含义运用经典热力学的原理，结合反映系统特征的模型，解决工业过程（特别是化工过程）中热力学性质的计算和预测、相平衡和化学平衡计算、能量的有效利用等实际问题。第十八页，共一百四十二页，2022年，8月28日2化工热力学的主要内容：原理-模型-应用3基本概念：1)systemandsurroundingclosed,openandisolated2)intensiveandextensivevariables3)statevariable4)equilibriumstate5)reversibleprocess6)themodynamicsprocessandcycle第十九页，共一百四十二页，2022年，8月28日4热力学性质的计算方法变量分析、普遍关系、引入模型、数学求解5纯物质的p-V-T相图（phasediagram）第二十页，共一百四十二页，2022年，8月28日6纯物质的p-V图7状态方程

流体p-V-T的解析表达式第二十一页，共一百四十二页，2022年，8月28日§2-4立方型状态方程

立方型方程可以化为V的三次方的形式。一般由斥力(repulsion)项和引力(attraction)项组成。一般情况下，prep>0,而patt<0第二十二页，共一百四十二页，2022年，8月28日典型的立方型方程，它们的常数可以通过普遍化关系式，从临界参数Tc,pc和偏心因子ω计算。特别是SRK和PR方程在工程上有广泛的应用。介绍几种重要的方程第二十三页，共一百四十二页，2022年，8月28日1vanderWaals（vdW）方程1）能同时表达汽液两相2）可计算出临界点3）准确度有限公式2-6第二十四页，共一百四十二页，2022年，8月28日由可解得（2-7）第二十五页，共一百四十二页，2022年，8月28日将a、b代入vdW方程，并用于临界点，得或（2-8）第二十六页，共一百四十二页，2022年，8月28日以Tc和pc表达的vdW常数为（2-9）（2-10）第二十七页，共一百四十二页，2022年，8月28日vdW方程形式简单，固定临界压缩因子0.375，计算容易。实际流体压缩因子数值在0.23~0.29之间。不足之处。vdW方程具有深远的理论意义，立方型状态方程，多数是基于vdW方程的改进。第二十八页，共一百四十二页，2022年，8月28日2Redlich-Kwong（RK）方程

斥力相与vdW相同，引力项与T是一个简单的T-0.5关系公式2-11第二十九页，共一百四十二页，2022年，8月28日公式2-12公式2-13第三十页，共一百四十二页，2022年，8月28日RK方程的Zc=1/3=0.333，仍偏大。RK方程较成功地用于气相p-V-T的计算，但液相的效果较差，不能预测纯流体的蒸汽压。第三十一页，共一百四十二页，2022年，8月28日3Soave（SRK）方程

1972年，Soave修正了RK方程

(2-14)(2-15)(2-18)第三十二页，共一百四十二页，2022年，8月28日(2-16)(2-17)临界等温线上，RK方程与SRK完全一样。因此SRK方程的Zc=1/3=0.333第三十三页，共一百四十二页，2022年，8月28日优点：1）较RK方程提高了表达纯物质汽液平衡的能力2）可用于混合物的汽液平衡计算，在工业上获得广泛应用。缺点：RK、SRK方程预测液相摩尔体积不够准确，Zc偏大第三十四页，共一百四十二页，2022年，8月28日4Peng-Robinson（PR）方程

（2-19）采用了类似于SRK方程中的a表达式（2-22）第三十五页，共一百四十二页，2022年，8月28日（2-20）（2-21）计算得临界压缩因子Zc=0.307

优点：PR方程预测液体摩尔体积的准确度较SRK有了明显改善。第三十六页，共一百四十二页，2022年，8月28日立方型状态方程的特点：1）形式简单，2）方程常数可以进行普遍化处理3）可以解得方程的体积根。4）由于内在缺陷，难以在大范围应用第三十七页，共一百四十二页，2022年，8月28日§2-5多常数状态方程立方型方程的发展是基于了vdW方程，而多常数状态方程是与virial方程相联系。

第三十八页，共一百四十二页，2022年，8月28日

多常数的高次型状态方程涉及更多的流体物性信息，适用范围更大，准确性更高，方程的预测效果更好。但计算量和复杂性增大。借助电算使其研究受到重视。第三十九页，共一百四十二页，2022年，8月28日1virial方程virial方程分为密度型(2-23)(2-24)压力型第四十页，共一百四十二页，2022年，8月28日B、C···（或B’、C’···）称作virial系数。任何状态方程可以通过级数展开，转化为Virial方程形式。

如对vdW方程展开成级数方程第四十一页，共一百四十二页，2022年，8月28日比较后即可将vdW方程和virial系数联系起来。在取无穷项的情况下，两者是等价的。第四十二页，共一百四十二页，2022年，8月28日1）virial系数的意义：微观上，反映了分子间的相互作用。第二virial系数B反映了两分子间的相互作用，第三virial系数C反映了三分子间的相互作用，…..宏观上，virial系数仅是温度的函数第四十三页，共一百四十二页，2022年，8月28日实际应用中常采用两项virial截断式高密度时高次相的影响非常敏感。第四十四页，共一百四十二页，2022年，8月28日2）第二virial系数的关联式①对应态关联式由Tsonopoulos提出，较多的应用于非极性、弱极性物质第四十五页，共一百四十二页，2022年，8月28日（2-26）（2-27）第四十六页，共一百四十二页，2022年，8月28日②从P-V-T数据确定第四十七页，共一百四十二页，2022年，8月28日由等温条件下的p-V-T数据，在密度不太高的条件下近似一条直线截距为第二virial系数B，斜率为第三virial系数C。对作图用外推至第四十八页，共一百四十二页，2022年，8月28日{B第四十九页，共一百四十二页，2022年，8月28日③利用Z~p图第二virial系数是与Z～p图上的等温线在p→0时的斜率有关。将V=ZRT/p代入式（2-23），得第五十页，共一百四十二页，2022年，8月28日p→0时，第三及以后各项为更高阶无穷小，所以（2-28）经微分处理得（2-29）第五十一页，共一百四十二页，2022年，8月28日随着温度的升高，Z～p图上的等温线在p→0时的斜率由负变为正，第二virial系数B只在某一特定温度下变为零，这一温度称为Boyle温度，用TB表示，即

第五十二页，共一百四十二页，2022年，8月28日目前高阶的virial系数的估算尚不成功。第五十三页，共一百四十二页，2022年，8月28日高次型状态方程与virial方程有一定的关系。较多见到的多常数高次型方程有BWR方程和马丁-侯方程（简称为MH方程）已经广泛地应用于化工及其它领域中，与立方型方程相比，高次型方程的准确性高，适用范围广，但计算量稍大。第五十四页，共一百四十二页，2022年，8月28日BWR方程的常数是从实验数据拟合得到的，也有相关的普遍化关联式。MH方程的常数能从纯物质的临界参数和蒸汽压数据计算出来。从形式上看，MH方程的数据规律性很好。第五十五页，共一百四十二页，2022年，8月28日2Benedict-Webb-Rubin（BWR）方程第五十六页，共一百四十二页，2022年，8月28日第一个能在高密度区表示流体p-V-T和计算汽液平衡的多常数方程，在工业上得到了一定的应用。第五十七页，共一百四十二页，2022年，8月28日BWR方程在应用中不断被改进，常数不断增加，准确性和使用范围也不断提高，但方程形式愈加复杂由于BWR方程的数学形式上的规律性不好，给数学推导、数值求根及方程的改进和发展带来一定的不便。第五十八页，共一百四十二页，2022年，8月28日3Martin-Hou（MH）方程

我国学者侯虞钧和美国的Martin教授在20世纪50年代初提出，数学形式整齐2-31第五十九页，共一百四十二页，2022年，8月28日温度函数很有规律第六十页，共一百四十二页，2022年，8月28日其中第六十一页，共一百四十二页，2022年，8月28日9个常数反映了较多的热力学性质的普遍化规律，只需输入纯物质的临界参数和一点的蒸汽压数据，就能从数学公式计算所有的常数。简便、可靠、适用范围广，可用于非极性至强极性化合物。是比较优秀的状态方程。第六十二页，共一百四十二页，2022年，8月28日MH方程已广泛用于流体p-V-T、汽液平衡、液液平衡、焓等热力学性质推算，并被用于大型合成氨装置的设计和过程模拟中。第六十三页，共一百四十二页，2022年，8月28日例1：P18用KR方程计算异丁烷：1）在420K和2MPa时的摩尔体积。2）380K时的饱和气、液相摩尔体积，已知蒸汽压2.25MPa。解：1）查附录Tc=408.1K,pc=3.648MPa,ω=0.1762）计算方程常数，写出方程形式第六十四页，共一百四十二页，2022年，8月28日第六十五页，共一百四十二页，2022年，8月28日3）图解求根由方程在图p-V上做出380K和420K两条等温线420K380K2.252第六十六页，共一百四十二页，2022年，8月28日4）迭代求解对于气体对于液体第六十七页，共一百四十二页，2022年，8月28日例2：欲在一7810cm3的钢瓶中装入1000g的丙烷，且在253.2℃（526.35K）下工作，若钢瓶的安全工作压力10MPa，问是否有危险？解：1）查临界参数及ωTc=369.85K,Pc=4.249MPa,ω=0.152

第六十八页，共一百四十二页，2022年，8月28日2）应用PR方程，由软件可计算得可以容纳的丙烷。即所以会有危险。第六十九页，共一百四十二页，2022年，8月28日§2-6对应态原理（CSP)对应态原理也是一种状态方程，以对比参数来表达状态方程。对比参数是指流体的真实值与临界值的比值，包括对比温度Tr、对比压力pr、对比体积Vr。第七十页，共一百四十二页，2022年，8月28日

对应态原理是预测流体性质最有效的方法之一。主要思路是从已知的参考流体的性质（或状态方程）获得研究流体的性质（或状态方程），其发展主要沿两条途径：一是多参数对应态原理，应用多的是三参数对应态原理一是形状因子对应态原理第七十一页，共一百四十二页，2022年，8月28日1二参数对应态原理vanderWaals首先提出，经过运算得vdW方程的对比形式为：（2-33）式中只含有纯数值和对比参数，即第七十二页，共一百四十二页，2022年，8月28日或表明在相同的对比温度和对比压力下，任何流体的对比体积（或压缩因子）是相同的。两参数对应态原理不够精确，只适用于简单的球形流体。第七十三页，共一百四十二页，2022年，8月28日2三参数对应态原理1）Zc作为第三参数Lydersen等引入Zc作为第三参数，压缩因子表示为第七十四页，共一百四十二页，2022年，8月28日2）偏心因子作为第三参数Pitzer研究了蒸汽压数据，发现简单流体在Tr=0.7时的对比蒸汽压近似等于0.1即而其它流体（除H2和He）的由此差别提出了偏心因子ω的概念第七十五页，共一百四十二页，2022年，8月28日简单流体ω=0，其它流体ω>0。偏心因子表达了一般流体与简单流体分子间相互作用的差异。第七十六页，共一百四十二页，2022年，8月28日第七十七页，共一百四十二页，2022年，8月28日1立方型状态方程立方型方程可以化为V的三次方的形式，可以得到解析的体积根。一般由斥力相和引力相组成1.1vanderWaals（vdW）方程

利用临界点的性质得到常数a,b的值小结第七十八页，共一百四十二页，2022年，8月28日1.2R-K方程改进了引力项，与T是一个简单的T-0.5关系第七十九页，共一百四十二页，2022年，8月28日1.3SRK方程常数a考虑了烃类在不同温度下的蒸汽压数据第八十页，共一百四十二页，2022年，8月28日1.4PR方程

拟合了蒸气压数据，预测液体摩尔体积的准确度较SRK有了明显改善第八十一页，共一百四十二页，2022年，8月28日2

virial方程两项截断式第八十二页，共一百四十二页，2022年，8月28日2.1virial系数的意义：微观上，反映了分子间的相互作用。第二virial系数B反映了两分子间的相互作用，第三virial系数C反映了三分子间的相互作用，…..宏观上，virial系数仅是温度的函数第八十三页，共一百四十二页，2022年，8月28日2.2第二virial系数的关联式①

对应态关联式由Tsonopoulos提出，较多的应用于非极性、弱极性物质第八十四页，共一百四十二页，2022年，8月28日②

从P-V-T数据确定由{B第八十五页，共一百四十二页，2022年，8月28日③利用Z~p图第八十六页，共一百四十二页，2022年，8月28日Boyle温度TB第八十七页，共一百四十二页，2022年，8月28日3Benedict-Webb-Rubin（BWR）方程第一个能在高密度区表示流体p-V-T和计算汽液平衡的多常数方程4Martin-Hou（MH）方程

数学形式整齐，温度函数有规律，是比较优秀的状态方程第八十八页，共一百四十二页，2022年，8月28日5对应态原理（CSP)对应态原理是以对比参数来表达的状态方程。对比参数是指流体的真实值与临界值的比值，包括对比温度Tr、对比压力pr、对比体积Vr。第八十九页，共一百四十二页，2022年，8月28日6偏心因子：表达了一般流体与简单流体分子间相互作用的差异。第九十页，共一百四十二页，2022年，8月28日三参数方程为：是简单流体的压缩因子表示代表研究流体相对于简单流体的偏差第九十一页，共一百四十二页，2022年，8月28日3）Lee—Kesler方程1975年，由Lee和Kesler提出的三参数对应态原理的解析形式。除简单流体外，选择正辛烷作参考流体(r)，其偏心因子ω(r)=0.3978。以得第九十二页，共一百四十二页，2022年，8月28日A）Z(0)、Z(r)分别代表简单流体和参考流体的压缩因子B）研究流体与参考流体的性质越接近，预测结果的准确性和可靠性越高。（2-37）第九十三页，共一百四十二页，2022年，8月28日C）在L-K方程中，简单流体和参考流体的状态方程均采用修正的BWR方程。简单流体的方程常数由简单流体的压缩因子和焓的数据拟合，参考流体的方程常数由正辛烷的数据得到。第九十四页，共一百四十二页，2022年，8月28日4）Teja方程1980年，Teja发展的三参数对应态原理采用了两个非球形的参考流体两个参考流体r1、r2可以采用不同的状态方程来描述。允许根据研究流体的性质对参考流体进行适当选择。第九十五页，共一百四十二页，2022年，8月28日3形状因子对应态原理基于保形溶液理论

f，h称为保形参数，与研究流体和参考流体的Tc、Vc之比有关第九十六页，共一百四十二页，2022年，8月28日

θ，Φ

称为形状因子，研究流体和参考流体的性质非常相似时可认为是近似的保形流体对，θ≈1，Φ≈1。一般情况θ，Φ是偏离1的。获得方程取决于两个关键因素第九十七页，共一百四十二页，2022年，8月28日A)形状因子θ，Φ

决定于研究流体和参考流体的性质。Leach等人以甲烷参考流体，针对烃类，能用于碳氢化合物的p-V-T和汽液平衡等性质的计算。第九十八页，共一百四十二页，2022年，8月28日B)参考流体的状态方程Z0

通常采用多常数的高次型状态方程。例：P222-4估计正丁烷在425.2K和4.4586MPa时的压缩因子第九十九页，共一百四十二页，2022年，8月28日解：1）查参数Tc=425.18K,pc=3.797MPa,ω=0.1933）查表

内差

2）计算第一百页，共一百四十二页，2022年，8月28日§2-7流体的饱和热力学性质常用的流体的饱和热力学性质主要有蒸汽压、汽化焓、汽化熵、饱和汽相摩尔体积、饱和液相摩尔体积。1饱和蒸汽压、汽化焓和汽化熵第一百零一页，共一百四十二页，2022年，8月28日1）饱和蒸汽压纯物质在一定温度下，能使汽液共存的压力为蒸汽压。p-T图上表达汽液平衡的蒸汽压曲线始于三相点而止于临界点。蒸汽压表达物性的唯一性，是温度的一元函数，其解析式为蒸汽压方程。第一百零二页，共一百四十二页，2022年，8月28日Clapeyron方程反映了蒸汽压关系（2-42）第一百零三页，共一百四十二页，2022年，8月28日仅是温度的函数假定为不随温度变化的常数B修正后得Antoine方程应用注意常数的使用条件则第一百零四页，共一百四十二页，2022年，8月28日

在缺乏蒸汽压数据或蒸汽压方程常数的条件下，可用经验方法估计第一百零五页，共一百四十二页，2022年，8月28日2）汽化焓汽液相平衡转化过程的潜热，仅是温度的函数，是重要的物性数据，焓值随温度升高而下降，达临界点时汽化焓为零。可以由Clapeyron方程计算，或由状态方程推算。常用Watson经验式第一百零六页，共一百四十二页，2022年，8月28日3）汽化熵汽液相平衡转化过程的熵变化。等于汽化焓除以汽化温度。2饱和液体摩尔体积

SRK、PR、BWR、MH-81等状态方程可用于气、液相性质的计算。但一般情况下，液相误差大于气相。第一百零七页，共一百四十二页，2022年，8月28日若只计算饱和液体体积，可用饱和液体摩尔体积方程。1）Rackett方程对大多数物质的计算误差为2%第一百零八页，共一百四十二页，2022年，8月28日2）修正的Rackett方程SpancerandDanner的修正式为：引入的Rackett常数ZRA需实验数据拟合，与ZC差别不大。对于存在缔合的物质，结果仍不满意。第一百零九页，共一百四十二页，2022年，8月28日

Campbell等将ZRA改为下列温度的函数：准确度有很大改善第一百一十页，共一百四十二页，2022年，8月28日3）Tait方程表达了等温线上液体的V-p关系。等温条件下，液体的摩尔体积随压力的增加而减小，只有在高压下才会明显。

p0、V0是给定温度下，某一已知的参考状态的压力和摩尔体积；D、E是两个与温度有关的常数。第一百一十一页，共一百四十二页，2022年，8月28日例：P252-61）饱和蒸汽压由Antoine方程计算第一百一十二页，共一百四十二页，2022年，8月28日2）饱和液相体积修正的Rackett方程第一百一十三页，共一百四十二页，2022年，8月28日3）饱和气相的摩尔体积采用virial截断式由普遍化关联式得到第二virial系数B解二次方程得体积根第一百一十四页，共一百四十二页，2022年，8月28日§2-8混合法则研究混合物性质时，常将混合物看成一个虚拟的纯物质，并具有虚拟的特征参数，将这些虚拟的特征参数代入纯物质的状态方程中就可以计算混合物的性质。混合物的虚拟参数强烈的依赖于混合物的组成。第一百一十五页，共一百四十二页，2022年，8月28日混合法则是指混合物的虚拟参数与混合物的组成和所含的纯物质的参数之间的关系式。通常在一定的理论指导下，引入适当的经验修正，再结合实验数据才能确定下来。混合物系统的符号和纯物质符号的规定见P27表2-1第一百一十六页，共一百四十二页，2022年，8月28日1Virial方程的混合法则第二Virial系数的混合法则为Bij由同温度下纯组分Virial系数Bi、Bj得到第一百一十七页，共一百四十二页，2022年，8月28日2立方型方程的混合法则两参数立方型方程中，b与分子的大小有关。a是分子间相互作用力的度量。第一百一十八页，共一百四十二页，2022年，8月28日RK方程中是相互作用参数，由实验数据拟合得到。近似认为SRK、PR方程中第一百一十九页，共一百四十二页，2022年，8月28日3BWR方程r数值见表2-24MH-81方程采用温度函数混合法则第一百二十页，共一百四十二页，2022年，8月28日k=3，4，5第一百二十一页，共一百四十二页，2022年，8月28日混合物状态方程的温度函数与纯物质相应的温度函数保持相同的符号。大多数情况下是二元相互作用参数，一般条件下，第一百二十二页，共一百四十二页，2022年，8月28日5修正的Rackett方程由纯物质的参数计算液体混合物的摩尔体积第一百二十三页，共一百四十二页，2022年，8月28日6对应态原理三参数对应态原理，常用临界参数混合法则。ω一般采用线性混合法则，将形状因子对应态原理推广到混合物，需要保形参数的混合法则。第一百二十四页，共一百四十二页，2022年，8月28日§2-9状态方程体积根的求解1状态方程体积根在p-V图上的几何形态一般以p为显函数的立方型状态方程可化为关于V的三次方程，如SRK方程第一百二十五页，共一百四十二页，2022年，8月2