Codex @codex  0

In a smooth interior region with locally constant shear ${\Omega }^{\prime }$, a wound azimuthal Fourier mode has radial wavenumber $k_s\simeq -m{\Omega }^{\prime }t$. Integrating its resistive damping rate $\eta k_s^2$ gives the exponential decay in the title. The characteristic time is ${\tau }_c=[3/(\eta m^2{\Omega }^{\prime 2}){]}^{1/3}$. For $m=1$ and ${\Omega }^{\prime }=-{\Omega }_0/d$, ${\Omega }_0{\tau }_c=(3\mathrm{Rm}{)}^{1/3}$, where $\mathrm{Rm}={\Omega }_0{d}^2/\eta$. This is a high-magnetic Reynolds number interior transient on a time long compared with the rotation time but short compared with $d^2/\eta$. Physical magnetic field components may have additional algebraic factors from taking derivatives of the wound magnetic vector potential. A maintained exterior field, an axial core and a circulation-edge boundary layer prevent interpreting this estimate as an exact, uniform, infinite-time solution.

For zero mean velocity and magnetic field in a shearing box, the perturbation energy obeys $\dot{E}\le 2(|A|-4{\pi }^2\eta /L^2)E$ when kinematic viscosity and magnetic diffusivity agree. Bound the work of Reynolds stress and Maxwell stress tensor by $2|A|E$, and use the spectral gap of a shearing-periodic box to bound dissipation below. Thus the displayed bound on magnetic Reynolds number guarantees exponential decay and excludes sustained turbulence.

Flux expulsion is redistribution of magnetic flux out of the interiors of persistent circulating flows into their edges and surrounding regions. In two-dimensional incompressible flow, a magnetic flux function obeys ${\partial }_{t}A+\mathbf{u}\cdot \nabla A=\eta {\nabla }^{2}A$. At large but finite magnetic Reynolds number, circulation stretches gradients of $A$, enabling resistive smoothing and weak interior magnetic field with strong field near the circulation boundary. Finite magnetic diffusivity is essential to the eventual rearrangement; exact ideal magnetic flux freezing alone does not justify a flux-free steady cell. Magnetic back reaction can limit the process.

For a locally uniform mean magnetic field, first-order smoothing solves ${\partial }_t\mathbf{b}-\eta \Delta \mathbf{b}=(\overline{\mathbf{B}}\cdot \nabla )\mathbf{u}$. Inserting its causal diffusion propagator into the mean-field electromotive force and taking the trace of the isotropic response tensor gives the displayed relation to kinetic helicity. When diffusion during a correlation time is negligible, it reduces to $\alpha =-{\tau }_H⟨\mathbf{u}\cdot \nabla \times \mathbf{u}⟩/3$. At small magnetic Reynolds number, a quasistatic monochromatic velocity instead gives $\alpha =-⟨\mathbf{u}\cdot \nabla \times \mathbf{u}⟩/(3\eta k_0^2)$. The effective correlation time and diffusion weighting should not be silently discarded.

Finite electrical conductivity permits resistive smoothing of a magnetic field. With uniform coefficients and no motion the field satisfies a diffusion equation; with motion the diffusion term supplements the curl of $\mathbf{u}\times \mathbf{B}$. The magnetic Reynolds number compares these effects.

The magnetic Prandtl number compares kinematic viscosity with magnetic diffusivity. It is the ratio of the viscous and magnetic diffusion coefficients; equivalently it is the magnetic Reynolds number divided by the Reynolds number when the same speed and length scale are used.

Suppose a steady solenoidal velocity satisfies $\Delta \mathbf{u}=-\mathbf{u}$, and spatial averaging permits integration by parts without boundary terms. For a uniform test magnetic field $\mathbf{B}$, the small-magnetic Reynolds number expansion gives ${\mathbf{b}}_1=(\mathbf{B}\cdot \nabla )\mathbf{u}$ and $\Delta {\mathbf{b}}_2=-\nabla \times (\mathbf{u}\times {\mathbf{b}}_1)$. Using the self-adjointness of the Laplacian inside the average,The second equality concerns the averaged mean-field electromotive force; it does not assert that the quadratic forcing itself is monochromatic.

// Target: fluid-mechanics.bigb

Write the physical magnetic field as ${\mathbf{B}}_{\mathrm{phys}}$ and use the constant mass density ${\rho }_0$ to introduce Alfvén velocity units, $\mathbf{B}={\mathbf{B}}_{\mathrm{phys}}/\sqrt{{\mu }_0{\rho }_0}$. If $P$ is the fluid's mechanical pressure, the magnetohydrodynamic total pressure per unit mass density is $p=P/{\rho }_0+B^2/2$. The ideal magnetohydrodynamic equations then reduce toThe first equation is momentum conservation: the magnetic pressure has been included in $p$, leaving magnetic tension explicitly. The second is the ideal magnetohydrodynamic induction equation, expressing magnetic flux freezing. The last two equations impose incompressibility and the solenoidal magnetic-field constraint.

The model assumes a smooth conducting-fluid continuum with constant mass density and constant magnetic permeability, nonrelativistic velocity, and scales on which a single-fluid description is appropriate. Negligible electrical resistivity removes magnetic diffusion; negligible viscosity removes viscous momentum diffusion. Equivalently, the magnetic Reynolds number and the Reynolds number are large on the scales being retained. The low-frequency approximation neglects the displacement current; Hall and other kinetic corrections to Ohm's law are absent. Incompressibility is an imposed approximation, commonly justified when the Mach number is small. There is no independently evolved thermodynamic pressure: $p$ enforces the incompressibility condition. We omit external body forces; a conservative one could instead be absorbed into $p$.

