Subcritical, subsonic case

General Soln Behavior for this case (different ICs than this run, so just to qualitatively compare):

*Note vr is actually negative for infall

Initial conditions are here: InitialConditions.pdf​

128³, fixed grid

At this resolution (dx=.03125), the critical radius (rcrit=.0158) is contained within the kernel (r_acc=.125), but below the grid scale (r_min=.027). So, the soln doesn't turn over the way it does in the above plots, at this resolution. Instead, the density and radial velocity are monotonic.

A density comparison of the new algorithm (left) vs. the Krumholz algorithm (right):

Density lineouts comparison-

Mach comparison-

Sink mass vs. time-