Les turbines aéronautiques doivent satisfaire à des normes d’émissions polluantes toujours en baisse. La qualité du mélange du carburant et de l’air dans la chambre de combustion est responsable de la formation de polluants nocifs pour l’environnement. La simulation aux grandes échelles (LES) permet d’étudier les mécanismes de mélanges turbulents de l’air et du carburant. La prise en compte de l’aspect liquide du carburant injecté devient nécessaire pour prédire correctement l’apparition de vapeur de carburant au sein du foyer. Le but de cette thèse est évaluer la fiabilité des simulations LES Euler-Euler dans une configuration complexe. Les processus d’injection, et d’évaporation du carburant liquide sont analysés et modélisés dans les simulations LES car ils pilotent la formation de vapeur de carburant. Les méthodes numériques pour résoudre les équations continues de la phase dispersée doivent permettre des simulations précises et robustes dans une configuration représentative d’une chambre de combustion. Les simulations présentées dans ces travaux reproduisent l’écoulement diphasique évaporant non-réactif du banc d’essai Mercato. Ce banc est équipé d’un système d’injection d’air vrillé et d’un atomiseur pressurisé-swirlé de kérosène typiques des foyers aéronautiques réels. Dans ces travaux, le modèle pour l’injection de liquide FIM-UR a été développé pour définir les conditions limites conduisant à un spray issu d’un atomiseur préssurisé-swirlé. Le kérosène employé dans les campagnes expérimentales est modélisé dans les simulations par un composé permettant d’obtenir des temps d’évaporation réalistes. Trois stratégies numériques ont été mises en place sur la configuration Mercato. Les comparaisons des résultats numériques aux mesures expérimentales ont permis d’évaluer la stratégie numérique conduisant à la meilleure précision. L’utilisation du schéma centré TTGC associé à un opérateur de viscosité artificielle localisée par un senseur adapté est optimale lorsque l’équation sur l’énergie décorrélée des gouttes est résolue. Cette stratégie permet de contrôler la localisation et les niveaux de viscosité par rapport à un schéma décentré. Les termes sources liés au mouvement mésoscopique permettent de redistribuer l’énergie dans les zones de compression ou de détente de la phase dispersée, et d’obtenir les bonnes répartitions des fluctuations dans la chambre de combustion. La stratégie retenue est comparée aux statistiques de la dynamique du spray résolu par une approche Lagrangienne employant la même injection monodispersse. Le méthode Euler-Euler conduit à la même précision de la dynamique de la phase dispersée que la méthode Euler-Lagrange. L’accès à l’évolution instationnaire de l’écoulement permet d’identifier les mêmes mécanismes de dispersion et de mélange dans les deux simulations. Des différences sur la répartition de diamètre moyen et de carburant dans la chambre ont été mis en évidence et reliés à la polydispersion locale qui n’est pas résolue dans l’approche Euler-Euler monodisperse et qui apparaît naturellement dans l’approche Euler-Lagrange malgré l’injection monodisperse. / Aeronautical gas turbines are facing growing demands on emission reductions. Indeed, the quality of the air-fuel mixture directly triggers the formation of pollutants degrading the environment. Large Eddy Simulation is an accurate numerical method to predict turbulent mixing in combustors. Adding the liquid phase of the fuel in these simulations also becomes necessary to properly predict the injection process and the vaporization of the fuel in the combustion chamber. The purpose of this dissertation is to evaluate the accuracy and reliability of Euler-Euler LES in a complex combustor configuration. The injection and vaporization processes of the fuel liquid phase are both modeled in the present LES as they drive the formation of the fuel gas phase. Moreover, the numerical methods that solve the continuous equations of the disperse phase must be accurate and robust in realistic combustor configurations. The simulations shown in the present study reproduce the non-reactive two-phase flow of the ONERA Mercato test bench. The experimental set-up is equipped with an air-swirler injection system and a pressure-swirled atomizer typical of actual turboengine combustors. In the present work the FIM-UR liquid injection model has been developed. It creates boundary conditions profiles for a liquid spray produced by a pressure-swirled atomizer. Kerosene used in the experiments is modeled in the present numerical simulations by a single species leading to a good estimate of the vaporization rate. Three numerical strategies have been tested on the Mercato configuration. Comparisons between experimental and LES results help defining the most accurate numerical strategy. The use of the centered numerical scheme TTGC stabilized by a localized artificial viscosity operator is best when the random uncorrelated energy of droplets is also resolved. Unlike an upwind numerical scheme, the selected strategy allows the user to control where and how much artificial viscosity is added. The source terms coming from the mesoscopic movement redistribute the energy in the compression or expanding zones of the disperse phase, and provide the proper distribution of fluctuations in the combustion chamber. The obtained strategy is compared with the statistics provided by a Lagrangian description of the liquid spray in the same mono-disperse injection. The Euler-Euler approach leads to the same accuracy in the same spray dynamics of the disperse phase as in the Euler-Lagrange method. Both unsteady flow simulations also provide the same dispersion and mixing processes in the Mercato set-up. Differences on the mean diameter and the fuel distribution in the combustion chamber are seen and related to the local poly-dispersion that cannot be resolved in the mono-disperse Euler-Euler approach and that naturally appear in the Euler-Lagrange method despite the mono-disperse injection.

