L'allumage et le réallumage de haute altitude présentent de grandes difficultés dans le cadre des chambres de combustion aéronautiques. Le succès d'un allumage dépend de multiples facteurs, des caractéristiques de l'allumeur à la taille des gouttes du spray en passant par le niveau de turbulence au point d'allumage. Déterminer la position optimale de l'allumeur ou le potentiel d'allumage d'une source d'énergie donnée à une position donnée sont ainsi des paramètres essentiels lors du design de chambre de combustion. Le but de ces travaux de thèse est d'étudier l'allumage forcé des chambres de combustion aéronautiques. Pour cela, des Simulation numériques aux Grandes Echelles (SGE) d'écoulements diphasiques réactifs sont utilisées et analysées. Afin de les valider, des données expérimentales issues du banc MERCATO installé à l'ONERA Fauga-Mauzac sont utilisées. Cela permet dans un premier temps de valider la méthodologie ainsi que les modèles utilisés pour les SGE diphasiques évaporantes avant leur utilisation dans d'autres conditions d'écoulement. Le cas diphasique réactif statistiquement stationnaire est ensuite comparé aux données disponibles pour évaluer les modèles en condition réactives. Ce cas est étudié plus en détail à travers l'analyse de caractéristiques de la flamme. Celle-ci semble être le théâtre de régimes de combustion très différents. On note aussi que la détermination de la méthode numérique la plus appropriée pour le calcul d'écoulements diphasiques n'est pas évidente. De plus, deux méthodes numériques différentes peuvent donner des résultats en bon accord avec l'expérience et pourtant avoir des modes de combustion différents. Les capacités de la SGE à correctement calculer un écoulement diphasique réactif étant validé, des SGE du phénomène transitoire d'allumage sont effectuées. La sensibilité observée expérimentalement de l'allumage aux conditions initiales, i.e. à l'instant de claquage, est retrouvé par les SGE. L'analyse met en évidence le rôle prépondérant de la dispersion du spray dans le développement initial du noyau de flamme. L'utilisation des SGE pour calculer les séquences d'allumage fournie de nombreuses informations sur le phénomène d'allumage, cependant d'un point de vue industriel, cela ne donne pas de résultat optimal, à moins de ne tester toutes les positions, ce qui rendrait le coût CPU déraisonnable. Des alternatives sont donc nécessaires et font l'objet de la dernière partie de ces travaux. On propose de dériver un critère local d'allumage, donnant la probabilité d'allumage à partir d'un écoulement diphasique (air et carburant) non réactif instationnaire. Ce modèle est basé sur des critères liés aux différentes phases menant à un allumage réussi, de la formation d'un premier noyau à la propagation de la flamme vers l'injecteur. Enfin, des comparaisons avec des données expérimentales sur des chambres aéronautiques sont présentées et sont en bon accord, indiquant que le critère d'allumage proposé, couplé avec une SGE d'écoulement diphasique non réactif, peut être utilisé pour optimiser la puissance et la position du système d'allumage. / Ignition and altitude reignition are critical issues for aeronautical combustion chambers. The success of ignition depends on multiple factors, from the characteristics of the igniter to the spray droplet size or the level of turbulence at the ignition site. Finding the optimal location of the igniter or the potential of ignition success of a given energy source at a given location are therefore parameters of primary importance in the design of combustion chambers. The purpose of this thesis is to study forced ignition of aeronautical combustion chambers. To do so, Large Eddy Simulations (LES) of two-phase reacting flows are performed and analyzed. First, the equations of the Eulerian formalism used to describe the dispersed phase are presented. To validate the successive LES, experimental data from the MERCATO bench installed at ONERA Fauga-Mauzac are used. It allows to validate the two-phase evaporating flow LES methodology and models prior to its use to other flow conditions. The statistically stationary two-phase flow reacting case is then compared to available data to evaluate the model in reacting conditions. This case is more deeply studied through the analysis of the characteristics of the flame. This last one appears to experience very different combustion regimes. It is also seen that the determination of the most appropriate methodology to compute two-phase flow flame is not obvious. Furthermore, two different methodologies may both agree with the data and still have different burning modes. The ability of the LES to correctly compute burning two-phase flow being validated, LES of the transient ignition phenomena are performed. The experimentally observed sensitivity of ignition to initial conditions, i.e. to sparking time, is recovered with LES. The analysis highlights the major role played by the spray dispersion in the development of the initial flame kernel. The use of LES to compute ignition sequences provides a lot of information about the ignition phenomena, however from an industrial point of view, it does not give an optimal result, unless all locations are tested, which brings the CPU cost to unreasonable values. Alternatives are hence needed and are the objective of the last part of this work. It is proposed to derive a local ignition criterion, giving the probability of ignition from the knowledge of the unsteady non-reacting two-phase (air and fuel) flow. This model is based on criteria for the phases of a successful ignition process, from the first kernel formation to the flame propagation towards the injector. Then, comparisons with experimental data on aeronautical chambers are done and show good agreement, indicating that the proposed ignition criterion, coupled to a Large Eddy Simulation of the stationary evaporating two-phase non-reacting flow, can be used to optimize the igniter location and power.

