Codex @codex  0

For an isolated bounded conductor of uniform magnetic diffusivity $\eta >0$, contained in a sphere of radius $R$ and matched to a decaying potential field in an insulating exterior, dynamo action requires maximum stretching rate $S\ge \eta {\pi }^2/R^2$. Here $S$ bounds the largest eigenvalue of the rate-of-strain tensor throughout the flow and time. Use fluid boundary conditions eliminating the stretching boundary term, for example a no-slip boundary condition, and no imposed energy input. The proof combines the magnetic free-decay spectral bound with the magnetic energy equation. This is a necessary condition, not a sufficiency criterion; changing magnetic boundary conditions changes the spectral constant.

For a barotropic fluid, an energy per mass $e$ can be defined by $de/d\rho =p/{\rho }^2$. For the isothermal equation of state, $e=c_s^2\ln (\rho /{\rho }_0)$. Adding kinetic, gravitational and magnetic energy gives the ideal magnetic energy flux $\mathbf{u}(E+p+B^2/(2{\mu }_0))-\mathbf{B}(\mathbf{u}\cdot \mathbf{B})/{\mu }_0$. This is a mechanical barotropic energy identity; it does not say that an isothermal gas is thermally isolated.

For fixed cloud shape, gravitational energy scales as $-C_{g}G{M}^2/L$ and magnetic energy as $C_B{\Phi }^2/({\mu }_{0}L)$. The magnetic-criticality threshold for negligible pressure is $M/\Phi >\sqrt{C_B/C_g}/\sqrt{{\mu }_{0}G}$. Ideal flux freezing keeps this ratio fixed during homologous collapse; geometry determines the numerical coefficient.

For smooth ideal magnetohydrodynamics of constant mass density, define $I_{\pm }={\int }_V|{\mathbf{z}}^{\pm }{|}^{2}dV$. Taking the dot product of the Elsässer variable equation with $2{\mathbf{z}}^{\pm }$ givesTherefore each integral is constant when its outward boundary flux vanishes, for example when both velocity and magnetic field are tangent to the fixed boundary, or for periodic boundary conditions. The kinetic energy plus magnetic energy is $\rho (I_{+}+I_{-})/4$, while the cross-helicity is $\sqrt{{\mu }_0\rho }(I_{+}-I_{-})/4$.

For a fixed Newtonian gravitational potential, total energy density is $\mathcal{E}=\rho (u^2/2+\Phi )+p/(\gamma -1)+B^2/(2{\mu }_0)$. Its flux combines advected kinetic/potential energy, specific enthalpy and the Poynting vector. Dotting the ideal magnetohydrodynamic momentum equation with velocity, adding the adiabatic internal energy equation and the magnetic energy equation cancels the magnetic work. The pressure terms combine into $-\nabla \cdot (p\mathbf{u})$, proving the conservation law. A time-dependent imposed potential instead contributes $\rho {\partial }_t\Phi$.

In vacuum-permeability magnetohydrodynamics, the magnetic energy in a volume is $E_B=\int B^2/(2{\mu }_0)dV$. Its volume density equals the magnetic pressure.

For a solenoidal vector field carrying current only in a conductor contained in a sphere of radius $R$, match the exterior field to a decaying potential field and include exterior magnetic energy. The smallest free-decay eigenvalue of the enclosing sphere is ${\pi }^2/R^2$, from its dipolar poloidal mode. The corresponding variational estimate gives the displayed bound. Continuity and the usual insulating-interface magnetic conditions are part of the admissible field class. This bound is the diffusive ingredient in Backus' necessary condition for dynamo action.

Under the energy-closed boundary assumptions of Backus' necessary condition for dynamo action, the magnetic energy equation gives ${\dot{E}}_B\le 2s(t)E_B$ after discarding nonnegative resistive dissipation. Integrating gives an upper bound by the time average of the largest spatial eigenvalue of the rate-of-strain tensor. Thus the field-amplitude exponent is at most its space-time supremum. The exponent of squared energy is twice the field-amplitude exponent.

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

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

Let $e_{ij}=({\partial }_i{u}_j+{\partial }_j{u}_i)/2$ be the rate-of-strain tensor of an incompressible flow, and let $s(t)$ be the supremum over the conductor of its largest eigenvalue. State Backus' necessary condition for dynamo action with its magnetic boundary conditions: an isolated bounded conductor of uniform positive magnetic diffusivity $\eta$, surrounded by an electrical insulator with a decaying potential exterior field, and no imposed magnetic field or boundary energy input. For definiteness take a no-slip boundary condition on the fluid, which eliminates the stretching surface term. If the conductor lies within a sphere of radius $R$ and $S={\sup }_{t}s(t)$, a necessary condition for a nondecaying dynamo isThe constant is the free-decay spectral bound for an insulating exterior, not a universal constant for every magnetic boundary condition. The condition is necessary, not sufficient, and involves maximum stretching rather than an rms velocity.

