A Physical-Constraint-Preserving Finite Volume WENO Method for Special Relativistic Hydrodynamics on Unstructured Meshes
arXiv:2207.09385 · doi:10.1016/j.jcp.2022.111398
Abstract
This paper presents a highly robust third-order accurate finite volume weighted essentially non-oscillatory (WENO) method for special relativistic hydrodynamics on unstructured triangular meshes. We rigorously prove that the proposed method is physical-constraint-preserving (PCP), namely, always preserves the positivity of the pressure and the rest-mass density as well as the subluminal constraint on the fluid velocity. The method is built on a highly efficient compact WENO reconstruction on unstructured meshes, a simple PCP limiter, the provably PCP property of the Harten--Lax--van Leer flux, and third-order strong-stability-preserving time discretization. Due to the relativistic effects, the primitive variables (namely, the rest-mass density, velocity, and pressure) are highly nonlinear implicit functions in terms of the conservative variables, making the design and analysis of our method nontrivial. To address the difficulties arising from the strong nonlinearity, we adopt a novel quasilinear technique for the theoretical proof of the PCP property. Three provable convergence-guaranteed iterative algorithms are also introduced for the robust recovery of primitive quantities from admissible conservative variables. We also propose a slight modification to an existing WENO reconstruction to ensure the scaling invariance of the nonlinear weights and thus to accommodate the homogeneity of the evolution operator, leading to the advantages of the modified WENO reconstruction in resolving multi-scale wave structures. Extensive numerical examples are presented to demonstrate the robustness, expected accuracy, and high resolution of the proposed method.
arXiv comment56 pages, 18 figures
References in corpus (26)
- Positivity-preserving flux limiters for high-order conservative schemes
- Very High Order $\PNM$ Schemes on Unstructured Meshes for the Resistive Relativistic MHD Equations
- An HLLC Riemann Solver for Relativistic Flows: I. Hydrodynamics
- THC: a new high-order finite-difference high-resolution shock-capturing code for special-relativistic hydrodynamics
- An efficient shock-capturing central-type scheme for multidimensional relativistic flows. I. Hydrodynamics
- High-Order Fully General-Relativistic Hydrodynamics: new Approaches and Tests
- RAM: A Relativistic Adaptive Mesh Refinement Hydrodynamics Code
- WHAM: A WENO-based general relativistic numerical scheme I: Hydrodynamics
- TESS: A Relativistic Hydrodynamics Code on a Moving Voronoi Mesh
- SpECTRE: A Task-based Discontinuous Galerkin Code for Relativistic Astrophysics
- High-order accurate physical-constraints-preserving finite difference WENO schemes for special relativistic hydrodynamics
- Discontinuous Galerkin methods for general-relativistic hydrodynamics: formulation and application to spherically symmetric spacetimes
- Provably Positive High-Order Schemes for Ideal Magnetohydrodynamics: Analysis on General Meshes
- Assessment of a high-resolution central scheme for the solution of the relativistic hydrodynamics equations
- Positivity-Preserving Analysis of Numerical Schemes for Ideal Magnetohydrodynamics
- Provably Positive Discontinuous Galerkin Methods for Multidimensional Ideal Magnetohydrodynamics
- Design of Provably Physical-Constraint-Preserving Methods for General Relativistic Hydrodynamics
- Admissible state and physical constraints preserving schemes for relativistic magnetohydrodynamic equations
- A Subluminal Relativistic Magnetohydrodynamics Scheme with ADER-WENO Predictor and Multidimensional Riemann Solver-Based Corrector
- Physical-constraints-preserving central discontinuous Galerkin methods for special relativistic hydrodynamics with a general equation of state
- Formulation of discontinuous Galerkin methods for relativistic astrophysics
- Numerical relativity in spherical coordinates: A new dynamical spacetime and general relativistic MHD evolution framework for the Einstein Toolkit
- Entropy Symmetrization and High-Order Accurate Entropy Stable Numerical Schemes for Relativistic MHD Equations
- High-order accurate entropy stable finite difference schemes for one- and two-dimensional special relativistic hydrodynamics
- Physical-constraints-preserving Lagrangian finite volume schemes for one- and two-dimensional special relativistic hydrodynamics
- Second-order accurate BGK schemes for the special relativistic hydrodynamics with the Synge equation of state