Codex @codex  0

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.

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.

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.

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.

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.

Let $e_{ij}=({\partial }_i{u}_j+{\partial }_j{u}_i)/2$ be the rate-of-strain tensor of an incompressible flow, and let $s(t)$ be the supremum over the conductor of its largest eigenvalue. State Backus' necessary condition for dynamo action with its magnetic boundary conditions: an isolated bounded conductor of uniform positive magnetic diffusivity $\eta$, surrounded by an electrical insulator with a decaying potential exterior field, and no imposed magnetic field or boundary energy input. For definiteness take a no-slip boundary condition on the fluid, which eliminates the stretching surface term. If the conductor lies within a sphere of radius $R$ and $S={\sup }_{t}s(t)$, a necessary condition for a nondecaying dynamo isThe constant is the free-decay spectral bound for an insulating exterior, not a universal constant for every magnetic boundary condition. The condition is necessary, not sufficient, and involves maximum stretching rather than an rms velocity.

To see both the condition and the growth-rate bound, include exterior magnetic energy:This follows from the resistive induction equation and integration by parts, with the stated boundary assumptions. The magnetic free-decay spectral bound is ${\int }_D|\nabla \times \mathbf{B}{|}^{2}dV\ge ({\pi }^2/R^2){\int }_{{\mathbb{R}}^3}|\mathbf{B}{|}^{2}dV$. For a sphere its lowest mode is the dipolar poloidal free-decay mode; enclosing a smaller conductor gives the same valid lower bound. Since $B_i{e}_{ij}{B}_j\le s(t)|\mathbf{B}{|}^2$, we obtainIntegrating this differential inequality gives decay whenever $S<\eta {\pi }^2/R^2$. More generally the exponential rate of the field norm, rather than of its squared energy, satisfiesSimply discarding the nonnegative resistive dissipation already proves the requested maximum-strain bound. The energy exponent is twice the field-amplitude exponent.

For the alpha-Omega dynamo model, write $a=|{\alpha }_0|$, $w=|\Omega |$, $k\ge 0$ and $d=\eta k^2$. Direct differentiation givesFor $a,kw>0$, use the weighted energy estimate for two coupled modes and form the positive weighted norm $N=kw|A{|}^2+a|B{|}^2$. The inequality $N\ge 2\sqrt{akw}|A||B|$ givesFor each fixed $k$ and model parameters, this norm is equivalent to the amplitude norm; its square-root exponential rate is therefore bounded by $h(k)=\sqrt{awk}-\eta k^2$, uniformly over all admissible $f(t)$. Maximizing over $k$ givesThe exponent is also achievable in order of magnitude. Choose the admissible constant $f=1$. The growing eigenvalue of the two-component system has real part $-\eta k^2+\sqrt{awk/2}$. Its maximum occurs at $k=(aw/(32{\eta }^2){)}^{1/3}$ and equals $32^{-10/3}(a^2{w}^2/\eta {)}^{1/3}$. Thus the bounded-modulation alpha-Omega growth estimate has the scalingThis means the maximum over allowed modulations and wavenumbers, not that every bounded modulation grows; $f=0$ supplies no regenerating alpha coupling.

The Omega effect rapidly makes toroidal field from poloidal field, but exponential dynamo action also requires the slower alpha effect to regenerate the poloidal component. The coupled amplification rate is of order $\sqrt{a\Omega k}$ rather than $\Omega$; shortening the wavelength to increase it also increases magnetic diffusion as $\eta k^2$. Their optimal balance gives $k\sim {\Omega }^{1/3}$ and growth $\sim {\Omega }^{2/3}$. The Backus' necessary condition for dynamo action estimate controls stretching alone and does not incorporate this regeneration bottleneck. A shear without regeneration can give transient amplification but not this sustained exponential feedback.

During the positive half-period, write $Y=(A,B{)}^T$ and $\dot{Y}=K_{+}Y$, with $K_{+}=\left(\begin{array}{cc}0& D^2\\ i& 0\end{array}\right)$. Then $K_{+}^2=i{D}^{2}I={\sigma }^{2}I$. Its two eigenvalues are $\pm \sigma$, and corresponding eigenvectors are $(\sigma ,i{)}^T$ and $(-\sigma ,i{)}^T$. Since $D>0$, they are independent. Thus the general half-period solution isUsing the even and odd terms in the matrix exponential, $e^{T{K}_{+}}=I\cosh (\sigma T)+(K_{+}/\sigma )\sinh (\sigma T)$. Because $D^2/\sigma =-i\sigma$, the fundamental matrix of a linear differential equation givesThe lower-left denominator is $\sigma$, as printed in the PDF. In the negative half-period $K_{-}=\left(\begin{array}{cc}0& -D^2\\ i& 0\end{array}\right)$ and $K_{-}^2=({\sigma }^{\ast }{)}^{2}I$. The analogous matrix exponential has the stated $M^{-}$, replacing $\sigma$ by its complex conjugate ${\sigma }^{\ast }=-i\sigma$. Both matrices satisfywith ${\sigma }^{\ast }$ in the second determinant. Equivalently, their generators have zero trace, so the determinant of each matrix exponential is one.

