---
title: Smagorinsky/Ladyzhenskaya LES Model
url: https://www.emergentmind.com/topics/smagorinsky-ladyzhenskaya-type-les-model
type: topic
---

# Smagorinsky/Ladyzhenskaya LES Model

A Smagorinsky/Ladyzhenskaya-type LES model is a large-eddy-simulation closure in which unresolved turbulence enters through a nonlinear effective viscosity depending on the local magnitude of the resolved velocity gradient or strain. In its classical LES form, the model closes the deviatoric subgrid-scale stress by an eddy viscosity proportional to \((C_s\Delta)^2|S|\); in its Ladyzhenskaya formulation, it appears as a nonlinear diffusion of \(p\)-Laplacian type, typically written as \(-\nabla\cdot(\bar\nu |\nabla u|^{p-2}\nabla u)\), with the Smagorinsky case corresponding to \(p=3\) [2508.07492]. Across the literature represented here, this model family appears in engineering LES, finite-element regularizations, stochastic PDEs, data assimilation, reduced-order modeling, and several solver-aware reformulations [2510.15819, 2411.10630].

## 1. Classical constitutive structure

The standard incompressible LES starting point is the filtered Navier–Stokes system
\[
\frac{\partial \bar v_i}{\partial t} + \bar v_j \frac{\partial \bar v_i}{\partial x_j}
= - \frac{1}{\rho} \frac{\partial \bar p}{\partial x_i}
+ \nu_N \frac{\partial^2 \bar v_i}{\partial x_j \partial x_j}
+ \frac{\partial}{\partial x_j} M_{ij},
\qquad
\frac{\partial \bar v_i}{\partial x_i}=0,
\]
with subgrid stress
\[
M_{ij} := \bar v_i \bar v_j - \widetilde{v_i v_j}.
\]
The classical eddy-viscosity ansatz is
\[
M_{ij} = \nu\, S_{ij},
\qquad
S_{ij} := \frac12\left(\frac{\partial \bar v_i}{\partial x_j} + \frac{\partial \bar v_j}{\partial x_i}\right),
\]
and the Smagorinsky choice is
\[
\nu = C_s^2 \delta^2 (2S_{ij}S_{ij})^{1/2}.
\]
In the more conventional deviatoric-stress form, this is
\[
\tau_{ij}^{\rm SGS}-\frac13\tau_{kk}^{\rm SGS}\delta_{ij}=-2\nu_t S_{ij},
\qquad
\nu_t=(C_s\Delta)^2|S|.
\]
This is the canonical Smagorinsky closure used as the reference point in several of the cited works [1901.01654, 2411.10630].

In compressible LES, the same closure is written in Favre-filtered variables and usually combined with a modeled isotropic SGS contribution and a turbulent heat flux. For compressible homogeneous isotropic turbulence, the SGS stress is written as
\[
\tau_{ij}-\frac{1}{3}\delta_{ij}\tau_{kk}
=
-2\overline{\rho}\,\nu_{\mathrm{sgs}}
\left(
\widetilde{S}_{ij}-\frac13\delta_{ij}\widetilde{S}_{kk}
\right),
\]
with
\[
\nu_{\mathrm{sgs}}=(C_s\Delta)^2|\widetilde{S}|,
\qquad
\tau_{kk}=2C_I\overline{\rho}\Delta^2|\widetilde{S}|^2,
\qquad
Q_j=-\frac{1}{Pr_t}\overline{\rho}\,\nu_{\mathrm{sgs}}\frac{\partial \widetilde{T}}{\partial x_j},
\]
where the cited implementation used \(C_s=0.16\), \(C_I=0.09\), \(Pr_t=0.71\), and \(\Delta=(\Delta x\Delta y\Delta z)^{1/3}\) [2309.05875]. This compressible formulation preserves the same constitutive idea—strain-dependent eddy viscosity—while extending it to density-weighted filtering and filtered energy transport.

