Asymptotic Preserving Linearly Implicit Additive IMEX-RK Finite Volume Schemes for Low Mach Number Isentropic Euler Equations

Saurav Samantaray

2024-11-06

Abstract:We consider the compressible Euler equations of gas dynamics with isentropic equation of state. Standard numerical schemes for the Euler equations suffer from stability and accuracy issues in the low Mach regime. These failures are attributed to the transitional behaviour of the governing equations from compressible to incompressible solution in the limit of vanishing Mach number. In this paper we introduce an extra flux term to the momentum flux. This extra term is recognised by looking at the constraints of the incompressible limit system. As a consequence the flux terms enable us to get a suitable splitting, so that an additive IMEX-RK scheme could be applied. Using an elliptic reformulation the scheme boils down to just solving a linear elliptic problem for the density and then explicit updates for the momentum. The IMEX schemes developed are shown to be formally asymptotically consistent with the low Mach number limit of the Euler equations and are shown to be linearly $L^2$ stable. A second order space time fully discrete scheme is obtained in the finite volume framework using a combination of Rusanov flux for the explicit part and simple central differences for the implicit part. Results of numerical case studies are reported which elucidate the theoretical assertions regarding the scheme and its robustness.

Numerical Analysis

What problem does this paper attempt to address?

This paper aims to solve the problems of standard numerical methods in terms of stability and accuracy in low - Mach number flows. Specifically:

1. **Stability problems**: In the low - Mach number limit, due to the limitation of the CFL condition, the time step of the standard explicit scheme needs to be proportional to the Mach number, which makes the solver almost impossible to evolve in practical applications, especially in the incompressible fluid region.

2. **Accuracy problems**: Numerical methods have difficulties in handling the transition behavior from compressible to incompressible, resulting in inaccurate solutions in the incompressible region. These inaccuracies mainly stem from the fact that numerical methods do not fully respect the transition characteristics of the equations, thus losing information in the process from continuous to discrete.

To solve these problems, the author introduced an additional flux term to improve the momentum flux and proposed a linearly implicit additive IMEX - RK finite - volume scheme. This scheme has the following characteristics:

- **Asymptotic Preserving (AP)**: The scheme can automatically be converted into a solver for incompressible fluids in the low - Mach number limit, and its stability requirements are independent of multi - scale or singular parameters.
- **Linearly implicit**: By introducing an additional flux term, the pressure term can be processed linearly and implicitly, avoiding the complexity of nonlinear solutions.
- **High - order accuracy**: By combining the Rusanov flux and the central difference, a second - order fully discrete finite - volume framework in space and time is constructed.

The paper also carried out numerical case studies to verify the correctness of the theoretical analysis and the robustness of the scheme.