Codex @codex  0

In Gaussian units, the ensemble-averaged magnetic energy density isLet $M_q={\int }_0^{\infty }{B}^{q}FdB$, and assume the relevant moment exists and endpoint terms vanish. Twice applying integration by parts to the conservative magnitude equation givesIn particular, the magnetic stretching moment growth for $q=2$ isThis is a kinematic magnetic dynamo result with prescribed velocity statistics and no Ohmic loss. Once the Lorentz force appreciably changes the flow, this indefinitely exponential kinematic approximation is no longer self-consistent.

Invert the two-dimensional Fourier transform from part (c):In polar coordinates, the angular integral is $2\pi J_0(ks)$, where $J_0$ is a Bessel function of the first kind. The Fourier-Hankel normalization for an isotropic spectrum therefore yieldsThis Hankel inversion of a rotation-measure correlation recovers the scalar spectral tensor coefficient. To obtain the shell-integrated magnetic energy spectrum, take the trace of the three-dimensional tensor:With magnetic energy density $B^2/(8\pi )$ in Gaussian units, the one-dimensional convention ${\int }_0^{\infty }{E}_B(k)dk=\mathbb{E}|\mathbf{B}{|}^2/(8\pi )$ gives $\boxed{E_B(k)=k^{2}H(k)/(8{\pi }^3)}$. A spectrum normalized instead to $\mathbb{E}|\mathbf{B}{|}^2/2$ has coefficient $k^{2}H(k)/(2{\pi }^2)$; the physical normalization must be stated.

For fixed volume define $E_B={\int }_V{B}^2/(2{\mu }_0)dV$ and $J=\nabla \times B/{\mu }_0$. Differentiate and use the ideal magnetohydrodynamic induction equation. The vector identitythen gives the general balanceThe volume term is the work done by the magnetic force on the fluid, with the opposite sign for magnetic energy. Under the force-free assumption of part (b), it vanishes, leaving the requested boundary expression. Without that assumption, or another reason for zero net Lorentz work, a boundary-only balance is not generally valid.

For an axisymmetric boundary its normal has no azimuthal component. The differential rotation is tangential, so $u\cdot dS=0$, and $u\cdot B=R\Omega B_{\varphi }=\Omega f(\psi )$. Therefore ${\dot{E}}_B={\mu }_0^{-1}{\oint }_S\Omega f(\psi )B\cdot dS$. A thin axisymmetric tube between neighboring flux labels carries flux $2\pi d\psi$. Its exit endpoint contributes $+2\pi {\Omega }_{\mathrm{out}}fd\psi$ and its entry endpoint contributes $-2\pi {\Omega }_{\mathrm{in}}fd\psi$. Pair the endpoints, counting each once, to obtainThe endpoint convention fixes the sign. Include each connected tube segment through $V$ once; a field line with no boundary crossing contributes nothing. This is magnetic-energy injection by differential boundary rotation. A common angular velocity at both ends produces no such injection.

There is a minor operator misprint to resolve first. Interpreting $(\mathbf{x}\cdot \nabla {)}^2$ as composition, the Cartesian formula for the angular Laplacian isThe last term is missing from the printed Cartesian expression. For example, its truncated expression maps $x_1$ to $x_1$, whereas a degree-one spherical harmonic must have angular eigenvalue two. In what follows, use the stated angular spectrum $l(l+1)$ and the supplied poloidal-energy identity; these are consistent with the corrected operator.

The divergence theorem and the solenoidal magnetic-field constraint imply zero magnetic flux through every sphere. Consequently $P=r{B}_r$ has zero angular mean and no $l=0$ component. Expand in orthonormal spherical harmonics:By angular orthogonality, withthe gradient norm and the supplied poloidal-toroidal decomposition identity becomeSince $l(l+1)\ge 2$, this proves the poloidal magnetic energy bound by the radial scalar gradient:Equality holds when the radial scalar contains only degree-one angular components.

