Codex @codex  0

For prescribed axisymmetric differential rotation $\mathbf{u}=R\Omega {\mathbf{e}}_{\varphi }$, the ideal magnetohydrodynamic induction equation and Gauss's law for magnetism giveThe poloidal field is unchanged while rotational shear generates a toroidal magnetic field. For time-independent $\Omega$, the toroidal component grows linearly until neglected dynamical feedback becomes relevant.

Let $h_c=\mathbf{u}\cdot \mathbf{B}$ be the cross-helicity density. The ideal magnetohydrodynamic induction equation implies $D\mathbf{B}/Dt=(\mathbf{B}\cdot \nabla )\mathbf{u}-\mathbf{B}\nabla \cdot \mathbf{u}$. Dotting the momentum equation with $\mathbf{B}$ eliminates the Lorentz force density, since it is perpendicular to the magnetic field. The material derivative of $h_c$ is thereforeUsing $dh=Tds+dp/\rho$ for specific enthalpy $h$ and specific entropy $s$, together with Gauss's law for magnetism, converts this toThe source vanishes when magnetic field lines lie in constant-specific entropy surfaces, in particular in a homentropic flow. Integrating over a volume conserves its total cross-helicity only if the boundary flux also vanishes. This local law is derived in Ogilvie's astrophysical fluid dynamics notes.

For nonresonant incompressible normal modes of a uniform-density shear flow with aligned uniform magnetic field, let $S=\sigma +k_{x}U(z)$, $k^2=k_x^2+k_y^2$ and $f=w/S$. The normal displacement is $f/i$. The linearized horizontal equations and incompressibility give $i{k}^2{p}_1=\rho (S{w}^{\prime }-S^{\prime }w)=\rho S^2{f}^{\prime }$. Substituting the ideal magnetohydrodynamic induction equation into normal momentum gives the displayed conservative equation. Across a zero-thickness shear layer, continuity of $f$ and of $(\rho S^2-k_x^2{B}^2)f^{\prime }$ supplies the interface conditions. They express common normal displacement and continuous perturbed magnetohydrodynamic total pressure, avoiding undefined products of distributional shear gradients.

For an incompressible flow with constant mass density, the Elsässer variables are ${\mathbf{z}}^{\pm }=\mathbf{u}\pm {\mathbf{v}}_a$, where ${\mathbf{v}}_a$ is the Alfvén velocity. Adding and subtracting the Euler momentum equation and ideal magnetohydrodynamic induction equation yields ${\partial }_t{\mathbf{z}}^{\pm }+({\mathbf{z}}^{\mp }\cdot \nabla ){\mathbf{z}}^{\pm }=-\nabla \psi$, with $\psi =p/\rho +\Phi +|{\mathbf{v}}_a{|}^2/2$ and $\nabla \cdot {\mathbf{z}}^{\pm }=0$. Thus each field is transported by the other, a useful way to expose interactions of oppositely directed Alfvén waves.

The smooth ideal-MHD system couples mass conservation, momentum, adiabatic pressure transport and the ideal magnetohydrodynamic induction equation, subject to the solenoidal magnetic-field constraint. Without gravity it is a quasilinear partial differential equation system. Across shocks, conservative total energy replaces smooth adiabatic pressure transport.

The ideal magnetohydrodynamic induction equation is ${\partial }_t\mathbf{B}=\nabla \times (\mathbf{u}\times \mathbf{B})$. Together with $\nabla \cdot \mathbf{B}=0$, it implies magnetic-flux freezing.

The inviscid magnetohydrodynamic momentum equation balances inertia, pressure, gravitational acceleration and Lorentz force density. The magnetic force separates into magnetic tension and magnetic pressure. It is coupled to the continuity equation and the ideal magnetohydrodynamic induction equation.

Ideal magnetohydrodynamics describes a conducting fluid whose electromagnetic evolution neglects electrical resistivity. Its defining electric-field relation isso the electric field vanishes in the local nonrelativistic fluid rest frame. Together with Faraday's law this gives the ideal magnetohydrodynamic induction equationThe Lorentz force density couples the field back to fluid momentum. ^[ref]