In this paper we present detailed Euler-Euler Large Eddy Simulations (LES) of dispersed bubbly flow in a rectangular bubble column. The motivation of this study is to investigate potential of this approach for the prediction of bubbly flows, in terms of mean quantities. The set of physical models describing the momentum exchange between the phases was chosen according to previous experiences of the authors. Experimental data, Euler-Lagrange LES and unsteady Euler-Euler Reynolds-Averaged Navier-Stokes model are used for comparison. It was found that the presented modelling combination provides good agreement with experimental data for the mean flow and liquid velocity fluctuations. The energy spectrum made from the resolved velocity from Euler-Euler LES is presented and discussed.

In this paper we present detailed Euler-Euler Large Eddy Simulations (LES) of dispersed bubbly flow in a rectangular bubble column. The motivation of this study is to investigate potential of this approach for the prediction of bubbly flows, in terms of mean quantities. The set of physical models describing the momentum exchange between the phases was chosen according to previous experiences of the authors. Experimental data, Euler-Lagrange LES and unsteady Euler-Euler Reynolds-Averaged Navier-Stokes model are used for comparison. It was found that the presented modelling combination provides good agreement with experimental data for the mean flow and liquid velocity fluctuations. The energy spectrum made from the resolved velocity from Euler-Euler LES is presented and discussed.

