Dynamo action

Dynamo action is growth or sustained maintenance of a magnetic field by conducting-fluid motion despite magnetic diffusion. In the kinematic problem the velocity is prescribed and the resistive induction equation is linear in the magnetic field. Magnetic feedback through the Lorentz force density matters once the field becomes dynamically important.

For an isolated bounded conductor of uniform magnetic diffusivity $\eta >0$, contained in a sphere of radius $R$ and matched to a decaying potential field in an insulating exterior, dynamo action requires maximum stretching rate $S\ge \eta {\pi }^2/R^2$. Here $S$ bounds the largest eigenvalue of the rate-of-strain tensor throughout the flow and time. Use fluid boundary conditions eliminating the stretching boundary term, for example a no-slip boundary condition, and no imposed energy input. The proof combines the magnetic free-decay spectral bound with the magnetic energy equation. This is a necessary condition, not a sufficiency criterion; changing magnetic boundary conditions changes the spectral constant.

Under the energy-closed boundary assumptions of Backus' necessary condition for dynamo action, the magnetic energy equation gives ${\dot{E}}_B\le 2s(t)E_B$ after discarding nonnegative resistive dissipation. Integrating gives an upper bound by the time average of the largest spatial eigenvalue of the rate-of-strain tensor. Thus the field-amplitude exponent is at most its space-time supremum. The exponent of squared energy is twice the field-amplitude exponent.

A mean-field dynamo evolves a slowly varying spatially averaged magnetic field. Correlations of fluctuating velocity and magnetic field supply a mean-field electromotive force. Its local linear response includes the alpha effect; gradients of the mean field can contribute additional transport terms. The distinction between a uniform test field used to measure a response coefficient and a spatially varying mean field used to obtain growth is essential.

An alpha-Omega mean-field dynamo couples toroidal-field generation by the Omega effect to poloidal-field regeneration by the alpha effect. A simple Fourier model is $\dot{A}={\alpha }_{0}f(t)B-\eta k^{2}A$, $\dot{B}=ik\Omega A-\eta k^{2}B$, where $A$ is a component of the vector potential for the poloidal field and $B$ is a toroidal-field amplitude. For bounded $|f|\le 1$, the bounded-modulation alpha-Omega growth estimate controls the maximum amplification.

In the two-component alpha-Omega dynamo with $|f(t)|\le 1$, apply the weighted energy estimate for two coupled modes with $a=|{\alpha }_0|$, $c=|k\Omega |$, and $d=\eta k^2$. The field-amplitude growth rate is at most $\sqrt{|{\alpha }_0\Omega k|}-\eta k^2$. Optimizing over $k$ gives the upper constant $3/2^{8/3}$ multiplying $({\alpha }_0^2{\Omega }^2/\eta {)}^{1/3}$. The admissible choice $f=1$ gives a lower constant $3/2^{10/3}$, so the maximum possible growth has the displayed scaling. An arbitrary chosen modulation need not grow. The statement concerns the model with freely adjustable wavenumber; physical scale-separation restrictions can constrain that optimization.

Differential rotation or shear stretches a poloidal magnetic field into a toroidal magnetic field. This is the Omega effect. By itself it need not regenerate the poloidal component; the alpha effect provides one possible feedback mechanism in an alpha-Omega dynamo. Without feedback, shear amplification can be transient or algebraic rather than sustained exponential dynamo action.

First-order smoothing neglects the fluctuating nonlinear velocity-magnetic-field product in the fluctuation resistive induction equation, but retains $⟨\mathbf{u}\times \mathbf{b}⟩$ in the mean-field electromotive force. For a constant test field it gives ${\partial }_t\mathbf{b}-\eta {\nabla }^2\mathbf{b}=({\mathbf{B}}_0\cdot \nabla )\mathbf{u}$. Its validity requires the omitted terms to be small compared with the retained forcing and diffusion terms.

For a velocity Fourier mode with wavenumber $k$, define $a=\eta k^2$, $d=a^2+{\Omega }^2$ and $\mathbf{C}=i({\mathbf{B}}_0\cdot \mathbf{k})\hat{\mathbf{u}}$. A cosine-forced mode has periodic response coefficients $(\mathbf{p},\mathbf{q})=(a\mathbf{C},\Omega \mathbf{C})/d$, while a sine-forced mode has $(-\Omega \mathbf{C},a\mathbf{C})/d$. Homogeneous transients decay as $e^{-at}$. The periodic response is a particular solution, not arbitrary initial data.

