ReTheta_t 相关性

R e θ t 既定的相关性由 Abu-Ghannam 和 Shaw [367] 提出。

在零压力梯度情况下,此相关性为:

Reθt=163+e(6.91-Tu)
(1513)

其中,Tu 为湍流强度。

由 Menter 等人 ([373]) 针对 Gamma ReTheta 转换模型提出的替换表达式为:

Reθt=803.73(Tu+0.6067)-1.027
(1514)

此表达式后经 Langtry [371] 修订,旨在改进较低自由流湍流强度时的结果,如下所示:

Reθt={(1173.51-589.428Tu+0.2196/Tu2);  Tu≤1.3331.5(Tu-0.5658)-0.671;  Tu>1.3
(1515)

下图中已绘制三种零压力梯度相关性。Simcenter STAR-CCM+ 中提供了所有这三种相关性,使用户能够灵活地自定义相关性,但Eqn. (1515) 是内部 "Suluksna-Juntasaro" 校准中默认使用的相关性。



Reθt vs. 自由流湍流强度

对于 Eqn. (1513) 至 Eqn. (1515) 中为零压力梯度呈现的三种相关性,最初均以包含压力梯度效应的形式呈现。有一点疑问:Reθt 相关性[378] 中需要压力梯度效应,因为它们可以引入压力梯度效应的双重计算。但是,它们以可选方式包括在 Simcenter STAR-CCM+ 中。如果存在压力梯度,相关性需要计算 Thwaites 参数 λθ,而在由 Menter 等人 [373] 提出的相关性中,需要计算加速度参数 K。

Thwaite 参数 λθ 定义为:

λθ=θ2ν∂U∂s 
(1516)

其中:

θ=νURe⁡θt
(1517)
U=|v¯|
(1518)
和
  • ν 为运动粘度。
  • v¯ 为平均速度矢量。

参数 K 定义为:

K=νU2∂U∂s 
(1519)

s 表示壁面平行方向,其中使用流向和壁面法向构造局部坐标系。

湍流强度从湍动能中得出,如下所示:

Tu=1002k/3U
(1520)

流向速度梯度计算方法如下:

 ∂U∂s=(∇v¯T⋅v¯)⋅v¯v¯⋅v¯
(1521)

其中,v¯T 表示 v¯ 的转置。

由于动量厚度 θ 必须从动量厚度雷诺数 Reθ 中衍生得出,因此相关性将成为迭代求解的隐式函数。

Abu-Ghannam 和 Shaw [367] 提出的非零压力梯度 (NZPG) 相关性为:

Reθt=163+e[F(λθ)(1-Tu6.91)]
(1522)

其中:

F(λθ)={6.91+12.75λθ+63.64λθ2;  λθ≤06.91+2.58λθ-12.27λθ2;  λθ>0
(1523)

为了确保数值稳定性,将应用以下限制:

-0.1≤λθ≤0.1Reθt≥20
(1524)

Menter 等人 ([373]) 提出的 NZPG 相关性为:

Reθt=803.73(Tu+0.6067)-1.027F(λθ,K)
(1525)
F(λθ,K)={1+Fλexp(-Tu/3);  λθ≤01+FK[1-exp(-Tu/1.5)]+0.556[1-e-23.9λθ]e-Tu/1.5;  λθ>0
(1526)

其中:

Fλ=10.32λθ+89.47λθ2+265.51λθ3
(1527)
FK=[0.0962(K∙106)+0.148(K∙106)2+0.0141(K∙106)3]
(1528)

为了确保数值稳定性,将应用以下限制:

-0.1≤λθ≤0.1-3x10-6≤K≤3x10-6Reθt≥20
(1529)

Langtry[371]提出的 NZPG 相关性为:

Reθt={(1173.51−589.428Tu+0.2196/Tu2)F(λθ)331.5(Tu−0.5658)−0.671F(λθ);Tu≤1.3;Tu>1.3
(1530)
F(λθ)={1+12.986λθ+123.66λθ2+405.689λθ3;  λθ≤01+0.275[1-exp(-35λθ)]exp(-2Tu);  λθ>0
(1531)

为了确保数值稳定性,将应用以下限制:

-0.1≤λθ≤0.1Tu≥0.027Reθt≥20
(1532)

下图对比显示了一系列自由流湍流强度的三种 NZPG 相关性。



Reθt vs. λθ 用于各种自由流湍流强度值