Provably convergent Newton-Raphson methods for recovering primitive variables with applications to physical-constraint-preserving Hermite WENO schemes for relativistic hydrodyn

05/24/2023
by   Chaoyi Cai, et al.
0

The relativistic hydrodynamics (RHD) equations have three crucial intrinsic physical constraints on the primitive variables: positivity of pressure and density, and subluminal fluid velocity. However, numerical simulations can violate these constraints, leading to nonphysical results or even simulation failure. Designing genuinely physical-constraint-preserving (PCP) schemes is very difficult, as the primitive variables cannot be explicitly reformulated using conservative variables due to relativistic effects. In this paper, we propose three efficient Newton–Raphson (NR) methods for robustly recovering primitive variables from conservative variables. Importantly, we rigorously prove that these NR methods are always convergent and PCP, meaning they preserve the physical constraints throughout the NR iterations. The discovery of these robust NR methods and their PCP convergence analyses are highly nontrivial and technical. As an application, we apply the proposed NR methods to design PCP finite volume Hermite weighted essentially non-oscillatory (HWENO) schemes for solving the RHD equations. Our PCP HWENO schemes incorporate high-order HWENO reconstruction, a PCP limiter, and strong-stability-preserving time discretization. We rigorously prove the PCP property of the fully discrete schemes using convex decomposition techniques. Moreover, we suggest the characteristic decomposition with rescaled eigenvectors and scale-invariant nonlinear weights to enhance the performance of the HWENO schemes in simulating large-scale RHD problems. Several demanding numerical tests are conducted to demonstrate the robustness, accuracy, and high resolution of the proposed PCP HWENO schemes and to validate the efficiency of our NR methods.

READ FULL TEXT

page 40

page 41

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...
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
02/09/2020

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

We propose and analyze a class of robust, uniformly high-order accurate ...
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
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
11/08/2021

Geometric Quasilinearization Framework for Analysis and Design of Bound-Preserving Schemes

Solutions to many partial differential equations satisfy certain bounds ...
research
01/13/2021

Reformulated dissipation for the free-stream preserving of the conservative finite difference schemes on curvilinear grids

In this paper, we develop a new free-stream preserving (FP) method for h...

Please sign up or login with your details

Forgot password? Click here to reset