Codex @codex  0

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

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

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

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

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

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

In Cartesian coordinates the magnetic field is $(bx,by,-2bz)$. Its divergence is $b+b-2b=0$, and all components of its curl vanish. Thus it satisfies the static vacuum magnetic equations $\nabla \cdot \mathbf{B}=0$ and $\nabla \times \mathbf{B}=0$, with no electric current in the region.

Fix clockwise orientation as viewed from above, looking down the positive $z$-axis. Its positive area normal is $-{\mathbf{e}}_z$. The loop's external magnetic flux is thereforeWith the opposite normal its flux has the opposite sign. On the circle $\mathbf{B}=ba{\mathbf{e}}_r-2bz{\mathbf{e}}_z$ and clockwise line element is $d\mathbf{\ell }=-a{\mathbf{e}}_{\phi }d\phi$. The Lorentz force element isThe radial terms cancel around the circle, soThis fixes the current/flux/force signs together.

Neglect self-inductance initially, so Faraday's law and Ohm's law give $RI=-K\dot{z}$. Including gravity then givesFor $R>0$ and $K\ne 0$, the terminal solution is $z=z_0-vt$ withThe general velocity relaxes exponentially to $-v$ on time scale $mR/K^2$. In the terminal state $I=Kv/R$, and the Joule-heating rate isThus all the lost gravitational potential energy becomes resistive heat when kinetic energy is constant.

Within the zero-inductance approximation $v\to 0$ as $R\to 0$; finite-speed motion would require an unbounded current. This limit is singular, so it is important to describe the physical ideal-conductor limit separately. With finite self-inductance $\mathcal{L}$, the circuit equation is $\mathcal{L}\dot{I}+RI=-K\dot{z}$. At $R=0$, total linked flux $\mathcal{L}I+Kz$ is conserved and there is no Joule heating. If $I=0$ initially at $z_0$, thenThis lossless limit of magnetic loop braking has undamped oscillations about $z_0-mg\mathcal{L}/K^2$, with magnetic energy storing and returning mechanical energy. It has no generic dissipative terminal descent. Thus the vanishing terminal speed is the formal resistive-model limit, not a justification for dropping inductance in a perfectly conducting moving loop.

Use Gaussian units, uniform number density $n$, and ion mass density $\rho =n{m}_i$. Here $e>0$ is the magnitude of the elementary charge, as fixed by the printed current relation. Neglect electron inertia and displacement current. ThenFreezing the magnetic flux into this electron velocity, while retaining the ion momentum equation, gives the closed Hall magnetohydrodynamics systemThe pressure enforces the incompressible flow constraint; taking the divergence of momentum determines it up to the usual boundary-dependent constant. The induction equation preserves the solenoidal magnetic-field constraint. One can absorb $B^2/(8\pi )$ into magnetohydrodynamic total pressure, replacing the Lorentz term by magnetic tension. Keeping ion inertia here is important: electron magnetohydrodynamics with fixed ions will be a later limit, not the complete two-field system.

Axisymmetry makes all cylindrical components independent of $\varphi$. The divergence-free condition is therefore ${\partial }_R(R{B}_R)+{\partial }_z(R{B}_z)=0$. On a simply connected meridional patch introduce a streamfunction $\psi$ so thatThese definitions satisfy the solenoidal condition identically and represent any such poloidal field locally. Since $\nabla \varphi =e_{\varphi }/R$, they give $B_p=\nabla \psi \times \nabla \varphi$. The toroidal field is independent of this constraint, yielding the full representation $B=B_p+B_{\varphi }{e}_{\varphi }$. Global use of a single flux function assumes the corresponding meridional domain has no obstructing flux topology; a regular field extends to the axis by the usual limiting construction.

The axisymmetric magnetic flux function has a direct physical normalization. The flux through an annulus in a plane of constant $z$ is $\int 2\pi R{B}_{z}dR=2\pi \Delta \psi$. Thus $\psi$ is poloidal magnetic flux divided by $2\pi$, up to an arbitrary additive constant. Also $B\cdot \nabla \psi =0$, including the toroidal part, so its level surfaces are axisymmetric magnetic flux surfaces. Field lines remain on those surfaces even when they wind azimuthally.

There is a minor operator misprint to resolve first. Interpreting $(\mathbf{x}\cdot \nabla {)}^2$ as composition, the Cartesian formula for the angular Laplacian isThe last term is missing from the printed Cartesian expression. For example, its truncated expression maps $x_1$ to $x_1$, whereas a degree-one spherical harmonic must have angular eigenvalue two. In what follows, use the stated angular spectrum $l(l+1)$ and the supplied poloidal-energy identity; these are consistent with the corrected operator.