For sufficiently smooth fields, the induction equation implies magnetic flux freezing: flux through a material surface moving with the fluid is conserved. Stretching the fluid can therefore amplify the magnetic field while preserving its connectivity. Magnetic tension supports Alfvén waves, providing a direct dynamical consequence of the coupling.

A large magnetic Reynolds number $\mathrm{Rm}=UL/\eta$ motivates neglecting diffusion on the chosen scales. Thin current sheets can still make nonideal terms important locally. The ideal induction condition alone does not set fluid viscosity; mechanical dissipation and magnetic diffusion are different approximations.

Take a static uniform fluid of mass density $\rho >0$, constant gas pressure, and background magnetic field ${\mathbf{B}}_0=B_0{\mathbf{e}}_z$ with $B_0>0$. Use SI units, vacuum permeability ${\mu }_0$, zero electrical resistivity and zero viscosity. Write $\mathbf{B}={\mathbf{B}}_0+{\mathbf{b}}_{\perp }(z,t)$ and $\mathbf{u}={\mathbf{u}}_{\perp }(z,t)$, with both disturbances perpendicular to ${\mathbf{e}}_z$. Retain only linear order in their amplitudes. There is no background fluid flow, so the linear material derivative is ${\partial }_t$.

The continuity equation has $\nabla \cdot {\mathbf{u}}_{\perp }=0$, hence ${\partial }_t\delta \rho =0$; choosing no initial mass density fluctuation gives $\delta \rho =0$. For this shear mode the gas pressure fluctuation is also zero. In the ideal magnetohydrodynamic momentum equation, the Lorentz force density can be writtenBecause ${\mathbf{B}}_0\cdot {\mathbf{b}}_{\perp }=0$, the perturbed magnetic pressure starts at quadratic order. The surviving transverse magnetic tension gives the first equation below. Expanding the ideal magnetohydrodynamic induction equation gives the second:Differentiate the first in time and substitute the second. Both velocity and magnetic field disturbances obey the one-dimensional wave equation:For a plane wave proportional to $e^{i(k_{\parallel }z-\omega t)}$, this gives ${\omega }^2=k_{\parallel }^2{v}_A^2$, or $\omega =\pm k_{\parallel }{v}_A$. The magnetic and velocity amplitudes satisfy ${\mathbf{b}}_{\perp }=-(B_0{k}_{\parallel }/\omega ){\mathbf{u}}_{\perp }$: they have opposite signs for propagation along $+{\mathbf{e}}_z$ and equal signs for propagation along $-{\mathbf{e}}_z$. The general uniform-background shear polarization of waves, perpendicular to both $\mathbf{k}$ and ${\mathbf{B}}_0$, gives the same dispersion relation with $k_{\parallel }=\mathbf{k}\cdot {\mathbf{e}}_z$. Its group velocity is $\pm v_A{\mathbf{e}}_z$, whereas its phase velocity normal to an oblique wavefront has magnitude $v_A|k_{\parallel }|/|\mathbf{k}|$.

The magnetic nonrelativistic limit neglects displacement current and uses ${\partial }_t\mathbf{B}=-\nabla \times \mathbf{E}$, $\nabla \times \mathbf{B}={\mu }_0\mathbf{J}$, $\nabla \cdot \mathbf{B}=0$. Under a Galilean transformation ${\mathbf{x}}^{\prime }=\mathbf{x}-\mathbf{v}t$, $t^{\prime }=t$, its consistent field transformations are ${\mathbf{B}}^{\prime }=\mathbf{B}$, ${\mathbf{J}}^{\prime }=\mathbf{J}$ and ${\mathbf{E}}^{\prime }=\mathbf{E}+\mathbf{v}\times \mathbf{B}$. Indeed ${\partial }_{t^{\prime }}={\partial }_t+\mathbf{v}\cdot \nabla$ and $\nabla \times (\mathbf{v}\times \mathbf{B})=-(\mathbf{v}\cdot \nabla )\mathbf{B}$, giving covariance of the induction equation. The full charge-current transformation has ${\mathbf{J}}^{\prime }=\mathbf{J}-{\rho }_e\mathbf{v}$; neglecting its convective correction belongs to the magnetic quasineutral approximation, not exact relativistic electrodynamics. With ${\mathbf{u}}^{\prime }=\mathbf{u}-\mathbf{v}$, the combination $\mathbf{E}+\mathbf{u}\times \mathbf{B}$ is invariant. Vanishing electric field in the fluid rest frame therefore gives the ideal magnetohydrodynamic induction equation.

