Développement numérique d'une méthode Moment of Fluid/ Level Set. Application à l'atomisation de jet et nappe liquides

With the increase in the passenger air travel, it has become necessary to design fuel efficient engines for long-haul and ultra-long-haul flights. Because of the continuous progress in the supercomputing power and advances in the computational fluid dynamics (CFD) methods, numerical simulations have been sought after as the choice for investigating the physical processes occurring inside the aircraft engines. Of the many processess, atomization of the liquid fuel, i.e., the process by which the injected liquid fuel breaks up into droplets remains imperative to be completely understood. Since atomization processgoverns the size of the fuel droplets produced, it has a direct influence on the evaporation rate, completeness of combustion, and even pollutant formation. However, due to the multiscale, multiphysical, and multiphase aspect of this process, it has become a challenge to numerically simulate it. Most often, the simulations run into under-resolution limitation when capturing the droplets. To mitigate this shortcoming and to make the simulations numerically tractable, this work presents two numerical methods of liquid/gas interface reconstruction to capture the liquid droplets {moment of fluid (MOF) method and hybrid moment of fluid{level set (HyMOFLS) method. These methods are coupled with consistent mass and momentum flux computation as well as ghost fluid method (GFM) for handling discontinuities in density and jump in pressure across the interface. The MOF method uses liquid volume fraction as well as liquid and gas phase centroids for interface reconstruction in each computational cell in the simulations. The advantage in using the phase centroids is that the neighbor cell data are not required in the interface reconstruction process resulting in a uniform treatment of interior and boundary cells in the computational domain. This method improves the liquid/gas interface orientation and reconstruction in the underresolved regions of the domain. The HyMOFLS method combines the MOF method and CLSVOF method such that MOF method is employed to capture under-resolved regions and CLSVOF method for resolved regions of the interface. The switch between the choice of these methods in the computational domain is made according to local mesh spacing and curvature of the interface. This method strikes a good balance between the reconstruction accuracy and modest computational cost requirement compared to MOF method. Hence, it is a natural choice for performing simulations of primary atomization at aircraft engine operating conditions. The HyMOFLS method is employed in to simulate primary atomization of liquid for two injection con_gurations used in aircraft engines under relevant operating conditions. First, a planar pre-filming Airblast atomization configuration is simulated using two gridresolutions and inlet velocity profile to investigate and analyze their effect on the atomization characteristics, i.e., droplet and ligament properties. Results suggest that the simulations are matching satisfactorily with the experiments and are of the same order of magnitude as the experimental data. Next, primary atomization of a turbulent liquid jet in gaseous crossow configuration is simulated under three di_erent (low, moderate,and high) density ratio operating conditions using three different mesh resolutions. The analyses of the results yielded that there is low probability of occurrence of secondary atomization due to aerodynamic forces, the jet bending and penetration decrease with increase in density ratio, and the wavelength of the instability waves formed on the windward side of the liquid jet decreases from low to moderate density ratio and increases from moderate to high density ratio.Avec l'intensification du trafic aérien, il est devenu nécessaire de concevoir des moteurs économes en carburant pour les vols longs et ultra-long-courriers. Avec l'augmentation de la puissance de calcul et les avancées des méthodes de calcul en dynamique des fluides (CFD), les simulations numériques sont essentielles pour étudier les processus physiques mis en jeu dans les moteurs d'avion. Parmi ces nombreux processus, l'atomisation du carburant liquide reste d'une très grande importance. En effet il détermine la taille des gouttelettes de carburant, et il a une influence directe sur le taux d'évaporation, l'optimisation de la combustion ainsi que sur la formation de polluants. Cependant, en raison de l'aspect multi-échelle, multi-physique et multiphasique de ce processus, il reste complexe de le simuler numériquement. En particulier, les simulations se heurtent à une limitation de la résolution lors de la capture des gouttelettes. Afin de repousser cette limitation, ce travail présente deux méthodes numériques de reconstruction de l'interface liquide/gaz pour capturer les gouttelettes de liquide : la méthode Moment Of Fluid (MOF) et une méthode Hybrid Moment Of Fluid-Level Set (HyMOFLS). La méthode MOF utilise la fraction de volume liquide ainsi que leur centroïde en phase liquide et gazeuse pour la reconstruction de l'interface dans chaque cellule. L'avantage de l'utilisation des centroïdes de phase est que les données des cellules voisines ne sont pas nécessaires dans le processus de reconstruction de l'interface, ce qui permet un traitement uniforme de toutes les cellules dans le domaine du calcul. De plus, la rupture numérique de structures sous-résolues est prise en compte, ce qui n'est pas le cas avec la méthode CLSVOF. La méthode HyMOFLS associe la méthode MOF et la méthode CLSVOF; la méthode MOF est utilisée dans les régions sous-résolues et la méthode CLSVOF pour les régions résolues de l'interface. Le passage entre ces méthodes se fait en fonction de la taille de mailles et de la courbure de l'interface. Cette méthode propose un compromis entre la précision de la reconstruction et le coût de calcul par rapport à la méthode MOF seule. Après de nombreuses validations, cette méthode est utilisée pour simuler deux configurations d'injection utilisées dans les moteurs d'avion dans des conditions de fonctionnement appropriées. Tout d'abord, une configuration d'atomisation d'un film plan (Airblast) est simulée pour différentes tailles de mailles et différents profils de vitesse d'entrée afin d'analyser leur influence sur l'atomisation, en particulier sur les caractéristiques des gouttelettes et des ligaments. On observe que les simulations se comparent de manière satisfaisante aux expériences et sont quantitativement du même ordre de grandeur que les données expérimentales. Ensuite, l'atomisation primaire d'un jet liquide turbulent dans une configuration Crossow est simulée pour trois rapports de masses volumiques différents (faible, modéré et élevé) et en utilisant trois tailles de mailles différentes. Les analyses des résultats ont montré une faible probabilité d'atomisation secondaire due aux forces aérodynamiques, ainsi qu'une diminution de la flexion et de la pénétration du jet avec l'augmentation du rapport de masse volumique, et enfin une diminution suivie d'une augmentation de la longueur d'onde des perturbations formées sur le dessus du jet liquide avec la croissance du rapport des masse volumique