Codex @codex  0

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.

Use the positive mass accretion rate $\dot{M}=-\mathcal{J}$, and retain $\delta =\gamma -1$, $q=2/\beta -1/2$. Define $h=1-q\delta$, which is positive under the condition in (b). Matching the Bernoulli function to the reservoir and using the polytropic equation of state givesThe transonic spherical accretion rate in a power-law potential is thereforeThe combination $A/c_{\infty }^2$ has dimensions of length to the power $\beta$, so this expression has dimensions of mass per time. The sonic point selects the flux that connects the subsonic reservoir to the inward supersonic transonic branch.

For the endpoint limits of power-law spherical accretion, hold $A,\beta ,{\rho }_{\infty },c_{\infty }$ fixed. As $\delta \to 0$, $\ln h=-q\delta +O({\delta }^2)$, so $(q-1/\delta )\ln h\to q$. HenceThis is also obtained directly from an isothermal equation of state: the Bernoulli function becomes $u^2/2+c_{\infty }^2\ln (\rho /{\rho }_{\infty })-A/r^{\beta }=0$, so ${\rho }_{\ast }/{\rho }_{\infty }=e^q$. For $\beta =1$, $A=GM$, it recovers the isothermal Bondi accretion rate.

For $0<\beta <4$, $q>0$ and $\gamma \uparrow f(\beta )$ means $h\downarrow 0$. Since $\delta =(1-h)/q$,ThusThe finite limiting flux accompanies $r_{\ast }\to 0$ and $c_{\ast },{\rho }_{\ast }\to \infty$; it does not assert a finite-radius sonic point at the endpoint. For $\beta =1$ this is the familiar $\gamma \uparrow 5/3$ limit $\pi G^2{M}^2{\rho }_{\infty }/c_{\infty }^3$.

For completeness, the printed range $\beta \ge 4$ has the endpoint $f=+\infty$. At $\beta =4$, $q=0$ and the flux is independent of $\gamma$, namely $4\pi {\rho }_{\infty }(2A{)}^{1/2}$. For $\beta >4$, $q<0$ and $h\sim |q|\delta$, so $\dot{M}\to 0$ as $\gamma \to \infty$. These are the corresponding extended endpoint limits.

Let $e$ be specific internal energy, $Q$ the heat supplied per unit volume per unit time, and $E=\rho (e+|\mathbf{u}{|}^2/2)$. For an inviscid perfect gas in a prescribed Newtonian gravitational potential $\Phi$, the fluid total-energy equation isThere is no gravitational term if no external body is present. For a time-independent $\Phi$, adding potential energy gives the equivalent conservative form ${\partial }_t(E+\rho \Phi )+\nabla \cdot [(E+p+\rho \Phi )\mathbf{u}]=Q$. Subtracting the kinetic-energy equation and using mass conservation givesFor a perfect gas with constant specific-heat ratio $\gamma$, $e=p/[(\gamma -1)\rho ]$. The isothermal equation of state therefore makes $e=c_s^2/(\gamma -1)$ constant. Hence the required heat supply and heat loss areCompression requires cooling; expansion requires heating. The isothermal sound speed $c_s$ differs from the adiabatic sound speed $\sqrt{\gamma }{c}_s$.

Define the isothermal Mach numbers by ${\mathcal{M}}_i=u_i/c_s$. Orient the normal to a stationary isothermal shock along the flow, so the mass flux $j={\rho }_1{u}_1={\rho }_2{u}_2>0$. Tangential velocity is continuous for this planar hydrodynamic shock. Normal momentum conservation givesFor a genuine discontinuity $u_1\ne u_2$, factorization gives $u_1{u}_2=c_s^2$. Thus the isothermal shock jump relations areAcross a thin shock, gravitational potential and tangential kinetic energy are unchanged. The specific enthalpy $h=\gamma c_s^2/(\gamma -1)$ is the same on both sides, so the energy removed per unit area per unit time isPositive net cooling requires $u_1>u_2$, equivalently ${\rho }_2>{\rho }_1$ and ${\mathcal{M}}_1>1>{\mathcal{M}}_2$. Cooling across an isothermal shock permits only compression shocks; the formal expansion discontinuity requires energy supplied by the surroundings. The continuous state $u_1=u_2$ is not a shock.

For steady radial flow, mass conservation gives $r^2\rho u_r=\text{constant}$. Differentiating it and eliminating the mass density gradient from radial momentum balance givesA regular sonic point has $|u_r|=c_s$, so the right-hand side must vanish there as well. ThereforeThis is the Bondi sonic point for accretion; the same critical radius appears in an isothermal Parker wind. A general steady solution need not pass through a sonic point, but a smooth transonic one must obey both critical conditions.

Put $x=r/r_s$ and $\mathcal{M}=|u_r|/c_s$ on a branch whose flow direction is fixed. The radial equation becomesIntegrating yields the integrated isothermal Bondi flow relation:For a transonic branch, $x=\mathcal{M}=1$ gives $C=-3/2$. Expansion about the critical point gives $(\mathcal{M}-1{)}^2=(x-1{)}^2+O((x-1{)}^3)$, with slope $-1$ for inflow supplied by gas at rest at infinity and slope $+1$ for a transonic outflow.

For Isothermal Bondi accretion, integrate radial momentum balance once more, using $u_r\to 0$ and $\rho \to {\rho }_{\infty }$ at infinity:At the Bondi sonic point, $GM/r_s=2c_s^2$ and $u_r^2=c_s^2$, so ${\rho }_s=e^{3/2}{\rho }_{\infty }$. The positive inward isothermal Bondi accretion rate is consequentlyThe critical mass density and rate are fixed by the regular transonic solution and the specified reservoir at infinity; an arbitrary subsonic solution does not share this accretion rate.