Skip to Main content Skip to Navigation
Journal articles

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
LPP - Laboratoire Paul Painlevé - UMR 8524, Inria Lille - Nord Europe
3 COFFEE - COmplex Flows For Energy and Environment
CRISAM - Inria Sophia Antipolis - Méditerranée , JAD - Laboratoire Jean Alexandre Dieudonné : UMR7351
4 ANGE - Numerical Analysis, Geophysics and Ecology
LJLL - Laboratoire Jacques-Louis Lions, Inria Paris-Rocquencourt
Abstract : 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.
Complete list of metadata

Cited literature [34 references]  Display  Hide  Download
Contributor : Caterina Calgaro Connect in order to contact the contributor
Submitted on : Sunday, January 13, 2013 - 4:28:10 PM
Last modification on : Thursday, August 4, 2022 - 4:54:52 PM
Long-term archiving on: : Saturday, April 1, 2017 - 4:18:25 AM


Files produced by the author(s)



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



Record views


Files downloads