Ce travail de thèse rentre dans le cadre du développement de modèles multi‐échelles de simulation de colonne d’absorption gaz‐liquide pour des applications de captage de CO2 en vue d’optimiser leur design. Il est le fruit d’une collaboration entre IFP Énergies nouvelles et l’Institut de Mécanique des Fluides de Toulouse. Les colonnes à garnissages représentent une technologie essentielle aux applications d’absorption gaz‐liquide. Dans les procédés de captage de CO2 aux amines, le solvant liquide s’écoule sur les parois du garnissage, idéalement sous forme de film ruisselant mouillant toute la surface disponible, le gaz, en régime turbulent, venant le cisailler à contre‐courant de manière à promouvoir un transfert de CO2 de la phase gaz vers la phase liquide. Un écoulement le plus homogène possible permet d’avoir les meilleures performances de transfert. Toutefois, l’expérience montre que des maldistributions, notamment de la phase liquide, peuvent apparaître même en cas de bonnes distributions en entrée de colonne. La distribution du liquide est régie par un phénomène de « dispersion » dont l’étude et la modélisation représentent le principal objet de cette thèse. Pour ce faire, ce travail de thèse s’appuie sur des travaux expérimentaux, réalisés sur une installation d’IFPEN à Lyon, et sur des travaux numériques réalisés dans l’équipe Interface de l’IMFT. Le premier axe de l’étude abordé a ainsi consisté en l’acquisition de données originales de distribution de liquide en partant d’une configuration d’alimentation sous forme d’un jet central en tête de colonne et ce pour deux types de garnissages métalliques : un garnissage structuré, le Mellapak 250.X et un garnissage vrac, l’IMTP. 40. La méthode de tomographie gamma a été mise en œuvre afin de mesurer l’atténuation d’un flux photonique par le liquide en mouvement ce qui permet d’établir des cartes de rétention locale de liquide sur une section de colonne. Les profils de rétention résultants ont été ensuite exploités afin de caractériser la dispersion de liquide dans le système pour des régimes d’écoulement allant des plus faibles aux plus fortes interactions gaz‐liquide. Ces résultats ont permis de développer un modèle simple d’advection diffusion faisant appel à un paramètre hydrodynamique clé qui est le « coefficient de dispersion », qui reproduit bien l’étalement du jet de liquide. Dans le cas du garnissage Mellapak 250.X nous avons pu mettre en évidence une dispersion qui est régie essentiellement par la géométrie du milieu. / This work is done within the framework of collaboration between IFPEN and l’Institut de Mécanique des Fluides de Toulouse (IMFT). It takes place in the general context of optimization of industrial gas‐liquid packed columns dedicated to CO2 capture using multi‐scale simulation approach of hydrodynamics and transfer. The main objective of this study is to develop predictive model for liquid dispersion in two‐phase flows in packed beds operating in the counter‐current mode. This model will help simulate the flow at the macro‐scale of a packed column. Packed columns are widely used for gas/liquid absorption processes since they generate subsequent exchange surface between phases with limited pressure drop. In the particular case of amine process, the two‐phase flow in the packing consists in thin trickling liquid films sheared by counter‐current gas flow that circulates in communicating tortuous channels. It is widely recognized that homogeneous flow ensures good separation performances. This is far from being the case at industrial scale since maldistributions, especially for liquid phase, do occur even with optimized liquid and gas distributors in the column inlets. Liquid distribution results from “dispersion” phenomenon which modelling is not fully handled in literature. Prediction of liquid dispersion has been the main objective of this Ph.D. research. This study is based upon two main approaches: an experimental study carried out using a 40 cm in diameter column operating in IFPEN (Lyon) and a numerical study carried out in Interface team in IMFT. The first part of this research focused on measuring spatial distribution of liquid hold‐up over different sections of the packed bed using a gamma ray tomography system. In order to focus into liquid dispersion phenomenon, a point source liquid feeding configuration was considered. Hold‐up maps obtained from photonic flux attenuation measurements were then used to characterize liquid dispersion from a source point for both packings (Mellapak 250.X and IMTP. 40) using a relatively simple advection‐diffusion model.

