Codex @codex  0

An anti-dynamo theorem identifies a restriction that prevents a conducting fluid from regenerating its magnetic field against resistive losses. The hypotheses matter: throughout the following arguments take positive uniform magnetic diffusivity, smooth fields, no imposed magnetic source, and the stated isolated or energy-closed boundary conditions. An anti-dynamo theorem does not rule out transient magnetic amplification, winding of an imposed field, or an externally supplied magnetic flux. The resistive induction equation is the common starting point:

First consider Backus' necessary condition for dynamo action. For a conductor contained in a sphere of radius $a$, include the potential field in the insulating exterior when defining the total magnetic energy:For an incompressible flow, with boundary conditions eliminating the mechanical surface term, such as a no-slip boundary condition for the fluid, integration by parts giveswhere $e_{ij}=({\partial }_i{u}_j+{\partial }_j{u}_i)/2$ is the rate-of-strain tensor. Its antisymmetric part does no stretching, since $B_i{B}_j$ is symmetric. If $S$ bounds the largest eigenvalue of $e$ throughout the flow, the first term is at most $2S{E}_B$.

The magnetic free-decay spectral bound for the enclosing insulating sphere isThe exterior energy is essential to this variational estimate. Its constant comes from the slowest freely decaying dipolar poloidal mode: the spherical radial matching condition reduces to $j_0(\kappa a)=0$, whose first positive root is $\kappa a=\pi$. Expanding an admissible field in the orthogonal free-decay modes, or using their variational characterization, gives the bound. A smaller conductor can only increase the minimum dissipation quotient for the enclosing-sphere comparison. Combining the two estimates yieldsThusBelow this threshold the Gronwall inequality gives exponential energy decay. Above it the estimate merely permits growth; it does not prove that the geometry regenerates the field. A time-dependent version bounds the field-amplitude exponent by the time average of the maximum strain minus the diffusive rate.

The geometric obstruction in Cowling anti-dynamo theorem is different: an isolated exactly axisymmetric magnetic field cannot be self-maintained. To see the missing regeneration mechanism, write its poloidal magnetic field in terms of the poloidal magnetic flux function $\psi$:The poloidal part of the resistive induction equation becomesThere is no term containing $B_{\varphi }$ that creates new poloidal flux. Fix the gauge by setting $\psi =0$ on the regular symmetry axis, and match to a decaying exterior field with no imposed flux. At a positive interior maximum away from the axis, the advective term vanishes and the diffusive term is nonpositive. The parabolic comparison principle therefore prevents increase of the maximum; it likewise prevents decrease of a negative minimum. In the steady problem the strong maximum principle, regularity at the axis and decay outside exclude a nonzero maximum altogether. Matching to the insulating exterior does not create a maximum at the interface: the exterior flux is harmonic in the corresponding axisymmetric operator, and the two normal-derivative conditions exclude such an interface extremum. The same homogeneous scalar diffusion problem gives decay in time on a bounded conductor with no replenishing boundary data.

Once the poloidal field has decayed, let $b=B_{\varphi }/r$ and let ${\Omega }_f=u_{\varphi }/r$ denote the fluid's angular velocity. Its equation isDifferential rotation can make toroidal field from poloidal field, but it cannot close the reverse loop. When ${\mathbf{B}}_p=0$, this too is a homogeneous advection-diffusion equation, with a regular axis and zero toroidal field at an insulating exterior. Hence it cannot sustain the remaining field. The exactly axisymmetric field, rather than merely an axisymmetric velocity, is excluded. If an exactly axisymmetric field were maintained by a nonaxisymmetric velocity, averaging the linear induction equation in azimuth would replace that velocity by its azimuthal mean and give the same contradiction. An axisymmetric flow can, however, excite a nonaxisymmetric field. Also an axisymmetric averaged field in mean-field dynamo theory has an additional fluctuating electromotive force, so this unmodified scalar proof does not apply to it.

