Positivity-preserving schemes for Euler equations: Sharp and practical CFL conditions

Caterina Calgaro (1, 2) , Emmanuel Creusé (1, 2) , Thierry Goudon (3) , Yohan Penel (1, 4)

1 SIMPAF - SImulations and Modeling for PArticles and Fluids
2 LPP - Laboratoire Paul Painlevé - UMR 8524
3 COFFEE - COmplex Flows For Energy and Environment
4 ANGE - Numerical Analysis, Geophysics and Ecology

When one solves PDEs modelling physical phenomena, it is of great importance to take physical constraints into account. More precisely, numerical schemes have to be designed such that discrete solutions satisfy the same constraints as exact solutions. For instance, the underlying physical assumptions for the Euler equations are the positivity of both den- sity and pressure variables. We consider in this paper an unstructured vertex-based tesselation in R2 . Given a MUSCL finite volume scheme and given a reconstruction method (including a limiting process), the point is to determine whether the overall scheme ensures the positivity. The present work is issued from seminal papers from Perthame and Shu (On positivity preserving finite volume schemes for Euler equations, Numer. Math. 73 (1996) 119-130) and Berthon (Robustness of MUSCL schemes for 2D unstructured meshes, J. Comput. Phys. 218 (2) (2006) 495-509). They proved in different frameworks that under assumptions on the cor- responding one-dimensional numerical flux, a suitable CFL condition guarantees that den- sity and pressure remain positive. We first analyse Berthon's method by presenting the ins and outs. We then propose a more general approach adding non geometric degrees of freedom. This approach includes an optimization procedure in order to make the CFL condition explicit and as less restric- tive as possible. The reconstruction method is handled independently by means of s- limiters and of an additional damping parameter. An algorithm is provided in order to specify the adjustments to make in a preexisting code based on a certain numerical flux. Numerical simulations are carried out to prove the accuracy of the method and its ability to deal with low densities and pressures.

Mots clés: Positivity-preserving MUSCL schemes; 2nd-order finite volume method; Euler equations; CFL condition; τ-limiters.; τ-limiters

Domaines: Equations aux dérivées partielles [math.AP]; Analyse numérique [math.NA]

Caterina Calgaro, Emmanuel Creusé, Thierry Goudon, Yohan Penel. Positivity-preserving schemes for Euler equations: Sharp and practical CFL conditions. Journal of Computational Physics, 2013, 234 (1), pp.417--438. ⟨10.1016/j.jcp.2012.09.040⟩. ⟨hal-00768479v2⟩