Codex @codex  0

Write $u=u_x$ and $B=B_z$. In this geometry the Lorentz force density is purely longitudinal and equals $-{\partial }_x(B^2/2{\mu }_0){\mathbf{e}}_x$. The ideal magnetohydrodynamic induction equation becomes $B_t+{\partial }_x(uB)=0$. Together with the continuity equation and adiabatic pressure evolution, the four primitive-variable equations areThey have the required quasilinear form ${\mathbf{U}}_t+A{\mathbf{U}}_x=0$ withHere the notation $u_x$ in derivatives means ${\partial }_{x}u$, whereas $u$ denotes the longitudinal velocity. Spatial dependence is one-dimensional; time dependence is allowed in the wave problem.

Put $c_f^2=\gamma p/\rho +B^2/({\mu }_0\rho )$, and assume $\rho >0$ and $c_f>0$. Expansion of the characteristic polynomial givesTwo independent eigenvectors for ${\lambda }_0$ can be chosen as $(1,0,0,0{)}^T$ and $(0,-B/{\mu }_0,0,1{)}^T$. For the two propagating branches, useful eigenvectors areThe repeated branch advects an entropy disturbance and a disturbance with unchanged total gas-plus-magnetic pressure. The $\pm$ branches are perpendicular fast magnetosonic waves, with speed $c_f$ relative to the fluid.

On a uniform background, a Fourier amplitude $\hat{\mathbf{U}}{e}^{i(kx-\omega t)}$ obeys $(kA-\omega I)\hat{\mathbf{U}}=0$. Thus the eigenvalues are the phase velocities $\omega /k$ for $k\ne 0$, and the four dispersion branches are $\omega =ku$ twice and $\omega =k(u\pm c_f)$. The reduced geometry retains neither an independent transverse velocity nor an additional propagating Alfvén wave.

To construct each finite-amplitude simple wave in magnetohydrodynamics, trace a state-space curve with tangent $r_{\pm }$, using $\rho$ as its parameter. Integrating the pressure and magnetic components givesThe constants $K$ and $b$ are the entropy and transverse magnetic flux freezing invariants along this state curve. HenceSince $A{\mathbf{U}}_{\rho }={\lambda }_{\pm }{\mathbf{U}}_{\rho }$, substituting this curve into the original system reduces it to ${\rho }_t+{\lambda }_{\pm }(\rho ){\rho }_x=0$. By the chain rule,Thus the characteristic speed itself satisfies the Inviscid Burgers equation. For initial speed $L(a)$, its smooth solution is ${\lambda }_{\pm }=L(a)$ along $x=a+tL(a)$, until the characteristic curves intersect. This derivation supplies genuine state curves for simple waves, not merely a linear dispersion relation.

Use the same length, thermal-time, velocity and temperature scales as in thermal convection. Write the imposed magnetic field as $B_0\hat{\mathbf{z}}$ and measure its perturbation in units of $B_0$. The physical linear magnetohydrodynamic momentum equation contains magnetic tension $(B_0^2/{\mu }_0){\partial }_z\mathbf{b}$; the gradient part of the Lorentz force density is absorbed into magnetohydrodynamic total pressure.

Divide the momentum equation by the viscous force scale ${\rho }_0\nu \kappa /d^3$. Its acceleration coefficient is $\kappa /\nu =1/\sigma$ and its magnetic coefficient isHere $Q$ is the Chandrasekhar number, comparing magnetic and viscous-diffusive effects, and $\zeta$ is the ratio of magnetic diffusivity to thermal diffusivity. The induction scale is $B_0\kappa /d^2$, so linearization givesComplete the linear system with $\dot{\theta }=w+{\nabla }^2\theta$, $\nabla \cdot \mathbf{u}=0$ and the solenoidal magnetic-field constraint $\nabla \cdot \mathbf{b}=0$. The omitted advective and perturbed-field tension terms are quadratic in the disturbances. Magnetic tension resists bending of the vertical field, while diffusion allows that bending to relax.

For a Cartesian change of basis represented by an orthogonal matrix $Q$, vectors transform as $E_i^{\prime }=Q_{ia}{E}_a$ and $B_i^{\prime }=Q_{ia}{B}_a$. Their squared norms are invariant and $Q_{ia}{Q}_{jb}{\delta }_{ab}={\delta }_{ij}$. Substitution into the given expression therefore givesthe transformation law of a Cartesian second-rank tensor. The magnetic field's extra axial sign under an improper physical reflection, if included, appears twice and cancels in its quadratic contribution.