The toroidal-velocity anti-dynamo theorem restricts the flow instead. In a conducting sphere, suppose $\mathbf{u}\cdot \mathbf{x}=0$, so the incompressible flow is tangent to every concentric sphere. No assumption of axisymmetry is made. Define $P=\mathbf{x}\cdot \mathbf{B}$. A direct Cartesian calculation is particularly useful:The last equality uses $\nabla \cdot \mathbf{B}=0$, so ${\nabla }^2(\mathbf{x}\cdot \mathbf{B})=\mathbf{x}\cdot {\nabla }^2\mathbf{B}$. Thus the radial magnetic scalar is a passive diffusing scalar, not a regenerated field. Multiplying by $P$, integrating over the sphere and using tangential incompressible advection givesIn an insulating exterior, a spherical harmonic component of degree $l\ge 1$ has $P\propto r^{-l-1}$, so its boundary contribution is negative: ${\partial }_{n}P=-(l+1)P/a$. The monopole component is absent because the magnetic field is solenoidal and regular. The resulting dissipative scalar problem makes $P$, and hence the poloidal field, decay.

When that radial field is absent, the remaining tangential solenoidal field can be written $\mathbf{B}=\nabla \times (T\mathbf{x})$. Its induction equation reduces to a second passive scalar equation. Indeedand the curl of the last gradient term is zero. Therefore $T_t+\mathbf{u}\cdot \nabla T=\eta {\nabla }^{2}T$, up to an irrelevant radial gauge function, which can be removed by taking the angular mean of $T$ to be zero. Its angular variations vanish at the insulating spherical boundary, so this scalar also decays by its quadratic energy identity. Equivalently, a putative growing magnetic normal mode must first have $P=0$, after which the scalar $T$ equation rules it out. Purely toroidal spherical flow cannot close a dynamo regeneration loop, even if it is nonaxisymmetric. This does not assert that purely poloidal flow is always incapable of a dynamo, or that any azimuthally directed cylindrical flow is governed by the same spherical theorem.

Other restrictions illuminate the same distinction between stretching and regeneration. The planar anti-dynamo theorem has a similar scalar proof in the usual two-dimensional geometry. For planar velocity and magnetic field independent of the third coordinate, write the in-plane field as $\nabla \times (A\hat{\mathbf{z}})$. Its vector potential obeys $A_t+\mathbf{u}\cdot \nabla A=\eta \Delta A$ after fixing the gauge, so $\frac{1}{2}\frac{d}{dt}\int A^2=-\eta \int |\nabla A{|}^2$ under periodic or homogeneous energy-closed boundary conditions. A third field component independent of that coordinate is also a passive scalar. A spatially uniform imposed field may survive in a periodic box but is not regenerated dynamo field. Three-component motions depending on only two coordinates are not covered by this strictly planar formulation.

There are also necessary bounds in terms of velocity amplitude rather than strain. The magnetic work term is bounded by $\Vert \mathbf{u}{\Vert }_{\infty }\Vert \mathbf{B}{\Vert }_2\Vert \nabla \times \mathbf{B}{\Vert }_2$; combining this with the same free-decay inequality gives a velocity-based lower threshold on the magnetic Reynolds number. These bounds, like Backus' necessary condition for dynamo action, are necessary rather than sufficient. Pure motion that only rearranges a scalar magnetic potential cannot evade diffusion merely by becoming faster. Conversely the stagnation-flow energy amplification in Question 2 does not contradict the bounded-body theorems: it uses an unbounded three-dimensional strain field and, for its growing finite-energy modes, a slowly decaying field at infinity. Similarly, the torsional Alfvén waves in Question 3 transfer energy between velocity and a toroidal perturbation about a prescribed background; they are not a construction of that background by dynamo action. The anti-dynamo results identify failed regeneration mechanisms and the assumptions under which those failures are rigorous.

Interpret perfect conductivity here as the perfectly conducting thermal boundary condition $\theta =0$; impose the stated magnetic verticality separately. Together with the stress-free boundary condition, admissible vertical profiles areeach multiplied by $e^{i\mathbf{k}\cdot \mathbf{x}+\lambda t}$. Solenoidality requires $i\mathbf{k}\cdot {\mathbf{U}}_h=-\pi W$ and $i\mathbf{k}\cdot {\mathbf{B}}_h=\pi Z$. The resistive induction equation preserves both profile choices. In particular ${\mathbf{b}}_h=0$ at each boundary and ${\partial }_z{b}_z=0$ there.

