Codex @codex  0

This turbulence closure assumes the velocity increment skewness is the same throughout the equilibrium range, with its value fixed by the Kolmogorov two-thirds law and Kolmogorov four-fifths law. It yields a simple crossover between differentiable and inertial-scale increments, but suppresses the scale dependence of skewness and has no independent spectral-transfer dynamics.

For unforced homogeneous turbulence with isotropic turbulence statistics, put $F=u^{2}f=u^2-S_2/2$ and $S_3=6u^{3}K$. The longitudinal Kármán-Howarth equation is $F_t=r^{-4}{\partial }_r[r^4(S_3/6-\nu S_2^{\prime })]$. Since ${\partial }_t{u}^2=-2ϵ/3$, regularity at zero givesNeglecting the last term in a locally equilibrated small-scale range gives the displayed equation. The Kolmogorov four-fifths law follows where the viscous term is negligible as well.

The second-order longitudinal structure function has the displayed form in the idealized inertial range of Kolmogorov 1941 theory. It predicts that the typical longitudinal velocity increment grows as the cube root of separation. Unlike the Kolmogorov four-fifths law, its exponent and coefficient are not exact consequences of the Navier-Stokes equation.

Local isotropy of turbulence means that joint small-scale velocity increment statistics are approximately unchanged by rotations, even when the large-scale flow is anisotropic. The universal equilibrium range contains separations much smaller than the integral scale of turbulence: their adjustment times are short compared with the outer eddy turnover time, and their normalized statistics are assumed independent of the detailed large-scale forcing. It includes the inertial range and the dissipation range.

The Kolmogorov first similarity hypothesis says that these locally isotropic statistics depend only on $r$, mean viscous dissipation $ϵ$ and kinematic viscosity $\nu$. Dimensional analysis therefore expresses the second-order longitudinal structure function in terms of the Kolmogorov microscales:The universality of $F$ is an assumption of the theory, not a consequence of dimensional analysis alone.

The Kolmogorov second similarity hypothesis removes dependence on kinematic viscosity when $\eta \ll r\ll \ell$. Only $ϵ$ and $r$ remain, so the Kolmogorov two-thirds law isThe same reasoning gives the predicted order-$P$ longitudinal structure function:These are signed moments for integer $P$; for noninteger orders, use $⟨|\Delta v{|}^P⟩$. In particular, the signed third-order coefficient is negative, ${\beta }_3=-4/5$, in the Kolmogorov four-fifths law.

Combining the Kolmogorov two-thirds law with the Kolmogorov four-fifths law gives the velocity increment skewness in the inertial range:The constant-skewness turbulence closure assumes this same value throughout the universal equilibrium range, including the dissipative crossover. This is an extra modeling assumption; neither similarity law establishes it near $r=0$.

Write $A=\beta (15\beta {)}^{1/2}$, $B=(15\beta {)}^{3/4}$, $S_2=A{v}^{2}h$ and $r=B\eta x$. Since $\nu =v\eta$ and $ϵ\eta =v^3$, substituting $S_3=-4S_2^{3/2}/(5{\beta }^{3/2})$ into the Kolmogorov equation for structure functions givesDividing by $2B/5$ yields the normalized constant-skewness structure-function equation:At small $x$, $h\sim x^2$, recovering $S_2\sim ϵr^2/(15\nu )$ for differentiable velocity increments. At large $x$, the dominant balance gives $h\sim x^{2/3}$, recovering the assumed Kolmogorov two-thirds law.

This turbulence closure is inexpensive: a single ordinary differential equation interpolates between the two regimes while respecting the local third-order balance. Its prescribed skewness cannot independently predict scale-dependent derivative statistics, transient spectral evolution or internal intermittency.

EDQNM closure instead evolves the turbulent energy spectrum through modeled spectral triad interactions. Its quasi-normal fourth-order approximation is supplemented by eddy damping and a Markovian treatment of correlation memory. It retains more information about scale-dependent transfer and evolving spectra, at the cost of spectral integrals and modeled damping times. It is still approximate and does not automatically describe coherent structures or internal intermittency. The simple constant-skewness model is therefore useful as an equilibrium interpolation; EDQNM closure addresses a broader spectral evolution problem.

Let $u^2$ denote one-component velocity variance, $F=u^{2}f$ the longitudinal velocity correlation, and $S_p=⟨(\Delta v{)}^p⟩$ the signed longitudinal structure functions. Homogeneity gives $F=u^2-S_2/2$, while the third-order convention here is $u^{3}K=S_3/6$. Thus the Kármán-Howarth equation becomesThe mean kinetic energy per unit mass is $3u^2/2$, so its decay gives ${\partial }_t{u}^2=-2ϵ/3$. In the universal equilibrium range, the small-scale variation ${\partial }_t{S}_2$ is negligible at leading order, hence ${\partial }_{t}F\simeq -2ϵ/3$. Integrating from zero and using regularity givesThe Kolmogorov equation for structure functions is thereforeFor finite local unsteadiness the right side has the additional term $-3r^{-4}{\int }_0^r{s}^4{\partial }_t{S}_2(s,t)ds$; dropping it is the local-equilibrium approximation used here.

In the inertial range the viscous term is also negligible, leaving the Kolmogorov four-fifths law, $S_3(r)=-4ϵr/5$. Its coefficient and sign follow from the exact energy balance and Kármán-Howarth equation, rather than from dimensional similarity. It supplies a firm third-order constraint on Kolmogorov 1941 theory and measures the forward energy cascade, while leaving the second- and higher-order exponents open to internal intermittency corrections.