La compréhension des écoulements biphasiques dans des milieux hétérogènes est déterminante pour un certain nombre de processus industriels. Le plus marquant étant la récupération assistée du pétrole piégé dans des roches poreuses,qui est un moteur pour la recherche dans ce domaine depuis plusieurs décennies.Pourtant, si l’imbibition spontanée de liquides dans un milieu poreux est relativement bien comprises, l’imbibition forcée d’un liquide dans un milieu poreux rempli d’un liquide plus visqueux n’a pas encore été décrite quantitativement. En combinant des expériences de microfluidique et d’imagerie confocale, nous avons étudié la morphologie et la dynamique d’imbibition forcée dans des milieux poreux réguliers. Nous avons mis en évidence la présence de quatre régimes d’imbibition,provenant de dynamiques d’invasion différentes à l’échelle du pore et ayant une signature claire sur les observables macroscopiques. L’étude de la dynamique locale associée à la microscopie confocale, qui donne accès à des informations tridimensionnelles,nous a permis de montrer que les transitions entre ces régimes étaient la conséquence de deux transitions de mouillage et d’une instabilité d’interface.De plus, contrairement à la majorité des études menées jusqu’à présent, nous nous sommes intéressés à l’évolution de ces motifs d’imbibition soumis à une injection continue sur des temps longs devant le temps de percolation. Cela a permis de mettre en évidence le vieillissement de ces motifs qui est le résultat d’un mûrissement d’origine capillaire. / Understanding two-phase flow in heterogeneous media is of great importance for a number of industrial processes. One of the most prominent examples is enhanced oil recovery which has driven fundamental and applied research in this fieldfor decades. However our understanding has remained extremely unbalanced. The case of spontaneous imbibition of a fluid in a porous medium is fairly well understood,whereas quantitative descriptions of forced imbibition in the presence of anunfavorable viscosity ratio is still lacking. Combining large-scale observations and confocal imaging of microfluidic experiments, we studied the morphology and dynamicsof forced imbibition in homogeneous porous media. We identify four classesof three-dimensional patterns resulting from different dynamics at the pore-scale, and having a clear signature on the macroscopic observables. By means of confocalmicroscopy allowing us to visualize three dimensional features of the local dynamics,we show that the transitions between the four imbibition scenarios result from two dynamical wetting transitions and one interfacial instability. In addition,unlike previous studies, we investigate the evolution of those patterns undercontinuous injection over long time scales. We evidence their aging according to acapillary-coarsening process.

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).

Lors de sa détente dans une turbine, la vapeur subit une chute d’enthalpie qui entraîne sa condensation spontanée sous forme d’un nuage de gouttelettes submicroniques. Ces gouttes vont se déposer sur les aubes aval et y former un film d’eau. Ce dernier est arraché sous l’effet de l’écoulement de vapeur environnant ce qui crée des gouttes de quelques dizaines de microns qui peuvent se redéposer sur les aubes aval. Ces phénomènes sont à l’origine de pertes, généralement regroupées sous le terme générique de ”pertes par humidité”, estimées grâce à la loi de Baumann. Le but de cette thèse est de mettre en place dans un code méridien des modèles permettant la prise en compte des phénomènes de condensation, de déposition et d’écoulement des films liquides afin de pouvoir estimer les pertes liées à chacun d’entre eux. Dans cette optique un modèle diphasique homogène permettant d’avoir accès à la fraction massique de liquide et au nombre de gouttes est implanté dans le code méridien. Ce modèle est validé sur un cas test expérimental de détente en tuyère et alimente le modèle de déposition. Les contributions de la diffusion, de la turbophorèse, de la thermophorèse, de la gravité et de l’inertie des gouttes à la déposition sont étudiées. Il apparaît nécessaire de prendre en compte la diffusion, la turbophorèse et la déposition inertielle sur les bords d’attaque des aubes. Un modèle permettant d’avoir accès à l’épaisseur et à la vitesse d’un film liquide soumis au cisaillement d’un écoulement environnant, à la force de frottement sur la paroi et aux effets de la rotation est mis en place. Ce modèle est validé par rapport à des résultats expérimentaux d’écoulement de film liquide sur une plaque plane dans des conditions proches de celles rencontrées en turbine à vapeur basse pression. Finalement, un calcul réalisé sur une géométrie réelle de turbine basse pression de 8 étages permet de démontrer l’applicabilité de la méthodologie mise en place sur un cas industriel. Les contributions des phénomènes étudiés aux pertes par humidité sont explicitées. / During the expansion in a turbine, the enthalpy of the steam fall. This fall leads the steam to cross the saturation line and brings about its spontaneous condensation and the appearance of fog droplets. The deposition of these droplets on downstream blades lies at the root of the creation of a liquid film. This film is torn off by steam flow and creates coarse water. These big droplets can impact downstream blades. Each of these phenomena induces some losses which are generally grouped in the ”wetness losses” estimated thanks to Baumann’s rule. The aim of this work is to develop, in a throughflow code, some models which allow to take into account the condensation, deposition and liquid film flow and the losses linked to each of these phenomena. An homogeneous two phases flow model is implemented in the throughflow code. Two transport equations on the mass fraction of liquid and on the droplets number are added to Euler’s equations written for gas phase. This model is validated on an experimental test case of expansion in nozzle and feeds the model of deposition. The influences of diffusion, turbophoresis, thermophoresis, gravity and of inertia of droplets on the deposition are studied. It appears to be necessary to take into account diffusion, turbophoresisand deposition due to inertia on the leading edges of blades. A model which allows to evaluate the thickness and the velocity of a liquid film submitted to aerodynamic shear, friction on blades and rotational effects has been developed. This model is validated on experimental results of a liquid film flow on a flat plate in some conditions representatives of those encountered in low pressure steam turbines. Finally a calculation realised on a real geometry of a 8-stages low pressure steam turbine has demonstrated that the methodology developed during this work can be used on an industrial test case. The contributions of the phenomena studied to wetness losses are evaluated.

