GLOBAL SIMULATIONS OF ACCRETION DISKS. I. CONVERGENCE AND COMPARISONS WITH LOCAL MODELS
- Department of Mathematics, University of Maryland, College Park, MD 20742-2421 (United States)
- Department of Astronomy and the Maryland Astronomy Center for Theory and Computation, University of Maryland, College Park, MD 20742-2421 (United States)
- Department of Astrophysical Sciences, Princeton University, Princeton, NJ 08544 (United States)
- JILA, University of Colorado at Boulder, Boulder, CO 80309-0440 (United States)
Grid-based magnetohydrodynamic (MHD) simulations have proven invaluable for the study of astrophysical accretion disks. However, the fact that angular momentum transport in disks is mediated by MHD turbulence (with structure down to very small scales) raises the concern that the properties of the modeled accretion disks are affected by the finite numerical resolution of the simulation. By implementing an orbital advection algorithm into the Athena code in cylindrical geometry, we have performed a set of global (but unstratified) Newtonian disk simulations extending up to resolutions previously unattained. We study the convergence of these models as a function of spatial resolution and initial magnetic field geometry. The usual viscosity parameter ({alpha}) or the ratio of thermal-to-magnetic pressure ({beta}) is found to be a poor diagnostic of convergence, whereas the average tilt angle of the magnetic field in the (r, {phi})-plane is a very good diagnostic of convergence. We suggest that this is related to the saturation of the MHD turbulence via parasitic modes of the magnetorotational instability. Even in the case of zero-net magnetic flux, we conclude that our highest resolution simulations (with 32 zones and 64 zones per vertical scale height) have achieved convergence. Our global simulations reach resolutions comparable to those used in local, shearing-box models of MHD disk turbulence. We find that the saturation predictors derived from local simulations correspond well to the instantaneous correlations between local flux and stress found in our global simulations. However, the conservation of magnetic flux implicit in local models is not realized in our global disks. Thus, the magnetic connectivity of an accretion disk represents physics that is truly global and cannot be captured in any ab initio local model.
- OSTI ID:
- 22020371
- Journal Information:
- Astrophysical Journal, Vol. 749, Issue 2; Other Information: Country of input: International Atomic Energy Agency (IAEA); ISSN 0004-637X
- Country of Publication:
- United States
- Language:
- English
Similar Records
TURBULENCE IN THE OUTER REGIONS OF PROTOPLANETARY DISKS. I. WEAK ACCRETION WITH NO VERTICAL MAGNETIC FLUX
Global MHD Simulations of Accretion Disks in Cataclysmic Variables (CVs). II. The Relative Importance of MRI and Spiral Shocks