混合物湍流

单个湍流模型用于计算所有相的湍流,方法为使用混合物特性和混合速度求解湍流物理量传输方程。

混合物物理量

混合物密度 ρm 定义为的所有其他 n 相物理量的线性插值(由体积分数 αi 加权):

1. EQUATION_DISPLAY
ρm=i=1nαiρi
(2445)

混合物动力粘度也是线性插值:

2. EQUATION_DISPLAY
μm=i=1nαiμi
(2446)

km εm 通过传输方程的解进行计算。

湍流动能是单位质量速度波动的能量,因此假设混合物 k 和相 k 相等:

3. EQUATION_DISPLAY
km=ki
(2447)

同样, εm 描述耗散率密度,并且假设对于混合物和相相等:

4. EQUATION_DISPLAY
εm=εi
(2448)

平均混合物速度被定义为:

5. EQUATION_DISPLAY
vm¯=i=1nαiρivii=1nαiρi
(2449)

混合物湍流涡粘度的计算假设各相的运动粘度相等:

6. EQUATION_DISPLAY
νmt=νit
(2450)

因此:

7. EQUATION_DISPLAY
μmt=ρmμitρi
(2451)

且相湍流粘度为:

8. EQUATION_DISPLAY
μit=ρiμmtρm
(2452)

K-Epsilon 混合物湍流

混合物湍动能 km 和混合物湍流耗散率 εm 的传输方程如下:

9. EQUATION_DISPLAY
t(ρmkm)+∇⋅(ρmkmvm¯)=∇⋅[(μm+μmtσk)km]+Pmkρm(εmε0)+Smk
(2453)
10. EQUATION_DISPLAY
t(ρmεm)+∇⋅(ρmεmvm¯)=∇⋅[(μm+μmtσε)εm]+[1tmeCε1PmεCε2f2ρm(εmtmeε0tm0)+Smε]
(2454)

其中:

  • tme=kmεm 为混合物大涡时间尺度
  • Cε1 Cε2 为模型系数
  • f2 为阻尼函数。

K-Omega 混合物湍流模型

混合物湍动能 km 和混合物比耗散率 ωm 的传输方程:

11. EQUATION_DISPLAY
t(ρmkm)+∇⋅(ρmkmvm¯)=∇⋅[(μm+σkμmt)km]+Pmkρmβ*fβ*(ωmkmω0k0)+Smk
(2455)
12. EQUATION_DISPLAY
t(ρmωm)+∇⋅(ρmωmvm¯)=∇⋅[(μm+σωμmt)ωm]+Pmωρmβfβ(ωm2ω02)+Smω
(2456)

雷诺应力传输混合物湍流模型

混合物耗散率方程

对于混合相 m ,耗散率方程类似于 K-Epsilon 混合物湍流模型:

13. EQUATION_DISPLAY
t(ρmεm)+∇⋅(ρmεmvm¯)=αm∇⋅[(μm+μmtσε)εm]+εmkm[Cε1(12tr(Pmε)+12Cε3tr(Gm))Cε2ρmεm]+ij(mijεj(ij)mjiεi)
(2457)

其中

  • Pmε 为混合物的湍流结果。
  • Gm 为混合物的浮力结果。

请参见雷诺应力传输方程