A plasma contains mobile charged species whose collective electromagnetic response matters. On scales large compared with the Debye length, positive and negative charge densities nearly cancel: this is quasineutrality, not the absence of electric currents or electric fields. For electrons of temperature $T_e$ and number density $n_e$, the characteristic plasma screening and oscillation scales areA collective plasma description requires many particles within a screening volume; long-wavelength, sufficiently slow behavior can often be described by magnetohydrodynamics with an appropriate fluid closure.

An individual particle feels the Lorentz force $q(\mathbf{E}+\mathbf{v}\times \mathbf{B})$. In a uniform magnetic field, its transverse motion gyrates while its parallel motion remains free, giving a helical path when both components are present. The nonrelativistic cyclotron frequency and gyroradius areWhen the system scale is much larger than $r_L$ and the relevant time scale longer than ${\Omega }_c^{-1}$, the field strongly organizes transport across it, provided collisions do not interrupt the gyromotion. A perpendicular electric field also gives the common guiding-centre drift ${\mathbf{v}}_E=\mathbf{E}\times \mathbf{B}/B^2$.

For a conducting fluid, the macroscopic Lorentz force density has the useful decompositionThese terms are magnetic pressure and magnetic tension. Magnetic tension resists bending, while magnetic-pressure gradients resist nonuniform compression. The plasma beta $\beta =2{\mu }_{0}p/B^2$ compares thermal and magnetic stresses.

With uniform magnetic diffusivity ${\eta }_m$, the resistive induction equation isIf the magnetic Reynolds number $R_m=UL/{\eta }_m\gg 1$, ideal magnetohydrodynamics applies approximately and magnetic flux freezing links the field geometry to fluid motion. Finite diffusivity permits field slippage and magnetic reconnection. Field tension supports Alfvén waves with $v_A=B/\sqrt{{\mu }_0\rho }$; coupled magnetic and gas compression also supports magnetosonic waves. These effects allow a field to transmit forces and angular momentum between distant parts of an astrophysical plasma.

Separate the velocity field and magnetic field into mean and fluctuating parts,Assume the Reynolds averaging operation commutes with derivatives and has the usual product rules. Averaging the resistive induction equation givesThe mean-field electromotive force is the contribution through which small-scale motions affect the large-scale magnetic field. If the mean field varies slowly compared with the fluctuation scales, its response can be expanded in the mean field and its spatial derivatives. With homogeneous isotropic statistics, the leading terms areThe alpha effect is the term parallel to the mean field. Since an electromotive response is a polar vector whereas a magnetic field is axial, isotropic $\alpha$ is a pseudoscalar; reflection-symmetric statistics force it to vanish. Nonzero kinetic helicity supplies the required handedness. The turbulent diffusivity term transports and smooths the mean field; constant positive $\beta$ adds $\beta {\nabla }^2\overline{\mathbf{B}}$ to its evolution.

One controlled calculation is first-order smoothing. For negligible mean flow, omit the fluctuating nonlinear product in the fluctuation induction equation, obtainingFor incompressible flow this forcing is $(\overline{\mathbf{B}}\cdot \nabla ){\mathbf{u}}^{\prime }-({\mathbf{u}}^{\prime }\cdot \nabla )\overline{\mathbf{B}}$. A short-memory approximation of duration $\tau$ replaces the inverse temporal response by multiplication by $\tau$. Substitution into the mean-field electromotive force then uses isotropy:ConsequentlyThe first sign follows by taking the trace of the alpha-response tensor: ${ϵ}_{ijk}{u}_j^{\prime }{\partial }_i{u}_k^{\prime }=-{\mathbf{u}}^{\prime }\cdot \nabla \times {\mathbf{u}}^{\prime }$. This establishes the helicity formula for isotropic first-order smoothing. Its controlled regimes include small magnetic Reynolds number or sufficiently short fluctuation correlation time. Outside such a closure, the coefficients involve response-weighted time correlations; an instantaneous helicity formula is not universal.

Work with the complex representation and take its real part to obtain the physical magnetic field. For constant magnetic diffusivity and the prescribed circular velocity $\mathbf{u}=s\Omega (s)\hat{\mathbf{\varphi }}$, the resistive induction equation and $\mathbf{B}=\nabla \times (A\hat{\mathbf{z}})$ giveTo see the last equation directly, $\mathbf{u}\times \mathbf{B}=-s\Omega B_s\hat{\mathbf{z}}=-\Omega A_{\varphi }\hat{\mathbf{z}}$. Taking its curl reproduces the induction equation if $A_t=-\Omega A_{\varphi }+\eta {\nabla }^{2}A$, up to a time-dependent additive constant removable by the gauge transformation. This is the scalar form of advection and magnetic diffusion in cylindrical coordinates.

Put $\psi =\varphi -\Omega (s)t$. The material derivative of $a{e}^{i\psi }$ is $a_t{e}^{i\psi }$, because ${\psi }_t+\Omega {\psi }_{\varphi }=0$. A radial derivative instead acts on both amplitude and phase:Also $A_{\varphi \varphi }=-a{e}^{i\psi }$. Substitution gives the shearing-coordinate magnetic flux equation, here for azimuthal mode number one:Writing it with $a_t$ on the left is important when considering zero magnetic diffusivity: one must not literally divide by $\eta =0$. In the ideal limit, $a_t=0$, so $A(s,\varphi ,t)=A(s,\varphi -\Omega (s)t,0)$. This is magnetic flux freezing: the scalar flux pattern is carried with each rotating fluid ring. It does not make the magnetic field stationary. Indeed,so differential rotation winds the field and can make its azimuthal component grow linearly in time. Rigid rotation, for which ${\Omega }^{\prime }=0$, only rotates the pattern.

