Les écoulements diphasiques en milieu poreux sont des phénomènes complexes et qui concernent de nombreux problèmes industriels. EDF travaille sur la faisabilité et la sécurité d'un stockage en couche géologique profonde de déchets nucléaires. Dans ce domaine la simulation des écoulements diphasiques en milieu poreux est particulièrement importante dans au moins trois domaines : tout d'abord lors de la phase de ventilation des galeries du stockage qui pourrait désaturer la roche présente et ainsi en modifier ses propriétés de rétention, mais également lors de la phase de resaturation des matériaux et enfin lors de l'arrivée de l'eau sur les parties métalliques contenues dans le stockage qui entraînera alors des phénomènes de corrosion et un dégagement d'hydrogène. Dans ce contexte, EDF souhaite se doter de méthodes numériques performantes et robustes ne nécessitant pas de conditions restrictives sur la forme des mailles. Ce travail s'inscrivant dans cette problématique, est consacré dans un premier temps au développement du schéma volumes finis SUSHI (Scheme Using Stabilization and Hybrid Interfaces) dans le code de mécanique d'EDF, Code_Aster afin de modéliser les écoulements diphasique en milieu poreux. Ce schéma a été développé en 2D et en 3D. Parallèlement une nouvelle formulation qui permet de traiter de manière uniforme les écoulements en milieu saturé et insaturé pour des problèmes miscibles et immiscibles est proposée. Différentes études modélisant des difficultés liées aux problématiques du stockage de déchets radioactifs en couches géologiques profondes ont été traitées. On peut citer l'étude d'un bi-matériau qui met en avant le ré-équilibrage capillaire d'un matériau par un autre possédant des propriétés et des conditions initiales en saturation très hétérogènes. On citera également l'étude de l'injection d'hydrogène dans un milieu initialement saturé en eau pure qui est tirée du benchmark " Ecoulement diphasique " proposé par le GNR MOMAS. Cette étude avait pour objectif de mettre en évidence le bon traitement de l'apparition d'une phase dans un milieu saturé et donc la pertinence de notre nouvelle formulation à traiter d'une manière unifié un problème d'écoulement saturé et un problème d'écoulement insaturé