Dû aux efforts apportés à la réduction des émissions de NOx dans des chambres de combustion aéronautiques il y a une tendance récente vers des systèmes à combustion pauvre. Cela résulte dans l'apparition de nouveaux types d'injecteur qui sont caractérisés par une complexité géométrique accrue et par des nouvelles stratégies pour l'injection du carburant liquide, comme des systèmes multi-point. Les deux éléments créent des exigences supplémentaires pour des outils de simulation numériques. La simulation à grandes échelles (SGE ou LES en anglais) est aujourd’hui considérée comme la méthode la plus prometteuse pour capturer les phénomènes d'écoulement complexes qui apparaissent dans une telle application. Dans le présent travail, deux sujets principaux sont abordés : Le premier est le traitement de la paroi ce qui nécessite une modélisation qui reste délicate en SGE, en particulier dans des géométries complexes. Une nouvelle méthode d'implementation pour des lois de paroi est proposée. Une étude dans une géométrie réaliste démontre que la nouvelle formulation donne de meilleurs résultats comparé à l’implémentation classique. Ensuite, la capacité d'une approche SGE typique (utilisant des lois de paroi) de prédire la perte de charge dans une géométrie représentative est analysée et des sources d'erreur sont identifiés. Le deuxième sujet est la simulation du carburant liquide dans une chambre de combustion. Avec des méthodes Eulériennes et Lagrangiennes, deux approches sont disponibles pour cette tâche. La méthode Eulérienne considère un spray de gouttelettes comme un milieu continu pour lequel on peut écrire des équations de transport. Dans la formulation Lagrangienne, des gouttes individuelles sont suivies ce qui mène à des équations simples. D’autre part, sur le plan numérique, le grand nombre de gouttes à traiter peut s’avérer délicat. La comparaison des deux méthodes sous conditions identiques (solveur gazeux, modèles physiques) est un aspect central du présent travail. Les phénomènes les plus importants dans ce contexte sont l'évaporation ainsi que le problème d'injection d'un jet liquide dans un écoulement gazeux transverse ce qui correspond à une version simplifiée d’un système multi-point. Le cas d'application final est la configuration d’un seul injecteur aéronautique, monté dans un banc d'essai expérimental. Ceci permet d'appliquer de manière simultanée tous les développements préliminaires de ce travail. L'écoulement considéré est non-réactif mais à part cela il correspond au régime ralenti d'un moteur d'avion. Dû aux conditions préchauffées, le spray issu du système d'injection multi-point s'évapore dans la chambre. Cet écoulement est simulé utilisant les approches Eulériennes et Lagrangiennes et les résultats sont comparés aux données expérimentales. / Due to efforts to reduce NOx emissions of aeronautical combustors, there is a recent trend towards lean combustion technologies. This results in novel injector designs, which are characterized by increased geometrical complexity and new injection strategies for the liquid fuel, such as multipoint systems. Both elements create additional challenges for numerical simulation tools. Large-Eddy simulation (LES) is regarded as the most promising method to capture complex flow phenomena in such an application. In the present work, two main areas of interest are considered: The first is wall modeling, which remains a challenging field in LES, in particular for complex geometries. A new implementation method for wall functions that uses a no-slip condition at the wall is proposed. It is shown that in a realistic burner geometry the new formulation yields improved results compared to a classical implementation. Furthermore, the capability of a typical LES with wall models to predict the pressure drop in a representative geometry is assessed and sources of error are identified. The second topic is the simulation of liquid fuel in a combustor. With Eulerian and Lagrangian methods, two different approaches are available for this task. The Eulerian approach considers a droplet spray as a continuum for which transport equations can be formulated. In the Lagrangian formulation, individual droplets are tracked, which leads to a simple formulation but can be challenging in terms of numerics due to the large number of particles to be treated. The comparison of these methods under identical conditions (gaseous flow solver, physical models) is a central aspect of the present work. The most important phenomena that are studied in view of the final application are evaporation and the problem of transverse liquid jets in a gaseous crossflow as a simplified representation of a multipoint system. The final application case is the configuration of a single aeronautical injector mounted in an experimental test bench. It allows to simultaneously apply all preliminary developments. The flow considered is non-reactive but otherwise corresponds to a partial load regime in an aeroengine Due to the pre-heated conditions, the spray issued by the multi-point injection undergoes evaporation. This flow is simulated using Eulerian and Lagrangian methods and the results are compared to experimental data.

