Codex @codex  0

EDQNM closure models the evolution of the turbulent energy spectrum through spectral triad interactions. Fourth-order correlations are approximated using products of second-order correlations, an eddy damping limits the resulting triple correlations, and their memory integrals are replaced by expressions using current spectra. The modeled damping times are additional physical input. This allows evolving spectra and scale-dependent transfer, but does not automatically recover coherent structures or internal intermittency. Bertoglio, Squires and Ferziger's simulation comparison explains the damping assumptions and tests their limitations.

Goldreich–Sridhar turbulence combines critical balance, balanced counterpropagating Alfvén waves, a local perpendicular cascade, and constant energy flux per unit mass $\mathcal{E}$. With ${\tau }_{\mathrm{nl}}\sim {\ell }_{\perp }/\delta u_{\ell }$ and $\mathcal{E}\sim \delta u_{\ell }^2/{\tau }_{\mathrm{nl}}$, one obtains $\delta u_{\ell }\sim (\mathcal{E}{\ell }_{\perp }{)}^{1/3}$. Defining the one-dimensional perpendicular turbulent energy spectrum by $\delta u_{\ell }^2\sim k_{\perp }E(k_{\perp })$ gives $E(k_{\perp })\sim {\mathcal{E}}^{2/3}{k}_{\perp }^{-5/3}$. Critical balance also gives $k_{\parallel }\sim {\mathcal{E}}^{1/3}{k}_{\perp }^{2/3}/v_A$. This scaling argument neglects scale-dependent alignment and intermittency corrections, which are additional physical effects rather than consequences of the stated assumptions.

The Iroshnikov-Kraichnan spectrum phenomenology assumes balanced, isotropic local Alfvénic turbulence with $\delta z_l\ll v_A$, decorrelated counterpropagating encounters, ${\tau }_A=l/v_A$ and ${\tau }_{\mathrm{nl}}=l/\delta z_l$. A weak wave cascade takes ${\tau }_{\mathrm{cas}}\sim {\tau }_{\mathrm{nl}}^2/{\tau }_A=v_{A}l/\delta z_l^2$. Constant energy flux gives $\delta z_l^4\sim ϵv_{A}l$ and the displayed shell-integrated turbulent energy spectrum. Isotropy is an additional assumption, not a consequence of a strong guide field.

The Loitsyansky integral is the second spatial moment of the trace velocity correlation tensor, with the displayed sign convention. If correlations decay sufficiently rapidly, it is the coefficient of $k^4/(24{\pi }^2)$ in the turbulent energy spectrum. In terms of the normalized longitudinal velocity correlation it is $8\pi u^2{\int }_0^{\infty }{r}^{4}f(r)dr$ when the boundary terms vanish; conventions omitting $8\pi$ also occur. Its evolution is $\dot{I}=8\pi [r^4{u}^{3}K{]}_{\infty }-12\nu L$, where $L$ is the Saffman integral. Thus it is invariant when $L=0$ and the longitudinal triple correlation is $o(r^{-4})$, not under arbitrary turbulent conditions.

Use incompressible ideal MHD at scales above the ion skin depth, a strong uniform mean field, and a statistically steady local energy cascade with mean flux $ϵ$ per unit mass. Assume balanced counterpropagating wave populations, comparable velocity and magnetic amplitudes in velocity units, and no additional scale-dependent alignment or intermittency correction. Write $\mathbf{b}=\delta \mathbf{B}/\sqrt{4\pi \rho }$ and introduce the fluctuation Elsässer variables ${\mathbf{z}}^{\pm }=\mathbf{u}\pm \mathbf{b}$. Their equations areOnly opposite populations distort one another. A single population has no leading nonlinear self-interaction, so the balance assumption is essential to a one-amplitude cascade estimate.

Let $\delta z_l$ denote the typical balanced amplitude and let $l_{\perp },l_{\parallel }$ be coherence lengths relative to the local mean field. The Alfvén wave time and eddy turnover time areFor $\chi \ll 1$, one counterpropagating encounter changes the packet fractionally by order $\chi$. Assume phases decorrelate between encounters. After $N$ encounters, random accumulation is of order $\sqrt{N}\chi$; an order-one change requires $N\sim {\chi }^{-2}$. The weak wave cascade time is thereforeWeakness by itself does not prove this random walk estimate; phase decorrelation and the interaction geometry are additional assumptions.