Le travail de cette thèse aborde le sujet de l’influence des particules non-pesantes et de taille macroscopique sur les écoulements en canal dans des conditions proches du seuil de la transition laminaire-turbulent. Les suspensions sont faiblement concentrées (fraction volumique φ = 5%). Le couplage hydrodynamique existant entre la phase dispersée et la phase continue est résolu numériquement par la Force-Coupling Method, et les particules sont suivies d’une façon lagrangienne. Dans un écoulement laminaire de Couette ou de Poiseuille plan, nous montrons que les contraintes induites par la phase solide augmentent avec l’inertie, et l’influence de la concentration est plus faible qu’en régime de Stokes. Les particules avancent avec un retard dans la direction de l’écoulement et migrent à travers les lignes de courant (effet Segré-Silberberg en Poiseuille). Les vitesses de migration et de glissement s’amplifient avec l’inertie et sont du même ordre de grandeur quand Rep = O(1). Quand les particules sont lâchées librement dans un écoulement de Poiseuille plan en-deça du seuil critique de transition à la turbulence, la suspension initiale- ment homogène (φ = 5%) devient stratifiée, après un temps d’écoulement de plusieurs dizaines d’unités de temps (rapport de la hauteur du canal sur la vitesse moyenne de l’écoulement). Après une centaine d’unités de temps, nous observons le développement d’une instabilité à l’interface entre la zone chargée en particules et la zone de fluide pur. Des motifs dunaires prennent place dans la direction de la vorticité. Ces motifs sont soutenus par des écoulements secondaires d’intensités faibles mais non-nulles. Dans un écoulement au-dessus du seuil de transition, nous avons étudié les profils des phases continues et dispersées et réalisé des visualisations 3D afin de comprendre pourquoi les particules macroscopiques diminuent le nombre de Reynolds critique de relaminarisation de l’écoulement. Nous observons que les particules provoquent une augmentation significative des fluctuations de vitesses dans les directions transverses et qu’elles modifient les structures rotationnelles de l’écoulement, qui deviennent plus petites, plus nombreuses et plus énergétiques (plus grandes vitesses de rotation). Le coefficient de frottement pariétal de l’écoulement de suspension en régime de transition est supérieur à celui de l’écoulement monophasique. Quand le nombre de Reynolds est diminué et que l’écoulement devient finalement laminaire, le coefficient de frottement pariétal rejoint la loi laminaire d’un écoulement monophasique, à condition de substituer la viscosité effective de la suspension à la viscosité du fluide dans l’expression du nombre de Reynolds. D’après nos résultats, la turbulence de l’écoulement de suspension est conservée jusqu’à des nombres de Reynolds bien inférieurs à celui de l’écoulement monophasique en canal, en accord avec les observations ex- périmentales de Matas, Morris et Guazzelli (PRL, 2003) pour une géométrie cylindrique. Par ailleurs, nous montrons que selon le sens de la transition, laminaire → turbulent ou turbulent → laminaire, le nombre de Reynolds critique de transition d’un régime à l’autre n’est pas le même. Nous n’avons pas observé d’influence significative de la concentration en ce qui concerne la valeur du nombre de Reynolds critique de relaminarisation pour les deux concentrations étudiées (φ = 2.5% et 5%). / This PhD addresses the influence of macroscopic and neutrally buoyant particles on the channel flows close to the laminar-turbulent transition regime. The suspension flow is moderately concentrated (solid volumetric concentration _ = 5%). The hydrodynamic coupling between the dispersed and carrier fluid is numerically resolved using the Force-Coupling Method approach. Particle trajectories are obtained by lagrangian tracking. In laminar wall-bounded flows, we show that the stress induced by the solid phase increases with inertia, and that the effect of the concentration is weaker than in the Stokes regime. The particles lag the flow and they migrate across the streamlines (Segré-Silberberg effect in Poiseuille flow). The migration and slip velocities are of the same order of magnitude for Rep = O(1). When the particles are freely suspended in a Poiseuille flow below the transition threshold, the initially homogeneous suspension (_ = 5%) becomes stratified after several ten time units (channel height/average flow velocity). After a hundred time units, the different rheological properties of the suspension segregated parts induce an instability yielding the formation of dune-like patterns, sustained by weak but finite secondary flows. In the fluctuating flow regime, we studied the profiles of the continuous and dispersed phase and realized 3D visualizations in order to understand why finite size particles delay the relaminarization threshold. The particles induce a significant increase of the velocity fluctuations in the transverse directions and they modify the rotational flow structures, which become smaller, more numerous and more energetic (larger rotation velocity). The wall-friction coefficient of the suspension flow in the transition regime is larger than the single-phase flow case. When the Reynolds number is decreased and the flow becomes laminar, the friction coefficient recovers the laminar law of a single phase flow provided that the fluid viscosity is replaced by the effective suspension viscosity in the Reynolds number definition. Our results clearly show that the two-phase channel flow turbulence is conserved down to a threshold well below the single phase flow limit, in agreement with the observations of Matas, Morris et Guazzelli (PRL, 2003) for a cylindrical geometry. In addition, we show that according to the transition direction, i.e. laminar 7! turbulent or turbulent 7! laminar, the switch from a regime to another does not occur at the same critical Reynolds number. Finally, in the limit of moderately concentrated (_ = 2.5−5%) suspension flow in a channel, the concentration has no significant influence on the critical Reynolds number.