De nos jours, la combustion représente encore un 90% de la production totale d'énergie au monde. La plupart des brûleurs de type industriel utilisent comme carburant des hydrocarbures en forme liquide. Cependant, un grand nombre d'études ont été dédiés aux flammes gazeuses et l'impact du spray liquide est encore loin d'être totalement compris. Le but de cet étude est l'amélioration de la modélisation des deux phénomènes principaux qui ont lieu entre l'atomisation du spray et la combustion, i.e. la dispersion des gouttes par la turbulence gazeuse et le procès d'évaporation dans le contexte de la Simulation Aux Grandes Echelles (SGE) des configurations complexes. Premièrement, l'approche Euler-Euler mésoscopique (Février et al. (2005)), basée sur une moyenne d'ensemble conditionnée et implémentée dans AVBP est amélioré. Le modèle de fermeture (Simonin et al. (2001); Kaufmann (2004)) pour les moments de deuxième ordre qui apparait dans les équations de transport résolues échoue quand appliqué à des configurations cisaillées (Riber (2007)). Plusieurs modèles proposés récemment par Masi (2010) et qui ont été valides a priori dans une configuration de nappe chargée de particules sont validés a posteriori dans la même configuration. Un analyse quantitative sur plusieurs cas avec diffèrent nombres de Stokes, nombres de Reynolds de la phase gazeuse et résolutions du maillage ont permit de retenir un modèle non-linéaire nommé 2EASM3, qui utilise le tenseur de déformations de la phase dispersée comme échelle de temps caractéristique. La deuxième partie a pour but l'amélioration du modèle d'évaporation implémenté dans AVBP. Ce modèle suppose une conduction infinie dans la phase liquide et symétrie sphérique dans la phase gazeuse ainsi que des lois simplifiées pour les propriétés thermodynamiques et de transport. Un nouveau modèle prenant en compte la dépendance de la viscosité du mélange gazeux avec la composition locale, et des nombres de Prandtl et Schmidt fixés par les valeurs à l'équilibre obtenus par moyen d'une simulation prenant en compte des lois complèxes pour les propriétés thermodynamiques et de transport est proposé. Cette nouvelle méthode produit des résultats en bon accord avec les mesures expérimentales pour l'évaporation d'une goutte isolé en une atmosphère d'azote au calme sans pourtant augmenter le cout du calcul. Finalement, l'impacte des nouveaux modèles est analysé dans une SGE de la configuration semi-industrielle MERCATO (García-Rosa (2008)). Bien que les données expérimentales ne soient pas suffisantes pour confirmer les résultats, les distributions de gouttes et de carburant gazeux sont significativement affectés par les modèles, ce qui pourrait avoir un impact directe sur le procès d'allumage. / Combustion still represents about 90% of the energy production in the world. Most industrial burners are fuelled with liquid hydrocarbons. However, most studies have been dedicated to gaseous ßames and the impact of liquid spray is still misunderstood. The purpose of this study is to improve the modelisation of two main phenomena occurring between atomization and combustion, i.e. the droplet dispersion in the turbulent gaseous flow and the evaporation process, in the context of Large Eddy Simulation (LES) of complex configurations. First, the mesoscopic Euler-Euler approach (Février et al. (2005)) based on a conditioned ensemble averaging and implemented in AVBP is improved. The closure model (Simonin et al. (2001), Kaufmann (2004)) for the second-order moments appearing in the transport equations solved fails in mean-sheared configurations (Riber (2007)). Several new models proposed by Masi (2010) and a priori tested in a particle-laden slab are tested a posteriori in the same configuration. A quantitative analysis based on several calculations varying the Stokes number, the gaseous Reynolds number and the grid resolution allows to retain a non-linear model using the particle rate-of-strain tensor as timescale and called 2EASM3. The second part consists in improving the evaporation model implemented in AVBP which assumes infinite conduction in the liquid and spherical symmetry in the gas phase along with simplified thermodynamics and transport properties calculation. A new model is proposed, where the dependence of gaseous mixture viscosity on local composition is accounted for, and the Prandtl and Schmidt numbers are fixed by a reference equilibrium calculation using complex thermodynamics and transport properties. This method shows good agreement with experimental measurements in the configuration of an isolated droplet evaporating in quiescent N2 without further increasing the computational cost. Finally, the impact of the new models is analysed in the LES of the MERCATO semi-industrial configuration (García-Rosa (2008)). Although the experimental data are not sufficient to confirm the results, both the droplet distribution and the fuel mass fraction are significantly affected, which would eventually affect the ignition process.