A finite-element realization used in reduced-basis LES replaces \(|S|\) by the Frobenius norm of the full gradient and defines a mesh-local viscosity,
\[
\nu_T(u) = C_S^2\sum_{K\in\mathcal{T}_h}h_K^2\big|\nabla u_{|_K}\big|\chi_K,
\]
which is still a Smagorinsky-type nonlinear viscosity, but intrinsically tied to the discrete mesh [2305.04598]. This suggests that, even within the classical family, the exact tensor norm and the relation between filter width and discretization are solver-dependent modeling choices rather than invariant features.

## 2. Ladyzhenskaya interpretation and generalized nonlinear viscosity

A central theme of the later mathematical literature is that Smagorinsky closure can be recast as a Ladyzhenskaya-type nonlinear constitutive law. In the data-assimilation formulation,
\[
\mathbf T_{\text{LES}}(\nabla u)
=
\left(\nu + \bar\nu |\nabla u|_F^{p-2}\right)\nabla u,
\]
so the LES model becomes
\[
\partial_t v + (v\cdot \nabla)v - \nu \Delta v
-\nabla\cdot\!\left(\bar\nu |\nabla v|_F^{p-2}\nabla v\right)
+ \nabla q = f,
\qquad
\nabla\cdot v=0.
\]
In this setting, the model is explicitly called “Smagorinsky/Ladyzhenskaya-type” because \(p=3\) yields the classical Smagorinsky nonlinear eddy viscosity [2508.07492]. The model is simultaneously a turbulence closure and a PDE regularization.

An even more general finite-element Ladyzhenskaya model writes the added nonlinear diffusion as
\[
-\nabla\cdot\left((C_S\delta)^r \|\nabla u\|_F^s \nabla u\right).
\]
The corresponding weak form contains
\[
(C_S\delta)^r(\|\nabla u\|_F^s \nabla u,\nabla v),
\]
and the paper states explicitly that \(r=2\) and \(s=1\) recover the classical Smagorinsky model [2510.15819]. In this sense, Smagorinsky appears as one distinguished member of a wider Ladyzhenskaya family parameterized by nonlinear diffusion exponents.

A related hybrid construction starts not from filtered Navier–Stokes, but from generalized Navier–Stokes with constitutive viscosity
\[
\mu_e=\mu_0\left(1+\lambda^2 |\mathbf D|^2\right)^q,
\qquad
\mathbf D=\frac12(\nabla\mathbf u+\nabla\mathbf u^T).
\]
For \(q=1/2\), the authors state that one recovers the Smagorinsky model by properly choosing \(\lambda\). After filtering, the remaining unclosed terms are approximated by Clark expansion rather than by a Boussinesq SGS ansatz [1509.04489]. This establishes a precise distinction between two traditions: in one, nonlinear viscosity is introduced after filtering as an SGS closure; in the other, it is already part of the constitutive PDE and filtering generates additional closure terms.

The relation is therefore structural rather than merely terminological. Smagorinsky closure is an LES eddy-viscosity model; Ladyzhenskaya regularization is a nonlinear-viscosity PDE model. The two coincide when the LES closure is written as a strain- or gradient-dependent viscosity and inserted directly into the momentum diffusion operator [2508.07492, 2510.15819].

## 3. Coefficients, dynamic procedures, and local variants

A persistent issue in the literature is the status of the Smagorinsky coefficient \(C_s\). In plane Couette LES with wall damping, the eddy viscosity is written as
\[
\nu_e=\{C_S\Delta(x_2)f_S(x_2)\}^2|S|,
\qquad
f_S(x_2)=1-\exp(-x_2^+/A^+),\quad A^+=25,
\]
with
\[
\Delta(x_2)=\{\Delta_1\Delta_2(x_2)\Delta_3\}^{1/3}.
\]
That study notes that the “typical value seems to be in the range of \(0.05\) to \(0.1\)” and uses \(C_S\) both as an LES parameter and as a bifurcation parameter for unstable periodic orbits [2107.14407]. In this wall-bounded context, the model is classical static Smagorinsky with van Driest damping.