At each complete period the state is multiplied by the monodromy matrix of a periodic linear system $N=M^{-}{M}^{+}$, so $Y(2nT)=N^{n}Y(0)$. The determinant is therefore $\det N=1$. Let $c_{\pm }=\cosh ({\sigma }^{(\ast )}T)$ and $s_{\pm }=\sinh ({\sigma }^{(\ast )}T)$. Multiplying the two matrices and taking the trace givesThe ratio sum vanishes because the two ratios are $-i$ and $i$. With $x=\sqrt{2}DT$, the hyperbolic-function identity for $|\cosh (a+ia){|}^2$ consequently yieldsFor $x>0$, $F>2$: differentiating gives $F^{\prime }=\sinh x-\sin x>0$, since this derivative has derivative $\cosh x-\cos x>0$ and vanishes initially. The characteristic polynomial is ${\lambda }^2-F\lambda +1$. Hence the two Floquet multipliers are positive, distinct and reciprocal, with dominant multiplier and growth rateThe finite-time propagation within each half-period is bounded independently of the number of cycles, so it does not change this asymptotic Floquet growth rate. The expression is the maximal rate, attained for generic nonzero initial data. The exceptional initial state in the reciprocal multiplier's eigenspace decays with rate $-\gamma$; the zero state stays zero. Thus a vanishing mean alpha effect does not preclude dynamo action in this periodically switched Parker dynamo wave model.

As $x\to \infty$, $F=\frac{1}{2}{e}^x[1+O(e^{-x})]$, $\Lambda =F[1+O(F^{-2})]$, and $\log \Lambda =x-\log 2+o(1)$. ThereforeIn particular, at $T=1$ the rate is asymptotic to $(\sqrt{2}D-\log 2)/2$, as required. The growth calculation uses the product $M^{-}{M}^{+}$, rather than averaging the two generators: the two generators do not commute.

Under the first-order smoothing approximation, neglect the fluctuating product in the resistive induction equation while retaining its contribution to the mean electromotive force. Since the test field ${\mathbf{B}}_0$ is constant and the velocity is a solenoidal vector field, the fluctuation equation isLet $a=\eta k^2>0$, $d=a^2+{\Omega }^2$, and ${\mathbf{C}}_j=i({\mathbf{B}}_0\cdot {\mathbf{k}}_j){\hat{\mathbf{u}}}_j$. For each Fourier mode, the diffusion equation becomes a scalar linear relaxation equation applied to each vector component. Take its long-time periodic response, with $\Omega \ne 0$. For the first mode,while for the second,Thus the complex response coefficients areThey giveAn arbitrary initial fluctuation also contains a homogeneous diffusive transient. In the forced modes this is $\mathrm{Re}{\sum }_j{\mathbf{r}}_j{e}^{i{\mathbf{k}}_j\cdot \mathbf{x}-at}$; more general initial modes also diffuse. The displayed harmonic expression is the periodic particular solution after those transients, not the most general solution at finite time. The constants $a,d$ and the factors of $i$ are unchanged by taking the final real part.

To justify the helicity relation, write ${\hat{\mathbf{u}}}_j={\mathbf{a}}_j+i{\mathbf{c}}_j$ with real vectors. Solenoidality gives ${\mathbf{k}}_j\cdot {\mathbf{a}}_j={\mathbf{k}}_j\cdot {\mathbf{c}}_j=0$. Their cross product is parallel to ${\mathbf{k}}_j$, andThis is the helicity vector of a solenoidal Fourier mode. In these conventions the spatially averaged kinetic helicity density of $\mathrm{Re}({\hat{\mathbf{u}}}_j{e}^{i{\mathbf{k}}_j\cdot \mathbf{x}})$ is $-H_j{k}^2$; the sign is set by the specified cross-product order.

For equal wavevectors, spatial averaging of two real harmonic fields givesTime averaging gives another factor $1/2$, because $⟨{\cos }^2\Omega t{⟩}_t=⟨{\sin }^2\Omega t{⟩}_t=1/2$ and the mixed temporal average is zero. Consequently each diagonal mode contributes to the mean-field electromotive forceDistinct modes with ${\mathbf{k}}_1\ne \pm {\mathbf{k}}_2$ have no spatially averaged cross term. The printed condition also permits the opposite-wavevector case. The possible cross terms still cancel after time averaging: the out-of-phase response coefficients are ${\mathbf{p}}_2=-\Omega {\mathbf{C}}_2/d$ and ${\mathbf{q}}_1=\Omega {\mathbf{C}}_1/d$, while integration by parts gives $⟨{\mathbf{u}}_1\times ({\mathbf{B}}_0\cdot \nabla ){\mathbf{u}}_2{⟩}_x=⟨{\mathbf{u}}_2\times ({\mathbf{B}}_0\cdot \nabla ){\mathbf{u}}_1{⟩}_x$ for the real spatial fields. This handles all the distinct wavevectors allowed by the question, assuming the usual periodic-cell or whole-space spatial average.

Thus the alpha tensor isBoth averaging factors matter. The alpha tensor is real and symmetric, and this contribution vanishes when the mode kinetic helicity densities vanish.

For the two perpendicular wavevectors and common $H_1=H_2=H$, put ${\alpha }_0=\eta H{k}^4/[2({\Omega }^2+{\eta }^2{k}^4)]$. ThenMean-field dynamo action is possible when ${\alpha }_0\ne 0$ and sufficiently long mean-field wavelengths are allowed. The missing $zz$ entry does not prevent an alpha-squared dynamo. Indeed a slowly varying mean magnetic field transverse to $z$, with dependence $e^{iKz+\lambda t}$, obeysThe two growth rates are $-\eta K^2\pm {\alpha }_{0}K$, so the larger is positive if $0<|K|<|{\alpha }_0|/\eta$. This supplies a brief constructive reason; scale separation also requires $|K|\ll k$. A spatially uniform test field itself has zero curl of its mean-field electromotive force and does not grow. If $H=0$, or if boundary conditions exclude all unstable long wavelengths, this particular alpha effect does not yield growth. The conclusion concerns the stated first-order smoothing approximation; no uncomputed turbulent-diffusion correction is assumed.