: Ce travail concerne l'écoulement inertiel en milieu poreux rencontré dans diversessituations telles que les écoulements autour des puits pour la récupération pétrolière, lesécoulements dans les réacteurs catalytiques, etc. En régime stationnaire, les différents modèlesmacroscopiques pour décrire ces écoulements inertiels (non-linéaires) demeurent encore sujetsà débat. Ces modèles consistent en une loi de Darcy corrigée de termes dont la dépendancevis à vis de la vitesse de filtration relève du régime d'écoulement. Dans ce travail, une attentionparticulière est portée tout d'abord à l'étude numérique (DNS), sur des structures modèles, de lalimite de stationnarité de l'écoulement monophasique newtonien qui correspond à la premièrebifurcation de Hopf, caractérisée par un nombre de Reynolds critique. La connaissance de cettelimite est cruciale puisqu'elle détermine le domaine de validité des modèles macroscopiquesstationnaires pertinents. Dans un deuxième temps, la dépendance de la déviation (inertielle) àla loi de Darcy par rapport aux propriétés de la structure poreuse (forme des grains, désordre)et à l'orientation de l'écoulement est étudiée dans le cas de structures 2D et 3D. Les propriétéseffectives de la structure à l'échelle macroscopique sont déterminées à partir de la résolutionnumérique des problèmes de fermeture associés au modèle macroscopique obtenu par prisede moyenne des équations de Navier-Stokes. Afin de déceler l'origine de cette déviation et sesdifférentes formes, l'évolution de la structure microscopique de l'écoulement en fonction dunombre de Reynolds est analysée. Plus particulièrement, le rôle des zones de recirculation, etles corrélations avec la courbure des lignes de courant multipliée par l’énergie cinétique localeet la variation de l’énergie cinétique le long de ces lignes sont étudiés. La dernière partie dutravail est consacrée à une étude numérique, toujours dans des situations modèles, de ladéviation à la loi de Darcy généralisée dans le cas de l'écoulement diphasique inertiel. / This work focuses on inertial flow in porous media encountered in differentindustrial situations such as flow around wells in oil recovery, flow in filters and in columns ofreactors for chemical engineering, etc. In stationary flow regime, the different macroscopicmodels describing inertial (non-linear) flow are still discussed. These models consist in theDarcy’s law with correction extra terms whose dependence upon the filtration velocity is afunction of the flow regime. In this work, a particular attention is attributed first to the numericalinvestigation (DNS), on model structures, of the limit of one phase Newtonian stationary flowwhich corresponds to the first Hopf bifurcation, characterized by a critical Reynolds number.The knowledge of this limit is crucial since it establishes the ranges of validity of the relevantmacroscopic stationary models. In a second step, the dependence of the deviation (inertial)from Darcy’s law on the properties of the porous structure (grains shape, disorder) and on theorientation of the flow is analyzed in 2D and 3D situations. The effective properties of thestructure and the flow at the macroscopic scale are obtained from the numerical resolution ofthe closure problems associated to the macroscopic model obtained from an up-scalingprocedure (volume averaging) of the Navier-Stokes equations. In order to identify the origin ofthe deviation and its different forms, the variation of the microscopic flow structure with theReynolds number is analyzed. More specifically, the role of the recirculation zones, and thecorrelations with flow streamlines curvature multiplied by the local kinetic energy and thevariation of the kinetic energy along these lines are studied. The last part of the work isdedicated to a numerical investigation of the deviation from the generalized Darcy’s law in thecase of two phase inertial flow.

