19–21 sept. 2022
Institut Henri Poincaré - Paris
Fuseau horaire Europe/Paris

Paola Allegrini (Institut de Mathématiques de Toulouse) - Un schéma numérique préservant l’asymptotique bas Mach pour les équations d’Euler

20 sept. 2022, 11:00
45m
Amphi Hermite (Institut Henri Poincaré - Paris)

Amphi Hermite

Institut Henri Poincaré - Paris

Description

Ce travail est réalisé en collaboration avec Marie-Hélène Vignal, Institut de Mathématiques de Toulouse, UT3-Paul Sabatier.

Dans cet exposé je m’intéresse au développement et à l’étude d’un schéma asymptotiquement préservant (AP) pour les équations d’Euler compressible dans la limite bas Mach.
Pour les écoulements hautement subsoniques, les ondes acoustiques étant très rapides par rapport à la vitesse du fluide, le gaz peut être considéré comme incompressible. D’un point de vue numérique, lorsque le nombre de Mach tend vers zéro, les schémas classiques explicites présentent deux inconvénients majeurs : ils perdent en consistance et imposent une contrainte sur le pas de temps très restrictive.
Nous avons développé un schéma asymptotiquement stable, c’est-à-dire soumis à une CFL indépendante du nombre de Mach, et asymptotiquement consistant, c’est-à-dire consistant avec le modèle incompressible à la limite bas Mach.
Ce type de schéma a été très étudié dans la littérature, notamment pour Euler isentropique [6, 5, 4, 7] mais aussi pour Euler complet [3, 2] avec diverses méthodes : schémas Lagrangiens, maillages décalés ou colocalisés.
Dans ce travail, nous proposons un schéma AP basé sur une discrétisation IMEX (Implicite-Explicite) pour la partie en temps et volumes finis colocalisés pour la partie spatiale.
Je présenterai notre schéma AP, son extension à l’ordre 2 et la procédure MOOD utilisée pour réduire les oscillations (problème classique de la montée en ordre). Enfin, je finirai mon exposé avec la présentation de quelques résultats préliminaires pour les équations de Navier-Stokes.

[1] P. Allegrini, M.-H. Vignal, Study of a low oscillatory second order all Mach IMEX scheme for the full Euler equations, soumis, 2022.

[2] W. Boscheri, G. Dimarco, R. Loubère, M. Tavelli, M.-H. Vignal, A second order all Mach number IMEX finite volume solver for the three dimensional Euler equations, J. Comp. Phys. 415, 2020.

[3]  M. Dumbser, V. Casulli, A conservative, weakly nonlinear semi-implicit finite volume scheme for the compressible Navier-Stokes equations with general equation of state, Applied Mathematics and Computation, vol. 272, issue P2, 479-497, 2016. 


[4]  G. Dimarco, R. Loubère, V. Michel-Dansac, M.-H. Vignal, Second order Implicit-Explicit Total Variation Diminishing schemes for the Euler system in the low Mach regime, J. Comput. Phys., 372, 178–201, 2018. 


[5]  G. Dimarco, R. Loubère, M.-H. Vignal, Study of a new asymptotic preserving scheme for the Euler system in the low Mach number limit, SIAM J. Sci. Comput., 39(5), A2099–A2128, 2017.

[6]  P. Degond, M. Tang, All speed scheme for the low Mach number limit of the isentropic Euler equations. Commun. Comput. Phys. 10(1), 1–31, 2011. 


[7]  P. Bruel, S. Delmas, J. Jung, V. Perrier, A low Mach correction able to deal with low Mach acoustics, J. Comp. Phys. 378, 723-759, 2019.

Documents de présentation