A subgrid-scale model represents the influence of unresolved motion on the quantities retained in a large eddy simulation. Its need follows directly from the nonlinear velocity product in the Navier-Stokes equations: averaging a product generally differs from multiplying the averages.
Let a bar denote a spatial averaging filter. With constant mass density , constant kinematic viscosity , and a time-independent translation-invariant filter on an unbounded or periodic domain, filtering the incompressible momentum equation giveswhere repeated indices use the Einstein summation convention, is pressure, and the subgrid-scale stress in this sign convention isThe assumed filter commutes with the displayed derivatives. The definition of is exact; a closure must approximate it using available resolved information.[ref]
A dynamic subgrid-scale eddy viscosity model by Massimo Germano and others (1990). Section2, equations3–7: filtering, filtered momentum and the unresolved stress
For an elementary example, average two equally weighted values of one velocity component, and . The mean is zero but the mean square is , so the corresponding stress component is . A resolved mean alone has lost the information needed to reconstruct that covariance.
The Smagorinsky model closes the deviatoric part with an eddy viscosity:where is the Kronecker delta, is the resolved strain-rate tensor, a dimensionless model coefficient and a characteristic filter width. The isotropic stress can be included in a modified pressure. This strain-dependent viscosity has its historical origin in Joseph Smagorinsky's atmospheric circulation modelling.[ref]
General Circulation Experiments with the Primitive Equations: I. The Basic Experiment by Joseph Smagorinsky. Abstract: lateral exchange coefficient based on horizontal strain and a length scale
[ref]A dynamic subgrid-scale eddy viscosity model by Massimo Germano and others (1990). Introduction, equations1–2: the standard Smagorinsky expression
For an incompressible flow, the model transfers resolved kinetic energy out at the rate per unit massA fixed nonnegative eddy viscosity therefore permits dissipation but excludes local turbulent backscatter. A dynamic subgrid-scale model can instead estimate its coefficient from resolved stresses at two filter levels. This addresses a modelling choice, rather than changing the exact definition of the unresolved stress.[ref]
A dynamic subgrid-scale eddy viscosity model by Massimo Germano and others (1990). Introduction and Section2: backscatter limitations and the dynamic two-filter construction
New to topics? Read the docs here!