蒸发和冷凝

Simcenter STAR-CCM+ 对液膜蒸发成气相和从气相冷凝成液膜的过程进行建模。

质量守恒

对多组分液膜的蒸发建模时,每个组分的蒸气压力及其蒸发率取决于混合物中不同组分的浓度。

每个组分 i 的组分质量通量在气体与液膜之间的交界面处守恒,因此:
1. EQUATION_DISPLAY
ρYi(v-h˙)-ρDidYidy=ρfYf,i(vf-h˙)-ρfDf,idYidy|f
(2770)

其中,在交界面处计算:

  • ρ ρf 分别为气体和液膜的密度
  • Yi Yf,i 分别为气体和液膜的质量分数
  • v vf 分别为气体和液膜的法向速度分量
  • Di Df,i 分别为气体和液膜的分子扩散系数
  • h˙ 为液膜厚度的变化率

总质量通量也守恒:

2. EQUATION_DISPLAY
ρ(v-h˙)=ρf(vf-h˙)
(2771)

Eqn. (2770)Eqn. (2771) 相结合将得到蒸汽质量分数 Yi 的以下微分方程:

3. EQUATION_DISPLAY
jNVρDidYi(y)dy=m˙v(jNVYi(y)1)
(2772)

其中, NV 为相互作用组分数,蒸发率定义如下:

4. EQUATION_DISPLAY
m˙v=-ρfh˙
(2773)

请注意,Eqn. (2770) 仅在低于饱和条件时才有效。

法向导数通过组分传递系数 st,i 和 Spalding 传递数 B 进行处理,因此:

5. EQUATION_DISPLAY
m˙v=-jNvst,j(Yc,j-Yj)1-jNvYjln(1+B)B
(2774)

Spalding 数 B 定义如下:

6. EQUATION_DISPLAY
B=jNvYj-jNvY1-jNvYj
(2775)

下标 c 表示网格单元值,并且可以假设 YYc

Eqn. (2774) 中的交界面气体质量分数 Yj Yf,j 通过背景摩尔权重 Wbg 相关,而交界面蒸汽质量分数 Yi

7. EQUATION_DISPLAY
Wbg=jNP,GXjWjjNP,GXj
(2776)
8. EQUATION_DISPLAY
YiXiWiXbgWbg+jNvXjWj
(2777)

每个组分的蒸发取决于其他组分,因为其蒸汽压力是混合物中这些组分浓度的函数。液体中组分 i 的蒸汽压力 pi 表示为:

9. EQUATION_DISPLAY
pi=γiXf,ipi*(T)
(2778)
其中:
  • γi 为活性系数,会考虑混合物中不同组分之间的相互作用。Simcenter STAR-CCM+ 提供两个模型来计算活性系数:Raoult 定律和修正 UNIFAC 模型。
  • pi* 为纯组分 i 的蒸汽压力。
  • Xf,i 为液体混合物中组分 i 的摩尔分数。

气体侧摩尔分数可以计算为:

10. EQUATION_DISPLAY
Xi=pip
(2779)

其中, p 为环境压力。

Raoult 定律
Raoult 定律假设液膜组分是具有类似分子结构的理想混合物。对于所有组分,Raoult 定律的活性系数约为 1。
修正 UNIFAC
修正 UNIFAC 模型是 UNIFAC(UNIQUAC 函数组活性系数)方法的自适应版本。UNIFAC 数据库 [631] 包含两个表,即由结构组列出的表面积和体积贡献以及不同组之间的能量相互作用参数。修正 UNIFAC 模型 [621] 将分子视为函数组的集合,并假设可以通过对组贡献求和来计算某些热力学属性。使用此方法,组分 i 的分子活性系数为 Eqn. (2778) 分为两部分:
11. EQUATION_DISPLAY
ln(γi)=ln(γiC)+ln(γiR)
(2780)
其中:
  • γiC 是表示由于分子大小差异而产生的贡献的组合项。
  • γiR 是表示由于分子相互作用的差异而产生的贡献的残差项。
残差项
混合物中组 m 中的 Xm 定义为:
12. EQUATION_DISPLAY
XmjNcXjvmjjNcXjkNg,jvkj
(2781)
其中:
  • Ng,j 为组分 j 中的组数。
  • Ng 为混合物中的总组数。
  • Nc 为混合物中的总组分数。
  • vki 为组分 j 中的结构组 k 数。
  • Xj 为混合物中组分 j 的摩尔分数。
混合物中组 m 的相对表面积定义为:
13. EQUATION_DISPLAY
Θm=XmQmNgnXnQn
(2782)

根据修正 UNIFAC 模型,组相互作用参数 Ψmn 计算如下:

14. EQUATION_DISPLAY
Ψmn=exp(amn+bmnT+cmnT2T)
(2783)
其中:
  • T 为温度。
  • amn bmn cmn 是从数据库中获取的组 mn 之间的组相互作用参数。
每个组 k 对活性系数的残差贡献计算如下:
15. EQUATION_DISPLAY
lnΓk=Qk{1lnmNgΘmΨmkmNgΘmΨkmmNgΘnΨnm}
(2784)
其中 ln(Γk) 通过假设分量以纯形式存在( Xj=1 )来计算。

对活性系数的残余贡献计算如下:

16. EQUATION_DISPLAY
ln(γiR)=kNg,ivki(ln(Γk)ln(Γki))
(2785)
组合项
分子尺寸差异产生的组合贡献可从以下来源获得:
17. EQUATION_DISPLAY
ln(γic)=1Vi+ln(Vi)z2qi[1ViFi+ln(ViFi)]
(2786)

组分 i 的表面积 qi 和体积 ri 估算如下:

18. EQUATION_DISPLAY
qi=kNg,iQkvki
(2787)
19. EQUATION_DISPLAY
ri=kNg,iRkvki
(2788)
其中 Rk Qk 表示从 UNIFAC 数据库中获取的组 k 体积和表面积贡献。以下参数定义为:
20. EQUATION_DISPLAY
Vi=rijNcXjrj
(2789)
21. EQUATION_DISPLAY
Vi=ri34jNcXjrj34
(2790)
22. EQUATION_DISPLAY
Fi=qijNcXjqj
(2791)

其中 z 为分子配位数,在文献中通常定义为常数值 10。

由于假设液膜较薄,因此认为气体侧的质量传递阻力明显大于液体侧的质量传递阻力。因此,液膜交界面摩尔分数浓度近似跟随关联的网格单元中心值。

使用活性系数、背景摩尔分数和液膜中的恒定组分分布,组分蒸发率可以表示如下:

23. EQUATION_DISPLAY
m˙v,i=Yim˙vρDidYidy
(2792)

这在所有条件下均有效。

除了流体动力学效应之外,蒸发率还可能受热效应影响。

能量守恒
交界面热通量守恒表示为:
24. EQUATION_DISPLAY
kdTdy-kfdTdy|fQv=0
(2793)

其中, k kf 表示气体和液膜导热率:

25. EQUATION_DISPLAY
Qv=jNvΔHivapmv,i˙
(2794)

这是交界面表面温度 Ts 的条件之一。

Eqn. (2792)Eqn. (2794) 相结合将得到总蒸发率表达式:

26. EQUATION_DISPLAY
m˙v=Q˙v+iNvΔHivapρDidYidyiNvΔHivapYi
(2795)

该表达式同时考虑了热效应,在所有条件下均有效。