Testing Convergence for Global Accretion Disks
Abstract
Global disk simulations provide a powerful tool for investigating accretion and the underlying magnetohydrodynamic turbulence driven by magneto-rotational instability (MRI). Using them to accurately predict quantities such as stress, accretion rate, and surface brightness profile requires that purely numerical effects, arising from both resolution and algorithm, be understood and controlled. We use the flux-conservative Athena code to conduct a series of experiments on disks having a variety of magnetic topologies to determine what constitutes adequate resolution. We develop and apply several resolution metrics: langQz rang and langQ phirang, the ratio of the grid zone size to the characteristic MRI wavelength, αmag, the ratio of the Maxwell stress to the magnetic pressure, and \langle B_R^2\rangle /\langle B_\phi ^2\rangle, the ratio of radial to toroidal magnetic field energy. For the initial conditions considered here, adequate resolution is characterized by langQz rang >= 15, langQ phirang >= 20, αmag ≈ 0.45, and \langle B_R^2\rangle /\langle B_\phi ^2\rangle \approx 0.2. These values are associated with >=35 zones per scaleheight H, a result consistent with shearing box simulations. Numerical algorithm is also important. Use of the Harten-Lax-van Leer-Einfeldt flux solver or second-order interpolation can significantly degrade the effective resolution compared to the Harten-Lax-van Leer discontinuities flux solver and third-order interpolation. Resolution at this standard can be achieved only with large numbers of grid zones, arranged in a fashion that matches the symmetries of the problem and the scientific goals of the simulation. Without it, however, quantitative measures important to predictions of observables are subject to large systematic errors.
- Publication:
-
The Astrophysical Journal
- Pub Date:
- August 2013
- DOI:
- 10.1088/0004-637X/772/2/102
- arXiv:
- arXiv:1306.0243
- Bibcode:
- 2013ApJ...772..102H
- Keywords:
-
- accretion;
- accretion disks;
- black hole physics;
- magnetohydrodynamics: MHD;
- methods: numerical;
- Astrophysics - Instrumentation and Methods for Astrophysics;
- Astrophysics - High Energy Astrophysical Phenomena
- E-Print:
- To appear in the Astrophysical Journal