Codex @codex  0

Let $h_c=\mathbf{u}\cdot \mathbf{B}$ be the cross-helicity density. The ideal magnetohydrodynamic induction equation implies $D\mathbf{B}/Dt=(\mathbf{B}\cdot \nabla )\mathbf{u}-\mathbf{B}\nabla \cdot \mathbf{u}$. Dotting the momentum equation with $\mathbf{B}$ eliminates the Lorentz force density, since it is perpendicular to the magnetic field. The material derivative of $h_c$ is thereforeUsing $dh=Tds+dp/\rho$ for specific enthalpy $h$ and specific entropy $s$, together with Gauss's law for magnetism, converts this toThe source vanishes when magnetic field lines lie in constant-specific entropy surfaces, in particular in a homentropic flow. Integrating over a volume conserves its total cross-helicity only if the boundary flux also vanishes. This local law is derived in Ogilvie's astrophysical fluid dynamics notes.

A cylindrical magnetostatic equilibrium is an axisymmetric vector field and pressure configuration independent of the axial coordinate. Its radial Lorentz force density balances the pressure gradient, as expressed by cylindrical magnetostatic pressure balance.

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.

The Lorentz force per unit area on an ideal current sheet is computed using the mean of its limiting magnetic fields. The sheet's own local magnetic field is equal and opposite on its two sides and cancels in this mean. Equivalently, integrate the Lorentz force density through a thin regularized layer: its self-field contribution averages to zero, leaving the external field acting on the surface current density.

Hartmann flow is a fully developed viscous conducting-fluid flow driven along a planar channel while a uniform magnetic field crosses the walls. Induction and the Lorentz force density couple its velocity and tangential induced field. For no-slip, normal-field walls separated by $2L$, the velocity is a hyperbolic cosine profile that approaches plane Poiseuille flow at zero imposed field and develops a uniform core with Hartmann layers at large Hartmann number.

The inviscid magnetohydrodynamic momentum equation balances inertia, pressure, gravitational acceleration and Lorentz force density. The magnetic force separates into magnetic tension and magnetic pressure. It is coupled to the continuity equation and the ideal magnetohydrodynamic induction equation.

In aligned steady ideal flow with $p=c_s^2\rho$ and constant isothermal sound speed, the Lorentz force density has no field-parallel component. Projecting momentum along the magnetic field makes $u^2/2+\Phi +c_s^2\ln (\rho /{\rho }_0)$ constant on each connected magnetic surface. The arbitrary reference density changes only the additive constant.

With the momentum-flux sign convention used for the Maxwell stress tensor, Maxwell's equations implywhere $\mathbf{g}$ is the electromagnetic momentum density and ${\mathbf{f}}_{\mathrm{L}}$ is the Lorentz force density exerted by the field on matter. Thus total field-plus-matter momentum is locally conserved.

In magnetically channelled accretion, a strong magnetic field guides an accretion flow along flux tubes. The component of the Lorentz force density parallel to the magnetic field vanishes, so the longitudinal flow can be treated hydrodynamically while the field controls the transverse geometry. A magnetic dipole field near a neutron star channels matter toward small polar caps.

Magnetic relaxation converts magnetic energy into material motion through the Lorentz force density, and dissipates the resulting kinetic energy through viscosity. In a perfect conductor, magnetic flux freezing preserves field-line connectivity throughout every smooth stage. A regular limiting state has zero velocity and is a magnetostatic equilibrium, but its magnetic field can contain current sheets. Conservation of one global magnetic helicity is weaker than preservation of the complete frozen field-line structure.

For constant mass density $\rho$, kinematic viscosity $\nu$ and magnetic permeability ${\mu }_0$, the incompressible momentum equation isIt is the Navier-Stokes equation with Lorentz force density. The magnetic force decomposes into magnetic tension and a magnetic pressure gradient. A uniform streamwise pressure gradient need not make ordinary pressure independent of transverse coordinates.

The leading force balance in a rapidly rotating conducting fluid when inertial and viscous forces are small, while Coriolis force, pressure, buoyancy and Lorentz force density remain significant. The Rossby number and Ekman number measure the omitted inertia and bulk viscosity. Rotation does not balance an arbitrary magnetic torque: impenetrable axisymmetric boundaries impose the Taylor constraint. Magnetic diffusion can remain in the induction equation even when viscosity is neglected in momentum balance.

