Two-Layer Models

A two-layer model solves for k but prescribes ε algebraically with distance from the wall in the viscosity dominated near-wall flow regions. In Simcenter STAR-CCM+, this approach is available with both the Standard and the Realizable K-Epsilon models.

For the two-layer models the dissipation rate near the wall is simply prescribed as:

Figure 1. EQUATION_DISPLAY
ε=k3/2lε
(1193)

where:

  • lε is a length scale function that is calculated depending on the Model Variant.

As per the suggestion of Jongen [307], the following wall-proximity indicator is used to combine the two-layer formulation with the full two-equation model:

Figure 2. EQUATION_DISPLAY
λ=12[1+tanh(Red-Rey*A)]
(1194)

where:

  • Red is the wall-distance Reynolds number given by Eqn. (1134).
  • Rey* defines the limit of applicability of the two-layer formulation and is a Model Coefficient.
  • A determines the width of the wall-proximity indicator and is defined such that the value of λ is within 1% of its far-field value for a given variation of ΔRey:
    Figure 3. EQUATION_DISPLAY
    A=|ΔRey|atanh0.98
    (1195)

    where ΔRey is a Model Coefficient.

The turbulent viscosity from the K-Epsilon model μt|kϵ is then blended with the two-layer value as follows:

Figure 4. EQUATION_DISPLAY
μt=λμt|kϵ+(1λ)μ(μtμ)2layer
(1196)

where:

Model Variants

In Simcenter STAR-CCM+ three variants of the two-layer formulation are available:

Model Variant lε (μtμ)2layer
Wolfstein [317]
Figure 5. EQUATION_DISPLAY
cld[1-exp (-Red2cl  )]
(1197)
Figure 6. EQUATION_DISPLAY
0.42RedCμ1/4[1-exp (-Red70 )]
(1198)
Norris-Reynolds [312]
Figure 7. EQUATION_DISPLAY
cldRedRed+5.3
(1199)
Figure 8. EQUATION_DISPLAY
0.42RedCμ1/4[1-exp (-Red50.5 )]
(1200)
Xu [318]

(specifically for natural convection)

Figure 9. EQUATION_DISPLAY
8.8d1+10/yv*+5.15×10-2yv*1v2_/k
(1201)
Figure 10. EQUATION_DISPLAY
0.544yv*1+5.025×104yv*1.65
(1202)

where:

Model Coefficients

Rey* ΔRey
60 10