Codex @codex  0

The Alfvén Mach number compares a characteristic flow speed with the Alfvén speed $v_A=B/\sqrt{{\mu }_0\rho }$. In an axisymmetric magnetohydrodynamic wind, the relevant poloidal quantity is $M_{Ap}=|{\mathbf{u}}_p|\sqrt{{\mu }_0\rho }/|{\mathbf{B}}_p|$. With magnetohydrodynamic mass loading $\rho {\mathbf{u}}_p=\kappa {\mathbf{B}}_p$, its square is ${\mu }_0{\kappa }^2/\rho$. The Alfvén surface is the locus $M_{Ap}=1$, and a smooth crossing requires the corresponding Alfvén-surface regularity condition for an axisymmetric wind. The total-speed and poloidal-speed definitions should be distinguished when there is substantial rotation or a strong toroidal magnetic field.

The wind integrals require a steady axisymmetric magnetohydrodynamic wind in ideal magnetohydrodynamics, with no viscous azimuthal torque or imposed toroidal loop voltage. Axisymmetry by itself is insufficient. Assume smooth connected poloidal flux surfaces with nonzero poloidal magnetic field; functions of the flux label below are understood on each such branch. In cylindrical coordinates, the given poloidal magnetic flux function obeysFor a steady ideal wind, the electric field is $\mathbf{E}=-\mathbf{u}\times \mathbf{B}$ and $\nabla \times \mathbf{E}=0$. Its toroidal component satisfies ${\partial }_z{E}_{\varphi }=0$ and ${\partial }_s(s{E}_{\varphi })=0$, so $E_{\varphi }=C/s$. Regularity at the axis, or equivalently the absence of an imposed loop voltage in an annular domain, gives $C=0$. Thus $({\mathbf{u}}_p\times {\mathbf{B}}_p{)}_{\varphi }=0$ and ${\mathbf{u}}_p=f{\mathbf{B}}_p$. The continuity equation and $\nabla \cdot {\mathbf{B}}_p=0$ now yieldConsequently $\kappa$ is the magnetohydrodynamic mass loading. Including the toroidal magnetic field, rewrite the full velocity asSince $\mathbf{u}\times \mathbf{B}=\omega \nabla \chi$, steady ideal induction gives $0=\nabla \omega \times \nabla \chi$. Hence $\omega =\omega (\chi )$, establishingThe invariant $\omega$ is the field-line angular velocity of an axisymmetric wind, not generally the fluid's angular velocity $\Omega$. When an ideal magnetic surface is anchored to a rigidly rotating star, $\omega$ equals the stellar angular velocity. With no poloidal outflow this reduces to Ferraro's law of isorotation.

To calculate the magnetic torque, use the azimuthal momentum equation. Axisymmetric pressure and Newtonian gravity have no azimuthal components; the magnetic pressure has none either. The azimuthal acceleration and magnetic tension contain the curvature terms characteristic of cylindrical coordinates:Multiplication by $s$ gives the Maxwell torque conservation in an axisymmetric wind equationUsing $u_{\varphi }=s\Omega$, $\rho {\mathbf{u}}_p=\kappa {\mathbf{B}}_p$ and ${\mathbf{B}}_p\cdot \nabla \kappa =0$ provesThis is the total angular momentum transport per unit poloidal magnetic flux. The first term is the material contribution and the second is the contribution of the Maxwell stress tensor. The specific magnetohydrodynamic angular-momentum invariant is $L=\ell /\kappa$ when $\kappa \ne 0$; its normalization differs from the paper's $\ell$.

Define the poloidal Alfvén Mach number bySubstituting $\Omega =\omega +\kappa B_{\varphi }/(\rho s)$ into $\ell$ givesAt a smooth crossing of the Alfvén surface, $M_A^2=1$ and the toroidal magnetic field is finite. The Alfvén-surface regularity condition for an axisymmetric wind is thereforeOne can also expose the apparent singularity explicitly:Finite passage through the Alfvén surface requires the numerator to vanish with the denominator. Equality of the numerator zeros is necessary; a globally smooth wind must also satisfy its other dynamical critical conditions. Here $s_A$ is the cylindrical distance from the rotation axis, not automatically the spherical Alfvén radius.

