An Asymptotic Preserving and Energy Stable Scheme for the Barotropic Euler System in the Incompressible Limit

K. R. Arun,Rahuldev Ghorai,Mainak Kar
2023-07-13
Abstract:An asymptotic preserving and energy stable scheme for the barotropic Euler system under the low Mach number scaling is designed and analysed. A velocity shift proportional to the pressure gradient is introduced in the convective fluxes, which leads to the dissipation of mechanical energy and the entropy stability at all Mach numbers. The resolution of the semi-implicit in time and upwind in space fully-discrete scheme involves two steps: solution of an elliptic problem for the density and an explicit evaluation for the velocity. The proposed scheme possess several physically relevant attributes, such as the positivity of density, the entropy stability and the consistency with the weak formulation of the continuous Euler system. The AP property of the scheme, i.e.\ the boundedness of the mesh parameters with respect to the Mach number and its consistency with the incompressible limit system, is shown rigorously. The results of extensive case studies are presented to substantiate the robustness and efficacy of the proposed scheme as well as the theoretical claims.
Numerical Analysis
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to design and analyze a numerical scheme that is both Asymptotic Preserving (AP) and energy - stable in the low - Mach - number limit for solving the barotropic Euler system. Specifically, the authors aim to overcome the serious problems encountered by traditional explicit time - stepping methods when dealing with low - Mach - number flows, such as high computational cost, numerical instability, and inaccurate solutions. ### Main problems 1. **Numerical stability in the low - Mach - number limit**: - In the case of low Mach numbers, the classical explicit time - stepping method, due to the CFL condition, results in a very small time step, thereby increasing the computational cost. - The low - Mach - number limit is a singular limit, and traditional explicit methods have difficulty handling the physical phenomena in this limit situation. 2. **Energy stability and entropy stability**: - Ensuring that the numerical scheme can maintain the stability of energy and entropy at all Mach numbers is crucial for guaranteeing the physical rationality and accuracy of the numerical solution. - Introduce a velocity shift proportional to the pressure gradient to ensure the dissipation of mechanical energy and entropy stability. 3. **Asymptotic fidelity (AP property)**: - The designed numerical scheme needs to be able to naturally transition to the incompressible limit system as the Mach number approaches zero and maintain the consistency and stability of the numerical solution. - It is necessary to prove the asymptotic fidelity of the scheme in the low - Mach - number limit, that is, the constraints on the time step and other discrete parameters do not depend on the Mach number. ### Solutions The authors propose a semi - implicit time - discretization method combined with an upwind finite - volume spatial - discretization. Specific steps include: - **Introducing velocity shift**: Introduce a velocity shift proportional to the pressure gradient in the convective flux to ensure the dissipation of mechanical energy and entropy stability. - **Semi - implicit time - discretization**: Overcome the stiffness - stability limitation through semi - implicit time - discretization, so that the time step can be larger, thereby reducing the computational cost. - **Solving elliptic equations**: Update the velocity by solving an elliptic equation for density to ensure the physical rationality of the scheme. - **Entropy stability**: Ensure that the numerical scheme satisfies the discrete entropy inequality, thereby guaranteeing the physical meaning and stability of the solution. ### Conclusion This paper successfully designs and analyzes a numerical scheme that is asymptotically faithful and energy - stable for solving the barotropic Euler system in the low - Mach - number limit. This scheme not only has a strict mathematical proof in theory but also demonstrates its robustness and high efficiency in multiple numerical experiments.