Positivity-preserving cell-centered Lagrangian schemes for multi-material compressible flows: From first-order to high-orders. Part I: The one-dimensional case

Abstract : One of the main issues in the field of numerical schemes is to ally robustness with accuracy. Considering gas dynamics, numerical approximations may generate negative density or pressure, which may lead to nonlinear instability and crash of the code. This phenomenon is even more critical using a Lagrangian formalism, the grid moving and being deformed during the calculation. Furthermore, most of the problems studied in this framework contain very intense rarefaction and shock waves. In this paper, the admissibility of numerical solutions obtained by high-order finite-volume-scheme-based methods, such as the discontinuous Galerkin (DG) method, the essentially non-oscillatory (ENO) and the weighted ENO (WENO) finite volume schemes, is addressed in the one-dimensional Lagrangian gas dynamics framework. After briefly recalling how to derive Lagrangian forms of the 1D gas dynamics system of equations, a discussion on positivity-preserving approximate Riemann solvers, ensuring first-order finite volume schemes to be positive, is then given. This study is conducted for both ideal gas and non ideal gas equations of state (EOS), such as the Jones-Wilkins-Lee (JWL) EOS or the Mie-Grüneisen (MG) EOS, and relies on two different techniques: either a particular definition of the local approximation of the acoustic impedances arising from the approximate Riemann solver, or an additional time step constraint relative to the cell volume variation. Then, making use of the work presented in [89, 90, 22], this positivity study is extended to high-orders of accuracy, where new time step constraints are obtained, and proper limitation is required. Through this new procedure, scheme robustness is highly improved and hence new problems can be tackled. Numerical results are provided to demonstrate the effectiveness of these methods. This paper is the first part of a series of two. The whole analysis presented here is extended to the two-dimensional case in [85], and proves to fit a wide range of numerical schemes in the literature, such as those presented in [19, 64, 15, 82, 84].
Type de document :
Article dans une revue
Journal of Computational Physics, Elsevier, 2016, 312, pp.385 - 415. 〈10.1016/j.jcp.2016.02.027〉
Liste complète des métadonnées

Littérature citée [86 références]  Voir  Masquer  Télécharger

https://hal.inria.fr/hal-01248910
Contributeur : François Vilar <>
Soumis le : mardi 29 décembre 2015 - 12:39:02
Dernière modification le : jeudi 26 juillet 2018 - 12:09:52

Fichier

lag_pos_1D.pdf
Fichiers produits par l'(les) auteur(s)

Identifiants

Collections

Citation

François Vilar, Chi-Wang Shu, Pierre-Henri Maire. Positivity-preserving cell-centered Lagrangian schemes for multi-material compressible flows: From first-order to high-orders. Part I: The one-dimensional case. Journal of Computational Physics, Elsevier, 2016, 312, pp.385 - 415. 〈10.1016/j.jcp.2016.02.027〉. 〈hal-01248910〉

Partager

Métriques

Consultations de la notice

127

Téléchargements de fichiers

48