Use the real magnetic field $\mathcal{B}$ in the coupled equations; its leading periodic approximation will be $\mathrm{Re}(\mathbf{B}{e}^{i\omega t})$. In the conducting incompressible flow, the magnetohydrodynamic momentum equation, resistive induction equation and solenoidal magnetic-field constraint areHere $p$ is pressure, $\mathbf{j}$ is electric current density, and $\eta$ and $\nu$ are magnetic diffusivity and kinematic viscosity. The moving-conductor Ohm law gives $\mathbf{E}=-\mathbf{u}\times \mathcal{B}+\eta \nabla \times \mathcal{B}$.

In the insulating exterior, neglecting displacement current, the magnetic field is curl-free away from the wire and solenoidal. The imposed electric current density suppliesThe exterior electric field also satisfies Faraday's law. At the rigid wall, impose the no-slip boundary condition $\mathbf{u}=0$. At finite conductivity, the interface conditions for electromagnetic fields give continuous normal magnetic field, continuous tangential magnetic field for equal permeability and no prescribed singular surface current, and continuous tangential electric field. The normal electric current density is zero at an insulating wall; it is automatic for the present two-dimensional fields with current only along $z$. The induced fields decay far from the wire and wall, and initial conditions, or selection of the periodic state after transients, complete the specification. In the limiting perfect conductor approximation a surface current can emerge, permitting a tangential-field jump.

The fully coupled solution need not remain monochromatic: a harmonic magnetic field produces both a mean and a twice-frequency Lorentz force density, and the resulting velocity can generate further harmonics. The following parts consistently use the leading magnetic response with fluid motion neglected.

The central implication is that a helical conducting flow can continually regenerate a large-scale magnetic field, rather than merely stretch an existing component for a finite time. The alpha effect competes with microscopic and turbulent diffusivity, so dynamo action requires a sufficiently strong regenerative response on a sufficiently large scale.

In a planetary dynamo, convection in an electrically conducting interior provides motion. Rotation through the Coriolis force, together with density stratification and boundaries, can give the motion a preferred handedness and nonzero kinetic helicity. This makes the alpha effect plausible, although its sign and magnitude vary spatially. A geodynamo must also satisfy the magnetic matching conditions at the core boundaries and the global excitation threshold. Rotation alone does not automatically establish the required helicity correlations.

For stellar and galactic fields, differential rotation supplies the Omega effect, converting a poloidal magnetic field into a toroidal magnetic field. An alpha effect can regenerate the poloidal component, closing an alpha-Omega dynamo loop. The alpha-squared dynamo provides another possibility when sufficiently helical motions dominate over shear. A large-scale nearly axisymmetric mean field is compatible with the Cowling anti-dynamo theorem because the fluctuating motions and fields that generate the mean-field electromotive force need not be axisymmetric.

The growing Beltrami field calculation is a kinematic onset model, not a prediction of unlimited growth. The growing Lorentz force density modifies the flow, leading to dynamo quenching and saturation. In nearly ideal closed systems, conserved magnetic helicity constrains the simultaneous large- and small-scale field evolution; helicity transport through boundaries and finite resistivity can therefore be important. Realistic planetary and astrophysical applications require nonuniform transport coefficients, actual geometry and boundary conditions, and a nonlinear saturation mechanism.

