Flux-Balanced Patankar-type Schemes for the Compressible Euler Equations

TL;DR

This paper explores flux-balanced Patankar-type schemes for the compressible Euler equations, emphasizing the importance of applying the Patankar-trick selectively to preserve physical quantities and improve numerical accuracy.

Contribution

It introduces a balanced approach to applying Patankar schemes, showing that applying it only to the density equation yields better preservation of contact discontinuities and accuracy.

Findings

Applying Patankar-trick only to density improves positivity preservation.

Balanced fluxes maintain contact discontinuities.

The proposed method outperforms traditional schemes in accuracy and efficiency.

Abstract

Positivity preservation of key physical quantities in the context of fluid flows, such as density and internal energy, is an essential property of a numerical scheme as otherwise the solution lacks physical relevance and has a not well-defined equation of state. One time integration technique that is capable of preserving the positivity of quantities for every time step size is the Patankar-trick and its variants. However, in the context of the Euler equations of gas dynamics, we wonder whether the Patankar-trick should be applied to the density and total energy equations or only to one of them. In this work, we discuss one drawback of the schemes when blindly applied to every positive conserved variable and additionally point out how to overcome the issue by balancing the involved numerical fluxes correctly. To illustrate our findings, we investigate modified Patankar--Runge--Kutta (MPRK) schemes in the context of the compressible Euler equations with and without stiff source terms. We discover that it is beneficial to only apply the Patankar-trick in the density equation and to balance the remaining numerical fluxes consistently rather than applying the trick also to the energy equation. This leads also to the preservation of contact discontinuities. We perform numerical experiments to demonstrate that the accuracy of the methods is maintained while the performance of our approach is superior to the traditional application of MPRK schemes.