The Scale-Adaptive Simulation (SAS) is an improved URANS formulation, which allows the resolution of the turbulent spectrum in unstable flow conditions. The SAS concept is based on the introduction of the von Karman length-scale into the turbulence scale equation. The information provided by the von Karman length-scale allows SAS models to dynamically adjust to resolved structures in a URANS simulation, which results in a LES-like behavior in unsteady regions of the flowfield. At the same time, the model provides standard RANS capabilities in stable flow regions.
In recent years the turbulent length-scale equation has been revisited by Menter and Egorov [130]. It was shown that the exact transport equation for the turbulent length-scale, as derived by Rotta [134], does actually introduce the second derivative of the velocity field and thereby into the turbulent scale equation. In Menter and Egorov [130], a two-equation turbulence model was presented using a formulation, which can be operated in RANS and SAS mode. While the further development of this model is still ongoing, it was considered desirable to investigate if the SAS term in the model could be transformed to existing two-equation models. The target two-equation model is the SST model and this leads to the formulation of the SST-SAS model.
The original version of the SST-SAS model (Menter and Egorov [131]) has undergone certain evolution and the latest model version has been presented in Egorov and Menter [197]. One model change is the use of the quadratic length scale ratio in the Equation 2–210 below, rather than the linear form of the original model version. The use of the quadratic length scale ratio is more consistent with the derivation of the model and no major differences to the original model version are expected. Another new model aspect is the explicitly calibrated high wave number damping to satisfy the requirement for an SAS model that a proper damping of the resolved turbulence at the high wave number end of the spectrum (resolution limit of the mesh) must be provided. In the following the latest model version of the SST-SAS model (Egorov and Menter [197]) will be discussed, which is also the default version in Ansys CFX.
The governing equations of the SST-SAS model differ from those of the SST RANS model [129] by the additional SAS source term in the transport equation for the turbulence eddy frequency :
(2–208) |
(2–209) |
where is the value for the regime of the SST model.
The additional source term reads for the latest model version Egorov and Menter [197]:
(2–210) |
This SAS source term originates from a term
in Rotta’s transport equation for the correlation-based length scale, see Menter and Egorov [130]. Because the integral is zero in homogeneous turbulence, it should in general be proportional to a measure related to inhomogeneity. The second velocity derivative was selected as this measure to ensure that the integral alone is modeled to zero at a constant shear rate, thus leading to and ultimately to in the SAS source term (Equation 2–210).
This model version (Egorov and Menter [197], Model
Version=2007
) is used as default. In order to recover the
original model formulation (Menter and Egorov [131]), the parameter Model Version
must be set to 2005
directly in the CCL:
FLUID MODELS: ... TURBULENCE MODEL: Model Version = 2005 Option = SAS SST END ... END
The model parameters in the SAS source term Equation 2–210 are
Here is the length scale of the modeled turbulence
(2–211) |
and the von Karman length scale given by
(2–212) |
is a three-dimensional generalization of the classic boundary layer definition
The first velocity derivative is represented in by , which is a scalar invariant of the strain rate tensor :
(2–213) |
Note, that the same also directly participates in (Equation 2–210) and in the turbulence production term .
The second velocity derivative is generalized to 3D using the magnitude of the velocity Laplacian:
(2–214) |
As a result, and are both equal to () in the logarithmic part of the boundary layer, where =0.41 is the von Karman constant.
Beside the use of the quadratic length scale ratio , the latest model version provides for the direct control of the high wave number damping. Two formulations are available:
The first formulation is the default and is realized by a lower constraint on the value in the following way:
(2–215)
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. Assuming the source term equilibrium (balance between production and destruction of the kinetic energy of turbulence) in both transport equations, one can derive the following relation between the equilibrium eddy viscosity , and :
(2–216)
This formula has a similar structure as the subgrid scale eddy viscosity in the LES model by Smagorinsky:
Therefore it is natural to adopt the Smagorinsky LES model as a reference, when formulating the high wave number damping limiter for the SST-SAS model. The limiter, imposed on the value, must prevent the SAS eddy viscosity from decreasing below the LES subgrid-scale eddy viscosity:
(2–217)
Substitution of and in condition above (Equation 2–217) results in the limiter value used in Equation 2–215. Similar to LES, the high wave number damping is a cumulative effect of the numerical dissipation and the SGS eddy viscosity. The model parameter has been calibrated using decaying isotropic turbulence. The default value of is 0.11, which provides nearly the same energy spectrum as the Smagorinsky LES model.
The second formulation limits the eddy viscosity directly:
(2–218)
The LES-WALE model is used for the calculation of the LES eddy viscosity because this model is suitable for transitional flows and does not need wall damping functions. This limiter can be turned on by setting the expert parameter
limit sas eddy viscosity = t
. The limiter (Equation 2–215 is then automatically turned off. The default value for the WALE model coefficient is 0.5 and can be specified by the expert parameterlimit sas eddy viscosity coef=0.5
.
Similar to the DES formulation, the SAS model also benefits from a switch in the numerical treatment between the steady and the unsteady regions. In DES, this is achieved by a blending function as proposed by Strelets [58], which allows the use of a second order upwind scheme with the CFX-Solver in RANS regions and a second order central scheme in unsteady regions. The blending functions are based on several parameters, including the mesh spacing and the ratio of vorticity and strain rate, as explained below.