Let $K$ be the kinetic energy, $M$ the magnetic energy and $E=K+M$:Dot the magnetohydrodynamic momentum equation with $\mathbf{v}$ and integrate. Incompressibility and the zero boundary velocity remove the advective and pressure work. The viscous term integrates to $-\rho \nu {\int }_D|\nabla \mathbf{v}{|}^2$. ThusThe ideal magnetohydrodynamic induction equation gives, with no boundary contribution,The Lorentz force density therefore transfers energy between the field and the fluid without creating total energy. Since $\nabla \cdot \mathbf{v}=0$ and $\mathbf{v}=0$ on the boundary, another integration by parts gives $\int |\nabla \mathbf{v}{|}^2=\int |{\mathbf{\omega }}_v{|}^2$, where ${\mathbf{\omega }}_v=\nabla \times \mathbf{v}$ is vorticity. The viscous magnetic-relaxation energy identity isThe conserved nonzero magnetic helicity prevents all the field energy from disappearing. Fix a bounded inverse-curl convention for the magnetic vector potential, so that $\Vert \mathbf{A}{\Vert }_2\le C_D\Vert \mathbf{B}{\Vert }_2$. The Cauchy-Schwarz inequality gives the magnetic-energy lower bound from helicity,Consequently $E$ decreases to a finite positive limit andUnder the intended regular long-time interpretation of smooth vorticity, the dissipation cannot keep producing peaks of fixed height and arbitrarily short duration. More precisely, if $d(t)=\Vert {\mathbf{\omega }}_v(t){\Vert }_2^2$ is uniformly continuous, uniformly continuous integrable dissipation tends to zero proves $d(t)\to 0$. The Poincaré inequality and the identity above then imply $\Vert \mathbf{v}(t){\Vert }_2\to 0$.

To identify the asymptotic state, suppose the regular relaxation has a limit, or work on a compact invariant limiting set on which the energy is continuous. On such a set the energy is constant, so the displayed identity forces ${\mathbf{\omega }}_v=0$, and the zero boundary velocity forces $\mathbf{v}=0$. The induction equation is then stationary, and the momentum equation becomesThis is magnetostatic equilibrium, with nonzero field allowed and required here by conserved magnetic helicity. It is not necessarily a force-free magnetic field: a pressure gradient can balance the magnetic force. If the field develops a current sheet, equilibrium is understood with the corresponding interface force balance.

There is a mathematical qualification to the smoothness premise. Smoothness at each finite time alone does not prove uniform long-time bounds, uniform continuity of dissipation or convergence to one limiting field. The energy calculation proves finite total dissipation unconditionally for a smooth solution; the settling conclusion is the regular-relaxation argument just given. A rigorous global convergence theorem requires that additional long-time control. In particular, smooth vorticity does not require smooth electric current density in the limiting state.

Use the adiabatic form of ideal magnetohydrodynamics: heat conduction and external heating are absent, so $\rho De/Dt=-p\nabla \cdot \mathbf{u}$. Infinite electrical conductivity and zero viscosity alone would not remove the thermal source terms on the cover page. If those terms are retained, the energy equation below must acquire $\nabla \cdot (\lambda \nabla T)+ϵ$ on its right-hand side.

Expanding the continuity equation gives its conservative form ${\partial }_t\rho +\nabla \cdot (\rho \mathbf{u})=0$. Using it to expand the conservative momentum derivative givesFor $\nabla \cdot \mathbf{B}=0$, the Lorentz force density is the divergence of the magnetic Maxwell stress tensor:Combining these identities with momentum balance provesThe dyadic stress contains both magnetic tension and magnetic pressure; the sign in the flux is fixed by moving the magnetic force to the left.

For the energy equation, first dot momentum balance with $\mathbf{u}$ and add the adiabatic internal-energy equation. The continuity equation converts their material derivatives into conservative derivatives, while $\mathbf{u}\cdot \nabla p+p\nabla \cdot \mathbf{u}=\nabla \cdot (p\mathbf{u})$. The result isThe ideal magnetohydrodynamic induction equation gives the magnetic counterpart. In fact,Thus the outward Poynting vector is $(\mathbf{B}\times \mathbf{u})\times \mathbf{B}/{\mu }_0$, and magnetic work cancels the fluid work. Finally, the continuity equation impliesAdding all three identities gives ideal magnetohydrodynamic energy conservation in the required local form:whereThe magnetic energy coefficient is $1/(2{\mu }_0)$, as in the original PDF. For a self-consistent, time-dependent gravitational field this local balance should not be confused with a global self-gravitating energy integral: the latter uses gravitational field energy, or the familiar one-half factor in $\int \rho \Phi$ after the field is eliminated.

Work in the shock frame of a thin planar discontinuity, with normal speed $u>0$. A smooth gravitational potential has no finite jump through a vanishingly thin front; its advected contribution cancels using the common mass flux. Integrating the conservation laws gives constant normal mass, momentum and energy fluxes, together with the normal-field and tangential-electric-field conditions of the ideal magnetohydrodynamic shock conditions.