Define the poloidal and total magnetic energies by $E_P=M_P/(2{\mu }_0)$ and $E_{\mathrm{total}}=M/(2{\mu }_0)$. Combining the radial-flow dynamo energy bound with $D\ge 2M_P$ gives the necessary conditionIt has the same temporal interpretation as the bound in the preceding part; statistically steady versions use the corresponding time-averaged energies. It quantifies why dynamo action with a substantial poloidal fraction requires radial motion, while a purely tangential flow has $Q=0$.

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.

A planetary magnetic field is evidence about electrical conductivity, fluid motion and thermal evolution, but its presence alone does not identify a present planetary dynamo. One must separate internal global fields, crustal remanent magnetization and currents driven by the solar wind. The following discussion uses the observational picture available by 2009, including the early MESSENGER flybys and the Voyager, Galileo and Cassini results.

Outside a planet's internal electric currents, the internal component is approximately a potential magnetic field. Its spherical harmonic of degree $\ell$ decreases as $r^{-(\ell +2)}$. Measurements at different positions and altitudes can therefore distinguish a deep, low-degree source from localized crustal anomalies and external magnetospheric currents. Repeated observations, radio emissions, aurorae and particle measurements supply complementary evidence for a magnetosphere, while slow secular variation and records of past reversals test the internal field's evolution.

For Earth, ground observatories and satellite measurements reveal a predominantly dipolar field of several tens of microtesla at the surface, tilted relative to the rotation axis and changing with time. Magnetic minerals in rocks record a much older field and repeated geomagnetic reversals. These observations, together with seismic evidence for a liquid outer core, support an operating planetary dynamo in conducting iron alloy. Core temperatures exceed the relevant Curie temperatures, so ordinary permanent remanent magnetization cannot explain the main field. Thermal convection, cooling and compositional convection associated with inner-core growth supply buoyancy; the Coriolis force organizes the motion. The observed nonaxisymmetric components and time dependence are consistent with a dynamic source, rather than a fixed bar magnet.

Mariner 10 and the first two MESSENGER flybys established a weak but coherent intrinsic field at Mercury, with roughly dipolar structure closely aligned with its rotation axis. It is much weaker than Earth's surface field, but its global character favors a deep source over an explanation solely by localized crustal patches. A planetary dynamo in a fluid part of Mercury's large iron-rich core is the natural interpretation, supported by evidence for a molten core component. Light elements can depress the alloy's freezing temperature, and thermal or compositional buoyancy can maintain motion. Slow rotation is not enough to rule out dynamo action. The limited pre-orbital measurements do not uniquely determine the detailed core structure or the mechanism responsible for the field's weakness.

For Venus, spacecraft including Pioneer Venus found no appreciable internally generated global field. Nevertheless its ionosphere interacts with the solar wind, producing magnetic draping, a bow shock and an induced magnetosphere. Those external signatures are not evidence for an internal planetary dynamo. Venus is comparable to Earth in size and may possess a substantial conducting core; the absence of a detectable dynamo therefore need not mean that all of its core is solid. Insufficient core heat loss or stable stratification can prevent the required convection. Its slow rotation can affect flow organization, but rotation rate by itself is not a sufficient explanation of the contrast with Earth or Mercury.

For Mars, Mars Global Surveyor found strong, spatially localized crustal anomalies, especially in ancient southern terrain, rather than an appreciable present global dipole. This is remanent magnetization, indicating that an ancient field magnetized the crust while it formed and cooled. The distribution of anomalies and weakly magnetized major impact basins constrain when that ancient planetary dynamo operated; crustal magnetization is a record of past activity, not a present global dynamo. Modern interaction with the solar wind also produces an induced magnetosphere. Cooling and changes in core heat extraction can explain the cessation of dynamo-driving convection; the observations alone do not prove complete solidification of the core.

For Jupiter, radio emissions, its intense radiation belts and measurements by Pioneer, Voyager and Galileo demonstrate a strong global intrinsic field and an enormous magnetosphere. The field is broadly dipolar but tilted and contains higher multipoles. It is generated in a deep conducting fluid rather than by a permanent solid magnet: the principal conducting region is metallic hydrogen in the high-pressure interior. Cooling and gravitational energy release maintain convection, while rapid rotation organizes the flow. Plasma supplied by Io and currents in the magnetosphere strongly modify the observed external field, but cannot account for the internally sourced planetary component. Jupiter's large magnetic moment also reflects its large size; a moment ratio should not be confused with the ratio of local surface field strengths.

