EDP Sciences
Free access
Volume 457, Number 2, October II 2006
Page(s) 371 - 384
Section Astrophysical processes
DOI http://dx.doi.org/10.1051/0004-6361:20065371

A&A 457, 371-384 (2006)
DOI: 10.1051/0004-6361:20065371

A high order Godunov scheme with constrained transport and adaptive mesh refinement for astrophysical magnetohydrodynamics

S. Fromang1, 2, P. Hennebelle3 and R. Teyssier4

1  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
2  Astronomy Unit, Queen Mary, University of London, Mile End Road, London E1 4NS, UK
3  Laboratoire de radioastronomie millimétrique, UMR 8112 du CNRS, École normale supérieure et Observatoire de Paris, 24 rue Lhomond, 75231 Paris Cedex 05, France
4  Service d'Astrophysique, CEA/DSM/DAPNIA/SAp, Centre d'Études de Saclay, L'orme des Merisiers, 91191 Gif-sur-Yvette Cedex, France

(Received 5 April 2006 / Accepted 30 June 2006 )

Aims. In this paper, we present a new method to perform numerical simulations of astrophysical MHD flows using the Adaptive Mesh Refinement framework and Constrained Transport.
Methods. The algorithm is based on a previous work in which the MUSCL-Hancock scheme was used to evolve the induction equation. In this paper, we detail the extension of this scheme to the full MHD equations and discuss its properties.
Results. Through a series of test problems, we illustrate the performances of this new code using two different MHD Riemann solvers (Lax-Friedrich and Roe) and the need of the Adaptive Mesh Refinement capabilities in some cases. Finally, we show its versatility by applying it to two completely different astrophysical situations well studied in the past years: the growth of the magnetorotational instability in the shearing box and the collapse of magnetized cloud cores.
Conclusions. We have implemented a new Godunov scheme to solve the ideal MHD equations in the AMR code RAMSES. We have shown that it results in a powerful tool that can be applied to a great variety of astrophysical problems, ranging from galaxies formation in the early universe to high resolution studies of molecular cloud collapse in our galaxy.

Key words: magnetohydrodynamics (MHD) -- methods: numerical

© ESO 2006

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.