Provably Physical-Constraint-Preserving Discontinuous Galerkin Methods for Multidimensional Relativistic MHD Equations

02/09/2020
by   Kailiang Wu, et al.
0

We propose and analyze a class of robust, uniformly high-order accurate discontinuous Galerkin (DG) schemes for multidimensional relativistic magnetohydrodynamics (RMHD) on general meshes. A distinct feature of the schemes is their physical-constraint-preserving (PCP) property, i.e., they are proven to preserve the subluminal constraint on the fluid velocity and the positivity of density, pressure, and specific internal energy. Developing PCP high-order schemes for RMHD is highly desirable but remains a challenging task, especially in the multidimensional cases, due to the inherent strong nonlinearity in the constraints and the effect of the magnetic divergence-free condition. Inspired by some crucial observations at the PDE level, we construct the provably PCP schemes by using the locally divergence-free DG schemes of the recently proposed symmetrizable RMHD equations as the base schemes, a limiting technique to enforce the PCP property of the DG solutions, and the strong-stability-preserving methods for time discretization. We rigorously prove the PCP property by using a novel "quasi-linearization" approach to handle the highly nonlinear physical constraints, technical splitting to offset the influence of divergence error, and sophisticated estimates to analyze the beneficial effect of the additional source term in the symmetrizable RMHD system. Several two-dimensional numerical examples are provided to confirm the PCP property and to demonstrate the accuracy, effectiveness and robustness of the proposed PCP schemes.

READ FULL TEXT

page 20

page 22

research
03/28/2022

Provably Positive Central DG Schemes via Geometric Quasilinearization for Ideal MHD Equations

In the numerical simulation of ideal MHD, keeping the pressure and densi...
research
02/07/2021

Minimum Principle on Specific Entropy and High-Order Accurate Invariant Region Preserving Numerical Methods for Relativistic Hydrodynamics

This paper explores Tadmor's minimum entropy principle for the relativis...
research
06/17/2022

Locally Structure-Preserving div-curl operators for high order Discontinuous Galerkin schemes

We propose a novel Structure-Preserving Discontinuous Galerkin (SPDG) op...
research
05/24/2023

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

This paper proposes and analyzes a novel efficient high-order finite vol...
research
07/19/2022

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

This paper designs and analyzes positivity-preserving well-balanced (WB)...
research
07/19/2022

A Physical-Constraint-Preserving Finite Volume WENO Method for Special Relativistic Hydrodynamics on Unstructured Meshes

This paper presents a highly robust third-order accurate finite volume w...

Please sign up or login with your details

Forgot password? Click here to reset