2024/09/09 by Saurav Samantaray, Samantaray, Saurav
Engineering · Mathematics · #35L45 #35L60 #35L65 #35L67 #65M06 #65M08 #Computational Fluid Dynamics and Aerodynamics #FOS: Mathematics #Gas Dynamics and Kinetic Theory #Navier-Stokes equation solutions #Numerical Analysis (math.NA)
paper · pdf · doi:10.48550/arxiv.2409.05854
openalex publication_date 2024/09/09 · openalex created_date 2024/10/22 · openalex updated_date 2026/07/28
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 L2 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.