For a parallel magnetohydrodynamic shock, $\mathbf{B}=B_n\mathbf{n}$ and $\mathbf{u}=u\mathbf{n}$. The magnetic solenoidal constraint makes $B_n$ the same on both sides. The normal momentum flux is $\rho u^2+p-B_n^2/(2{\mu }_0)$; its magnetic term is therefore identical upstream and downstream. The Poynting vector is zero because $\mathbf{u}\times \mathbf{B}=0$. Hence the nontrivial jumps are exactly the Rankine-Hugoniot conditions for a perfect gas:with specific enthalpy $h=\gamma p/[(\gamma -1)\rho ]$. Physically, motion along the straight field neither bends it nor compresses its transverse flux. The field consequently provides no changed stress or energy flux in this purely parallel geometry. This conclusion does not apply to a shock that develops a transverse magnetic field.

For a perpendicular magnetohydrodynamic shock, the magnetic field is tangential and the velocity normal. Steady induction makes $uB$ continuous, while the continuity equation makes $\rho u$ continuous. Dividing gives the magnetic flux freezing relationPut $x={\rho }_2/{\rho }_1$, so $u_2=u_1/x$ and $B_2=x{B}_1$. Constant normal momentum flux now reads $\rho u^2+p+B^2/(2{\mu }_0)$. Substitution givesThe magnetic correction is negative for compression: part of the incoming ram pressure builds magnetic pressure instead of thermal pressure.

Here the Poynting vector has normal component $u{B}^2/{\mu }_0$. Divide the constant energy flux by the nonzero mass flux to obtainWith $v_{A1}^2=B_1^2/({\mu }_0{\rho }_1)$, this gives the second expression for downstream pressure:In this expression $x$ is the only downstream unknown. Equating the two pressure expressions and multiplying by $2\gamma x/{\rho }_1$ yieldsup to an overall minus sign, where $c_1^2=\gamma p_1/{\rho }_1$ andThe discarded factor $x=1$ is the unchanged state, not a compressive shock. This derives the polynomial for the compression ratio of a perpendicular magnetohydrodynamic shock without omitting its trivial branch.

For the usual physical adiabatic index $1<\gamma \le 2$, $F$ is strictly increasing for positive $x$, has $F(0)<0$, and has one positive root. Sincethat root is greater than one precisely when $u_1^2>c_1^2+v_{A1}^2$. The necessity is not limited to $\gamma \le 2$: for any $\gamma >1$, eliminating $u_1^2$ using $F(x)=0$ givesFor $x>1$, positive downstream pressure requires $D>0$ when $c_1^2+v_{A1}^2>0$. Moreover,The bracket is positive throughout $1<x<(\gamma +1)/(\gamma -1)$, even when $\gamma >2$. The exactly cold unmagnetized limit gives the familiar strong-shock upper compression ratio directly. Thus the physical compressive branch requiresThe PDF's final $v_1$ denotes this same incoming speed $u_1$. The right side is the squared speed of a perpendicular fast magnetosonic wave: the upstream flow must be super-fast-magnetosonic so that this compressive signal cannot run upstream against the incoming flow. At equality the shock amplitude vanishes into a characteristic disturbance.

Take $u<0$ for inward radial velocity, and assume a constant specific-heat ratio $\gamma >1$. The isentropic flow has $p=K{\rho }^{\gamma }$ and adiabatic sound speed $c^2=\gamma K{\rho }^{\gamma -1}$. Steady spherical mass conservation and the radial Euler equation for fluid motion giveEliminating the mass density derivative gives $(u-c^2/u)u^{\prime }=2c^2/r-GM/r^2$, orDifferentiating $c^2\propto {\rho }^{\gamma -1}$ gives the corresponding sound speed equation:At a smooth sonic point, $u_s^2=c_s^2$. A finite derivative requires the numerator of the flow equation to vanish too. ThusThis regularity requirement concerns a smooth transonic branch, rather than an arbitrary singular solution reaching unit Mach number.

Integrating the momentum equation using the specific enthalpy $c^2/(\gamma -1)$ yields the Bernoulli equationFor a finite positive mass accretion rate, its large-radius solution hasThese expressions satisfy both the flow equations asymptotically and show $u\to 0$ and $\rho \to {\rho }_{\infty }$. In particular $u^{\prime }/u\sim -2/r$ agrees with the first boxed equation. Combining the Bernoulli equation with sonic regularity additionally gives $c_s^2=2c_{\infty }^2/(5-3\gamma )$, so a finite positive sonic point occurs in the range identified below.

