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)		1123 - 1132
Section		Astrophysical processes
DOI		http://dx.doi.org/10.1051/0004-6361:20077943

A&A 476, 1123-1132 (2007)
DOI: 10.1051/0004-6361:20077943

II. The effect of transport coefficients

S. Fromang¹, J. Papaloizou¹, G. Lesur², and T. Heinemann¹

¹ Department of Applied Mathematics and Theoretical Physics, University of Cambridge, Centre for Mathematical Sciences, Wilberforce Road, Cambridge, CB3 0WA, UK
² Laboratoire d'astrophysique, UJF CNRS, BP 53, 38041 Grenoble Cedex 9, France
e-mail: S.Fromang@damtp.cam.ac.uk

(Received 24 May 2007 / Accepted 30 September 2007)

Abstract
Aims.We study the influence of the choice of transport coefficients (viscosity and resistivity) on MHD turbulence driven by the magnetorotational instability (MRI) in accretion disks.
Methods.We follow the methodology described in Paper I: we adopt an unstratified shearing box model and focus on the case where the net vertical magnetic flux threading the box vanishes. For the most part we use the operator split code ZEUS, including explicit transport coefficients in the calculations. However, we also compare our results with those obtained using other algorithms (NIRVANA, the PENCIL code and a spectral code) to demonstrate both the convergence of our results and their independence of the numerical scheme.
Results.We find that small scale dissipation affects the saturated state of MHD turbulence. In agreement with recent similar numerical simulations done in the presence of a net vertical magnetic flux, we find that turbulent activity (measured by the rate of angular momentum transport) is an increasing function of the magnetic Prandtl number Pm for all values of the Reynolds number Re that we investigated. We also found that turbulence disappears when the Prandtl number falls below a critical value $\rm Pm_{c}$ that is apparently a decreasing function of Re. For the limited region of parameter space that can be probed with current computational resources, we always obtained $\rm Pm_{c}>1$ .
Conclusions.We conclude that the magnitudes of the transport coefficients are important in determining the properties of MHD turbulence in numerical simulations in the shearing box with zero net flux, at least for Reynolds numbers and magnetic Prandtl numbers that are such that transport is not dominated by numerical effects and thus can be probed using current computational resources.

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.