多组分总体沸腾

多组分沸腾质量传递模型用于计算相交界面处的摩尔通量和温度。 计算相之间的质量传递、动量传递和能量传递时需要这些术语。

考虑固定坐标参考坐标系,并将 vi 定义为此坐标系中组分 i 的速度。 组分 i 的摩尔通量表示为:

1. EQUATION_DISPLAY
Ni=civi[molm2s]
(2074)

对所有组分求和以获得总摩尔通量:

2. EQUATION_DISPLAY
Nt=i=1nNi=ctv
(2075)

其中,摩尔平均速度 v 定义为:

3. EQUATION_DISPLAY
v=i=1nxivi
(2076)

现在,可以定义扩散通量(即相对于混合物通量的组分 i 的通量)。 选择 v 作为混合物参考速度,因此将摩尔扩散通量定义为:

4. EQUATION_DISPLAY
Ji=ci(vi-v)
(2077)

且:

5. EQUATION_DISPLAY
i=1nJi=0
(2078)

摩尔通量 Ni 通过以下公式与摩尔扩散通量相关:

6. EQUATION_DISPLAY
Ni=Ji+civ=Ji+xiNt
(2079)

薄膜模型

薄膜模型的假设汇总如下:

  • 质量传递的所有阻抗都位于与相边界相邻的薄膜中。
  • 此薄膜中的质量传递通过稳态分子扩散单独进行。
  • 通过薄膜的质量传递发生在与交界面垂直的方向(即忽略沿交界面的成分梯度)。
  • 薄膜外没有成分梯度。

扩散通量方程

在先前的假设下,可以使用有限差分计算跨薄膜的组分浓度梯度,扩散通量的方程变为:

7. EQUATION_DISPLAY
Ji=Ni-xiNt=ctj=1n1kij*Δxj
(2080)

其中,将有限通量质量传递系数矩阵表示为:

8. EQUATION_DISPLAY
[k*]=[k][Ξ]
(2081)

[k] [Ξ] 分别为低通量质量传递系数矩阵和有限通量校正系数矩阵。

零通量限制中的质量传递系数

低通量质量传递系数矩阵 [k] 定义为:

9. EQUATION_DISPLAY
[k]=1l[B]1
(2082)
10. EQUATION_DISPLAY
Bii=xiDin+k=1kinxkDij
(2083)
11. EQUATION_DISPLAY
Bij=xi(1Dij1Din)
(2084)

Eqn. (2082) 中显示的长度尺度 l 衍生自薄膜理论,是薄膜层的厚度。 此物理量可以由舍伍德数计算得出,因为:

12. EQUATION_DISPLAY
Sh=kdD=dl
(2085)

从中可得到 l=dSh

生成的液膜厚度是流动条件、几何和流体属性的函数。 在 Simcenter STAR-CCM+ 中,只能为每个相间相互作用指定一个特征长度 d 和一个舍伍德数 Sh 。 因此,为所有传递的组分计算并应用相同的 l 值。

质量传递校正系数矩阵

校正系数的形式取决于采用的质量传递模型。 对于两个液膜模型,遵循从蒸汽到液体的正通量约定,蒸汽和液体的校正矩阵分别为:

13. EQUATION_DISPLAY
[Ξ]V=[Φ]V[exp[Φ]V[I]]1
(2086)

14. EQUATION_DISPLAY
[Ξ]L=[Φ]Lexp[Φ]L[exp[Φ]L[I]]1
(2087)

其中:

15. EQUATION_DISPLAY
Φii=Nictκin+k=1kinNkctκik
(2088)
16. EQUATION_DISPLAY
Φij=Ni(1ctκij1ctκin)
(2089)

交界面能量平衡

施加跨汽-液交界面的能量通量的连续性:

17. EQUATION_DISPLAY
EV=EL=E
(2090)

其中, EV EL 为交界面处能量通量的法向分量。 Eqn. (2090) 可以重新写为:

18. EQUATION_DISPLAY
qV+i=1nNiVH¯iV(TI)=qL+i=1nNiLH¯iL(TI)
(2091)

其中热通量由以下公式给出:

19. EQUATION_DISPLAY
qV=hV*(TV-TI)
(2092)
20. EQUATION_DISPLAY
qL=hL*(TI-TL)
(2093)

根据薄膜模型,蒸汽相中的有限通量热传递系数计算如下:

21. EQUATION_DISPLAY
hV*=hV(ΦHV)/(expΦHV-1)=ΞHhV
(2094)

其中, hV 为零通量热传递系数, ΞH=(ΦHV)/(expΦHV-1) 为 Ackermann 校正因子,且

22. EQUATION_DISPLAY
ΦH=i=1nNiWicpi/h
(2095)

同样,对于液体:

23. EQUATION_DISPLAY
hL*=hL(ΦHL)expΦHL/(expΦHL-1)=ΞHLhL
(2096)

其中, ΞHL=(ΦHL)expΦHL/(expΦHL-1)

由于液相热传递系数的值较高, hL 的高通量校正预计接近于 1。

交界面处的组分平衡

假设两相之间交界面任一侧的相成分彼此处于即时局部平衡中以及与交界面温度和系统压力处于即时局部平衡中:

24. EQUATION_DISPLAY
yiI=Ki0(P,TI)xiI
(2097)

其中假定理想的混合物行为。 摩尔分数还遵循约束:

25. EQUATION_DISPLAY
xiI=1
(2098)
26. EQUATION_DISPLAY
yiI=1
(2099)

