Spalart–Allmaras turbulence model

This article has multiple issues. Please help or discuss these issues on the talk page. (Learn how and when to remove these template messages) This article relies largely or entirely on a single source. Relevant discussion may be found on the talk page. Please help by introducing citations to additional sources. Find sources:  –  ·  ·  · scholar · JSTOR (April 2017) This section does not cite any sources. Please help by adding citations to reliable sources. Unsourced material may be challenged and removed. (April 2017) (Learn how and when to remove this template message) This article may be too technical for most readers to understand. Please to make it understandable to non-experts, without removing the technical details. (May 2011) (Learn how and when to remove this template message) (Learn how and when to remove this template message)

In physics, the Spalart–Allmaras model is a one-equation model that solves a modelled transport equation for the kinematic eddy turbulent viscosity. The Spalart–Allmaras model was designed specifically for aerospace applications involving wall-bounded flows and has been shown to give good results for boundary layers subjected to adverse pressure gradients. It is also gaining popularity in turbomachinery applications.

In its original form, the model is effectively a low-Reynolds number model, requiring the viscosity-affected region of the boundary layer to be properly resolved ( y+ ~1 meshes). The Spalart–Allmaras model was developed for aerodynamic flows. It is not calibrated for general industrial flows, and does produce relatively larger errors for some free shear flows, especially plane and round jet flows. In addition, it cannot be relied on to predict the decay of homogeneous, isotropic turbulence.

It solves a transport equation for a viscosity-like variable ${\displaystyle {\tilde {\nu }}}$. This may be referred to as the Spalart–Allmaras variable.

Original model[]

The turbulent eddy viscosity is given by

${\displaystyle \nu _{t}={\tilde {\nu }}f_{v1},\quad f_{v1}={\frac {\chi ^{3}}{\chi ^{3}+C_{v1}^{3}}},\quad \chi :={\frac {\tilde {\nu }}{\nu }}}$

${\displaystyle {\frac {\partial {\tilde {\nu }}}{\partial t}}+u_{j}{\frac {\partial {\tilde {\nu }}}{\partial x_{j}}}=C_{b1}[1-f_{t2}]{\tilde {S}}{\tilde {\nu }}+{\frac {1}{\sigma }}\{\nabla \cdot [(\nu +{\tilde {\nu }})\nabla {\tilde {\nu }}]+C_{b2}|\nabla {\tilde {\nu }}|^{2}\}-\left[C_{w1}f_{w}-{\frac {C_{b1}}{\kappa ^{2}}}f_{t2}\right]\left({\frac {\tilde {\nu }}{d}}\right)^{2}+f_{t1}\Delta U^{2}}$

${\displaystyle {\tilde {S}}\equiv S+{\frac {\tilde {\nu }}{\kappa ^{2}d^{2}}}f_{v2},\quad f_{v2}=1-{\frac {\chi }{1+\chi f_{v1}}}}$

${\displaystyle f_{w}=g\left[{\frac {1+C_{w3}^{6}}{g^{6}+C_{w3}^{6}}}\right]^{1/6},\quad g=r+C_{w2}(r^{6}-r),\quad r\equiv {\frac {\tilde {\nu }}{{\tilde {S}}\kappa ^{2}d^{2}}}}$

${\displaystyle f_{t1}=C_{t1}g_{t}\exp \left(-C_{t2}{\frac {\omega _{t}^{2}}{\Delta U^{2}}}[d^{2}+g_{t}^{2}d_{t}^{2}]\right)}$

${\displaystyle f_{t2}=C_{t3}\exp \left(-C_{t4}\chi ^{2}\right)}$

${\displaystyle S={\sqrt {2\Omega _{ij}\Omega _{ij}}}}$

The rotation tensor is given by

${\displaystyle \Omega _{ij}={\frac {1}{2}}(\partial u_{i}/\partial x_{j}-\partial u_{j}/\partial x_{i})}$

where d is the distance from the closest surface and ${\displaystyle \Delta U^{2}}$ is the norm of the difference between the velocity at the trip (usually zero) and that at the field point we are considering.

The constants are

${\displaystyle {\begin{matrix}\sigma &=&2/3\\C_{b1}&=&0.1355\\C_{b2}&=&0.622\\\kappa &=&0.41\\C_{w1}&=&C_{b1}/\kappa ^{2}+(1+C_{b2})/\sigma \\C_{w2}&=&0.3\\C_{w3}&=&2\\C_{v1}&=&7.1\\C_{t1}&=&1\\C_{t2}&=&2\\C_{t3}&=&1.1\\C_{t4}&=&2\end{matrix}}}$

Modifications to original model[]

According to Spalart it is safer to use the following values for the last two constants:

${\displaystyle {\begin{matrix}C_{t3}&=&1.2\\C_{t4}&=&0.5\end{matrix}}}$

Other models related to the S-A model:

DES (1999) [1]

DDES (2006)

Model for compressible flows[]

There are two approaches to adapting the model for compressible flows. In the first approach, the turbulent dynamic viscosity is computed from

${\displaystyle \mu _{t}=\rho {\tilde {\nu }}f_{v1}}$

where ${\displaystyle \rho }$ is the local density. The convective terms in the equation for ${\displaystyle {\tilde {\nu }}}$ are modified to

${\displaystyle {\frac {\partial {\tilde {\nu }}}{\partial t}}+{\frac {\partial }{\partial x_{j}}}({\tilde {\nu }}u_{j})={\mbox{RHS}}}$

where the right hand side (RHS) is the same as in the original model.

Boundary conditions[]

Walls: ${\displaystyle {\tilde {\nu }}=0}$

Freestream:

Ideally ${\displaystyle {\tilde {\nu }}=0}$, but some solvers can have problems with a zero value, in which case ${\displaystyle {\tilde {\nu }}\leq {\frac {\nu }{2}}}$ can be used.

This is if the trip term is used to "start up" the model. A convenient option is to set ${\displaystyle {\tilde {\nu }}=5{\nu }}$ in the freestream. The model then provides "Fully Turbulent" behavior, i.e., it becomes turbulent in any region that contains shear.

Outlet: convective outlet.

References[]

• Spalart, P. R. and Allmaras, S. R., 1992, "A One-Equation Turbulence Model for Aerodynamic Flows" AIAA Paper 92-0439

External links[]

• This article was based on the Spalart-Allmaras model article in CFD-Wiki
• What Are the Spalart-Allmaras Turbulence Models? from kxcad.net
• The Spalart-Allmaras Turbulence Model at NASA's Langley Research Center Turbulence Modelling Resource site