Misc,

VADER: A Flexible, Robust, Open-Source Code for Simulating Viscous Thin Accretion Disks

, and .
(Jun 25, 2014)

Abstract

The evolution of thin axisymmetric viscous accretion disks is a classic problem in astrophysics. While such models provide only approximations to the true processes of instability-driven mass and angular momentum transport, their simplicity makes them invaluable tools for both semi-analytic modeling and simulations of long-term evolution where two- or three-dimensional calculations are too computationally costly. Despite the utility of these models, there is no publicly-available framework for simulating them. Here we describe a highly flexible, general numerical method for simulating viscous thin disks with arbitrary rotation curves, viscosities, boundary conditions, grid spacings, equations of state, and rates of gain or loss of mass (e.g., through winds) and energy (e.g., through radiation). Our method is based on a conservative, finite-volume, second-order accurate discretization of the equations, which we solve using an unconditionally-stable implicit scheme. We implement Anderson acceleration to speed convergence of the scheme, and show that this leads to factor of \~5 speed gains over non-accelerated methods in realistic problems. We have implemented our method in the new code Viscous Accretion Disk Evolution Resource (VADER), which is freely available for download from <a href="https://bitbucket.org/krumholz/vader/">this https URL</a> under the terms of the GNU General Public License.

Tags

Users

  • @ericblackman

Comments and Reviews