For the small-radius inward free fall branch, neglecting enthalpy against $GM/r$ in the Bernoulli equation predicts $u^2\sim 2GM/r$. The continuity equation then fixes $\rho \propto r^{-3/2}$, and consequently $c^2\propto r^{-3(\gamma -1)/2}$. The neglected-to-leading ratio isIt tends to zero precisely whenThis makes the assumed supersonic flow and leading free fall self-consistent. At equality, the ratio is of constant order rather than vanishing; above it, the proposed asymptotic balance fails. The isothermal limit $\gamma =1$ also permits free fall, with logarithmic rather than power-law enthalpy.

For the stellar funnel, write the magnetic scalar potential as $\psi =\mu \cos \theta /r^2$, where $R=r\sin \theta$ and $z=r\cos \theta$. Its magnetic dipole field has $B_r=-2\mu \cos \theta /r^3$ and $B_{\theta }=-\mu \sin \theta /r^3$, so a magnetic-field-line equation isThe edge of the polar flux tube meets the stellar surface at ${\theta }_{\ast }\simeq a/R_{\ast }$. Thus, near the pole, ${\theta }^2\simeq {\theta }_{\ast }^{2}r/R_{\ast }$, and its transverse radius is $R\simeq r\theta$. This proves the dipolar flux-tube area lawEquivalently magnetic flux conservation gives $A|\mathbf{B}|=\mathrm{constant}$ and $|\mathbf{B}|\propto r^{-3}$. The stipulated small cross section keeps the tube close to the polar axis, where its direction is approximately radial; the approximation is not valid all the way to arbitrarily large $r$.

Outside the star the prescribed field is a gradient, hence $\nabla \times \mathbf{B}=0$ and $\mathbf{j}=0$. Its Lorentz force density is therefore $\mathbf{j}\times \mathbf{B}=0$. The magnetic-pressure and magnetic-tension parts cancel in this exact current-free dipole. Field-aligned velocity also has $\mathbf{u}\times \mathbf{B}=0$, consistent with a steady prescribed field. The following calculation uses the narrow polar funnel approximation.

Now steady mass conservation reads $\rho uA=\mathrm{constant}$. With $A\propto r^3$, replace the spherical factor 2 by 3:The Bernoulli equation is unchanged because the magnetic field supplies no work or force. On a gravity-dominated free fall branch, $u^2\sim 2GM/r$, but $\rho \propto r^{-5/2}$ and $c^2\propto r^{-5(\gamma -1)/2}$. HenceThe same bound follows from regular transonic accretion in a power-law tube: $u_s^2=c_s^2=GM/(3r_s)$ and the reservoir-matched Bernoulli equation gives $c_s^2=2c_{\infty }^2/(7-5\gamma )$. The stellar surface cuts off the mathematical small-$r$ limit; the criterion refers to a funnel that becomes gravity dominated and crosses its sonic point before reaching that surface.

Write $m_0={\mu }_0\rho$ and treat the imposed mean magnetic field as uniform on the forcing scale. In first-order smoothing, retain terms linear in the fluctuating fields in their equations, while keeping their quadratic correlation when calculating the mean-field electromotive force. The inertial fluctuation, nonlinear induction term and quadratic fluctuating Lorentz force density are neglected in the steady response. The linear magnetic force isAbsorb the gradient into the total pressure. The resulting equations arewith $\nabla \cdot \mathbf{u}=\nabla \cdot \mathbf{b}=0$. Taking the divergence of the first equation gives ${\nabla }^{2}P=0$ for solenoidal forcing. In a periodic or homogeneous Fourier calculation its nonconstant part therefore vanishes.

For each Fourier mode with wavevector $\mathbf{q}$, $|\mathbf{q}|=k$, define $\ell =\mathbf{B}\cdot \mathbf{q}$ andSolving the two coupled linear equations gives the monochromatic coupled magnetic and velocity responseThe magnetic coupling suppresses modes with large projection of their wavevector along the imposed field.