Des études internes à Liebherr-Aerospace France, concernant la conception de nouveaux compresseurs lubrifiés par gaz réfrigérant, ont montré que dans des conditions de fonctionnement spécifiques, un mélange de vapeur et de liquide apparaît au sein du compresseur, au lieu d'une phase vapeur seule. De ce fait, le comportement des paliers à feuilles lubrifiés au gaz réfrigérant est étudié, y compris la possibilité d'un écoulement diphasique du lubrifiant. L'étude porte sur le comportement du lubrifiant uniquement, dans des conditions de fonctionnement qui sont celles des paliers à feuilles. L'approche Thermo-Hydrodynamique décrit les caractéristiques du lubrifiant telles que la pression, la densité, la viscosité et la température. Dans ce modèle, une équation de Reynolds généralisée pour écoulement turbulent, une équation d'état non-linéaire pour écoulement diphasique et une équation de l'énergie tridimensionnelle pour film-mince et écoulement turbulent sont utilisées. Les paramètres globaux du palier sont calculés en régime permanent. / Internal experiments at Liebherr-Aerospace FRANCE on new refrigerant-lubricated compressor designs have shown that under specific operating conditions, a mixture of vapor and liquid appears in the compressor, instead of a single-phase vapor flow. Therefore, refrigerant-lubricated foil bearings behavior is studied, including the likelihood of two-phase flow in the lubricant. We focus on the lubricant behavior only, in the operating conditions of foil bearings. The Thermo-Hydrodynamic approach describes lubricant characteristics such as pressure, density, viscosity, and temperature. It involves the use of a generalized Reynolds equation for turbulent flow, a nonlinear cubic equation of state for two-phase flow and a 3D turbulent thin-film energy equation. Journal bearing global parameters are calculated for steady-state conditions.