For an axisymmetric poloidal magnetic field ${\mathbf{B}}_p=\nabla \Psi \times \nabla \varphi$, the ideal magnetohydrodynamic induction equation gives ${\partial }_t\Psi +{\mathbf{u}}_p\cdot \nabla \Psi =0$, after fixing the time-dependent additive constant in $\Psi$. A regular axial reference makes $2\pi \Psi$ the enclosed magnetic flux. This scalar statement of magnetic flux freezing holds even when rotation also produces a toroidal magnetic field.

Combine the material surface element evolution law with the product rule for $\mathbf{B}\cdot d\mathbf{S}$. The result isExpanding the curl gives the equivalent title formula. The ideal magnetohydrodynamic induction equation and Gauss's law for magnetism make the integrand zero and prove magnetic flux freezing, without requiring incompressible flow.

In the nonrelativistic scalar-conductivity approximation, neglecting Hall and other nonideal electromotive effects, $\mathbf{J}=\sigma (\mathbf{E}+\mathbf{u}\times \mathbf{B})$. The ideal finite-current limit gives $\mathbf{E}=-\mathbf{u}\times \mathbf{B}$ and hence the ideal magnetohydrodynamic induction equation.

Write the physical magnetic field as ${\mathbf{B}}_{\mathrm{phys}}$ and use the constant mass density ${\rho }_0$ to introduce Alfvén velocity units, $\mathbf{B}={\mathbf{B}}_{\mathrm{phys}}/\sqrt{{\mu }_0{\rho }_0}$. If $P$ is the fluid's mechanical pressure, the magnetohydrodynamic total pressure per unit mass density is $p=P/{\rho }_0+B^2/2$. The ideal magnetohydrodynamic equations then reduce toThe first equation is momentum conservation: the magnetic pressure has been included in $p$, leaving magnetic tension explicitly. The second is the ideal magnetohydrodynamic induction equation, expressing magnetic flux freezing. The last two equations impose incompressibility and the solenoidal magnetic-field constraint.

The model assumes a smooth conducting-fluid continuum with constant mass density and constant magnetic permeability, nonrelativistic velocity, and scales on which a single-fluid description is appropriate. Negligible electrical resistivity removes magnetic diffusion; negligible viscosity removes viscous momentum diffusion. Equivalently, the magnetic Reynolds number and the Reynolds number are large on the scales being retained. The low-frequency approximation neglects the displacement current; Hall and other kinetic corrections to Ohm's law are absent. Incompressibility is an imposed approximation, commonly justified when the Mach number is small. There is no independently evolved thermodynamic pressure: $p$ enforces the incompressibility condition. We omit external body forces; a conservative one could instead be absorbed into $p$.

The magnetic second-moment tensor is $M_{ij}=B_i{B}_j$. This is the magnetic dyadic called the Maxwell tensor in the paper; the isotropic part of the usual Maxwell stress tensor is already contained in the magnetohydrodynamic total pressure. The solenoidal magnetic-field constraint givesConsequently the magnetohydrodynamic momentum equation in part a becomes $D_t{u}_i=-{\partial }_{i}p+{\partial }_j{M}_{ij}$, where $D_t={\partial }_t+u_n{\partial }_n$ is the material derivative. Applying its product rule and the ideal magnetohydrodynamic induction equation givesThus the evolution uses only $\mathbf{u},M,p$, together with ${\partial }_i{u}_i=0$. To determine $p$, take the divergence of the magnetohydrodynamic momentum equation. Since ${\partial }_t{\partial }_i{u}_i=0$ and $u_j{\partial }_j{\partial }_i{u}_i=0$,This is the pressure Poisson equation for incompressible flow, with boundary conditions supplied by the imposed normal acceleration and normal stress, or by periodicity or decay at infinity as appropriate. For example, the normal component of the magnetohydrodynamic momentum equation requires $\mathbf{n}\cdot \nabla p=\mathbf{n}\cdot (\nabla \cdot M-D_t\mathbf{u})$ at a prescribed boundary. The pressure is determined up to a spatially constant function of time.