For the specified interior angular velocity, ${\Omega }^{\prime }=-{\Omega }_0/d$ and ${\Omega }^{\prime \prime }=0$ away from the circulation edge. Let $T={\Omega }_{0}t$, $x=s/d$ and $\mathrm{Rm}={\Omega }_0{d}^2/\eta$. The amplitude equation becomesFor a smooth initial amplitude on the cell scale, at $x$ bounded away from zero and the edge, the $-T^{2}a$ term dominates when $T\gg 1$. Its accumulated damping is order one at $T=O({\mathrm{Rm}}^{1/3})$, whereas the integrated terms proportional to $T$ are only $O({\mathrm{Rm}}^{-1/3})$ and the terms independent of $T$ are $O({\mathrm{Rm}}^{-2/3})$. The resulting cubic-time resistive damping in a differentially rotating cell isThe magnetic field has this leading exponential decay envelope. Its components also contain the algebraic winding factor visible in $B_{\varphi }$, so the claim concerns the cubic exponent, not an exact amplitude with no prefactor. Equivalently the radial wavenumber grows as $|k_s|\simeq |{\Omega }^{\prime }|t$ and the integrated resistive damping is ${\int }_0^t\eta k_s^{2}dt=\eta ({\Omega }^{\prime }{)}^2{t}^3/3$.

This is phase mixing in magnetic flux expulsion. The circulation creates ever finer gradients, allowing weak magnetic diffusion to act much sooner than its unsheared time $t_{\eta }=d^2/\eta$:Rapid winding followed by resistive smoothing expels the imposed field from the cell interior. The maintained exterior magnetic field is redistributed into a circulation-edge boundary layer. The result is an interior high-magnetic Reynolds number transient, not a uniform assertion at the axis and the nonsmooth edge, nor the exact $t\to \infty$ solution with a continuously imposed far field. The diffusion length at ${\tau }_c$ is $\sqrt{\eta {\tau }_c}/d=3^{1/6}{\mathrm{Rm}}^{-1/3}$, identifying where a local bulk approximation can fail. The linear profile also needs smoothing at a regular rotation axis and at $s=d$ if used as a fully smooth global flow.

In photospheric magnetoconvection, circulating and diverging convection transports weak magnetic flux towards converging downflow regions. Solar granulation concentrates field in intergranular lanes, while larger-scale supergranulation helps form the solar magnetic network. Flux expulsion explains why a highly conducting photosphere can have relatively weak field in cell interiors and intermittent concentrations at their edges: high electrical conductivity delays diffusion on the original scale, but does not prevent diffusion on the fine scales created by advection. A long-lived circulation must last long enough for the expulsion time; rapidly changing cells need not reach this idealized state. Once the field becomes strong, the Lorentz force modifies the convection, so this prescribed-flow calculation alone neither determines the final field strength nor describes strong-field sunspots.

Take the historical scope to be the state of solar dynamo theory around 2001. A successful explanation must account for organized patterns, not merely the existence of a magnetic field. For the Sun, the main constraints are the roughly eleven-year sunspot cycle and twenty-two-year magnetic solar cycle; opposite leading polarities in the two hemispheres and reversal between successive cycles according to Hale's polarity law; the equatorward drift of the activity belts in the solar butterfly diagram; the systematic bipolar-region tilt of Joy's law; and reversal of the large-scale polar field near activity maximum. The field is intermittent, with strong sunspots coexisting with weaker surface magnetic flux, and cycle amplitudes vary, including extended solar grand minima such as the Maunder minimum. These are separate constraints on regeneration, migration, emergence, saturation and long-term variability.

The fundamental success is a physically viable source of continuing magnetic energy. In a highly conducting star the resistive induction equation permits stretching and folding to transfer kinetic energy into magnetic energy; magnetic diffusion allows changes of topology and opposes amplification. A kinematic magnetic dynamo can have growing solutions even though the velocity is prescribed, demonstrating the regenerative mechanism. The Cowling anti-dynamo theorem rules out a self-sustaining entirely axisymmetric magnetic configuration; an approximately axisymmetric large-scale field can nevertheless be maintained by nonaxisymmetric convection and its correlations. A relic field undergoing passive diffusion does not by itself explain repeated organized reversals and the systematic relation of magnetic activity to rotation and convection.

A tractable description is a mean-field dynamo. Split velocity and field into means and fluctuations. Averaging the induction equation givesThe mean-field electromotive force is where unresolved correlations enter. A simple local isotropic closure is $\mathcal{E}=\alpha \overline{\mathbf{B}}-{\eta }_t\nabla \times \overline{\mathbf{B}}$. The alpha effect couples the toroidal magnetic field back to a poloidal magnetic field, while the Omega effect produces a toroidal magnetic field from differential rotation:In the short-correlation, locally isotropic version of first-order smoothing, $\alpha \simeq -\tau ⟨{\mathbf{u}}^{\prime }\cdot \nabla \times {\mathbf{u}}^{\prime }⟩/3$ and ${\eta }_t\simeq \tau ⟨|{\mathbf{u}}^{\prime }{|}^2⟩/3$. Rotation and stratification can produce kinetic helicity with opposite signs in the two hemispheres, allowing a hemispherically organized alpha effect. These formulae explain a possible mechanism, but their controlled assumptions are not automatically satisfied by stellar turbulence at very large magnetic Reynolds number.

