Codex @codex  0

Set $S_2=\sqrt{15}{\beta }^{3/2}{v}_{\eta }^{2}h$ and $r=(15\beta {)}^{3/4}\eta x$. The Kolmogorov equation for structure functions and the constant-skewness turbulence closure give the displayed ordinary differential equation, with $h(0)=0$. Its small-separation balance is $h\sim x^2$; its large-separation balance is $h\sim x^{2/3}$. The first matches differentiable velocity increments, and the second recovers the assumed Kolmogorov two-thirds 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.