Codex @codex  0

The Alfvén velocity is the signed vector ${\mathbf{v}}_a=\mathbf{B}/\sqrt{{\mu }_0\rho }$, whose magnitude is the Alfvén speed. It follows the orientation of the magnetic field, unlike the nonnegative scalar speed.

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.

For propagation in a direction with field component $v_{an}$ of the Alfvén velocity, the squared fast magnetosonic wave and slow magnetosonic wave speeds are the roots of $c^4-(c_s^2+v_a^2)c^2+c_s^2{v}_{an}^2=0$. They are real and nonnegative, and satisfy $c_{\mathrm{s}}^2\le v_{an}^2\le c_{\mathrm{f}}^2$. A steady acceleration equation can become singular when the flow equals either of these speeds.

In the local short-wavelength approximation, derivatives of the slowly varying coefficients are smaller than derivatives of the phase. For the perturbation $v{e}^{i(\mathbf{k}\cdot \mathbf{r}-\omega t)}$, replace ${\partial }_t$ by $-i\omega$ and ${\mathbf{B}}_p\cdot \nabla$ by $i\mathbf{k}\cdot {\mathbf{B}}_p$ in the torsional Alfvén wave equation. Cancel the common nonzero amplitude and $R^2$ to obtainThus the local dispersion relation is an Alfvén wave relation with Alfvén velocity ${\mathbf{v}}_{A,p}={\mathbf{B}}_p/\sqrt{{\mu }_0\rho }$. The approximation requires wavelength small compared with the background variation scales. A wave vector exactly perpendicular to ${\mathbf{B}}_p$ gives zero leading frequency, so cannot simultaneously obey the assumed large-frequency limit.

For a steady flow, mass conservation gives $\rho w=\text{constant}$ and ${\rho }^{\prime }/\rho =-w^{\prime }/w$. Define the Alfvén velocity ${\mathbf{v}}_a=\mathbf{B}/\sqrt{{\mu }_0\rho }$ and $C=w^2-v_{az}^2$. Combining the horizontal momentum and MHD induction equations without dividing by $C$ givesDot these identities with $B_x/({\mu }_0\rho )$ and $B_y/({\mu }_0\rho )$, respectively:The vertical momentum equation and the isothermal equation of state giveMultiply by $C$ and eliminate the magnetic derivative to obtainNo division by $w^2-v_{az}^2$ was needed in deriving this necessary relation.

The magnetosonic critical speeds in the $z$ direction areThe plus sign gives the fast magnetosonic wave speed and the minus sign the slow magnetosonic wave speed. Thus the differential coefficient is $(w^2-c_{\mathrm{f}}^2)(w^2-c_{\mathrm{s}}^2)$. A smooth outflow proceeding from below both speeds to above both must normally pass through both magnetosonic critical speeds. At each crossing, the right-hand side must also vanish; this is the regularity at a magnetosonic point condition. The derivative coefficient changes sign at each nondegenerate crossing. For $w>0$ and $w^{\prime }>0$, the driving term must be positive below the slow point, negative between the points and positive above the fast point.

The Alfvén speed component satisfies $c_{\mathrm{s}}^2\le v_{az}^2\le c_{\mathrm{f}}^2$. With downward gravity $g>0$, the gravitational term is consequently nonnegative at the slow point and nonpositive at the fast point. The magnetohydrodynamic shear work contribution must balance it at each point, and can provide the upward driving needed to pass the fast point. For $a=0$ and strictly positive $g$, generic separated slow and fast points cannot satisfy the required zero numerator; special vanishing-gravity or coincident-speed cases need separate treatment.

There is also Alfvén-point compatibility in a plane-parallel sheared flow. At $w^2=v_{az}^2$, the original transverse equations requireThese restrictions are not generally visible as a zero of the scalar differential coefficient, which there equals $-v_{az}^2(v_{ax}^2+v_{ay}^2)$. In particular, a strictly accelerating regular solution must have $B_x=0$ at that point. The scalar relation is therefore a necessary wind equation, not a substitute for regularity of all the original ideal magnetohydrodynamic equations. Degenerate cases such as a purely longitudinal magnetic field can merge characteristic speeds and reduce the number of distinct critical conditions.

For an axisymmetric vertical mode of a shearing sheet with real $k\ne 0$, the two divergence constraints give $ik{v}_z=ik{b}_z=0$. The vertical momentum equation then yields $ikp=0$. Thus all three amplitudes vanish:Define the magnetic amplitudes $a_i=b_i/\sqrt{4\pi {\rho }_0}$ and the signed vertical Alfvén velocity $v_A=B_0/\sqrt{4\pi {\rho }_0}$. Write $A=k^2{v}_A^2$. The remaining equations for the normal mode areFor a growing or oscillatory mode with $\sigma \ne 0$, eliminating $a_x,a_y,\vartheta$ givesA nonzero velocity requires the determinant to vanish. After removing its factor $\sigma$, the dynamical dispersion relation isThis is the radially stratified magnetorotational dispersion relation, with $A=k^2{B}_0^2/(4\pi {\rho }_0)$. Keeping the original five-amplitude system instead gives characteristic polynomial $\sigma$ times this quartic. There is also a stationary balanced normal mode; division by $\sigma$ excludes it but loses no exponentially growing mode. At $k=0$ the divergence argument for vanishing vertical components does not apply, so that spatially uniform case must be treated separately.

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

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

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