This thesis addresses the experimental investigation of the cavitation phenomenon and its main consequences on the normal operation of a safety relief valve (SRV). More particularly, limitation of the mass flux discharged and alteration of the hydraulic fluid forces behavior is of main interest for the proper design and sizing of such devices. In nuclear or thermal engineering systems, the use of SRVs is mandatory since it represents the ultimate protection device before an accident occurs, caused by a sudden pressurization of the system. A careful design and sizing of the SRV is therefore essential. The complete understanding of the physics taking place in the flow through the valve is required to guaranty and optimize the security of the protected process.<p><p>In order to investigate the above effects of cavitation in a SRV, two different orifice sized valves (API 2J3 type and a transparent model based on an API 1 1/2G3 type) are tested in two different experimental facilities expressly built for this purpose. Instead of using a spring, the design of both valves allows the adjustment of the disc at any desired lift. Hence the static behavior of the valves is investigated. Both facilities, operating at different magnitude scales, allow the study of single phase and cavitating flow conditions required to properly determine the most important hydraulic characteristics, and access on any potential scaling effect between both sized SRVs. Experimental techniques used for the determination of the hydraulic characteristics include temperature, flow rate, fluid forces and pressure measurements both upstream and downstream the test sections. <p><p>Results show a similar influence of cavitation on the flow characteristics of both valves, minimizing any potential scaling effect. The liquid pressure recovery factor FL, which is normally used to identify a choked flow condition in a control valve, is experimentally determined for the first time in a SRV. The existence of a local minimum located at small openings of the lift indicates a change on the flow characteristics of both valves, which is related to the location of the minimum cross section of the flow that does not remain constant for every lift position. An extended experimental campaign is performed to analyse the effect of the blowdown ring adjustment located around the nozzle of the API 2J3 valve. Results confirm that the position of the ring has an important contribution for the hydraulic forces acting on the valve disc. <p><p>In the second part of the research, precise optical diagnostic techniques are successfully applied in the transparent valve to locally characterize the flow topology in a SRV experiencing cavitation. These results are innovative and enrich the experimental database available in the literature for the characterization and understanding of the flow physics in such devices. In a first configuration, high speed visualization is applied to observe qualitatively the flow pattern and the inception of liquid vaporization. Particle tracking results suggest that vapor bubbles are formed in the core of vortices detached from the shear layers attached to the valve. These rotational structures promote lower pressure regions allowing the liquid to vaporize. In the second configuration, particle image velocimetry is applied to extract the velocity field in both single phase and cavitating flow conditions. Results of PIV confirm the existence of a submerged jet just downstream the minimum section. This jet is characterized by two non-symmetric shear layers at its sides. Under cavitation conditions, PIV results confirm that vapor bubbles are formed preferentially inside the jet shear layers. The phenomenon of mass flux limitation caused by cavitation is reproduced at small openings of the valve and interaction with the flow topology is highlighted. It is observed that limitation of the flow occurs when the vena contracta is shifted towards the minimum geometrical section of the flow. Finally, instabilities of the flow downstream the critical section are investigated in the frequency domain by means of time resolved data. Results suggest that vortex shedding mechanism is dominated by a constant Strouhal number which is slightly affected by the valve opening. <p><p>In the last part of the research, the methodology used in water is extended and applied to cryogenic liquids. Two different geometries are investigated experimentally and numerically using water and liquid nitrogen as working fluids. Results suggest that both the flow coefficient (determined at single flow conditions), and the liquid recovery factor (used to identify choked flows), are independent on the fluid properties and therefore, an hydraulic similarity relation can be proposed.<p><p>This research project was carried out at the von Karman Institute for Fluid Dynamics (VKI), in Belgium, in close collaboration and with the funding of Centre Technique des Industries Mécaniques (CETIM) in France. / Doctorat en Sciences de l'ingénieur / info:eu-repo/semantics/nonPublished

