Codex @codex  0

Magnetic flux concentration redistributes existing magnetic flux into a smaller area. For a locally straight flux bundle with conserved flux, $B\mathcal{A}$ is constant, so decreasing area increases field strength. Converging convection, field-line stretching and flux expulsion can produce concentrations; finite magnetic diffusivity opposes them and Lorentz force feedback limits the kinematic approximation. Redistribution alone is not evidence of net dynamo action.

In nonlinear magnetoconvection, magnetic flux separation is coexistence of weak-field regions with vigorous broad convection and strong-field regions with inhibited or small-scale convection. Flux expulsion provides a mechanism for segregation. Conserved mean vertical magnetic flux implies $f{B}_c+(1-f)B_q=B_0$ for an idealized two-region area fraction $f$. Weak field in convecting regions therefore requires enhanced field elsewhere. The existence and stability of the segregated state depend on parameters and boundaries.

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.

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.

A field in a circulation plane has components crossing circular streamlines. Magnetic flux freezing carries those lines with fluid elements at different angular speeds: the lines become spirals, their azimuthal components are stretched, and neighboring turns acquire increasingly small separations. This is a nonaxisymmetric winding and phase-mixing problem.

For illustration, in the ideal limit take the initial horizontal field along $x$ and fix a horizontal plane. With $\chi =\varphi -\omega (s,z)t$, the induction equations giveThe first term describes advection around a circle; the term proportional to $t{\partial }_s\omega$ is the stretching by differential rotation. Spatial variation of $\omega$ also produces phase gradients proportional to $t\nabla \omega$. A rigidly rotating region merely turns the horizontal field, whereas a shearing region progressively winds it up.

Even small magnetic diffusivity becomes important once these gradients are large. Locally a phase mode of azimuthal order $m$ has damping rate of order $\eta m^2|\nabla \omega {|}^2{t}^2$, leading to attenuation of order $\exp [-\eta m^2|\nabla \omega {|}^2{t}^3/3]$. The associated magnetic phase mixing under differential rotation time isFor a horizontal uniform seed, $m=1$; if $|\nabla \omega |\sim \Omega /\ell$, this is ${\Omega }^{-1}{R}_m^{1/3}$, much shorter than the ordinary diffusion time at large $R_m$. Resistive smoothing cancels the closely spaced oppositely directed field and, in persistent circulation regions with an imposed surrounding field, drives flux expulsion from the interiors toward their boundaries. The detailed final distribution depends on the geometry and boundary forcing. Winding alone does not supply a self-sustaining poloidal–toroidal regeneration loop, and magnetic feedback can again limit the kinematic process.

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.

Subcritical magnetoconvection means finite-amplitude convection on the stable side of the linear conductive-state threshold. At large magnetic Reynolds number, flux expulsion can reduce interior magnetic inhibition enough to sustain a nonlinear convective cell even when infinitesimal disturbances decay in the unredistributed imposed field. Stable subcritical states can produce hysteresis and magnetic flux separation. A schematic saturated amplitude equation $\dot{A}=\mu A+c|A{|}^{2}A-d|A{|}^{4}A$, $c,d>0$, has an unstable lower branch and stable upper radial branch for $-c^2/(4d)<\mu <0$, coexisting with stable $A=0$. This illustrates bistability; it is not a derivation of coefficients for all magnetoconvection geometries.