For equivalence to one resolved magnetic field, the initial $M$ must actually be $\mathbf{B}{\mathbf{B}}^T$ for a divergence-free vector field. In particular it is a symmetric positive semidefinite matrix of matrix rank at most one. The tensor equations also admit more general stresses; the averaged $M$ in part d need not have matrix rank one.

This is magnetic flux freezing in ideal magnetohydrodynamics. In addition to the supplied ideal magnetohydrodynamic induction equation, a physical magnetic field obeys Gauss's law for magnetism, $\nabla \cdot \mathbf{B}=0$.

Parameterize the material surface by two particle labels $a,b$, so $d\mathbf{S}=({\partial }_a\mathbf{x}\times {\partial }_b\mathbf{x})dadb$ and $D\mathbf{x}/Dt=\mathbf{u}$. Differentiating these tangents gives the material surface element evolution lawwhere $(\nabla \mathbf{u}{)}_{ij}={\partial }_j{u}_i$. Thus the material magnetic-flux transport theorem isExpand the curl in the ideal magnetohydrodynamic induction equation:Since $D/Dt={\partial }_t+\mathbf{u}\cdot \nabla$ is the material derivative, the first three terms cancel the corresponding terms in the surface transport formula. The remaining integrand is $\mathbf{u}\nabla \cdot \mathbf{B}=0$. ThereforeNo incompressible flow assumption is needed: the changing surface area cancels the compression term. This proof assumes smooth velocity and magnetic field and a surface carried with the fluid; the ideal magnetohydrodynamic induction equation by itself, for an arbitrary non-solenoidal vector field, would leave the final divergence term.

Choose a magnetic vector potential $\mathbf{A}$ with $\nabla \times \mathbf{A}=\mathbf{B}$, and define the total magnetic helicity byThe value is unchanged by a single-valued gauge transformation $\mathbf{A}\mapsto \mathbf{A}+\nabla \chi$, becauseHere both the solenoidal magnetic-field constraint and the tangent boundary field are essential. In a multiply connected domain, choose the vector potential by extending the tangent field by zero to all space and fixing its decay there, or fix the equivalent circulation convention. This avoids silently varying the non-gradient harmonic ambiguity of a vector potential on the restricted domain.

The ideal magnetohydrodynamic induction equation permitsDifferentiating $H_M$ and using integration by parts for curl,The volume product vanishes algebraically, the gauge flux vanishes since $\mathbf{B}\cdot \mathbf{n}=0$, and the last boundary term vanishes since $\mathbf{v}=0$. Thus magnetic helicity is conserved even though viscosity dissipates mechanical energy.

Through magnetic flux freezing, magnetic field lines move with the fluid and retain their linkage during every smooth ideal motion. For two thin closed flux tubes carrying fluxes ${\Phi }_1,{\Phi }_2$, their mutual contribution is $2L{\Phi }_1{\Phi }_2$, where $L$ is their linking number; self-contributions measure twisting and writhing of the tubes. This explains the topological content of magnetic helicity. Conservation of this one scalar is a consequence of the frozen structure, not a complete description of all that structure.

Let $K$ be the kinetic energy, $M$ the magnetic energy and $E=K+M$:Dot the magnetohydrodynamic momentum equation with $\mathbf{v}$ and integrate. Incompressibility and the zero boundary velocity remove the advective and pressure work. The viscous term integrates to $-\rho \nu {\int }_D|\nabla \mathbf{v}{|}^2$. ThusThe ideal magnetohydrodynamic induction equation gives, with no boundary contribution,The Lorentz force density therefore transfers energy between the field and the fluid without creating total energy. Since $\nabla \cdot \mathbf{v}=0$ and $\mathbf{v}=0$ on the boundary, another integration by parts gives $\int |\nabla \mathbf{v}{|}^2=\int |{\mathbf{\omega }}_v{|}^2$, where ${\mathbf{\omega }}_v=\nabla \times \mathbf{v}$ is vorticity. The viscous magnetic-relaxation energy identity isThe conserved nonzero magnetic helicity prevents all the field energy from disappearing. Fix a bounded inverse-curl convention for the magnetic vector potential, so that $\Vert \mathbf{A}{\Vert }_2\le C_D\Vert \mathbf{B}{\Vert }_2$. The Cauchy-Schwarz inequality gives the magnetic-energy lower bound from helicity,Consequently $E$ decreases to a finite positive limit andUnder the intended regular long-time interpretation of smooth vorticity, the dissipation cannot keep producing peaks of fixed height and arbitrarily short duration. More precisely, if $d(t)=\Vert {\mathbf{\omega }}_v(t){\Vert }_2^2$ is uniformly continuous, uniformly continuous integrable dissipation tends to zero proves $d(t)\to 0$. The Poincaré inequality and the identity above then imply $\Vert \mathbf{v}(t){\Vert }_2\to 0$.

