Correlations for ReTheta_t

A well-established correlation for R e θ t was developed by Abu-Ghannam and Shaw [374].

In the case of a zero pressure gradient, this correlation is:

Figure 1. EQUATION_DISPLAY
Reθt=163+e(6.91-Tu)
(1528)

where Tu is the turbulence intensity.

An alternative expression, developed for the Gamma ReTheta transition model by Menter et al. [380] is:

Figure 2. EQUATION_DISPLAY
Reθt=803.73(Tu+0.6067)-1.027
(1529)

and this was later revised by Langtry [378] with a view to improving the results at lower free-stream turbulence intensities, as follows:

Figure 3. EQUATION_DISPLAY
Reθt={(1173.51-589.428Tu+0.2196/Tu2);  Tu1.3331.5(Tu-0.5658)-0.671;  Tu>1.3
(1530)

The three zero-pressure gradient correlations are plotted in the following figure. All three are provided in Simcenter STAR-CCM+ to allow the user flexibility to customize the correlations, but Eqn. (1530) is the default correlation used in the internal “Suluksna-Juntasaro” calibration.



Reθt vs. Free-Stream Turbulence Intensity

The three correlations presented for zero pressure gradient in Eqn. (1528) through Eqn. (1530) were originally presented in a form that include the effects of a pressure gradient. There is some doubt that the pressure gradient effects are required in the Reθt correlation [385], since they arguably introduce a double-accounting of the pressure gradient effects, but they are nevertheless optionally included in Simcenter STAR-CCM+. In the presence of a pressure gradient, the correlations require the computation of the Thwaites parameter, λθ, and in the case of the Menter et al. [380] correlation, the acceleration parameter, K.

Thwaite’s parameter λθ is defined as:

Figure 4. EQUATION_DISPLAY
λθ=θ2νUs
(1531)

where:

Figure 5. EQUATION_DISPLAY
θ=νUReθt
(1532)
Figure 6. EQUATION_DISPLAY
U=|v¯|
(1533)
and:
  • ν is the kinematic viscosity.
  • v¯ is the mean velocity vector.

The parameter K is defined as:

Figure 7. EQUATION_DISPLAY
K=νU2Us
(1534)

s denotes the wall-parallel direction, where the local coordinate system is constructed using the flow direction and the wall-normal direction.

The turbulence intensity is obtained from the turbulent kinetic energy as follows:

Figure 8. EQUATION_DISPLAY
Tu=1002k/3U
(1535)

The streamwise velocity gradient is calculated as:

Figure 9. EQUATION_DISPLAY
Us=(v¯Tv¯)v¯v¯v¯
(1536)

where v¯T denotes the transpose of v¯.

Since the momentum thickness θ must be derived from the momentum thickness Reynolds number Reθ, the correlation becomes an implicit function that is solved iteratively.

The non-zero pressure gradient (NZPG) correlation of Abu-Ghannam and Shaw [374] is:

Figure 10. EQUATION_DISPLAY
Reθt=163+e[F(λθ)(1-Tu6.91)]
(1537)

where:

Figure 11. EQUATION_DISPLAY
F(λθ)={6.91+12.75λθ+63.64λθ2;  λθ06.91+2.58λθ-12.27λθ2;  λθ>0
(1538)

For numerical robustness, the following limits are applied:

Figure 12. EQUATION_DISPLAY
-0.1λθ0.1Reθt20
(1539)

The NZPG correlation of Menter et al. [380] is:

Figure 13. EQUATION_DISPLAY
Reθt=803.73(Tu+0.6067)-1.027F(λθ,K)
(1540)
Figure 14. EQUATION_DISPLAY
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
(1541)

where:

Figure 15. EQUATION_DISPLAY
Fλ=10.32λθ+89.47λθ2+265.51λθ3
(1542)
Figure 16. EQUATION_DISPLAY
FK=[0.0962(K106)+0.148(K106)2+0.0141(K106)3]
(1543)

For numerical robustness, the following limits are applied:

Figure 17. EQUATION_DISPLAY
-0.1λθ0.1-3x10-6K3x10-6Reθt20
(1544)

The NZPG correlation of Langtry [378] is:

Figure 18. EQUATION_DISPLAY
Reθt={(1173.51589.428Tu+0.2196/Tu2)F(λθ)331.5(Tu0.5658)0.671F(λθ);Tu1.3;Tu>1.3
(1545)
Figure 19. EQUATION_DISPLAY
F(λθ)={1+12.986λθ+123.66λθ2+405.689λθ3;  λθ01+0.275[1-exp(-35λθ)]exp(-2Tu);  λθ>0
(1546)

For numerical robustness, the following limits are applied:

Figure 20. EQUATION_DISPLAY
-0.1λθ0.1Tu0.027Reθt20
(1547)

A comparison of the three NZPG correlations is shown in the following figure for a range of free-stream turbulence intensities.



Reθt vs. λθ for various values of Free-Stream Turbulence Intensity