Isothermal Bondi accretion

Isothermal Bondi accretion has $p=c_s^2\rho$ with constant sound speed. Its regular transonic solution crosses the Bondi sonic point and approaches free fall near the accretor.

For steady Isothermal Bondi accretion from gas at rest at infinity, define $\lambda =\dot{M}{c}_s^3/(4\pi G^2{M}^2{\rho }_{\infty })$, $x=r/r_s$ and $y=u/c_s>0$, where $r_s=GM/(2c_s^2)$. The Bernoulli function and mass conservation give $y^2/2-\ln y=2\ln x+2/x-\ln (4\lambda )$. The left side has minimum $1/2$ at $y=1$, and the radius-dependent part on the right has minimum two at $x=1$. Existence over all positive radii therefore requires the displayed bound. Equality permits the smooth accelerating-inward transonic branch and gives the usual isothermal Bondi accretion rate. Smaller positive rates have an everywhere-subsonic reservoir-connected branch, and zero rate gives a hydrostatic atmosphere. An absorbing inner boundary inside the sonic radius selects the critical inflow; steady isothermal equations and the outer density alone do not specify that inner condition.

The regular spherical isothermal inflow from mass density ${\rho }_{\infty }$ at rest at infinity has sonic mass density $e^{3/2}{\rho }_{\infty }$. Combining it with the Bondi sonic point and mass conservation gives $\dot{M}=\pi e^{3/2}{G}^2{M}^2{\rho }_{\infty }/c_s^3$.

For inward speed $u>0$, steady spherical mass conservation and radial momentum balance give

With $x=r/r_s$ and $\mathcal{M}=|u_r|/c_s$, the radial isothermal equation integrates to ${\mathcal{M}}^2/2-\ln \mathcal{M}=2/x+2\ln x+C$. A regular transonic branch has $C=-3/2$, and the reservoir-fed accretion branch crosses $(x,\mathcal{M})=(1,1)$ with slope $-1$.

The isothermal Bondi sonic point is $r_s=GM/(2c_s^2)$, where the regular transonic solution has $u=c_s$.

At the critical point, both the coefficient of $du/dr$ and the remaining right-hand side of the Isothermal Bondi equation vanish, allowing a finite derivative through the sonic transition.