Dynamic Smagorinsky models replace fixed \(C_s\) by coefficients determined from resolved fields at grid and test-filter levels. In the thermally driven cavity study, the dynamic model is described as allowing the Smagorinsky constant “to vary in space and time depending on an algebraic identity between the subgrid-scale stresses at two different filtered levels and the resolved turbulent stresses.” On the same coarse mesh, the dynamic model captured \(93\%\) of the reference peak vertical velocity fluctuation, whereas the standard model predicted about \(67\%\) near the wall; the dynamic simulation cost about \(33\%\) more CPU time [2110.00389]. In compressible homogeneous isotropic turbulence, dynamic Smagorinsky and Vreman were reported to be closer to DNS than static Smagorinsky for enstrophy and dilatation trends, while all models reproduced the overall kinetic-energy decay [2309.05875].

The local applicability of the dynamic procedure is not straightforward. A detailed analysis of the local DSM showed that, without averaging in homogeneous directions, the dynamic coefficient has a singular leading-order form and can develop exceedingly large positive values. The proposed remedy is the dynamic gradient Smagorinsky model (DGSM), which replaces the strain tensor by the full resolved velocity gradient tensor. In its basic form,
\[
\tau_{ij}-\frac13\tau_{kk}\delta_{ij}
=
-2 C_{\text{GS}}\overline{\Delta}^2
|\nabla \overline{\mathbf u}|\,
(\partial_j \overline u_i),
\]
with a corresponding local Germano coefficient \(C_{\text{GS}}\). The paper shows that the leading-order singularity of local DSM is removed, and the DGSM remained stable in decaying isotropic turbulence, temporal mixing layers, and turbulent channel flow without spatial or temporal averaging, while local DSM became unstable in the reported tests [2110.05585]. This does not eliminate the broader problem of coefficient sensitivity, but it shifts the model from a singular local dynamic law to a bounded one.

These developments show that “the Smagorinsky constant” is not a single question. It is alternately treated as a fixed empirical constant, a dynamic field, a bifurcation or homotopy parameter, or a learned constitutive function. The family remains unified by the nonlinear-viscosity structure, not by a universal prescription for \(C_s\).

## 4. Anisotropy, wall treatment, and discretization-aware reformulations

A major reformulation is the Lattice Eddy Simulation (LAES) approach. LAES starts from a node-interaction picture on a CFD mesh and derives extra terms
\[
T_u=(C_s\Delta x)^2 S \frac{\partial^2 u}{\partial x^2}
+(C_s\Delta y)^2 S \frac{\partial^2 u}{\partial y^2}
+(C_s\Delta z)^2 S \frac{\partial^2 u}{\partial z^2},
\]
with analogous expressions for \(T_v\) and \(T_w\). On isotropic meshes, \(\Delta x=\Delta y=\Delta z=\Delta\), these reduce to
\[
T_u=(C_s\Delta)^2 S \nabla^2 u,
\qquad
\nu_t=(C_s\Delta)^2 S,
\]
which the authors state is identical to classical Smagorinsky LES [2211.12810]. On anisotropic meshes, however, LAES is not standard Smagorinsky LES: it retains directional widths explicitly and yields what the paper describes as an “asymmetric subgrid stress matrix” if interpreted in classical LES terms.

This anisotropic form is used to argue against the usual need for wall damping or dynamic adjustment. LAES sets \(C_s=0.08\) in all simulations and attributes improved wall-normal fluctuations to the wall-normal filter width being set to \(\Delta y\). Since \(\Delta y\) is small near the wall on channel meshes, wall-normal turbulent diffusion is automatically reduced. The paper therefore claims that LAES “need not use any ad hoc damping functions or any dynamic filtering procedure for wall turbulence,” but the same text also states that the universality claim is only empirically suggested by moderate-\(Re_\tau\) channel-flow evidence and is not broadly established [2211.12810].

