「第5章 相平衡【复习准备】」.doc

第5章 相平衡重点难点：相平衡判据及相律，相平衡计算，活度系数模型中方程参数的确定。汽液平衡不同计算方法的特点与实际选择。1) 相平衡基础(1) 相平衡的判据平衡判据的确定应以热力学第二定律为依据。利用熵增原理可以来判断过程自发进行的方向和限度。由于温度、压力为实际过程容易测量与控制的参数，故采用Gibbs自由能的平衡判据在相平衡和化学平衡中应用最为广泛。在恒温恒压下，一切能自发进行的过程都将引起体系Gibbs自由能的减少。当体系达到平衡态时，体系的Gibbs自由能为最小。即恒温恒压下的封闭体系，平衡的判据可表达为 (dG)T,p=0 (5-1)根据式(5-1)，可推出相平衡的条件为：“各相的温度相等、压力相等、各组分在各相的化学位相等”。数学表达式为： (5-2)上式为以化学位表示的相平衡判据式。化学位在计算上不方便，故在解决相平衡问题时用组分逸度代替化学位，可得： (i=1,2,3,N) (5-3)此为以组分逸度表示的相平衡判据式，表明在一定T、p下处于平衡状态的多相多组分系统中，任一组分i在各相中的分逸度必定相等。 由于逸度比化学位在计算上更为方便，故式(5-3)是解决相平衡问题最实用的公式，是所有多组分体系相平衡计算的出发点。(2) 相律 相律是多组分多相平衡体系所遵循的最普遍规律，说明确定多组分多相平衡状态所需要的独立强度量的数目。Gibbs推导的相律表示如下： F(自由度)= C(组分数) - (相数)+2 (5-5)处于相平衡状态下，各相的温度、压力以及各相的组分的组成均已被确定，但描述体系的平衡状态无须使用全部的变量，只要由相律求得的自由度数的变量即可。例如二元系汽液平衡，组分数C= 2，相数=2，自由度F=2 。即只要在变量(温度T、压力p、气相组成yi、液相组成xi)中指定两个变量，则体系的状态便可确定。2) 互溶系统的汽液平衡计算通式相平衡的判据应用于汽液平衡(VLE)时，即 (i=1,2,3,N) (5-4)式中表示混合物中组分i的逸度，上标V指汽相，L指液相。上式既是汽液平衡的准则，又是汽液平衡计算的基本公式。N个组分共有N个式子，从数学的角度来看，求解条件已得到满足。然而具体应用时，需要建立混合物中组分的逸度与体系的温度、压力以及汽液相组成的关系，即汽液平衡模型。汽液平衡计算中常用的热力学处理方法(常用的汽液平衡模型)有两种：(1) 状态方程法(EOS法) 用一个能同时适用于汽液两相的状态方程来计算汽、液相的组分逸度系数。当汽液两相平衡时，有 (i=1,2,3,N) (5-5)上式即为多元系统相平衡计算的通式(EOS法)。式中、分别为汽、液相中组分i的逸度系数，可以用一个适合于汽、液两相的状态方程和混合法则来计算。 EOS法在高压、中压及低压下均适用，不必作不同程度的简化。近年来，EOS法已成为令人瞩目的方法，特别是高压汽液平衡，甚至超临界体系均可以方便地用EOS法进行汽液平衡的计算。(2) 活度系数法(i法)由状态方程计算汽相逸度系数，活度系数方程校正液相的非理想性。取以Lewis-Randall 规则为基础的标准态。当汽液两相平衡时，有 (i1,2,3,N) (5-6)上式即为多元系统相平衡计算的通式(活度系数法)。式中： (i=1,2,3,N) (5-7)这就是中低压下常用的汽液平衡计算通式。活度系数法不适用于高压汽液平衡的计算，因为活度系数方程中没有考虑p对i的影响。通常，需针对具体的汽液平衡系统，对式(5-7)做相应的简化处理。 在远离临界点区域的中低压范围内，可以忽略压力对液相的、的影响，Poynting因子近似为1，即。则可得中低压下汽液平衡的关系式为： (i1,2,3,N) (5-8) 低压下，汽相可视为理想气体，而液相为真实溶液，此时，则汽液平衡的关系式为： (i=1,2,3,N) (5-9)上式即为低压下非理想溶液的汽液平衡关系式，这是我们要讲的重点内容。 在中等压力下，汽相可以视为理想溶液，有，则汽液平衡的关系式为： (5-10) 中低压下体系中各组分性质相似，汽、液两相均可视为理想溶液，有，则汽液平衡的关系式为： 或 (5-11)式中fiV和fi L分别表示在体系温度、压力下纯气体i和纯液体i的逸度。此时，汽液平衡比(常数)Ki (5-12)可看出，对理想溶液，Ki仅为温度和压力的函数。对低分子的烃类组成的体系在远离临界点的压力下，其性质接近理想溶液，此时Ki 可用式(5-12)表示。De Priester早已制出了某些低分子烃类的p-T-K列线图。 当汽相为理想气体，液相为理想溶液完全的理想系，则 yi p= xi pis (i=1,2,N) (5-13) 注意：比较EOS法和i法的各自的优缺点。根据体系及条件选择合适的方法进行计算。3）汽液平衡的计算汽液平衡的计算是要得出体系处于平衡时的压力、温度及汽、液相组成之间的关系，典型的汽液平衡计算包括泡点、露点计算及闪蒸计算。 对N元体系，处于汽液平衡状态，总的变量数为2N个(温度T、压力p、N-1个液相摩尔分数xi，N-1个汽相摩尔分率yi)，其自由度F=2N-N=N，故对于N元体系，解决汽液平衡问题，首先应指定N个独立变量。当N个变量一经指定，则其余N个变量可通过列出N个组分的平衡关系式联立求解而得。(1) 中压下汽液平衡的计算多组分体系汽液平衡问题的求解，如对N元体系的汽液平衡，必须解N个方程。 (i1,2,3,N) (5-8)对于每个组分必须计算的热力学函数如下：热力学函数决定该函数的变量计算所需要的参数及方程pisTpis =f (T)中的所有常数 Antoine方程T, p, y1,yN-1Tci , pci, Vci. Zci,i求Bij virial方程isTpis, Tci, pci,i求Bi virial方程iT, x1,xN-1ViL, gij-gii, gji-gjj Wilson方程对中压下泡点、露点的计算，可以通过联解式(5-8)与组成归一化关系式求得。泡点计算：联解式(5-8)与露点计算：联解式(5-8)与用式(5-8)进行中压汽液平衡的计算，一般需要双重循环，适用于计算机求解。泡点温度、泡点压力、露点温度、露点压力的计算框图可见教材。(2) 活度系数方程中模型参数的确定(确定端值)实验数据确定活度系数的目的在于求活度系数方程式中的常数。根据相平衡方程式可以求得平衡组成，但在已知平衡组成的情况下也可以反过来求解活度系数，即通过实验测得各个相的组成时，根据相平衡关系式可以求出活度系数i。对真实溶液求活度系数i的目的是为了建立活度系数模型，包括两个步骤：第一，选择合适的活度系数模型；第二，确定模型中的参数。实验测定通常都是在低压或中压下进行的。在中、低压下，汽液平衡关系式为：将上式取对数，得 式中、可由状态方程计算，采用virial方程两项式计算即可。利用VLE关系式计算活度系数时，除了从实验测得的数据以外，还需确定第二virial系数。低压下，汽相可视为理想气体，则： (5-14)这是一个经常使用的公式，为了计算活度系数i 的值，必须要有T、p、xi、yi的实验数据。由实验数据测得的活度系数，通常用一个方程来关联，那么对一体系用什么方程关联比较合适呢？我们知道各种关联液相活度系数与组成之间的方程都是根据不同的溶液模型，通过超额Gibbs自由能与组成之间的关系而导出的，因此在选择合适的关联式时，可根据溶液超额Gibbs自由能与组成之间的关系形式作出初步的判断。 从汽液平衡组成求模型参数(多点法) 对仅含两个参数的关联式，原则上只要一个正确的实验点的数据。但这样的方法很可能由于该实验点的偶然误差而导致计算结果有很大的偏差，因此，一般需采用较多的实验点作为计算的依据，取平均值，这样才能获得比较可靠的计算结果。 对二元溶液，由汽液平衡数据算出不同液相组成时的活度系数1、2，再计算方程参数。如果有多个实验点，取平均值。 由无限稀释活度系数求模型参数 无限稀释活度系数是指混合物中的组分i在无限稀释条件下的活度系数，若用表示，则有对Margules 方程： 对van Laar 方程： 对Wilson 方程 可以看出，与活度系数关联式中的模型参数有着简单的关系，只要已知，的数据，就可算出活度系数关联式中方程的参数。可以通过实验确定无限稀释活度系数。 根据共沸组成求模型参数（单点法） 我们知道，一些非理想性较大的系统可以表现出共(恒)沸现象。混合物的共沸数据(二元系统有Taz、paz、xaz=yaz)反映了系统的非理想性，是汽液平衡数据的重要特殊点，故在测定上也尤为仔细，准确度较高，经常被用于求解活度系数的模型参数。特别是常减压条件下的共沸点数据已有较多的积累。若二元系在一定温度、压力下能产生共沸，则可根据共沸数据较方便地求模型参数。 设汽相为理想气体，在共沸条件下，xi=yi，则 (5-15)由此，可以根据共沸点的汽液平衡数据计算出共沸点的活度系数值，再结合具体的活度系数模型解出两个模型参数。(3) 如何判断体系能否形成共(恒)沸物？可以利用相对挥发度12是否等于1来判断体系有无共沸物生成。 对二元系 (1-16)对低压体系，汽相可以看成理想气体，则 (5-17)若存在共沸点，则12 =1。 计算x1=0时(12) x1=0和x1=1时(12) x1=1的值，判断 若(12) x1=01，而(12) x1=1 1 或 (12) x1=0 1，而(12) x1=11 则存在共沸点判断12 是否等于1：因为12是x1的连续函数，如果12从小于1变化到大于1，或者从大于1变化到小于1，中间必然经过12 =1的点，即在某一组成时一定存在共沸点。(4) 汽液平衡数据的热力学一致性检验汽液平衡实验数据是否可靠，可用Gibbs-Duhem方程来检验，若实验数据与Gibbs-Duhem方程符合，则称此汽液平衡数据具有热力学一致性。当温度、压力变化时，Gibbs-Duhem方程变为 (5-18)上式即为变温、变压下的Gibbs-Duhem方程的形式，是进行热力学一致性检验的基本工具，可直接用于汽液平衡数据的热力学一致性检验。在恒温、恒压下，Gibbs-Duhem方程变为对