La régulation des polluants a mené à la création de nouveaux systèmes de combustion. Le carburant étant stocké sous forme liquide, sa transformation jusqu’à sa combustion est complexe. La capacité de la Simulation aux grandes échelles à simuler des écoulements turbulents réactifs a été montrée sur des cas académiques comme sur des configurations industrielles, tout en prenant en compte les phénomènes multiphysiques intervenant dans ces configurations, mais les études sur la structure de flamme diphasique sont encore trop peu nombreuses. La présence de deux solveurs pour la simulation d’une phase liquide étant disponible dans le code AVBP, leur utilisation permet une comparaison et une compréhension des phénomènes en jeu combinant dispersion, évaporation, et combustion. La première partie de l’étude relate la validation du modèle d’injection FIM-UR. Ce modèle est capable de reconstruire les profils de vitesses et de granulométrie à l’injecteur sans avoir à simuler les phénomènes d’atomisation primaire et secondaire. Une validation en régime turbulent avait déjà été réalisée, et on propose ici de valider le modèle dans un cas laminaire. Des comparaisons entre simulations monodisperses et polydisperse et des expériences sont effectuées. La simulation monodisperse Lagrangienne donne une bonne structure globale mais la simulation polydisperse Lagrangienne permet de retrouver le comportement au centre du cône avec la présence des petites gouttes et à la périphérie du cône par la présence des grosses gouttes. De plus, des améliorations sont apportées au modèle pour le formalisme Eulérien et montrent de bons résultats. La partie suivante s’intéresse à caractériser un spray polydisperse par une distribution monodisperse. En effet, au cas où une approche polydisperse n’est pas possible, le choix du diamètre moyen à prendre pour une simulation monodisperse est délicat. On propose donc d’analyser le comportement d’un spray polydisperse en le comparant à ceux de sprays monodisperses. Deux configurations académiques sont choisies : des cas de Turbulence Homogène Isotrope chargée en particules pour étudier la dynamique, et des calculs d’évaporation 0D. Trois paramètres sont étudiés pour la dynamique : la concentration préférentielle (ou ségrégation), la traînée moyenne et la traînée réduite moyenne. Cette dernière et la ségrégation de la distribution polydisperse semblent affectées par les tailles de goutte les plus faibles, et la concentration préférentielle apparait alors comme la moyenne des ségrégations des classes qui la composent pondérées par l’inverse du nombre de Stokes associé à chacune de ces classes. La traînée moyenne de la simulation polydisperse possède un comportement proche des diamètres moyens D10 et D20. Ces analyses nous poussent donc à choisir le D10 pour caractériser la dynamique d’un spray polydisperse. Les calculs d’évaporation 0D ne permettent pas dans un premier temps de caractériser efficacement la masse évaporée d’un spray polydisperse par celle d’un spray monodisperse équivalent, mais la définition de nouveaux diamètres issus de la littérature des lits fluidisés comme le D50% le permet, ce qui le place autour du D32. On propose donc de caractériser l’évaporation d’un spray polydisperse par ce diamètre. Enfin, la dernière partie étudie la structure de flamme diphasique dans la chambre MERCATO, à l’aide du formalisme Lagrangien, monodisperse et polydisperse, mais aussi en utilisant le formalisme Eulérien. La validation du modèle FIM-UR du premier chapitre et ses améliorations sont utilisées pour représenter les conditions d’injection liquide. En plus d’un calcul polydisperse, deux simulations monodisperses Lagrangiennes sont réalisées en prenant les diamètres moyens D10 et D32, suite à la partie précédente. Des comparaisons qualitatives et des validations sont réalisées, en comparant des profils de vitesses gazeuses axiale et fluctuante et vitesse axiale liquide issus de l’expérience. / Regulations on pollutants have led to the creation of new combustion systems. Giving that fuel is stored in a liquid form, its evolution until combustion is complex. The ability of Large Eddy Simulation has been demonstrated on academic cases, as well as on industrial configurations, by taking into account the multi-physics phenomena, but there is a lack of studies about two-phase flow flame structures. Two solvers for the simulation of two-phase flows are available in the AVBP code, hence both simulations are performed to compare and increase understanding of the phenomena involved such as dispersion, evaporation and combustion. The first part of the study focuses on the validation of the FIM-UR injection model. This model is able to build velocity and droplet profiles at the injector, without simulating primary and secondary break up. A validation in a turbulent case has already been done, and this study validates the model in a laminar case. Comparisons between monodisperse and polydisperse simulations, and experiments are performed. The monodisperse Lagrangian simulation shows good results but the polydisperse simulation is able to represent profiles in the center of the cone by small droplets and at the peripheral part of the cone, by big ones. Moreover, improvements in the Eulerian model exhibit good results. The next section tries to evaluate the impact of polydispersion. Indeed, when a polydisperse approach is not available, choosing the mean diameter can be tricky. A comparison between the behavior of polydisperse spray and monodisperse sprays ones is realised. Two academic cases are studied: Homogeneous Isotropic Turbulence with particles to analyze the dynamics, and 0D evaporation cases. For the dynamics, preferential concentration, mean drag and reduced mean drag are studied. The latter and preferential concentration are affected by small droplets, and the preferential concentration of a polydisperse spray is equivalent to the average of preferential concentration of classes, extracted from the polydisperse distribution, weighted by the inverse of the Stokes number of each class. The mean drag behaves like the D10 and D20 mean drags. This analysis allows us to choose the D10 to characterize a polydisperse distribution for the dynamics. Zero-D evaporation simulations cannot characterize the polydisperse spray evaporated mass by the evaporated mass of monodisperses sprays. New definitions of diameters from fluidized bed literature enable the use of D50%, which is close to D32. We propose to use this diameter to characterize the evaporation of a polydisperse spray. Finally, the last section studies the structure of two-phase flames in the MERCATO bench, using the Lagrangian formalism, monodisperse and polydisperse but also using the Eulerian formalism. The validation of FIM-UR model and improvements from the first section are used to represent liquid injection conditions. A polydisperse simulation is realized and two monodisperse simulations are computed using mean diameters D10 and D32, thanks to the previous section. Qualitative comparisons and validations are realized, comparing gaseous velocity profiles and liquid velocity profiles. Good agreements are found and the mean diameter D32 seems to be close to the polydisperse spray. A comparison between mean flames is done with an Abel transform of the flame from the experiments. The flame has an "M shape", anchored by small recirculation zones out of the swirler, and by a point at the tip of the central recirculation zone. Then, the impact of droplet distributions is analyzed. Even if few bigger droplets from the polydisperse distribution are convected in the hot gases due to bigger particular time and evaporation time, two-phase flow flame structures are equivalent. Different combustion regimes appeared with premixed flames and pockets of fuel burning in the hot gases.

