La conception de nouveaux moteurs impose de respecter des normes de sécurité concernant les performances d'allumage et de ré-allumage en conditions critiques. Des campagnes d'essais étant onéreuses, les industriels cherchent donc à disposer d'outils numériques fiables. Afin d'améliorer la simulation des écoulements, le caractère multicomposant du carburant doit être pris en compte. L'objectif de cette thèse est d'étudier l'influence de l'évaporation d'un brouillard de gouttes sur un écoulement réactif. Pour cela, une étude de la propagation d'une flamme laminaire 1D est réalisée à l'aide d'un code de calcul multiphysique (CEDRE). Un train continu de gouttes monodisperse est injecté, les gouttes étant mono ou bicomposant. L'influence de la dynamique d'évaporation sur la combustion est étudiée. Deux cinétiques chimiques réduites multicomposant sont comparées. La composition, le diamètre et la richesse initiale des gouttes ont un impact sur la structure de flamme, la vitesse de flamme et la composition des gaz brûlés. Ensuite, l'effet de l'évaporation est étudié en phase d'allumage pour un brouillard de gouttes polydisperses monocomposant avec un modèle de noyau d'allumage local. L’écoulement instationnaire non-réactif dans un secteur de chambre industriel (MERCATO) est calculé avec une approche LES. Le caractère instationnaire, voire périodique, de la phase dispersée est mis en évidence en certains points de l'écoulement. Les résultats, associés au modèle d'allumage et à des critères, sont utilisées pour réaliser une carte de probabilité d'allumage. Des essais de calcul d'allumage complet de la chambre sont réalisés. Les résultats indiquent une surestimation des termes sources liés à l'évaporation de la phase dispersée et à la combustion. / The design of new aircraft engines needs in particular to comply with safety standards for the performance of stabilized combustion and ignition or re-ignition under critical conditions. Experimental campaigns are expensive, so numerical tools are needed. To improve the accuracy of the models used to simulate flow, the multicomponent nature of the fuel must be taken into account, whether it is kerosene or alternative fuel. The objective of this thesis is to study the influence of a droplet mist vaporization on a reactive flow. For this, an academic study of the propagation of a 1D laminar flame is performed using a CFD code {CEDRE). A continuous stream of monodisperse droplets is injected, the droplets being mono or bicomponent. The influence of the dynamics of evaporation on combustion is particularly studied. Two reduced multicomponent chemical kinetics are compared. The composition, the diameter and the initial equivalent ratio of droplets have an impact on the structure of the flame, the flame speed and composition of the burnt gases. A local ignition kernel model is applied to study the influence ofvaporization on ignition in the case of monocomponent, polydisperse droplets. Experimental data are available for a monosector combustion chamber (MERCATO) so the non-reactive unsteady flow is simulated with a LES approach. The unsteady, sometimes periodic, nature of the dispersed phase is highlighted in some points of the flow. A ignition model is applied to instantaneous flow fields and criteria are analysed to realise an ignition probability map which validates the approach. Finally, ignition of a combustion chamber is tested. The results point out an overestimation of source terms related to the evaporation of the dispersed phase and combustion.

