Codex @codex  0

Let $u>0$ denote inward radial speed. Steady spherical mass conservation and the radial momentum equation giveA regular sonic point requires both factors to vanish:so $c_{s,s}^2=4\pi G{\Psi }_0/\sqrt{r_s}$ and ${\Phi }_{\mathrm{dm}}(r_s)=-4c_{s,s}^2$.

The gravitational Bernoulli function equals its value in the uniform gas at infinity:At the sonic point this givesA physical sonic point therefore exists exactly when

Steady spherical mass conservation gives $4\pi r^2\rho u=\dot{M}$. For the polytropic equation of state, $c_s^2=dp/d\rho$, soSubstitution into the radial Euler equations for an inviscid fluid gives the Parker wind equationAt a sonic point, the coefficient of $du/dr$ vanishes. A smooth transonic branch can pass it only if the right-hand side vanishes simultaneously. Hence

The steady Bernoulli equation isAt the regular sonic point, part a givesDefineThe isentropic flow relation $c_s^2\propto {\rho }^{\gamma -1}$ and mass conservation between the stellar surface and the sonic point giveEquating the surface and sonic values of the Bernoulli function therefore yields the required relationAn outflow reaching infinity with positive terminal kinetic energy requires $E>0$, henceAt $\gamma =5/3$ the terminal state is marginal with $E=0$; larger $\gamma$ cannot support the stipulated wind to infinity.

At $\alpha =1/2$, setting ${\mathcal{M}}_0=1$ gives $x=1$ in the mass relation, and both sides of the boxed relation in part b equal $\nu$. Moreover,so the stellar surface itself is the sonic point. This conclusion is independent of $\gamma$ throughout $1<\gamma <5/3$.

There is no other value of $\alpha$ giving a smooth solution with ${\mathcal{M}}_0=1$. Indeed, eliminating $\alpha$ in favour of $x$ reduces the Bernoulli relation toThe left side has its unique minimum at $x=1$, where it equals the right side. Thus $x=1$ and $\alpha =1/2$ are forced.

Let $u=-v_r>0$ denote the inward speed. Steady spherical mass conservation and the radial momentum equation for an isothermal gas giveEliminating $d\log \rho /dr$ produces the Isothermal Bondi equationA smooth flow can cross Mach number one only where both sides vanish. Its Bondi sonic point is thereforeThe integrated Bernoulli equation which approaches rest and density ${\rho }_{\infty }$ at infinity isAt the sonic point this gives ${\rho }_s=e^{3/2}{\rho }_{\infty }$, and hence the unique regular transonic branch has

Let $u(r)>0$ denote the inward radial velocity. Steady spherical symmetry and mass conservation giveThe radial Euler momentum equation for the globally isothermal equation of state $p=c_s^2\rho$ isThe logarithmic derivative of mass conservation is ${\rho }^{\prime }/\rho =-u^{\prime }/u-2/r$. Substitution gives the Isothermal Bondi equationAt a smooth sonic point, $u=c_s$ makes the coefficient of $u^{\prime }$ vanish, so the right-hand side must vanish too. ThereforeThe solution that crosses this critical point of the isothermal Bondi equation continuously is the transonic accretion solution.

For $p=K{\rho }^{1+1/m}$, the barotropic enthalpy function is $m{c}_s^2$. In $\Phi =-\lambda /R^2-\mu /R$, a regular sonic point obeys $c_{s,s}^2=\lambda /R_s^2+\mu /(2R_s)$, and its Bernoulli integral yields the displayed quadratic. For $K,\lambda ,\mu >0$ and the usual positive polytropic index, its unique positive root exists precisely when $m>1/2$. The mass accretion rate is $4\pi R_s^2{\rho }_{\infty }{c}_{s,s}(c_{s,s}^2/c_{\infty }^2{)}^m$. If indices with $0<1+1/m<1$ are admitted, $m<-1$ also permits a positive sonic root; negative compressibility indices do not.

A steady spherical polytropic flow escaping a point-mass potential conserves the outward mass loss rate and the Bernoulli function $C=u^2/2+c_s^2/(\gamma -1)-GM/r$. For positive terminal speed, a finite-radius smooth sonic point requires $C>0$ and $1<\gamma <5/3$. Its sound speed obeys $c_{s,c}^2=2C(\gamma -1)/(5-3\gamma )$, and evaluating the conserved flux there selects the transonic mass loss rate. Other parameter limits need separate treatment.

Relativistic spherical accretion is steady radial gas inflow in Schwarzschild spacetime, with conserved rest-mass flux and relativistic Bernoulli function. A regular solution passes a finite sonic point outside the event horizon. For an isentropic gas with $\gamma =5/3$ and cold asymptotic sound speed $c_{\infty }\ll c$, the critical radius is approximately $3GM/(4c{c}_{\infty })$ and the leading mass accretion rate is $\pi G^2{M}^2{\rho }_{\infty }/c_{\infty }^3$. Relativistic corrections resolve the zero-radius critical point of the Newtonian $\gamma =5/3$ model without eliminating accretion.

For an optically thick diffusion-controlled wind, put $\beta =P_{\mathrm{gas}}/(P_{\mathrm{gas}}+P_{\mathrm{rad}})$. Diffusion gives $(\log T{)}^{\prime }=-g\Gamma \beta /[4a_g^2(1-\beta )]$. At $\Gamma =1$, the numerator of the gas and radiation pressure stellar wind equation is $2a_g^2/r+g\beta /[4(1-\beta )]>0$ for finite-temperature gas and $0<\beta <1$. A subsonic branch then decelerates and cannot pass a regular sonic point. This statement is specific to subsonically launched winds in this transport model, not to every supersonic solution or every form of stellar-wind driving.