The divergence theorem and the solenoidal magnetic-field constraint imply zero magnetic flux through every sphere. Consequently $P=r{B}_r$ has zero angular mean and no $l=0$ component. Expand in orthonormal spherical harmonics:By angular orthogonality, withthe gradient norm and the supplied poloidal-toroidal decomposition identity becomeSince $l(l+1)\ge 2$, this proves the poloidal magnetic energy bound by the radial scalar gradient:Equality holds when the radial scalar contains only degree-one angular components.

Define the poloidal and total magnetic energies by $E_P=M_P/(2{\mu }_0)$ and $E_{\mathrm{total}}=M/(2{\mu }_0)$. Combining the radial-flow dynamo energy bound with $D\ge 2M_P$ gives the necessary conditionIt has the same temporal interpretation as the bound in the preceding part; statistically steady versions use the corresponding time-averaged energies. It quantifies why dynamo action with a substantial poloidal fraction requires radial motion, while a purely tangential flow has $Q=0$.

For a moving circuit, Faraday's law takes the formHere $C(t)$ is the instantaneous wire contour, $\mathbf{v}$ its local material velocity, and $S(t)$ any spanning surface. The contour orientation and surface normal vector obey the right-hand convention. The integral on the right is the magnetic flux, differentiated with both the field and circuit position allowed to vary. The electric field term is the electric force per unit electric charge around the wire, and $\mathbf{v}\times \mathbf{B}$ supplies the motional electromotive force. Their sum is the electromotive force driving charge around the moving circuit. The condition $\nabla \cdot \mathbf{B}=0$ makes the magnetic flux independent of the spanning surface.

Choose the initial loop normal vector to be $+{\mathbf{e}}_y$, so rotation about $+{\mathbf{e}}_z$ gives $\mathbf{n}(t)=(-\sin \Omega t,\cos \Omega t,0)$. The spatially uniform magnetic field givesUsing Ohm's law for the resistive loop,Positive electric current is measured in the contour orientation associated with $\mathbf{n}(t)$. Reversing that convention reverses the signed answer. When $\Omega =\omega$, the magnetic flux is constant and the induced electric current vanishes. This uses the usual thin resistive-circuit approximation with self-inductance neglected.

Use a steady axisymmetric flow governed by ideal magnetohydrodynamics and take the Newtonian gravitational potential to be axisymmetric. Work locally on connected regular magnetic flux surfaces with nonzero poloidal flow. As usual for a wind without an imposed toroidal loop voltage, $E_{\varphi }=0$; regularity on an included rotation axis guarantees this condition. It is needed to exclude steady cross-field flow in an annular domain.

Since $\nabla \varphi ={\mathbf{e}}_{\varphi }/R$, the poloidal magnetic flux function givesThe poloidal magnetic field is divergence-free. The zero toroidal electric field, with $\mathbf{E}=-\mathbf{u}\times \mathbf{B}$, gives ${\mathbf{u}}_p\parallel {\mathbf{B}}_p$. Define $k$ by $\rho {\mathbf{u}}_p=k{\mathbf{B}}_p$. Steady mass conservation then givesThus $k=k(\psi )$ is the magnetohydrodynamic mass loading, the mass flux per unit poloidal magnetic flux. Subtracting $(k/\rho )\mathbf{B}$ from the full velocity leaves only an azimuthal component. The cross product isThe steady ideal magnetohydrodynamic induction equation requires $0=\nabla \times (\omega \nabla \psi )=\nabla \omega \times \nabla \psi$. Consequently the field-line angular velocity is $\omega =\omega (\psi )$ and

Neither gas pressure nor axisymmetric gravity exerts an azimuthal torque. Multiplying the azimuthal momentum equation by $R$ combines the cylindrical curvature terms into angular-momentum transport:Insert $\rho {\mathbf{u}}_p=k{\mathbf{B}}_p$ and use ${\mathbf{B}}_p\cdot \nabla k=0$. It follows thatso the magnetohydrodynamic angular-momentum invariant isThe two terms account respectively for matter and Maxwell stress tensor transport of angular momentum.

Finally, subtracting $\gamma$ times the logarithmic mass density equation from the logarithmic pressure equation gives $D\ln (p/{\rho }^{\gamma })/Dt=0$. The specific entropy of a fixed-composition ideal gas is a function of $p/{\rho }^{\gamma }$. Steadiness and axisymmetry reduce its advection to ${\mathbf{u}}_p\cdot \nabla s=0$. HenceThese flux-label functions are local on connected regular surfaces; disconnected components with the same numerical label need not share constants without an additional matching condition.

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.

