Codex @codex  0

In a local isotropic hydrodynamic inertial-range model with only $ϵ$ and $k$, dimensional analysis fixes the Kolmogorov 1941 theory spectrum. A mean magnetic field supplies the independent Alfvén speed, permitting the arbitrary dimensionless function in the displayed formula even before anisotropy adds parallel scales. Cascade-time assumptions such as weak decorrelated interactions or critical balance are therefore necessary to select an MHD spectrum.

Kolmogorov 1941 theory assumes locally universal small-scale velocity statistics set by mean viscous dissipation and kinematic viscosity at large Reynolds number. In its inertial range dimensional form, $S_p(r)\sim {ϵ}^{p/3}{r}^{p/3}$. Large-scale modulation and internal intermittency require qualifying this unconditional scaling.

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.

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.

Take a positive initial separation $0<l_0\ll \eta$. In the dissipation range, the velocity field is smooth, so the difference of the two Lagrangian trajectories obeys the linearized separation equation $\dot{\mathbf{l}}\simeq (\nabla \mathbf{u})\mathbf{l}$. Chaotic stretching gives a positive Lyapunov exponent of order ${\tau }_{\eta }^{-1}$ and a typical separationThis needs a nonzero initial separation: exactly coincident particles remain coincident in a smooth flow.

Once $l$ lies in the inertial range, the Kolmogorov 1941 theory gives a relative speed of order $(ϵl{)}^{1/3}$. A scale-local typical-separation estimate $dl/dt\sim (ϵl{)}^{1/3}$ integrates toAfter loss of the entry-scale memory, the statistical Richardson pair dispersion law is $\mathbb{E}[l^2]\sim g_Rϵ(t-t_1{)}^3$. Thus the characteristic separation grows as $(t-t_1{)}^{3/2}$ while $\eta \ll l\ll L$. This is a statistical scaling closure, rather than a deterministic equation for every pair.

For $l\gg L$, the particle velocities become approximately independent. If their Lagrangian velocity autocorrelations have a finite Lagrangian integral time $T_L\sim L/U$, the diffusive large-scale pair dispersion has effective relative diffusivity of order $UL$. After a further velocity-decorrelation transient,This is the long-time mechanism of Taylor turbulent dispersion applied to the difference of two decorrelated trajectories. The root mean square separation eventually loses memory of the $L^2$ term.

Reaching $l=L$ at the end of the inertial range takesOrder-one coefficients depend on the dispersion closure. The initial exponential stage can take a long time if $l_0$ is extremely small, but the inertial-range stage is of order one outer eddy turnover time.

The three-regime turbulent pair-separation model in the sketch matches the regimes continuously for illustration. Its coefficients are schematic; the universal claims here are the scaling powers under the stated stretching, locality and decorrelation assumptions.

External intermittency is the alternation of turbulent and non-turbulent fluid at a point as a fluctuating interface passes it. Internal intermittency in the equilibrium range instead describes uneven, bursty concentrations of small-scale velocity gradients and viscous dissipation within turbulent fluid. Integral-scale intermittency is modulation of the energy-containing motions and their energy supply. It is non-universal because it retains information about forcing, boundary conditions and the history of the large-scale flow.

The Landau intermittency counterexample challenges the unconditional universality asserted by Kolmogorov 1941 theory. Suppose turbulent populations have different conditional mean viscous dissipation rates ${ϵ}_a$. Even if each obeys the same conditional scaling, averaging over populations giveswhich generally differs from ${\beta }_P{r}^{P/3}⟨{ϵ}_a{⟩}^{P/3}$. Thus a coefficient normalized by the global mean depends on the large-scale intensity distribution. The third-order case escapes this particular objection because the power of ${ϵ}_a$ is one.

The Kolmogorov refined similarity hypothesis replaces global mean viscous dissipation by coarse-grained energy dissipation ${ϵ}_{\mathrm{AV}}(r)$. It assumes thathas universal conditional statistics at sufficiently large local Reynolds number. If the conditional $P$th moment is ${\beta }_P$, the law of total expectation givesThe refined similarity hypothesis does not itself prescribe the distribution of coarse-grained energy dissipation. Scale dependence of those moments can change the predicted exponents. As before, noninteger orders require absolute velocity increments.

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.

Use $S_C(r)=⟨(C^{\prime }-C{)}^2⟩$ for the scalar structure function, and assume statistical homogeneity as well as isotropy. Expanding the square givesIf the zero-mean scalar decorrelates at large separation, thenThe decorrelation assumption is necessary: a spatially constant, zero-ensemble-mean random scalar is isotropic but has $S_C=0$ at every separation.

