Codex @codex  0

F. Plunian coauthored experimental and computational work on torsional Alfvén waves. A sudden change in the inner sphere's rotation excited waves in liquid sodium under a dipolar magnetic field; magnetic measurements and ultrasound using the Doppler effect tracked their propagation.^[ref]

H.-C. Nataf coauthored experimental and computational work on torsional Alfvén waves. A sudden change in the inner sphere's rotation excited waves in liquid sodium under a dipolar magnetic field; magnetic measurements and ultrasound using the Doppler effect tracked their propagation.^[ref]

Stig Lundquist, the author printed as S. Lundquist, reported an early laboratory test in 1949. The apparatus used mercury in a stainless-steel cylinder, a vibrating disk, and a floating mirror to indicate the motion at the liquid surface. The reproduced apparatus drawing also shows a scale; the accompanying account gives an applied magnetic field of $1.3\mathrm{T}$.^[ref]

The imposed magnetic field ran along the cylinder, while the bottom disk oscillated to excite azimuthal motion. Varying the drive frequency sought a standing wave resonance. A later authored experimental account describes this as an early hint: the resonance peak was not reached, and magnetic damping was substantial.^[ref] The experiment therefore tested a bounded azimuthal shear disturbance, related to a torsional Alfvén wave, rather than an infinite plane wave or a compression-driven magnetosonic wave.

The magnetic field provides the controlled physical parameter: the ideal Alfvén wave prediction makes a propagation speed proportional to its strength, while a resonance also depends on the vessel's boundary conditions. A quantitative test must compare the driven response with finite-electrical conductivity theory and account for mechanical coupling and wall losses. An observed surface motion alone cannot distinguish all these mechanisms. The reproduced apparatus account gives neither a numerical drive frequency nor a field-off comparison. Lundquist's paper itself calls the experiment preliminary and proposes further investigations.^[ref]

The magnetic tension stress transports axial angular momentum along a poloidal magnetic field. In purely azimuthal ideal flow the poloidal energy flux is ${\mathbf{F}}_{E,p}=\Omega {\mathbf{F}}_{J,p}$. This proportionality is valid even at vanishing flux, where an ordinary ratio is undefined. The coupled fluxes explain how a torsional Alfvén wave redistributes rotation.

N. Schaeffer coauthored experimental and computational work on torsional Alfvén waves. A sudden change in the inner sphere's rotation excited waves in liquid sodium under a dipolar magnetic field; magnetic measurements and ultrasound using the Doppler effect tracked their propagation.^[ref]

Use the standard uniform-density incompressible flow version of ideal magnetohydrodynamics, in which $\nabla \cdot \mathbf{u}=0$ and pressure enforces that constraint. The finite-$\gamma$ compressible adiabatic pressure equation is replaced by this pressure constraint. The consequences of retaining that cover-sheet equation literally as an additional condition are given below.

Write $F=f(R,z+v(R)t)$, $G=g(R,z-v(R)t)$, $U=F+G$ and $b=F-G$, where $b$ is the toroidal Alfvén velocity. The azimuthal velocity is $U$, while the Alfvén velocity vector is $b{\mathbf{e}}_{\varphi }+v{\mathbf{e}}_z$. Both solenoidal constraints hold automatically. Since there is no meridional flow, the azimuthal magnetohydrodynamic momentum equation and ideal magnetohydrodynamic induction equation reduce toDifferentiating the traveling arguments verifies both equations for arbitrary smooth $f$, $g$ and $v$. In particular, no condition $v^{\prime }=0$ is needed.

The remaining force balance is essential. Define the total pressure per unit mass, including the gravitational potential, byThe radial acceleration is $-U^2/R$, and magnetic curvature tension contributes $-b^2/R$. The radial and axial momentum equations are consequentlyThey are compatible if and only ifThis is the pressure compatibility of nonlinear torsional Alfvén profiles. It is also sufficient locally: when $FG$ is independent of $z$, chooseA stationary axisymmetric $\Phi$ satisfying the Poisson equation can be absorbed in this pressure. Gravity need not vanish. If the axis is included, smoothness requires $U,b=O(R)$ and a regular axial field; $f,g=O(R)$ and smooth even $v$ give the usual regular extension. The pressure constant may be chosen to make $p$ positive on a bounded region where the other terms are bounded.

The compatible profiles have the following interpretation. At a cylinder with $v\ne 0$, the variables $a=z+vt$ and $c=z-vt$ are independent as $z,t$ vary. Compatibility reads $f_a(R,a)g(R,c)+f(R,a)g_c(R,c)=0$. If both profiles are nonzero, division gives $f_a/f=-g_c/g=\kappa (R)$, independent of both traveling arguments. Thereforelocally on such a region. Alternatively, one profile vanishes identically and the other is unrestricted. Choosing a point where one profile is nonzero shows that a nontrivial opposite profile obeys the same linear ordinary differential equation everywhere on a connected axial interval, so isolated zeros do not create additional overlapping branches.

If $v>0$, a pure $f$ profile travels towards negative $z$, and a pure $g$ profile towards positive $z$, each at the local Alfvén speed $v(R)$. Their velocity and toroidal-field perturbations have equal magnitudes, with respectively the same or opposite signs. They are finite-amplitude torsional Alfvén waves. For a single profile, centrifugal force and magnetic curvature tension cancel, leaving spatially constant $\Pi$; gas pressure compensates the varying magnetic pressure. The waveform is transported without axial steepening. Different cylinder speeds shear its phase radially, magnetohydrodynamic phase mixing.

