2004/06/30 by Chi-kwan Chan, Chi‐kwan Chan, Dimitrios Psaltis +1
Physics and Astronomy · #Accretion (finance) #Accretion disc #Advection #Astrophysical Phenomena and Observations #Astrophysics and Star Formation Studies #Boundary (topology) #Boundary value problem #Focus (optics) #Numerical analysis #Pulsars and Gravitational Waves Research #Spectral method #Stability (learning theory) #astro-ph
paper · pdf · doi:10.1086/430511
published as Astrophys.J. 628 (2005) 353-367 · 15 pages, 11 figures. To appear in the ApJ. Substantially revised to include the energy equation and an application to fallback disks. High resolution plots and animations of the simulations are available at http://www.physics.arizona.edu/~chan/research/astro-ph/0406073/index.html
arxiv created 2005/03/26 · openalex publication_date 2005/07/18 · openalex created_date 2016/06/24 · openalex updated_date 2026/08/05
We present a numerical method for studying the normal modes of accretion flows around black holes. In this first paper, we focus on two-dimensional, viscous, hydrodynamic disks, for which the linear modes have been calculated analytically in previous investigations. We use pseudospectral methods and low-storage Runge-Kutta methods to solve the continuity equation, the Navier-Stokes equation, and the energy equation. We devise a number of test problems to verify the implementation. These tests demonstrate the ability of spectral methods to handle accurately advection problems and to reproduce correctly the stability criteria for differentially rotating hydrodynamic flows. They also show that our implementation is able to handle sound waves correctly with nonreflective boundary conditions, to recover the standard solution for a viscous-spreading ring, and to produce correctly the Shakura-Sunyaev steady-disk solution. Finally, we have applied our algorithm to the problem of a nonaxisymmetric viscous-spreading ring and verify that such configuration is unstable to nonaxisymmetric perturbations.