We present a numerical scheme for viscoelastic flow based on a second-order central differencing method recently introduced in the context of incompressible newtonian flow; the incompressibility constraint is treated with the projection method. The result is a simple and efficient scheme that is readily adaptable to a wide class of differential constitutive equations and flow geometries. We imp...

In this paper, we discuss a two-dimensional numerical study of isothermal high Mach number Bondi– Hoyle–Lyttleton flow around a small accretor. The flow is found to be unstable at high Mach numbers, with the instability appearing even for a larger accretor. The instability appears to be the unstable radial mode of the accretion column predicted by earlier analytic work.