气液平衡的计算方法

合肥学院HefeiUniversity《化工热力学》过程论文题目：气液平衡的计算方法系别：化学与材料工程系专业：化学工程与工艺学号：1303021001姓名：教师：于晓飞气液平衡的计算方法摘要:气液平衡计算是化学过程中一项十分重要的计算。气液平衡的计算方法有几种，活度系数法，状态方程法（EOS法），GEMC和GDI方法计算流体气液相平衡。在气液平衡的计算中有三种泡点计算、露点计算和闪蒸计算,这里我们对闪蒸计算不做研究。关键词：气液平衡计算方法GEMCGDI正文:气液平衡计算的基本公式及计算类型：相平衡的判据应用于气液平衡，即为：广=fL（i=l，2,3,...,N）式中,f为混合物中组分i的逸度；上标V指的iii是气相；上标L指的是液相。上式既是气液平衡的准则，有事气液平衡计算的基本公式。具体应用时,需要建立混合物中组分的逸度fV、fL与体系的温度、压力以及气液相平衡组成ii关系.1.1活度系数法根据溶液热力学力论，将液相中组分的逸度与组分的活度系数相联系,简称活度系数法.对液相，由活度与活度系数的定义式得出广=YXf0式中，f0为标准态的逸度，以取Lewis-Randall定则为基准的标准态，即纯液体ii在体系的温度下的逸度。f0=fL=pS"Sexp'怜⑪i式中，指数项expJPi-dp称为Poynting因子，其意义是压力对f0影响的校正。psRT ii对气相将fL与fV表达式带入式中，得iii=l，2，.，N)y=YXps0Sexpi=l，2，.，N)iiiiiRTps式中,y和「别为汽、液相中组分i的摩尔分数心为气相混合物中组分i在体系温度T，体系压力p下的逸度系数；Y为液相中组分i的活度系数；ps为ii纯组分i在体系温度T时的饱和蒸气压；0s为纯组分i在体系温度T与其饱和蒸气压pS时的逸度系数；Vl为纯组分i在体系温度T时液相的摩尔体积.ii1.2GEMC方法计算原理GEMC方法可同时在两个模拟盒子中进行蒙特卡罗(MC)模拟，二者相对独立,但保持热力学相关，即满足相平衡条件(压力、温度和化学势相等),其温度T、总体积V和两个盒子中的总粒子数N保持不变。为达到相平衡，在模拟过程中需要进行3种不同的蒙特卡罗移动,按不同的接受概率接受，以满足相平衡条件：两个模拟盒的粒子分别在盒内自由移动,包括粒子的平动、转动等,以达到盒内平衡，其接受概率为P=min[l,exp(-AU/kT)](1)在保持总体积不变的条件下,在两个模拟盒子间进行体积的涨落,以达到两模拟盒压力相等，其接受概率为P=min1,exp(-(P=min1,exp(-(AU+WKT+Nln守lV-AV+N2ln-^V2)保持粒子总数不变的条件下,在两个模拟盒子间进行粒子交换，以达到两模拟盒化学势相等,接受概率为P=min1,NV12exp(-(P=min1,NV12exp(-(AU+AU)/KT)123)式中,U、V和N分别为两盒子的能量、体积和粒子数,T为体系温度，k为玻尔兹曼常数.1.3GDI方法计算原理流体的蒸发焓与饱和蒸气压可通过Clapeyron方程关联起来dp_AH

~dT~TAV式中，p、T分别为体系压力和温度，AHv为蒸发焓，AVv为两相的体积差•对含气体的相平衡，(4)式一般变形为(5)dlnpAH(5)_pTAVV式右边气液两相的蒸发焓及体积差均可由分子模拟求得，进而式转变为压力对温度的常微分方程,可通过数值法求解。因此在已知体系气液共存线上一个点作为积分起点的情况下，整个气液共存线可通过分子模拟并结合热力学积分计算得到。本工作采用三种分子的沸点作为积分参考点,液体和气体分别采用分子动力学和蒙特卡罗方法模拟,通过预测-矫正法计算积分，最终计算气液共存线。1。4Q函数法(间接法)汽液平衡时，按判据式(1-6.13),f(v)_f(L)(k=1,…,K),如气相米用逸度因

kk子、液相采用第I种活度因子分别计算气液相的非理想性，得

TOC\o"1-5"\h\z\o"CurrentDocument"pyp=p*p*xyexp[V*l（p一p*）/RT],k=1,2,,Kkkkkkk,I m,k k（2-2.1）整理上式可得系统总压p,p=Spy=Ep*p*xyexp[V*L（p一p*）/RT]/pk kkkk,I m,k k kk=1 k=12—2。2）式中y用式（1-7.46、7.47）代入,得k,IEp*p*xexp[V*l（p一p*Ep*p*xexp[V*l（p一p*）/RT]kk~~k rn^k kexp<Q+-KBx2-2。3）一艺xj=1jWj丿x[j,K]x[k,K]一VE心'

RT]〔也丿Lkx[k,K]jiSxj7 ' j丿x[j,K]Sp、-艺xjSx注意当i=K，式中对x的偏导数全为零。式（2-2。3）的意义在于：如果暂时不考K虑p*、p*、V*l、p、He和Ve，则式中除了Q以外，其它的变量就是已输入的kk m,k kmmT、p、x。而Q函数正是T、p、x的函数，式（2-2。3）实质上是一个Q函数的偏微分方程，只要有足够数量的一系列T、p、x的实验数据，原则上可以解得Q=Q（T,p,x）•有了Q,可用式（1-7。46、47）计算y，代入式（2—2。1）即可求k，I得y.至于那些暂时放在一边的变量：其中p*、p*和V*l是纯组分性质，与混合kk m,k物无关。p决定于气相组成y,可利用上次迭代的y值计算，但还需要使用合适k的状态方程，从这个意义上说，T、p、x推算y并不是完全的无模型，但当压力不太高时,气相非理想性远没有液相的那样强烈，在压力较低时，采用截止到第二维里系数的维里方程足以估算这种非理想性，甚至可以令p=1，也不致带严重k误差。至于He和Ve，后者很小，常可忽略,前者对于恒温数据不起作用，对于mm恒压数据，实践证明略去后影响不大。总之，这一方法基本上不使用模型,或者严格地说，不使用液相活度因子模型，而它是气液平衡计算中最关键的模型。式（2-2.3）原则可以求解,但实践上却有很大困难,因为导数出现在exp中,是一个超越型的偏微分方程，没有解析解,只能通过数值方法求解。国内外学者已发展

