On the Existence of Pressure-Equilibrium-Preserving Numerical Fluxes for Supercritical Fluids
Abstract
In this work we propose a new existence theorem for numerical-flux functions for supercritical fluids that are pressure-equilibrium-preserving (PEP). In particular, we characterize the existence of consistent numerical-flux functions that satisfy an algebraic PEP property. Our theory links the existence of PEP schemes to geometric properties of the equation of state describing the thermodynamics of the fluid. When these geometric properties fail for a pair of states on the same isobar, no PEP schemes of the considered form can exist on a domain containing those states. In our analysis the equation of state itself can be fully general only needing to satisfy some fundamental thermodynamic principles. Recently, PEP compatibility conditions for general equations of state have been derived [1] that rely on the existence of certain thermodynamic derivatives which are not guaranteed to be defined under fundamental thermodynamic principles. The geometric conditions in our theory do not depend on the existence of these derivatives and recover the recent compatibility conditions in the case that these derivatives are defined. Finally, using numerical experiments we demonstrate for two supercritical fluids that our existence conditions are restrictive and thus that no PEP schemes of the form we consider exist on the domain we specify for these fluids. Using our geometric perspective, we also shed light on mechanisms by which PEP schemes can develop numerical issues, which we also demonstrate using numerical experiments.
Related
- Analysis, thermodynamics, and a numerical solver for a pressure-temperature equilibrium closure of the four-equation model
- Structure-Preserving Numerical Schemes for a Time-Fractional Porous Medium Equation with Nonlocal Pressure
- Uniqueness and structure-preserving numerical simulations for an Ultrahyperbolic Schrödinger equation
- Stability of ideal Nash equilibrium points in generalized multi-criteria games, interval-valued games, and existence in single-leader--multi-follower problems
- A unified framework for original energy dissipation of equilibrium-preserving exponential Runge--Kutta methods for L2 gradient flows
BibTeX
@misc{klein2026existence,
title = {{On the Existence of Pressure-Equilibrium-Preserving Numerical Fluxes for Supercritical Fluids}},
author = {Robin Ben Klein},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2609.17492},
eprint = {2609.17492},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.17492},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - On the Existence of Pressure-Equilibrium-Preserving Numerical Fluxes for Supercritical Fluids
AU - Robin Ben Klein
PY - 2026
DA - 2026/10/06
PB - arXiv
DO - 10.48550/arxiv.2609.17492
UR - https://arxiv.org/abs/2609.17492
AB - In this work we propose a new existence theorem for numerical-flux functions for supercritical fluids that are pressure-equilibrium-preserving (PEP). In particular, we characterize the existence of consistent numerical-flux functions that satisfy an algebraic PEP property. Our theory links the existence of PEP schemes to geometric properties of the equation of state describing the thermodynamics of the fluid. When these geometric properties fail for a pair of states on the same isobar, no PEP schemes of the considered form can exist on a domain containing those states. In our analysis the equation of state itself can be fully general only needing to satisfy some fundamental thermodynamic principles. Recently, PEP compatibility conditions for general equations of state have been derived [1] that rely on the existence of certain thermodynamic derivatives which are not guaranteed to be defined under fundamental thermodynamic principles. The geometric conditions in our theory do not depend on the existence of these derivatives and recover the recent compatibility conditions in the case that these derivatives are defined. Finally, using numerical experiments we demonstrate for two supercritical fluids that our existence conditions are restrictive and thus that no PEP schemes of the form we consider exist on the domain we specify for these fluids. Using our geometric perspective, we also shed light on mechanisms by which PEP schemes can develop numerical issues, which we also demonstrate using numerical experiments.
ER -