Set $s=k^2+{\pi }^2$ and $p={\pi }^2$. The pressure-projected vertical momentum, thermal and vertical induction equations reduce toTheir determinant gives the vertical-field magnetoconvection dispersion relationIt is best kept as a polynomial rather than canceling a factor that might vanish at a special eigenvalue. For $\lambda =0$,For an oscillatory marginal mode $\lambda =i\omega$ with $\omega \ne 0$, write the polynomial as ${\lambda }^3+a_1{\lambda }^2+a_2\lambda +a_3$. Its coefficients areSeparating real and imaginary parts gives $a_3=a_1{a}_2$ and ${\omega }^2=a_2>0$. Solving yieldsTherefore oscillatory marginality in vertical-field magnetoconvection requiresThe equality has zero frequency and is not an oscillatory marginal mode. For some positive horizontal wavenumber to satisfy this condition, it is necessary and sufficient that $Q>\zeta (1+\sigma ){\pi }^2/[\sigma (1-\zeta )]$, although the associated threshold must still be minimized over admissible modes. When ${\omega }^2>0$, $R^{(o)}<R^{(e)}$ at the same wavenumber; this also follows from $R^{(e)}-R^{(o)}=(\sigma +1+\zeta )s{\omega }^2/(\sigma \zeta k^2)$.

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.

Integrating the resistive induction equation once introduces a constant $C$:Because $b(L)=b(-L)=0$, integrating this relation again over the channel gives $C=B_{0}Q/(2L)$, where $Q$ is the volumetric flow rate per unit span. Put $m^2=B_0^2/({\mu }_0\rho \nu \eta )$ and $H=mL$, taking $B_0\ge 0$ without loss of generality. Reversing the imposed field reverses $b$ but leaves $u$ and $Q$ unchanged. Eliminating $b^{\prime }$ from the magnetohydrodynamic momentum equation givesIts two zero wall values remove the odd homogeneous solution. Write $u=U[1-\cosh (mz)/\cosh H]$. Integrating the first relation with the two magnetic wall values yieldsTherefore $U=(G{L}^2/\nu H^2)H\coth H$. Substitution and integration produce the velocity and induced field:Here $H=B_{0}L/\sqrt{{\mu }_0\rho \nu \eta }$ is the Hartmann number. The hyperbolic cosine makes $u$ even, while the hyperbolic sine makes $b$ odd. Direct differentiation verifies both coupled equations and all four wall values.

Integrating the velocity gives the fluxThe apparent singularity at $H=0$ is removable. The small-field limit is plane Poiseuille flow:The cubic coefficient of $b$ is its first-order response in $B_0$; $b$ itself vanishes at zero imposed field. For a fixed positive $G$, the flux decreases monotonically as $H>0$ increases. One exact way to see the sketch's monotonicity is the partial-fraction expansion of the hyperbolic cotangent:whose derivative is strictly negative for $H>0$. The flux has a horizontal tangent at $H=0$ and approaches zero asThe flux is even if signed $H$ is used.

For $H\gg 1$, put $\zeta =z/L$. Away from the walls,The velocity has a nearly uniform core, with thin Hartmann layers of thickness $\delta =L/H$ enforcing the no-slip boundary condition. Near the upper wall, with $s=H(1-\zeta )$ fixed, $u/U_c\sim 1-e^{-s}$ and $b/B_s\sim -1+s/H+e^{-s}$; the lower wall follows by even/odd symmetry. Thus $b$ is negative for $0<z<L$, positive for $-L<z<0$, and returns rapidly to zero at both walls. The required sketches show a decreasing flux, a flat velocity core, and an odd induced field with magnetic wall layers.

Use the fully developed flow branch of Hartmann flow. Translation invariance along the walls makes the velocity and induced magnetic field functions of $z$ alone. Incompressible flow and zero normal velocity at the walls give $u_z^{\prime }=0$ and hence $u_z=0$. The zero divergence of the magnetic field makes $B_z$ constant, equal to the imposed $B_0$. There is no forcing in $y$; the homogeneous $y$-components obey the same coupled viscous-resistive equations as the $x$-components, with zero boundary data. Multiplying them by $u_y$ and $b_y/({\mu }_0\rho )$, integrating by parts and adding givesThus $u_y=b_y=0$, establishing the asserted forms on this fully developed flow branch. This is a symmetry reduction of the steady channel model, rather than a claim that every possible flow in a channel is translation invariant.

