Simulations from Waters et al. (2019)
The top panels show the inflow (blue lines) and outflow (red lines) rates as a function of radius. The difference between these curves is the net accretion rate, shown as dashed black lines; it shows barely noticeable fluctuations away from the Bondi rate as the simulation evolves.
The middle panels show colormaps of pressure (in units of
p_0, the reference pressure at the outer boundary) on a linear scale over the full domain. The bottom panels zoom into the region within the Bondi radius and show instead log-scale plots of the azimuthal component of angular momentum in units of
R_s c. This region looks nearly identical to the spherically symmetric Bondi solution when looking at density, radial velocity, or pressure; the sonic surface is unchanged and appears in the bottom left corner as a black contour.
The initial conditions for all of these runs is the analytic Bondi solution. Time is shown in the upper right corner of each plot in units of the dynamical time
Constant entropy BC runs
Constant pressure BC runs
These simulations were run on the Institutional Computing clusters at Los Alamos National Laboratory under allocation award w19_rhdccasims (PI: T. Waters). We thank the high performance computing team for their continued efforts, as well as the developers of the public GRMHD code Athena++.