For the Iroshnikov-Kraichnan spectrum, impose isotropy, $l_{\parallel }\sim l_{\perp }\sim l$, in addition to $\delta z_l\ll v_A$. Then ${\tau }_A=l/v_A\ll {\tau }_{\mathrm{nl}}=l/\delta z_l$ andHence $\delta z_l\sim (ϵv_{A}l{)}^{1/4}$. A local shell-integrated turbulent energy spectrum has $\delta z_l^2\sim kE(k)$, $k\sim l^{-1}$, soThe isotropy is imposed phenomenologically; a strong guide field actually distinguishes directions and does not force this isotropic picture. The spectrum is not an automatically valid asymptotic result for every strongly magnetized cascade.

For the anisotropic weak Alfvénic cascade, keep $l_{\parallel }$ approximately fixed while cascading to smaller $l_{\perp }$, and require $\chi \ll 1$. This geometry follows from resonant three-wave Alfvén interactions: frequency conservation together with parallel-wavenumber conservation makes one member of an oppositely propagating resonant triad have zero parallel wavenumber, so the propagating modes retain their parallel wavenumber in the leading resonant transfer. Thus there is no leading parallel cascade in this ideal weak-wave ordering. Using the fixed parallel scale,Therefore $\delta z_l\sim (ϵv_A/l_{\parallel }{)}^{1/4}{l}_{\perp }^{1/2}$. For the one-dimensional perpendicular spectrum, $\delta z_l^2\sim k_{\perp }{E}_{\perp }(k_{\perp })$ givesUnlike the isotropic estimate, the parallel length does not decrease with the perpendicular length. Since $\chi =\delta z_l{l}_{\parallel }/(v_A{l}_{\perp })\propto l_{\perp }^{-1/2}$, the interactions become stronger down the cascade. If ${\chi }_0\ll 1$ at outer perpendicular scale $L_{\perp }$, the weak approximation breaks at $l_{\perp ,\ast }\sim L_{\perp }{\chi }_0^2$. This is the weak-to-strong transition of an Alfvénic cascade.

For Goldreich–Sridhar turbulence, abandon the weak-collision estimate and impose critical balance, ${\tau }_A\sim {\tau }_{\mathrm{nl}}$. The local cascade is strong, with transfer time of order one eddy turnover time. ThenThe parallel scale is set dynamically, not held fixed:At smaller scales $l_{\parallel }/l_{\perp }$ grows, so the eddies are increasingly elongated along the local magnetic field. Strong turbulence does not mean $\delta z_l\sim v_A$: small transverse amplitudes can still be strong when the eddies are sufficiently anisotropic to satisfy critical balance. Imbalance, alignment or nonlocal interactions would require different cascade reasoning.

Finally, in the usual isotropic hydrodynamic Kolmogorov 1941 theory inertial-range argument, the only relevant dimensional quantities are $ϵ$ and $k$. With $[ϵ]=L^2{T}^{-3}$ and $[E]=L^3{T}^{-2}$, writing $E\propto {ϵ}^a{k}^b$ forces $a=2/3$, $b=-5/3$. This uses the physical assumptions that outer-scale, viscous and intermittency effects do not supply another leading dependence; dimensions alone are not a proof of those assumptions.

A magnetic guide field supplies the independent speed $v_A$, and anisotropy adds independent parallel lengths. Even before adding anisotropy, dimensions permitwith arbitrary dimensionless $\Phi$. For example $\Phi (x)\propto x^{1/2}$ reproduces the isotropic Iroshnikov-Kraichnan spectrum estimate. This dimensional freedom of an MHD spectrum explains why propagation, decorrelation, anisotropy and cascade-time assumptions must select among the spectra: purely dimensional considerations do not select a unique MHD exponent.

Assume statistically steady, three-dimensional turbulence at large Reynolds number, with local statistical homogeneity and statistical isotropy at scales small compared with the forcing scale $L$. The Kolmogorov 1941 theory additionally assumes a local energy cascade, constant mean energy flux per unit mass $ϵ$, and negligible viscous dissipation within the inertial range. Its dimensional form neglects corrections from internal intermittency.