Les écoulements diphasiques dans les centrales de type réacteur à eau pressurisée appartiennent à des régimes très variés allant du faible nombre de Mach jusqu'aux ondes de chocs. Calculer des solutions approchées précises de ces écoulements peut s'avérer délicat dans certains régimes. On s'intéresse dans cette thèse à la conception et à l'étude de méthodes numériques robustes et stables à grand pas de temps, capables de calculer des solutions approchées précises quel que soit le régime d'écoulement, y compris sur maillage grossier. Une stratégie pour construire de tels schémas consiste à : utiliser un schéma semi implicite basé sur un splitting d’opérateurs pour séparer la résolution approchée des phénomènes rapides de celles des phénomènes lents ; corriger les flux numériques afin d’améliorer la précision du schéma dans certains régimes. Deux approches sont utilisées pour analyser la capacité du schéma numérique à gérer plusieurs régimes d'écoulement. L’approche des schémas asymptotic preserving est utilisée pour traiter le système de la dynamique des gaz avec termes sources raides. On utilise ensuite la notion de schéma tout-régime pour le système de la dynamique des gaz et les systèmes diphasiques homogénéisés HRM et HEM à bas nombre de Mach. Des propriétés garantissant la stabilité et la robustesse des schémas ont été obtenues, et en particulier des inégalités d'entropie discrètes. L'implémentation de ces méthodes a permis de mener des expériences numériques en 1D et 2D sur maillage non structuré qui confirment le gain en précision et en temps de calcul des schémas asymptotic preserving et tout-régime ainsi construits par rapport à des schémas numériques classiques. / Two-phase flows in Pressurized Water Reactors belong to a wide range of Mach number flows. Computing accurate approximate solutions of those flows may be challenging from a numerical point of view as classical finite volume methods are too diffusive in the low Mach regime. In this thesis, we are interested in designing and studying some robust numerical schemes that are stable for large time steps and accurate even on coarse meshes for a wide range of flow regimes. An important feature is the strategy to construct those schemes. We use a mixed implicit-explicit strategy based on an operator splitting to solve fast and slow phenomena separately. Then, we introduce a modification of a Suliciu type relaxation scheme to improve the accuracy of the numerical scheme in some regime of interest. Two approaches have been used to assess the ability of our numerical schemes to deal with a wide range of flow regimes. The first approach, based on the asymptotic preserving property, has been used for the gas dynamics equations with stiff source terms. The second approach, based on the all-regime property, has been used for the gas dynamics equations and the homogeneous two-phase flows models HRM and HEM in the low Mach regime. We obtained some robustness and stability properties for our numerical schemes. In particular, some discrete entropy inequalities are shown. Numerical evidences, in 1D and in 2D on unstructured meshes, assess the gain in term of accuracy and CPU time of those asymptotic preserving and all-regime numerical schemes in comparison with classical finite volume methods.

