Past exam of the mathematics course of the University of Cambridge / 2012 / iii / Paper 62 (Astrophysical Fluid Dynamics)

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.

Let $u=u_r<0$ for accretion. The steady spherical continuity equation and radial Euler equations for an inviscid fluid giveFor a barotropic fluid, $dp/dr=c_s^{2}d\rho /dr$. Dividing continuity equation by $r^2\rho u$ gives ${\rho }^{\prime }/\rho =-2/r-u^{\prime }/u$. Eliminating this derivative yieldsAt a smooth transonic critical point, the coefficient of $du/dr$ vanishes. A finite derivative then requires the right side to vanish simultaneously:These are the sonic and regularity conditions. Otherwise the derivative is singular and the proposed smooth transonic passage fails. The local slope must also be a real root compatible with the desired branch. For example, defining $C=\rho d(c_s^2)/d\rho$ at the point and differentiating both sides givesThis explains why simultaneous vanishing is a necessary regularity condition, rather than an automatic proof that every potential admits a critical crossing.

Define the barotropic enthalpy function by $dh=dp/\rho$. The momentum equation integrates to the Bernoulli equationChanging the reference mass density in $h$ changes only the constant $B$. For an isothermal closure this logarithmic function is a barotropic barotropic pressure potential; it should not be confused with the constant thermodynamic specific enthalpy of an ideal gas held at fixed temperature.

For isothermal transonic accretion in the Paczyński-Wiita potential, assume $M>0$ and $c_s>0$, and define $r_B=GM/(2c_s^2)$ and $\beta =16c_s^2/c^2$. The critical equation iswhose two roots are $r_G+\frac{1}{2}{r}_B(1\pm \sqrt{1+\beta })$. Only the plus root lies outside $r_G$. ThereforeFor this isothermal case $C=0$, and the slope relation reduces toThe positive slope is the inward-accretion solution connecting a small inward speed at large radius to the supersonic inward branch. Thus the exterior critical point is a genuine nondegenerate transonic point.

Choose the reference mass density as ${\rho }_0$, so $h=c_s^2\log (\rho /{\rho }_0)$ and $B=0$ from the conditions at infinity. At $u_c=-c_s$,Substitution into $\dot{M}=4\pi r_c^2{\rho }_c{c}_s$ givesThis is the rate selected by the smooth transonic accretion solution. Arbitrary static or subsonic formal solutions are not assigned this rate merely by specifying conditions at infinity. As $r_G/r_B\to 0$, the result tends to the isothermal Bondi accretion rate $\pi e^{3/2}{\rho }_0(GM{)}^2/c_s^3$, a useful normalization check.

The flux representation gives $B_R=-{\psi }_z/R$ and $B_z={\psi }_R/R$. Therefore $\nabla \cdot \mathbf{B}=0$ and $\mathbf{B}\cdot \nabla \psi =0$. With magnetohydrodynamic mass loading $k=k(\psi )$,Thus the steady continuity equation is satisfied.

The cross product in the ideal magnetohydrodynamic induction equation simplifies toHence steady induction requires $\mathbf{B}\cdot \nabla (v/R)=0$. On a regular connected poloidal flux surface, this integrates locally toThis is the field-line angular velocity of an axisymmetric wind. It need not equal the gas angular velocity when $B_{\varphi }\ne 0$, because $u_{\varphi }=R\omega +k{B}_{\varphi }/\rho$. The regular-surface qualification matters: if $\psi$ is constant and the field is purely toroidal, induction alone cannot turn an arbitrary $v(R,z)/R$ into a function of $\psi$.

Put $\mathbf{F}=(\nabla \times \mathbf{B})\times \mathbf{B}$ and $E=|\mathbf{u}{|}^2/2+\Phi$. Dot the pressure-free momentum equation with $\mathbf{u}$. Since $\mathbf{B}\cdot \mathbf{F}=0$, axisymmetry givesThe azimuthal momentum equation in cylindrical coordinates similarly givesSubtracting $\omega$ times the second relation from the first, and using $\mathbf{B}\cdot \nabla \omega =0$, proves that $E-\omega R{u}_{\varphi }$ is constant along a flowing regular poloidal flux surface. For nonzero mass loading its corotating energy invariant of an axisymmetric magnetic wind isThis is the invariant labeled $ϵ$ in the question. It is a corotating mechanical energy, distinct from the total magnetohydrodynamic Bernoulli invariant, which includes electromagnetic energy flux. In standard notation it equals total specific energy minus $\omega$ times the total specific angular-momentum invariant.

