Morphology and Mach Number Dependence of Subsonic Bondi-Hoyle Accretion

We carry out three-dimensional computations of the accretion rate onto an object (of size R sink and mass m) as it moves through a uniform medium at a subsonic speed v infinity. The object is treated as a fully absorbing boundary (e.g., a black hole). In contrast to early conjectures, we show that for an accretor with R sink << R A = 2 Gm / v infinity 2 in a gaseous medium with adiabatic index gamma = 5/3, the accretion rate is independent of Mach number and is determined only by m and the gas entropy. Our numerical simulations are conducted using two different numerical schemes via the Athena++ and Arepo hydrodynamics solvers, which reach nearly identical steady-state solutions. We find that pressure gradients generated by the isentropic compression of the flow near the accretor are sufficient to suspend much of the surrounding gas in a near-hydrostatic equilibrium, just as predicted from the spherical Bondi-Hoyle calculation. Indeed, the accretion rates for steady flow match the Bondi-Hoyle rate, and are indicative of isentropic flow for subsonic motion where no shocks occur. We also find that the accretion drag may be predicted using the Safronov number, Theta = R A /R sink, and is much less than the dynamical friction for sufficiently small accretors (R sink << R A ).