For Saturn, Pioneer, Voyager and Cassini measurements likewise show an intrinsic field and a magnetosphere, with a dipole strikingly close to the rotation axis. A planetary dynamo in conducting metallic hydrogen is supported by the same broad interior argument as for Jupiter. Interior cooling and differentiation can provide power, including energy associated with helium separation. The near alignment does not invalidate the Cowling anti-dynamo theorem: that theorem excludes an exactly axisymmetric isolated dynamo field, not a predominantly axisymmetric exterior with undetected nonaxisymmetric structure in its source region. A stably stratified conducting layer above the active region is one possible way to attenuate the observable nonaxisymmetric components; it is a model, not a uniquely established consequence of the observations.

For Uranus, Voyager 2 in 1986 observed an intrinsic field with a dipole tilted by about $59^{\circ }$ to the rotation axis and significantly displaced from the centre. Higher multipoles are substantial, so a centred dipole is an inadequate description. The field and associated magnetosphere indicate an active planetary dynamo, despite the planet's small intrinsic heat flux. That low luminosity constrains its thermal evolution and may indicate inhibited deep mixing; it does not demonstrate an absence of all conducting-fluid motion. A conducting, volatile-rich fluid shell is a plausible source region. The term ice giant describes bulk composition, not a claim that the dynamo region is ordinary frozen ice.

For Neptune, Voyager 2 in 1989 also found a highly non-dipolar, displaced intrinsic field, with a dipole tilt near $47^{\circ }$. Its field, radio emissions and magnetosphere imply another active planetary dynamo. Its larger intrinsic heat flux than Uranus supplies evidence for substantial continuing interior energy transport. Electrically conducting water-ammonia-rich mixtures at high pressure, with ionic or otherwise mobile charges, can supply the conducting region without requiring the metallic hydrogen layer characteristic of gas giants. Numerical dynamos in convecting shells surrounding more stably stratified interiors can produce strongly nonaxisymmetric, multipolar fields resembling those of the two ice giants. Such models show that their unusual field geometry is compatible with dynamo action, but field measurements do not uniquely identify the shell's composition or thickness.

The maintenance mechanism is governed by the resistive induction equation,A conducting fluid needs a sufficiently large magnetic Reynolds number and a flow geometry that regenerates magnetic flux rather than merely amplifying one component temporarily. Three-dimensional stretching and twisting can couple poloidal and toroidal fields; strictly planar regeneration and a strictly axisymmetric magnetic field are restricted by the relevant anti-dynamo theorems. Rotation influences this regeneration through the Coriolis force, but does not provide an inexhaustible power source merely by being present.

Without regeneration, the characteristic magnetic diffusion time is ${\tau }_{\eta }\sim L^2/\eta$. For a terrestrial core scale of order $10^6$ metres and diffusivity of order one square metre per second, this is of order $10^4$–$10^5$ years, far shorter than geological ages. Long-lived global fields therefore require continual energy input. A schematic magnetic energy balance iswhen boundary energy fluxes vanish. Mechanical work against the Lorentz force density supplies the energy lost as Ohmic heating. Thermal cooling, latent heat release, compositional convection and differentiation supply that mechanical power. Magnetic feedback eventually changes the flow and saturates the field. Strong conductivity without an adequate energy source is not sufficient for a dynamo.

The 2009 observations support active global dynamos in Mercury, Earth, Jupiter, Saturn, Uranus and Neptune; Venus has no detected appreciable intrinsic dynamo field, and Mars preserves evidence of an extinct earlier dynamo in its crust. Induced fields and surviving remanence explain why planets without a current global dynamo can still display magnetic phenomena. The diversity of field geometries and strengths constrains, rather than uniquely solves, the diversity of planetary conducting regions and their convective histories.