Near the footpoint, negligible $k{B}_{\varphi }/\rho$ gives $u_{\varphi }\simeq R\omega$, while the poloidal velocity is parallel to the poloidal magnetic field. The gas approximately corotates and slides along a rotating field line as a bead on a wire, with effective potentialA decreasing effective potential converts corotating potential energy into poloidal kinetic energy. This is the mechanism of magnetocentrifugal acceleration; local launching does not alone guarantee global escape.

For an attractive point mass, take $\Phi =-GM/(R^2+z^2{)}^{1/2}$ and ${\omega }^2=GM/R_0^3$. Parameterize the outward straight line by $R=R_0+s\sin \alpha$, $z=s\cos \alpha$. The first derivative of ${\Phi }_{\mathrm{eff}}$ vanishes at $s=0$ because the footpoint is in circular radial balance. The meridional Hessian there is ${\omega }^2\mathrm{diag}(-3,1)$, soThus the thirty-degree magnetocentrifugal launching criterion in its strict negative-curvature form isAbove this angle the equilibrium is unstable to an outward displacement; below it the near-footpoint potential initially rises. The source's proportionality constant must be negative to describe attraction.

The marginal straight-line launch at thirty degrees has zero quadratic curvature. For precisely the straight line and point-mass potential here, let $x=s/R_0$ and $a=\sin \alpha$. Direct expansion givesAt $\alpha =\pi /6$, this is $-7x^3/16+O(x^4)$, so the outward direction is downhill at cubic order. Consequently the strict inequality is the nondegenerate quadratic criterion; it should not be interpreted as excluding every one-sided nonlinear marginal launch at equality. An exactly resting bead remains an equilibrium until displaced.

In the inertial frame the steady convective acceleration is $-R{\Omega }^2{\mathbf{e}}_R$. For a barotropic fluid, $\nabla p/\rho =\nabla h$, where $h$ is the barotropic enthalpy function. The equilibrium momentum equation therefore giveson each connected fluid region. Uniform rotation is essential: it makes the centrifugal acceleration derivable from the centrifugal potential $-{\Omega }^2{R}^2/2$.

The Cowling approximation sets the gravitational-potential perturbation to zero while retaining background gravity. The barotropic pressure force linearizes as $\delta (\nabla p/\rho )=\nabla \delta h$, with $\delta h=p^{\prime }/\rho =c_s^2{\rho }^{\prime }/\rho$. Denote its mode amplitude by $W$. The advective time derivative on a scalar mode is ${\partial }_t+\Omega {\partial }_{\varphi }$, so the stated positive-frequency exponential produces $i\sigma$ with $\sigma =\omega +m\Omega$. Linearizing the radial centrifugal term produces $-2\Omega v_{\varphi }$, while linearizing the azimuthal convective term produces $2\Omega v_R$. ThusThese are the pressure-gradient and Coriolis acceleration terms in the corotating-frame form. Independently, linearizing continuity equation givesThe equilibrium mass density is axisymmetric, so no additional azimuthal background-density derivative appears.

For $\sigma \ne 0$ and $D=4{\Omega }^2-{\sigma }^2\ne 0$, invert the horizontal momentum system:Inserting these into the continuity equation, dividing by $i\sigma$ and multiplying by $D$ givesThe two terms proportional to $\rho W_R$ cancel; this cancellation leaves precisely the background-density derivative in the last term. Since $-D/{\sigma }^2=1-4{\Omega }^2/{\sigma }^2$, this is exactlyThe singular frequency cases require the original velocity equations rather than this inversion.

In the low-frequency anelastic approximation for a rotating barotropic star, neglect the left side. Spherical background mass density satisfies ${\rho }_R/R={\rho }_z/z$, understood by its smooth limiting form on the equatorial plane. For $W=z{R}^m$,The remaining equation reduces toWriting $q=2\Omega /\sigma$, the bracket factors as $(m+1-q)(q+1)$. The regular desired branch is thereforeThe other algebraic factor corresponds to $\sigma =-2\Omega$, where the eliminated horizontal system is singular; it is not classified by this pressure-equation inversion. A constant-density background makes the bulk factor identically zero but does not invalidate the displayed solution.

For the regular branch, an explicit velocity check is particularly informative. Up to a common mode normalization,It satisfies $R{v}_R+z{v}_z=0$ and $\nabla \cdot \mathbf{v}=0$. Thus $\nabla \cdot (\rho \mathbf{v})=0$ for every spherical mass density profile, directly verifying the anelastic continuity condition. The motion is tangential to spherical shells and is the sectoral inertial mode of a slowly rotating barotropic star. For single-valued azimuthal modes, $m$ is a positive integer. With the source's exponential convention,For $m\ge 2$ the pattern is retrograde relative to the star but prograde in the inertial frame. This is a leading slow-rotation anelastic mode, not an exact solution of the compressible equations whose left-hand term was discarded.