An alpha-Omega dynamo naturally supports oscillation and migration rather than just monotone amplification. For example, in the northern hemisphere choose local right-handed directions $x$ equatorwards, $y$ azimuthally and $z$ radially outwards. With $U_y=Sz$ and poloidal potential $A\hat{\mathbf{y}}$, the local equations areFor $e^{pt+ikx}$ the local alpha-Omega dynamo wave dispersion isThe growing Parker dynamo wave has phase velocity $-\mathrm{Im}p/k$, equatorward for $\alpha S<0$. Thus a migration direction, an excitation threshold and a finite oscillation period arise from the feedback between the alpha effect and shear, rather than being independently imposed. The local example illustrates the mechanism; it does not determine the global stellar period, parity or latitude range. Dipole and quadrupole parity in a mean-field dynamo depend on the spatial coefficients, coupling between hemispheres and boundary conditions.

A major difficulty is matching that mechanism to the measured internal rotation. Helioseismology constrains differential rotation and the tachocline, giving a plausible strong-shear region near the base of the convection zone. At low northern latitudes a positive radial shear combined with the usual positive convective $\alpha$ would give poleward propagation in the simple local wave model, rather than the equatorward solar butterfly diagram. A different sign or location of the alpha effect, an interface solar dynamo, and magnetic transport must therefore be examined rather than assuming that any alpha-Omega dynamo predicts the observed migration. The stable region below the convection zone may store strong toroidal magnetic flux; magnetic buoyancy can then produce rising tubes and bipolar active regions. Flux-tube emergence also has to account for their low latitudes and Joy's law tilts, not simply for a large interior field.

The Babcock-Leighton mechanism supplies a more directly observable route to regeneration. Tilted bipolar regions disperse; cancellation of opposite polarities and poleward transport of surviving flux alter the polar poloidal magnetic field. Differential rotation subsequently winds that field into a toroidal magnetic field. A flux-transport solar dynamo combines this surface source with meridional circulation in a star and turbulent magnetic diffusivity. Models with poleward surface transport and an equatorward deep return flow can reproduce an equatorward activity belt together with poleward surface-flux migration and an appropriate phase relation of the polar field. This was already an explicit model construction in the 1999 flux-transport calculation. Its success remains conditional on the assumed circulation, transport coefficients and emergence prescription; the internal return flow was not then measured well enough to make these quantities unique. The source requires emerged active regions and therefore also leaves the recovery from a state with almost no spots as a nontrivial problem.

Saturation is another limitation of kinematic models. Exponential growth cannot persist indefinitely: the Lorentz force alters the shear and convection, and changes the correlations that generate the alpha effect. Algebraic dynamo quenching can produce bounded cycles but fitting a quenching coefficient is not a first-principles amplitude prediction. Near conservation of magnetic helicity in a highly conducting fluid adds constraints on the growth and saturation of a large-scale field; boundaries and helicity transport matter. Magnetic buoyancy, intermittency and the interaction of mean and fluctuating fields likewise make a local isotropic closure incomplete. Fluctuating regeneration and flow, or nonlinear interactions, can produce cycle irregularity and solar grand minima, but reproducing one irregular time series does not establish a unique explanation of the Maunder minimum.

Other solar-type stars provide an important independent test. Long-term chromosphere observations already showed both cyclic and irregular activity, with rapidly rotating young stars tending to be more active and older slower rotators more often showing smooth cycles; see the original stellar activity survey. These qualitative trends support a rotational-convective origin and stellar rotation-activity feedback through wind-driven magnetic braking of a solar-type star. The relevant control need not be rotation alone: the stellar activity Rossby number, comparing rotation and the convection time, expresses why stellar structure matters. A predictive theory should explain the diversity of cycle periods and geometries, not assign a separate arbitrary $\alpha$, diffusivity and circulation to every star. Stellar activity measurements also trace magnetic heating or spots, not automatically the full internal magnetic geometry.

By 2001, dynamo theory explained credible mechanisms for regeneration, polarity cycling and migrating activity, and suitably constructed models reproduced important solar patterns. It had not supplied a unique, self-consistent quantitative explanation of the full solar and stellar activity phenomenology. The gap lay in determining transport and regeneration from the actual stellar flows, nonlinear saturation and emergence, and robust predictions across different stars. Three-dimensional calculations supported magnetic generation, but numerical diffusivities and the restricted range of resolved scales limited direct extrapolation to stellar conditions. Agreement of a prescribed-coefficient mean-field dynamo with selected observations is meaningful evidence of feasibility, not a demonstration that its detailed mechanism or parameters have been uniquely established.

Magnetic flux concentration and flux expulsion describe redistribution of an existing field, not necessarily creation of net flux. For a weak field in a prescribed conducting flow, the resistive induction equation isIts ideal part combines transport, stretching and compression:Thus a converging flow can intensify a component, and stretching can increase magnetic energy even if the signed net flux is zero. At large magnetic Reynolds number $\mathrm{Rm}=UL/\eta$, diffusion is initially slow on the eddy scale, so magnetic flux freezing carries field lines with the flow. These statements do not imply that an ideal closed eddy can simply lose its initially enclosed flux.

The missing step is the production of small scales by differential circulation. In a two-dimensional incompressible eddy, the Cartesian magnetic flux function obeysDifferent streamlines have different circulation rates. An initially smooth field is wound into fine folds; diffusion and reconnection then smooth the variation of $A$ inside the eddy. As $A$ becomes nearly uniform there, its gradient and hence the interior magnetic field become small. The surviving flux lies mainly near the circulation boundary and adjacent separatrices, where stronger gradients remain. This is flux expulsion. With exactly zero diffusivity, field topology and the flux of material loops remain frozen, so a relaxed flux-free interior is not obtained by pure advection alone. Rigid rotation also fails to generate the differential winding responsible for this process.

The phase mixing in magnetic flux expulsion estimate makes the enhanced role of small diffusion explicit. If the angular-velocity gradient is of order $U/L^2$, a wound radial wavenumber grows as $k_r\sim Ut/L^2$. The accumulated damping exponent is of orderIt becomes order one atThis is much shorter than the large-scale diffusion time $(L/U)\mathrm{Rm}$ at large $\mathrm{Rm}$, but longer than one turnover time. It is a kinematic differential-shear estimate, not a universal time for every turbulent or rapidly evolving convection cell.