Write the Maxwell stress tensor as $T=\mathbf{EE}+\mathbf{BB}-\frac{1}{2}(E^2+B^2)I$. Its divergence isThe identity $(\mathbf{A}\cdot \nabla )\mathbf{A}-\frac{1}{2}\nabla A^2=-\mathbf{A}\times (\nabla \times \mathbf{A})$ follows by contracting two Levi-Civita symbols, or directly by differentiating components. Using Maxwell's equations consequently yieldsHence the local conservation of electromagnetic momentum isThe two subtracted terms are the Lorentz force density, while $\mathbf{E}\times \mathbf{B}$ is the electromagnetic momentum density in the units used here.

Write $\mathbf{u}={\partial }_t\mathbf{\xi }$. In the uniform equilibrium, the linearized ideal magnetohydrodynamic equations reduce toIntegrating from the undisplaced reference state, and using the uniform background magnetic field, givesThese are the perturbations induced by the fluid displacement; independent time-independent changes to the reference state are excluded. The linear Lorentz force density is $(\nabla \times \delta \mathbf{B})\times \mathbf{B}/{\mu }_0$. Its magnetic pressure and magnetic tension parts giveFor a Fourier mode, put $s=\mathbf{k}\cdot \mathbf{\xi }$, $q=\mathbf{B}\cdot \mathbf{\xi }$ and $K=\mathbf{k}\cdot \mathbf{B}$. Then $\delta p=-i\gamma ps$ and $\delta \mathbf{B}=i(K\mathbf{\xi }-\mathbf{B}s)$. Substitution givesIf $s=q=0$, the fluid displacement is perpendicular to both the wave vector and the magnetic field. Only magnetic tension restores it, andThis is the Alfvén wave. The Alfvén velocity is the field-directed propagation vector: $\omega =\pm \mathbf{k}\cdot {\mathbf{v}}_a$. The signed phase velocity normal to the wavefront is $\pm v_a\cos \theta \hat{\mathbf{k}}$; the group velocity is $\pm {\mathbf{v}}_a$. The distinction matters for oblique propagation. For generic directions the Alfvén wave has one transverse polarization, with $\delta \rho =\delta p=0$.

To obtain the other magnetohydrodynamic waves, take the scalar products with $\mathbf{k}$ and $\mathbf{B}$:The determinant condition for a nonzero pair $(s,q)$ isThus the fast magnetosonic wave and slow magnetosonic wave haveTheir fluid displacements lie in the plane of $\mathbf{k}$ and $\mathbf{B}$ and are generally compressive. In the fast magnetosonic wave, gas and magnetic pressure provide the stronger restoring combination; in the slow magnetosonic wave, their perturbations oppose one another. The phase speeds depend on direction. They satisfy $v_{\mathrm{sl}}\le \min (c_s,v_a)$ and $v_{\mathrm{f}}\ge \max (c_s,v_a)$. For parallel propagation the two speeds are $c_s$ and $v_a$, with a degeneracy between a transverse branch and the Alfvén wave. For perpendicular propagation $v_{\mathrm{f}}^2=c_s^2+v_a^2$ and $v_{\mathrm{sl}}=0$; the latter is a nonpropagating limiting disturbance. These special directions require interpreting the polarizations by continuity rather than assuming three distinct nonzero frequencies.

In a steady flow the ideal magnetohydrodynamic induction equation gives $\nabla \times (\mathbf{u}\times \mathbf{B})=0$. Since both vectors have only poloidal magnetic field components, write $\mathbf{u}\times \mathbf{B}=F_{\varphi }{\mathbf{e}}_{\varphi }$. Axisymmetry impliesFor a field regular on the axis, $C=0$; equivalently one may impose zero toroidal electromotive force. Under this physical condition, $\mathbf{u}$ is parallel to $\mathbf{B}$. Write $\rho \mathbf{u}=k\mathbf{B}$. Continuity and $\nabla \cdot \mathbf{B}=0$ then give $\mathbf{B}\cdot \nabla k=0$. The poloidal magnetic flux function has $\mathbf{B}\cdot \nabla \psi =0$, so locally on each connected magnetic surface,This regular-axis alignment of steady poloidal ideal flow needs its zero-circulation condition: axisymmetry alone does not imply this alignment on a domain excluding the axis. An explicit counterexample is available in a cylindrical annulus. Take nonzero constants $U,b$, positive ${\rho }_c,c_s$, andHere $\psi =-bz$, $\nabla \cdot \mathbf{B}=\nabla \times \mathbf{B}=0$ and $\mathbf{u}\times \mathbf{B}=Ub{\mathbf{e}}_{\varphi }/R$ has zero curl. The flow satisfies continuity, $p^{\prime }=-\rho {\Phi }^{\prime }$, and ${\nabla }^2\Phi =4\pi G\rho$. Its acceleration and Lorentz force density vanish, so it solves the steady equations while $\mathbf{u}$ and $\mathbf{B}$ are perpendicular. The nonzero toroidal circulation and the axial singularity explain why this counterexample is excluded from a regular polar accretion column. The remaining derivation uses the regular, aligned branch.