At separation $l$, let $\delta u_l$ be a typical velocity increment. The eddy turnover time is ${\tau }_l\sim l/\delta u_l$. Constant transfer rate then gives $ϵ\sim \delta u_l^2/{\tau }_l\sim \delta u_l^3/l$, henceHere $E(k)$ is the one-dimensional turbulent energy spectrum, obtained from $\delta u_l^2\sim kE(k)$ with $k\sim l^{-1}$. The scaling holds only forThe Kolmogorov length scale follows by setting the scale-dependent Reynolds number $\delta u_{l}l/\nu$ to one. The associated Kolmogorov microscales give ${\tau }_{\eta }=(\nu /ϵ{)}^{1/2}$. At larger scales the forcing matters; at smaller scales kinematic viscosity matters.

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.

Use the Goldreich–Sridhar turbulence model in its balanced form. Assume a uniform strong guide field; anisotropic fluctuations $k_{\parallel }\ll k_{\perp }$; comparable counterpropagating Alfvén wave amplitudes; a statistically steady inertial range with negligible forcing and dissipation within it; transfer mainly between comparable perpendicular scales; and a scale-independent energy flux per unit mass $\mathcal{E}$. The estimate also assumes no extra scale-dependent alignment factor in the nonlinear interaction and neglects intermittency corrections.

Let ${\ell }_{\perp }$ and ${\ell }_{\parallel }$ be an eddy's perpendicular and parallel scales, and $\delta u_{\ell }\sim \delta B_{\ell }/\sqrt{4\pi {\rho }_0}$ its Alfvénic velocity amplitude. The nonlinear interaction is supplied by the oppositely propagating Elsässer variable, so comparable populations giveThe Alfvén propagation time is ${\tau }_A\sim {\ell }_{\parallel }/v_A$. For strong turbulence, critical balance makes these times comparable; the decorrelation and cascade time is therefore of order ${\tau }_{\mathrm{nl}}$, not a much longer time obtained by adding many weak encounters.

A constant perpendicular energy flux givesDefine the one-dimensional perpendicular turbulent energy spectrum so that $\int E(k_{\perp })d{k}_{\perp }$ measures fluctuation energy per unit mass. Locality of the spectrum then givesHencewhere $C$ is a dimensionless normalization constant not determined by the scaling argument. The magnetic and kinetic spectra have the same leading scaling under the assumed Alfvénic amplitude relation.

Critical balance also determines the scale-dependent anisotropy:Thus the eddies become relatively more elongated along the guide field at smaller perpendicular scales. The $-5/3$ exponent follows under the listed cascade assumptions; it is not established for every imbalanced or aligned MHD state solely by invoking critical balance.

Let $C(r)=⟨\mathbf{u}(\mathbf{x})\cdot \mathbf{u}(\mathbf{x}+\mathbf{r})⟩$, with zero mean velocity, and normalize the turbulent energy spectrum by ${\int }_0^{\infty }E(k)dk=⟨|\mathbf{u}{|}^2⟩/2$. Use the Fourier transform pairIsotropic turbulence makes the transform radial, while the energy normalization gives $E(k)=2\pi k^2\hat{C}(k)$. ThusAssume, for example, $\int (1+r^4)|C(r)|d^{3}r<\infty$, sufficient to control the remainder. Expanding the radial kernel gives $\sin (kr)/(kr)=1-k^2{r}^2/6+O(k^4{r}^4)$. With the Saffman integral $L=\int C(r)d^{3}r$ and Loitsyansky integral $I=-\int r^{2}C(r)d^{3}r$, this provesThe original PDF defines the two distinct integrals; the TeX aid accidentally repeats the definition of $I$ in place of $L$.

For a growing volume $V$, let ${\mathbf{P}}_V={\int }_V\mathbf{u}dV$ be its momentum per unit density. Statistical homogeneity givesFor volumes whose boundary-to-volume ratio tends to zero and integrable $C$, divide by $|V|$ and use dominated convergence:If the volume contains many weakly dependent eddies with uncompensated random impulses, the central limit theorem suggests root-mean-square total momentum proportional to $|V{|}^{1/2}$; its variance is extensive and $L$ is nonzero. This is a physical expectation, not a theorem for all turbulent fields: cancellation between cells can eliminate the extensive term.