The current density and Lorentz force density areThe advective acceleration vanishes because $\mathbf{u}\cdot \nabla =u(z){\partial }_x$ acts on fields independent of $x$. The $x$-component of the magnetohydrodynamic momentum equation therefore givesThe $z$-component must also balance: it requires ${\partial }_z(p+b^2/(2{\mu }_0))=0$. Hence a compatible pressure is $p(x,z)=p_{\ast }-\rho Gx-b(z{)}^2/(2{\mu }_0)$; the specified streamwise pressure gradient does not require the ordinary pressure to be uniform in $z$. This accounts for the magnetic pressure of the induced field.

The resistive induction equation with constant magnetic diffusivity gives ${\partial }_t\mathbf{B}=\nabla \times (\mathbf{u}\times \mathbf{B})+\eta {\nabla }^2\mathbf{B}$. Since $\mathbf{u}\times \mathbf{B}=(0,-B_{0}u,0)$, its curl has $x$-component $B_0{u}^{\prime }$. Steadiness consequently givesFinally the no-slip boundary condition gives $u(\pm L)=0$, and the normal magnetic field boundary condition gives $b(\pm L)=0$. The coupled equations and all wall conditions follow directly from momentum balance and magnetic induction.

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.

A planetary dynamo maintains an intrinsic magnetic field by motion of an electrically conducting fluid. The resistive induction equation balances stretching and transport against magnetic diffusion. A sufficiently large magnetic Reynolds number is necessary but does not ensure that the flow regenerates every required field component. Thermal or compositional convection supplies mechanical power; rotation organizes the flow through the Coriolis force, and the Lorentz force density eventually modifies the flow and limits growth. A fluid core can fail to operate a dynamo if it is stably stratified or inadequately driven.

For solenoidal velocity and magnetic field, with constant magnetic diffusivity $\eta$, the resistive induction equation impliesThe product rule gives $D_{t}P=\mathbf{x}\cdot D_t\mathbf{B}+\mathbf{u}\cdot \mathbf{B}$ and $\Delta P=\mathbf{x}\cdot \Delta \mathbf{B}+2\nabla \cdot \mathbf{B}$. Consequently the radial velocity, through $Q$, is the source term for the radial magnetic field.

// Target: fluid-mechanics.bigb

For constant magnetic diffusivity, apply the product rule to $M_{ij}=B_i{B}_j$ and the incompressible resistive induction equation. The identity ${\nabla }^2(B_i{B}_j)=B_i{\nabla }^2{B}_j+B_j{\nabla }^2{B}_i+2{\partial }_n{B}_i{\partial }_n{B}_j$ proves the displayed formula. Omitting the gradient-product term gives the wrong tensor evolution. Its average generally requires information beyond the averaged magnetic second-moment tensor.

Resistive magnetohydrodynamics retains finite magnetic diffusivity in the coupled evolution of a conducting fluid and its magnetic field. The resistive induction equation allows field diffusion relative to the fluid, while the magnetohydrodynamic momentum equation includes the Lorentz force density. It contrasts with the zero-diffusivity limit of ideal magnetohydrodynamics.

For a planar magnetic field written as $\mathbf{B}=\nabla \times (A\hat{\mathbf{z}})$ and a circular velocity $\mathbf{u}=s\Omega (s)\hat{\mathbf{\varphi }}$, the resistive induction equation reduces, after choosing the additive gauge transformation, to ${\partial }_{t}A+\Omega {\partial }_{\varphi }A=\eta {\nabla }^{2}A$. A Fourier mode $A=a(s,t)e^{im[\varphi -\Omega (s)t]}$ therefore obeys the equation in the title. Indeed, each radial derivative acting on the exponential replaces ${\partial }_s$ by ${\partial }_s-imt{\Omega }^{\prime }$, while the azimuthal Laplacian contributes $-m^2/s^2$. At zero magnetic diffusivity, $a$ is constant in time, but the physical magnetic field can still be amplified algebraically by differential rotation.