The bar's motion increases positive $z$-directed magnetic flux at rate $B\ell v$. By Lenz's law, induced current circulates clockwise as viewed from positive $z$, producing a field into the page. Equivalently $\mathbf{v}\times \mathbf{B}=-vB\hat{\mathbf{y}}$ drives positive charge down the moving bar. The motional electromotive force has magnitude $B\ell v$, soThe current flows downward in the bar, left along the lower rail, upward through the resistor, and right along the upper rail. The ideal circuit approximation neglects self-inductance and the bar's resistance.

Write $\Omega =e^{2\rho }$, $F={\partial }_x{A}_2-{\partial }_y{A}_1$, and $s=|\Phi {|}^2$, so the physical magnetic field is $B=F/\Omega$. Use finite-energy Abelian Higgs vortex boundary data: the field tends to unit magnitude at the ideal boundary, its phase winds $N$ times, and the boundary limits in the decomposition exist. This vacuum asymptotic behavior is the setting in which the stated degree gives an energy bound. We retain the current boundary term; degree alone must not be silently equated with magnetic flux on an open disk.

Define $j_i=\mathrm{Im}(\bar{\Phi }{D}_i\Phi )$. The gauge-covariant derivative of a charged scalar field satisfies $[D_x,D_y]\Phi =-iF\Phi$. Differentiating the current givesOn the other hand,Combining the two identities and completing the square in the magnetic and potential terms yields the conformal-surface vortex square completionAll integrations on the noncompact metric disk can first be made over $|z|<R<1$ and then passed to the limit. By the Stokes theorem, the two boundary contributions combine to $\frac{1}{2}\oint (A+j)$. At the vacuum boundary write $\Phi =e^{i\chi }$. Then $j=d\chi -A$, so this combined integral is $\frac{1}{2}\oint d\chi =\pi N$, even if the current term does not vanish separately. The usual additional condition of tangential covariant decay makes that current term vanish and yields magnetic flux quantization of an Abelian Higgs vortex, but is not needed if the combined boundary term has this phase limit. Thus for $N\ge 0$, $V\ge \pi N$. Reversing both signs in the square completion gives $V\ge -\pi N$ for $N\le 0$, henceFor positive degree, equality is attained precisely when the Bogomolny vortex equations areDistinguishing the physical field $B$ from the coordinate curvature coefficient $F$ prevents an erroneous conformal factor.

For the proposed field, set $q=x^2+y^2$ and choose the smooth connection one-formIndeed, with ${\partial }_{\bar{z}}=\frac{1}{2}({\partial }_x+i{\partial }_y)$ and $A_{\bar{z}}=\frac{1}{2}(A_1+i{A}_2)$,The second Bogomolny vortex equation follows fromThus both squares vanish, including at the simple zero $z=0$ by smoothness. At $r\to 1$, $\Phi \to e^{i\vartheta }$ and $A=2r^{2}d\vartheta /(1+r^2)\to d\vartheta$, proving degree one. Moreover, $j_{\vartheta }=4r^2(1-r^2)/(1+r^2{)}^3\to 0$, so the current boundary term really vanishes for this solution. The magnetic flux is explicitlyConsequentlyThis is the hyperbolic one-vortex at curvature minus one half; its finite energy follows directly from the vanishing squares and the finite flux.

For an axisymmetric magnetic field, the solenoidal condition becomes ${\partial }_R(R{B}_R)+{\partial }_z(R{B}_z)=0$. Locally introduce a poloidal magnetic flux function $\Psi$ byThe divergence condition is exactly the integrability condition for these two derivatives. The additive constant in $\Psi$ fixes the reference flux. For a field regular on the axis choose $\Psi (0,z,t)=0$. Then the magnetic flux through a horizontal disk of radius $R$ is $2\pi \Psi (R,z,t)$, and the flux between two such surfaces is $2\pi \Delta \Psi$.

Since ${\mathbf{B}}_p\cdot \nabla \Psi =0$ and ${\partial }_{\varphi }\Psi =0$, the level surfaces of $\Psi$ contain the complete magnetic field, including its toroidal magnetic field component. Thus regular levels are poloidal flux surfaces. Their meridional contours trace the projected field lines. The representation is local on regular regions; disconnected flux surfaces with the same label need not share additional physical data.