The subject of this PhD thesis is the development of an efficient system that can simultaneously pump and separate a gas-liquid mixture, in particular an oil-air mixture. Two-phase flows are encountered in many applications (petroleum extraction, flow in nuclear power plant pumps, pulp and paper processing, etc.) but this study is mainly focused on lubrication systems of aircraft gas turbine engines.<p><p>The pump and separator system (PASS) for two-phase flows developed in this PhD thesis aims to perform three functions simultaneously:<p>• Send back the oil to the tank (oil pumping)<p>• Separate the air from the oil (de-aeration)<p>• Separate the oil from the air (de-oiling) and release the sealing air into the atmosphere (venting). <p>Particular care is given to the liquid flow rate lost at the gas outlet of the system.<p>Consequently, it could replace the scavenge pumps and oil-air separators existing in present lubrication systems. This modification provides several advantages: simplification of the lubrication circuit, reduction of oil consumption and of the size of the lubrication system.<p><p>This research is divided into three axes: the theoretical study of the important physical mechanisms taking place inside the two-phase flow pump and separator system, the experimental development, tests and optimization of different PASS prototypes, and also the numerical simulations of the two-phase flow inside these prototypes. Although the experiments were the central pillar of this research, the three axes were closely imbricated.<p><p>The PASS design includes three main components:<p>• An inlet chamber with one or several tangential inlets giving a natural centrifugation to the flow,<p>• An impeller (forced centrifugation) with an axial and a radial part followed by a volute chamber,<p>• A metallic foam that lets pass micron and sub-micron droplets and which is followed by an axial vent port.<p><p>The centrifugation causes the liquid (oil) to move radially outwards in an annular body (a liquid ring) generating pressure. The thickness of this liquid ring inside the impeller is mainly determined by the pressure coefficient (related to the back-pressure and the rotational speed). When the back-pressure increases, the thickness of the liquid ring increases too. An advantage of the PASS is that it does not impose any relation between the liquid head and the liquid flow rate, contrary to common centrifugal pump. It self-regulates the radial position of the gas-liquid interface to sustain the operating conditions.<p><p>The de-aeration efficiency mainly depends on the pressure coefficient (for a constant liquid viscosity or temperature) or on the thickness of the liquid ring. The pressure gradient which appears in the liquid rotating in an annular body acts like a dam for the gas phase. Indeed, the gas movement is mainly determined by the pressure field (buoyancy) while the liquid distribution is dominated by centrifugal and Coriolis forces. Buoyancy tends to accumulate the gas phase near low pressure areas (PASS hub, suction side of the blades, clearances between closed impeller and casing).<p><p>The first oil-air PASS prototype produces high viscous losses due to the high peripheral velocity and liquid viscosity. Therefore, the pumping efficiency is poor compared to common impeller pumps. However, the pumping is not the key function of the PASS and a power consumption below 5 kW is acceptable for the application considered in this work. For applications that require lower power consumptions, a reduction of the rotational speed must be considered.<p><p>Thus, the rotational speed and the impeller diameter are two major constraints for the PASS design which determine the de-aeration and pumping efficiencies. The impeller diameter also influences the size of passage sections for the air flow. The air velocity must be kept as low as possible because the entrainment of droplets increases when the air velocity rises (drag forces on droplets). Indeed, this large influence of the air flow rate on the oil consumption (de-oiling efficiency) was demonstrated by a theoretical analysis, the experiments and the CFD simulations. The production of droplets in the inlet pipes when the two-phase flow is annular is a key phenomenon regarding the oil consumption.<p><p>In addition to the air flow rate, other variables also influence the oil consumption:<p>• Air-oil temperature: when the temperature rises, the oil consumption increases because the surface tension and the oil density are reduced. Moreover, as the air density also decreases, the air velocity rises.<p>• Oil flow rate: the oil consumption rises more or less linearly with the oil flow rate. However, the influence of the oil flow rate on the inlet droplet size is uncertain.<p>• Rotational speed: the rotational speed has obviously a strong impact on the oil consumption without metallic foam. However, experiments showed that the metallic foam efficiency is almost independent on the rotational speed. Therefore, the oil consumption with the Retimet foam does not depend on the PASS rotational speed.<p>• Altitude or air density: the oil consumption decreases when the air density is reduced because the drag forces on droplets also decrease.<p>The gas density (altitude) is also supposed to influence the de-aeration efficiency but this could not be tested or simulated in this work (the de-aeration efficiency gets probably better when decreasing the gas density because the buoyancy forces increase).<p><p>Theory, experiments and numerical simulations also allowed the prediction of performance of the first oil-air prototype for real in-flight operating conditions. Two problems have been identified: the de-aeration efficiency at MTO and cruise ratings and the oil leak throughout the vent in cold start and windmilling. To solve them, some modifications of the lubrication system have been suggested. With these modifications, the oil-air PASS should become very efficient and attractive for engine manufacturers. / Doctorat en Sciences de l'ingénieur / info:eu-repo/semantics/nonPublished

