4.14.1. Shielding Function

The starting point for the development of the Shielded Detached Eddy Simulation (SDES) shielding function is the existing Delayed Detached Eddy Simulation (DDES) models (see Detached Eddy Simulation (DES)). In DDES, the primary goal is to model the attached and mildly separated boundary layers in RANS mode and then switch to LES mode for separated (detached) shear layers. Starting from a given RANS model, a DDES formulation is achieved in two-equation models by a modification of the sink term in the -equation. In the -equation of the SST model, the modification is reformulated as follows:

(4–288)

where

(4–289)

In the previous equation, is the turbulence length scale, is the maximum edge length of the local cell, is the DDES shielding function and is a coefficient. When equals zero, the RANS model is recovered. This happens when one or both of the following conditions are satisfied:

(4–290)

and/or

(4–291)

The function is designed to protect attached boundary layers from the switch to the LES regime. The limit at which the boundary layer becomes affected depends on the ratio of the grid spacing to the boundary layer thickness (for example, =/). Without shielding (that is, when =0) the limit is approximately =1, so that at <1 the solution in the boundary layer becomes affected by the DES sink term. With the shielding function as given in Detached Eddy Simulation (DES), the limit for is reduced to 0.2-0.3, depending on the pressure gradient. Violating these limits has severe consequences, as the RANS solution will become compromised, and you might not be aware of it. DES / DDES methods were originally developed for aircraft simulations, where the boundary layer is typically thin relative to the dimensions of the aircraft / wing. Under those conditions (and with diligent mesh generation practices), the previously described limits for can typically be ensured. However, it is harder to ensure these limits for general industrial applications, due to lower Reynolds numbers, less control over the grid spacing, and more complex interactions between components.

The first goal of the development of the SDES method was therefore the reformulation of the function to achieve more reliable shielding, meaning the ability to maintain RANS boundary layers even under much more severe mesh refinement in attached boundary layers than is possible for DDES. For this purpose, a new unpublished and Ansys proprietary shielding function called was developed, which provides substantially safer shielding than the functions currently used. The resulting sink term is equivalent to the one from the DDES model:

(4–292)

where

(4–293)

Figure 4.7: Eddy Viscosity Profiles shows the eddy viscosities for self-similar zero-pressure-gradient boundary layer simulations under mesh refinement, where =/ (note the different range of the refinement ratio for the two models). The SDES shielding (on the right) is clearly superior to the DDES shielding, and maintains the RANS boundary layer even under severe mesh refinement.

Figure 4.7: Eddy Viscosity Profiles

Eddy Viscosity Profiles