Smagorinsky-Lilly Subgrid Scale Model

The Smagorinsky-Lilly SGS model is based on the Prandtl’s mixing length model and assumes that a kinematic SGS viscosity can be expressed in terms of the length scale and the strain rate magnitude of the resolved flow.

(1)

τ_{i j}^{'} = ρ \tilde{u_{i} u_{j}} - ρ \tilde{u_{i}} \tilde{u_{j}} = - 2 μ_{s} \tilde{S_{i j}}

where

$μ_{s} = ρ {(C_{S} Δ)}^{2} \sqrt{2 \tilde{S_{i j}} \tilde{S_{i j}}}$ is the subgrid turbulent viscosity,
Model constant 0.17 < $C_{s}$ < 0.21 (Lilly, 1996)
$\tilde{S_{i j}} = \frac{1}{2} (\frac{\partial \tilde{u_{i}}}{\partial x_{j}} + \frac{\partial \tilde{u_{j}}}{\partial x_{i}})$ is the filtered (resolved) strain rate tensor

A major drawback of this model is that the model constant ( $C_{s}$ ) does not vary in space and time. Furthermore, this model has no correct wall behavior and it is too dissipative for laminar turbulent transition cases. These limitations led to the development of a dynamic model for which the model constant is allowed to vary depending on the grid resolution and flow regime.