Thèse soutenue

Simulation numérique tridimensionnelle de l'aérodynamique interne non réactive des moteurs à allumage commandé par une méthode hybride RANS / LES

FR  |  
EN
Auteur / Autrice : Al Hassan Afailal
Direction : Rémi Manceau
Type : Thèse de doctorat
Discipline(s) : Mécanique des fluides et énergétique
Date : Soutenance le 16/12/2020
Etablissement(s) : Pau
Ecole(s) doctorale(s) : École doctorale sciences exactes et leurs applications (Pau, Pyrénées Atlantiques ; 1995-)
Partenaire(s) de recherche : Laboratoire : Laboratoire de mathématiques et de leurs applications (Pau) - Laboratoire de Mathématiques et de leurs Applications [Pau]
Equipe de recherche : IFP Energies Nouvelles
Jury : Président / Présidente : Jacques Borée

Résumé

FR  |  
EN

L'aérodynamique interne est un élément fondamental pour améliorer la combustion dans les moteurs à allumage commandé. Une meilleure maitrise des écoulements internes est permise grâce aux outils de simulation CFD qui sont de plus en plus utilisés dans le processus de développement des moteurs à allumage commandé. Cette thèse avait pour objectif d’étendre l'approche hybride RANS/LES-temporelle dite HTLES, initialement dédiée pour des écoulements statistiquement stationnaires, aux écoulements moteurs avec des parois mobiles et des modes opératoires cycliques, puis de la valider dans des configurations représentatives des écoulements moteurs. Cette approche vise à modéliser les régions proches parois par approche statistique RANS et tend continûment vers la LES temporelle loin des parois si la discrétisation spatiale et temporelle est suffisamment résolue. Le formalisme temporel permet une hybridation RANS/LES consistante dans un écoulement statistiquement stationnaire, les deux méthodes se basant sur des opérateurs temporels (respectivement la moyenne temporelle et le filtrage temporel). Une première amélioration de l’approche HTLES a été proposée en ajoutant une fonction de protection qui impose le mode RANS dans la région proche paroi, indépendamment de la discrétisation locale (spatiale et temporelle). Dans les écoulements cycliques, l’approche HTLES modélise les échelles turbulentes non-résolues en se basant sur des moyennes de phase des grandeurs résolues qui sont inconnues lors de la simulation. La moyenne glissante exponentielle (EWA) a été utilisée afin d’approximer ces moyennes de phase. Une formule pour définir la largeur de la moyenne glissante a été proposée de sorte que les fluctuations turbulentes (hautes fréquences) soient filtrées des quantités résolues, tout en conservant les composantes cycliques (basses fréquences). Cette approche a été implémentée dans le code de calcul industriel CONVERGE CFD. Elle a d'abord été validée dans deux configurations stationnaires : un canal plan infini et un banc volute. A cet effet, les résultats ont été comparés aux données de référence et aux résultats RANS et LES. Dans les régions proches parois où le maillage est sous résolu pour la LES, EWA-HTLES a mieux prédit l’écoulement grâce à l'utilisation du mode RANS, permettant une meilleure prédiction des pertes de charge. La résolution des grandes échelles dans la région centrale a permis d'obtenir des prédictions aussi précises qu’une simulation LES en termes de vitesses moyennes et des fluctuations. La validation de l'EWA-HTLES a également été effectuée dans deux configurations moteurs : le tumble compressé et le moteur Darmstadt, tous deux présentant des caractéristiques aérodynamiques typiques aux moteurs à allumage commandé telles que la génération et la compression du mouvement de tumble et la variabilité cyclique. Pour chaque configuration, un nombre total de 40 cycles consécutifs simulés à l'aide de EWA-HTLES a été utilisé pour calculer les deux premiers moments statistiques. Les résultats ont été comparés aux données de la PIV, et aux résultats donnés par les simulations RANS et LES. Les résultats ont montré que le modèle développé arrive à contrôler correctement la transition RANS-LES dans des configurations complexes avec des conditions d'écoulement non stationnaires et des déformations géométriques importantes, assurant le mode RANS aux parois et la LES au centre du cylindre. La résolution des grandes échelles a permis une bonne prédiction des phénomènes instationnaires, particulièrement l'évolution des caractéristiques du mouvement de tumble et des phénomènes associés aux variabilités cycliques, tels que l'augmentation locale de vitesses fluctuantes. Les résultats de l'EWA-HTLES sont similaires à ceux prédits par la LES et meilleurs que ceux donnés par les simulations RANS. Ces résultats montrent des perspectives encourageantes pour l'application de cette méthode dans de nombreuses configurations industrielles.