4.12.2. Transport Equations for the SST-SAS Model

From a fundamental standpoint, all the currently used two-equation models suffer from lack of an underlying exact transport equation, which can serve as a guide for the model development on a term by term basis. The reason for this deficiency lies in the fact that the exact equation for does not describe the large scales, but the dissipative scales. The goal of a two-equation model is however the modeling of the influence of the large scale motions on the mean flow. Due to the lack of an exact equation, the - and the -equations are modeled in analogy with the equation for the turbulent kinetic energy, , using purely heuristic arguments. A more consistent approach for formulating a scale-equation has been developed by Rotta (1968, 1972). Instead of using purely heuristic and dimensional arguments, Rotta formulated an exact transport equation for turbulent kinetic energy times length scale, . Rotta’s equation represents the large scales of turbulence and can therefore serve as a basis for term-by-term modeling. The transport equations for the SST-SAS model as implemented in Ansys Fluent are based on transforming Rotta’s approach to - (SST) and defined as:

(4–263)

(4–264)

For the detailed derivation see Egorov and Menter  [431]. The transport equations of the SST-SAS model Equation 4–263 and Equation 4–264 differ from those of the SST RANS model  [428] by the additional SAS source term in the transport equation for the turbulence eddy frequency , Equation 4–264. In Equation 4–264 is the value for the - regime of the SST model.

The additional source term reads (see Egorov and Menter  [431] for details):

(4–265)

This SAS source term originates from a second order derivative term in Rotta’s transport equations  [430]. The model parameters in the SAS source term Equation 4–265 are

Here is the length scale of the modeled turbulence

(4–266)

and the von Kármán length scale is a three-dimensional generalization of the classic boundary layer definition

(4–267)

The first velocity derivative is represented in by , which is equal to , a scalar invariant of the strain rate tensor :

(4–268)

Note, that the same also directly participates in (Equation 4–265) and in the turbulence production term . The second velocity derivative is generalized to 3D using the magnitude of the velocity Laplacian:

(4–269)

So defined and are both equal to in the logarithmic part of the boundary layer, where is the von Kármán constant.

The model also provides a direct control of the high wave number damping. This is realized by a lower constraint on the value in the following way:

(4–270)

This limiter is proportional to the mesh cell size , which is calculated as the cubic root of the control volume size . The purpose of this limiter is to control damping of the finest resolved turbulent fluctuations. The structure of the limiter is derived from analyzing the equilibrium eddy viscosity of the SST-SAS model  [431].