混合物湍流

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

混合物物理量

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

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

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

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

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

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

3. EQUATION_DISPLAY
k m = k i
(2447)

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

4. EQUATION_DISPLAY
ε m = ε i
(2448)

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

5. EQUATION_DISPLAY
v m ¯ = i = 1 n α i ρ i v i i = 1 n α i ρ i
(2449)

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

6. EQUATION_DISPLAY
ν m t = ν i t
(2450)

因此:

7. EQUATION_DISPLAY
μ m t = ρ m μ i t ρ i
(2451)

且相湍流粘度为:

8. EQUATION_DISPLAY
μ i t = ρ i μ m t ρ m
(2452)

K-Epsilon 混合物湍流

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

9. EQUATION_DISPLAY
t ( ρ m k m ) + ∇⋅ ( ρ m k m v m ¯ ) = ∇⋅ [ ( μ m + μ m t σ k ) k m ] + P m k ρ m ( ε m ε 0 ) + S m k
(2453)
10. EQUATION_DISPLAY
t ( ρ m ε m ) + ∇⋅ ( ρ m ε m v m ¯ ) = ∇⋅ [ ( μ m + μ m t σ ε ) ε m ] + [ 1 t m e C ε 1 P m ε C ε 2 f 2 ρ m ( ε m t m e ε 0 t m 0 ) + S m ε ]
(2454)

其中:

  • t m e = k m ε m 为混合物大涡时间尺度
  • C ε 1 C ε 2 为模型系数
  • f 2 为阻尼函数。

K-Omega 混合物湍流模型

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

11. EQUATION_DISPLAY
t ( ρ m k m ) + ∇⋅ ( ρ m k m v m ¯ ) = ∇⋅ [ ( μ m + σ k μ m t ) k m ] + P m k ρ m β * f β * ( ω m k m ω 0 k 0 ) + S m k
(2455)
12. EQUATION_DISPLAY
t ( ρ m ω m ) + ∇⋅ ( ρ m ω m v m ¯ ) = ∇⋅ [ ( μ m + σ ω μ m t ) ω m ] + P m ω ρ m β f β ( ω m 2 ω 0 2 ) + S m ω
(2456)

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

混合物耗散率方程

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

13. EQUATION_DISPLAY
t ( ρ m ε m ) + ∇⋅ ( ρ m ε m v m ¯ ) = α m ∇⋅ [ ( μ m + μ m t σ ε ) ε m ] + ε m k m [ C ε 1 ( 1 2 t r ( P m ε ) + 1 2 C ε 3 t r ( G m ) ) C ε 2 ρ m ε m ] + i j ( m i j ε j ( i j ) m j i ε i )
(2457)

其中

  • P m ε 为混合物的湍流结果。
  • G m 为混合物的浮力结果。

请参见雷诺应力传输方程