化工热力学第五章教案

1授课内容第五章非均相系统的热力学性质计算§5-1引言1非均相系统热力学性质计算的内容1)确定平衡状态要进行非均系统热力学性质计算，必须首先确定相平衡，所以混合物的相平衡计算十分重要。确定了平衡状态后,各相的性质计算就属于均相性质的范畴。2）相平衡的类型不同的非均相系统包含了不同的相平衡，典型的有汽液平衡(VLE)、液液平衡(LLE)、固液平衡(SLE)等，虽然类型不同，但计算原理和方法类似。重点讲解汽液平衡的计算。3)计算互成平衡的各个相的性质对于混合物，相平衡关系主要是指T，p和各相的组成之间的关系，作为非均相系统的性质，还包括互成平衡的各相的其它热力学性质。它们的计算需要将混合物的相平衡准则与反映混合物特性的模型（状态方程+混合法则或活度系数模型）结合起来。2Gibbs-Duhem方程在热力学计算中的作用Gibbs-Duhem方程是混合物中各组分的偏摩尔性质的约束关系，不仅在检验偏摩尔性质的模型时非常有用，而且，由于有些偏摩尔性质（如lnγi等）与混合物的相平衡数据相联系，所以，Gibbs-Duhem方程在相平衡数据的检验和推算中也有重要作用。3本章的主要内容1）混合物的相图和相平衡计算2）汽液平衡数据的一致性检验3）热力学性质的推算和预测2§5-2混合物的汽-液平衡(本章的重点)N个组分的两相系统，在一定T,p下达到汽液平衡混合物的汽液平衡系统该系统的基本强度性质是T,p，汽相组成和液相组成共有2+(N-1)+(N-1)=2N个。由相律知，N元的两相平衡系统的自由度是f=N–2+2=N，若给定N个独立变量，其余N个强度性质就能确定下来，这是汽液平衡计算的主要任务。完成了汽液平衡计算，该非均相系统中的任何一个相的其它热力学性质就容易得到了1混合物的汽-液相图相律提供了确定系统所需要的强度性质数目，对于二元汽、液相混合物，其基本的强度性质是(T，p，x1，y1)，系统的自由度为f=2–M+2=4-M（M是相的数目），系统的最小相数为M=1，故最大的自由度数是f=3，表明最多需要3个强度性质来确定系统。这样，二元汽-液相图就需要表达成三维立体曲面形式。但在等温条件或等压条件，系统状态可以表示在二维平面上，在汽液共存时，M=2，f=1，故汽液平衡（VLE）关系就能表示成曲线。1）等压二元相图在固定压力条件下，单相区的状态可以表示在温度~组成的平面上，汽液平衡关系可以表示成温度~组成（T~x1和T~y1）的等压二元相图，还可表示为x~y曲线。见P107图5-2pTV相L相x1,x2,···,xN-1y1,y2,···,yN-13图5-2等压二元系统的相图2）等温二元相图在固定温度条件下，单相区的状态可以表示在压力~组成的平面上，汽液平衡关系可以表示成压力~组成（p~x1和p~y1）的等温二元相图。在实际应用中，还可以表示成二元汽液平衡关系曲线x~y图。见P108图5-3图5-3等温二元系统的相图3）理想系统泡点线在等温二元系统相图中，连接1sp和2sp的斜虚线代表了理想系统（即汽相是理想气体混合物，液相是理想溶液）的泡点线，理想系统的泡点线方程为4）真实液相互溶系统相对理想系统偏差的分类1211222121()sssssppypypxpxpppx44.1）一般正偏差系统泡点线位于理想系统的泡点线上方，但不产生极大值，称之为一般正偏差系统；4.2）一般负偏差系统泡点线位于理想系统的泡点线下方而又不产生极小值时，称为一般负偏差系统；4.3）有共沸点的系统若泡点线产生了极值点，称为共沸点。其沸点的温度和压力分别称为共沸温度（Taz）和共沸压力(paz)。在共沸点，泡点线与露点线相切，汽、液相组成相等，并称为共沸组成，即共沸点分为最高压力共沸点和最低压力共沸点。对于p-x-y图上的最高压力共沸点，一般会表现为T-x-y图上的最低温度共沸点。01x1y1p-x12sp01p-x12sp0110p-x1p-y1LV1siiippx1siiippxazaziixy0T-x1T-y1LV1min11sazazaziippxppxy1sp1sp5p-x-y图上的最低压力共沸点，一般也会表现为T-x-y图上的最高温度共沸点。5）二元部分互溶系统的等压相图有些混合物，汽液平衡系统中的液相出现了部分互溶（即分层液相）的情况，此时，系统实际上是汽-液-液三相平衡（VLLE）。由于汽液液平衡时M=3，在等温或等压条件下f=0，相图上的汽-液-液平衡关系是一个固定的三相点。图5-4二元部分互溶系统的等压相图2汽液平衡的准则和计算方法1）汽液平衡的准则N元系统的汽液平衡准则可以表示如下：2）汽液平衡的计算方法2.1）EOS法：若汽、液相的组分逸度系数可以用一个同时适合于汽、液两相的状态方程及其混合法则来计算，这种方法称为状态方程法，或简称EOS法。01x1y10T-x1T-y1LV110p-x1(1,2,,)vliiffiNmin11sazazaziippxppxy6EOS法对于状态方程的要求是很高的。2.2）EOS+γ法若液相中组分的逸度用活度系数计算，即（采用了对称归一化的活度系数），则汽液平衡准则为若采用不对称归一化的活度系数，即则汽液平衡关系为这种用状态方程和活度系数两个模型来处理汽液平衡的方法称为状态方程+活度系数法，或简称EOS+γ法。模型的选择主要由系统特征决定2.3）汽液平衡常数汽液平衡还可用汽液平衡常数Ki来表示，定义为,一般是T，p的函数3汽液平衡计算类型N元汽液平衡系统的自由度是N，故必须指定N个强度性质作为独立变量，才能确定汽液平衡。汽液平衡计算的目的是从指定的N个独立变量，确定其余基本的从属变量。指定N个独立变量的方案不同，构成了不同的汽液平衡计算类型。常见的汽液平衡计算类型见P111表5-2(1,2,,)vliiffiNvvlliiiiiifpyfpx(12)vliiiiyxiN则有，，，lliiiiffx(1,2,,)vliiffiN(1,2,,)vliiiiipyfxiN*,liSolventiiifHx*,(1,2,,)viiiSolventiipyHxiNiiiyKx()(1,2,,)(12)liiiviiliiiiviiyKEOSiNxyfKEOSiNxp或（），，，7表5-2常见的汽液平衡计算类型计算类型独立变量待确定的基本的从属变量1等温泡点计算T,x1,x2···xN-1p,y1,y2···yN-12等压泡点计算p,x1,x2···xN-1T,y1,y2···yN-13等温露点计算T,y1,y2···yN-1p,x1,x2···xN-14等压露点计算p,y1,y2···yN-1T,x1,x2···xN-15闪蒸计算T,p,z1,z2,···,zN-1x1,x2···xN;y1,y2···yN和η3.1）泡点计算第一、二类型是泡点计算，即确定某一组成的液体混合物在一定压力下的沸点（泡点温度）或一定温度下的蒸汽压（泡点压力），以及平衡汽相组成。要确定的N+1个基本的从属变量可以从N个平衡准则方程和汽相组成的归一化方程联列求解出来：3.2）露点计算第三、四类型是露点计算，即确定某一组成的汽体混合物在一温度下的露点压力或一定压力下的露点温度，以及平衡液相组成。待确定的N+1个基本从属变量也是从平衡准则和液相组成的归一化方程联列求解得到。3.3）闪蒸计算第五类型是闪蒸计算。闪蒸的名词来源于液体流过阀门等装置，由于压力突然降低而引起急骤蒸发，产生部分汽化，形成互成平衡的汽、液两相（也可以是汽相产生部分冷凝）。在一定温度、压力条件下，总组成为zi的混合物分为相互成平衡的汽、液两相，1(1,2,,)1vliiNiiffiNy1(1,2,,)1vliiNiiffiNxVVL8闪蒸计算的目是就的确定汽、液相组成（xi,yi）及汽相分率闪蒸计算输入了N+1个强度性质(增加一个强度性质是为了计算方便)，输出结果中，除了两相组成之外，还有另外一个性质即汽相分率η，共有2N+1个未知量，它们是从汽液平衡准则和物料平衡方程及归一化方程组成的2N+1个方程组联列求解比较泡点计算、露点计算和闪蒸计算可知，在泡点时，液相组成等于总组成，汽相分率等于0；在露点时，汽相组成等于总组成，汽相分率等于1；闪蒸时，汽、液相组与总组成均不相等，汽相分率在0和1之间4状态方程法（EOS法）计算混合物的汽液平衡状态方程法计算混合物汽液平衡的主要步骤如下：(1)选定的一个能适用于汽、液两相的状态方程，并结合混合法则推导出组分逸度系数的表达式（能用于汽、液两相的组分逸度的计算）；(2)由纯组分的有关参数（如临界性质、偏心因子等）得到各纯组分的状态方程常数，并得到相互作用参数（在混合法则中）；(3)由迭代法求解汽液平衡准则或组成归一化方程组。通常要借助计算机完成5关于相互作用参数通过计算目标函数的极小值得到相互作用参数11(1,2,,)()(1)(1,2,,)()11汽液平衡准则物料平衡或（归一化方程）vliiiiiNNiiiiffiNzxyiNxyexp11111explnlnln逸度：vNPNNPNvliiiiljijijiijyOBffx2exp11组成：NPNicalijijOByyexp11exp13压力或沸点：或NPNPcalcaljjjjpOBTTp423组合：（是加权因子）OBOBOB96状态方程+活度系数法（EOS+法）计算混合物的汽液平衡EOS+法分别采用两个模型来分别计算汽相和液相的组分逸度。若采用对称归一化定义的活度系数，则平衡准则可以转化为1）系统条件下纯组分逸度的计算根据等温条件下纯组分逸度随压力的变化其中称为Poynting因子2）平衡条件3）简化式3.1）对于低压下的理想溶液（理想系统）汽液平衡准则可以简化成3.2）对于中等压力条件下的体系近似取Poynting因子常用的EOS+法的相平衡准则为3.3）常、减压条件下的汽液平衡这是一类最常见的汽液平衡，常将汽相作为理想气体，液相作为非理想溶液处理，即则平衡准则为4）计算方法(1,2,,)vliiiiiypfxiNlnTfVpRT()lnd2得siplllsiiiislipfVVpppfRTRT()exp若不是很高时，有所以sslsssiiiiilslssiiiiiipfppVppfppRTlsiiiVppRT()expvsiiiiiiyppxiN(1,2,,)1i(1,2,,)siiiypxpiN1ivsiiiiiypxpiN(1,2,,)11viisiiiiyppxiN(1,2,,)11vii1i10常、减压汽液平衡计算中只涉及一个活度系数模型和一个计算蒸汽压的方程，不再需要状态方程。应用的计算软件中，包括了该类型汽液平衡计算的内容，其中所有的活度系数模型是Wilson方程，饱和蒸汽压用Antoine方程计算，计算类型分为四种。等压泡点的计算框图见P120图5-8若进行等温泡点的计算，可利用解析式求出气相组成和泡点压力7低压气体在液体中的溶解度气体在液体中的溶解度属于汽液平衡的一种特殊情况。由于在溶液状态下，混合物中的轻组分不一定能以纯的液态存在（混合物温度可能超过气体组分的临界温度），故将这种溶解平衡称为气液平衡（GLE）。轻组分处于超临界状态，故在液相的溶解度很低（常称为溶质），此时，溶质组分采用不对称归一化定义的活度系数更合理，所以溶质组分（1）的汽液平衡准则为溶剂组分（2）没有超临界，仍采用对称归一化定义的活度系数，其汽液平衡准则当系统的压力较低时，气相