Past exam of the mathematics course of the University of Cambridge / 2005 / iii / Paper 76 / 2

In mean-field electrodynamics, use Reynolds averaging to write ${\mathbf{U}}_{\mathrm{tot}}=\overline{\mathbf{U}}+\mathbf{u}$ and ${\mathbf{B}}_{\mathrm{tot}}=\overline{\mathbf{B}}+\mathbf{b}$, with zero-mean fluctuations. Assume the averaging operation commutes with the derivatives and that resolved quantities vary little across an averaging region. Averaging the induction equation then givesThis is the mean-field electromotive force, a vector correlation with units of velocity times magnetic field. The closure problem is to express this correlation in terms of the resolved field.

In the kinematic regime the flow statistics are prescribed, so the response is linear in the mean magnetic field, after independently excited fluctuation transients have decayed. If the mean-field scale $L$ is much greater than the turbulent scale $\ell$, and its evolution is slow compared with the response memory, expand that response in the local mean magnetic field and its gradients:Without scale separation one generally has a spatially nonlocal and temporally retarded response rather than constant local coefficients.

For homogeneous turbulence the coefficients are independent of position. Invariance under proper rotations makes the rank-two response proportional to ${\delta }_{ij}$, and the isotropic rank-three gradient response proportional to the alternating tensor. Choosing the diffusivity sign convention givesThe alpha effect can regenerate a mean magnetic field, while $\beta$ supplies turbulent magnetic diffusivity. For constant positive $\beta$, $\nabla \times (-\beta \nabla \times \overline{\mathbf{B}})=\beta {\nabla }^2\overline{\mathbf{B}}$ because the magnetic field is solenoidal. Thus the total diffusivity is $\eta +\beta$.

The electromotive vector is a polar vector, whereas the magnetic field is an axial vector. Therefore $\alpha$ is a pseudoscalar and changes sign under reflection; $\beta$ is an ordinary scalar. A nonzero isotropic alpha effect requires broken reflection symmetry, such as nonzero kinetic helicity. If “isotropic” includes reflection symmetry as well as rotations, $\alpha$ must vanish. The dimensions are $[\alpha ]=L/T$ and $[\beta ]=L^2/T$.

Take zero mean flow locally and incompressible fluctuations. Subtracting the averaged induction equation givesThe first-order smoothing approximation omits the second term on the right when solving for $\mathbf{b}$, but retains $⟨\mathbf{u}\times \mathbf{b}⟩$ in the mean equation. For a locally uniform, slowly changing mean magnetic field the forcing becomes $(\overline{\mathbf{B}}\cdot \nabla )\mathbf{u}$.

There are two controlled ways to make the omitted fluctuation term small. At small eddy magnetic Reynolds number $R_m=u_{\mathrm{rms}}\ell /\eta$, diffusion gives $b/B\sim R_m$, so the omitted $ub/\ell$ is smaller than the retained $uB/\ell$ by $R_m$. Alternatively, if $u_{\mathrm{rms}}\tau /\ell \ll 1$ for the fluctuation response time $\tau$, then $b/B\sim u_{\mathrm{rms}}\tau /\ell$ and the short-memory response is perturbative. Mean-field scale separation is additionally required. Merely calling the turbulence weak or assuming that the mean magnetic field is small does not justify the omission.

With homogeneous initial transients removed, the causal linear solution isInsert this into the electromotive correlation. Isotropy gives ${\alpha }_{ij}=\alpha {\delta }_{ij}$, so taking the trace yields the helicity formula for isotropic first-order smoothing:Using the alternating-tensor sign givesThe minus sign follows from ${ϵ}_{imn}=-{ϵ}_{min}$. This derivation specifies the diffusion and correlation weighting, rather than asserting that alpha is proportional to instantaneous helicity for arbitrary turbulence.

If $\eta \tau /{\ell }^2\ll 1$, diffusion over a correlation time is negligible. Define the helicity correlation time byFor an ordinary decaying correlation with nonzero mean kinetic helicity, this gives the familiar estimateIt is the kinetic helicity, not magnetic helicity. Opposite flow handedness reverses the alpha effect. The time-integral formula remains meaningful if the equal-time helicity vanishes or its correlation changes sign, when replacing it by a single positive correlation time may fail. In the same weak-diffusion short-memory approximation the gradient forcing $-(\mathbf{u}\cdot \nabla )\overline{\mathbf{B}}$ gives $\beta =\frac{1}{3}{\int }_0^{\infty }⟨\mathbf{u}(t)\cdot \mathbf{u}(t-\tau )⟩d\tau \simeq {\tau }_U{u}_{\mathrm{rms}}^2/3$; the two correlation times need not be equal.

