Numerical Models of Viscous Accretion Flows Near Black Holes
  We report on a numerical study of viscous fluid accretion onto a black hole. The flow is axisymmetric and uses a pseudo-Newtonian potential to model relativistic effects near the event horizon. The numerical method is a variant of the ZEUS code. As a test of our numerical scheme, we are able to reproduce results from earlier, similar work by Igumenshchev and Abramowicz and Stone et al. We consider models in which mass is injected onto the grid as well as models in which an initial equilibrium torus is accreted. In each model we measure three ``eigenvalues'' of the flow: the accretion rate of mass, angular momentum, and energy. We find that the eigenvalues are sensitive to r_{in}, the location of the inner radial boundary. Only when the flow is always supersonic on the inner boundary are the eigenvalues insensitive to small changes in r_{in}. We also report on the sensitivity of the results to other numerical parameters. 
Jonathan C. McKinney and Charles F. Gammie