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

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.

42 pages, 22 figures, accepted for publication in J. Comp. Phys

References in corpus (4)

• Athena: A New Code for Astrophysical MHD
• 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

Cited by in corpus (30)

• The Athena++ Adaptive Mesh Refinement Framework: Design and Magnetohydrodynamic Solvers
• Dedalus: A Flexible Framework for Numerical Simulations with Spectral Methods
• The Arepo public code release
• Beyond second-order convergence in simulations of magnetised binary neutron stars with realistic microphysics
• GRAthena++: puncture evolutions on vertex-centered oct-tree AMR
• Resistive relativistic MHD simulations of astrophysical jets
• The Decoupling of Binaries from Their Circumbinary Disks
• Systematic construction of upwind constrained transport schemes for MHD
• Dissipative Magnetohydrodynamics for Non-Resistive Relativistic Plasmas
• Quokka: A code for two-moment AMR radiation hydrodynamics on GPUs
• Local Simulations of Heating Torques on a Luminous Body in an Accretion Disk
• A constrained transport divergence-free finite element method for Incompressible MHD equations
• Improved Coupling of Hydrodynamics and Nuclear Reactions via Spectral Deferred Corrections
• Addition of tabulated equation of state and neutrino leakage support to IllinoisGRMHD
• A $4^{\mathrm{th}}$-order accurate finite volume method for ideal classical and special relativistic MHD based on pointwise reconstructions
• A high-order weighted finite difference scheme with a multi-state approximate Riemann solver for divergence-free magnetohydrodynamic simulations
• Long-term GRMHD simulation of magnetic field in isolated neutron stars
• A Multistate Low-dissipation Advection Upstream Splitting Method for Ideal Magnetohydrodynamics
• WENO-Wombat: Scalable Fifth-Order Constrained-Transport Magnetohydrodynamics for Astrophysical Applications
• RAM: Rapid Advection Algorithm on Arbitrary Meshes
• High-order finite-volume integration schemes for subsonic magnetohydrodynamics
• Formation of Orion Fingers
• On collective nature of nonlinear torsional Alfvén waves
• Digital Infrastructure in Astrophysics
• General relativistic hydrodynamics code for dynamical spacetimes with curvilinear coordinates, tabulated equations of state, and neutrino physics
• A Co-Scaling Grid for Athena++
• Towards Higher Order Accuracy in Self-Gravitating Hydrodynamics
• Extremely High-Order Convergence in Simulations of Relativistic Stars
• Numerical Methods for Simulating Star Formation
• Novel approach to general curvilinear coordinates for plasma fluid applications