In the small-$R_m$, slowly varying quasistatic regime, diffusion is instead essential. The response givesFor velocity concentrated at wavenumber $k_0$, this becomes $\alpha =-⟨\mathbf{u}\cdot \mathbf{\varpi }⟩/(3\eta k_0^2)$. Its effective memory is the magnetic diffusion time on the eddy scale. The controlled low-$R_m$ setting is developed in Moffatt's analysis of turbulent dynamo action; the causal-response calculation above also covers short-memory turbulence without assuming that diffusion is always negligible.

For constant isotropic $\alpha ,\beta$ and zero mean flow, the mean induction equation isLet $\overline{\mathbf{B}}=\hat{\mathbf{B}}{e}^{\sigma t+i\mathbf{k}\cdot \mathbf{x}}$, with $\mathbf{k}\cdot \hat{\mathbf{B}}=0$. On this two-dimensional transverse space, $i\mathbf{k}\times$ has eigenvalues $\pm k$. For example, if $\mathbf{k}=k{\mathbf{e}}_z$, the vectors ${\mathbf{e}}_x\pm i{\mathbf{e}}_y$ have these curl eigenvalues. The homogeneous alpha-squared dynamo growth criterion follows fromChoose the circular polarization of the mode to match the sign of $\alpha$. This is an alpha-squared dynamo: the alpha effect supplies both parts of the regeneration loop without differential rotation.

If arbitrarily long wavelengths are permitted, every nonzero alpha has a growing band $0<k<|\alpha |/{\eta }_T$. Maximizing the quadratic givesprovided $k_{\ast }\ell \ll 1$ so the mean-field approximation is valid. A uniform field has $k=0$, zero curl and zero growth; a constant alpha cannot amplify that uniform mode. In a periodic domain of size $L$, the lowest nonzero wavenumber is $2\pi /L$, giving $|\alpha |L/{\eta }_T>2\pi$. Other boundaries and geometries give different critical constants, but the general scale criterion is $|\alpha |L/{\eta }_T$ sufficiently large. Nonzero kinetic helicity alone is therefore not a sufficient finite-domain dynamo criterion: its induced alpha must overcome diffusion on an available mean-field scale.

The geodynamo operates in the electrically conducting liquid of Earth's outer core. buoyancy from thermal and compositional convection supplies motion, while the Coriolis force organizes rotating flow. Rotation together with radial stratification permits kinetic helicity, typically with opposite signs in the two hemispheres. In mean-field language this allows an alpha effect that converts poloidal and toroidal components into one another. A large-scale alpha-squared dynamo is therefore a plausible regeneration mechanism even without dominant differential rotation.

The growth condition compares this regeneration with the sum of molecular and turbulent magnetic diffusion on the core scale. A spherical model must impose the appropriate magnetic matching at the core boundaries and permit spatially varying, generally anisotropic alpha and diffusivity tensors. Rotation selects an axis, and stratification and boundaries select positions, so homogeneous isotropic turbulence is only a local idealization. Such a model can favor a predominantly axial dipole, but the preferred parity and threshold must be found from its global eigenproblem rather than deduced from the plane-wave formula alone. If mean shear is significant, alpha–Omega regeneration can also contribute.

The first-order-smoothing helicity formula is a useful guide to the sign and magnitude of alpha, not an automatically controlled numerical closure for the core. The important convective motions need not have small eddy magnetic Reynolds number or short enough correlation time; scale separation can also be limited. Furthermore, the Lorentz force modifies the flow and its correlations once the field grows. Saturation therefore requires nonlinear feedback, including dynamo quenching and magnetic-helicity constraints where relevant.

The observed field varies and occasionally reverses. A constant-coefficient growing kinematic eigenmode does not itself explain these reversals or select a saturated amplitude. Spatial and temporal variations, competing modes, and nonlinear flow–field interactions must be included. mean-field dynamo theory supplies a physical regeneration mechanism and a testable threshold, while the geodynamo requires a rotating, bounded, nonlinear realization of that mechanism.