Abstract: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.

Conic Programming for Circuit Equations with Rational Current Controlled Resistors.

On the Convergence of the Accelerated Riccati Iteration Method.

SOCP Relaxations of Optimal Power Flow Problem Considering Current Margins in Radial Networks

Stochastic Linear-Quadratic Control with Conic Control Constraints on an Infinite Time Horizon

Conic approximation to nonconvex quadratic programming with convex quadratic constraints

Lanczos-type Variants of the COCR Method for Complex Nonsymmetric Linear Systems

Kkt Solution And Conic Relaxation For Solving Quadratically Constrained Quadratic Programming Problems

A Conic Model for Electrolyzer Scheduling

A Nonmonotone Trust-Region Method of Conic Model for Unconstrained Optimization

A Study on the Strong Duality of Second-Order Conic Relaxation of AC Optimal Power Flow in Radial Networks

Fixed-Order H2-Conic Control

A Constrained Optimization Approach for Lcp

A Study on the Strong Duality of Conic Relaxation of AC Optimal Power Flow in Radial Networks

Parametric Solutions to the Regulator-Conjugate Matrix Equations

A Method for Controlling Parallel and Series RLC Circuits with Time-Varying Resistance, Inductance and Capacitance

Conic Descent Redux for Memory-Efficient Optimization

Quasi-Minimal Residual Variants of the COCG and COCR Methods for Complex Symmetric Linear Systems in Electromagnetic Simulations

Control of Complex Nonlinear Dynamic Rational Systems.

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

Optimal placement of capacitors in a radial network using conic and mixed integer linear programming

Further Development in Convex Conic Reformulation of Geometric Nonconvex Conic Optimization Problems