Dynamic Smagorinsky Model
Dynamic Smagorinsky: Theoretical Foundations
Overview
Professor, what aspect of the Smagorinsky model did the dynamic Smagorinsky model improve?
The biggest problem with the Smagorinsky model is that it uses a fixed value for the model constant $C_s$, even though it depends on the flow. The dynamic Smagorinsky model (Germano et al., 1991; Lilly, 1992) is a method to determine $C_s$ locally and dynamically during the computation. It utilizes a mathematical identity called the Germano Identity.
Germano Identity
What is the Germano Identity?
It is an identity concerning two-stage filtering: the grid filter $\bar{\phantom{u}}$ (width $\Delta$) and the test filter $\hat{\phantom{u}}$ (width $\hat{\Delta} = 2\Delta$).
This $L_{ij}$ (Leonard stress tensor) can be calculated directly from known resolved-scale quantities. On the other hand, from the consistency condition with the model expression for the SGS stress,
Dynamic Calculation of $C_s^2$
How is $C_s^2$ determined?
Use Lilly's (1992) least squares method.
$\langle \cdot \rangle$ represents averaging in spatial directions (homogeneous directions) or Lagrangian averaging. Without this averaging, $C_s^2$ can become locally negative or oscillate violently, making the computation unstable.
What is the physical meaning of $C_s^2$ becoming negative?
A negative $C_s^2$ means backscatter, i.e., reverse energy transport from the SGS scale to the resolved scale. This is a physically possible phenomenon, but numerically it can cause instability. In implementation, the lower limit of $C_s^2$ is clipped to zero, or negative values are suppressed using Lagrangian averaging.
Advantages of the Dynamic Model
What are the specific advantages of the dynamic model?
| Advantage | Explanation |
|---|---|
| $C_s \to 0$ at walls | Automatically gives correct wall behavior without Van Driest damping |
| Reproduction of transitional flow | $C_s \to 0$ in laminar regions, avoiding excessive dissipation |
| Adaptation to different Re numbers | No need for prior tuning of constants |
| Partial reproduction of backscatter | Physically allows $C_s^2 < 0$ |
The Revolutionary Nature of Germano's Idea: "Let the Flow Itself Determine the Constant"
The revolutionary aspect of the dynamic procedure proposed by Massimo Germano in 1991 lies in the idea that "the SGS constant does not need to be determined by humans." By solving the Germano identity using a test filter and determining the locally optimal $C_s$ from the flow field itselfโthis idea fundamentally changed the conventional thinking of "calibrating model parameters through experiments." However, the initial paper had the problem of the constant becoming negative and causing divergence in the computation. The following year, Lilly (1992) stabilized it using averaging via the least squares method, giving it a practical form.
Computational Methods for Dynamic Smagorinsky
Test Filter Implementation
How exactly is the test filter implemented?
For unstructured grids, the average (volume-weighted) of cell center values from neighboring cells is often used as the test filter. For structured grids, top-hat filters or Gaussian filters can be used.
| Filter Type | Implementation Method | Accuracy |
|---|---|---|
| Top-hat (Box) | Simple average of neighboring cells | Low (problematic on non-uniform grids) |
| Volume-weighted average | $\hat{\phi}_P = \sum_f V_f \phi_f / \sum_f V_f$ | Medium |
| Gaussian | Gaussian weights based on distance | High (suited for structured grids) |
Stabilization Techniques
How do you deal with the problem of $C_s^2$ becoming negative and causing instability?
There are mainly three approaches.
1. Clipping: Clip to zero when $C_s^2 < 0$. Simplest but loses physical backscatter.
2. Spatial averaging: Average $\langle L_{ij}M_{ij}\rangle$ in homogeneous directions (e.g., spanwise direction). Suitable for channel flow.
3. Lagrangian averaging (Meneveau et al. 1996): Time average along fluid particle paths. Applicable to non-homogeneous flows.
The Lagrangian dynamic model is the most versatile and is implemented in OpenFOAM as dynLagrangian. It is also available in Fluent as the Dynamic Smagorinsky-Lilly model.
Solver Settings
Please tell me how to set it up in each solver.
The dynamic Smagorinsky is theoretically the most elegant method among LES SGS models, and its major strength is that it requires no constant tuning. However, the trade-off is implementation complexity and cost (test filter computation).
The Problem: Is "Twice" the Right Size for the Test Filter Width?
The test filter width $\hat{\Delta}$ in the dynamic procedure is conventionally set to $\hat{\Delta} = 2\Delta$, but many people are at a loss for a clear answer when asked why it's twice. Theoretically, it's "the ratio that yields the maximum information within the range where inertial subrange scale similarity holds," but in reality, it's a rather empirical choice. Studies have tried $\hat{\Delta} = 3\Delta$ or $4\Delta$, and cases of grid dependence have been reported. Instead of "just setting it to twice," we recommend conducting sensitivity tests.
Dynamic Smagorinsky in Practice
Applicability
In what situations should the dynamic Smagorinsky model be used?
| Suitable Applications | Reason |
|---|---|
| Transitional flow (flow including laminarโturbulent transition) | $C_s$ automatically becomes zero in laminar regions |
| Industrial LES with complex geometry | No tuning required, versatile |
| LES resolving near-wall regions | Correct wall behavior without Van Driest damping |
| Systematic studies at different Re numbers | Constants adjust automatically |
How much does the computational cost increase compared to the standard Smagorinsky model?
Related Topics
Experience the theory firsthand with the interactive simulator for this field
All Simulators