Entropy conservative discretization of compressible Euler equations with an arbitrary equation of state

Alessandro Aiello,Carlo De Michele,Gennaro Coppola
2024-08-02
Abstract:This study proposes a novel spatial discretization procedure for the compressible Euler equations which guarantees entropy conservation at a discrete level when an arbitrary equation of state is assumed. The proposed method, based on a locally-conservative discretization, guarantees also the spatial conservation of mass, momentum, and total energy and is kinetic energy-preserving. In order to achieve the entropy-conservation property for an arbitrary non-ideal gas, a general strategy is adopted based on the manipulation of discrete balance equations through the imposition of global entropy conservation and the use of a summation by parts rule. The procedure, which is extended to an arbitrary order of accuracy, conducts to a general form of the internal-energy numerical flux which results in a nonlinear function of thermodynamic and dynamic variables and still admits the mass flux as a residual degree of freedom. The effectiveness of the novel entropy-conservative formulation is demonstrated through numerical tests making use of some of the most popular cubic equations of state.
Fluid Dynamics,Numerical Analysis
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is to develop a new spatial discretization method for solving the compressible Euler equations to ensure entropy conservation at the discrete level under any Equation of State (EoS). Specifically, the research objectives include: 1. **Entropy Conservation**: Ensure that in numerical simulations, for smooth, non - viscous flows, entropy can be discretely conserved. 2. **Any Equation of State**: The proposed method is applicable not only to ideal gases but also to any non - ideal gases (such as real gases), thus expanding the application range of existing methods. 3. **Structure - Preserving Properties**: In addition to entropy conservation, this method also guarantees the spatial conservation of mass, momentum and total energy, and is Kinetic Energy Preserving (KEP). ### Research Background and Motivation Traditional numerical discretization methods are prone to strong nonlinear instabilities when dealing with compressible flows at high Reynolds numbers, especially when using non - dissipative schemes to discretize the convection terms. These instabilities usually stem from the failure of numerical methods to reproduce some structural properties or invariants of the continuous system at the discrete level. For example, the compressible Euler equations have many induced symmetries, and these symmetries are usually not accurately reproduced by the discrete system. To improve the fidelity of numerical simulations and suppress instabilities, in recent years, researchers have begun to focus on Structure - Preserving or Physically Compatible Methods. Such methods aim to discretely reproduce the structural properties of the continuous system on any finite grid. ### Main Contributions This paper proposes a new discretization method, and its main features are as follows: - **Entropy Conservation**: By enforcing global entropy conservation, the numerical fluxes that satisfy entropy conservation are derived. - **Any Equation of State**: The method is applicable to any Equation of State, especially the Equation of State of real gases, which is of great significance in engineering applications. - **High - Order Accuracy**: The method can be extended to any order of accuracy, and a detailed derivation is carried out at the second - order accuracy. - **Local Conservativeness**: It ensures the spatial conservation of mass, momentum and total energy, and is kinetic energy - preserving. ### Mathematical Derivation and Verification The author, through a series of mathematical derivations, shows how to design numerical fluxes through the global entropy conservation constraint and proves that these fluxes can maintain entropy conservation under any Equation of State. In addition, the article also verifies the effectiveness of the proposed method through numerical tests, using some popular cubic equations of state for testing. In conclusion, this research provides a new and efficient numerical method for dealing with compressible flows of real gases, which can maintain important physical properties at the discrete level, thereby improving the accuracy and stability of numerical simulations.