Cette thèse présente une contribution à l’étude des transferts d’eau au sein des piles à combustible de type PEMFC, un aspect clé de cette technologie. Une approche de simulation numérique est développée en couplant un modèle de type réseau de pores dans la couche de diffusion (DM), une approche mixte continue –réseau de pore dans la couche microporeuse (MPL) et une modélisation par compartiments dans la couche active. L’approche développée prend en compte les transferts couplés de chaleur et d’eau via notamment la modélisation des phénomènes de changement de phase dans la DM et la MPL (évaporation et condensation). Dans une première partie, nous étudions le cas où l’eau migre dans l’assemblage MPL-DM directement en phase liquide. L’impact de la variation de pression dans la phase gazeuse sur la distribution de la phase liquide est étudié. L’épaisseur optimale de la MPL est également étudiée. Dans une seconde partie, nous étudions des situations où l’eau se forme par condensation dans la couche de diffusion. Nous étudions tout d’abord l’impact des propriétés de la couche de diffusion et de la MPL sur le diagramme de condensation. Ensuite nous analysons l’impact de la formation de l’eau liquide sur la distribution de courant locale. Enfin, l’impact de la mouillabilité sur les figures de condensation est étudié. Cette dernière étude est vue comme un premier pas vers l’étude des mécanismes de dégradation dans le régime de condensation. / This thesis is a contribution to the study water transfers within PEM fuel cell, a key element of this technology. A numerical simulation tool is developed coupling a pore network model in the gas diffusion layer (DM), a mixed continuum – pore network approach in the microporous layer (MPL) and a model by compartments in the catalyst layer. The developed approach takes into account the coupled heat and water transfers through the modeling of phase change phenomena (evaporation – condensation) in the DM and in the MPL. In the first part, we study the case where water migrates into the MPL-DM assembly directly in liquid phase. The impact of gas pressure variation on liquid phase distribution is studied. The optimal thickness of MPL is studied too. In the second part we study situations where liquid water essentially formed by condensation in the diffusion layer. We first study the impact of DM and MPL properties on the condensation diagram. Then we analyze the impact of liquid water formation on the local current density distribution. Finally the impact of wettability modifications on the liquid water patterns is studied. This last study is considered as a first step toward the study of degradation mechanisms in the condensation regime.

Un grand effort a été réalisé ces dernières années dans la compréhension du champ d'écoulement et du transfert de chaleur dans les canaux de refroidissement internes présents dans les pales de turbine. En effet, des systèmes de refroidissement avancés ont non seulement conduit à l'augmentation de l'efficacité de la turbine à gaz en augmentant la température d'entrée de la turbine au-dessus de la température de fusion du matériau, mais également en augmentant la durée de vie de la turbine. Pour permettre de tels progrès, des techniques expérimentales et numériques modernes ont été largement appliquées afin d'interpréter et d'optimiser l'aérodynamique et le transfert de chaleur dans les canaux de refroidissement internes. Cependant, les données disponibles sont limitées dans le cas des canaux de refroidissement internes dans les aubes de rotor de turbine. Les gradients de rotation et de température introduisent des forces de flottabilité de type Coriolis et centripète dans le référentiel rotatif, modifiant de manière significative l'aérothermodynamique par rapport aux passages stationnaires. Dans le cas des pales de rotor de turbine, la plupart des investigations sont soit basées sur des mesures ponctuelles, soit sont contraintes à des régimes de rotation faibles. L'objectif principal de ce travail est d'étudier le débit détaillé et le transfert de chaleur d'un canal de refroidissement interne à des conditions de fonctionnement dimensionnelles sans moteur représentatives. Ce travail introduit une section d'essai en laboratoire qui exploite des canaux à nervures sur un large éventail de nombres de Reynolds, de rotation et de flottabilité. Dans le présent travail, le nombre de Reynolds va de 15,000 à 55,000, le nombre de rotation maximum est égal à 0.77 et le nombre maximal de flottabilité est égal à 0.77. La nouvelle installation expérimentale consiste en une conception polyvalente qui permet l'interchangeabilité de la géométrie testée, de sorte que les canaux de différents rapports d'aspect et les géométries de nervure peut être facilement installé. La particle image velocimetry et la thermographie à cristaux liquides sont effectuées pour fournir des mesures précises de vitesse et de transfert de chaleur dans les mêmes conditions opératoires, ce qui conduit à un ensemble de données expérimentales unique. De plus, des simulations à grands virages sont réalisées pour donner une image de l'ensemble du champ d'écoulement et compléter les observations expérimentales. En outre, l'approche numérique vise à fournir une méthodologie robuste qui est capable de fournir des prédictions haute-fidélité de la performance des canaux de refroidissement internes. / A great effort has been carried out over the recent years in the understanding of the flow field and heat transfer in the internal cooling channels present in turbine blades. Indeed, advanced cooling schemes have not only lead to the increase of the gas turbine efficiency by increasing the Turbine Inlet Temperature above the material melting temperature, but also the increase of the turbine lifespan. To allow such progresses, modern experimental and numerical techniques have been widely applied in order to interpret and optimize the aerodynamics and heat transfer in internal cooling channels. However, the available data is limited in the case of internal cooling channels in turbine rotor blades. Rotation and temperature gradients introduce Coriolis and centripetal buoyancy forces in the rotating frame of reference, modifying significantly the aerothermodynamics from that of the stationary passages. In the case of turbine rotor blades, most of the investigations are either based on point-wise measurements or are constraint to low rotational regimes. The main objective of this work is to study the detailed flow and heat transfer of an internal cooling channel at representative engine dimensionless operating conditions. This work introduces a laboratory test section that operates ribbed channels over a wide range of Reynolds, Rotation and Buoyancy numbers. In the present work, the Reynolds number ranges from 15,000 to 55,000, the maximum Rotation number is equal to 0.77, and the maximum Buoyancy number is equal to 0.77. The new experimental facility consists in a versatile design that allows the interchangeability of the tested geometry, so that channels of different aspect ratios and rib geometries can be easily fitted. Particle Image Velocimetry and Liquid Crystal Thermography are performed to provide accurate velocity and heat transfer measurements under the same operating conditions, which lead to a unique experimental data set. Moreover, Large Eddy Simulations are carried out to give a picture of the entire flow field and complement the experimental observations. Additionally, the numerical approach intends to provide a robust methodology that is able to provide high fidelity predictions of the performance of internal cooling channels.