A related concentration mechanism occurs at converging flow boundaries. For a locally vertical bundle, conservation of flux gives $B\mathcal{A}\simeq \mathrm{constant}$, so reducing its cross-section increases $B$. In a simple converging boundary region with strain rate $U/L$, balancing strain and transverse diffusion givesThis is the thickness estimate for that local strain-diffusion model; it need not be the thickness or energy scaling of every closed-eddy flux-expulsion problem. Stronger field eventually exerts Lorentz force feedback and changes the flow. The equipartition magnetic field $B_{\mathrm{eq}}=\sqrt{{\mu }_0\rho }U$ is a useful indication that the kinematic assumption can fail, not a universal upper bound on concentrated fields.

At the Sun's photosphere, hot upflows form solar granules, and their spreading horizontal motion transports weak vertical field toward cooler intergranular lanes and downflows. A thin-layer illustration, neglecting vertical transport terms, is ${\partial }_t{B}_z+{\nabla }_h\cdot ({\mathbf{u}}_h{B}_z)=0$. It gives $D_h{B}_z/Dt=-B_z{\nabla }_h\cdot {\mathbf{u}}_h$: horizontal convergence intensifies the vertical flux density. This surface transport is related to the redistribution by convective eddies, but is not identical to the closed-streamline mechanism in the preceding planar example. Larger-scale supergranulation also gathers magnetic elements toward its boundaries, helping organize the solar magnetic network.

Concentration can be further intensified by compressible thermal dynamics. In the thin-tube picture of convective collapse, an unstable downward flow drains gas, reduces internal pressure and allows contraction at approximately conserved flux. Lateral balance $p_i+B^2/(2{\mu }_0)\simeq p_e$ then permits field stronger than the local kinetic equipartition estimate. This evacuation mechanism is excluded by a strict Boussinesq approximation. It is not a description of every evolving surface concentration: gas pressure, dynamic stresses and magnetic curvature can all matter when thin-tube equilibrium fails. The resulting strong, intermittent magnetic elements are associated with intergranular structure; their radiative appearance also depends on thermal structure and opacity. The organized fields of sunspots require a broader account of flux supply and magnetic feedback, rather than attribution to weak-field expulsion alone.

Thus convection redistributes and concentrates flux, finite diffusion enables classical expulsion, and nonlinear magnetic and thermal feedback determines the strong-field outcome. An expulsion argument alone neither supplies the original solar field nor establishes a self-sustaining solar dynamo.

Choose the planar magnetic field convention $\mathbf{B}=\nabla \times (A{\mathbf{e}}_z)$, so $B_s=s^{-1}{A}_{\varphi }$ and $B_{\varphi }=-A_s$. The resistive induction equation, with constant magnetic diffusivity and a suitable time-dependent gauge constant removed, reduces to advection and diffusion of this Cartesian magnetic flux function:For the specified differential rotation, its mth azimuthal Fourier mode obeysThe physical field is the real part of the complex mode. Substitute $a=s^{m}g(t)\exp [-im{\Omega }_{0}t-i{s}^{2}f(t)]$. A direct radial differentiation givesOn the left, rigid rotation cancels, leaving $\dot{g}/g-i{s}^2\dot{f}+im\omega s^2$. Equating the constant and quadratic coefficients yieldsThe initial condition fixes $f(0)=0$ and $g(0)=1$, since the constant C was already factored out. Set ${\mu }^2=4i\eta m\omega$. Then the exact quadratic-shear magnetic flux solution follows by differentiation:Indeed $\dot{f}=m\omega {\mathrm{sech}}^2(\mu t)$ and $4i\eta f^2=m\omega {\tanh }^2(\mu t)$. Also $\dot{g}/g=-(m+1)\mu \tanh (\mu t)=-4i\eta (m+1)f$. For an unspecified initial normalization the second formula has a multiplicative constant $C^{\prime }$, but here $C^{\prime }=1$.

For $\eta >0$ and $\omega >0$, choose $\mu =(1+i)\sqrt{2\eta m\omega }$. The asymptotic condition is $\mathrm{Re}(\mu )t\gg 1$: then $f$ tends to the finite constant $m\omega /\mu$ and $g\sim 2^{m+1}\exp [-(m+1)\mu t]$. At any fixed s, differentiating A to obtain the magnetic field multiplies it by factors involving $m/s$ and $sf$, which have no additional exponential time dependence. Therefore its envelope decays asUsing the paper's magnetic Reynolds number definition gives the precise scalingThus the decay time scales as ${\mathrm{Rm}}^{1/2}/{\Omega }_0$ for positive ${\Omega }_0$ and fixed m. The mechanism is magnetic phase mixing under differential rotation: winding produces progressively fine radial structure on which magnetic diffusion acts efficiently. For negative $\omega$, choose the square root with positive real part, replacing $\omega$ by $|\omega |$ in the decay time. If $\omega =0$, the initial harmonic mode simply rotates, with $f=0$, $g=1$; there is no shear-enhanced decay. The printed positive-root and Reynolds-number expressions implicitly use $\omega >0$.

Take zero mean flow locally and incompressible fluctuations. Subtracting the averaged induction equation givesThe first-order smoothing approximation omits the second term on the right when solving for $\mathbf{b}$, but retains $⟨\mathbf{u}\times \mathbf{b}⟩$ in the mean equation. For a locally uniform, slowly changing mean magnetic field the forcing becomes $(\overline{\mathbf{B}}\cdot \nabla )\mathbf{u}$.

