MHD simulations of the magnetorotational instability in a shearing box with zero net flux

Free access article

Issue		A&A Volume 476, Number 3, December IV 2007


Page(s)		1113 - 1122
Section		Astrophysical processes
DOI		http://dx.doi.org/10.1051/0004-6361:20077942

A&A 476, 1113-1122 (2007)
DOI: 10.1051/0004-6361:20077942

I. The issue of convergence

S. Fromang and J. Papaloizou

Department of Applied Mathematics and Theoretical Physics, University of Cambridge, Centre for Mathematical Sciences, Wilberforce Road, Cambridge, CB3 0WA, UK
e-mail: S.Fromang@damtp.cam.ac.uk

(Received 24 May 2007 / Accepted 30 September 2007)

Abstract
Aims.We study the properties of MHD turbulence driven by the magnetorotational instability (MRI) in accretion disks. To do this we perform a series of numerical simulations for which the resolution is gradually increased.
Methods.We adopt the local shearing box model and focus on the special case for which the initial magnetic flux threading the disk vanishes. We employ the finite difference code ZEUS to evolve the ideal MHD equations.
Results.Performing a set of numerical simulations in a fixed computational domain with increasing resolution, we demonstrate that turbulent activity decreases as resolution increases. The highest resolution considered is 256 grid cells per scale height. We quantify the turbulent activity by measuring the rate of angular momentum transport through evaluating the standard $\alpha$ parameter. We find $\alpha$ = 0.004 when (N_x,N_y,N_z) = (64,100,64), $\alpha$ = 0.002 when (N_x,N_y,N_z) = (128,200,128) and $\alpha$ = 0.001 when (N_x,N_y,N_z) = (256,400,256). This steady decline is an indication that numerical dissipation, occurring at the grid scale is an important determinant of the saturated form of the MHD turbulence. Analysing the results in Fourier space, we demonstrate that this is due to the MRI forcing significant flow energy all the way down to the grid dissipation scale. We also use our results to study the properties of the numerical dissipation in ZEUS. Its amplitude is characterised by the magnitude of an effective magnetic Reynolds number $\rm Re_{M}$ which increases from 10⁴ to 10⁵ as the number of grid points is increased from 64 to 256 per scale height.
Conclusions.The simulations we have carried out do not produce results that are independent of the numerical dissipation scale, even at the highest resolution studied. Thus it is important to use physical dissipation, both viscous and resistive, and to quantify contributions from numerical effects, when performing numerical simulations of MHD turbulence with zero net flux in accretion disks at the resolutions normally considered.

Key words: accretion, accretion disks -- magnetohydrodynamics (MHD) -- methods: numerical

What is OpenURL?

The OpenURL standard is a protocol for transmission of metadata describing the resource that you wish to access. An OpenURL link contains article metadata and directs it to the OpenURL server of your choice. The OpenURL server can provide access to the resource and also offer complementary services (specific search engine, export of references...). The OpenURL link can be generated by different means.

If your librarian has set up your subscription with an OpenURL resolver, OpenURL links appear automatically on the abstract pages.
You can define your own OpenURL resolver with your EDPS Account. In this case your choice will be given priority over that of your library.
You can use an add-on for your browser (Firefox or I.E.) to display OpenURL links on a page (see http://www.openly.com/openurlref/). You should disable this module if you wish to use the OpenURL server that you or your library have defined.