A staggered grid, Lagrangian-Eulerian remap code for 3-D MHD simulations

T D Arber, A W Longbottom, C L Gerrard, A M Milne

Research output: Contribution to journalArticlepeer-review

254 Citations (Scopus)


In this paper an approach to multidimensional magnetohydrodynamics (MI-ID) which correctly handles shocks but does not use an approximate Riemann solver is proposed. This approach is simple and is based on control volume averaging with a staggered grid. The method builds on the older and often overlooked technique of on each step taking a fully 3-D Lagrangian step and then conservatively remapping onto the original grid. At the remap step gradient limiters are applied so that the scheme is monotonicity-preserving. For Euler's equations this technique, combined with an appropriately staggered grid and Wilkins artificial viscosity, can give results comparable to those from approximate Riemann solvers. We show how this can be extended to include a magnetic field, maintaining the divergence-free condition and pressure positivity and then present numerical test results. Where possible a comparison with other shock capturing techniques is presented and the advantages and disadvantages of the proposed scheme are clearly explained. (C) 2001 Academic Press.

Original languageEnglish
Pages (from-to)151-181
Number of pages31
JournalJournal of Computational Physics
Publication statusPublished - 20 Jul 2001


  • MHD
  • shocks
  • Lagrangian
  • FLOW
  • FLUX


Dive into the research topics of 'A staggered grid, Lagrangian-Eulerian remap code for 3-D MHD simulations'. Together they form a unique fingerprint.

Cite this