To calculate the mean-field electromotive force, group each pair of opposite wavevectors into a real mode and average its spatial phase. Put $\mathbf{n}=\mathbf{q}/k$, and define its forcing correlationsFor any solenoidal real mode $\mathbf{v}=\mathbf{a}\cos (\mathbf{q}\cdot \mathbf{x})+\mathbf{c}\sin (\mathbf{q}\cdot \mathbf{x})$, its two amplitude vectors are perpendicular to $\mathbf{q}$. Direct multiplication givesAlso, integration by parts gives $⟨(L\mathbf{g})\times (L\mathbf{f})⟩=-{\ell }^2⟨\mathbf{f}\times \mathbf{g}⟩$, where $L=\mathbf{B}\cdot \nabla$. Expanding the response products therefore yieldsSumming over independent real modes, or integrating over their angular spectrum, is the required general answer for $⟨\mathbf{u}\times \mathbf{b}⟩$.

At zero field, $\mathcal{E}(0)={\sum }_{\mathbf{q}}{\mathbf{C}}_{\mathbf{q}}/(\nu \eta k^4)$. Thus the stated zero-field condition removes the unweighted sum of the cross correlations. For anisotropic correlated forcing it does not necessarily remove their differently weighted sum at nonzero field. The displayed general formula retains that allowed contribution. It disappears if cross forcing is absent mode by mode, or under the jointly isotropic statistics assumed for the final part.

The helicity content can also be expressed directly using the actual fluctuating fields. Define their kinetic helicity $H_u(\mathbf{q})=⟨{\mathbf{u}}_{\mathbf{q}}\cdot \nabla \times {\mathbf{u}}_{\mathbf{q}}⟩$ and current helicity $H_b(\mathbf{q})=⟨{\mathbf{b}}_{\mathbf{q}}\cdot \nabla \times {\mathbf{b}}_{\mathbf{q}}⟩$. Substitute $\mathbf{f}=\nu k^2\mathbf{u}-L\mathbf{b}/m_0$ and $\mathbf{g}=\eta k^2\mathbf{b}-L\mathbf{u}$ into their cross correlation. This gives the exact identityIt shows the competing kinetic and magnetic contributions without introducing a spurious singularity when the factor on the left vanishes. For a monochromatic field, ${\mathbf{a}}_b=\nabla \times \mathbf{b}/k^2$ is a magnetic vector potential, so $H_b=k^2⟨{\mathbf{a}}_b\cdot \mathbf{b}⟩$: current helicity is $k^2$ times magnetic helicity. In the isotropic weak-field limit, with zero-field fluctuations ${\mathbf{u}}_0,{\mathbf{b}}_0$, the scalar alpha effect isThis identifies both signs and their different viscous and magnetic relaxation times.

For jointly isotropic forcing let $h_f,h_g$ now denote the total forcing helicities, and set $S_h=\nu h_g/m_0-\eta h_f$. The cross-forcing contribution is even in $\mathbf{B}$, whereas isotropy requires any vector depending only on $\mathbf{B}$ to be parallel to it and odd under reversal. Hence that contribution vanishes. A natural symmetric alpha tensor isOnly its longitudinal coefficient contributes to $\mathcal{E}={\alpha }_{\parallel }\mathbf{B}$. Introduce $h=|B|/[k\sqrt{m_0\nu \eta }]$ and $A_0=\nu \eta k^4$. Angular integration givesTherefore the strong-field quenching of an isotropic electromotive force isprovided $S_h\ne 0$. If the two helicity contributions cancel exactly, this electromotive force vanishes.

There is a tensor qualification to the printed scaling claim. The natural transverse coefficient instead contains $\frac{1}{2}{\int }_0^1(1-{\mu }^2)(1+h^2{\mu }^2{)}^{-2}d\mu \sim \pi /(8h)$ and scales as $|B{|}^{-1}$. The relation ${\mathcal{E}}_i={\alpha }_{ij}{B}_j$ by itself leaves transverse components undetermined. One may choose the effective representative ${\alpha }_{ij}={\alpha }_{\parallel }{\delta }_{ij}$, for which the requested $|B{|}^{-3}$ scaling holds. It is the observable longitudinal response that is fixed, rather than every component of a unique field-dependent tensor.

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 $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.