A fourth-order accurate finite volume method for ideal MHD via upwind constrained transport

Kyle Gerard Felker, James McLellan Stone

Research output: Contribution to journalArticle

10 Scopus citations

Abstract

We present a fourth-order accurate finite volume method for the solution of ideal magnetohydrodynamics (MHD). The numerical method combines high-order quadrature rules in the solution of semi-discrete formulations of hyperbolic conservation laws with the upwind constrained transport (UCT) framework to ensure that the divergence-free constraint of the magnetic field is satisfied. A novel implementation of UCT that uses the piecewise parabolic method (PPM) for the reconstruction of magnetic fields at cell corners in 2D is introduced. The resulting scheme can be expressed as the extension of the second-order accurate constrained transport (CT) Godunov-type scheme that is currently used in the Athena astrophysics code. After validating the base algorithm on a series of hydrodynamics test problems, we present the results of multidimensional MHD test problems which demonstrate formal fourth-order convergence for smooth problems, robustness for discontinuous problems, and improved accuracy relative to the second-order scheme.

Original languageEnglish (US)
Pages (from-to)1365-1400
Number of pages36
JournalJournal of Computational Physics
Volume375
DOIs
StatePublished - Dec 15 2018

All Science Journal Classification (ASJC) codes

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

Keywords

  • Constrained transport
  • High-order finite volume method
  • Magnetohydrodynamics
  • Numerical methods

Fingerprint Dive into the research topics of 'A fourth-order accurate finite volume method for ideal MHD via upwind constrained transport'. Together they form a unique fingerprint.

  • Cite this