An unsplit Godunov method for ideal MHD via constrained transport

Thomas A. Gardiner, James McLellan Stone

Research output: Contribution to journalArticlepeer-review

415 Scopus citations

Abstract

We describe a single step, second-order accurate Godunov scheme for ideal MHD based on combining the piecewise parabolic method (PPM) for performing spatial reconstruction, the corner transport upwind (CTU) method of Colella for multidimensional integration, and the constrained transport (CT) algorithm for preserving the divergence-free constraint on the magnetic field. We adopt the most compact form of CT, which requires the field be represented by area-averages at cell faces. We demonstrate that the fluxes of the area-averaged field used by CT can be made consistent with the fluxes of the volume-averaged field returned by a Riemann solver if they obey certain simple relationships. We use these relationships to derive new algorithms for constructing the CT fluxes at grid cell corners which reduce exactly to the equivalent one-dimensional solver for plane-parallel, grid-aligned flow. We show that the PPM reconstruction algorithm must include multidimensional terms for MHD, and we describe a number of important extensions that must be made to CTU in order for it to be used for MHD with CT. We present the results of a variety of test problems to demonstrate the method is accurate and robust.

Original languageEnglish (US)
Pages (from-to)509-539
Number of pages31
JournalJournal of Computational Physics
Volume205
Issue number2
DOIs
StatePublished - May 20 2005

All Science Journal Classification (ASJC) codes

  • Numerical Analysis
  • Modeling and Simulation
  • Physics and Astronomy (miscellaneous)
  • General Physics and Astronomy
  • Computer Science Applications
  • Computational Mathematics
  • Applied Mathematics

Fingerprint

Dive into the research topics of 'An unsplit Godunov method for ideal MHD via constrained transport'. Together they form a unique fingerprint.

Cite this