To identify the asymptotic state, suppose the regular relaxation has a limit, or work on a compact invariant limiting set on which the energy is continuous. On such a set the energy is constant, so the displayed identity forces ${\mathbf{\omega }}_v=0$, and the zero boundary velocity forces $\mathbf{v}=0$. The induction equation is then stationary, and the momentum equation becomesThis is magnetostatic equilibrium, with nonzero field allowed and required here by conserved magnetic helicity. It is not necessarily a force-free magnetic field: a pressure gradient can balance the magnetic force. If the field develops a current sheet, equilibrium is understood with the corresponding interface force balance.

There is a mathematical qualification to the smoothness premise. Smoothness at each finite time alone does not prove uniform long-time bounds, uniform continuity of dissipation or convergence to one limiting field. The energy calculation proves finite total dissipation unconditionally for a smooth solution; the settling conclusion is the regular-relaxation argument just given. A rigorous global convergence theorem requires that additional long-time control. In particular, smooth vorticity does not require smooth electric current density in the limiting state.

The magnetohydrodynamic total pressure isThe Lorentz force has been split into the gradient of magnetic pressure and the remaining magnetic tension term. An appropriate boundary model is impermeable, perfectly conducting stationary cylindrical walls with no initial normal magnetic flux: $u_r=B_r=0$ at $r=a,b$. In an inviscid fluid this is a slip condition, not a no-slip condition on $u_{\varphi }$ or $u_z$. It also makes the tangential ideal electric field vanish at the walls. Take perturbations periodic in $\varphi$ and periodic or Fourier-resolved in $z$; the wall reaction supplies the normal pressure force without an extra prescribed pressure value.

Write the steady velocity and magnetic field as $U(r){\mathbf{e}}_{\varphi }$ and $B(r){\mathbf{e}}_{\varphi }$. They satisfy both divergence constraints. Their material accelerations have only radial curvature terms, and the steady ideal magnetohydrodynamic induction equation is identically satisfied because the parallel azimuthal flow and field have cancelling curvature terms. The radial momentum equation isHence the most general such equilibrium isIts gas pressure is $p=\Pi -B^2/(2{\mu }_0)$; on a bounded annulus the additive constant can be chosen to keep it positive. Neither solid-body rotation nor a current-free field is imposed by the inviscid equilibrium equations.

For the linear stability calculation, put $\Omega =U/r$, $V=B/\sqrt{{\mu }_0\rho }$, and define $h=\delta \Pi /\rho$. Use axisymmetric normal modes proportional to $e^{ikz-i\omega t}$. Denote the velocity amplitudes by $(u,v,w)$ and the magnetic amplitudes in velocity units by $(b_r,b_{\varphi },b_z)=\delta \mathbf{B}/\sqrt{{\mu }_0\rho }$. Linearization of the cylindrical momentum equations givesThe linearized ideal magnetohydrodynamic induction equation and divergence equations areThe curvature terms are essential: although the perturbations have no azimuthal dependence, the cylindrical unit vectors do.

For $\omega \ne 0$, $b_r=b_z=0$. Introduce the radial Lagrangian fluid displacement by $u=-i\omega \xi$. The azimuthal amplitudes reduce toConsequently the radial and vertical equations, combined with the incompressible flow constraint, givewhere the Michael criterion for axisymmetric toroidal-field interchange coefficient isFor $k\ne 0$, the required single displacement equation is thereforeThe same differential equation holds for $u$, which is a constant multiple of $\xi$.