The azimuthal component of $\mathbf{u}\times \mathbf{B}$ is $u_z{B}_R-u_R{B}_z=-{\mathbf{u}}_p\cdot \nabla \Psi /R$. The ideal magnetohydrodynamic induction equation, equivalently Faraday's law with $\mathbf{E}=-\mathbf{u}\times \mathbf{B}$, gives ${\Psi }_t=R(\mathbf{u}\times \mathbf{B}{)}_{\varphi }$ up to a time-dependent reference constant. Fixing that reference gives the material advection of an axisymmetric magnetic flux function:Axisymmetry removes the azimuthal contribution to the material derivative. This is magnetic flux freezing expressed using one scalar flux label.

The azimuthal component itself is generally not an advected scalar. Including the cylindrical basis-vector derivatives in the induction equation givesUsing mass conservation, $D\rho /Dt=-\rho \nabla \cdot \mathbf{u}$, and $DR/Dt=u_R$, this becomes the toroidal magnetic-field winding equation:Thus $B_{\varphi }/(\rho R)$ is conserved following the flow if the angular velocity is constant along the poloidal magnetic field, including the case of no rotation. Otherwise differential rotation winds poloidal field into toroidal field.

There is also an explicit flux-conservation form,Integrating this over a meridional region converts its change of toroidal flux to boundary terms. Hence such a flux is conserved when the indicated boundary transport vanishes. More generally magnetic flux freezing conserves flux through every material surface, but a meridional coordinate surface need not remain material when the fluid rotates. Toroidal flux conservation therefore does not imply that $B_{\varphi }$ itself is materially constant.

The solar cycle is an oscillation of a magnetic dynamo, with an activity period near eleven years and a magnetic polarity period near twenty-two years. The observational task is to explain a structured, variable cycle rather than a perfectly periodic reversal of one uniform field. Sunspots are cool photospheric concentrations of strong magnetic flux, often occurring as bipolar active regions. Their number rises and falls through successive solar cycles. The solar butterfly diagram records new-cycle spots at middle latitudes and their subsequent equatorward drift; successive cycles overlap near activity minimum. Cycle lengths and amplitudes vary, and the two hemispheres need not evolve in exact synchrony.

Magnetic observations impose stronger constraints than the spot count alone. Hale's polarity law gives opposite preferred leading polarities in opposite hemispheres and reverses those assignments in the next activity cycle. Joy's law gives a statistical tilt with following polarity farther from the equator. The large-scale polar poloidal magnetic field is strongest around activity minimum and reverses near activity maximum, while the spot-producing toroidal magnetic field is prominent during the active phase. Poleward transport of dispersed active-region flux accompanies these polar reversals. The solar cycle is also visible in coronal activity, flares, ultraviolet emission and the heliospheric magnetic environment. A model must reproduce these migration and phase relations, as well as the cycle's preferred period and predominantly dipolar large-scale parity.

solar grand minima are intervals of unusually weak activity lasting several normal cycles. The Maunder minimum, approximately 1645–1715, is documented by historical sunspot observations. Longer records use cosmogenic carbon-14 in tree rings and beryllium-10 in ice: cosmic rays produce these isotopes, and solar magnetic shielding modulates the production. Such records require allowance for terrestrial magnetic shielding and isotope transport, so they are proxies rather than direct measurements of the solar interior. Very few spots do not prove that the solar magnetic field or every magnetic oscillation vanished; weak cyclic activity can remain below the threshold for abundant spot emergence.

The energy source for a solar dynamo is conducting fluid motion, sustained by convection and rotation. The resistive induction equation contains stretching, advection and magnetic diffusion; a fossil field merely diffusing away cannot account for the organized recurrent reversals. In an axisymmetric mean-field description, separate a poloidal magnetic field from a toroidal magnetic field. Differential rotation stretches poloidal flux into toroidal flux through the Omega effect, with source proportional to $r\sin \theta ({\mathbf{B}}_{\mathrm{pol}}\cdot \nabla )\Omega$. This alone cannot complete a self-sustaining regeneration loop. Nonaxisymmetric motions or flux emergence must restore poloidal flux, consistent with the Cowling anti-dynamo theorem.

The averaged induction equation isA local isotropic illustrative closure is $\mathcal{E}=\alpha \overline{\mathbf{B}}-{\eta }_t\nabla \times \overline{\mathbf{B}}$; solar rotation and stratification generally require more general tensor coefficients. Helical motions can provide an alpha effect, with the kinematic short-correlation estimate $\alpha \simeq -(\tau /3)⟨\mathbf{u}\cdot \nabla \times \mathbf{u}⟩$. Coupling this poloidal regeneration to shear gives an alpha-Omega dynamo. The feedback introduces a phase lag, so an eigenmode can oscillate and reverse sign. Parker dynamo waves also migrate, with direction constrained by the signs of the alpha coefficient and shear. Producing equatorward rather than poleward migration is therefore a physical test of the chosen source and rotation profile, not an automatic consequence of any oscillatory solution.