To see both the condition and the growth-rate bound, include exterior magnetic energy:This follows from the resistive induction equation and integration by parts, with the stated boundary assumptions. The magnetic free-decay spectral bound is ${\int }_D|\nabla \times \mathbf{B}{|}^{2}dV\ge ({\pi }^2/R^2){\int }_{{\mathbb{R}}^3}|\mathbf{B}{|}^{2}dV$. For a sphere its lowest mode is the dipolar poloidal free-decay mode; enclosing a smaller conductor gives the same valid lower bound. Since $B_i{e}_{ij}{B}_j\le s(t)|\mathbf{B}{|}^2$, we obtainIntegrating this differential inequality gives decay whenever $S<\eta {\pi }^2/R^2$. More generally the exponential rate of the field norm, rather than of its squared energy, satisfiesSimply discarding the nonnegative resistive dissipation already proves the requested maximum-strain bound. The energy exponent is twice the field-amplitude exponent.

For the alpha-Omega dynamo model, write $a=|{\alpha }_0|$, $w=|\Omega |$, $k\ge 0$ and $d=\eta k^2$. Direct differentiation givesFor $a,kw>0$, use the weighted energy estimate for two coupled modes and form the positive weighted norm $N=kw|A{|}^2+a|B{|}^2$. The inequality $N\ge 2\sqrt{akw}|A||B|$ givesFor each fixed $k$ and model parameters, this norm is equivalent to the amplitude norm; its square-root exponential rate is therefore bounded by $h(k)=\sqrt{awk}-\eta k^2$, uniformly over all admissible $f(t)$. Maximizing over $k$ givesThe exponent is also achievable in order of magnitude. Choose the admissible constant $f=1$. The growing eigenvalue of the two-component system has real part $-\eta k^2+\sqrt{awk/2}$. Its maximum occurs at $k=(aw/(32{\eta }^2){)}^{1/3}$ and equals $32^{-10/3}(a^2{w}^2/\eta {)}^{1/3}$. Thus the bounded-modulation alpha-Omega growth estimate has the scalingThis means the maximum over allowed modulations and wavenumbers, not that every bounded modulation grows; $f=0$ supplies no regenerating alpha coupling.

The Omega effect rapidly makes toroidal field from poloidal field, but exponential dynamo action also requires the slower alpha effect to regenerate the poloidal component. The coupled amplification rate is of order $\sqrt{a\Omega k}$ rather than $\Omega$; shortening the wavelength to increase it also increases magnetic diffusion as $\eta k^2$. Their optimal balance gives $k\sim {\Omega }^{1/3}$ and growth $\sim {\Omega }^{2/3}$. The Backus' necessary condition for dynamo action estimate controls stretching alone and does not incorporate this regeneration bottleneck. A shear without regeneration can give transient amplification but not this sustained exponential feedback.

Use Faraday's law ${\partial }_t\mathbf{B}=-\nabla \times \mathbf{E}$ and the solenoidal magnetic-field constraint $\nabla \cdot \mathbf{B}=0$. In the nonrelativistic, single-fluid approximation, neglect Hall and other nonideal electromotive terms and use the moving-conductor moving-conductor Ohm law, $\mathbf{E}+\mathbf{u}\times \mathbf{B}=\mathbf{J}/\sigma$. Infinite electrical conductivity with finite current gives $\mathbf{E}=-\mathbf{u}\times \mathbf{B}$, henceThis is the ideal magnetohydrodynamic induction equation. For comparison, neglecting displacement current in Ampère-Maxwell equation gives $\mathbf{J}=\nabla \times \mathbf{B}/{\mu }_0$; with uniform finite electrical conductivity it produces magnetic diffusion ${\eta }_m{\nabla }^2\mathbf{B}$, where the magnetic diffusivity is ${\eta }_m=1/({\mu }_0\sigma )$. The ideal approximation requires a large magnetic Reynolds number $UL/{\eta }_m$. Dropping displacement current is useful for this finite-conductivity comparison, but Faraday's law and the ideal Ohm relation already suffice for the ideal induction equation.

Expanding the curl and using the solenoidal magnetic-field constraint gives the material derivative formCombine this with mass conservation, $D\rho /Dt=-\rho \nabla \cdot \mathbf{u}$, to obtainNow parametrize a material curve by a fixed label $a$: $\mathbf{X}(a,t)$ obeys ${\partial }_t\mathbf{X}=\mathbf{u}(\mathbf{X},t)$. Differentiating with respect to $a$ shows that its tangent $\mathbf{\ell }={\partial }_a\mathbf{X}$ evolves by $D\mathbf{\ell }/Dt=(\mathbf{\ell }\cdot \nabla )\mathbf{u}$. This is exactly the same linear ordinary differential equation as for $\mathbf{B}/\rho$. Initially parallel tangents remain parallel by uniqueness, with a label-dependent proportionality factor constant along each particle trajectory. Thus magnetic field lines are transported as material curves, wherever the field and fluid flow are smooth and the field is nonzero. This is the field-line form of magnetic flux freezing.

