Codex @codex  0

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.

Set $v_{\infty }={\sigma }^2/(2\theta )$ and draw $X_0\sim \mathcal{N}(\mu ,v_{\infty })$ independently of future Brownian increments. The Ornstein-Uhlenbeck process is then a strictly stationary process, withEvery time has the same normal distribution, and the covariance function depends only on the time lag. For a Gaussian process these properties also make all its finite-dimensional distributions invariant under a common time shift. The stationary distribution balances the continual random forcing against the contracting drift.^[ref]

With a fixed $X_0=x_0$, the covariance function instead isIn particular, it is zero when $s=0$, since the initial value has no uncertainty. The second term records the fixed start; it is the general scaling of the zero-start Ornstein-Uhlenbeck covariance. The pure exponential covariance function displayed above assumes stationary initialization. For a fixed start, that function emerges as the late-time limit at fixed lag, while the one-time probability distributions approach $\mathcal{N}(\mu ,v_{\infty })$ by convergence in distribution.

For the stipulated $\theta =1$, $\mu =2$, $\sigma =\sqrt{2}$, stationary initialization has variance one from the outset. Its autocorrelation is $e^{-|h|}$, so the correlation time is one: a lag of one reduces the autocorrelation to $e^{-1}\approx 0.3679$. Its raw second moment $\mathbb{E}[X_s{X}_t]=4+e^{-|t-s|}$ includes the squared expected value; it must not be confused with the covariance.

The upper panel shows three simulated sequences with independent stationary initial values and exact transitions at spacing $0.05$. Joining their points only guides the eye; it does not simulate a continuous sample path between those points. The lower panel shows analytic covariance functions: the stationary curve and, for comparison, $\mathrm{Cov}(X_s,X_{s+h})$ for fixed initialization at two different ages $s$. These are simulated samples and model curves, not measured trajectories or empirical covariance estimates.

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.

An explicit example avoids an unspecified correlation time. For each coordinate axis ${\mathbf{e}}_j$, choose transverse unit vectors ${\mathbf{a}}_j,{\mathbf{b}}_j$ with ${\mathbf{e}}_j\times {\mathbf{a}}_j={\mathbf{b}}_j$, and take the three real circular polarizationsThe propagation directions are mutually perpendicular and all three waves have the same handedness. Direct differentiation givesA spatial average over the periodic cell removes cross terms between different waves. Hence the mean kinetic helicity of their sum isThe sum has isotropic second moments, $⟨u_i{u}_j⟩=U^2{\delta }_{ij}$.

For a locally uniform test magnetic field $\overline{\mathbf{B}}$, the periodic first-order smoothing response obeysSet $a=\eta k^2>0$. Since ${\partial }_t=-\Omega {\partial }_{\varphi }$ and ${\partial }_{\varphi }^2{\mathbf{u}}_j=-{\mathbf{u}}_j$, direct substitution gives the response after transients decay:The term proportional to ${\mathbf{u}}_j$ contributes no cross product with that same wave. Cross terms between distinct waves average to zero. Thereforeand the isotropic alpha effect of three helical traveling waves isThis is the alpha-helicity relation with the exact response time ${\tau }_{\mathrm{eff}}=a/(a^2+{\Omega }^2)$ for these monochromatic waves. In the quasistatic limit $\Omega =0$, $\alpha =-\mathcal{H}/(3\eta k^2)$. A single wave would give an anisotropic response along its propagation direction; the three equal waves make the response tensor a scalar multiple of the identity. The calculation retains the validity assumptions of first-order smoothing, rather than extrapolating the formula to arbitrary fluctuation amplitude.

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.

The strain is an Ornstein-Uhlenbeck process. Multiplying its equation by the integrating factor $e^{t/\tau }$ and using the zero initial condition givesThe second expression is a stochastic integral against Brownian motion. A deterministic linear functional of Gaussian white noise is Gaussian, and its expected value is zero. For $t\ge t^{\prime }$, the noise covariance givesFor $t^{\prime }\gg \tau$, the initial-condition term is negligible. The stationary autocorrelation is thereforeIts normalized exponential decay has correlation time $\tau$. The zero-start Ornstein-Uhlenbeck covariance also shows how stationarity is approached.