An inconsistency exists in bubble force models used in the standard Euler-Euler simulations. The bubble force models are typically developed by assuming that the forces act on the bubbles' centers of mass. However, in the standard Euler-Euler model, each bubble force is a function of the local gas volume fraction because the phase-averaging method is used. This inconsistency can lead to gas over-concentration in the center or near the wall of a channel when the bubble diameter is larger than the computational cell size. Besides, a mesh-independent solution may not exist in such cases. In addition, the bubble deformation is not fully considered in the standard Euler-Euler model. In this thesis, a particle-center-averaging method is used to represent the bubble forces as forces that act on the bubbles' centers of mass. A particle-center-averaged Euler-Euler approach for bubbly flow simulations is developed by combining the particle-center-averaged Euler-Euler framework with a Gaussian convolution method. The convolution method is used to convert the phase-averaged and the particle-center-averaged quantities. The remediation of the inconsistency in the standard Euler-Euler model by the particle-center-averaging method is demonstrated using a simplified two-dimensional test case. Bubbly flows in different vertical pipes are used to validate the particle-center-averaged Euler-Euler approach. The bubbly flow simulation results for the particle-center-averaged Euler-Euler model and the standard Euler-Euler model are compared with experimental data. For monodisperse simulations, the particle-center-averaging method alleviates the over-predictions of the gas volume fraction peaks for wall-peaking cases and for finely dispersed flow case. Whereas, no improvement is found in the simulated gas volume fraction for center-peaking cases because the over-prediction caused by the inconsistency has been smoothed by the turbulent dispersion. Moreover, the axial gas and liquid velocities simulated with both Euler-Euler models are similar, which proves that the closure models for bubble forces and turbulence are correctly applied in the particle-center-averaged Euler-Euler model. For fixed polydisperse simulations, the particle-center-averaging method can also alleviate the over-prediction of the gas volume fraction peak in the center or near the wall of a pipe. The axial gas velocities simulated with both Euler-Euler models are about the same. Comparisons are also made for the simulation results of bubbly flows in a cylindrical bubble column and the experimental data. The gas volume fractions and the axial gas velocities simulated with both Euler-Euler models almost coincide with each other, which indicates that the sink and source terms for the continuity equations and the degassing boundary are set correctly in the particle-center-averaged Euler-Euler model. An oblate ellipsoidal bubble shape is considered in the particle-center-averaged Euler-Euler simulations by an anisotropic diffusion. The influence of bubble shape on the simulation results of bubbly pipe flows is investigated. The results show that considering the oblate ellipsoidal bubble shape in simulations can further alleviate the over-predictions of the gas volume fraction peaks for wall peaking cases, but it has little influence on the gas volume fractions of center-peaking cases and the axial gas velocities.