In particular $L=0$ for an ensemble built from locally momentum-compensated eddies with zero net impulse, or when the complete incompressible velocity correlation tensor decays rapidly enough to have no longitudinal $r^{-3}$ tail. Merely imposing zero total momentum in one finite periodic box is insufficient to prove $L=0$ in an infinite-volume limit: the finite-box constraint fixes one global mode. Rapid decay of the trace $C$ alone is also insufficient; the individual tensor components can retain algebraic tails whose trace cancels.

In this question $\Delta v$ is the velocity increment in the direction of the separation, so $S_2(r)=⟨(\Delta v{)}^2⟩$ is a longitudinal velocity structure function. This interpretation is necessary for the kernel in part (iii); the full vector increment has a different kernel.

A velocity component that varies very little over distance $r$ is almost the same at the two sample points, so its common part cancels from their difference. A Fourier component acquires the difference factor $e^{i\mathbf{k}\cdot \mathbf{r}}-1$. Its squared magnitude is $2[1-\cos (\mathbf{k}\cdot \mathbf{r})]$, of order $(kr{)}^2$ when $kr\ll 1$ and of order one after averaging rapid phases when $kr\gg 1$. Thus a longitudinal structure function attenuates large eddies and retains smaller eddies, but is a smooth filter rather than an exact scale cutoff.

The heuristic cutoff $k_c=\pi /r$ marks wavelengths of order $2r$. If motions below that cutoff are treated as perfectly coherent and those above it as uncorrelated at the two points, isotropy allocates one third of their velocity variance to the longitudinal component. Since the turbulent energy spectrum integrates to one half of total velocity variance, the uncorrelated high-wavenumber component contributesorThe factor $3/4$ follows from the component and energy normalizations. The approximation discards the finite gradient of eddies larger than $r$; part (iii) quantifies that missing contribution.

Let $E(k)$ be the one-dimensional turbulent energy spectrum, normalized so that $\frac{1}{2}⟨|{\mathbf{u}}^{\prime }{|}^2⟩={\int }_0^{\infty }E(k)dk$. Large energy-containing eddies receive energy from the forcing. In the inertial range, nonlinear transfer passes this energy towards smaller scales at a mean rate $ϵ$ per unit mass. Dimensional analysis of Kolmogorov 1941 theory givesThe dissipation range begins near the inverse of the Kolmogorov length scale, where kinematic viscosity removes the transferred kinetic energy. The large-scale spectrum depends on the forcing; the rising low-$k$ branch in this sketch is illustrative, rather than a universal power law.

The tube has neither homogeneous nor isotropic large-scale turbulence: gravity selects a vertical direction, its walls constrain the motion, and its density gradient varies with height. The width $W$ bounds the transverse size of the energy-containing eddies. At sufficiently high Reynolds number, much smaller eddies may become approximately locally isotropic and exhibit an inertial range, even though the full tube flow does not satisfy the homogeneous-isotropic assumptions.

To obtain the Reynolds stress, write $u_i=U_i+u_i^{\prime }$ and $p=P+p^{\prime }$, with $⟨u_i^{\prime }⟩=0$. Averaging the incompressible Navier-Stokes equations uses $⟨u_i{u}_j⟩=U_i{U}_j+⟨u_i^{\prime }{u}_j^{\prime }⟩$ and gives the Reynolds-averaged momentum equationHere the final term is the mean buoyancy force, with the hydrostatic reference removed. Thus averaging the quadratic momentum flux introduces the covariance $⟨u_i^{\prime }{u}_j^{\prime }⟩$; multiplied by ${\rho }_0$, it has the dimensions of a stress. Its divergence transfers momentum between the mean flow and the fluctuations.

An eddy viscosity model represents this Reynolds stress byThe isotropic term can be absorbed into the mean pressure, while the deviatoric term gives a diffusion of mean momentum with turbulent viscosity ${\nu }_T\sim u_{\ast }\ell$, where $u_{\ast }$ and $\ell$ are representative turbulent speed and mixing length. Unlike kinematic viscosity, this coefficient depends on the flow. The closure assumes a local relation between stress and mean strain; anisotropic, buoyancy-driven large eddies may require a more detailed Reynolds stress model.