Les simulations des écoulements diphasiques à l’échelle réelle de l’application nécessitent des modèles pour les termes non fermés des équations macroscopiques. Des simulations numériques directes à particule résolue utilisant la méthode de pénalisation visqueuse ont été réalisées afin de mesurer les interactions entre des particules de différentes formes (sphérique et ellipsoïdale) et le fluide porteur à différents régimes d'écoulement (de stokes à l'inertiel). Deux méthodes ont été développées durant cette thèse afin d'extraire les forces hydrodynamiques ainsi que le transfert de chaleur sur les frontières immergées représentant les particules. Plusieurs validations ont été conduites pour différentes configurations de particules : de la simulation d’une particule isolée à un réseau aléatoire de sphères en passant par réseau cubique face centrée de sphères. Une corrélation du nombre de Nusselt est proposée pour un sphéroïde allongé plongé dans un écoulement uniforme. / The simulations of multiphase flows at real application scale need models for unclosed terms in macroscopic equations. Particle-Resolved Direct Numerical Simulations using Viscous Penalty Method have been carried out to quantify the interactions between particles of different shapes (spheres, ellipsoids) and the carrier fluid at different regimes (from Stokes to inertial). Two methods have been developed to extract hydrodynamic forcesand heat transfers on immersed boundaries representing the particles. Validations have been conducted for various configuration of particles: from an isolated sphere and spheroid to Face-Centered Cubic to a random arrangement of spheres. A correlation of the Nusselt number for an isolated prolate spheroid past by a uniform flow is proposed.

Une étude numérique et des investigations expérimentales sont menées sur une bulle de Taylor ascendante au passage d’une contraction ou d’un élargissement brusque dans un liquide Newtonien stagnant. Le code CFD est issu du logiciel libre Gerris utilisant la méthode VOF (Volume Of Fluid) pour représenter l’interface liquide-gaz. La méthode numérique est vérifiée en utilisant les résultats existants sur des bulles de Taylor ascendantes dans des colonnes droites. L’étude expérimentale permet d’étudier une bulle d’azote dans un mélange eau-glycérol pour différentes concentrations. Le rapport des diamètres de la conduite varie de 0,69 à 1,72. Les images de la bulle montante dans les singularités sont capturées par une caméra haute vitesse. Les variations de vitesse de la bulle et de l’épaisseur du film liquide sont étudiées. Les résultats montrent que l’élargissement brusque avec un rapport de diamètres plus élevé entraîne plus de perturbations sur la bulle avec de forts effets sur sa queue. Les queues de bulles instables sont coupées en petites bulles dans certains cas et une carte d’écoulement permet de le prédire. En outre, il a été observé que les variations de la forme de la bulle dépendent de la longueur des bulles. En ce qui concerne la bulle passant par une contraction, le phénomène de blocage a été observé et une carte de prédiction a été proposée. Finalement, cette étude, basée sur un intervalle assez large du nombre d’Eötvös et du rapport de diamètres, propose de nouvelles connaissances pour mieux comprendre l’ascension d’une bulle de Taylor dans les singularités. / Numerical simulations and experimental investigations of rising individual Taylor bubble through a vertical sudden expansion and contraction in stagnant Newtonian liquid is presented. The CFD procedure is based on the open source package Gerris which adopts the volume-of-fluid (VOF) method to represent the gas/liquid interface. The numerical method is verified using the existing results of single Taylor bubbles rising in straight columns. The experiments investigate a nitrogen bubble rising in a water-glycerol mixture for different concentrations. The pipe diameter ratio ranges from 0.69 to 1.72. The images of the bubble rising through the singularities are captured by a high-speed camera. Our investigations focus on the transient process of the bubbles passing the singularity. The variations of the bubble velocity and the liquid film thickness are investigated. The results show that the greater expansion ratios yield more perturbations on the bubbles and have strong effects on the tail of the bubble. The unstable bubble tails are cut off into smaller bubbles in some of the test cases and a bubble break-up regime map obtained by simulations has been proposed. The bubble shape variations depend also on the length of the bubbles. For a bubble passing through a contraction, the blocking phenomenon has been observed and a map has been proposed. Finally, this study, based on a large range of Eötvös numbers and expansion/contraction ratios, provides new insights to better understand the effect of singularities on rising Taylor bubbles.