Le risque individuel de rupture des anévrismes cérébraux est un enjeu majeur dans la prise en charge clinique des anévrismes asymptomatiques. La rupture anévrismale se produit lorsque la contrainte intra-pariétale dépasse la contrainte à rupture du matériau composant la paroi. Notre étude a pour objectif d'être un pas vers une nouvelle mesure biomécanique du risque individuel de rupture des anévrismes cérébraux. Dans un premier temps, une étude expéri- mentale fut menée pour caractériser le comportement biomécanique de la paroi anévrismale sur 16 échantillons d'anévrismes prélevés chirurgicalement. L'expérimentation sur les échan tillons de poche anévrismale a permis de dégager trois grandes classes de tissus pour chaque sexe (homme et femme) : souple, rigide et intermédiaire. Il apparaît que tous les anévrismes non rompus appartiennent à la catégorie rigide ou intermédiaire et que tous les anévrismes rompus correspondent à la catégorie souple. Ceci permet de mettre en évidence une corrélation entre le risque de rupture et les propriétés du matériau composant la paroi anévrismale. Dans un deuxième temps, des simulations d'interaction fluide/structure (FSI) ont été réalisées pour comparer les déformations d'un anévrisme " patient spécifique " constitué d'un matériau dégradé et non dégradé. Les résultats montrent que les propriétés du matériaux ont un impact majeur sur l'ampleur de la variation de volume anévrismale diastolosystolique. Les changements en terme de variations de volume en fonction des caractéristiques du tissu sont potentiellement visualisables à l'aide de l'imagerie médicale. Une analyse des incertitudes des paramètres est aussi présentée et montre la robustesse des résultats aux incertitudes des données d'entrée. Il a ensuite été démontré sur 12 cas " patient-spécifique " d'anévrismes différents (forme, taille, localisation et conditions aux limites différentes) qu'il existe toujours une différence significative en terme de variation de volume au cours du cycle cardiaque entre un anévrisme dont la paroi est composé d'un matériau rigide et d'un matériau souple. Cette étude suggère donc que la variation de volume anévrismale pourrait être utilisée comme une base pour une évaluation individuelle du risque de rupture des anévrismes cérébraux. / The individual risk of rupture of cerebral aneurysm is a major stake in the clinical treatment. The aneurismal rupture occurs when the intra-parietal stress exceeds the rupture stress of the material of the aneurismal wall. The goal of our study is to be a step toward a new biomechanical measure of an individual risk of rupture of cerebral aneurysm. First, an experimental study was performed to characterize the biomechanical behavior of the aneurismal wall on 16 samples of aneurysms removed by neurosurgery. The experimentation on the samples allowed us to reach three main categories of tissues for each sex (female and male): soft, intermediate and stiff. All the unruptured aneurysms belong to the stiff category or the intermediate category and all the ruptured aneurysms belong to the soft category. This is allowed us to give prominence to the correlation between the risk of rupture and the properties of the material of the aneurismal wall. Then, Fluid/Structure interaction computations (FSI) were performed to compare the strain of a “patient-specific” aneurysm composed of a degraded and undegraded material. The results show that the properties of the material have a major impact on the scope of the aneurismal volume variation over the cardiac cycle. The volume variation changes depending on the properties of the tissue are potentially viewable by medical imaging. A study of the uncertainties of the parameters is also proposed and shows the robustness of the results. We also demonstrated on 12 cases of “patient-specific” aneurysms that a significant difference stiff exists in terms of volume variation over the cardiac cycle between an aneurysm composed of a stiff and a soft material. This study suggests that the aneurismal volume variation could be used as a basis for an evaluation of the individual risk of rupture of cerebral aneurysms.

