气液平衡的计算方法.docx

合肥学院

HefｅiＵniversｉty

《化工热力学》过程论文

题目:

气液平衡得计算方法

系别：

化学与材料工程系

专业:

化学工程与工艺

学号:

1303021０01

姓名:

于晓飞

教师:

高大明

气液平衡得计算方法

摘要:气液平衡计算就就是化学过程中一项十分重要得计算。气液平衡得计算方法有几种,活度系数法,状态方程法(EOS法)，ＧEMＣ和ＧDI方法计算流体气液相平衡。在气液平衡得计算中有三种泡点计算?、露点计算和闪蒸计算,这里我们对闪蒸计算不做研究。

关键词:气液平衡计算方法ＧEＭＣGDＩ

正文：

气液平衡计算得基本公式及计算类型:相平衡得判据应用于气液平衡,即为:（i=１,2,３,…,Ｎ)式中,为混合物中组分ｉ得逸度;上标V指得就就是气相;上标L指得就就是液相。

上式既就就是气液平衡得准则,有事气液平衡计算得基本公式。具体应用时，需要建立混合物中组分得逸度、与体系得温度、压力以及气液相平衡组成关系。

1、1活度系数法

根据溶液热力学力论,将液相中组分得逸度与组分得活度系数相联系，简称活度系数法。对液相,由活度与活度系数得定义式得出

=

式中,为标准态得逸度,以取Lewis-Ｒａnｄａll定则为基准得标准态,即纯液体i在体系得温度下得逸度。

==

式中，指数项称为Ｐoyntinｇ因子,其意义就就是压力对影响得校正。

对气相

将与表达式带入式中,得

p=（i＝1,2,…,N)

式中，和分别为汽、液相中组分i得摩尔分数;为气相混合物中组分ｉ在体系温度T,体系压力p下得逸度系数；为液相中组分i得活度系数;为纯组分ｉ在体系温度T时得饱和蒸气压;为纯组分i在体系温度T与其饱和蒸气压时得逸度系数;为纯组分i在体系温度T时液相得摩尔体积。

1、2GEMＣ方法计算原理

GEMＣ方法可同时在两个模拟盒子中进行蒙特卡罗(MＣ)模拟，二者相对独立,但保持热力学相关,即满足相平衡条件(压力、温度和化学势相等）,其温度T、总体积V和两个盒子中得总粒子数N保持不变、为达到相平衡,在模拟过程中需要进行3种不同得蒙特卡罗移动,按不同得接受概率接受，以满足相平衡条件：

1)两个模拟盒得粒子分别在盒内自由移动,包括粒子得平动、转动等，以达到盒内平衡,其接受概率为

Ｐ=min[1，exｐ(－ΔＵ／kＴ)](1)

2)在保持总体积不变得条件下,在两个模拟盒子间进行体积得涨落,以达到两模拟盒压力相等,其接受概率为

P＝min(2）

3)保持粒子总数不变得条件下,在两个模拟盒

子间进行粒子交换,以达到两模拟盒化学势相等,接受概率为

P=miｎ(３）

式中,U、V和N分别为两盒子得能量、体积和粒子

数,T为体系温度，k为玻尔兹曼常数、

1、3?GＤＩ方法计算原理

流体得蒸发焓与饱和蒸气压可通过Ｃlapｅyron方程关联起来:

(4)

式中,p、Ｔ分别为体系压力和温度，ΔHv为蒸发焓,ΔVv为两相得体积差、对含气体得相平衡，（４)式一般变形为

(5)

(５)式右边气液两相得蒸发焓及体积差均可由分子模拟求得,进而式转变为压力对温度得常微分方程,可通过数值法求解、因此在已知体系气液共存线上一个点作为积分起点得情况下，整个气液共存线可通过分子模拟并结合热力学积分计算得到、本工作采用三种分子得沸点作为积分参考点,液体和气体分别采用分子动力学和蒙特卡罗方法模拟，通过预测－矫正法计算积分,最终计算气液共存线。

1、4Q函数法（间接法)

汽液平衡时,按判据式(1-６、1３）,(k=１,…,Ｋ),如气相采用逸度因子、液相采用第I种活度因子分别计算气液相得非理想性,得

, ? ?(2－2、1)

整理上式可得系统总压p，

?? ??(2－2、２）

式中用式(1-７、46、７、4７)代入,得

? ?（2-2、3)

注意当ｉ=K,式中对ｘK得偏导数全为零。式（２－2、3)得意义在于:如果暂时不考虑、、、、和，则式中除了Q以外，其她得变量就就就是已输入得T、p、x。而Q函数正就就是T、ｐ、x得函数,式(2－２、3)实质上就就是一个Q函数得偏微分方程，只要有足够数量得一系列T、p、ｘ得实验数据,原则上可以解得Ｑ=Q(Ｔ,p，x)。有了Ｑ，可用式(１-7、46、47)计算?k,I,代入式(２-2、1)即可求得ｙ。至于那些暂时放在一边得变量:其中、和就就是纯组分性质,与混合物