A Parker dynamo wave is a traveling or oscillatory mode of a mean-field dynamo coupling poloidal and toroidal magnetic-field components through an alpha effect and shear. A reduced complex-amplitude model can represent fluctuating alpha by a time-dependent coupling. Its periodically switched growth must be calculated from the monodromy matrix of a periodic linear system, rather than by averaging noncommuting generators.

For $\dot{A}=D^{2}f(t)B$ and $\dot{B}=iA$, with $f$ alternating between $+1$ and $-1$ for intervals of length $T$, the monodromy matrix of a periodic linear system has determinant one and trace $\cosh (\sqrt{2}DT)+\cos (\sqrt{2}DT)>2$. Its reciprocal Floquet multipliers give positive dominant Floquet growth rate for $D,T>0$, despite the zero mean of $f$. An exceptional initial state in the stable eigenspace decays.

The alpha effect is the part of a mean-field electromotive force linear in the mean magnetic field itself, $\mathcal{E}=\mathbf{\alpha }\overline{\mathbf{B}}+\cdots$. The alpha tensor can be anisotropic; it need not be a scalar multiple of the identity. Its curl can couple transverse mean-field components and produce an alpha-squared dynamo.

An alpha-squared dynamo generates a mean magnetic field using an alpha effect without requiring an additional large-scale shear coupling. For $\mathbf{\alpha }={\alpha }_0\mathrm{diag}(1,1,0)$ and a transverse mean field varying as $e^{iKz}$, the growth rates are $-\eta K^2\pm {\alpha }_{0}K$. A nonzero ${\alpha }_0$ therefore allows growth for sufficiently long wavelengths, provided the domain and scale separation permit them.

An alpha-squared dynamo may use a tensorial alpha effect rather than an isotropic coefficient. For constant alpha tensor $\mathbf{\alpha }$ and magnetic diffusivity $\eta >0$, a steady nonzero Fourier mode with wavevector $\mathbf{K}$ satisfies $\eta K^2\mathbf{B}=i\mathbf{K}\times (\mathbf{\alpha }\mathbf{B})$ and $\mathbf{K}\cdot \mathbf{B}=0$. Anisotropy can change both the critical alpha magnitude and the preferred wavevector direction.

For $\mathbf{\alpha }=\mathrm{diag}({\alpha }_0,{\alpha }_0,{\alpha }_1)$, set $\gamma ={\alpha }_0/\eta$, $\delta ={\alpha }_1/{\alpha }_0$, and $q^2=k^2+l^2$. A steady nonzero Fourier mode obeys the displayed condition when the denominator is nonzero. At fixed $q>0$ and $0<\delta <1$, minimizing over $m^2\ge 0$ gives $m_{\ast }^2=(1-2\delta )q^2$ and ${\gamma }_{\min }^2=4(1-\delta )q^2$ for $\delta <1/2$; for $\delta \ge 1/2$ it gives $m_{\ast }^2=0$ and ${\gamma }_{\min }^2=q^2/\delta$. This is a boundary-constrained threshold of an anisotropic alpha-squared dynamo.

The alpha tensor is the linear map from a uniform test magnetic field to the corresponding mean-field electromotive force. Its components satisfy ${\mathcal{E}}_i={\alpha }_{ij}{B}_{0j}$. Anisotropic helical flows can produce a real symmetric rank-two response, such as ${\alpha }_0\mathrm{diag}(1,1,0)$.

The mean-field electromotive force is the vector correlation $\mathcal{E}=⟨\mathbf{u}\times \mathbf{b}⟩$. Its curl enters the mean resistive induction equation. It is a local vector with units of velocity times magnetic field, not the scalar circuit voltage called electromotive force. A space-time average of periodic modes is often used to measure its transport coefficients.

A space-time average over periodic modes integrates over the spatial periodic cell and one temporal period. For real spatial Fourier modes sharing a wavevector, $⟨\mathrm{Re}(\mathbf{v}{e}^{i\mathbf{k}\cdot \mathbf{x}})\times \mathrm{Re}(\mathbf{w}{e}^{i\mathbf{k}\cdot \mathbf{x}}){⟩}_x=\mathrm{Re}({\mathbf{v}}^{\ast }\times \mathbf{w})/2$. A matching sine or cosine temporal factor supplies another $1/2$. Opposite wavevectors can have nonzero spatial cross averages, so their temporal coefficients must be checked rather than discarded merely because the wavevectors are distinct.