There are two controlled ways to make the omitted fluctuation term small. At small eddy magnetic Reynolds number $R_m=u_{\mathrm{rms}}\ell /\eta$, diffusion gives $b/B\sim R_m$, so the omitted $ub/\ell$ is smaller than the retained $uB/\ell$ by $R_m$. Alternatively, if $u_{\mathrm{rms}}\tau /\ell \ll 1$ for the fluctuation response time $\tau$, then $b/B\sim u_{\mathrm{rms}}\tau /\ell$ and the short-memory response is perturbative. Mean-field scale separation is additionally required. Merely calling the turbulence weak or assuming that the mean magnetic field is small does not justify the omission.

With homogeneous initial transients removed, the causal linear solution isInsert this into the electromotive correlation. Isotropy gives ${\alpha }_{ij}=\alpha {\delta }_{ij}$, so taking the trace yields the helicity formula for isotropic first-order smoothing:Using the alternating-tensor sign givesThe minus sign follows from ${ϵ}_{imn}=-{ϵ}_{min}$. This derivation specifies the diffusion and correlation weighting, rather than asserting that alpha is proportional to instantaneous helicity for arbitrary turbulence.

If $\eta \tau /{\ell }^2\ll 1$, diffusion over a correlation time is negligible. Define the helicity correlation time byFor an ordinary decaying correlation with nonzero mean kinetic helicity, this gives the familiar estimateIt is the kinetic helicity, not magnetic helicity. Opposite flow handedness reverses the alpha effect. The time-integral formula remains meaningful if the equal-time helicity vanishes or its correlation changes sign, when replacing it by a single positive correlation time may fail. In the same weak-diffusion short-memory approximation the gradient forcing $-(\mathbf{u}\cdot \nabla )\overline{\mathbf{B}}$ gives $\beta =\frac{1}{3}{\int }_0^{\infty }⟨\mathbf{u}(t)\cdot \mathbf{u}(t-\tau )⟩d\tau \simeq {\tau }_U{u}_{\mathrm{rms}}^2/3$; the two correlation times need not be equal.

In the small-$R_m$, slowly varying quasistatic regime, diffusion is instead essential. The response givesFor velocity concentrated at wavenumber $k_0$, this becomes $\alpha =-⟨\mathbf{u}\cdot \mathbf{\varpi }⟩/(3\eta k_0^2)$. Its effective memory is the magnetic diffusion time on the eddy scale. The controlled low-$R_m$ setting is developed in Moffatt's analysis of turbulent dynamo action; the causal-response calculation above also covers short-memory turbulence without assuming that diffusion is always negligible.

The geodynamo operates in the electrically conducting liquid of Earth's outer core. buoyancy from thermal and compositional convection supplies motion, while the Coriolis force organizes rotating flow. Rotation together with radial stratification permits kinetic helicity, typically with opposite signs in the two hemispheres. In mean-field language this allows an alpha effect that converts poloidal and toroidal components into one another. A large-scale alpha-squared dynamo is therefore a plausible regeneration mechanism even without dominant differential rotation.

The growth condition compares this regeneration with the sum of molecular and turbulent magnetic diffusion on the core scale. A spherical model must impose the appropriate magnetic matching at the core boundaries and permit spatially varying, generally anisotropic alpha and diffusivity tensors. Rotation selects an axis, and stratification and boundaries select positions, so homogeneous isotropic turbulence is only a local idealization. Such a model can favor a predominantly axial dipole, but the preferred parity and threshold must be found from its global eigenproblem rather than deduced from the plane-wave formula alone. If mean shear is significant, alpha–Omega regeneration can also contribute.

The first-order-smoothing helicity formula is a useful guide to the sign and magnitude of alpha, not an automatically controlled numerical closure for the core. The important convective motions need not have small eddy magnetic Reynolds number or short enough correlation time; scale separation can also be limited. Furthermore, the Lorentz force modifies the flow and its correlations once the field grows. Saturation therefore requires nonlinear feedback, including dynamo quenching and magnetic-helicity constraints where relevant.

The observed field varies and occasionally reverses. A constant-coefficient growing kinematic eigenmode does not itself explain these reversals or select a saturated amplitude. Spatial and temporal variations, competing modes, and nonlinear flow–field interactions must be included. mean-field dynamo theory supplies a physical regeneration mechanism and a testable threshold, while the geodynamo requires a rotating, bounded, nonlinear realization of that mechanism.

Let $⟨\cdot ⟩$ denote the spatial average. Assume periodicity, or an equivalent homogeneous averaging limit in which total derivatives average to zero and integration by parts is valid. Fix the induced field to have zero mean. Expand the solution of the resistive induction equation in small magnetic Reynolds number:The uniform test magnetic field commutes with differentiation. Since $\Delta \mathbf{u}=-\mathbf{u}$, the first equation has solution ${\mathbf{b}}_1=(\mathbf{B}\cdot \nabla )\mathbf{u}$.

The quadratic second-order forcing need not be monochromatic, so one must not simply replace its inverse Laplacian by multiplication by $-1$. Instead use the self-adjointness of the Laplacian under the average:This proves the monochromatic small-Reynolds-number mean electromotive force expansion

For the symmetry of the first-order monochromatic alpha tensor, writeIts antisymmetric part is determined by contraction with the Levi-Civita symbol. The contraction identity givesThe second term vanishes by incompressibility; the first is $⟨{\partial }_j(u_j{u}_p)⟩=0$. Since every three-dimensional antisymmetric second-rank tensor is equivalent to its contracted axial vector, ${\alpha }^{(1)}$ is symmetric. No Fourier series or Fourier transform is required.

