RANS 湍流传热

对于 RANS 湍流模型,能量方程中平均热通量的定义基于 Boussinesq 逼近法。对于 SKE LRe 模型,Simcenter STAR-CCM+ 提供温度通量模型,该模型将 Boussinesq 近似法替换为湍流热通量的代数公式。

默认情况下,平均热通量 在能量方程 Eqn. (1145) 中假设与湍流涡粘度成比例,即:q¯

1. EQUATION_DISPLAY
q¯=-(κ+μtCpPrt)T¯
(1343)

其中:

  • κ 为流体的导热率。
  • μt 为湍流涡粘度,由各自湍流模型给出。对于 RST 模型, 已重新定义,由 Eqn. (1340) 给出。μt
  • Cp 为比热。
  • Prt 为湍流普朗特数。
  • T¯ 为平均温度。

但是,对于浮力主导的情况或非常靠近壁面的位置,此假设不适用。由 Kenjeres 等人 [346] 提出的补救方法为用湍流热通量本身的代数公式替代 Boussinesq 近似。此公式是雷诺应力各向异性和温度变化的函数,需要为其求解额外的传输方程。代数温度通量模型的性能与近壁湍流行为的正确近似有很大关系,因此需要低雷诺数湍流模型。

对于温度通量模型,热通量定义为:

2. EQUATION_DISPLAY
q¯=-κT¯-ρCpvθ¯
(1344)

其中,ρ 为密度。

湍流热通量 的代数公式为:vθ_

3. EQUATION_DISPLAY
vθ_=-Ctu0kε(Ctu1RT¯+Ctu2v¯vθ¯+Ctu3βθ2_g)+Ctu4Avθ¯
(1345)

其中:

  • Ctu0、、、 和 为模型系数Ctu1Ctu2Ctu3Ctu4
  • k 为湍动能。
  • ε 为湍流耗散率。
  • REqn. (1308) 给出的雷诺应力张量。
  • β 为热膨胀系数。
  • g 为重力矢量。

A 为雷诺应力各向异性张量,定义为:

4. EQUATION_DISPLAY
A=1kR23I
(1346)

其中,I 为单位张量。

温度变化 通过求解附加传输方程来计算:θ2¯

5. EQUATION_DISPLAY
t(ρθ2¯)+∇⋅(ρθ2¯v¯)=∇⋅[(μ+μtσθ2)θ2¯]+Gθρεθ
(1347)

其中:

温度变化 的结果定义为:Gθ

6. EQUATION_DISPLAY
Gθ=-ρvθ¯T¯
(1348)

温度变化耗散率 从热时间尺度 的定义求得,即:εθTθ

7. EQUATION_DISPLAY
Tθ=θ2_2εθ
(1349)

且假设恒定的湍流-热时间尺度比:

8. EQUATION_DISPLAY
R=Tθ/T
(1350)

其中:

湍流热通量 还用于在 和 的传输方程中设置浮力结果项 ,即:vθ_kεGb

9. EQUATION_DISPLAY
Gb=-ρβ gvθ¯
(1351)

模型系数

R σθ2
0.5 1.0
Ctu0 Ctu1 Ctu2 Ctu3 Ctu4
0.12 1.0 0.6 0.6 1.5