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 and introduce the fluctuation Elsässer variables . Their equations are
Only 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 denote the typical balanced amplitude and let be coherence lengths relative to the local mean field. The Alfvén wave time and eddy turnover time are
For , one counterpropagating encounter changes the packet fractionally by order . Assume phases decorrelate between encounters. After encounters, random accumulation is of order ; an order-one change requires . The weak wave cascade time is therefore
Weakness 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, , in addition to . Then and
Hence . A local shell-integrated turbulent energy spectrum has , , so
The 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 approximately fixed while cascading to smaller , and require . 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 . For the one-dimensional perpendicular spectrum, gives
Unlike the isotropic estimate, the parallel length does not decrease with the perpendicular length. Since , the interactions become stronger down the cascade. If at outer perpendicular scale , the weak approximation breaks at . 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, . The local cascade is strong, with transfer time of order one eddy turnover time. Then
The parallel scale is set dynamically, not held fixed:
At smaller scales grows, so the eddies are increasingly elongated along the local magnetic field. Strong turbulence does not mean : 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 . With and , writing forces , . 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 , and anisotropy adds independent parallel lengths. Even before adding anisotropy, dimensions permit
with arbitrary dimensionless . For example 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 . 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 , let be a typical velocity increment. The eddy turnover time is . Constant transfer rate then gives , hence
Here is the one-dimensional turbulent energy spectrum, obtained from with . The scaling holds only for
The Kolmogorov length scale follows by setting the scale-dependent Reynolds number to one. The associated Kolmogorov microscales give . At larger scales the forcing matters; at smaller scales kinematic viscosity matters.
Take a positive initial separation . In the dissipation range, the velocity field is smooth, so the difference of the two Lagrangian trajectories obeys the linearized separation equation . Chaotic stretching gives a positive Lyapunov exponent of order and a typical separation
This needs a nonzero initial separation: exactly coincident particles remain coincident in a smooth flow.
Once lies in the inertial range, the Kolmogorov 1941 theory gives a relative speed of order . A scale-local typical-separation estimate integrates to
After loss of the entry-scale memory, the statistical Richardson pair dispersion law is . Thus the characteristic separation grows as while . This is a statistical scaling closure, rather than a deterministic equation for every pair.
For , the particle velocities become approximately independent. If their Lagrangian velocity autocorrelations have a finite Lagrangian integral time , the diffusive large-scale pair dispersion has effective relative diffusivity of order . 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 term.
Reaching at the end of the inertial range takes
Order-one coefficients depend on the dispersion closure. The initial exponential stage can take a long time if is extremely small, but the inertial-range stage is of order one outer eddy turnover time.
Figure 1.
Exponential, Richardson and diffusive regimes of turbulent pair separation
. Details.
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.
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 , mean viscous dissipation and kinematic viscosity . Dimensional analysis therefore expresses the second-order longitudinal structure function in terms of the Kolmogorov microscales:
The universality of is an assumption of the theory, not a consequence of dimensional analysis alone.
The Kolmogorov second similarity hypothesis removes dependence on kinematic viscosity when . Only and remain, so the Kolmogorov two-thirds law is
The same reasoning gives the predicted order- longitudinal structure function:
These are signed moments for integer ; for noninteger orders, use . In particular, the signed third-order coefficient is negative, , in the Kolmogorov four-fifths law.
The zeroth law of turbulence states that at fixed outer velocity and integral scale of turbulence , the mean viscous dissipation per unit mass approaches a nonzero finite value as the Reynolds number becomes large. In an established forward energy cascade, this gives , with of order one. The small-scale velocity gradients become large enough to offset decreasing kinematic viscosity; this is the turbulent dissipation anomaly.
At the Kolmogorov microscales, the local Reynolds number is of order one, , and the viscous dissipation is . Solving these balances gives
Substituting the outer-scale estimate, and omitting order-one constants, gives
The associated eddy turnover time is .
A passive scalar is transported by the velocity field without feeding back on it:
Take its mean to be zero, and define the scalar dissipation rate by the half-scalar variance convention . In a statistically equilibrated energy cascade of scalar fluctuations, this is also the flux of half-scalar variance toward small scales.
The Kolmogorov two-thirds law gives a typical velocity increment and hence eddy turnover time . Assuming local scalar transfer on this same time scale, . The Obukhov-Corrsin theory therefore gives the scalar analogue:
This applies to separations at which both direct forcing and molecular diffusion are negligible and the transporting velocity increments lie in the inertial range. For very large or very small ratios , the scalar and velocity cutoff scales differ, so this common range must be checked. Defining as dissipation of the full scalar variance instead simply changes the convention for .
For positive initial separation below the Kolmogorov length scale, a phenomenological pair-dispersion model has smooth-flow exponential stretching, inertial-range Richardson pair dispersion, and diffusive large-scale pair dispersion beyond the integral scale of turbulence. A positive Lyapunov exponent is needed for the first regime, scale locality and initial-memory loss for the second, and integrable relative-velocity autocorrelation for the third. The inertial-range crossing time is of order , one outer eddy turnover time.