A Subluminal Relativistic Magnetohydrodynamics Scheme with ADER-WENO Predictor and Multidimensional Riemann Solver-Based Corrector

Abstract

The relativistic magnetohydrodynamics (RMHD) set of equations has recently seen increased use in astrophysical computations. Even so, RMHD codes remain fragile. The reconstruction can sometimes yield superluminal velocities in certain parts of the mesh. In this paper we present a reconstruction strategy that overcomes this problem by making a single conservative to primitive transformation per cell followed by higher order WENO reconstruction on a carefully chosen set of primitives that guarantee subluminal reconstruction of the flow variables. For temporal evolution via a predictor step we also present second, third and fourth order accurate ADER methods that keep the velocity subluminal during the predictor step. The RMHD system also requires the magnetic field to be evolved in a divergence-free fashion. In the treatment of classical numerical MHD the analogous issue has seen much recent progress with the advent of multidimensional Riemann solvers. By developing multidimensional Riemann solvers for RMHD, we show that similar advances extend to RMHD. As a result, the face-centered magnetic fields can be evolved much more accurately using the edge-centered electric fields in the corrector step. Those edge-centered electric fields come from a multidimensional Riemann solver for RMHD which we present in this paper. In this paper we also develop several new test problems for RMHD. We show that RMHD vortices can be designed that propagate on the computational mesh as self-preserving structures. These RMHD vortex test problems provide a means to do truly multidimensional accuracy testing for RMHD codes. Several other stringent test problems are presented. We show the importance of resolution in certain test problems. Our tests include a demonstration that RMHD vortices are stable when they interact with shocks.

67 pages, 9 figures, accepted for publication in Journal of Computational Physics

References in corpus (9)

• 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
• Efficient, High Accuracy ADER-WENO Schemes for Hydrodynamics and Divergence-Free Magnetohydrodynamics
• WhiskyMHD: a new numerical code for general relativistic magnetohydrodynamics
• An HLLC Solver for Relativistic Flows -- II. Magnetohydrodynamics
• IllinoisGRMHD: An Open-Source, User-Friendly GRMHD Code for Dynamical Spacetimes
• WHAM: A WENO-based general relativistic numerical scheme I: Hydrodynamics
• Equation of State in Numerical Relativistic Hydrodynamics
• A High-Order Relativistic Two-Fluid Electrodynamic Scheme with Consistent Reconstruction of Electromagnetic Fields and a Multidimensional Riemann Solver for Electromagnetism

Cited by in corpus (22)

• Discontinuous Galerkin algorithms for fully kinetic plasmas
• SpECTRE: A Task-based Discontinuous Galerkin Code for Relativistic Astrophysics
• ADER discontinuous Galerkin schemes for general-relativistic ideal magnetohydrodynamics
• A High-Order Relativistic Two-Fluid Electrodynamic Scheme with Consistent Reconstruction of Electromagnetic Fields and a Multidimensional Riemann Solver for Electromagnetism
• High order ADER schemes and GLM curl cleaning for a first order hyperbolic formulation of compressible flow with surface tension
• Admissible state and physical constraints preserving schemes for relativistic magnetohydrodynamic equations
• Space-time adaptive ADER discontinuous Galerkin schemes for nonlinear hyperelasticity with material failure
• Runge-Kutta discontinuous Galerkin methods for the special relativistic magnetohydrodynamics
• High-order accurate entropy stable nodal discontinuous Galerkin schemes for the ideal special relativistic magnetohydrodynamics
• High-order accurate entropy stable adaptive moving mesh finite difference schemes for special relativistic (magneto)hydrodynamics
• A Physical-Constraint-Preserving Finite Volume WENO Method for Special Relativistic Hydrodynamics on Unstructured Meshes
• On physical-constraints-preserving schemes for special relativistic magnetohydrodynamics with a general equation of state
• Efficient Alternative Finite Difference WENO Schemes for Hyperbolic Conservation Laws
• Implicit gradients based conservative numerical scheme for compressible flows
• A comparison of approximate non-linear Riemann solvers for Relativistic MHD
• High-order finite-difference entropy stable schemes for two-fluid relativistic plasma flow equations
• Physical Constraint Preserving Higher Order Finite Volume Schemes for Divergence-Free Astrophysical MHD and RMHD
• Provably Physical-Constraint-Preserving Discontinuous Galerkin Methods for Multidimensional Relativistic MHD Equations
• An analytical solution of the isentropic vortex problem in the special relativistic magnetohydrodynamics
• An Alternative Finite Difference WENO-like Scheme with Physical Constraint Preservation for Divergence-Preserving Hyperbolic Systems
• Implicit multiderivative collocation solvers for linear partial differential equations with discontinuous Galerkin spatial discretizations
• Minimum Principle on Specific Entropy and High-Order Accurate Invariant Region Preserving Numerical Methods for Relativistic Hydrodynamics