For a tube carrying outward mass flux $d\dot{M}=\kappa {\mathbf{B}}_p\cdot d\mathbf{S}$, the outward angular momentum flux isThus the stellar torque is ${\dot{J}}_{\star }=-\int \omega s_A^{2}d\dot{M}$, with the integral over the escaping wind. Magnetic stresses give each mass element an effective angular-momentum lever arm equal to the cylindrical Alfvén crossing radius. In a strongly sub-Alfvénic region the fluid approximately corotates with the magnetic surface; beyond that region it can carry the accumulated angular momentum away. For an outward field and positive $\omega$, a super-Alfvénic expanding tube with $s>s_A$ has $B_{\varphi }<0$, so the magnetic stress transports positive angular momentum outwards. If $s_A$ greatly exceeds the stellar surface lever arm, a modest mass-loss rate can cause substantial wind-driven magnetic braking of a solar-type star. An unmagnetized outflow would instead carry approximately the material angular momentum fixed at its launch radius.

All invariants in this part refer to the selected magnetic field line. Its squared poloidal Alfvén number isAt the Alfvén surface crossing it equals unity, so ${\rho }_a={\mu }_0{k}^2$ and $M_{Ap}^2=1/y$. Subtract the two azimuthal integrals in part (a) to findFinite azimuthal fields and velocity at $y=1$ require the Alfvén-surface regularity condition for an axisymmetric wind $\ell =\omega R_a^2$. HenceHere $v_{\varphi }$ denotes the azimuthal velocity relative to field-line rotation. For gravity of a point mass at $z=0$, $\Phi =-GM/R=-GM/(R_{a}x)$.

Insert these quantities in the corotating energy integral with $w=0$, and divide by the gravitational energy scale $GM/R_a$. This gives the cold radial magnetohydrodynamic wind integral,withBoth are dimensionless. The apparent $0/0$ at $x=y=1$ must be evaluated on the smooth solution branch, not assigned an arbitrary value.

A constant-energy curve obeys $f_x+f_{y}dy/dx=0$. At a nondegenerate magnetosonic point, the vanishing coefficient of the mass density derivative requires simultaneous vanishing of its numerator:The curve must admit a real crossing slope. Expanding about such a point gives $f_{yy}{m}^2+2f_{xy}m+f_{xx}=0$, $m=dy/dx$, so a regular saddle point of a scalar function crossing has $f_{xy}^2-f_{xx}{f}_{yy}>0$. These are the form of the slow and fast regularity conditions in a warm wind, with the enthalpy contribution retained in $f$.

In the strictly cold approximation used to obtain the displayed $f$, the slow speed is zero. There is a cold-limit degeneracy of a slow magnetosonic point: a finite-speed interior slow crossing need not survive. Indeed,so this cold equation has no ordinary sub-Alfvénic slow critical point. A small retained enthalpy is needed to locate that point, which can approach the launching boundary as the sound speed tends to zero. The fast critical point still obeys the simultaneous conditions above wherever the reduced expression is nonsingular.

Introduce the squared poloidal Alfvén numberSolving the two linear azimuthal invariants givesAt an Alfvén surface, $M_A^2=1$. Smooth passage through the apparent singularity requires both numerators to vanish at the same cylindrical radius $r_A$, giving the Alfvén-surface regularity condition for an axisymmetric windThe finite limiting values of $u_{\varphi }$ and $B_{\varphi }$ then follow by l'Hopital's rule from the local variation of $M_A$ and $r$ along the field line; the algebraic invariants alone fix their combination rather than each value separately at the critical point.