Past exam of the mathematics course of the University of Cambridge / 2007 / iii / Paper 75 / 3

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.

At leading anisotropic order, perpendicular Alfvénic turbulence evolves independently of the compressive fluctuations. The slow modes and the entropy mode are transported by that turbulence. The fast magnetosonic frequency is of order $k_{\perp }(v_A^2+c_s^2{)}^{1/2}$, much larger than the reduced Alfvén frequency; the low-frequency reduced magnetohydrodynamics ordering filters out that fast-wave sector.

The passive evolution can be made explicit for a uniform ideal adiabatic plasma. Let $c_s^2=\gamma p_0/{\rho }_0$, $\delta =\delta \rho /{\rho }_0$, $b=\delta B_{\parallel }/B_0$, and define the normalized entropy perturbationUse $D={\partial }_t+{\mathbf{u}}_{\perp }\cdot {\nabla }_{\perp }$ and ${\nabla }_{\parallel }={\partial }_z+(\delta {\mathbf{B}}_{\perp }/B_0)\cdot {\nabla }_{\perp }$. Perpendicular force balance makes the leading total-pressure perturbation vanish:The ideal entropy equation gives $Ds=0$. The leading continuity and parallel-induction equations have the same compressive divergence term; subtracting them gives $D(\delta -b)=-{\nabla }_{\parallel }{u}_{\parallel }$. Parallel momentum gives $D{u}_{\parallel }=-{\nabla }_{\parallel }\delta p/{\rho }_0$. PutThe two compressive equations therefore becomeLet $c_T=v_A{c}_s/(v_A^2+c_s^2{)}^{1/2}$ be the tube speed and $h=c_s^{2}r/c_T$. ThenThus the slow and entropy fluctuations in reduced magnetohydrodynamics have characteristic amplitudes $f^{\pm }=u_{\parallel }\pm h$ obeyingThe two slow-wave populations propagate in opposite directions along the bent field lines while being mixed across them by ${\mathbf{u}}_{\perp }$. The entropy perturbation is carried by the fluid without a restoring-wave frequency. A pure entropy perturbation has $r=b=u_{\parallel }=0$ and $\delta =-s/\gamma$, so its density variation is accompanied by a temperature change with no leading pressure variation.

The coefficients and advecting fields in these three scalar equations come from the Alfvénic cascade; the scalars do not enter the leading perpendicular Alfvén equations. They therefore develop small perpendicular scales through passive mixing. If a scalar $f$ has constant variance flux ${\mathcal{F}}_f$ and is mixed on the Goldreich–Sridhar time ${\tau }_{\ell }\sim {\mathcal{E}}^{-1/3}{\ell }_{\perp }^{2/3}$, thenThe flux and amplitude need not equal the Alfvénic energy flux. This passive-spectrum statement uses the same locality and stationary inertial-range assumptions as part (a), and ignores additional damping. If a passive mode is absent initially and is not forced, its homogeneous leading equation keeps it absent; the Alfvénic cascade does not have to generate every compressive mode.

In the ideal, leading anisotropic system, there are five separately conserved quadratic cascade quantities: the two oppositely propagating Alfvén populations, the two slow-wave populations, and the entropy fluctuation. With the amplitudes of part (b), one convenient normalization iswhere ${\mathbf{z}}_{\perp }^{\pm }={\mathbf{u}}_{\perp }\pm \delta {\mathbf{B}}_{\perp }/\sqrt{4\pi {\rho }_0}$. Constant background factors convert the velocity variances into modal energies and the entropy variance into an entropy-related free-energy contribution. These are the five quadratic cascade invariants of anisotropic magnetohydrodynamic turbulence.

To verify conservation rather than just count polarizations, the perpendicular Elsässer variables obeyDot each equation with $2{\mathbf{z}}_{\perp }^{\pm }$ and integrate. The propagation and nonlinear terms become surface fluxes of $|{\mathbf{z}}_{\perp }^{\pm }{|}^2$, and the pressure term becomes a surface flux because the corresponding field is solenoidal. For periodic boundaries, or boundaries with vanishing invariant flux, all these terms vanish and $d{I}_A^{\pm }/dt=0$. Their sum is proportional to perpendicular kinetic plus magnetic energy, while their difference is proportional to perpendicular cross-helicity; the two Elsässer energy invariants contain both pieces of information.

Likewise the slow-mode equations giveHere ${\mathbf{u}}_{\perp }$ is divergence-free and ${\nabla }_{\parallel }$ differentiates along the divergence-free vector $\hat{\mathbf{z}}+\delta {\mathbf{B}}_{\perp }/B_0$, so both integrals are boundary fluxes. Finally $Ds=0$ gives $d{I}_e/dt=-\int {\mathbf{u}}_{\perp }\cdot {\nabla }_{\perp }{s}^{2}dV=0$. This proves all five quadratic conservation laws under the same boundary assumptions.

The two Alfvén cascades interact dynamically: each population strains the other and is needed for its nonlinear transfer between scales. Their integrated invariants are nevertheless conserved separately, so this interaction does not transfer one population's invariant into the other. The slow and entropy cascades are driven by the Alfvénic mixing and field-line motion, with no feedback on the perpendicular cascade at the retained order. Their injection rates set independent fluxes; those fluxes need not be equal, even when the resulting perpendicular spectra share the $-5/3$ scaling.

The count is a count of independent quadratic modal cascade channels. Ideal passive advection also preserves higher moments, such as $\int F(s)dV$ for arbitrary smooth $F$, and analogous integrals for the slow characteristic scalars. These are not additional quadratic modal energies. Forcing and dissipative terms make the corresponding invariants source-and-sink balances; a stationary inertial range conserves their flux through scales rather than their total content in a continually forced domain.