Actuellement, les codes de calcul de composants thermohydrauliques de réacteurs nucléaires du CEA et d'EDF utilisent des modèles physiques diphasiques de mélange à 3 ou 4 équations. Or, il existe un fort besoin industriel pour des modèles physiques plus sophistiqués tels que le modèle diphasique à 6 équations voire des modèles multichamps. Par ailleurs, le code système CATHARE du CEA, un des codes systèmes les plus utilisés aujourd'hui sur le plan international, utilise un modèle physique diphasique à 6 équations et un schéma numérique semi-implicite de type ICE à maillages décalés. Le schéma de CATHARE est connu pour sa robustesse dans une large gamme de configurations d'écoulement. En s'inspirant de l'expérience de CATHARE, on propose de mettre en œuvre un schéma volumes finis colocalisés de type « pressure based ». Le but est d'obtenir un décentrement des flux qui assure la robustesse du schéma tout en gardant une bonne précision. De plus, le fait de pouvoir utiliser des maillages colocalisés (structurés ou non-structurés) permet de traiter de différentes configurations complexes de cœur de réacteurs et de réaliser des calculs fin d'inter assemblage. Le schéma doit conserver exactement la masse et l'énergie et la solution numérique doit converger lorsque l'on raffine le maillage. Le schéma doit être capable de traiter des cas d'apparition et de disparition des phases, par exemple le cas de la colonne bouillante où il y a changement de phase dû au transfert de chaleur, des cas de tuyères avec changement de phase dû à un élargissement ou un rétrécissement brusque, ou de séparation de phase par gravité. En outre le schéma doit être capable de traiter des configurations de calcul à faible nombre de Mach, par exemple le cas du renoyage d'un cœur de réacteur. L'objet de la thèse consistera à développer une méthode Volumes Finis co-localisés (dans l'esprit de Ghidaglia et al.) et la direction d'investigation s'inspirera des travaux de Jeong et al. qui a conduit au code CUPID. Références. Ghidaglia, J. M., Kumbaro, A., & Le Coq, G. (2001). On the numerical solution to two fluid models via a cell centered finite volume method. European Journal of Mechanics-B/Fluids, 20(6), 841-867. Jeong, J. J., Yoon, H. Y., Cho, H. K., Kim, J., & Park, I. K. (2008). A semi-implicit numerical scheme for a transient two-fluid three-field model on an unstructured grid. International Communications in Heat and Mass Transfer, 35(5), 597-605. / Currently, there is a strong industrial need for sophisticated physical models such as the two-phase model with 6 equations or multi-field models for the thermo-hydraulic calculation of nuclear reactor components. In addition, the system code CATHARE of CEA, one of the most utilized system codes at the international level, employs a two-phase 6 equations and a semi-implicit numerical scheme of ICE type on staggered grids. The code CATHARE is known for its robustness in a wide range of flow configurations. Drawing on the experience of CATHARE, we propose to implement a collocated finite volume pressure based scheme. The aim is to obtain a decentering of flux which ensures the robustness of the scheme while keeping good accuracy. In addition, being able to use collocated grids (structured or unstructured) can handle different complex configurations. The scheme should conserve exactly the mass and energy and the numerical solution needs to converge when the mesh is refined. The scheme should be able to handle cases of phase appearance and disappearance, for example in the case of boiling column where phase change is due to heat transfer, the case of nozzles with phase change due to a widening or abrupt narrowing, or phase separation by gravity. In addition the scheme should be capable of calculation configurations at low Mach number. The purpose is to develop a co-located Finite Volume method (in the spirit of Ghidaglia et al.) And direction of investigation is build on work of Jeong et al. which led to CUPID code. References: Ghidaglia, J. M., Kumbaro, A., & Le Coq, G. (2001). On the numerical solution to two fluid models via a cell centered finite volume method. European Journal of Mechanics-B/Fluids, 20(6), 841-867. Jeong, J. J., Yoon, H. Y., Cho, H. K., Kim, J., & Park, I. K. (2008). A semi-implicit numerical scheme for a transient two-fluid three-field model on an unstructured grid. International Communications in Heat and Mass Transfer, 35(5), 597-605.

Ce travail de thèse est consacré à la modélisation et simulation numérique d'écoulements diphasiques liquide-gaz mettant en jeu des transferts de chaleur. La simulation de configurations où la prise en compte des effets de compressibilité de la phase gazeuse est indispensable (micropompes, microactionneurs, etc...) a nécessité l'utilisation d'un modèle original, considérant le liquide incompressible et le gaz compressible sous l'hypothèse faible Mach. Lors de cette thèse, ce modèle a été implémenté dans un code diphasique prenant en compte l'interface à l'aide d'une méthode de front-tracking. Des cas tests ont été développés spécifiquement afin de vérifier la conservation de l'énergie pour des configurations de complexité croissante. Les résultats des cas tests ont permis de mettre en évidence la difficulté à assurer la conservation de l'énergie lorsque l'interface n'est pas discontinue mais lissée, comme c'est le cas dans la méthode de front-tracking standard. Une méthode de traitement d'interface hybride a été proposée, rétablissant le caractère discontinu de l'interface avec la reconstruction d'une fonction indicatrice de phase échelon, tandis que le déplacement de l'interface est assuré d'un pas de temps à l'autre à l'aide du front-tracking. Les résultats obtenus avec cette nouvelle méthode hybride sont très satisfaisants, la méthode hybride permettant d'assurer la conservation de l'énergie et de la masse avec précision dans les simulations. / This thesis is devoted to the modeling and the numerical simulation of liquid-gas flows in non isothermal micro-cavities or micro-channels. The objective is to describe two-phase flows in which compressibility of the gaseous phase plays a key role (as for instance in micropumps, microactuators, etc...). An original model is developed, considering in the same computational domain, an incompressible liquid and a compressible gas under the low Mach approximation. This model has been implemented in a code using the front-tracking method for the interface description. In order to check the proper satisfaction of the energy balance, specific test cases have been developed considering several configurations of increasing complexity. It has been shown from these test cases that energy conservation can hardly be satisfied when the interface is described by the means of a smooth function, which is done in the front-tracking method. An hybrid method has been proposed, restoring the discontinous nature of the interface. It makes use of a step function combined with the front tracking method. Results obtained with this new hybrid method show that mass conservation and energy balance are very properly enforced during the computations. This thesis is devoted to the modelling and numerical simulation of liqui-gas flow envolving heat transfer. Simulation of configurations where it is essential to take into account the compressibility nature of the gaseous phase(for instance micropumps, microactuators, etc...) require the use of an original model, considering an incompressible liquid and a compressible gas under the low Mach assuption. During this thesis, the model has been implemented in a multiphase flow code using the front-tracking method to handle the interface. Test cases have been developped specifically to check the energy conservation for differents configurations of an increasing complexity. Numerical results highlighted difficulties encountered to ensure the energy conservation while using smooth description of the interface, as it is the case in the standard front-tracking method. An hybrid method has been proposed, restoring the discontinous character of the interface by reconstructing a step maker function, whereas the front displacement from a time step to the next, is still handled with the front-tracking. Results obtained using this new hybride method are very satisfactory, the hybrid method allowing the code to ensure accurately the energy and mass conservation during the computations.

Cette étude concerne l’étude des écoulements diphasiques dans le cadre du refroidissement des composants électroniques en systèmes embarqués. L’étude bibliographique a permis de sélectionner le refroidissement par spray comme technique prometteuse pour dissiper des flux de chaleur au-delà de 100 W/cm2. Une étude hydraulique, utilisant de l’eau et du HFE7100 comme fluides de refroidissement, nous a permis de valider des modèles permettant de déterminer la taille et la vitesse de gouttes provenant d’une sélection de buses de spray. Pour la partie thermique, nous avons conçu une section d’essais (évaporateur) permettant de pulvériser en spray afin d’étudier le refroidissement avec une boucle fermée diphasique.Vu la complexité du système de spray influencé par plusieurs paramètres et phénomènes physiques, nous avons isolé le phénomène d’ébullition nucléée dans une configuration en ébullition nucléée avec un élément chauffant identique à celui employé avec le refroidissement par spray. Pour améliorer les échanges thermiques, 6 surfaces avec différentes structurations (macroscopiques, microscopiques et hybrides) ont été sélectionnées. Les résultats de tests avec ces surfaces ont été comparés avec une surface lisse tant pour le refroidissement par spray que pour le refroidissement en vase.D’une part, avec un refroidissement par spray, les surfaces macrostructurées nous ont permis de dissiper des puissances thermiques de l’ordre de 140 W/cm2 avec d’importants coefficients de transfert thermique. D’autre part, avec un système de refroidissement par immersion, une des surfaces hybrides a montré être la plus performante.Les résultats reportés dans cette thèse ont permis d’approfondir la compréhension des mécanismes de transfert de chaleur en refroidissement par spray. De même, ils ouvrent la voie à l’étude des améliorations et optimisations du système permettant de l’employer en systèmes embarqués. / This dissertation concerns the study of two-phase flow cooling of electronic components in embedded systems. From a literature review, Spray Cooling was selected as a promising technique for dissipating heat fluxes above 100 W/cm2. A hydraulic study, using water and HFE7100 as coolants, has validate models for determining the size and speed of drops from a selection of spray nozzles. Regarding the thermal study, we have designed a test section (evaporator) to study cooling in a two-phase closed loop system.Given the complexity of Spray Cooling systems, which are influenced by several parameters and involve several physical phenomena, the nucleate boiling phenomenon has been isolated in a Pool Boling system with an identical heating element as Spray Cooling experiment. To improve heat exchange, 6 surfaces with different structures (macroscopic, microscopic and hybrid) were selected. The boiling test results with these surfaces have been compared with a smooth surface for both Spray Cooling and Pool Boiling.On one hand, in Spray Cooling tests, the macrostructured surfaces dissipated heat flux up to 140 W/cm2 with significant heat transfer coefficients. On the other hand, in the Pool Boling system, one of the hybrid surfaces has shown to be the most efficient.The results reported in this dissertation contributes on the understanding of the boiling mechanisms of heat transfer in Spray Cooling. Likewise, they open the way to the study of improvements and optimizations of the system for its use in embedded systems.

Lors de sa rentrée dans l’atmosphère d’une planète, un engin spatial subit un échauffement important dû aux frottements des gaz atmosphériques sur la paroi. Cette élévation de température conduit à une dégradation physico-chimique du bouclier thermique de l’objet constitué de matériaux composites. Un composite est constitué de divers matériaux qui s’ablatent différemment. Dans cette thèse, nous nous intéressons essentiellement à la fusion d’un matériau durant sa phase de rentrée atmosphérique. Nous sommes donc en présence de trois phases : solide, liquide et gaz. Pour simuler ce phénomène, des méthodes numériques robustes ont été mises au point pour calculer l’écoulement diphasique compressible autour de l’objet. Le couplage entre le solide et l’écoulement fluide a aussi été étudié. Les méthodes numériques développées durant cette thèse sont basées sur une approche volumes finis. Une stratégie de décomposition d’opérateurs est utilisée pour résoudre le modèle diphasique à cinq équations avec les termes de dissipation modélisant l’écoulement fluide. L’idée principale de cette décomposition d’opérateurs est de séparer les phénomènes acoustiques et dissipatifs des phénomènes de transport. Un traitement implicite de l’étape acoustique est réalisé tandis que l’étape de transport est résolue explicitement. Le schéma semi-implicite global est alors très robuste, conservatif et préserve les discontinuités de contact. Les conditions d’interface entre les domaines fluide et solide sont déduites des bilans de masse et d’énergie à la paroi. Le front de fusion est suivi explicitement grâce à une formulation ALE des équations. La robustesse de l’approche et l’apport de la formulation semi-implicite sont finalement démontrés grâce à des expériences numériques mono et bidimensionnelles sur maillages curvilignes mobiles. / During atmospheric re-entry phase, a spacecraft undergoes a sudden increase of the temperature due to the friction of atmospheric gases. This rise drives to a physical-chemical degradation of the thermal protective system of the object made of composite material. A composite is made of several materials with ablates differently. In this thesis, we mainly focus on the melting of an object during its re-entry phase. Therefore there are three phases: solid, liquid and gas phases. In order to simulate this phenomenon, robust numerical methods have been developed to compute a compressible multiphase flow. The coupling strategy between the solid and the fluid have also been studied. Solvers developed in the present work are based on Finite Volume Method. A splitting strategy is used to compute compressible two-phase flows using the five-equation model with viscous and heat conduction effects. The main idea of the splitting is to separate the acoustic and dissipative phenomena from the transport one. An implicit treatment of the acoustic step is performed while the transport step is solved explicitly. The overall scheme resulting from this splitting operator strategy is very robust, conservative, and preserves contact discontinuities. The boundary interface condition between the solid and the multiphase flow is enforced by mass and energy balances at the wall. The melting front is tracked explicitly using an ALE formulation of the equations. The robustness of the approach and the interest of the semi-implicit formulation are demonstrated through numerical simulations in one and two dimensions on moving curvilinear grids.

L'allumage d'un moteur fusée cryotechnique (carburants liquides) est une phase critique. La moindre anomalie dans la procédure d'allumage peut conduire à la destruction du lanceur. L'objectif de cette thèse est de développer une méthodologie s'appuyant sur la simulation aux grandes échelles (LES) pour étudier les phénomènes physiques impliqués dans un tel allumage. L'intérêt de la méthode LES est de pouvoir capturer les couplages instationnaires entre la turbulence, les processus diphasiques et la cinétique chimique. L'outil numérique est tout d'abord validé sur des cas académiques et expérimentaux, puis appliqué à un moteur fusée réel. Une approche graduelle est employée : les différents cas de validation présentent une complexité croissante, permettant d'isoler les processus physiques principaux. Ce travail de recherche montre que l'approche de la simulation aux grandes échelles, dans un contexte de calcul massivement parallèle, peut être utilisée pour étudier la séquence complète d'allumage dans un moteur fusée réel.

De nombreuses applications industrielles mettent en jeu des écoulements gaz-particules, comme les turbines aéronautiques et les réacteurs à lit fluidifié de l'industrie chimique. La prédiction des propriétés de la phase dispersée, est essentielle à l'amélioration et la conception des dispositifs conformément aux nouvelles normes européennes des émissions polluantes. L'objectif de cette thèse est de développer le formalisme Euler-Lagrange dans un solveur parallèle et non-structure pour la simulation aux grandes échelles pour ce type d'écoulements. Ce travail est motivé par l'augmentation rapide de la puissance de calcul des machines massivement parallèles qui ouvre une nouvelle voie pour des simulations qui étaient prohibitives il y a une décennie. Une attention particulière a été portée aux structures de données afin de conserver une certaine simplicité et la portabilité du code sur des différentes architectures. Les développements sont valides pour deux configurations : un cas académique de turbulence homogène isotrope décroissante et un calcul polydisperse d'un jet recirculant charge en particules. L'équilibrage de charges de particules est mis en évidence comme une solution prometteuse pour les simulations diphasiques Lagrangiennes afin d'améliorer les performances des calculs lorsque le deséquilibrage est trop important.

La pulvérisation du liquide est un mécanisme très important lors de la combustion dans les moteurs de fusée et de type diesel. C'est dans ce cadre que s'inscrivent les travaux réalisés dans cette thèse. L'objectif est de modéliser les mécanismes physiques qui pilotent l'atomisation assistée par air. Dans cette configuration, de nombreux paramètres peuvent influer sur la fragmentation du liquide tels que les vitesses initiales du liquide ou du gaz, les masses volumiques des deux phases ou la tension de surface par exemple ... Ce phénomène est donc particulièrement complexe et la résolution exacte de chaque fragmentation est très difficile. Dans ce travail, une approche stochastique est proposée afin de modéliser les différentes étapes d'atomisation assistée par air. Celle-ci permet de décrire la phase d'atomisation primaire et les statistiques du liquide non fragmenté au voisinage de l'injecteur. A partir de ces résultats, des fragments liquides sont formés et transportés dans l'écoulement de façon à décrire l'intégralité du spray. Le mouvement de la phase gazeuse est résolu par simulation des grandes échelles. Le caractère turbulent et instationnaire de l'écoulement est ainsi pris en compte, ce qui permet le suivi lagrangien des éléments liquides issus de la phase d'atomisation primaire. Leurs collisions ainsi que les fragmentations ou coalescences qui en résultent sont également modélisés par l'intermédiaire d'une approche stochastique. Ces différentes étapes permettent de simuler l'atomisation assistée par air et de déterminer la longueur du cœur liquide, l'angle du spray, le diamètre de Sauter moyen et de comparer ces résultats à des données expérimentales.