A $4^{\mathrm{th}}$-order accurate finite volume method for ideal classical and special relativistic MHD based on pointwise reconstructions

Abstract

We present a novel implementation of a genuinely $4^{\mathrm{th}}$-order accurate finite volume scheme for multidimensional classical and special relativistic magnetohydrodynamics (MHD) based on the constrained transport (CT) formalism. The scheme introduces several novel aspects when compared to its predecessors yielding a more efficient computational tool. Among the most relevant ones, our scheme exploits pointwise to pointwise reconstructions (rather than one-dimensional finite volume ones), employs the generic upwind constrained transport averaging and sophisticated limiting strategies that include both a discontinuity detector and an order reduction procedure. Selected numerical benchmarks demonstrate the accuracy and robustness of the method.

30 pages, 14 figures

References in corpus (8)

• PLUTO: a Numerical Code for Computational Astrophysics
• A High Order Godunov Scheme with Constrained Transport and Adaptive Mesh Refinement for Astrophysical MHD
• An Unsplit Godunov Method for Ideal MHD via Constrained Transport in Three Dimensions
• ECHO: an Eulerian Conservative High Order scheme for general relativistic magnetohydrodynamics and magnetodynamics
• An HLLC Solver for Relativistic Flows -- II. Magnetohydrodynamics
• A Simple and Accurate Riemann Solver for Isothermal MHD
• Systematic construction of upwind constrained transport schemes for MHD
• A comparison of approximate non-linear Riemann solvers for Relativistic MHD

Cited by in corpus (5)

• Relativistic reconnection with effective resistivity: I. Dynamics and reconnection rate
• EuroHPC SPACE CoE: Redesigning Scalable Parallel Astrophysical Codes for Exascale
• PyPLUTO: a data analysis Python package for the PLUTO code
• High-order finite-volume integration schemes for subsonic magnetohydrodynamics
• The Yin-Yang Magnetic Flux Eruption (Yin-Yang-MFE) Code: A Global Corona Magnetohydrodynamic Code with the Yin-Yang grid