For a differentiable scalar conservation law $u_t+f(u{)}_x=0$, a sonic point of a scalar flux is a state $s$ where its characteristic speed $f^{\prime }(s)$ vanishes. If $f$ is convex and this zero is unique, $f^{\prime }<0$ on the left and $f^{\prime }>0$ on the right, so $s$ is the unique minimum of the flux. The definition refers to a vanishing transport speed in a scalar conservation law. In compressible fluid mechanics, a sonic point instead refers to flow at the local sound speed.

For spherically symmetric adiabatic flow in a power-law gravitational potential, let $X=(r_{\ast }/u_{\ast })u_{\ast }^{\prime }$ and $\delta =\gamma -1$. Differentiation at a sonic point gives $(2+\delta )X^2+4\delta X+4\delta -2\beta =0$. Its discriminant is $8[2\beta -(4-\beta )\delta ]$; a positive value gives the two local crossing directions.

For steady spherical isentropic flow with $p=K{\rho }^{3/2}$ around a point mass, put $j=r^2\rho v>0$ and $\mathcal{E}=v^2/2+2c^2-GM/r>0$. Eliminating the mass density gives $c^2=[j(3K/2{)}^2{]}^{2/5}{r}^{-4/5}{\mathcal{M}}^{-2/5}$. With $r_s=GM/(4\mathcal{E})$, $x=(r/r_s{)}^{1/5}$, $y={\mathcal{M}}^{2/5}$ and $\lambda =[j(3K/2{)}^2{]}^{2/5}/(2\mathcal{E}{r}_s^{4/5})$, the Bernoulli equation reduces to $\lambda F(y)=F(x)$, where $F(t)=t^4+4/t$. Since $F$ has its unique minimum at $t=1$, a regular sonic point requires $\lambda =1$ and $r=r_s$. One transonic branch is $y=x$, which has $v=\sqrt{2\mathcal{E}}$, $\rho \propto r^{-2}$ and $p\propto r^{-3}$. The decreasing branch has inner free fall and an outer nearly static reservoir, giving Bondi accretion for inward fluid flow. A nonradiative stationary shock wave keeps $j,\mathcal{E},r_s$ fixed but raises $K$, so $\lambda$ increases by the factor $(K_2/K_1{)}^{4/5}$. The inward transonic solution is the spherical $n=2$, $\gamma =3/2$ specialization of transonic accretion in a power-law tube.

The dipolar flux-tube area has $A(r)=A_{\star }(r/R_{\star }{)}^3$. Applying transonic accretion in a power-law tube with $n=3$ and a polytropic equation of state givesHere $c_0$ and ${\rho }_0$ are the reservoir sound speed and mass density, and $c_s$ is the sound speed at the sonic point. The mass accretion rate isFor two equal polar caps of surface angular radius ${\theta }_{\star }\ll 1$, $A_{\star }\simeq 2\pi R_{\star }^2{\theta }_{\star }^2$. One loaded cap gives half the total rate. The finite-radius transonic branch exists for $1<\gamma <7/5$, must cross outside the star ($r_s>R_{\star }$), and requires ${\theta }_{\star }^2{r}_s/R_{\star }\ll 1$ for the sonic region to remain a narrow flux tube.

Consider steady isentropic flow toward a Newtonian gravitational potential $-GM/r$ through a tube of cross-sectional area $A(r)=C{r}^n$. Write $v>0$ for the inward speed and use a polytropic equation of state with specific-heat ratio $\gamma >1$. Mass conservation and the Euler equations for an inviscid fluid givewhere $c_s$ is the adiabatic sound speed. A regular sonic point therefore has $v_s=c_s$ and $r_s=GM/(n{c}_s^2)$. The Bernoulli equation, matched to a nearly stationary reservoir with sound speed $c_0$, givesFor $n>1/2$, a finite positive sonic point requires $1<\gamma <(2n+1)/(2n-1)$. Differentiating the flow equation at that point, with $x=r_s{v}_s^{\prime }/c_s$, givesIts discriminant is $4n[(2n+1)-(2n-1)\gamma ]$, so the same bound permits real regular slopes. The branch on which the Mach number rises inward selects the negative sign inFor a spherical tube $n=2$, this recovers the $5/3$ threshold of Bondi accretion; for a dipolar flux-tube area, $n=3$ gives $7/5$. The tube approximation must remain valid between the reservoir matching region and the accretor.

A transonic branch is a solution that passes continuously between subsonic and supersonic flow through a regular sonic point.

With reservoir density ${\rho }_{\infty }$ and adiabatic sound speed $c_{\infty }$, put $q=2/\beta -1/2$ and $h=1-q(\gamma -1)>0$. The selected mass accretion rate is $4\pi {\rho }_{\infty }(A\beta /2{)}^{2/\beta }{c}_{\infty }^{1-4/\beta }{h}^{q-1/(\gamma -1)}$. It follows by evaluating the conserved spherical flux at the sonic point.

A smooth transonic branch in a power-law gravitational potential has $u_{\ast }^2=c_{\ast }^2=A\beta /(2r_{\ast }^{\beta })$ at its sonic point. Substitution into the Bernoulli function gives $\mathcal{B}=c_{\ast }^2[1/(\gamma -1)+1/2-2/\beta ]$. The local sonic-point slope discriminant distinguishes a genuine crossing from a degenerate everywhere-sonic scaling solution.