This research project was carried out at the von Karman Institute for Fluid Dynamics (VKI), in Belgium, in collaboration and with the funding of Centre Technique des Industries Mécaniques (CETIM) in France.<p>The flow of a mixture of two fluids in pipes can be frequently encountered in nuclear, chemical or mechanical engineering, where gas-liquid eactors, boilers, condensers, evaporators and combustion systems can be used. The presence of section changes or more generally geometrical singularities in pipes may affect significantly the behavior of twophase flow and subsequently the resulting pressure drop and mass flow rate. Therefore, it is an important subject of investigation in particular when the application concerns industrial safety valves.<p>This thesis is intended to provide a thorough research on two-phase (air-water) flow phenomena under various circumstances. The project is split in the following steps. At first, experiments are carried out in simple geometries such as smooth and sudden divergence and convergence singularities. Two experimental facilities are built; one in smaller scale in von Karman Institute and one in larger scale in CETIM. During the first part of the study, relatively simple geometrical discontinuities are investigated. The characterization and modeling of contraction and expansion nozzles (sudden and smooth change of section) is carried out. The pressure evolution is measured and pressure drop correlations are deduced. Flow visualization is also performed with a high-speed camera; the different flow patterns are identified and flow regime maps are established for a specific configuration.<p>A dual optical probe is used to determine the void fraction, bubble size and velocity upstream and downstream the singularities.<p>In the second part of the project, a more complex device, i.e. a Safety Relief Valve (SRV), mainly used in nuclear and chemistry industry, is thoroughly studied. A transparent model of a specific type of safety valve (1 1/2" G 3") is built and investigated in terms of pressure evolution. Additionally, flow rate measurements for several volumetric qualities and valve openings are carried out for air, water and two-phase mixtures. Full optical access allowed identification of the structure of the flow. The results are compared with measurements performed at the original industrial valve. Flowforce analysis is performed revealing that compressible and incompressible flowforces in SRV are inversed above a certain value of valve lift. This value varies with critical pressure ratio, therefore is directly linked to the position at which chocked flow occurs during air valve operation. In two-phase flow, for volumetric quality of air=20%, pure compressible flow behavior, in terms of flowforce, is remarked at full lift. Numerical simulations with commercial CFD code are carried out for air and water in axisymmetric 2D model of the valve in order to verify experimental findings.<p>The subject of modeling the discharge through a throttling device in two-phase flow is an important industrial problem. The proper design and sizing of this apparatus is a crucial issue which would prevent its wrong function or accidental operation failure that could cause a hazardous situation. So far reliability of existing models predicting the pressure drop and flow discharge in two-phase flow through the valve for various flow conditions is questionable. Nowadays, a common practice is widely adopted (standard ISO 4126-10 (2010), API RP 520 (2000)); the Homogeneous Equilibrium Method with the so-called !-method, although it still needs further validation. Additionally, based on !-methodology, Homogeneous Non-Equilibrium model has been proposed by Diener and Schmidt (2004) (HNE-DS), introducing a boiling delay coefficient. The accuracy of the aforementioned models is checked against experimental data both for transparent model and industrial SRV. The HNE-DS methodology is proved to be the most precise among the others. Finally, after application of HNE-DS method for air-water flow with cavitation, it is concluded that the behavior of flashing liquid is simulated in such case. Hence, for the specific tested conditions, this type of flow can be modeled with modified method of Diener and Schmidt (CF-HNE-DS) although further validation of this observation is required. / Doctorat en Sciences de l'ingénieur / info:eu-repo/semantics/nonPublished

