2.2.2. Large-Eddy Simulation Approach

In the LES approach, the filtered flow quantities are solved by the transport equations on the CFD mesh. For example,   now refers to the filtered flow velocity. The undetermined stress term   in Equation 2–3 is called the sub-grid-scale (SGS) stress tensor, which is essential in LES modeling. Ansys Forte offers two types of LES models, the Smagorinsky model and the dynamic structure model.

The Smagorinsky model [89] is based on a viscosity assumption and accounts for the dissipative nature of turbulent flows, that is, kinetic energy is dissipated from the large-scale to the small-scale motions. It relates the SGS stress to the strain rate of the filtered flow field:

(2–30)

The SGS viscosity, , is modeled as:

(2–31)

in which is the local CFD mesh size, is the magnitude of the filtered strain rate tensor (see Equation 2–25 , but note that the overbar   now means filtering), and is a model constant, 0.17. The sub-grid kinetic energy, , is modeled as:

(2–32)

where is a model constant, 0.101. It is noted that both and depend on the CFD mesh size, .

The sub-grid flux term   in the species transport Equation 2–1 is modeled as:

(2–33)

in which is the SGS turbulent diffusivity. Similarly, the sub-grid flux term H in the energy equation (Equation 2–5 ) is modeled as:

(2–34)

in which is the SGS thermal conductivity and is related to the SGS thermal diffusivity and heat capacity by . The SGS mass and thermal diffusivity are related to the SGS viscosity by:

(2–35)

(2–36)

where and are the turbulent Schmidt and Prandtl numbers, respectively. The viscosity-based assumption makes Smagorinsky model effective in dissipating the kinetic energy from large-scale to small-scale motions. However, such dissipation can be excessive, and energy-containing flow structures may not be adequately resolved unless a very fine mesh is used. Although, its numerical stability in complex flows is an advantage.

The dynamic structure model [71] is a non-viscosity and similarity-based model, relating the SGS stress to the “Leonard stress” derived from a larger filter:

(2–37)

where L is the Leonard stress tensor, , and . An explicit filtering operation (denoted as  ) needs to be performed to compute the Leonard stress tensor, and the explicit filter size is chosen as twice of the local CFD mesh size. Essentially, the first term on the right-hand-side of Equation 2–37 can be viewed as a similarity model relating the SGS stress to the Leonard stress, with a dynamic scaling factor of representing the ratio of sub-grid kinetic energy to sub-filter kinetic energy. The second term is an optional added viscosity-based model [96] . It is suggested to be activated for fuel injection and spray simulations for numerical stability purposes. The viscosity is modeled as

(2–38)

where and . The symbols   and   are injection velocity and nozzle diameter in sprays.

To use the dynamic structure model, the sub-grid flux terms and can still adopt the viscosity-based models shown in Equation 2–33 and Equation 2–34 . The SGS viscosity is calculated as:

(2–39)

where . A transport equation for the sub-grid kinetic energy, , needs to be solved to provide closure

(2–40)

where . The dissipation rate of the sub-grid kinetic energy has been modeled as .

The dynamic structure model offers good prediction of the sub-grid stresses [81] , and has been validated in a number of benchmark flow problems [96] , [98] . Note that in view of our best practice recommendations in engine flow simulations, the Smagorinsky model is recommended as the default option due to its superior numerical stability in dealing with flows near compelx boundaries.