Function for Hybrid Upwind Scheme

Recognizing the superiority of an upwind scheme for RANS computations and the benefits of central differencing for LES, Travin and others [365] proposed the following blending coefficient for DES:

Figure 1. EQUATION_DISPLAY
σHU=tanh{[min(A,10)]3}
(1450)

where:

Figure 2. EQUATION_DISPLAY
A=max(lratiomax(g,10-20)-0.5,0)
(1451)
Figure 3. EQUATION_DISPLAY
lratio=CDESΔρCμ1.5max(Sˆ,CdesT/τ)μ+μt
(1452)
Figure 4. EQUATION_DISPLAY
g=tanh{[min(B,10)]4}
(1453)
Figure 5. EQUATION_DISPLAY
B=2Wmax(S,W)max(Sˆ2,10-20)
(1454)
Figure 6. EQUATION_DISPLAY
Sˆ2=W2+S22
(1455)

with the modulus of the mean vorticity tensor W and the modulus of the mean strain rate tensor S defined according to Eqn. (1131) and Eqn. (1129) respectively.

The local convective scale, τ is defined as:

Figure 7. EQUATION_DISPLAY
τ=Δ|v|
(1456)

CdesT is a Model Coefficient and CDES is given in Spalart-Allmaras DES—Model Coefficients.

Model Coefficient

CdesT
1