Implications for Large-Eddy Simulation · Arun, Kamal, Colonius & Johnson · J. Fluid Mech.

In 3D turbulence, the energy cascade exhibits a spectral bump exceeding Kolmogorov's scaling near the viscous cutoff — the bottleneck effect.
Both share a common structural origin: small-scale shear layers.
Stokes Flow Regularization generates dynamic LES closures. A mixed model (eddy viscosity + nonlinear gradient) mitigates the artificial bottleneck by capturing local residual stress structure.
A normality-based decomposition of filtered velocity gradients reveals that shear layers drive strain self-amplification, vortex stretching, and the backscatter causing the bottleneck.
Rather than the conventional strain + vorticity split, the VGT is decomposed in a principal reference frame into three physically distinct constituents:
Symmetric & normal — axis-aligned stretching and compression
Non-normal — captures shear layers, the key small-scale structure
Antisymmetric & normal — captures vortex tubes
Using Gaussian filtering, interscale energy transfer is expanded in multiscale velocity gradients and decomposed into three mechanisms:

~60% of forward cascade; normal straining dominates
~40% of forward cascade; shear vorticity dominates
Negligible net in inertial range, but drives backscatter in subinertial range

The normality-based decomposition unambiguously distinguishes tube-like structures (rigid rotation) from sheet-like structures (shear vorticity) — a distinction obscured by symmetry-based analysis.
Associated with rigid rotation. Filtering preserves tube-like character. Dominant in unfiltered DNS at low .
Associated with shear vorticity. Widths of . Softened by filtering; dominant in small-scale energy transfer.

The mixed model reproduces filtered DNS partitioning across all resolved scales. The eddy viscosity model deviates near the filter width — its statistics collapse onto a DNS at (effective ).

The eddy viscosity model produces excessive backscatter from the strain–vorticity covariance term — specifically from shear layer contributions — in the inertial range.
This artificial bottleneck has the same shear-layer origin as the true DNS bottleneck, but appears at the wrong scales because the model retains unfiltered-like shear layer structure.
Adding a nonlinear gradient component (mixed model) explicitly resolves scale-local terms, eliminating the spurious backscatter and restoring correct cascade efficiency.
Straining of multiscale shear layers accounts for >50% of energy transfer; shear layer stretching drives scale-local vortex stretching.
The bottleneck effect — true and artificial — is almost entirely due to normal straining interacting with shear vorticity at small scales.
Explicitly capturing scale-local terms via a nonlinear gradient component reproduces filtered DNS statistics and eliminates the artificial bottleneck.
Normality-based analysis provides a principled, quantitative tool for evaluating and designing LES closures across flow regimes.
Multiscale Velocity Gradients & Energy Transfer