Abstract:This paper develops novel and robust central discontinuous Galerkin (CDG) schemes of arbitrarily high-order accuracy for special relativistic magnetohydrodynamics (RMHD) with a general equation of state (EOS). These schemes are provably bound-preserving (BP), meaning they consistently preserve the upper bound for subluminal fluid velocity and the positivity of density and pressure, while also (locally) maintaining the divergence-free (DF) constraint for the magnetic field. For 1D RMHD, the standard CDG method is exactly DF, and its BP property is proven under a condition achievable by the BP limiter. For 2D RMHD, we design provably BP and locally DF CDG schemes based on the suitable discretization of a modified RMHD system, which is the relativistic analogue of Godunov's symmetrizable form of the non-relativistic MHD system [S. K. Godunov, Numerical Methods for Mechanics of Continuum Medium , 1 (1972) 26–34]. A key novelty in our schemes is the meticulous discretization of additional source terms in the modified RMHD equations, so as to precisely counteract the influence of divergence errors on the BP property across overlapping meshes. Notably, we provide rigorous proofs of the BP property for our CDG schemes and first establish the theoretical connection between BP and discrete DF properties on overlapping meshes for RMHD. Owing to the absence of explicit expressions for primitive variables in terms of conserved variables, the constraints of physical bounds are strongly nonlinear, making the BP proofs highly nontrivial. We overcome these challenges through technical estimates within the geometric quasilinearization (GQL) framework [K. Wu & C.-W. Shu, SIAM Review , 65 (2023) 1031–1073], which equivalently converts the nonlinear constraints into linear ones. Furthermore, we introduce a new 2D cell average decomposition on overlapping meshes, which relaxes the theoretical BP CFL constraint and reduces the number of internal nodes, thereby enhancing the efficiency of the 2D BP CDG method. Finally, we implement the proposed CDG schemes for extensive RMHD problems with various EOSs, demonstrating their robustness and effectiveness in challenging scenarios like ultra-relativistic blasts and jets in strongly magnetized environments.

Positivity-preserving DG and central DG methods for ideal MHD equations

Structure-Preserving Oscillation-Eliminating Discontinuous Galerkin Schemes for Ideal MHD Equations: Locally Divergence-Free and Positivity-Preserving

Positivity-preserving third order DG schemes for Poisson–Nernst–Planck equations

High-order Accurate Entropy Stable Nodal Discontinuous Galerkin Schemes for the Ideal Special Relativistic Magnetohydrodynamics.

A new discretely divergence-free positivity-preserving high-order finite volume method for ideal MHD equations

An Entropy-Stable Discontinuous Galerkin Discretization of the Ideal Multi-Ion Magnetohydrodynamics System

A Divergence-Free and $H(div)$-Conforming Embedded-Hybridized DG Method for the Incompressible Resistive MHD equations

PHYSICAL-CONSTRAINT-PRESERVING CENTRAL DISCONTINUOUS GALERKIN METHODS FOR SPECIAL RELATIVISTIC HYDRODYNAMICS WITH A GENERAL EQUATION OF STATE

Provably Positivity-Preserving Constrained Transport (PPCT) Second-Order Scheme for Ideal Magnetohydrodynamics

Uniformly High-Order Structure-Preserving Discontinuous Galerkin Methods for Euler Equations with Gravitation: Positivity and Well-Balancedness

Positivity-Preserving Well-Balanced Central Discontinuous Galerkin Schemes for the Euler Equations under Gravitational Fields

Well-balanced Positivity-Preserving High-Order Discontinuous Galerkin Methods for Euler Equations with Gravitation

An entropy stable nodal discontinuous Galerkin method for the resistive MHD equations: Continuous analysis and GLM divergence cleaning

Positivity-preserving and energy-dissipating discontinuous Galerkin methods for nonlinear nonlocal Fokker-Planck equations

Entropy-Stable Gauss Collocation Methods for Ideal Magneto-Hydrodynamics

GQL-Based Bound-Preserving and Locally Divergence-Free Central Discontinuous Galerkin Schemes for Relativistic Magnetohydrodynamics

A high-order positivity-preserving single-stage single-step method for the ideal magnetohydrodynamic equations

A Locally Divergence-Free Oscillation-Eliminating Discontinuous Galerkin Method for Ideal Magnetohydrodynamic Equations

Positivity-preserving and entropy-bounded discontinuous Galerkin method for the chemically reacting, compressible Euler equations. Part II: The multidimensional case