For a narrow stream tube, integrating continuity gives constant $A\rho |\mathbf{u}|$. On a nonzero-flow tube the magnetohydrodynamic mass loading $k$ is constant, so $A\rho |\mathbf{u}|=A|k||\mathbf{B}|$. Therefore $A|\mathbf{B}|$ is constant and $A\propto |\mathbf{B}{|}^{-1}$. This also follows directly from conserved magnetic flux through a flux tube.

For constant isothermal sound speed, $\nabla p/\rho =c_s^2\nabla \ln (\rho /{\rho }_0)$, with arbitrary reference density ${\rho }_0$. Project momentum along $\mathbf{B}$. The Lorentz force density has no component in that direction, and alignment implies $\mathbf{B}\cdot (\mathbf{u}\cdot \nabla )\mathbf{u}=\mathbf{B}\cdot \nabla (|\mathbf{u}{|}^2/2)$. Thus the isothermal magnetic Bernoulli integral isChanging ${\rho }_0$ only shifts the Bernoulli function by a constant. Since $k$ is constant along the magnetic field,Differentiating the isothermal magnetic Bernoulli integral therefore gives

In the slender polar accretion column, the dipolar flux-tube area is $A(z)=A_{\ast }(z/R_{\ast }{)}^3$. The decreasing axial field is a leading approximation, not an exactly solenoidal field $(0,0,B_z(z))$ throughout a cylinder. Indeed $\nabla \cdot \mathbf{B}={\partial }_z{B}_z\ne 0$ for that literal field. Near the axis a small radial component $B_R\simeq -(R/2){\partial }_z{B}_z=3R{B}_z/(2z)$ supplies the required radial divergence. Its magnitude is smaller by $R/z$, even though its divergence is leading order. Keeping this expanding flux tube while neglecting transverse forces yields the intended one-dimensional model.

Use the positive inward speed $v=-u_z$; the signed sonic velocity is $u_z=-c_s$. Conservation of mass and the Bernoulli equation giveEliminating ${\rho }^{\prime }$ yields the isothermal dipolar accretion equation,A smooth transonic branch requires both sides to vanish at its sonic point, soThis is the critical point of the flow equation; a generic subsonic solution need not cross it. A crossing in the exterior column requires $z_s\ge R_{\ast }$, and a sonic point strictly outside the star requires $z_s>R_{\ast }$. Differentiating the equation at the crossing gives $2(v_s^{\prime }{)}^2=3c_s^2/z_s^2$. For inward accretion that accelerates toward the star, $v_s^{\prime }=-\sqrt{3/2}{c}_s/z_s$.

The reservoir boundary condition fixes $ϵ=c_s^2\ln ({\rho }_{\infty }/{\rho }_0)$. At the sonic point, ${\Phi }_s=-3c_s^2$ and $v_s^2/2=c_s^2/2$, givingFinally use the sonic mass flux through the dipolar flux-tube area:This rate belongs to the smooth transonic branch of the idealized column, with $A_{\ast }$ the supplied total loaded area. The reservoir condition alone does not force every steady solution onto this branch. A real narrow dipolar column must match an outer flow where the slender approximation ceases to apply.