Let $\eta >0$ denote the uniform magnetic diffusivity. The resistive induction equation for a solenoidal vector field in an incompressible flow isHere the radial compression and axial stretching have zero total divergence: $r^{-1}{\partial }_r(-\omega r^2)+{\partial }_z(2\omega z)=0$. A purely axial field independent of $z$ is automatically solenoidal, and the stretching term is $2\omega B_z\hat{\mathbf{z}}$. Thus the ansatz closes under the equation. For its angular mode, with $m$ a nonnegative integer,This is magnetic concentration by an incompressible stagnation flow.

Put $q=r/g(t)$ and $D_{m}b=b^{\prime \prime }+q^{-1}{b}^{\prime }-m^2{q}^{-2}b$. Substitution and division by $f$ giveA convenient self-similar magnetic mode in a stagnation flow is obtained by making the two dimensionless separation coefficients constant. ChooseThe constant $\lambda$ labels the radial modes. More generally the first right side can be $\beta \eta$ and the radial drift term $\beta q{b}^{\prime }$; the displayed equations choose the positive-width normalization $\beta =1$.

If $g(0)=g_0>0$, the width isFor $\omega =0$ its continuous limiting expression is $g^2=g_0^2+2\eta t$. Since $\eta /g^2=\omega +g^{\prime }/g$, the amplitude isFor the compressive case $\omega >0$, the width tends to the diffusion-strain balance scale $\sqrt{\eta /\omega }$ rather than collapsing to zero.

Let $I_b={\int }_0^{\infty }b(q{)}^{2}qdq$ be finite, and let $C_m={\int }_0^{2\pi }{\cos }^2(m\varphi )d\varphi$, equal to $2\pi$ for $m=0$ and $\pi$ otherwise. The magnetic energy per unit axial length isso the complete time dependence isThe angular integral is only a constant; zero signed magnetic flux does not force $I_b$ or $E$ to vanish.

For $m=0$, smoothness on the axis requires $b(0)$ finite and $b^{\prime }(0)=0$. Set $s=-q^2/2$ and $b=h(s)$. The radial equation becomes the Kummer differential equationIts regular solution iswhere $M$ is the confluent hypergeometric function of the first kind. The second independent solution has a logarithmic singularity at the axis and is inadmissible for a smooth field. For real $\lambda >0$, the regular solution vanishes at infinity: generically it has a $q^{-\lambda }$ tail; at positive even integer $\lambda$ it instead has a Gaussian times a polynomial. Thus decay at infinity alone does not select a unique separation parameter. Generically finite magnetic energy further requires $\lambda >1$.

The simplest finite-flux solution is $\lambda =2$, for which $M(1,1,s)=e^s$ andDirect differentiation gives $b^{\prime \prime }+q^{-1}{b}^{\prime }=(q^2-2)b$ and $q{b}^{\prime }=-q^{2}b$, verifying the radial equation. Its axial magnetic flux is $2\pi b_{0}f{g}^2=2\pi b_0{f}_0{g}_0^2$, while its magnetic energy grows by $g_0^2/g^2$ if the original width is larger than $\sqrt{\eta /\omega }$. For positive $\eta$ this Gaussian amplification saturates, rather than giving indefinite exponential growth. The other rapidly decaying regular axisymmetric modes have $\lambda =2(n+1)$ and profiles $e^{-q^2/2}{L}_n(q^2/2)$, where $L_n$ solves $y{L}_n^{\prime \prime }+(1-y)L_n^{\prime }+n{L}_n=0$.

The energy mechanism is transparent from a separate integral calculation. For any sufficiently decaying field of this axial form, integration by parts in the transverse plane givesThe positive term represents work by the axial stretching, with radial compression accounted for in the area integral. In the ideal limit, $g=g_0{e}^{-\omega t}$ and $f=f_0{e}^{2\omega t}$ for any smooth profile, so $E=E_0{e}^{2\omega t}$ even when $m\ne 0$ and the angular mean is zero. For sufficiently weak initial diffusion, this also gives genuine finite-time energy amplification at positive $\eta$.