Dans les travaux de recherche développés durant cette thèse, une première approche desécoulements magnétohydrodynamiques (MHD) multiphasiques est proposée. Cette approche seconcentre sur les phénomènes liés à l’interaction entre une dynamique interfaciale et un écoulementde coeur MHD. Le couplage induit entre rhéologie de surface etMHDde sous-phase est particulièrementillustré par le développement d’un viscosimètre annulaire surfacique, dédié à l’étude desmétaux liquides progressivement oxydés. En premier lieu sont introduits les éléments théoriquespropres à laMHDet à la rhéologie de surface. La modélisation de leur couplage fait intervenir deuxparamètres interfaciaux : les viscosités surfaciques dilatationnelle et de cisaillement. L’influence respectivede ces deux paramètres sur l’écoulement MHD de sous-phase est étudiée analytiquementet numériquement dans le cas d’un écoulement (stratifié) annulaire MHD permanent. Dans la configurationretenue, un champ magnétique uniforme vertical est imposé, perpendiculairement à lasurface liquide graduellement oxydée. Le rôle décisif des contraintes visqueuses interfaciales concernantl’(in)activation des couches de Hartmann est démontré, conduisant à des topologies atypiquesd’écoulement MHD. Le viscosimètre annulaire MHD est ensuite proposé en tant que méthodeexpérimentale originale, permettant la mesure sélective des viscosités surfaciques de fluides électroconducteurs.Les premières campagnes expérimentales aboutissent à une estimation de la viscositéde cisaillement interfacial d’un alliage métallique (GaInSn) liquide à température ambiante. Finalement,une ouverture sur l’écoulement MHD à proximité d’une inclusion gazeuse sphérique rigideest discutée en annexe de ce projet, en lien avec des conditions mécaniques variables à l’interfaceliquide/gaz pilotées par la rhéologie de surface, constituant une première approche vers la descriptiondes écoulements MHD dispersés. / In this thesis work, a first approach towards the description of magnetohydrodynamic(MHD) multiphase flows is proposed, based on the investigation of the role of surface rheologyin permanent regime. The study of the coupling between bulk MHD and surface rheology is particularlymotivated by the development of an annular surface viscometer devoted to liquid metalstopped with an oxidation layer. First, theoretical foundations of MHD and of surface rheologyare separately introduced. The modelling brings out a strong coupling between bulk and surfacevelocities, the latter being dependent on two interfacial parameters: the surface shear and dilatationalviscosities. Their respective influence is analytically and numerically investigated in the caseof a permanent annular (stratified) MHD end-driven flow. In the considered geometry, a vertical(uniform) magnetic field is imposed, perpendicular to the gradually oxidising liquid surface. Thecontribution of planar surface viscous stresses to the possible electrical activation of Hartmannlayers is demonstrated, leading to a variety of atypical MHD flow patterns. The annular MHDviscometer is then developed as a first device able to perform selective measurement of the surfaceviscosities of electroconductive fluids. First experiments lead to an estimated value of the surfaceshear viscosity for a liquid alloy (GaInSn) at room temperature. Some hints are finally given to investigatethe study case of aMHDflow past a rigid fluid sphere, with varying interfacial conditionsgoverned by surface rheology, as an additionalwork constituting a first step towards the descriptionof dispersed MHD flows.

