Thèse soutenue

Modèles de turbulence auto-adaptatifs pour la simulation des écoulements compressibles instationnaires
FR  |  
EN
Accès à la thèse
Auteur / Autrice : Grégoire Pont
Direction : Paola CinnellaJean-Christophe RobinetPierre Brenner
Type : Thèse de doctorat
Discipline(s) : Mécanique-matériaux
Date : Soutenance le 08/04/2015
Etablissement(s) : Paris, ENSAM
Ecole(s) doctorale(s) : École doctorale Sciences des métiers de l'ingénieur (Paris)
Partenaire(s) de recherche : Laboratoire : Laboratoire de Dynamique des Fluides (Paris) - Laboratoire de Dynamique des Fluides
Jury : Président / Présidente : Jean-Pierre Croisille
Examinateurs / Examinatrices : Paola Cinnella, Jean-Christophe Robinet, Pierre Brenner, Ivan Mary
Rapporteurs / Rapporteuses : Philippe Spalart, Maria Vittoria Salvetti

Résumé

FR  |  
EN

Cette thèse est principalement dédiée à la simulation des écoulements massivement décollés dans le domaine spatial. Nous avons restreint notre étude aux écoulements d'arrière-corps, pour lesquels ces décollements sont imposés par des changements brutaux de la géométrie. Dans le domaine spatial, le caractère fortement compressible des écoulements rencontrés impose l'utilisation de schémas numériques robustes. D'un autre coté, la simulation fine de la turbulence impose des schémas d'ordre élevé et peu dissipatifs. Ces deux spécifications, apparemment contradictoires, doivent pourtant coexister au sein d'une même simulation. Les modèles de turbulence ainsi que les schémas de discrétisation sont indissociables et leur couplage doit impérativement être considéré. Les schémas numériques doivent garder leur précision formelle dans des géométries complexes et des maillages très irréguliers imposés par le contexte industriel. Cette étude analyse le schéma de discrétisation utilisé dans le code de calcul FLUSEPA développé par Airbus Defence & Space. Ce schéma est robuste et précis pour des écoulements avec chocs et il présente une faible sensibilité au maillage (l'ordre 3 étant conservé même sur des maillages fortement perturbés). Malheureusement, le schéma possède une trop faible résolvabilité liée à un niveau de dissipation trop élevé pour envisager des simulations hybrides RANS/LES. Pour pallier à cet inconvénient, nous nous sommes penchés vers une solution basée sur un recentrage conditionnel et local : dans les zones dominées par des structures tourbillonnaires, une fonction analytique assure un recentrage local lorsque la stabilité numérique le permet. Cette condition de stabilité assure le couplage entre le schéma et le modèle. De cette manière, les viscosités laminaire et tourbillonnaire sont les seules à jouer un rôle dans les régions dominées par la vorticité et servent aussi à stabiliser le schéma numérique. Cette étude présente de plus une comparaison qualitative et quantitative de plusieurs modèles hybrides RANS/LES, à égalité de maillage et de schéma utilisés Pour cela, un certain nombre d'améliorations (notamment de leur capacité à résoudre les instabilités de Kelvin-Helmohlotz sans retard), proposées dans la littérature ou bien introduites dans cette thèse, sont prises en compte. Les applications numériques étudiées concernent des géométries allant de la marche descendante au lanceur spatial complet à échelle réduite.