Write $S_2(r)=⟨(\Delta v{)}^2⟩$ and normalize the turbulent energy spectrum by ${\int }_0^{\infty }E(k)dk=⟨|\mathbf{u}{|}^2⟩/2$. In isotropic incompressible turbulence the longitudinal velocity correlation isCombining $S_2=2[Q_{LL}(0)-Q_{LL}(r)]$ with $Q_{LL}(0)=2\int Edk/3$ gives the exact spectral filterThe expression has a smooth removable limit at $x=0$. Its asymptotic forms areAn eddy much larger than $r$ contributes nearly the same velocity at the two points and mostly cancels in the increment. An eddy much smaller than $r$ contributes decorrelated velocities and does not cancel. Replacing the smooth transition of $H$ by a step at $kr=\pi$ therefore gives the common estimate $3S_2/4\simeq {\int }_{\pi /r}^{\infty }E(k)dk$. The cutoff corresponds roughly to a half wavelength equal to the point separation; its exact numerical location is conventional.

Keeping the leading nonzero low-wavenumber term gives the improved estimateIt remains imperfect because $kr$ is not small throughout the whole low-wavenumber interval and the exact kernel oscillates above the transition. The data sheet's alternative coefficient $1/{\pi }^2$ is a nearby approximate matching prescription; $1/10$ is the exact Taylor coefficient. Neither piecewise filter is an exact identity.

Higher-order longitudinal structure functions also retain large-scale information. Decompose $\mathbf{u}=\mathbf{U}+{\mathbf{u}}_s$, where $\mathbf{U}$ varies on a length much larger than $r$. Its longitudinal increment is $\Delta U_L\simeq r{\partial }_L{U}_L$, which is small but not zero. Expanding $(\Delta U_L+\Delta u_{s,L}{)}^p$ gives its own moments and mixed moments with the small-scale increment. Rare large gradients or correlated amplitude modulation can matter particularly strongly at large $p$.

Thus a velocity increment is not an exact scale-local bandpass filter. Its moments may retain large-scale strain and forcing information, undermining an unqualified claim of universal small-scale statistics determined by one local parameter. This observation alone does not refute an asymptotic inertial range theory: for a three-dimensional $k^{-5/3}$ spectrum, the low-wavenumber correction is mainly supplied by wavenumbers of order $1/r$, while an isolated integral-scale smooth contribution behaves as $r^2$ and is subleading to $r^{2/3}$. The issue is whether large-scale contamination and mixed moments actually become negligible in the regime and moment order being used.

The turbulent energy spectrum $E(k)$ distributes kinetic energy over wavenumber shells, normalized by ${\int }_0^{\infty }E(k)dk=⟨|\mathbf{u}{|}^2⟩/2$. For isotropic incompressible flow its tensor form is ${\Phi }_{ij}=E(k)({\delta }_{ij}-k_i{k}_j/k^2)/(4\pi k^2)$. Its low-wavenumber behavior records large-eddy momentum statistics.

The velocity correlation tensor is $Q_{ij}(\mathbf{r})=⟨u_i(\mathbf{x})u_j(\mathbf{x}+\mathbf{r})⟩$. For isotropic turbulence it splits into longitudinal and transverse parts. Incompressibility constrains them, while the Fourier transform gives the turbulent energy spectrum. Its trace is the scalar velocity-dot-product correlation.

For balanced anisotropic Alfvénic turbulence, suppose ${\tau }_A=l_{\parallel }/v_A\ll {\tau }_{\mathrm{nl}}=l_{\perp }/\delta z$ and resonant interactions preserve the parallel scale during the perpendicular cascade. Random encounter accumulation gives ${\tau }_{\mathrm{cas}}=v_A{l}_{\perp }^2/(\delta z^2{l}_{\parallel })$. Constant energy flux then gives $\delta z\propto l_{\perp }^{1/2}$ and the displayed perpendicular turbulent energy spectrum. Unlike isotropic phenomenology, the fixed parallel scale is independent of the decreasing perpendicular scale.