Dans l’industrie chimique, la majorité des réacteurs est équipée d’un dispositif « évent de sécurité » permettant d’éviter leur éclatement en cas de surpression accidentelle conséquente à un emballement. La nature physique du rejet (gazeux ou diphasique gaz-liquide) influe fortement sur la taille requise de l’évent. L’objectif de notre étude est de connaitre mieux le comportement de l’écoulement (désengagement) en cas d’un emballement de réaction et de prédire la nature du rejet (mono ou diphasique) notamment à l’échelle industrielle.Une étude expérimentale d’une réaction d'estérification fut réalisée dans un réacteur en verre de 0,5 l afin de visualiser la nature de l’écoulement et de mesurer la fraction de vide moyenne (ᾱ) lors de l’emballement de réaction. Une étude paramétrique fut effectuée de façon à identifier, pour le régime d'écoulement et la transition, les paramètres clés : la viscosité et l’agitation. Un calorimètre pseudo-adiabatique (VSP2) fut utilisé pour étudier la thermodynamique de l’emballement d’estérification. Les bilans massique et énergétique ont été utilisés pour calculer la vitesse superficielle de la vapeur (jg,max). En outre, le jg, max fut calculé aussi à partir de différentes corrélations, issues d’études en colonne à bulles. Une comparaison entre les résultats obtenus a montré l’inadéquation de ces corrélations pour un système réactif. Une carte de régime pour un système réactif a été construite pour la première fois (jg,max versus ᾱ, combinée à l'observation du régime d'écoulement dans le réacteur), avec une méthodologie pour la prédiction de la nature du rejet en cas d’emballement. / In chemical industry, most reactors are equipped with an emergency relief vent to prevent bursting in case of accidental overpressure due to a runaway reaction scenario. The physical nature of the vent release (gas phase or gas-liquid) strongly influences the necessary vent size. The objective is to enhance the knowledge on the flow behavior (disengagement) during a runaway reaction and to be able to predict the nature of the vent flow (1 or 2-phase) namely at industrial scale.Experiments of esterification were done in a 0.5 l glass reactor in order to visualize the flow pattern (hydrodynamic) and to measure the average void fraction (ᾱ) during the runaway. A parametric study was carried out to identify the key parameters onto the flow pattern and transition, they were: the viscosity and the stirring. A pseudo-adiabatic calorimeter (VSP2) was used to study the runaway P and T profile and the kinetic of the esterification. The thermodynamic data and the mass & heat balances were used to calculate the superficial vapor velocity (jg,max). Furthermore, this jg,max was also calculated from different correlations, obtained from bubble column research. A comparison showed the incorrectness of using these correlations for a reactive system. A flow pattern map (jg,max versus ᾱ combined with the observation of the flow pattern,) for a reactive system was built up for the first time together with a methodology to predict the vent flow nature.

Cette thèse expérimentale porte sur des impacts de gouttes en caléfaction, aussi appelée régime Leidenfrost. Dans ce cas, la goutte est isolée thermiquement et mécaniquement du substrat surchauffé par une fine couche de vapeur. On s'affranchit ainsi de la friction visqueuse. Les substrats présentent des textures micrométriques localisées. On retrouve un régime autosimilaire d'étalement aux temps courts. On caractérise des régimes de recouvrement d'un défaut ponctuel par la goutte. Ces régimes sont dictés par l'épaisseur de lamelle par rapport à celle du défaut. Les défauts génèrent des excroissances dans leur sillage dont la dynamique peut être approchée selon deux modèles inertio-capillaires ; l'un valable aux temps courts, l'autre aux temps plus longs. En présence de plusieurs défauts, on fragmente la lamelle depuis plusieurs sites selon le ratio entre largeur des défauts et épaisseur de la lamelle. On simule par un algorithme de pavage le motif final sur lequel le fluide se concentre à la fin de la fragmentation. Grâce à des rugosités plus complexes on peut canaliser l'étalement de la goutte. On exacerbe alors l'étalement maximal dans l'axe des canaux. On peut aussi inhiber l'étalement par des textures circulaires. Les rugosités affectent le temps avant rebond de la goutte. On exhibe une dépendance générale unique entre temps avant rebond et étalement maximal. La dépendance est valable lorsque les textures exacerbent l’étalement, aussi bien que lorsqu’elles l’inhibent / The presented work deals with drop impacts in Leidenfrost regime. In such a case, the drop is thermally and mechanically isolated from the overheated substrate by a thin vapor layer. Viscous friction can therefore be neglected. The substrates are shaped with localized micrometric textures. We rediscover a self-similar spreading regime at short times. We characterize covering regimes of a single defect by the drop. Those regimes are driven by the ratio between lamella thickness and defect thickness. Defects give rise to excrescences in their wake, whose dynamic can be approached by two inertial-capillary recoil models. One is valid at short times and the other at longer times. In the case of several defects, we break up the lamella from different sites according to the ratio between defect width and lamella thickness. We numerically predict, with a tessellation algorithm, the pattern on which the fluid is localized at the end of the fragmentation. Through more complex textures, we can channel the drop spreading. The spreading is increased in the directions of the channels. One can also inhibit the spreading with circular textures. The textures affect the time before drop rebound. We exhibit a general and unique dependency between time before rebound and maximal spreading. This dependency is valid when textures increase the spreading as well as when they inhibit it