A different reformulation arises from time-discretization-aware LES. For forward Euler, the time difference satisfies the exact filter-swap identity
\[
\partial_t^\tau = \overline{(\cdot)}^\tau \circ \partial_t,
\]
which yields an exact space-time residual decomposition with a new temporal term
\[
r^\tau_{\mathrm{time}}(u)=\overline{f(u)}^\tau-f(u).
\]
Its leading-order form is
\[
r^\tau_{\mathrm{time}}(u)
=
\frac{\tau}{2}\partial_t f(u)+O(\tau^2)
=
-\frac{\tau}{2}(f'(u))^2\partial_x u + O(\tau^2),
\]
a Lax–Wendroff-type diffusion. The resulting augmented Smagorinsky closure is
\[
m^{h,\tau}(\cdot)
=
-\, \theta_h^2 h^2\, |\partial_x^h \cdot|\, \partial_x^h \cdot
-\, \theta_\tau\, \tau\, (f'(\cdot))^2\, \partial_x^h \cdot .
\]
In Burgers tests, this time-aware augmentation remained accurate at coarse time steps where space-only Smagorinsky closures degraded [2606.17759]. This suggests that “Smagorinsky-type” dissipation can be made discretization-aware in both space and time, not only in the filter-width sense.

## 5. Mathematical analysis, stochastic formulations, and data assimilation

The Ladyzhenskaya interpretation makes the model accessible to monotone-operator and SPDE analysis. For stochastic Ladyzhenskaya–Smagorinsky equations with damping,
\[
du(t)+\big[\mu A(u(t))+B(u(t))+\beta C(u(t))\big]dt
=
f(t)\,dt+\sum_{k=1}^\infty \sigma_k(t,u(t))\,dW_k(t),
\]
with
\[
A(u)=-P\,\mathrm{div}\Big((1+|\nabla u|^2)^{\frac{p-2}{2}}\nabla u\Big),
\qquad
C(u)=P(|u|^{r-2}u),
\]
the paper proves local monotonicity for \(p>\frac d2+1\), \(r\ge2\), and global monotonicity for \(p>2\), \(r\ge4\), followed by existence and pathwise uniqueness of strong solutions and a small-time large deviation principle [2106.10861]. The Smagorinsky case is \(p=3\). This places the model in a rigorous SPDE framework rather than only an LES-engineering one.

For the 3D periodic stochastic Ladyzhenskaya–Smagorinsky equations,
\[
du +  \left( u \cdot \nabla u -  \nu \Delta u - \nabla \cdot \left( \bar{\nu} \,  |\nabla u|^{r-2} \nabla u \right) + \nabla  p \right) dt
=
f\,dt + g(t,u)\,dW,
\]
the long-time averaged dissipation rate
\[
\varepsilon=\varepsilon_0+\varepsilon_M
\]
is shown to satisfy an upper bound of Kolmogorov order \(U^3/L\), and the paper emphasizes that it remains finite as \(\nu\to0\). The conclusion drawn there is that, in a periodic domain without boundary layers, the model “does not over-dissipate” [2510.12004]. This sharply qualifies a common criticism of Smagorinsky closure: the cited result does not deny over-dissipation in wall-bounded flows, but it does separate that phenomenon from the periodic bulk case.

In continuous data assimilation, the forecast model
\[
\partial_t v + (v\cdot \nabla)v - \nu \Delta v
-\nabla\cdot\!\left(\bar\nu |\nabla v|_F^{p-2}\nabla v\right)
+ \nabla q
=
f + \mu I_h(u)-\mu I_h(v)
\]
is nudged toward NSE observations \(u\). In 2D, the assimilated system is proved globally well posed for \(p\ge\frac52\) under \(c_0\mu h^2\le \nu\), and the error satisfies
\[
\|u-v\|^2 \le C_{u,\Omega}\,\bar\nu(1-e^{-t}) + \|u(0)-v(0)\|\,e^{-t},
\]
so that
\[
\|u-v\|_{L^2} = \mathcal O(\bar\nu^{1/2})
\quad\text{as }t\to\infty.
\]
The paper interprets the residual floor as the irreducible model-observation mismatch caused by the extra LES dissipation [2508.07492]. In this formulation, Smagorinsky/Ladyzhenskaya-type LES is not only a turbulence model but also a biased forecast model whose structural error can be quantified.

At the fully discrete level, the EMAC-Ladyzhenskaya formulation combines nonlinear viscosity with the EMAC convection form,
\[
c(u,v,w)=2(D(u)v,w)+((\nabla\cdot u)v,w),
\]
which conserves energy, linear momentum, and angular momentum under weak incompressibility enforcement. The resulting EMAC-LM scheme retains the same asymptotic order as skew-symmetric discretizations but has a Gronwall factor without explicit Reynolds-number dependence, and the reported step-flow benchmarks showed better long-time preservation of momentum and angular momentum than the skew form [2510.15819].

## 6. Applications, limitations, and current directions

The cited applications span incompressible channel flow, plane Couette flow, atmospheric boundary layers, buoyancy-driven cavity flow, compressible homogeneous isotropic turbulence, supersonic jets, reduced-order models, SPH reformulations, and machine-learned non-Boussinesq closures. This breadth shows that Smagorinsky/Ladyzhenskaya-type modeling functions less as a single model than as a constitutive backbone reused across numerics and flow classes [2211.12810, 2309.05875].

At the same time, several limitations recur. Static Smagorinsky is repeatedly described as over-dissipative, especially near walls or in mean-strain-dominated regions [2107.14407, 2212.12353]. In coarse LES of a heated cavity at \(Ra=10^9\), the standard model predicted first moments reasonably well but underpredicted near-wall fluctuation levels relative to the dynamic model [2110.00389]. In compressible homogeneous isotropic turbulence, all tested SGS models captured the overall decay trends, but the same assessment reported underprediction of peak enstrophy and dilatation and overprediction of peak temperature variance, underscoring the difficulty of compressible turbulence dynamics for classical eddy-viscosity closures [2309.05875].

Another recurring theme is that discretization can rival or dominate SGS modeling. In low-order LES of supersonic jets, the authors concluded that the characteristics of numerical discretization can be as important as the effects of the SGS models, and that static Smagorinsky, dynamic Smagorinsky, and Vreman produced similar behavior on the refined grid [2212.12353]. This suggests that model assessment cannot be separated from numerical scheme order, built-in dissipation, and filter-width definition.

Several recent directions explicitly move beyond the strict Boussinesq assumption. In APG turbulent boundary layers, a numerically consistent machine-learned SGS stress
\[
\boldsymbol{\tau}_{\text{SGS}}
=
C_1(\mathbf{\Pi};\boldsymbol{\theta})\, \Delta^2 \lVert\mathbf{S}\rVert \mathbf{S}
+
C_2(\mathbf{\Pi};\boldsymbol{\theta})\, \Delta^2(\mathbf{S}\mathbf{R}-\mathbf{R}\mathbf{S})
\]
was proposed to overcome “the limitations of linear eddy-viscosity closures in complex flows.” The reported a posteriori tests improved mean velocity and wall-shear stress relative to the Dynamic Smagorinsky model and achieved monotonic convergence with grid refinement [2601.20265]. This does not discard Smagorinsky entirely—the first tensor basis term is still Smagorinsky-like—but it treats linear stress-strain proportionality as insufficient in non-equilibrium wall turbulence.

Other extensions change the constitutive law rather than the tensor basis. A PDE-constrained optimization framework on the 1D Kuramoto–Sivashinsky equation replaces the fixed Smagorinsky law
\[
\nu(s)=C_s^2\delta^2|s|
\]
by an optimized constitutive function \(\nu(|s|)\). The identified optimal law is nonlinear, can become negative for small strain, and reduced trajectory divergence relative to the standard Smagorinsky model, while also illustrating the inherent limitations of the eddy-viscosity paradigm [1901.01654]. This suggests that the best model inside the Smagorinsky/Ladyzhenskaya class need not be a positive linear function of \(|S|\).

Finally, some reformulations shift the modeling target altogether. In SPH, a rigorous filtering derivation shows that properly smoothed SPH solves the same filtered equations as explicit LES and could, in principle, use Smagorinsky or dynamic Smagorinsky closures, although the paper argues that approximate deconvolution is more natural in that context [1807.11244]. In reduced-order modeling, Smagorinsky viscosity has been embedded into reduced basis formulations and monitored through a Kolmogorov-spectrum indicator rather than only through residual norms [2305.04598].

Taken together, these works portray the Smagorinsky/Ladyzhenskaya-type LES model as a durable but non-final closure class. Its enduring feature is the nonlinear viscosity law; its unresolved questions concern coefficient universality, tensor structure, wall asymptotics, stochastic consistency, solver dependence, and the extent to which eddy viscosity alone can represent backscatter, anisotropy, and non-equilibrium transfer.

Source: https://www.emergentmind.com/topics/smagorinsky-ladyzhenskaya-type-les-model