Spectral Methods for Time-Dependent Studies of Accretion Flows. II. Two-Dimensional Hydrodynamic Disks with Self-Gravity

Abstract

Spectral methods are well suited for solving hydrodynamic problems in which the self-gravity of the flow needs to be considered. Because Poisson's equation is linear, the numerical solution for the gravitational potential for each individual mode of the density can be pre-computed, thus reducing substantially the computational cost of the method. In this second paper, we describe two different approaches to computing the gravitational field of a two-dimensional flow with pseudo-spectral methods. For situations in which the density profile is independent of the third coordinate (i.e., an infinite cylinder), we use a standard Poisson solver in spectral space. On the other hand, for situations in which the density profile is a delta function along the third coordinate (i.e., an infinitesimally thin disk), or any other function known a priori, we perform a direct integration of Poisson's equation using a Green's functions approach. We devise a number of test problems to verify the implementations of these two methods. Finally, we use our method to study the stability of polytropic, self-gravitating disks. We find that, when the polytropic index Gamma is <= 4/3, Toomre's criterion correctly describes the stability of the disk. However, when Gamma > 4/3 and for large values of the polytropic constant K, the numerical solutions are always stable, even when the linear criterion predicts the contrary. We show that, in the latter case, the minimum wavelength of the unstable modes is larger than the extent of the unstable region and hence the local linear analysis is inapplicable.

13 pages, 9 figures. To appear in the ApJ. High resolution plots and animations of the simulations are available at http://www.physics.arizona.edu/~chan/research/astro-ph/0512448/index.html

References in corpus (7)

• On the Dynamical Foundations of Alpha Disks
• The Thermal Regulation of Gravitational Instabilities in Protoplanetary Disks II. Extended Simulations with Varied Cooling Rates
• "Mariage des Maillages": A new numerical approach for 3D relativistic core collapse simulations
• A multi-domain spectral method for scalar and vectorial Poisson equations with non-compact sources
• Thermal stability of self-gravitating, optically thin accretion disks
• Self-regulated Accretion Disks
• Spectral Methods for Time-Dependent Studies of Accretion Flows. I. Two-dimensional, Viscous, Hydrodynamic Disks

Cited by in corpus (8)

• On the time variability of geometrically-thin black hole accretion disks I : the search for modes in simulated disks
• MHD Simulations of Accretion onto Sgr A*: Quiescent Fluctuations, Outbursts, and Quasi-Periodicity
• Studies of Thermally Unstable Accretion Disks around Black Holes with Adaptive Pseudo-Spectral Domain Decomposition Method I. Limit-Cycle Behavior in the Case of Moderate Viscosity
• A Fast Potential and Self-Gravity Solver for Non-Axisymmetric Disks
• Self-Gravitational Force Calculation of Infinitesimally Thin Gaseous Disks on Nested Grids
• Multi-scale simulations of black hole accretion in barred galaxies: Self-gravitating disk models
• Spectral Methods for Time-Dependent Studies of Accretion Flows. III. Three-Dimensional, Self-gravitating, Magnetohydrodynamic Disks
• Self-gravity in thin protoplanetary discs: 1. The smoothing-length approximation versus the exact self-gravity kernel