Skip to main content

Spalart-Allmaras Detached Eddy Simulation (DES)

Properties

  • One equation model based on a modified turbulence viscosity, ν~\tilde{\nu}

Model equations

The ν~\tilde{\nu} transport equation is given by:

\Ddt{\rho \tilde{\nu}} = \div \left( \rho D_\tilde{\nu} \tilde{\nu} \right) + \frac{C_{b2}}{\sigma_{\nu_t}} \rho \mag{\grad \tilde{\nu}}^2 + C_{b1} \rho \tilde{S} \tilde{\nu} \left(1 - f_{t2}\right) - \left(C_{w1} f_w - \frac{C_{b1}}{\kappa^2} f_{t2}\right) \rho \frac{\tilde{\nu}^2}{\tilde{d}^2} + S_\tilde{\nu}

where the length scale d~\tilde{d} is defined by:

d~=min(ΨCDESΔ,y) \tilde{d} = \min \left( \Psi C_{DES} \Delta, y \right)

and Ψ\Psi is the low Reynolds number correction function:

Ψ2=min[102,11Cb1Cw1κ2fw[ft2+(1ft2)fv2]fv1max(1010,1ft2)] \Psi^2 = \min \left[ 10^2, \frac{1 - \frac{1 - C_{b1}}{C_{w1} \kappa^2 f_w^{*}} \left[ f_{t2} + \left(1 - f_{t2}\right) f_{v2} \right]}{f_{v1} \max \left(10^{-10}, 1-f_{t2} \right)} \right]

Default model coefficients

σνt\sigma_{\nu_t}Cb1C_{b1}Cb2C_{b2}Cw1C_{w1}Cw2C_{w2}Cw3C_{w3}
2/30.13550.622Cb1κ2+1+Cb2σνt\frac{C_{b1}}{\kappa^2} + \frac{1 + C_{b2}}{\sigma_{\nu_t}} 0.32
Cv1C_{v1}CsC_{s}CDESC_{\mathit{DES}}CkC_{k}Ct3C_{t3}Ct4C_{t4}fwf_w^{*}
7.1 | 0.3 | 0.65 | 0.07 | 1.2 | 0.5 | 0.424

Initialisation

Usage

The model is specified using:

LES
{
turbulence on;
LESModel SpalartAllmarasDES;

// Optional entries
SpalartAllmarasDESCoeffs
{
// Apply low-Reynolds number correction; default = yes
lowReCorrection yes;
}
}

Boundary conditions

Inlet

Outlet

Walls

  • wall functions

Further information

Source code:

References:

  • Principle reference: Spalart et al. 67
  • Low Reynolds number correction: Spalart et al. 68