Cette étude porte sur l'utilisation conjointe (appelée ” couplage ”) de modèle de turbulence à grandes échelles LES (Large Eddy Simulation) et du modèle multiphasique VOF (Volume of Fluid). Cette utilisation conjointe est nécessaire dans de nombreuses applications industrielles comme celles de l'automobile où l'on recherche par exemple à évaluer les prestations diphasiques de l'auvent liées au phénomène d'entraînement et d'ingestion d'une nappe d'eau par le HVAC (système d'air conditionné). Cependant, l'utilisation conjointe de ces méthodes nécessite un traitement particulier de la turbulence proche de la surface liquide afin de reproduire convenablement la quantité de mouvement transmise depuis la phase gazeuse, motrice, jusque dans la phase liquide.<p><p>Basée sur une approche numérique et expérimentale, notre étude est articulée autour de trois axes. Tout d'abord, le développement de techniques de mesures spécifiques pour l'étude expérimentale de notre problématique: le LeDaR pour mesurer les déformées d'une interface et la PIV d'interface afin d'accéder aux champs de vitesse et de turbulence dans chacune des deux phases. Le second axe est la constitution d'une base de données expérimentales sur une configuration de type jet impactant sur une surface liquide représentative des phénomènes rencontrés dans l'auvent. Enfin, le troisième axe de travail est l'évaluation des modèles existants dans le code Ansys Fluent et à partir de cette analyse le développement et la validation de modèles de couplage LES-VOF.<p><p>L'évaluation des modèles développés a permis de valider une stratégie de calcul adaptée aux simulations de l'entraînement d'une surface d'eau par un écoulement d'air turbulent. / Doctorat en Sciences de l'ingénieur / info:eu-repo/semantics/nonPublished

Ce travail de thèse porte sur la modélisation et la simulation numérique de la collision et l'agglomération de particules colloïdales dans un écoulement fluide turbulent par une nouvelle méthode. Ces particules sont sensibles dans une même mesure aux effets brownien et turbulent. La première partie du travail concerne la modélisation du phénomène physique,allant du transport des particules jusqu'à la modélisation des forces d'adhésion physico-chimiques en passant par l'étape cruciale qui est la détection des interactions entre les particules (collisions). Cette détection des collisions est dans un premier temps étudiée par rapport aux algorithmes classiques existants dans la littérature. Bien que très efficaces dans le cadre de particules soumises à l'agitation turbulente, les conclusions de cette partie exposent les limites des méthodes existantes en termes de coûts numériques, pour le traitement d'un ensemble de colloïdes soumis au mouvement brownien. La seconde partie du travail oriente alors les travaux vers une vision novatrice du phénomène physique considéré. Le caractère diffusif aléatoire est alors considéré d'un point de vu stochastique, comme un processus conditionné dans l'espace et dans le temps. Ainsi, une nouvelle méthode de détection et de traitement des collisions de particules soumises exclusivement à un mouvement diffusif est présentée et validée, exposant un gain considérable en termes de coûts numériques. Le potentiel de cette nouvelle approche est validé et ouvre de nombreuses pistes de réflexion dans l'utilisation des méthodes stochastiques appliqués à la représentation de la physique / Ph.D thesis focuses on modeling and numerical simulation of collision and agglomeration of colloidal particles in a turbulent flow by using a new method. These particles are affected by both Brownian and turbulent effects. The first part of the work deals with current models of the physical phenomenon, from the transport of single particles to a model for physico-chemical adhesive forces, and points out the critical step which is the detection of interactions between particles (collisions). This detection is initially studied by applying classical algorithms existing in the literature. Although they are very efficient in the context of particles subject to turbulent agitation, first conclusions show the limitations of these existing methods in terms of numerical costs, considering the treatment of colloids subject to the Brownian motion. The second part of this work proposes a new vision of the physical phenomenon focusing on the random diffusive behaviour. This issue is adressed from a stochastic point of view as a process conditionned in space and time. Thus, a new method for the detection and treatment of collisions is presented and validated, which represents considerable gain in terms of numerical cost. The potential of this new approach is validated and opens new opportunities for the use of stochastic methods applied to the representation of physics

Les accidents graves résultant de la fusion d’un coeur de réacteur nucléaire doivent être anticipés pour améliorer l’efficacité de leur mitigation. De tels accidents sont survenus à TMI-2 (1979) et à Fukushima (2011). Suite à un accident de perte de refroidissement, l’échauffement du coeur et l’oxydation de la gaine de combustible suivie d’un renoyage (injection d’eau) peuvent entraîner l’effondrement des barres de combustible et la formation d’un lit de débris dans le coeur. La vapeur produite lors du renoyage peut activer l’oxydation exothermique du Zircaloy, entraînant la fusion partielle des matériaux. Cette évolution engendre des zones à porosité réduite limitant la pénétration de l’eau et/ou des zones imperméables. Dans cette situation, l’efficacité de l’injection d’eau dans le coeur pour arrêter la progression de la dégradation et empêcher la fusion du coeur du réacteur peut être considérablement réduite. Dans ce cadre, l’IRSN a lancé le programme PEARL visant à étudier la thermohydraulique du renoyage des lits de débris chauds entourés d’une zone plus perméable simulant la présence de zones intactes ou moins endommagées dans le coeur. Dans cette thèse, les expériences PEARL ont été modélisées et simulées avec ICARE/CATHARE pour évaluer l’évolution d’un renoyage d’un lit de débris surchauffé entouré d’un bypass de perméabilité plus grande. La thermohydraulique du processus a été analysée et l’effet de différents paramètres (géométrie, conditions aux limites) sur le comportement de renoyage a été évalué. Sous certaines conditions, l’entraînement de l’eau dans le bypass a été identifié et évalué. Un modèle analytique a été mis au point ensuite pour étudier de façon approfondie le renoyage d’un milieu poreux hétérogène surchauffé composé de deux lits de débris de perméabilité et de porosité différentes et pour décrire l’entraînement de l’eau dans le bypass. Ce modèle calcule les principales variables caractérisant le processus de renoyage, telles que la vitesse du front de trempe, le taux de conversion eau-vapeur et le débit d’eau entraîné dans le bypass.Il fournit de bons résultats qualitatifs et quantitatifs concernant la redistribution du débit d’eau par rapport aux résultats expérimentaux. Ce modèle a plusieurs avantages. Il est écrit sous une forme plutôt générale incluant les termes de correction de Forchheimer et les termes croisés non nuls dans l’équation de Darcy-Forchheimer généralisée. Les différentes options des équations de quantité de mouvement proposées, y compris les changements dans les corrélations et les lois de frottement interfacial, peuvent être testées facilement. La comparaison des calculs avec les résultats expérimentaux indique qu’il est nécessaire d’inclure une loi de frottement interfacial pour obtenir de bonnes prédictions. L’extrapolation à l’échelle du réacteur est simple et des calculs ont été effectués pour évaluer l’impact des paramètres géométriques du lit de débris (granulométrie, porosité, dimensions) ainsi que les conditions thermiques et hydrauliques (température, pression, débit d’injection). Ainsi, le modèle est très utile pour estimer le temps de trempe total et latempérature maximale qui pourraient être atteinte dans le lit de débris à grande échelle. Cela permet d’évaluer la probabilité de réussite du renoyage d’un lit de débris chauds formé lors d’un scénario accidentel hypothétique. / Severe accidents arising from the fusion of a nuclear reactor core must be anticipated to enhance the efficiency of their mitigation. Such accidents have occurred at TMI-2 (1979) and Fukushima (2011). Following a loss of coolant accident, core heating and oxidation of the fuel cladding followed by reflooding (injection of water) may lead to the collapse of fuel rods and formation of porous debris bed in the core. Steam produced upon reflooding may activate the exothermic oxidation of Zircaloy leading to partial melting of materials. Such evolution generates zones with reduced porosity limiting coolant penetration and/or impermeable blocked zones. In this situation, the efficiency of injecting water into the core to stop the progress of degradation and prevent the reactor core melting may be significantly reduced. In this scope, IRSN launched PEARL program to investigate the thermal hydraulics of reflooding of hot debris beds surrounded by a more permeable zone simulating the presence of intact or less damaged zones in the core. The PEARL experiments were modeled and simulated using ICARE/CATHARE code to assess the evolution of a bottom reflooding of a superheated debris bed surrounded by a bypass of larger permeability. The thermal hydraulics of the quenching process has been analyzed and the effect of each of the initial conditions on the reflooding behavior was assessed. The effect of pressure was investigated and related to the entrainment of injected water at quench front level into the bypass. An analytical model was then developed to investigate thoroughly the reflooding of a superheated heterogeneous porous medium, composed of two layers of contrasting permeability and porosity, and to describe the water entrainment in the bypass. This model computes the main variables characterizing the reflooding process such as quench front velocity, water-to-steam conversion ratio, and the flow rate of water entrained in the bypass. It provides good qualitative and quantitative results for the two-phase flow redistribution as compared to experimental results. This model has several advantages. It is written in a rather general form including the Forchheimer correction terms and non-zero cross-terms in the generalized Darcy-Forchheimer momentum equation. Variations of proposed momentum equations including changes in correlations andinterfacial friction laws can be tested easily and efficiently. Comparison of the calculations against experimental results indicated that it is necessary to include an interfacial friction law to obtain good predictions. This model allows performing fast evaluations of the efficiency of cooling bycomputing the fraction of the injected flow rate that participates in cooling. Upscaling to the reactor scale is straightforward and calculations were performed to assess the impact of geometric parameters of the debris bed (particle size, porosity, dimensions) as well as thermal hydraulic conditions (temperature, pressure, injection flow rate) on the reflooding process. Thus the model is very useful to estimate the total quenching time and the maximum temperature that could be reached by the hot debris bed at large scales. This allows assessing the probability of a successful quenching of a hot debris bed formed during a hypothetical accidental scenario.