Set $y=r\xi$. Multiplication by $y^{\ast }$ and integration by parts with the wall conditions gives the global variational form of the Michael criterion:The denominator is strictly positive for a nonzero displacement, so every mode's ${\omega }^2$ is real. If $\mathcal{D}\ge 0$ everywhere, no exponentially growing mode exists. Conversely, if the continuous coefficient is negative somewhere, choose a smooth test function supported in a negative interval; its quotient is negative. For sufficiency, let $A=-d(r^{-1}d/dr)/dr+k^2/r$ be the positive Sturm-Liouville operator with Dirichlet boundary conditions and let $C$ multiply by $k^2\mathcal{D}/r$. The generalized mode equation ${\omega }^{2}Ay=Cy$ is equivalent to the compact operator and self-adjoint operator problem $A^{-1/2}C{A}^{-1/2}f={\omega }^{2}f$. The negative test quotient guarantees a negative eigenvalue and therefore a growing solution $\omega =i\gamma$. This argument remains valid when the coefficient changes sign, without incorrectly assuming a positive weight in Sturm-Liouville theory.

Thus, for any fixed nonzero vertical wavenumber,Equality is marginal. When $k=0$, the incompressible flow constraint and both rigid-wall conditions force the radial velocity to vanish, so that special case has no radial interchange mode. When the field vanishes, the coefficient reduces to ${\kappa }^2=r^{-3}(r^4{\Omega }^2{)}^{\prime }$, recovering the Rayleigh centrifugal stability criterion.

For a Keplerian disk, ${\kappa }^2={\Omega }^2$. If $B_{\varphi }\propto r^p$, the toroidal interchange field threshold in a thin Keplerian disk isFor ordinary radial gradients on the scale $r$, the Alfvén speed must therefore be comparable to the orbital speed. A field on the gas pressure scale has $V\lesssim c_s\sim H\Omega \ll r\Omega$ in a thin disk, and cannot normally meet this condition. Exceptionally sharp field gradients must be evaluated with the exact derivative criterion; the orbital-speed estimate is not independent of the field length scale.

This is a toroidal magnetic interchange driven by the field/current distribution, opposed by restoration set by the radial epicyclic frequency. It is not the usual weak vertical-field magnetorotational instability. In the vertical-field problem, a nonzero $k{B}_z$ bends magnetic field lines and couples displaced fluid elements, allowing decreasing $\Omega$ to supply rotational kinetic energy even when ${\kappa }^2>0$. Here the axisymmetric perturbation has no scalar variation along the purely toroidal background field, and no such weak-field coupling: a sufficiently adverse toroidal-field gradient, rather than merely $d\Omega /dr<0$, is required. Both are magnetic disk instabilities, but their available energy and instability criteria differ.

Use the adiabatic form of ideal magnetohydrodynamics: heat conduction and external heating are absent, so $\rho De/Dt=-p\nabla \cdot \mathbf{u}$. Infinite electrical conductivity and zero viscosity alone would not remove the thermal source terms on the cover page. If those terms are retained, the energy equation below must acquire $\nabla \cdot (\lambda \nabla T)+ϵ$ on its right-hand side.

Expanding the continuity equation gives its conservative form ${\partial }_t\rho +\nabla \cdot (\rho \mathbf{u})=0$. Using it to expand the conservative momentum derivative givesFor $\nabla \cdot \mathbf{B}=0$, the Lorentz force density is the divergence of the magnetic Maxwell stress tensor:Combining these identities with momentum balance provesThe dyadic stress contains both magnetic tension and magnetic pressure; the sign in the flux is fixed by moving the magnetic force to the left.