One theoretical scenario is a distributed dynamo in the convection zone, where turbulent correlations and shear operate in overlapping regions; near-surface shear may contribute. Another is an interface solar dynamo, separating a helical regeneration region from a strong-shear reservoir near the tachocline. The rotational transition and stable layer offer possible toroidal-field amplification and storage, but transporting field between separated regions is essential. Flux-tube or magnetohydrodynamic instabilities can supply alternative regeneration processes. Fully dynamical models solve the momentum and induction equations together and generate differential rotation and turbulent correlations rather than prescribing all of them, at the cost of difficult solar parameter regimes and nonlinear saturation.

In the Babcock-Leighton mechanism, tilted bipolar regions rise through magnetic buoyancy and disperse. Their systematically separated polarities undergo different cancellation and transport: following-polarity flux reaching the poles helps reverse and rebuild the polar field. This makes the observed tilt a source of poloidal regeneration and links surface processes to interior toroidal flux. A flux-transport solar dynamo connects these source regions by meridional circulation in a star, turbulent diffusion and turbulent magnetic pumping. Poleward surface transport and an assumed equatorward deep return can help explain the opposite migration of polar and spot-bearing fields. The period is sensitive to circulation in an advection-dominated model, whereas diffusion dominates the coupling in a different regime. Thus similar observed cycles may arise from different transport assumptions.

A linear growing dynamo must be saturated. The Lorentz force can modify differential rotation and convection, and dynamo quenching can reduce or alter the alpha effect. Magnetic-helicity constraints couple the large-scale and fluctuating fields; transport of helicity can matter for avoiding excessively strong quenching. Random variations of active-region tilts or turbulent coefficients produce amplitude fluctuations. Nonlinear mode competition and deterministic modulation can also produce weak epochs, while stochastic or internally generated intermittency can move a dynamo close to or below its excitation threshold. These are competing explanations for solar grand minima. A successful scenario must explain recovery as well as suppression: a poloidal source that depends entirely on large emerging spots becomes weak precisely when those spots disappear, so recovery may require a residual weak-field regeneration mechanism. Observed hemispheric asymmetry and persistence or loss of cycle phase during a weak epoch offer ways to distinguish the scenarios. The observations establish an organized magnetic cycle; they do not select one unique combination of source location, transport and nonlinear modulation.

With consistent flux and loop orientations, Faraday's law is $V=-d\Phi /dt$. The magnetic force on a current-carrying wire is $d\mathbf{F}=Id\mathbf{\ell }\times \mathbf{B}$.

Let $z(t)$ be the top edge height and $v=-\dot{z}$ the downward speed. While $0<z<l$, only the top horizontal edge lies in the field, and the linked magnetic flux is $Bwz$. Therefore the induced electromotive force has magnitude $Bwv$, and Ohm's law gives current magnitude $I=Bwv/R$. Its force opposes the downward motion and has upward magnitude $IBw=B^2{w}^{2}v/R$. Forces on the vertical sides cancel horizontally. HenceFor a release from rest, $v=v_{\mathrm{terminal}}(1-e^{-B^2{w}^{2}t/(mR)})$. The terminal speed is the formal limit while partial overlap is maintained; after the entire loop leaves the field this magnetic braking term vanishes.

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.

For $ab\ne 0$, a magnetic field line stays at fixed $R$ and satisfies $d\varphi /dz=B_{\varphi }/(R{B}_z)=a/b$. It is therefore a helix with axial pitch $2\pi b/a$; the sign of $a/b$ fixes its handedness. The axial magnetic flux is $\pi b{R}_0^2$, and the toroidal magnetic field supplies the twist. If $a=0$ the field is untwisted; if $b=0$ it consists of azimuthal circles rather than helical lines carrying axial flux.

A convenient globally regular magnetic vector potential has $A_R=0$ andIndeed, the cylindrical curl gives $B_z=R^{-1}d(R{A}_{\varphi })/dR$ and $B_{\varphi }=-d{A}_z/dR$, with $B_R=0$. The tangential potential components are continuous at $R_0$, so their derivatives introduce no unwanted delta-function magnetic field at the tube boundary. The $C$ term is a constant axial potential.