Ce travail concerne l’étude d’écoulements incompressibles soumis à la rotation avec un solveur haute précision dans des géométries semi-complexes. La technique numérique mise en œuvre combine des schémas compacts, une méthode de projection multi domaine directe et un traitement efficace de la singularité à l’axe basé sur des conditions de parité dans l’espace de Fourier. Le solveur a été parallélisé avec une approche hybride MPI-OpenMP pour réduire les temps de calcul. Dans un premier temps, les précisions spatiales et temporelles de la méthode numérique et la scalabilité du solveur ont été vérifiées. La capacité du solveur à traiter des écoulements plus complexes a été évaluée en considérant des écoulements de type éclatement tourbillonnaire et un écoulement turbulent en conduite cylindrique. Dans un second temps, plusieurs écoulements typiques des machines tournantes ont été étudiés. Le premier écoulement est un écoulement turbulent incompressible isotherme dans un étage simplifié d’un compresseur haute pression d’une turbine à gaz. Les simulations menées ont mises en évidence l’effet de la rotation sur l’écoulement, notamment sur les instabilités se développant le long des parois et sur les différentes structures cohérentes. Le second cas traité est un écoulement turbulent de jet impactant un disque en rotation avec un fort confinement et transfert thermique. Une attention particulière a été portée sur les champs hydrodynamiques et thermique le long du rotor. Enfin, une étude préliminaire d’un jet turbulent impactant un disque fixe d’épaisseur non nulle dans une configuration moins contrainte avec prise en compte du couplage conduction-convection a été réalisée. / This work deals with the study of rotating incompressible flows with a high accurate solver in semi complex geometries. The numerical method used in this work combines compact schemes, a direct multidomain projection method and an efficient axis treatment based on parity conditions in Fourier space. The use of cylindrical coordinates introduces this mathematical singularity. In order to reduce the calculation time, the solver was parallelized with an hybrid MPI-OpenMP parallelization. First, the spatial and temporal accuracies of the numerical method and the scalability of the solver were checked. Then, the capability of the algorithm to deal with more complex flows was verified. Vortex breakdown flows and turbulent pipe flow were studied. In the second step, typical flows of turbomachineries and rotating systems were considered. The first flow was an incompressible isothermal turbulent flow in a high pressure compressor of gas turbine. The different simulations highlighted the rotation effects on the flows, especially on the instabilities appearing along the walls and the coherent structures. The second considered flow was a turbulent impinging jet on a rotating disk with heat transfer in a small aspect ratio cavity. The hydrodynamic fields and heat transfer near the rotor were analyzed in detail. Finally, a preliminary investigation of an impinging jet on a non-rotating disk in a larger aspect ratio cavity with a coupling between conduction and convection transfer was carried out.