To derive the next identity, set $\mathbf{C}=\mathbf{u}\times (\mathbf{B}\cdot \nabla )\mathbf{u}$. Expanding the curl and cross product in components,The average of the second term is zero: integration by parts makes it $-⟨C_i{\partial }_j{u}_j⟩$. Integrating the first term by parts givesExtracting the coefficient of $B_j$ now yieldsInterchanging $i$ and $j$ reverses the cross product, proving the antisymmetry of the second-order monochromatic alpha tensor:The second coefficient of the alpha tensor can therefore contribute turbulent magnetic pumping, while its symmetric part vanishes at this order.

Magnetoconvection is the competition between thermally driven convection and the force and transport of a conducting fluid's magnetic field. An imposed magnetic field resists motions that bend its field lines through magnetic tension, while magnetic diffusion permits slippage. The competition changes the onset threshold, preferred spatial scale, temporal behaviour and finite-amplitude organization. The geometry matters: a disturbance nearly invariant along a horizontal imposed magnetic field can avoid much of its tension, whereas an impermeable layer with a vertical imposed magnetic field cannot generally avoid bending it.

A useful quantitative example is a nonrotating Boussinesq layer of depth $d$, heated from below with temperature drop $\Delta T>0$ and threaded by a vertical magnetic field $B_0\hat{\mathbf{z}}$. Let $\nu$, $\kappa$, $\eta$ denote kinematic viscosity, thermal diffusivity and magnetic diffusivity. In thermal-diffusion units define the Rayleigh number, Prandtl number, magnetic-to-thermal diffusivity ratio and Chandrasekhar number byUse stress-free impermeable boundaries, fixed temperature perturbation zero, and magnetic boundary conditions keeping the field vertical. More explicitly, the normal modes below satisfy $w={\partial }_z^{2}w=\theta =b_x=b_y={\partial }_z{b}_z=0$ at $z=0,1$. For a horizontal wavenumber $a>0$, take $w,\theta \propto \sin \pi z$, $b_z\propto \cos \pi z$, and put $s=a^2+{\pi }^2$. The linear temperature and resistive induction equation give $(\lambda +s)\theta =w$ and $(\lambda +\zeta s)b_z={\partial }_{z}w$. Applying the solenoidal projection to the magnetohydrodynamic momentum equation removes pressure; the magnetic coupling is $Q\zeta {\partial }_z\mathbf{b}$, and the vertical projected buoyancy has coefficient $R{a}^2/s$. Substituting the temperature and magnetic responses gives the vertical-field magnetoconvection dispersion relationThese stated boundary conditions make the single vertical mode exact in this model; other boundaries change numerical thresholds.

For steady marginal convection, set $\lambda =0$. The threshold curve isThe imposed magnetic field adds a positive penalty. Differentiating with respect to $x=a^2$ shows that the minimizing wavenumber obeysAt $Q=0$ this gives $a_c=\pi /\sqrt{2}$ and $R_c=27{\pi }^4/4$. With increasing $Q$, the root moves to larger $x$ because the left-hand side is strictly increasing for $x\ge {\pi }^2/2$. The strong-field wavenumber selection in magnetoconvection limit isThus a strong vertical imposed field raises the steady threshold and selects narrower convection cells. At small horizontal wavenumber the incompressibility constraint makes a given vertical motion require large horizontal displacement and magnetic bending. Increasing horizontal wavenumber reduces this magnetic penalty, but increases viscous and thermal diffusion. Their balance selects the narrow cells. This statement is specific to the vertical-field geometry: wavenumber selection in oblique-field magnetoconvection can instead favour nearly field-invariant motions.

Oscillatory marginality in vertical-field magnetoconvection occurs because magnetic and temperature perturbations need not relax at the same rate. To exhibit the phase balance, substitute $\lambda =i\omega$ with $\omega \ne 0$ and write the neutral condition asThe imaginary part must vanish because $R$ is real. Dividing that imaginary part by $\omega s$ givessoA nonzero oscillation consequently requires $0<\zeta <1$ and a sufficiently strong imposed magnetic field. In this regime thermal diffusion can erase the temperature anomaly faster than magnetic diffusion relaxes bent field lines. The magnetic restoring force can outlast the buoyant anomaly, reverse the motion, and maintain the phase lag needed for oscillatory instability. The real part gives the oscillatory neutral curveOnly its portions with ${\omega }^2>0$ are physical. The first onset is the lower of the minimized steady and admissible oscillatory thresholds; the existence of an oscillatory neutral solution at a given wavenumber alone does not establish that it is first. For example, at fixed $P>0$ and $0<\zeta <1$, both minimizing wavenumbers grow like $Q^{1/6}$ for large $Q$, while $R_s\sim {\pi }^{2}Q$ and $R_o\sim D{\pi }^{2}Q$ with $D=\zeta (P+\zeta )/(1+P)<1$. At those oscillatory minima, the positive order-$Q$ contribution to ${\omega }^2$ dominates its order-$Q^{2/3}$ negative contribution, so the oscillatory branch is genuinely admissible and becomes the lower threshold. An overstability therefore precedes steady onset in this limit. More generally nonlinear convection can oscillate even when linear onset is steady; the linear calculation is not a classification of every time-dependent state.

The nonlinear regime introduces a different feedback. The magnetic Reynolds number $\mathrm{Rm}=Ud/\eta$ measures induction relative to diffusion for a finite convective velocity $U$. At large but finite $\mathrm{Rm}$, persistent convection rolls can redistribute magnetic flux. In a two-dimensional incompressible roll, a flux function $\mathcal{A}$ obeysCirculation repeatedly stretches the flux pattern, building gradients on which even small magnetic diffusivity acts efficiently. Over the flux-rearrangement time this permits flux expulsion: the interior of a vigorous roll becomes weakly magnetized, while magnetic flux is concentrated near its edges. The weak-field interior then suffers less magnetic inhibition and can support broader, stronger convection than the linear narrow-cell prediction for the original uniform magnetic field. Exact zero resistivity would preserve material flux and cannot by itself justify the relaxed flux-expelled state.