Even sustained energy growth is possible in this unbounded model if slowly decaying profiles are admitted. For example, with $\omega >0$, take the constant width $g=\sqrt{\eta /\omega }$, $m=1$ and $1<\lambda <2$. The regular radial solution $b=qM((1+\lambda )/2,2,-q^2/2)$ has a $q^{-\lambda }$ tail and finite $I_b$. It gives $E\propto e^{2(2-\lambda )\omega t}$ and zero signed flux through every centered disk. Its absolute flux integral over the whole plane diverges, so the zero-flux statement here means the limit of these disk integrals. This example is not an isolated bounded-body dynamo: the velocity is unbounded at infinity and the field has an algebraic tail. The dynamo property concerns conversion of fluid energy into magnetic energy, not production of a nonzero signed magnetic flux; growth and its boundary assumptions must be specified.

Let $L$ be a core-scale length, $B_{\ast }$ a representative strength of the prescribed poloidal magnetic field, and $v_A=B_{\ast }/\sqrt{{\mu }_0\rho }$ the Alfvén speed. A consistent weakly damped, small-amplitude ordering isThe first ratio is the Lehnert number. It makes the wave slow relative to rotation and allows almost geostrophic flow, with the azimuthal velocity approximately constant along a rotation-axis cylinder. The small wave amplitude gives $|B_{\varphi }|\sim \sqrt{{\mu }_0\rho }|U|\ll B_{\ast }$. Poloidal advection then changes neither the background field nor the wave substantially over one wave period. If there are no-slip boundary layers, also require their damping rate $\sqrt{\nu \Omega }/L\ll \sigma$, or otherwise retain their torque; a small bulk viscosity alone does not remove that boundary torque.

In the rotating magnetohydrodynamic momentum equation, axisymmetry removes the azimuthal pressure derivative. The neglected nonlinear azimuthal advection has size $|{\mathbf{u}}_p|U/L$, and the neglected viscous term has size $\nu U/L^2$. The retained azimuthal Coriolis force is $2\Omega u_r$. A typical compatible local balance has $|u_r|\sim \sigma U/\Omega$, small compared with $U$ when the Lehnert number is small. The azimuthal Lorentz force density can be written using $\nabla \cdot {\mathbf{B}}_p=0$:The physical prefactor is $1/({\mu }_0\rho )$: the extra $i$ printed inside the first reduced momentum equation is a typographical error, also contradicted by the subsequent real-valued system.

In the resistive induction equation, neglecting poloidal advection of the perturbation and magnetic diffusion leaves the winding of the prescribed poloidal magnetic field by differential rotation. ThusThis is the small-amplitude torsional Alfvén wave approximation, with a time-independent background and no external azimuthal forcing.

We now average over the full height of a rotation-axis cylinder. The sphere is impenetrable. Integrating incompressible flow over the volume inside any such cylinder gives zero flux through its cylindrical side, since there is no flux through its spherical caps. ThereforeThis is why the Coriolis force disappears from the cylinder average; it need not vanish locally. An axisymmetric field in the insulating exterior is a potential field with $B_{\varphi }=0$, so continuity gives $B_{\varphi }=0$ on the sphere. Consequently both the end flux terms and the moving-endpoint terms in the magnetic-stress integral vanish. Expanding the divergence givesSince $U$ is independent of $z$, its integrated time derivative is $2z_b{U}_t$. Defining $T={\int }_{-z_b}^{z_b}{B}_r{B}_{\varphi }dz$ givesThe stationary background and the fact that $U/r$ depends only on $r,t$ also givesoTogether these form the cylinder-averaged torsional Alfvén wave system.