Cette thèse est dédiée à l'étude expérimentale et numérique des phénomènes de condensation dans un pressuriseur de chaufferie nucléaire embarquée conçu par l'entreprise DCNS. Au sein du système est observé un écoulement diphasique avec changement d'état, descendant dans un tube vertical partiellement immergé dans du liquide. Le liquide sous-refroidi s'écoule dans le tube sous l'effet de la gravité en entrainant de la vapeur saturée qui se condense à son contact. Dans un premier temps, l'étude des phénomènes de condensation est réalisée sur un banc expérimental dont la conception, le dimensionnement, la réalisation et l'exploitation ont été effectués au cours de cette thèse. Le dispositif est une version simplifiée à l'échelle 1 du système réel. Cette étude expérimentale a permis d'identifier les régimes d'écoulements diphasiques présents dans le système et d'analyser l'impact de trois paramètres : le débit de liquide injecté, le sous-refroidissement de ce liquide à l'entrée du tube et la pression en sortie de tube. Par ailleurs, l'étude propose un modèle déterminant les conditions pour laquelle le tube d'essais est complètement noyé, ce qui minimise les phénomènes de condensation. Dans un second temps, une simulation numérique de l'écoulement diphasique est réalisée par CFD sur ANSYS Fluent à l'aide du modèle " Volume Of Fluid ". La prise en compte de la condensation de la vapeur au contact du liquide sous-refroidi a nécessité le développement, grâce aux fonctions utilisateur du logiciel, d'un modèle spécifique s'appuyant sur un modèle de la littérature. Les résultats obtenus ont permis de reproduire qualitativement le déversement du liquide observé expérimentalement en entrée de tube ainsi que la condensation de la phase vapeur au contact de la phase liquide. / This thesis focuses on condensation phenomena occurring in a patented system designed by DCNS Company. The system ensures pressure regulation of nuclear boiler rooms embarked in naval vessels. Within the system occurs a downward liquid/vapor flow with phase change in a partially submerged vertical pipe. Subcooled liquid flows into the pipe and induce a suction of saturated vapor which condense in contact with liquid phase. This study aim to understand the physical phenomena occurring in the system and to simulate these phenomena. First, the condensation phenomena are studied with a test loop which conception, dimensioning, realization and exploitation were realized during this study. The experimental study identified the flow regimes occurring in the system and analyzed the impact of three parameters: the liquid flow rate, the liquid subcooling and the pressure at the system outlet. This study also proposed a simple model which determine the necessary conditions to completely drown the pipe, which minimize condensation phenomena. Thereafter, the numerical simulation of the two-phase flow was realized by CFD on ANSYS Fluent with the "Volume of Fluid" model. The condensation simulation was achieved by applying a condensation model used in literature. This model was added to the numerical code with User-Defined Functions (UDF). The simulation reproduced the weir regime flow at the pipe entrance and the vapor condensation inside the pipe.

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.

Le travail présenté dans ce mémoire concerne l'étude du mouvement d'inclusions solides en suspension dans un écoulement turbulent de canal isotherme et anisotherme par simulation numérique directe. Grâce à cet outil de simulation, nous avons pu analyser l'influence de l'inertie et du croisement de trajectoires sur les caractéristiques dynamiques et thermiques de la phase dispersée ainsi que sur celles du fluide vu au sein d'une turbulence non homogène. Cela pourra notamment faciliter le développement futur des modélisations euléro-lagrangienne et euléro-eulérienne dans ce type de turbulence. Nous avons essayé d'examiner au mieux la conséquence de ces effets sur les statistiques de la phase dispersée telles que la concentration, la moyenne et l'écart type de la vitesse, les covariances fluide-particules, les corrélations triples de la vitesse des particules. Parallèlement à cela, les caractéristiques du fluide vu par les particules ont été étudiées, nous nous sommes intéressés à la vitesse de dérive, aux tensions de Reynolds du fluide vu, et la décorrélation temporelle des fluctuations de la vitesse du fluide vu. En ce qui concerne la partie thermique, nous présentons les statistiques thermiques de la phase dispersée et du fluide vu par les particules au sein de l'écoulement anisotherme vertical descendant. Tous ces aspects sont développés en gardant à l'esprit le cadre général de cette étude, c'est-à-dire comprendre le comportement thermique de la phase dispersée et fournir des informations concernant des grandeurs importantes intervenant dans la modélisation de tels écoulements.