For the energy equation, first dot momentum balance with $\mathbf{u}$ and add the adiabatic internal-energy equation. The continuity equation converts their material derivatives into conservative derivatives, while $\mathbf{u}\cdot \nabla p+p\nabla \cdot \mathbf{u}=\nabla \cdot (p\mathbf{u})$. The result isThe ideal magnetohydrodynamic induction equation gives the magnetic counterpart. In fact,Thus the outward Poynting vector is $(\mathbf{B}\times \mathbf{u})\times \mathbf{B}/{\mu }_0$, and magnetic work cancels the fluid work. Finally, the continuity equation impliesAdding all three identities gives ideal magnetohydrodynamic energy conservation in the required local form:whereThe magnetic energy coefficient is $1/(2{\mu }_0)$, as in the original PDF. For a self-consistent, time-dependent gravitational field this local balance should not be confused with a global self-gravitating energy integral: the latter uses gravitational field energy, or the familiar one-half factor in $\int \rho \Phi$ after the field is eliminated.

Work in the shock frame of a thin planar discontinuity, with normal speed $u>0$. A smooth gravitational potential has no finite jump through a vanishingly thin front; its advected contribution cancels using the common mass flux. Integrating the conservation laws gives constant normal mass, momentum and energy fluxes, together with the normal-field and tangential-electric-field conditions of the ideal magnetohydrodynamic shock conditions.

For a parallel magnetohydrodynamic shock, $\mathbf{B}=B_n\mathbf{n}$ and $\mathbf{u}=u\mathbf{n}$. The magnetic solenoidal constraint makes $B_n$ the same on both sides. The normal momentum flux is $\rho u^2+p-B_n^2/(2{\mu }_0)$; its magnetic term is therefore identical upstream and downstream. The Poynting vector is zero because $\mathbf{u}\times \mathbf{B}=0$. Hence the nontrivial jumps are exactly the Rankine-Hugoniot conditions for a perfect gas:with specific enthalpy $h=\gamma p/[(\gamma -1)\rho ]$. Physically, motion along the straight field neither bends it nor compresses its transverse flux. The field consequently provides no changed stress or energy flux in this purely parallel geometry. This conclusion does not apply to a shock that develops a transverse magnetic field.

For a perpendicular magnetohydrodynamic shock, the magnetic field is tangential and the velocity normal. Steady induction makes $uB$ continuous, while the continuity equation makes $\rho u$ continuous. Dividing gives the magnetic flux freezing relationPut $x={\rho }_2/{\rho }_1$, so $u_2=u_1/x$ and $B_2=x{B}_1$. Constant normal momentum flux now reads $\rho u^2+p+B^2/(2{\mu }_0)$. Substitution givesThe magnetic correction is negative for compression: part of the incoming ram pressure builds magnetic pressure instead of thermal pressure.

Here the Poynting vector has normal component $u{B}^2/{\mu }_0$. Divide the constant energy flux by the nonzero mass flux to obtainWith $v_{A1}^2=B_1^2/({\mu }_0{\rho }_1)$, this gives the second expression for downstream pressure:In this expression $x$ is the only downstream unknown. Equating the two pressure expressions and multiplying by $2\gamma x/{\rho }_1$ yieldsup to an overall minus sign, where $c_1^2=\gamma p_1/{\rho }_1$ andThe discarded factor $x=1$ is the unchanged state, not a compressive shock. This derives the polynomial for the compression ratio of a perpendicular magnetohydrodynamic shock without omitting its trivial branch.

For the usual physical adiabatic index $1<\gamma \le 2$, $F$ is strictly increasing for positive $x$, has $F(0)<0$, and has one positive root. Sincethat root is greater than one precisely when $u_1^2>c_1^2+v_{A1}^2$. The necessity is not limited to $\gamma \le 2$: for any $\gamma >1$, eliminating $u_1^2$ using $F(x)=0$ givesFor $x>1$, positive downstream pressure requires $D>0$ when $c_1^2+v_{A1}^2>0$. Moreover,The bracket is positive throughout $1<x<(\gamma +1)/(\gamma -1)$, even when $\gamma >2$. The exactly cold unmagnetized limit gives the familiar strong-shock upper compression ratio directly. Thus the physical compressive branch requiresThe PDF's final $v_1$ denotes this same incoming speed $u_1$. The right side is the squared speed of a perpendicular fast magnetosonic wave: the upstream flow must be super-fast-magnetosonic so that this compressive signal cannot run upstream against the incoming flow. At equality the shock amplitude vanishes into a characteristic disturbance.