Use the fully developed flow branch of Hartmann flow. Translation invariance along the walls makes the velocity and induced magnetic field functions of $z$ alone. Incompressible flow and zero normal velocity at the walls give $u_z^{\prime }=0$ and hence $u_z=0$. The zero divergence of the magnetic field makes $B_z$ constant, equal to the imposed $B_0$. There is no forcing in $y$; the homogeneous $y$-components obey the same coupled viscous-resistive equations as the $x$-components, with zero boundary data. Multiplying them by $u_y$ and $b_y/({\mu }_0\rho )$, integrating by parts and adding givesThus $u_y=b_y=0$, establishing the asserted forms on this fully developed flow branch. This is a symmetry reduction of the steady channel model, rather than a claim that every possible flow in a channel is translation invariant.

The current density and Lorentz force density areThe advective acceleration vanishes because $\mathbf{u}\cdot \nabla =u(z){\partial }_x$ acts on fields independent of $x$. The $x$-component of the magnetohydrodynamic momentum equation therefore givesThe $z$-component must also balance: it requires ${\partial }_z(p+b^2/(2{\mu }_0))=0$. Hence a compatible pressure is $p(x,z)=p_{\ast }-\rho Gx-b(z{)}^2/(2{\mu }_0)$; the specified streamwise pressure gradient does not require the ordinary pressure to be uniform in $z$. This accounts for the magnetic pressure of the induced field.

The resistive induction equation with constant magnetic diffusivity gives ${\partial }_t\mathbf{B}=\nabla \times (\mathbf{u}\times \mathbf{B})+\eta {\nabla }^2\mathbf{B}$. Since $\mathbf{u}\times \mathbf{B}=(0,-B_{0}u,0)$, its curl has $x$-component $B_0{u}^{\prime }$. Steadiness consequently givesFinally the no-slip boundary condition gives $u(\pm L)=0$, and the normal magnetic field boundary condition gives $b(\pm L)=0$. The coupled equations and all wall conditions follow directly from momentum balance and magnetic induction.

Let $D/Dt={\partial }_t+\mathbf{u}\cdot \nabla$ be the material derivative, and write $h$ for specific enthalpy. The ideal magnetohydrodynamic induction equation and Gauss's law for magnetism giveDotting the momentum equation with $\mathbf{B}$ eliminates the Lorentz force density. Therefore the cross-helicity density satisfiesThe first law of thermodynamics gives $dh=Tds+dp/\rho$ for specific entropy $s$, without requiring uniform entropy. ConsequentlySince $\nabla \cdot \mathbf{B}=0$, the first term on the right is a divergence. This proves the cross-helicity conservation law, with fluxFor the perfect gas used here, $h=\gamma p/[(\gamma -1)\rho ]$.

The remaining magnetic field has $\mathbf{B}=B_{\varphi }(R){\mathbf{e}}_{\varphi }+B_z(R){\mathbf{e}}_z$. The radial Lorentz force density separates into a magnetic pressure gradient and inward magnetic tension:Here $(\mathbf{B}\cdot \nabla )\mathbf{B}$ has radial component $-B_{\varphi }^2/R$ because ${\partial }_{\varphi }{\mathbf{e}}_{\varphi }=-{\mathbf{e}}_R$. Thus magnetostatic equilibrium givesIntegrating the axial component of the magnetostatic Ampère-Maxwell equation gives $I(R)=cR{B}_{\varphi }/2$, since regularity removes the integration constant. Substituting $B_{\varphi }=2I/(cR)$ combines the toroidal terms intoThis is the cylindrical magnetostatic pressure balance equation in terms of enclosed electric current.

In the narrow polar flux tube, the magnetic field lines are almost radial and the transverse width is much smaller than the radial scale. The Lorentz force density has no component along the magnetic field, so the longitudinal magnetically channelled accretion is hydrodynamic. Let $v(r)>0$ denote inward speed. Mass conservation, the Euler equations for an inviscid fluid, and the polytropic equation of state giveUsing the adiabatic sound speed $c^2=\gamma K{\rho }^{\gamma -1}$ and $A\propto r^3$,The Bernoulli equation iswhere the right-hand side comes from matching to a nearly stationary reservoir with negligible gravitational potential. The fixed magnetic field provides transverse confinement; neglecting magnetic forces in this one-dimensional equation concerns their longitudinal component.

Let $\mathbf{\sigma }$ be the stated Maxwell stress tensor. Taking its divergence and usinggivesThe Maxwell equations in vacuum with sources areSubstitution, together with antisymmetry of the cross product, yieldsTherefore the required vector is the electromagnetic momentum densityandThis is local conservation of electromagnetic momentum. The tensor component ${\sigma }_{ij}$ is the flux of $i$-momentum in the $j$ direction in the convention of the question, $\mathbf{g}$ is field momentum per unit volume, and $\rho \mathbf{E}+\mathbf{J}\times \mathbf{B}$ is the Lorentz force density transferred to matter.

