RNG K-Epsilon 湍流

RNG K-Epsilon 模型是双方程湍流模型,它可对湍动能 k 和湍流耗散率 ε 的输运方程进行求解,从而封闭雷诺平均纳维-斯托克斯方程。

Yakhot 和其他人 [818] 将一个名为“重新归一化组”(RNG) 理论的统计方法应用于纳维-斯托克斯方程。RNG 理论适用于以下事实,即不同长度尺度的涡对湍流有贡献。它在计算耗散时从全局角度考虑了这些不同的尺度,而不是依赖一个湍流尺度。

最初提出的 RNG 模型未显式考虑可压缩性或浮力效应。但是,Simcenter STAR-CCM+ 中的标准 K-Epsilon 模型对这些效应进行了建模。

湍流粘度关系

湍流涡粘度 μt 计算如下:

1. EQUATION_DISPLAY
μt=ρCμkT
(4060)

其中:

湍流时间尺度 T 计算如下:

T T 带可实现尺度选项
2. EQUATION_DISPLAY
max(Te,Ctνε)
(4061)
3. EQUATION_DISPLAY
max(min(Te,CTCμfμS),Ctνε)
(4062)

其中:

传输方程

动能 k 和湍流耗散率 ε 的传输方程为:

4. EQUATION_DISPLAY
t(ρk)+∇⋅(ρkv¯)=∇⋅[(μ+μtσk)k]+Pkρ(εε0)+Sk
(4063)
5. EQUATION_DISPLAY
t(ρε)+∇⋅(ρεv¯)=∇⋅[(μ+μtσε)ε]+1TeCε1PεCε2ρ(εTeε0T0)Cε4ρε∇⋅v¯+SRNG+Sε
(4064)

其中:

  • v¯ 为平均速度。
  • μ 为动力粘度。
  • σk σε Cε1 Cε2 Cε4 模型系数
  • ε0 为源项中抵消湍流衰减 [316] 的环境湍流值。可以施加环境源项还会导致 Eqn. (1170) 给出的单位时间尺度 T0 的定义如下。
  • Pk Pε 结果项
  • Sk 和 为用户指定的源项。 Sε

SRNG Eqn. (4064) 中,表示平均流变形对湍流的影响,并定义如下:

6. EQUATION_DISPLAY
SRNG=Cμη3(1η/η0)1+βη3ρε2k
(4065)

其中:

结果项

结果项 Pk Pε 定义如下:

7. EQUATION_DISPLAY
Pk=Gk+Gnl+GbϒM
(4066)
8. EQUATION_DISPLAY
Pε=Gk+Gnl+Cε3Gb
(4067)

其中:

模型系数

CM (Sarkar) Ct CT Cε1 Cε2 Cε4 Cμ β σε σk
2 1 0.6 1.42 1.68 0.387 0.085 0.012 0.719 0.719
Cε3

可用文献未明确指定此系数。默认情况下,此系数根据 [305] 计算为:

9. EQUATION_DISPLAY
Cε3=tanh|vb||ub|
(4068)

其中, vb ub 为与重力矢量 g 平行和垂直的速度分量。

此公式会在自然对流边界层外将该系数设为零。

另外, Cε3 可以在任何位置用作常数,或根据浮力结果项 Gb 来指定,如下所示:

10. EQUATION_DISPLAY
Cε3={1   for Gb00   for Gb<0
(4069)