CE TRAVAIL PORTE SUR LA MODELISATION ET LA SIMULATION D'ECOULEMENTS REACTIFS TURBULENTS A CHIMIE NON INFINIMENT RAPIDE PAR UNE FONCTION DENSITE DE PROBABILITE CONJOINTE. LA RESOLUTION DE L'EQUATION D'EVOLUTION TEMPORELLE DE LA PDF CONJOINTE ANALYTIQUEMENT OU NUMERIQUEMENT EST TRES COMPLEXE. EN FAIT, LE PROBLEME EST SIMPLIFIE EN INTEGRANT CETTE EQUATION DANS L'ESPACE DES VITESSES. ON OBTIENT AINSI UNE APPROCHE HYBRIDE PROBABILISTE EULERIENNE LAGRANGIENNE. UNE METHODE DE MONTE-CARLO EST UTILISEE ENSUITE, POUR SIMULER CETTE EQUATION. DANS CETTE EQUATION, LE PHENOMENE DE MELANGE A PETITE ECHELLE DOIT ETRE REPRESENTE PAR UN MODELE. LA COMPARAISON D'UN NOUVEAU MODELE DE MELANGE AVEC DES MODELES CLASSIQUES D'ECHANGE AVEC LA MOYENNE A ETE FAITE DANS LE CAS D'UNE COUCHE DE MELANGE THERMIQUE ET DANS LE CAS D'UNE COUCHE DE MELANGE REACTIVE. LES RESULTATS DES CALCULS SONT ENCOURAGEANTS. LA POSSIBILITE DE L'APPLIQUER A DES CONFIGURATIONS INDUSTRIELLES COMPLEXES 2D AXISYMETRIQUE (LE BLUFF-BODY) ET 3D (CHAMBRE DE COMBUSTION DEXTRE) EST DEMONTREE

[Résumé en français] La simulation numérique de réacteurs industriels promet de devenir un outil précieux aidant à prospecter de nouvelles technologies de mélange, à dimensionner et à optimiser les conditions de fonctionnement du réacteur. Ces calculs deviennent difficiles dans le cas de réactions rapides (polymérisations anioniques, précipitations) où existe un fort couplage entre la réaction chimique et le mélange. Ce travail de thèse a pour but d'évaluer et de valider la simulation d'écoulements réactifs pour des applications de génie chimique avec une méthode que nous avons choisie pour sa capacité à tenir compte des interactions entre processus de mélange turbulent et réactions chimiques rapides. La méthode sélectionnée est une méthode basée sur la résolution numérique des équations de Navier-Stokes en milieu turbulent et l'équation de transport de la fonction de densité de probabilité jointe de la composition (PDF) à l'aide d'un algorithme de Monte-Carlo. L'évaluation s'est faite sur un réacteur de précipitation semi-batch de type Taylor-Couette. Chaque étape du calcul a été validée expérimentalement. Au terme de ce travail, il a été montré que l'emploi de ce type d'approche pour le calcul d'un réacteur industriel est hasardeux et délicat. Il demande une validation de l'écoulement, le calage de 2 constantes de modélisation et la connaissance précise des cinétiques de réaction. Enfin, l'utilisation d'un algorithme de type Monte-Carlo engendre un bruit numérique qui selon le type d'écoulement peut nuire à la fiabilité des résultats. Le cas traité n'a pas donné des résultats totalement conformes à l'expérience. La méthode numérique utilisée n'est pas pour autant remise en question car le cas traité était par nature délicat

Les limitations récentes imposées par ICAO-CAEP, qui règlent les émissions de NOx, mènent à l'implémentation du concept de combustion lean dans les moteurs aéronautiques. Du point de vue du design, il faudrait étudier de façon plus approfondie la combustion lean et donc la Computational Fluid Dynamics(CFD) peut être un outil essentiel à ce but. Plusieurs phénomènes sont impliqués et différentes stratégies de modélisation, avec des différences en termes de coûts de calcul, sont disponibles. Néanmoins, jusqu'à présent, peu d'outils numériques peuvent prendre en compte les effets de la préparation du combustible liquide dans les calculs réactifs. Les conditions limites d'atomisation sont normalement déterminées par des approches corrélatives qui ne couvrent pas toutes les conditions de fonctionnement et les caractéristiques géométriques des brûleurs aéronautiques. Cependant, comme on peut lire dans la première partie de ma thèse, où plusieurs études de cas de littérature sont analysées, l'impact de la préparation du combustible liquide peut être extrêmement important. Les considérations basées sur des approches corrélatives ne sont pas fiables. Des méthodes prédictives focalisées sur l'atomisation du combustible sont nécessaires. Cette activité de recherche vise donc à développer un outil numérique général, utilisable dans le domaine industriel et capable de modéliser la phase liquide de son injection jusqu'à la génération d'un spray dispersé. Le modèle ELSA (Eulerian Lagrangian Spray Atomization), basé sur une approche eulérienne dans la région dense et une lagrangienne dans la zone diluée, a été choisi à ce but. Le solveur traite le combustible liquide jusqu'à la génération d'une phase dispersée et prend en compte le processus d'atomisation par l'introduction de la densité d'interface liquide-gaz. Néanmoins, si on applique cette méthode dans un environnement réactif fortement turbulent comme un brûleur aéronautique on peut rencontrer plusieurs limites. Par conséquent, on a dédié une attention particulière tout d'abord à l'étude du terme de flux turbulent à l'intérieur de l'équation de la fraction volumique liquide. Cette quantité est d'une importance primaire pour un écoulement turbulent, avec des vitesses de glissement élevées entre ses phases. Une nouvelle fermeture de second ordre pour cette variable est proposée et validée sur un cas de jet en crossflow. Ensuite, pour gérer un environnement réactif, un nouveau modèle d'évaporation est intégré dans le code et évalué par rapport aux résultats expérimentaux. Enfin, une autre méthode de dériver la distribution de la taille des gouttes dans le contexte ELSA pour l'injection lagrangienne est présentée et validée avec des simulations DNS. Pour conclure, ce travail introduit une nouvelle méthode pour une description unifiée de la combustion et de l'atomisation dans les calculs CFD. L'approche proposée devrait conduire à une description complète de l'évolution du combustible et à une caractérisation plus pertinente de l'écoulement réactif. Plusieurs aspects qui sont également mis en évidence sont améliorables et peuvent offrir des suggestions pour d'ultérieurs travaux.