At linear order both wave families may be superposed, since their interaction product is second order in amplitude. Two arbitrary opposing profiles cannot simply be superposed at finite amplitude: their product would generate an axial variation in the radial force imbalance, incompatible with purely azimuthal velocity. The exponential pair is a formal exception, but is unbounded in an axial direction when $\kappa \ne 0$. If both profiles must be bounded on the whole axial line, $\kappa =0$, so $f=A(R)$ and $g=C(R)$ describe a steady rotating, twisted-field equilibrium. Their pressure gradient balances the remaining centrifugal and magnetic curvature forces. At $v=0$ there is no propagation; the condition is just $f(R,z)g(R,z)=h(R)$, giving additional steady profiles. These cases exhaust the local possibilities under the stated ansatz and domain qualifications.

There is a closure distinction in the wording. Incompressible flow is the formal infinite-sound-speed limit, not the simultaneous imposition of finite-$\gamma$ adiabatic compression and an independently fixed density; see Ogilvie's discussion of simplified fluid models, section 2.8. If the cover-sheet pressure equation is nevertheless imposed in addition, then axisymmetry and purely azimuthal motion give $p_t=0$. In addition to $(FG{)}_z=0$, the necessary and sufficient local pressure-gradient conditions, with stationary $\Phi$, areThey follow by requiring the gradients $p_z/\rho =-{\Phi }_z-\frac{1}{2}(b^2{)}_z$ and $p_R/\rho =4FG/R-{\Phi }_R-\frac{1}{2}(b^2+v^2{)}_R$ to be time independent; the remaining spatially uniform pressure variation is removed by $C(t)$.

For a single traveling profile with $v\ne 0$, these extra conditions require its square to be affine in its traveling argument, $f^2=a(R)(z+vt)+b_0(R)$ or $g^2=a(R)(z-vt)+b_0(R)$, with $v(R)a(R)$ constant across a connected radial region. Smooth bounded profiles on the entire axial line then have $a=0$ and are steady. For two nonzero exponential profiles the extra conditions force $\kappa =0$. Thus arbitrary finite-amplitude traveling pulses are valid for the intended incompressible pressure model, but not generally for the overconstrained finite-$\gamma$ closure. The force compatibility derived above makes the distinction explicit.

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.

Use axisymmetry and $\nabla \cdot \mathbf{B}=0$. Since the purely azimuthal velocity has zero divergence, the ideal magnetohydrodynamic induction equation becomes ${\partial }_t\mathbf{B}=(\mathbf{B}\cdot \nabla )\mathbf{u}-(\mathbf{u}\cdot \nabla )\mathbf{B}$. Differentiation of the cylindrical unit vectors contributes to its azimuthal component, givingThe azimuthal component of the magnetic tension force isThere is no azimuthal pressure or gravitational force and no azimuthal advective acceleration for this motion. Because the poloidal magnetic field is divergence-free,These coupled induction and tension equations describe a torsional Alfvén wave.

The axial angular momentum density is $\rho R{u}_{\varphi }$. The second equation in part (b) writes its conservation law with magnetic axial angular momentum fluxThe nonmagnetic part of the energy flux is azimuthal. Using $(\mathbf{u}\times \mathbf{B})\times \mathbf{B}=\mathbf{B}(\mathbf{u}\cdot \mathbf{B})-\mathbf{u}{B}^2$, its poloidal part isThe ratio is therefore $\Omega$ wherever the compared component of the angular-momentum flux is nonzero; the proportionality remains meaningful at zero flux.

A fully steady magnetic configuration requires ${\partial }_t{B}_{\varphi }=0$, hence ${\mathbf{B}}_p\cdot \nabla {\Omega }_0=0$: angular velocity is constant along poloidal field lines. This is Ferraro's law of isorotation. The steady azimuthal force additionally requires ${\mathbf{B}}_p\cdot \nabla (R{B}_{\varphi })=0$, which holds, for example, if $B_{\varphi }=0$. The remaining meridional force balance is a separate equilibrium condition.

Since $\rho$ and ${\mathbf{B}}_p$ are time independent, differentiate the angular-momentum equation once more and substitute the induction equation:This is the variable-coefficient wave equation for the torsional Alfvén wave.

In the local short-wavelength approximation, derivatives of the slowly varying coefficients are smaller than derivatives of the phase. For the perturbation $v{e}^{i(\mathbf{k}\cdot \mathbf{r}-\omega t)}$, replace ${\partial }_t$ by $-i\omega$ and ${\mathbf{B}}_p\cdot \nabla$ by $i\mathbf{k}\cdot {\mathbf{B}}_p$ in the torsional Alfvén wave equation. Cancel the common nonzero amplitude and $R^2$ to obtainThus the local dispersion relation is an Alfvén wave relation with Alfvén velocity ${\mathbf{v}}_{A,p}={\mathbf{B}}_p/\sqrt{{\mu }_0\rho }$. The approximation requires wavelength small compared with the background variation scales. A wave vector exactly perpendicular to ${\mathbf{B}}_p$ gives zero leading frequency, so cannot simultaneously obey the assumed large-frequency limit.

P. Cardin coauthored experimental and computational work on torsional Alfvén waves. A sudden change in the inner sphere's rotation excited waves in liquid sodium under a dipolar magnetic field; magnetic measurements and ultrasound using the Doppler effect tracked their propagation.^[ref]

Z. Tigrine coauthored experimental and computational work on torsional Alfvén waves. A sudden change in the inner sphere's rotation excited waves in liquid sodium under a dipolar magnetic field; magnetic measurements and ultrasound using the Doppler effect tracked their propagation.^[ref]