On Local Minimum Entropy Principle of High-Order Schemes for Relativistic Euler Equations

This paper establishes the minimum entropy principle (MEP) for the relativistic Euler equations with a broad class of equations of state (EOSs) and addresses the challenge of preserving the local version of the discovered MEP in high-order numerical schemes. At the continuous level, we find out a family of entropy pairs for the relativistic Euler equations and provide rigorous analysis to prove the strict convexity of entropy under a necessary and sufficient condition. At the numerical level, we develop a rigorous framework for designing provably entropy-preserving high-order schemes that ensure both physical admissibility and the discovered MEP. The relativistic effects, coupled with the abstract and general EOS formulation, introduce significant challenges not encountered in the nonrelativistic case or with the ideal EOS. In particular, entropy is a highly nonlinear and implicit function of the conservative variables, making it particularly difficult to enforce entropy preservation. To address these challenges, we establish a series of auxiliary theories via highly technical inequalities. Another key innovation is the use of geometric quasi-linearization (GQL), which reformulates the nonlinear constraints into equivalent linear ones by introducing additional free parameters. These advancements form the foundation of our entropy-preserving analysis. We propose novel, robust, locally entropy-preserving high-order frameworks. A central challenge is accurately estimating the local minimum of entropy, particularly in the presence of shock waves at unknown locations. To address this, we introduce two new approaches for estimating local lower bounds of specific entropy, which prove effective for both smooth and discontinuous problems. Numerical experiments demonstrate that our entropy-preserving methods maintain high-order accuracy while effectively suppressing spurious oscillations.