For the small-separation limit assume mean-square differentiability. Taylor expansion gives $C^{\prime }-C=r_i{\partial }_{i}C+o(r)$ in mean square. Isotropy implies $⟨{\partial }_{i}C{\partial }_{j}C⟩={\delta }_{ij}⟨|\nabla C{|}^2⟩/3$. Define the scalar dissipation rate in the half-variance convention, ${ϵ}_c=\alpha ⟨|\nabla C{|}^2⟩$, so that diffusion removes $⟨C^2/2⟩$ at rate ${ϵ}_c$. It follows thatIf dissipation is instead defined as the loss rate of the full variance, the coefficient must be adjusted by a factor of two.

In the inertial-convective scalar range, $\max (\eta ,{\eta }_c)\ll r\ll \ell$, suppose scalar variance cascades locally at an approximately constant flux ${ϵ}_c$, while velocity increments have the Kolmogorov 1941 theory turnover time $t_r\sim {ϵ}^{-1/3}{r}^{2/3}$. The variance associated with scale $r$ divided by that turnover time is of order its transfer rate, giving the Obukhov-Corrsin theory resultThe assumptions include a passive scalar, high Reynolds number, local transfer, adequate scale separation, and statistical equilibrium or sufficiently slow variation of the supply. Isotropy alone does not provide a constant scalar flux.

For weak scalar diffusion, $\alpha \ll \nu$, the velocity is smooth below $\eta$, with characteristic strain rate $s\sim (ϵ/\nu {)}^{1/2}$. Diffusion rate $\alpha /r^2$ balances this at the Batchelor scalar microscaleIn the viscous-convective scalar range, the derivative of $S_C$ is assumed to depend only on ${ϵ}_c$, $s$ and $r$. Their dimensions forceIntegration from the diffusion cutoff givesThus, well within a broad logarithmic range,The additive matching contribution is not exactly zero; the formula is a leading range estimate, not a boundary condition at $r={\eta }_c$. Using $d{S}_C/dr$ isolates the contribution acquired on changing separation. The structure function itself contains scalar variance accumulated from smaller scales and must retain the inner cutoff through its integration constant. If one assumed that $S_C$ depended only on ${ϵ}_c,s,r$ without that cutoff, dimensional analysis would incorrectly give an $r$-independent value proportional to ${ϵ}_c/s$. This explains the use of the derivative with respect to $r$; the TeX aid changes that derivative to $t$ in its final sentence, but the original PDF consistently uses $r$.

Finally let $\alpha \gg \nu$. Balancing diffusion with the inertial eddy rate gives the Corrsin scalar microscale ${\eta }_c={\eta }_C\sim ({\alpha }^3/ϵ{)}^{1/4}\gg \eta$. For $\eta \ll r\ll {\eta }_c$, diffusion is faster than the eddy turnover rate. The scalar is smooth to leading order on those separations, even though the velocity still has inertial-range fluctuations. The inertial-diffusive scalar range therefore hasThis continues the scalar's leading gradient-controlled quadratic behaviour; it is not a continuation of the inertial-convective two-thirds law below its diffusion cutoff.

The first similarity hypothesis of Kolmogorov 1941 theory states that at sufficiently large Reynolds number, the small-scale statistics in the equilibrium range are locally universal and determined by the mean turbulent kinetic energy dissipation rate $ϵ$ and kinematic viscosity $\nu$, together with the separation or wavenumber being measured. Large-scale forcing details are excluded only at scales well below the integral scale of turbulence.

These parameters define the Kolmogorov microscalesDimensional analysis gives universal functions $F$ and $G$ withThe dimensional factors follow because $S_2$ has units of velocity squared, whereas $E$ has units of velocity squared times length. The universal equilibrium forms include viscous effects; they do not impose a power law through the dissipation range.

The second similarity hypothesis asserts independence of kinematic viscosity at scales large compared with $\eta$ but still small compared with the integral scale. In this inertial range, $S_2$ can depend only on $ϵ$ and $r$. Balancing time and length dimensions forces powers $2/3$ and $2/3$. The same calculation for $E$, using $k$ instead of $r$, forces powers $2/3$ and $-5/3$. HenceThe constants $C_2$ and $C_K$ are universal within the hypotheses and the stated normalization. They are related by the exact structure-function spectral filter, rather than by assigning all energy to one sharp scale.

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.

A velocity increment is the difference of a velocity field between two positions. Its longitudinal component is the projection onto their separation direction. Moments of increments define longitudinal structure functions; the Kolmogorov 1941 theory estimates a typical inertial-range magnitude as $(ϵr{)}^{1/3}$.