Asymptotic preserving finite volume method for the compressible Euler equations: analysis via dissipative measure-valued solutions

K.R. Arun,Amogh Krishnamurthy,Mária Lukáčová-Medvid'ová
2024-07-18
Abstract:We propose and analyze a new asymptotic preserving (AP) finite volume scheme for the multidimensional compressible barotropic Euler equations to simulate low Mach number flows. The proposed scheme uses a stabilized upwind numerical flux, with the stabilization term being proportional to the stiff pressure gradient, and we prove its conditional energy stability and consistency. Utilizing the concept of dissipative measure-valued (DMV) solutions, we rigorously illustrate the AP properties of the scheme for well-prepared initial data. In particular, we prove that the numerical solutions will converge weakly to a DMV solution of the compressible Euler equations as the mesh parameter vanishes, while the Mach number is fixed. The DMV solutions then converge to a classical solution of the incompressible Euler system as the Mach number goes to zero. Conversely, we show that if the mesh parameter is kept fixed, we obtain an energy stable and consistent finite-volume scheme approximating the incompressible Euler equations as the Mach number goes to zero. The numerical solutions generated by this scheme then converge weakly to a DMV solution of the incompressible Euler system as the mesh parameter vanishes. Invoking the weak-strong uniqueness principle, we conclude that the DMV solution and classical solution of the incompressible Euler system coincide, proving the AP property of the scheme. We also present an extensive numerical case study in order to illustrate the theoretical convergences, wherein we utilize the techniques of K-convergence.
Numerical Analysis
What problem does this paper attempt to address?
The core problem that this paper attempts to solve is to develop and analyze a new Asymptotic Preserving (AP) finite - volume method for solving multi - dimensional compressible isentropic Euler equations to simulate low - Mach - number fluid flows. Specifically, the research aims to prove that the proposed numerical scheme has conditional energy stability and consistency, and strictly demonstrate the AP property of this scheme by using the concept of Dissipative Measure - Valued (DMV) solutions. ### Main problem decomposition 1. **Low - Mach - number limit problem**: - When the Mach number \( \varepsilon \to 0 \), the solutions of the compressible Euler equations will converge to the solutions of the incompressible Euler system. This process is called the zero - Mach - number limit or the low - Mach - number limit. - Through the framework of DMV solutions, the paper proves that under appropriate initial conditions, the numerical solutions will weakly converge to the classical solutions of the incompressible Euler system. 2. **Asymptotic preservation (AP property) of the numerical method**: - A new finite - volume method of upwind - stable scheme is proposed, in which the stabilizing term is proportional to the pressure gradient. - It is proved that with fixed grid parameters, this numerical method can approximate the incompressible Euler equations as the Mach number approaches zero; at the same time, with a fixed Mach number, as the grid parameters approach zero, it can approximate the compressible Euler equations. 3. **Energy stability and consistency**: - The energy stability and consistency of this numerical method are proved to ensure its reliability in numerical simulations. ### Formula summary - **Pressure potential energy**: \[ \psi_\gamma(\varrho) = \frac{\varrho^\gamma}{\gamma - 1}, \quad \gamma > 1 \] - **Relative internal energy**: \[ \Pi_\gamma(\varrho) = \psi_\gamma(\varrho) - \psi_\gamma(1) - \psi'_\gamma(1)(\varrho - 1) \] - **Total energy identity**: \[ \partial_t \left( \frac{1}{2} \varrho |u|^2 + \frac{1}{\varepsilon^2} \Pi_\gamma(\varrho) \right) + \nabla_x \cdot \left( \left( \frac{1}{2} \varrho |u|^2 + \frac{1}{\varepsilon^2} (\psi_\gamma(\varrho) - \psi'_\gamma(1)\varrho + p) \right) u \right) = 0 \] - **Modified momentum equation**: \[ \partial_t (\varrho u) + \nabla_x \cdot (\varrho u \otimes (u - \delta u)) + \frac{1}{\varepsilon^2} \nabla_x p = 0 \] ### Conclusion Through strict mathematical analysis and numerical experiments, this paper demonstrates the effectiveness and accuracy of the proposed Asymptotic Preserving finite - volume method in the low - Mach - number limit. By introducing the concept of DMV solutions, it not only proves the AP property of the numerical method, but also provides a new theoretical tool for understanding and dealing with low - Mach - number fluid flows.