Dynamic Subgrid Scale Model

Recognizing C s variations in space and time, Germano et al. (1991) proposed the dynamic model to compute the value of C s rather than specifying it explicitly.

It is implemented by utilizing two filters: a cutoff filter Δ and a test (coarse) cutoff filter Δ .

The subgrid stress tensor τ i j ' with the cutoff filter ( Δ ) is: τ i j ' = ρ u i u j ˜ ρ u i ˜ u j ˜ = 2 ρ ( C S Δ ) 2 | S ˜ | S i j ˜ . Where | S ˜ | = 2 S i j ˜ S i j ˜ is the strain rate magnitude.

Figure 1 shows a resolved turbulence region utilizing Large Eddy Simulation (LES) and a modeled region assuming the subgrid tensor τ i j ' .
Figure 1. Energy Spectrum for LES Using the Cutoff Filter Width ( Δ )


The test subgrid stress tensor T i j with the coarse filter ( Δ ) can be written as

T i j = ρ u i u j ˜ ˜ ρ u i ˜ ˜ u j ˜ ˜ = 2 ρ ( C S Δ ) 2 | S ˜ ˜ | S i j ˜ ˜

where
  • | S ˜ ˜ | = 2 S i j ˜ ˜ S i j ˜ ˜ is the coarse filtered strain rate magnitude.
  • S i j ˜ ^ = 1 2 ( u i ˜ ^ x j + u j ˜ ^ x i ) is the filtered strain rate tensor, using the coarse cutoff filter.
Figure 2 shows a resolved LES region and a corresponding subgrid modelled region (T) when the coarse filter is employed.
Figure 2. Energy Spectrum for LES Using the Test Cutoff Filter Width ( Δ )


Because of the coarse filtering, the test (coarse) subgrid stress tensor T i j should be a summation of the coarse filtered subgrid stress tensor τ i j ' ˜ and the Leonard stress tensor L i j .

L i j = T i j τ i j ' ˜ =   ρ u i u j ˜ ˜ ρ u i ˜ ˜ u j ˜ ˜ ρ u i u j ˜ ˜ + ρ u i ˜ u j ˜ ˜ = ρ u i ˜ u j ˜ ˜ ρ u i ˜ ˜ u j ˜ ˜

where
  • τ i j ' ˜ is the subgrid tensor for the cutoff filter (or grid filtered), then test filtered.

τ i j ' ˜ = ρ u i u j ˜ ˜ ρ u i ˜ u j ˜ ˜ = 2 ρ ( C S Δ ) 2 | S ˜ | S i j ˜ ˜

  • L i j is the Leonard subgrid stress tensor, representing the contribution to the subgrid stresses by turbulence length scales smaller than the test filter but larger than the cutoff filter.

The Leonard subgrid stress tensor can be arranged as

L i j = 2 ρ ( C S Δ ^ ) 2 | S ˜ ˜ | S i j ˜ ˜ + 2 ρ ( C S Δ ) 2 | S ˜ | S i j ˜ ˜ = 2 ρ C S 2 Δ 2 ( | S ˜ | S i j ˜ ˜ α 2 | S ˜ ˜ | S i j ˜ ˜ )

where α = Δ ^ / Δ .

The Leonard subgrid stress tensor can be rewritten as

L i j = 2 ρ C S 2 Δ 2 ( | S ˜ | S i j ˜ ˜ α 2 | S ˜ ˜ | S i j ˜ ˜ ) = C S 2 M i j

where M i j = 2 ρ Δ 2 ( | S ˜ | S i j ˜ ˜ α 2 | S ˜ ˜ | S i j ˜ ˜ ) .

Since the above equation is overdetermined a minimum least square error method is used to determine the coefficient C s .

C s 2 = L i j M i j M i j M i j

In order to avoid numerical instabilities associated with the above equation, as the numerator could become negative, averaging of the error in the minimization is employed.
C s 2 = L i j M i j M i j M i j