注意,仅当至少一个汽-液平衡系数 Ki0 小于 1 且至少一个 Ki0 大于 1 时,才可以求解这些交界面状态关系。

点求解器的方程

考虑交界面划分的两个混合物 x y 。 上标 L V 分别用于液体和蒸汽中的总物理量, I 表示交界面物理量。 交界面被视为没有质量传递阻抗的表面,平衡条件在任何时候都占优。 假设交界面处或薄膜中没有发生化学反应。

当交界面不固定时,通量 Ni 参考交界面,而不参考基准固定参考坐标系(有关更多详细信息,请参见 [554])。 假设在与交界面垂直的方向发生质量传递:

27. EQUATION_DISPLAY
Nix=Niy=Nii=1,2,...,n
(2100)
28. EQUATION_DISPLAY
i=1n
(2101)

其中, Nix Ni 的法向分量,从总体相 x 指向交界面 I ,而 Niy 从界面 I 指向总体相 y

一般情况下,液体中可能存在 m 组分,气体中可能存在 n ,它们并非都易挥发或可冷凝。 由于易挥发和可冷凝组分之间存在 k 耦合,因此 m-k 组分不易挥发, n-k 组分不可冷凝。

因此,要计算的变量为 m 交界面摩尔分数 xiI n 交界面摩尔分数 yiI k 交界面摩尔通量 Ni 和交界面温度 TI

因此,变量总数为 m+n+k+1 ,模型的方程为:

  • 基于交界面摩尔通量(液体侧)的 m-1 个方程
29. EQUATION_DISPLAY
Ni-xiINt=ctLi=1n1kil*L(xjI-xj)i=1,k+1,...,m1
(2102)
  • 基于交界面摩尔通量(蒸汽侧)的 n-1 个方程

30. EQUATION_DISPLAY
Ni-yiINt=ctVi=1n1kij*V(yjyjI)i=1,k+1,...,n1
(2103)
  • 基于交界面处能量平衡的 1 个方程
31. EQUATION_DISPLAY
hV*(TVTI)+i=1nNiVH¯iVWi(TI)=hL*(TITL)+i=1mNiLH¯iLWi(TI)
(2104)
  • 基于交界面处平衡关系的 k 个方程
32. EQUATION_DISPLAY
yiI=KixiIi=1,2,...,k
(2105)
  • 基于交界面(液体侧)处的摩尔分数约束的 1 个方程
33. EQUATION_DISPLAY
i=1mxiI=1
(2106)
  • 基于交界面(蒸汽侧)处的摩尔分数约束的 1 个方程
34. EQUATION_DISPLAY
i=1nyiI=1
(2107)

因此,方程数为 m+n+k+1 ,即变量数。 如果组分 i 不易挥发或不可冷凝,则 Eqn. (2102)Eqn. (2103) 中的 Ni=0 。 如果液体是纯物质 ( m=1 ),则不始终满足 Eqn. (2102)Eqn. (2105) 形式的方程。

同样,如果蒸汽是纯物质 ( n=1 ),则不始终满足 Eqn. (2103)Eqn. (2106) 形式的方程。

Eqn. (2106)Eqn. (2102) 是使用牛顿求解器为计算域的每个网格单元求解的非线性代数方程组。

稀释离散的简化模型

当 Maxwell-Stefan 扩散模型不活动时,进行稀释离散的假设,即 xi0 (对于 i=1,2,...,n1 )和 xn1 。 在这种假设下,仅考虑每个离散组分与列表中最后几个组分(溶剂)之间的相互作用。 忽略离散组分之间的相互作用。 因此,Eqn. (2102)Eqn. (2103) 分别简化为:

35. EQUATION_DISPLAY
Ni-xiINt=ctLκinLΦiiLeΦiiLeΦiiL1(xiI-xi)
(2108)

36. EQUATION_DISPLAY
Ni-yiINt=ctVκinVΦiiVeΦiiL1(yiyiI)
(2109)

其中, κin=Dinl Φii=Ntctκin

其他方程不受稀释离散假设的影响。

简称

下表中定义了多组分沸腾质量传递模型的方程中使用的项。

ci

的摩尔密度 i=ρi/Mi

[molm3]

ct

混合物摩尔密度 =i=1nci

[molm3]

cp

比热

[JkgK]

d 系统的特征长度尺度 m

|D|

Fick 扩散系数的矩阵

[m2s]

|D|

Maxwell-Stefan 扩散系数的矩阵

[m2s]

h

热传递系数

[Wm2K]

H

比焓

[Jkg]

|I| 单位张量

Ji

摩尔扩散通量 =ci(vi-v)

[molm2s]

l 液膜厚度 m

Mi

组分 i 的摩尔质量

[kgmol]

Ni

组分 i 的摩尔通量

[molm2s]

Sh 舍伍德数
v

摩尔平均参考速度 =i=1nxivi

[ms]

xi , yi 组分 i 的总摩尔分数(液体侧、蒸汽侧)
xiI , yiI 组分 i 的交界面摩尔分数(液体侧、蒸汽侧)

Wi

i 的分子量

[gmol]

μi

i 的摩尔化学势

[Jmol]

ρi

i的质量密度 =ciMi

[kgm3]

ρt

混合物质量密度 =i=1nρi

[kgm3]

ωi i的质量分数 =ρi/ρt
[Φ] 质量传递率因子的矩阵
[Ξ] 高通量校正因子的矩阵