--- title: Smagorinsky Model in LES Turbulence url: https://www.emergentmind.com/topics/smagorinsky-model type: topic --- # Smagorinsky Model in LES Turbulence The Smagorinsky model is a large-eddy simulation (LES) subgrid-scale closure in which unresolved stresses are represented by an algebraic eddy viscosity proportional to a filter length squared and to the magnitude of the resolved deformation rate. In filtered notation, its classical form is \[ \tau^{sgs}_{ij}-\tfrac13\tau^{sgs}_{kk}\delta_{ij}=-2\,\nu_t\,\overline S_{ij}, \qquad \nu_t=(C_S\Delta)^2\,|\overline S|, \] while, in one periodic-box formulation, the same closure appears through the nonlinear diffusion term \[ -\nabla\!\cdot\Bigl((C_s\delta)^2\,|\nabla u|\,\nabla u\Bigr) \] in the momentum equation [2110.05585][1601.06635]. The model is one of the most studied sub-grid scale representations in LES, and its simplicity has made it a reference point for analyses of dissipation, wall effects, dynamic coefficient procedures, reduced-order approximations, and structurally enriched SGS closures [2409.06047]. ## 1. Classical closure and parameterization In LES, filtering the incompressible Navier–Stokes equations introduces the subgrid-scale stress \[ \tau_{ij}=\overline{u_i u_j}-\overline u_i\,\overline u_j, \] which must be modeled. The classical Smagorinsky closure adopts the Boussinesq hypothesis and assumes that the deviatoric SGS stress is aligned with the resolved rate-of-strain tensor, \[ \overline S_{ij}=\tfrac12\bigl(\partial_i\overline u_j+\partial_j\overline u_i\bigr), \qquad |\overline S|=\sqrt{2\,\overline S_{ij}\,\overline S_{ij}}, \] through the eddy viscosity \[ \nu_t=(C_S\Delta)^2\,|\overline S| . \] Equivalent formulations appear in incompressible, compressible, and periodic-box settings, sometimes written with \(|\nabla u|\) rather than \(|S|\), but the constitutive idea is the same: unresolved motions are modeled by an effective viscosity that grows with local resolved deformation [2308.07978][1601.06635]. The filter width \(\Delta\) is typically tied to the discretization. Examples in the literature include \(\Delta=2\pi/N\) on an \(N^3\) pseudo-spectral grid, \(\Delta^3=\Delta_1\Delta_2\Delta_3\) on a rectangular grid, and \(\Delta=\sqrt{\Delta x^2+\Delta y^2+\Delta z^2}\) in HIGRAD/FIRETEC [2308.07978][2110.05585][2409.06047]. Representative coefficient choices reported in the literature include \(C_s\approx0.1\), \(C_s=0.16\) for homogeneous isotropic turbulence, \(C_S\approx0.17\) in a canopy LES implementation, and values in the range \(\approx0.1\text{–}0.2\) in many flows [1601.06635][2308.07978][2409.06047][2411.05640]. A defining feature of the classical model is that it is purely dissipative. The modeled SGS dissipation in one standard formulation is \[ \varepsilon_m=C_S^2\,\Delta^2\,|\overline S|^3, \] so increasing \(C_S\) increases dissipation and damps small scales more strongly [2110.05585]. This also means that the classical eddy-viscosity model cannot represent backscatter, i.e. energy transfer from unresolved fluctuations back into the resolved mean or resolved scales [2201.08914][2308.07978]. ## 2. Energy balance and dissipation scaling For the periodic Smagorinsky model on \(\Omega=(0,L_\Omega)^3\), \[ u_t+u\cdot\nabla u-\nu\Delta u-\nabla\!\cdot\Bigl((C_s\delta)^2|\nabla u|\nabla u\Bigr)+\nabla p=f, \qquad \nabla\cdot u=0, \] the exact energy equality for strong solutions is \[ \frac{1}{2|\Omega|}\frac{d}{dt}\|u(t)\|_{L^2}^2 +\frac{\nu}{|\Omega|}\|\nabla u\|_{L^2}^2 +\frac{(C_s\delta)^2}{|\Omega|}\|\nabla u\|_{L^3}^3 = \frac{1}{|\Omega|}\int_\Omega f\cdot u\,dx . \] The instantaneous total model dissipation is therefore \[ \varepsilon_S(u)=\frac{\nu}{|\Omega|}\|\nabla u\|_{L^2}^2 +\frac{(C_s\delta)^2}{|\Omega|}\|\nabla u\|_{L^3}^3, \] with long-time average \[ \langle\varepsilon_S\rangle = \limsup_{T\to\infty}\frac1T\int_0^T\varepsilon_S(u(t))\,dt \] [1601.06635]. Using the large-scale quantities \[ F=\Bigl(\tfrac1{|\Omega|}\|f\|_{L^2}^2\Bigr)^{1/2},\qquad U=\Bigl\langle \tfrac1{|\Omega|}\|u\|_{L^2}^2\Bigr\rangle^{1/2}, \qquad L=\min\Bigl\{|\Omega|^{1/3},\,F/\|\nabla f\|_{L^\infty},\,F/\|\nabla f\|_{L^2},\,F/\|\nabla f\|_{L^3}\Bigr\}, \] with \(Re=UL/\nu\), the periodic-box analysis yields the bound \[ \langle\varepsilon_S\rangle \le 3\,\frac{U^3}{L} +\frac{9}{8}\,\frac{\nu U^3}{L^3} +3\,C_s^2\Bigl(\frac{\delta}{L}\Bigr)^2\frac{U^3}{L}, \] or equivalently \[ \langle\varepsilon_S\rangle \le 3\,\frac{U^3}{L} \Bigl(1+\frac{3}{8\,Re}+C_s^2\Bigl(\frac{\delta}{L}\Bigr)^2\Bigr). \] As \(\delta\to0\) or \(C_s\to0\), the last term vanishes; as \(Re\to\infty\), the \(Re^{-1}\) term vanishes [1601.06635]. This result is central in the dissipation theory of the model. For body-force-driven turbulence under periodic boundary conditions, it gives the same \(U^3/L\) scaling as the incompressible Navier–Stokes equations up to a vanishing \(Re^{-1}\) correction and a vanishing \((\delta/L)^2\) correction. The stated conclusion is that, without boundary layers, the unmodified Smagorinsky model does not over-dissipate [1601.06635]. ## 3. Boundary layers, mesh effects, and the over-dissipation problem The Smagorinsky model is often reported to over-damp or over-dissipate in wall-bounded and shear-driven flows. In lid-driven shear flow on an under-resolved mesh \(h\) satisfying \[ O(Re^{-1})0\) and \(\delta>0\), consistent with numerical evidence on the effects of model viscosity without wall damping [1804.00745]. That analysis also states that, for the lid-driven problem, the turbulent boundary layer is a more important length scale than the Kolmogorov microscale. Shear flows display two small scales, \[ \eta\sim Re^{-3/4}L \qquad\text{and}\qquad \delta_{bl}\sim Re^{-1}L, \] and the controlling scale for energy dissipation in the lid-driven problem is \(\delta_{bl}\sim Re^{-1}L\), not \(\eta\) [1804.00745]. This is consistent with the periodic-box conclusion that the dominant over-dissipation mechanism is associated with boundary effects rather than with the turbulent cascade itself [1601.06635]. Several remedies have therefore been developed. One is damping of the SGS viscosity near walls. With a damping function \(\beta(x)\), the model becomes \[ u_t+u\cdot\nabla u-\nu\Delta u+\nabla p-\nabla\cdot\bigl[\beta(x)(C_s\delta)^2|\nabla u|\nabla u\bigr]=0, \qquad \nabla\cdot u=0, \] and the time-averaged dissipation is \[ \langle\varepsilon_{SMD}(u)\rangle = \limsup_{T\to\infty}\frac1T\int_0^T \frac1{|\Omega|}\int_\Omega \nu|\nabla u|^2+(C_s\delta)^2\beta(x)|\nabla u|^3\,dx\,dt . \] For a modified \(C^1\) damping function \(\beta_d\) with algebraic wall contact order \(\alpha\ge2\), the reported bound is \[ \langle\varepsilon_{SMD}(u)\rangle \le \bigl[C_1+C_2(C_s\delta/L)^2\bigr]\frac{U^3}{L}, \] which restores \(U^3/L\) scaling independently of \(R\) as \(R\to\infty\) [1608.02655]. Another remedy is spatially adaptive damping. One recent adaptive formulation replaces the constant coefficient by \[ C_S(x)=C_{S0}\exp[-\beta\,\Phi(x)], \qquad \Phi(x)=\text{distance from }x\text{ to }\partial\Omega, \] so that \[ \nu_t(x)=[C_S(x)\Delta]^2\,|S(\bar u)|. \] That work states existence and uniqueness of weak solutions and an energy dissipation bound for the adaptive model, with the explicit aim of reducing over-dissipation near boundaries [2411.05640]. ## 4. Dynamic, corrected, and structurally enriched variants The classical constant-coefficient model has generated a broad family of variants designed to reduce excessive dissipation, recover backscatter-like effects, or relax the Boussinesq alignment assumption. | Variant | Defining modification | Stated purpose or effect | |---|---|---| | Dynamic Smagorinsky | \(C_s^2=\langle L_{ij}M_{ij}\rangle/\langle M_{ij}M_{ij}\rangle\) | Compute the coefficient on-the-fly from two filter levels | | Dynamic gradient Smagorinsky | Replace \(\overline S_{ij}\) by \(G_{ij}=\partial_j\overline u_i\) | Remove the local dynamic singularity without averaging | | Corrected Smagorinsky model | Add \(-\frac{C_s^4\delta^2}{\mu^2}\Delta w_t\) | Incorporate energy flow from unresolved fluctuations back into the mean | | Dynamic tensor-coefficient model | Replace scalar coefficient by tensor \(C_{ij}\) | Relax local stress–strain alignment | | Helicity SGS model | Add helicity-gradient/vorticity coupling to Smagorinsky viscosity | Incorporate structural information beyond turbulence intensity | The dynamic Smagorinsky procedure of Germano et al. compares SGS stresses at a grid filter \(\overline{(\cdot)}\) of width \(\Delta\) and a test filter of width \(\widehat\Delta=2\Delta\). In one formulation, \[ C_s^2=\frac{\langle L_{ij}M_{ij}\rangle}{\langle M_{ij}M_{ij}\rangle}, \] and in rotating stratified turbulence the coefficient was stabilized by whole-domain averaging in the absence of walls [2308.07978]. A related issue is that the widely used dynamic Smagorinsky model becomes singular if used locally without averaging. The local dynamic gradient Smagorinsky model addresses this by replacing the rate-of-strain tensor by the full resolved velocity gradient, \[ G_{ij}=\partial_j\overline u_i, \qquad \tau_{ij}-\tfrac13\tau_{kk}\delta_{ij} = -2\,C_{\rm GS}\,\Delta^2\,|G|\,G_{ij}, \] which removes the singularity because the denominator now decays only when the full \(G_{ij}\) vanishes [2110.05585]. The corrected Smagorinsky model is motivated by the fact that the classical model is an eddy-viscosity approximation to a resolved mean velocity and cannot represent a flow of energy from unresolved fluctuations to the resolved mean velocity. Its governing equation is \[ w_t-\frac{C_s^4\delta^2}{\mu^2}\Delta w_t +w\cdot\nabla w -\nu\Delta w +\nabla q -\nabla\!\cdot\!\bigl((C_s\delta)^2|\nabla w|\,\nabla w\bigr) =f, \qquad \nabla\cdot w=0, \] and the added Voigt-type term is designed to model the omitted \(\frac{d}{dt}\|u'\|^2\) contribution from the Reynolds-stress energy balance [2201.08914]. The dynamic tensor-coefficient Smagorinsky model replaces the scalar eddy-viscosity coefficient by a tensorial coefficient. In one reduced formulation, \[ \tau_{ij}^{sgs}-\tfrac13\delta_{ij}\tau^{sgs}_{kk} = -\bigl(C_{ik}S_{kj}+C_{jk}S_{ki}\bigr)\,|S|\,\Delta^2, \] with four independent dynamic coefficients after imposing tracelessness constraints. That model is reported to vanish in laminar flow and at solid boundaries, and to have the correct asymptotic behavior in the near-wall region of a turbulent boundary layer [2202.05502]. A different structural enrichment incorporates helicity. In LES notation, the deviatoric SGS stress is modeled as \[ \tau^{SGS}_{ij}\big|_D = -2\,\nu_S\,\overline S_{ij} + \eta_S\,[\partial_iH_S\,\overline\omega_{*j} +\partial_jH_S\,\overline\omega_{*i} -\tfrac23\delta_{ij}(\overline\omega_*\cdot\nabla H_S)], \] with \[ \nu_S=(C_s\Delta)^2|\overline S|, \qquad \eta_S=C_h\,\Delta^2/|\overline S|. \] The stated motivation is that the standard Smagorinsky framework uses only information about turbulence intensity through turbulent energy, whereas helicity supplies structural or geometrical information of turbulence [2603.26559]. ## 5. Functional analysis, well-posedness, and regularity The Smagorinsky model has also become a substantial object of nonlinear analysis. For the unsteady rotational Smagorinsky model, the classical viscous term is replaced by a weighted curl operator, \[ \partial_t v+\omega\times v+\curl\!\bigl(C_\alpha d(x)^\alpha |\omega|\,\omega\bigr)+\nabla\bar q=f, \qquad \omega=\curl v, \qquad \Div v=0, \] with \(d(x)=\mathrm{dist}(x,\partial\Omega)\). In the natural case \(p=3\), the appropriate space is \[ V=W^{1,3}_{0,\sigma}(\Omega,d^\alpha), \] and existence of weak solutions follows from Bochner pseudo-monotone evolution equations. The threshold \(\alpha<2\) is identified as intrinsic, and \(p=3\) is the minimal exponent making the wall-degenerate convective term manageable in the pseudo-monotone setting [2107.00236]. Recent variational refinements formulate the SGS stress as \[ \tau_{SGS}(u)=-2\,\nu_t(u)\,S(u), \qquad \nu_t(u)=(C_s\Delta)^2|S(u)|+\epsilon\,\mathcal F(u), \] where \(\epsilon>0\) is a small regularization parameter and \(\mathcal F(u)\) ensures \(\nu_t\ge\nu_0>0\). In that framework, existence of weak solutions is obtained on bounded Lipschitz domains for \(u_0\in H\) and \(f\in L^2(0,T;H^{-1}(\Omega))\), while uniqueness is asserted under the additional regularity \(u\in L^4(0,T;W^{1,4}(\Omega)^n)\) [2411.10630]. Analyses with dynamic or evolving boundaries extend the same program. One formulation couples the bulk Smagorinsky system \[ \partial_tu+(u\cdot\nabla)u-\nu\Delta u+\nabla p-\nabla\cdot\tau(u)=f, \qquad \tau(u)=2C_s^2\delta^2|S(u)|S(u), \] to a Wentzell-type tangential boundary evolution and proves existence, uniqueness, an energy estimate, and higher-order Sobolev regularity [2411.06230]. Related corrected-model analyses in time-dependent domains derive \(H^s\)-bounds, \(L^2\)-error estimates relative to Navier–Stokes, and long-time convergence statements under assumptions on forcing and domain smoothness [2412.16230]. A further high-order extension introduces both a mild inverse-Laplacian temporal correction and an \(s\)-order dissipation term, \[ \partial_t w - C_s\,\delta_s\,\mu\,(-\Delta)^{-2}\partial_t w +(w\cdot\nabla)w -\nu\Delta w +\nabla q -\nabla\cdot\bigl[(C\delta)^s|\nabla w|\nabla w\bigr] =f, \] and proves enhanced regularity and error estimates in higher Sobolev spaces [2411.05249]. A common theme across these analyses is that the nonlinear SGS operator is treated through monotonicity, coercivity, compactness, or pseudo-monotonicity. This has turned the Smagorinsky closure from a purely empirical LES device into a model with a substantial functional-analytic theory. ## 6. Discretization, reduced models, and representative applications The model remains prominent in applied LES because its algebraic form is easy to embed in existing solvers. In HIGRAD/FIRETEC, the Smagorinsky closure \[ \tau^r_{ij}=-2\,\rho_g\,\nu_r\,\overline S_{ij}, \qquad \nu_r=(C_S\Delta)^2|\overline S|, \] was implemented with fixed \(C_S=0.17\) and \(\Delta\) equal to the local cell diagonal. In neutral, no-fire simulations of flow over a homogeneous forest edge, the modeled fraction \(\langle K_S\rangle_t\) was \(<5\%\) of the total Smagorinsky TKE profile, peaks in modeled TKE occurred at shear interfaces, and resolved TKE profiles, mean velocity profiles, sweep/ejection ratios, and spectral slopes were comparable to those of FIRETEC’s original Linn 1.5-order TKE closure. Against field data at mid-canopy, the reported root-mean-square errors were \(0.35\) for Linn and \(0.57\) for Smagorinsky [2409.06047]. In decaying rotating stratified turbulence, LES with constant-coefficient Smagorinsky, dynamic Smagorinsky, and dynamic Clark closures were compared against filtered DNS on a \(192^3\) grid. The reported result was that all SGS model predictions were very similar, with the classical Smagorinsky model displaying the highest deviation compared to DNS. Dynamic Smagorinsky improved the fit up to \(\kappa\approx10\text{–}20\), but both Smagorinsky variants over-dissipated the smallest LES scales and failed to reproduce inverse-cascade behavior because \(C_s^2>0\) everywhere, so no energy flows from subgrid to resolved scales even when the DNS does [2308.07978]. The model has also become a target for advanced discretizations. A certified reduced-basis Smagorinsky model for steady flows used the Empirical Interpolation Method to approximate the nonlinear eddy diffusion term and reported speed-up factors \(\simeq170\text{–}1\,460\times\) for a backward-facing step and \(\simeq1\,300\text{–}3\,200\times\) for a lid-driven cavity in FreeFem++ [1709.00243]. In the virtual element setting, a general-order discretization with the local eddy viscosity \[ \nu_S(\mathbf u)(\mathbf x) = C_S^2\sum_{T\in\mathcal T_h}h_T^2\,|\nabla\mathbf u|_T|_{\ell^2}\,\chi_T(\mathbf x), \qquad C_S=0.1, \] was tested up to \(Re=10\,000\). The reported conclusion was that isotropic refinement with hanging nodes enhanced accuracy relative to anisotropic meshes, used fewer degrees of freedom than uniform meshes, and yielded the most stable Newton convergence behavior [2510.03563]. Across these settings, a consistent pattern emerges. The classical Smagorinsky model remains attractive because it is local, algebraic, and computationally light; it captures many mean and resolved statistics adequately in homogeneous or wall-free settings; but it also inherits structural limitations tied to positivity of eddy viscosity, local stress–strain alignment, and sensitivity to boundary layers. Much of the subsequent literature can be read as a sequence of attempts to preserve its simplicity while correcting those three liabilities through damping, dynamic coefficient evaluation, modified dissipation, tensorial coefficients, or additional structural information such as helicity [1601.06635][2110.05585][2202.05502][2603.26559]. Source: https://www.emergentmind.com/topics/smagorinsky-model