Codex @codex  0

Use the aligned poloidal flow from part (i), and define the magnetohydrodynamic mass loading $\eta =\rho \alpha$, so $\rho {\mathbf{u}}_p=\eta {\mathbf{B}}_p$. Axisymmetric mass conservation and $\nabla \cdot {\mathbf{B}}_p=0$ give ${\mathbf{B}}_p\cdot \nabla \eta =0$. Introduce the poloidal magnetic flux function $\Psi$ by ${\mathbf{B}}_p=R^{-1}\nabla \Psi \times {\mathbf{e}}_{\varphi }$, and write $u_{\varphi }=R\Omega$. ThenThe steady induction equation says the curl of this expression is zero. Hence $\nabla {\Omega }_F\times \nabla \Psi =0$, wherethe field-line angular velocity ${\Omega }_F$ is constant on each connected magnetic surface and along its field lines.

To control the difference between matter rotation and field-line rotation, use the azimuthal component of the ideal magnetohydrodynamic momentum equation. The inertial and magnetic terms are respectively $\rho R^{-1}{\mathbf{u}}_p\cdot \nabla (R{u}_{\varphi })$ and $({\mu }_{0}R{)}^{-1}{\mathbf{B}}_p\cdot \nabla (R{B}_{\varphi })$; pressure and gravity have no azimuthal derivative. Multiplying by $R$ and using constant mass loading givesThus the magnetohydrodynamic angular-momentum invariant $L=R^2\Omega -R{B}_{\varphi }/({\mu }_0\eta )$ is constant along the line when $\eta \ne 0$. Let $v_{Ap}=|{\mathbf{B}}_p|/\sqrt{{\mu }_0\rho }$ and define $M_{Ap}^2=|{\mathbf{u}}_p{|}^2/v_{Ap}^2={\mu }_0{\eta }^2/\rho$. Eliminating $B_{\varphi }$ between the two invariants givesFor a highly sub-Alfvénic poloidal flow, $M_{Ap}^2\ll 1$, this is $\Omega \simeq {\Omega }_F$ along a nonsingular field-line segment. This proves sub-Alfvénic isorotation in a steady axisymmetric wind, with the correction explicitly displayed. If the poloidal flow is exactly zero, the induction equation directly gives ${\mathbf{B}}_p\cdot \nabla \Omega =0$, namely Ferraro's law of isorotation. The approximate claim presumes regular invariants: small $M_{Ap}$ alone is not a uniform relative bound if $L/R^2$ is allowed to become arbitrarily large.