For the stated magnetic vector potential, part a givesApplying the same identity once more and using the Ampère-Maxwell equation in the magnetostatic limit givesExpanding the Lorentz force density ${\mathbf{F}}_m=\mathbf{J}\times \mathbf{B}$, using $\nabla \varphi ={\mathbf{e}}_{\varphi }/r$ and the vector triple-product identity, separates its poloidal and azimuthal parts:where

In an axisymmetric magnetostatic state, pressure and gravitational forces are poloidal, so the azimuthal component of the Lorentz force density must vanish. Hence $f=0$ andThe two gradients are locally parallel, so $L\alpha$ is constant on each regular level surface of $\beta$. Thereforefor an arbitrary flux function $F$.

A magnetostatic force-free magnetic field is in equilibrium when its electric current exerts no Lorentz force density:Thus $\mathbf{J}$ is pointwise parallel to $\mathbf{B}$, so for some scalar $\alpha (\mathbf{r})$,Using the magnetostatic Ampère's law, $\nabla \times \mathbf{B}=(4\pi /c)\mathbf{J}$, givesTaking the divergence and using $\nabla \cdot \mathbf{B}=0$ yieldsThe force-free parameter is therefore constant along each magnetic field line.

In a magnetostatic equilibrium, inertia is absent and all forces balance. A force-free magnetic field is the limiting case in which the Lorentz force density vanishes:Thus the current density is parallel to the magnetic field, so for some scalar force-free parameter $\alpha (\mathbf{r})$,Using Ampère's law in magnetostatics givesTaking the divergence and using both the divergence of a curl is zero and $\nabla \cdot \mathbf{B}=0$ yieldssoThe force-free parameter is therefore constant along every magnetic field line.

A planetary dynamo maintains an intrinsic magnetic field by motion of an electrically conducting fluid. The resistive induction equation balances stretching and transport against magnetic diffusion. A sufficiently large magnetic Reynolds number is necessary but does not ensure that the flow regenerates every required field component. Thermal or compositional convection supplies mechanical power; rotation organizes the flow through the Coriolis force, and the Lorentz force density eventually modifies the flow and limits growth. A fluid core can fail to operate a dynamo if it is stably stratified or inadequately driven.

Resistive magnetohydrodynamics retains finite magnetic diffusivity in the coupled evolution of a conducting fluid and its magnetic field. The resistive induction equation allows field diffusion relative to the fluid, while the magnetohydrodynamic momentum equation includes the Lorentz force density. It contrasts with the zero-diffusivity limit of ideal magnetohydrodynamics.

In inviscid magnetostrophic balance within an impenetrable axisymmetric container, with no azimuthal buoyancy or external mechanical torque, the cylindrically integrated Lorentz force density must vanish on every cylinder coaxial with rotation. The azimuthal pressure derivative integrates to zero, and the cylindrical Coriolis integral is proportional to net radial flux, zero by incompressibility and impermeability. Multiplying by cylinder area and lever arm gives the equivalent zero-torque form. This rotating-fluid condition is distinct from magnetic-helicity relaxation to a force-free plasma state. Its original primary source is Taylor's 1963 rotating-fluid paper.

For constant mass density $\rho$, kinematic viscosity $\nu >0$, ideal induction, zero boundary velocity and tangent magnetic field, put $E={\int }_D[\rho |\mathbf{v}{|}^2/2+|\mathbf{B}{|}^2/(2{\mu }_0)]dV$. Multiplying the momentum equation by $\mathbf{v}$ and integrating cancels pressure work and advective flux. The magnetic-energy derivative is minus the integral of Lorentz force density dotted with velocity, canceling the same work in the kinetic-energy equation. Integration by parts then gives the displayed identity. For incompressible flow with zero boundary velocity, $\int |\nabla \mathbf{v}{|}^2=\int |\nabla \times \mathbf{v}{|}^2$.

In a Z-pinch, an axial electric current produces a toroidal magnetic field whose inward magnetic tension confines a plasma. A cylindrical magnetostatic equilibrium balances the resulting Lorentz force density against a pressure gradient.