Lors de la phase de rentrée atmosphérique d'une sonde spatiale, la paroi du corps est le siège de phénomènes physico-chimiques complexes. Nous nous intéressons dans cette thèse au cas où le matériau solide de l'objet de vol comporte plusieurs constituants s'ablatant de façon différentielle. En particulier, l'un de ces constituants subit un changement de phase donnant lieu à l'apparition d'une phase liquide. Nous sommes en présence de trois phases : solide, liquide et gaz. Les travaux effectués dans cette thèse correspondent au développement de méthodes numériques en 2D capables de modéliser les différentes interfaces en présence ainsi que l'évolution des fluides ou des matériaux séparés par celle-ci. L'enjeu principal de la thèse est de proposer des méthodes et des algorithmes de couplage pour l'écoulement diphasique, la thermique multimatériaux et les changements de phase (fusion et sublimation) / During atmospheric reentry phase of a spacecraft, its body surface is the seat of complex physico-chemical phenomena. We focus in this thesis on the case where the wall of the flying object has several components ablating differentially. In particular, one of those components undergoes a phase change giving the rise to the introduction of a liquid phase. We have three phases in the domain: solid, liquid and gas phases.The work done in this thesis corresponds to the development of 2D numerical methods which can modelize the different interfaces. The main issue of this thesis is to propose methods and algorithms for coupling the two-phase flow, multi-material heat problems and phase changes (melting and sublimation).

Dans cette thèse, nous nous sommes intéressés à la mesure et à la modélisation des propriétés d'écoulements diphasiques (perméabilités relatives) dans des milieux poreux carbonatés représentatifs de la zone de transition de réservoirs pétroliers carbonatés. Notre étude expérimentale a porté sur deux roches carbonatées de même minéralogie mais de propriétés pétrophysiques différentes (unimodale et bimodale). Les fluides utilisés sont une saumure, une huile minérale pour les expériences en mouillabilité franche à l'eau et une huile brute pour les expériences en mouillabilité altérée (obtenue par macération). Durant les expériences en « écoulements permanents », la perte de charge est mesurée continûment et le champ de saturation locale est mesuré en fin d'expérience par atténuation gamma. Les principaux résultats expérimentaux de cette étude sont les suivantes :- La mouillabilité dépend de la saturation initiale en huile et donc de la hauteur de la zone de transition. Elle a une influence sur les perméabilités relatives à l'huile ET à l'eau.- L'hystérésis des perméabilités relatives à l'eau est d'autant plus prononcée que la saturation initiale en huile augmente. Une observation similaire est faite sur la courbure des perméabilités relatives à l'huile.- Les évolutions des Sorw en fonction des Soi ne vérifient pas toutes, une loi monotone de type Land. Leur évolution dépend aussi de la structure poreuse.- Les deux types de roches carbonatés présentent une altération de mouillabilité différente, la structure vacuolaire favorisant l'altération de mouillabilité. Nos données expérimentales ont été interprétées en utilisant différents modèles d'hystérésis de la littérature. Il apparait clairement des différences entre les prédictions numériques et nos données expérimentales car ces modèles ne prennent pas en compte l'évolution de la mouillabilité. Nous avons proposé un modèle de calcul des perméabilités relatives qui décrit nos données expérimentales de manière satisfaisante. / In this work, we take an interest in measuring and modeling the multiphase flow properties (relative permeabilities) in different carbonate rock types that are representative of the transition zone of carbonate reservoirs. Our experimental study focused on two different carbonate rock types, with almost the same mineralogy, but different petrophysical properties (bimodal and unimodal). The system of fluids used is composed of brine, a mineral oil (for experiments at water-wet conditions), and a crude oil (for experiments including wettability alteration after an ageing process at high temperature). During the experiments performed by the steady-state core flooding method, the pressure drop is continuously measured along the core, and the saturation profiles are obtained at the end of the flooding, by gamma ray measurements. The main results of this study are as follows:- Wettability depends on the initial oil saturation and hence the heights above the oil-water contact. It affects the oil relative permeabilities AND water relative permeabilities- The hysteresis observed on water relative permeability is even more pronounced that the initial oil saturation increases. A similar observation is made on the curvatures of the oil relative permeability.- The evolution of Sorw as function of Soi does not always follow a classical Land correlation. Their evolution also depends on the pore structure.- The two rock types exhibit different responses to wettability alteration. The vuggy structure might promote oil-wetness.Our experimental data were interpreted using different hysteretic models of the literature. We observed significant differences between the results predicted by the models and our experimental observations which are attributed mainly to the failures to take into account the change in wettability. We proposed an improved hysteresis model, which described our experimental data satisfactorily.