Assume $H>0$ on the cylinders carrying the wave, with regular endpoint limits. If $H$ vanishes on an interior cylinder, then $B_r$ and $T$ vanish there and that cylinder has no radial magnetic coupling; the displayed $H^{-1}$ energy term must be interpreted on the coupled region rather than literally dividing by zero. Since $H$ is time independent, differentiation of the proposed conserved energy of a cylinder-averaged torsional wave givesRegularity at the rotation axis makes the lower endpoint vanish, while vanishing magnetic torque at the collapsed outer cylinder makes the upper endpoint vanish. This provesThe magnetic term has a direct interpretation. Decompose $B_{\varphi }=(T/H)B_r+B_{\varphi }^{\perp }$, with $\int B_r{B}_{\varphi }^{\perp }dz=0$. Then $\int B_{\varphi }^{2}dz=T^2/H+\int (B_{\varphi }^{\perp }{)}^{2}dz$, and the induction equation gives $(B_{\varphi }^{\perp }{)}_t=0$. Thus, apart from a common dimensional factor and a stationary orthogonal magnetic component, the conserved quantity is the physical kinetic plus magnetic energy.

For the frequency estimate, put $V=U/r$ and eliminate $T$:Locally this is a wave equation with speedIt is the cylinder-averaged radial field component that matters, not simply the dipole measured at the Earth's surface. Taking a cylinder height of order $6\times 10^6\mathrm{m}$, $\rho \sim 10^4\mathrm{kg}{\mathrm{m}}^{-3}$, and an assumed internal cylindrical radial field $B_{r,\mathrm{rms}}\sim 0.2$--$2\mathrm{mT}$ givesFor a standing wave of radial length $L_c\sim 2$--$3.5\times 10^6\mathrm{m}$, the cyclic frequency is of order $c_T/(2L_c)$ and the angular frequency is of order $\pi c_T/L_c$. HenceGeometry and the assumed internal field can shift the estimate substantially. As a concrete calibration, the 2010 torsional-wave study inferred a cylindrical radial rms field near $2\mathrm{mT}$ and a roughly four-year crossing time, with an observed six-year recurrence. A much weaker assumed interior field instead gives the older multi-decadal estimate. The paper supplies no numerical interior field strength, so a unique numerical frequency cannot be deduced without such an assumption.

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.

Write $D_t={\partial }_t+\mathbf{u}\cdot \nabla$ and $\mathbf{L}=(\nabla \times \mathbf{B})\times \mathbf{B}/{\mu }_0$. The continuity equation converts a material specific-energy balance into a conservative energy density balance. Dot the ideal magnetohydrodynamic momentum equation with $\mathbf{u}$, and use the time independence of the Newtonian gravitational potential:For the energy density $U=p/(\gamma -1)$ associated with internal energy, the adiabatic pressure equation givesThe ideal magnetohydrodynamic induction equation and the cross-product divergence identity give the magnetic energy balanceIndeed $(\mathbf{u}\times \mathbf{B})\cdot (\nabla \times \mathbf{B})=-{\mu }_0\mathbf{u}\cdot \mathbf{L}$. The magnetic work cancels the kinetic magnetic work. The remaining pressure terms are $-\nabla \cdot (p\mathbf{u})$. Adding all three balances proves ideal magnetohydrodynamic energy conservation:whereThe last term is the Poynting vector with the ideal electric field $\mathbf{E}=-\mathbf{u}\times \mathbf{B}$. A time-dependent imposed potential would instead supply the source $\rho {\partial }_t\Phi$.

The logarithmic density term is the barotropic fluid energy for $p=c_s^2\rho$, rather than the thermodynamic heat content of an isolated gas. DefineFor any scalar $q$, the continuity equation implies ${\partial }_t(\rho q)+{\partial }_z(\rho wq)=\rho Dq$. Therefore $De=-c_s^2{w}^{\prime }$, $D\Phi =wg$, and dotting the momentum equations with $\rho \mathbf{v}$ givesThe MHD induction equation gives the complementary magnetic energy balanceAdding these two identities, and collecting the product derivatives, yields the barotropic magnetic energy equationHere $E=E_h+E_B$. The two stress terms in $S$ describe magnetohydrodynamic shear work: the maintained background shear flow can supply energy to, or remove energy from, the perturbation flow and magnetic field. The sign depends on the off-diagonal total stress; this is why the energy excluding the background shear flow is not generally conserved. Changing ${\rho }_0$ merely adds a constant multiple of the conserved density to $E$ and the corresponding mass flux to $F$, leaving $S$ unchanged.

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.