The weak interior field also follows from a steady leading-order balance for a coherent two-dimensional roll. Let $\psi$ be its streamfunction and suppose its interior streamlines are nested closed curves with $\nabla \psi \ne 0$. In the small-diffusivity interior, $\mathbf{u}\cdot \nabla {\mathcal{A}}_0=0$ implies ${\mathcal{A}}_0=F(\psi )$. Integrate the full steady flux equation over the region enclosed by one streamline. Because $\mathbf{u}\cdot \mathbf{n}=0$ there, the advection integral vanishes; the divergence theorem then gives $\oint {\partial }_n\mathcal{A}ds=0$. In a regular interior expansion this requires $F^{\prime }(\psi )\oint {\partial }_n\psi ds=0$. The normal derivative of $\psi$ has one nonzero sign around a regular nested streamline, so $F^{\prime }(\psi )=0$. Hence the leading interior flux function is constant, and its spatial gradient, which gives the two-dimensional magnetic field, vanishes. The mismatch with external flux is accommodated by magnetic boundary layers. This argument assumes the coherent steady roll and a regular interior; it does not prove that every turbulent flow expels flux.

The mean imposed vertical magnetic flux is conserved for horizontally periodic boundaries. In a schematic two-region picture, let a fraction $f$ have weak field $B_c$ in vigorous convection and the remainder have field $B_q$. ThenThe intensified field outside the active cells inhibits surrounding convection. This gives magnetic flux separation: vigorous weak-field patches coexist with strong-field regions whose motion is small or occurs on a finer scale. The redistribution is a consequence of nonlinear induction and back reaction, rather than disappearance of the net imposed flux.

This feedback can permit subcritical magnetoconvection. Infinitesimal perturbations decay below the conductive state's linear threshold because they encounter the full imposed field. A finite-amplitude roll can nevertheless persist by first expelling enough flux to reduce its own magnetic restraint. A sufficient finite disturbance can then switch between the stable conductive state and a convective state, producing hysteresis. A schematic amplitude equation makes the distinction between local linear stability and finite-amplitude persistence explicit:Here $\mu$ measures distance from linear onset; this is an illustrative normal form, not a claim that these coefficients have been derived for every magnetic layer. Writing $r=|A{|}^2$, the nonzero equilibria satisfyThey appear at a saddle-node bifurcation at $\mu =-c^2/(4d)$. For $-c^2/(4d)<\mu <0$, the zero state is stable, the lower radial-amplitude branch is unstable, and the upper branch is stable because its radial derivative has sign $c-2d{r}_{+}<0$. Thus finite-amplitude convection can coexist with a linearly stable conductive state below the onset threshold. Large $\mathrm{Rm}$ facilitates the flux-expulsion mechanism but is not a universal sufficient condition: geometry, the imposed magnetic field, the capacity to form coherent recirculating cells, and magnetic back reaction all matter. The broad nonlinear cells in weak-field regions can coexist with much narrower cells in magnetically dominated regions, linking the nonlinear flux distribution back to scale selection.

Let ${\mathbf{x}}^{\prime }=\mathbf{x}-\mathbf{v}t$, $t^{\prime }=t$ for a constant relative velocity $\mathbf{v}$. The Galilean transformation gives ${\nabla }^{\prime }=\nabla$ and ${\partial }_{t^{\prime }}={\partial }_t+\mathbf{v}\cdot \nabla$. In the magnetic Galilean limit of the pre-Maxwell equations, takeThe Ampère's law without displacement current and the solenoidal constraint are immediately unchanged. For Faraday's law, the constant-velocity vector identity is $\nabla \times (\mathbf{v}\times \mathbf{B})=\mathbf{v}\nabla \cdot \mathbf{B}-(\mathbf{v}\cdot \nabla )\mathbf{B}$. ThusThese are transformations of the magnetic nonrelativistic approximation. The exact charge-current transformation also contains the convective term $-{\rho }_e\mathbf{v}$ in ${\mathbf{J}}^{\prime }$; ignoring it requires the quasineutral magnetic limit. They are not exact transformations of the full Maxwell equations.

Since ${\mathbf{u}}^{\prime }=\mathbf{u}-\mathbf{v}$, the combination ${\mathbf{E}}^{\prime }+{\mathbf{u}}^{\prime }\times {\mathbf{B}}^{\prime }$ equals $\mathbf{E}+\mathbf{u}\times \mathbf{B}$. For an ideally conducting fluid the rest-frame electric field vanishes, so the ideal moving-conductor Ohm law is $\mathbf{E}+\mathbf{u}\times \mathbf{B}=0$. Substitution in Faraday's law givesThe ideal magnetohydrodynamic induction equation requires nonrelativistic fluid motion, negligible displacement current, and negligible nonideal terms in the electric-field relation. In particular, magnetic diffusion must be slow compared with advection: the magnetic Reynolds number $\mathrm{Rm}=UL/{\eta }_m\gg 1$, where ${\eta }_m=({\mu }_0\sigma {)}^{-1}$ is the magnetic diffusivity. The single-fluid moving-conductor Ohm law must apply on the scales considered, with corrections to that relation negligible; viscosity is not itself excluded by this induction equation.

For a smooth ideal evolution, magnetic flux freezing states that the magnetic flux through every material surface is conserved and that magnetic field lines are carried with the fluid. Their connectivity is preserved: ideal smooth advection does not reconnect them. These implications are stated here without proof, as requested.