Numerical Computation of Two Dimensional Wind Accretion of Isothermal   Gas
  A new numerical algorithm for calculating isothermal wind accretion flows has been developed and is applied here to the analysis of the hydrodynamics of two-dimensional plane symmetric accretion flows in wind-fed sources. Polar coordinates are used to ensure fine resolution near the object. It is found that a thin accretion column is formed which shows wave-like oscillations. Small accretion disks are formed temporarily around the object. Mass accretion rate and angular momentum accretion rate exhibit quasi-periodic oscillations. The amplitudes of the oscillations depend on the size of the inner boundary, the number of grid points and the method of calculation. For a smaller size of the inner boundary, finer grids and more accurate numerical schemes, the amplitudes of the oscillation become larger. 
Eiji Shima, Takuya Matsuda, Ulrich Anzer, Gerhard Borner and Henri   M.J. Boffin