Le sujet de cette thèse est la modélisation de la combustion turbulente en milieu diphasique telle qu'elle a lieu par exemple dans les moteurs Diesel. Certaines expériences notamment celles présentées dans cette thèse suggèrent que quand le combustible nécessaire à la combustion est sous forme d'un nuage de gouttes, les fluctuations en particulier de richesse dans le milieu gazeux sont plus fortes. Cela doit être pris en compte dans la modélisation et pourrait être dû à la présence des gouttes. Celles-ci en s'évaporant, jouent le rôle de sources de combustible ponctuelles, aléatoirement réparties dans l'espace, ce qui crée plus de fluctuations à l'échelle de l'espace intergoutte. Ces fluctuations supplémentaires sont prises en compte ici par un groupe de termes lié à l'évaporation, dans l'équation de la variance de la variable de mélange, qui sont fermés en considérant les équations locales d'un milieu diphasique. L'importance d'une nouvelle variable décrivant l'état du mélange à la surface des gouttes au moment de l'évaporation d'une particule fluide a été montrée. Une équation fermée a été écrite pour cette variable. Ces deux nouvelles caractéristiques permettent de mieux décrire l'état du mélange. Pour représenter la combustion nous avons d'abord utilisé un modèle qui suppose que la chimie est infiniment rapide. L'état du milieu thermochimique est alors complètement dépendant de l'évaporation et du mélange sans influence de la cinétique chimique. Le modèle MIL de Gonzalez et Borghi est ensuite utilisé pour prendre en compte certains effets chimiques. Il suppose que la chimie est infiniment brusque, c'est-à-dire à température d'activation très élevée. Ces deux modèles ont été testés pour deux configurations expérimentales, une flamme de type jet et une bombe simulant la combustion dans les moteurs Diesel. Les taux de fluctuation trouvés lors de ces simulations sont assez élevés et mettent en avant la nécessité de bien estimer la dissipation des fluctuations scalaires. Les résultats trouvées en utilisant ces modèles sont comparables à ceux mesurés, notamment dans le cas de la bombe de type Diesel, le modèle MIL permet de retrouver un temps d'auto-allumage et une longueur de la flamme réaliste.

Les simulations et l'étude expérimentale présentée dans ce mémoire ont pour but d'améliorer les diagnostics optiques afin de comprendre la physique de la pulvérisation et l'évaporation des gouttes. L'anémométrie Phase Doppler permet déjà de connaître la corrélation taille-vitesse dans un milieu relativement dilué. L'objectif des simulations numériques expliquées dans ce mémoire est de présenter la faisabilité d'étendre les mesures PDA à la mesure de d'indice de réfraction des gouttes. L'extension proposée, appelée l'anémométrie Phase Doppler à point sde gloire est basée sur l'analyse de la lumière diffusée par les points de gloires de la réflexion et de la réfraction d'une goutte se séplaçant dans un volume de mesure type PDA. Une des autres difficulté concernant l'étude de l'atomisation est la forte densité des sprys à la sortie des injecteurs de types automobile. la vélocimétrie laser par corrélation permet de connaître la vitesse des structures par l amesure d'un temps de vol. Dans le cadre de ce travail de thèse, cette technique a été validée et implantée au laboratoire.

On s'intéresse dans ce travail à la simulation et l'approximation d'écoulements diphasiques avec transition de phase. Il s'agit de modéliser la dynamique d'une bulle de cavitation. Le modèle repose sur les équations d'Euler en coordonnées sphériques, l'interface gaz-liquide étant indiquée par une fonction de couleur. Dans une première partie, aucun transfert de masse n'est supposé. Un schéma Lagrange-projection est d'abord proposé : seule l'interface est déplacée à la vitesse du fluide. La projection repose sur un algorithme de suppression-découpage qui assure que les volumes ne deviennent pas négatifs. Le second chapitre traite du terme source géométrique. On construit un schéma équilibre articulé autour du schéma VFRoe-ncv pour lequel on propose une correction entropique non paramétrique. Un méthode d'ordre élevé de type Galerkin discontinu est ensuite étudiée dans le cadre de la magnétohydrodynamique. L'intégration en temps est réalisée par une méthode de type Adams-Bashforth, qui s'avère bien adaptée aux algorithmes de pas de temps local. La deuxième partie de la thèse traite de la modélisation du changement de phase liquide-vapeur. L'inf-convolution et la transformée de Legendre définissent une structure naturelle pour la construction de loi de pression de mélange. En particulier on montre que la construction de Maxwell est équivalente à la construction de l'enveloppe convexe de l'énergie de van der Waals. Un algorithme de transformée de Legendre rapide, adapté à la prise en compte correcte des bords du domaine de calcul de loi d'état, est développé. La méthode est appliquée à la construction de lois de pression tabulées de mélanges binaires miscibles ou immiscibles de gaz raides.