The contribution compares results of Large Eddy Simulations of the cavitating flow in a model oil hydraulic spool valve using an Euler-Euler and a one-way coupled Euler- Lagrange model. The impact of the choice of the empirical constants in the Kunz cavitation model is demonstrated. Provided these are chosen appropriately the approach can yield reasonable agreement with the corresponding experiment. The one-way Euler-Lagrange model yields less agreement. It is demonstrated that this is due to the lack of realistic volumetric coupling, rarely accounted for in this type of method. First results of such an algorithm are presented featuring substantially more realism.

Les turbines aéronautiques doivent satisfaire à des normes d'émissions polluantes toujours en baisse. La qualité du mélange du carburant et de l'air dans la chambre de combustion est responsable de la formation de polluants nocifs pour l'environnement. La simulation aux grandes échelles (LES) permet d'étudier les mécanismes de mélanges turbulents de l'air et du carburant. La prise en compte de l'aspect liquide du carburant injecté devient nécessaire pour prédire correctement l'apparition de vapeur de carburant au sein du foyer. Le but de cette thèse est évaluer la fiabilité des simulations LES Euler-Euler dans une configuration complexe. Les processus d'injection, et d'évaporation du carburant liquide sont analysés et modélisés dans les simulations LES car ils pilotent la formation de vapeur de carburant. Les méthodes numériques pour résoudre les équations continues de la phase dispersée doivent permettre des simulations précises et robustes dans une configuration représentative d'une chambre de combustion. Les simulations présentées dans ces travaux reproduisent l'écoulement diphasique évaporant non-réactif du banc d'essai Mercato. Ce banc est équipé d'un système d'injection d'air vrillé et d'un atomiseur pressurisé-swirlé de kérosène typiques des foyers aéronautiques réels. Dans ces travaux, le modèle pour l'injection de liquide FIM-UR a été développé pour définir les conditions limites conduisant à un spray issu d'un atomiseur préssurisé-swirlé. Le kérosène employé dans les campagnes expérimentales est modélisé dans les simulations par un composé permettant d'obtenir des temps d'évaporation réalistes. Trois stratégies numériques ont été mises en place sur la configuration MERCATO. Les comparaisons des résultats numériques aux mesures expérimentales ont permis d'évaluer la stratégie numérique conduisant à la meilleure précision. L'utilisation du schéma centré TTGC associé à un opérateur de viscosité artificielle localisée par un senseur adapté est optimale lorsque l'équation sur l'énergie décorrélée des gouttes est résolue. Cette stratégie permet de contrôler la localisation et les niveaux de viscosité par rapport à un schéma décentré. Les termes sources liés au mouvement mésoscopique permettent de redistribuer l'énergie dans les zones de compression ou de détente de la phase dispersée, et d'obtenir les bonnes répartitions des fluctuations dans la chambre de combustion. La stratégie retenue est comparée aux statistiques de la dynamique du spray résolu par une approche Lagrangienne employant la même injection monodispersse. Le méthode Euler-Euler conduit à la même précision de la dynamique de la phase dispersée que la méthode Euler-Lagrange. L'accès à l'évolution instationnaire de l'écoulement permet d'identifier les mêmes mécanismes de dispersion et de mélange dans les deux simulations. Des différences sur la répartition de diamètre moyen et de carburant dans la chambre ont été mis en évidence et reliés à la polydispersion locale qui n'est pas résolue dans l'approche Euler-Euler monodisperse et qui apparaît naturellement dans l'approche Euler-Lagrange malgré l'injection monodisperse.