For the flux statement, take a material surface $\mathbf{X}(a,b,t)$ and let ${\mathbf{t}}_a={\partial }_a\mathbf{X}$, ${\mathbf{t}}_b={\partial }_b\mathbf{X}$. Both tangents obey $D{t}_i/Dt=({\partial }_j{u}_i)t_j$. Its oriented material surface element is $\mathbf{A}={\mathbf{t}}_a\times {\mathbf{t}}_{b}dadb$. Differentiating the cross product, rather than assuming its transport rule, givesEquivalently, $D\mathbf{A}/Dt=(\nabla \cdot \mathbf{u})\mathbf{A}-(\nabla \mathbf{u}{)}^T\mathbf{A}$, with $(\nabla \mathbf{u}{)}_{ij}={\partial }_j{u}_i$. Contracting this derived rule with the induction equation gives a pointwise cancellation:Integrating over the fixed material labels therefore proves conservation of flux through an open material surface:The surface need not be closed; its boundary is carried with the fluid. The result follows from material transport, rather than from the zero flux through a closed surface.

For a homologously shrinking cloud, write $V\sim L^3$ and keep its shape factors fixed. Conserved mass gives $\rho \sim M/L^3$; conserved magnetic flux gives $B\sim \Phi /L^2$. Thus the gravitational and magnetic energies scale aswhere $C_g,C_B>0$ are dimensionless geometry factors. Both grow in magnitude as $L^{-1}$, so collapse cannot reduce magnetic support relative to gravity while the mass-to-flux ratio is frozen. With negligible gas pressure, contraction lowers the combined potential energy only when its coefficient of $L^{-1}$ is negative. Consequently a necessary critical mass-to-flux ratio condition isHere $\Phi$ denotes the magnitude of the conserved threading flux. The numerical coefficient depends on geometry and boundary conditions; the scaling argument does not determine it or make the condition sufficient in the presence of other support.

For adiabatic pressure support during gravitational collapse, $p\propto {\rho }^{\gamma }\propto L^{-3\gamma }$. The pressure-support scale is $pV\propto L^{3-3\gamma }$, so relative to either gravity or magnetic energy,Pressure becomes more important as $L$ decreases if $\gamma >4/3$, equally important in scaling if $\gamma =4/3$, and less important if $\gamma <4/3$. In particular, a monatomic perfect gas with $\gamma =5/3$ becomes increasingly pressure supported. For isothermal pressure support during gravitational collapse, the isothermal equation of state gives $p\propto \rho \propto L^{-3}$, so $pV$ is constant and $pV/|W_g|\propto L$: isothermal pressure becomes less important during collapse. The same comparisons hold against magnetic support because its energy has the same $L^{-1}$ scaling as gravity.

Let the constant mass density be $\rho >0$ and put ${\mathbf{v}}_a=\mathbf{B}/\sqrt{{\mu }_0\rho }$, the Alfvén velocity. Since both the velocity and the magnetic field have zero divergence, the magnetic tension and magnetic pressure decomposition of the Lorentz force givesThe Newtonian gravitational potential $\Phi$ remains in $\psi$; uniform mass density does not justify dropping a prescribed gravitational acceleration. Adding and subtracting the two equations proves the Elsässer variable equationsHere ${\mathbf{z}}^{\pm }=\mathbf{u}\pm {\mathbf{v}}_a$; each Elsässer variable is transported by the other.

Define $I_{\pm }={\int }_V|{\mathbf{z}}^{\pm }{|}^{2}dV$. Taking the dot product with $2{\mathbf{z}}^{\pm }$ gives the Elsässer energy invariant balanceThe divergence theorem proves $d{I}_{\pm }/dt=0$ whenever the net boundary flux vanishes. In particular, $\mathbf{u}\cdot \mathbf{n}=\mathbf{B}\cdot \mathbf{n}=0$ makes ${\mathbf{z}}^{\pm }\cdot \mathbf{n}=0$, so both fluxes vanish pointwise. Periodic boundary conditions or decay at infinity are also sufficient. These boundary conditions matter: fixed volume alone gives no conservation.

The kinetic energy plus magnetic energy isExpanding $|\mathbf{u}\pm {\mathbf{v}}_a{|}^2$ and using the cross-helicity definition yieldsBoth energy and cross-helicity therefore follow from the two Elsässer energy invariants. The quantity $E$ here is exactly the requested kinetic energy plus magnetic energy, without an additional gravitational term.