了多种方法，根据所采用数值方法的不同,可以分为几种类型：最早是Barker]提出的方法，其核心是选择一个过量函数模型代入式（2—2.3）,利用一系列T、p、x的实验数据，拟合得到模型参数和Q函数.第二种是Mixon等]发展的有限差分法，它以差分来逼近式（2-2。3）中的导数，然后利用Newton法迭代求得离散格点上的Q值•这种方法不依赖于任何过量函数模型，是严格的无模型法•它对二元系的计算非常成功，得到广泛的应用。但用于三元系时,收敛速度极慢，且求解过程不稳定。第三种是样条函数法，包括适用于二元系的三次样条函数法和适用于任意组分数的曲面样条函数法。特别是曲面样条函数法,它不仅能方便地用于二元系和三元系，也能成功地应用于多元系，更重要的是不同组分数的计算方法可以统一在一个框架下.大量实例计算表明，没有收敛的困难,不受多元系Q函数曲面类型的限制。下面介绍计算简单的Barker法和适用于多元系的曲面样条函数法，至于其它方法,感兴趣的读者可以参考相关的文献和著作。1.5Barker法Barker法的核心是选择一个合适的过量函数模型,其中包括若干待定模型参数,代入式（2-2.3）后，利用一系列T、p、x的实验数据,拟合得到模型参数，这就得到Q函数。原则上，第4章中介绍的各种过量函数或活度因子模型均可以作为Barker法的候选模型•但这些模型都是针对特定的对象而建立起来的，都有一定的适用范围，这就使得Barker法的准确度受到所选模型可靠性和适用性的限制,对多元系问题会更突出。解决的办法是尽可能选用灵活性大的经验模型，例如对于二元系,可以采用如下的Redlich—Kister型的经验Q函数模型，Q=xxEna（x一x）j1 2j=oj1 2（2—2。4）相对应的活Iny =x21,I2相对应的活Iny =x21,I2度因子为]NA（x一x）j+2xxj=0j1 2 1 2-2。5）NjA（x一x）jjj=0j1 2lny =x22,I1nA（x-x）j12xxj=0j2 1 12NjA（x一x）j-Jj=0j2 12-2.6）根据T、p、x实验数据的多少和计算精度的要求，可以选择不同的N值。式(2-2.4)是一个关于A的线性方程，如果已知不同组成下的Q函数值，j可以非常方便地采用最小二乘法关联得到N+1个A。具体计算时，我们可以采用j如下的迭代过程:(1)假设气相为理想气体、液相为理想溶液,计算气相组成的初值y。；⑵由状态方程计算气相逸度因子申*和申，由式(2—2.7)计算各组分的液kkexp[V*lexp[V*l(p-p*)/RT]}m,k kY=pyP/k,I kk(2—2.7)由式(2—2。8)计算各实验点的过量吉氏函数Q;Q=fxInykk,Ik=1(2—2.8)由计算得到的过量吉氏函数Q关联式(2—2.4)中的未知参数A(j=0，…，N)；j由式(2—2.5、2。6)计算各组分的活度因子丫；k,I由式(2-2.9)计算各组分新的气相组成yi;y=p*p*xYexp|V*l(p-p*)/RT]/pp，k=1, ,Kkkkkk,I m,k k k(2-2.9)比较yi和y。,如果两者不相等，则令yo=yi,转步骤(2)，进行新一轮循环迭代，直至达到规定的计算进度。参考文献|1]郭天民，多元汽液平衡和精馏，化学工业出版社，1983|2]朱自强，姚善泾，金彰礼，流体相平衡原理及其应用，浙江大学出版社，1990|3]PrausnitzJ.M.，LichthenthalerR。N.，deAzevedoE.G.,Molecularthermodynamicsoffluid—phaseequilibria，3rded.,PrenticeHallPTR,1999[4]AssaelM.J.，TruslerJ.P。M.，TsolakisT。F.，Thermophysicalpropertiesoffluids,ImperialCollegePress，London,1996[5]PrausnitzJ.M。,AndersonT。F。,Computercalculationsforthemulticomponentvapor—liquidandliquid—liquidequilibria,Prentice-HallInc。,EnglewoodCliffs.1980[6]MalanowskiS。，FluidPhaseEquilibria，8，197(1982)[7]MalanowskiS.,FluidPhaseEquilibria,9,311(1982)[6]WeirR。D.，deLoosThW.,Measurementofthethermodynamicpropertiesofmultiplephases，Elsevier，Amsterdam,2005HalaE。,PickJ。,FriedV。，VilimO.，Vapour—liquidequilibria,2ndEd。,Pergamon,Oxford，1967BarkerJ。A。，Austr.J.Chem.,6,207(1953)AbbottM。M.，vanNessH.C。,FluidPha