Dû aux efforts apportés à la réduction des émissions de NOx dans des chambres de combustion aéronautiques il y a une tendance récente vers des systèmes à combustion pauvre. Cela résulte dans l'apparition de nouveaux types d'injecteur qui sont caractérisés par une complexité géométrique accrue et par des nouvelles stratégies pour l'injection du carburant liquide, comme des systèmes multi-point. Les deux éléments créent des exigences supplémentaires pour des outils de simulation numériques. La simulation à grandes échelles (SGE ou LES en anglais) est aujourd'hui considérée comme la méthode la plus prometteuse pour capturer des phénomènes d'écoulement complexes qui apparaissent dans une telle application. Dans le présent travail, deux sujets principaux sont abordés: Le premier est le traitement de la paroi ce qui nécessite une modélisation qui reste délicate en SGE, en particulier dans des géométries complexes. Une nouvelle méthode d'implémentation pour des lois de paroi est proposée. Une étude dans une géométrie réaliste démontre que la nouvelle formulation donne de meilleurs résultats comparé à l'implémentation classique. Ensuite, la capacité d'une approche SGE typique (utilisant des lois de paroi) de prédire la perte de charge dans une géométrie représentative est analysée et des sources d'erreur sont identifiées. Le deuxième sujet est la simulation du carburant liquide dans une chambre de combustion. Avec des méthodes Eulériennes et Lagrangiennes, deux approches sont disponibles pour cette tâche. La méthode Eulérienne considère un spray de gouttelettes comme un milieu continu pour lequel on peut écrire des équations de transport. Dans la formulation Lagrangienne, des gouttes individuelles sont suivies ce qui mène à des équations simples. D'autre part, sur le plan numérique, le grand nombre de gouttes à traiter peut s'avérer délicat. La comparaison des deux méthodes sous conditions identiques (solveur gazeux, modèles physiques) est un aspect central du présent travail. Les phénomènes les plus importants dans ce contexte sont l'évaporation ainsi que le problème d'injection d'un jet liquide dans un écoulement gazeux transverse ce qui correspond à une version simplifiée d'un système multi-point. Le cas d'application final est la configuration d'un seul injecteur aéronautique, monté dans un banc d'essai expérimental. Ceci permet d'appliquer de manière simultanée tous les développements préliminaires de ce travail. L'écoulement considéré est non-réactif mais à part cela il correspond au régime ralenti d'un moteur d'avion. Dû aux conditions préchauffées, le spray issu du sstème d'injection multi-point s'évapore dans la chambre. Cet écoulement est simulé, utilisant les approaches Eulériennes et Lagrangiennes et les résultats sont comparés aux données expérimentales.

Les systèmes fixes par aspersion d'eau sont souvent cités comme un moyen pour améliorer le niveau de sécurité dans les tunnels routiers. Cependant, leurs effets sur un incendie en tunnel et les interactions possibles avec le système de ventilation sont encore mal connus. Le travail présenté ici vise à étudier l'action de ces systèmes par une approche phénoménologique de l'évaporation et l'utilisation d'un modèle CFD tridimensionnel. L'emploi d'un code CFD incluant des modèles diphasiques repose sur la simulation de deux situations expérimentales : une expérience académique de spray évaporant et des essais d'incendie avec aspersion dans un tunnel à échelle réduite. Les résultats obtenus sont relativement satisfaisants, en utilisant une modélisation Euler-Lagrange pour l'écoulement diphasique. L'influence du système d'aspersion sur la propagation de la couche de fumée est plus particulièrement étudiée. Cependant, ce type de modélisation nécessite des données d'entrée d'un niveau difficilement accessible sans avoir recours à une expérience de caractérisation des buses. Une approche plus phénoménologique, à intégrer dans un modèle unidimensionnel, est envisagée comme alternative.

Dans ce mémoire, on s'intéresse à la simulation des écoulements liquide-vapeur en transition de phase. Pour décrire ces écoulements, une approche bifluide moyennée à deux pressions indépendantes est retenue. Cette description du mélange liquide-vapeur s'appuie sur le modèle à sept équations de Baer et Nunziato. On étudie les aptitudes de cette modélisation à simuler les transitions de phase apparaissant en ingénierie nucléaire.<br /><br />Dans un premier temps, on élabore un cadre thermodynamique théorique pour décrire les écoulements liquide-vapeur. Dans ce cadre, on réalise la fermeture du modèle de Baer et Nunziato. De nouvelles modélisations sont proposées pour les termes d'interaction entre les phases. Ces nouvelles modélisations dotent le modèle bifluide à deux pressions d'une inégalité d'entropie. On étudie ensuite les propriétés mathématiques de ce modèle. Sa partie convective hyperbolique se présente sous une forme non-conservative. On étudie tout d'abord la définition de ses solutions faibles. Divers régimes d'écoulement sont alors mis à jour pour le mélange diphasique. Ces différents régimes d'écoulement présentent des analogies avec le comportement fluvial et torrentiel des écoulements en rivière. Les stabilités linéaire et non-linéaire de l'équilibre liquide-vapeur sont ensuite établies. Pour affiner notre description des interactions diphasiques, on étudie pour finir l'implémentation d'un modèle de turbulence, ainsi que l'implémentation d'une procédure de reconstruction pour la densité d'aire interfaciale.<br /><br />On s'intéresse ensuite à la simulation de ce modèle. Suivant une approche à pas fractionnaires, une méthode numérique est élaborée dans un formalisme Volumes Finis. Pour réaliser l'approximation de la partie convective, diverses adaptations non-conservatives de solveurs de Riemann standard sont tout d'abord proposées. A l'inverse du cadre non-conservatif classique, l'ensemble de ces schémas converge vers une unique solution. Un nouveau schéma de relaxation est ensuite proposé pour approcher la dynamique des transferts interfaciaux. L'ensemble de la méthode numérique se caractérise alors par la préservation des équilibres liquide-vapeur. Dans un premier temps, cette méthode numérique est employée à la comparaison des différentes modélisations bifluides à une et deux pressions. On l'applique ensuite à la simulation des écoulements liquide-vapeur dans les circuits hydrauliques des réacteurs à eau sous pression en configuration accidentelle.

En simulation de réservoir, les variations spatiales des propriétés des roches sont données par un modèle géologique défini sur une grille pouvant comporter plusieurs centaines de millions de mailles. Simuler des écoulements polyphasiques sur ce type de grille peut être très coûteux en temps de calculs. Les ingénieurs de gisement procèdent souvent à une mise à l'échelle de ces propriétés sur une grille plus grossière, utilisée ensuite pour la simulation. Cette façon de procéder ne permet pas toujours de bien prendre en compte l'influence sur l'écoulement des hétérogénéités présentes à l'échelle fine. Plusieurs méthodes multi-échelles ont donc été proposées ces dernières années afin d'obtenir un meilleur compromis entre les temps de calculs et la précision des solutions obtenues. Actuellement, la plupart de ces méthodes permettent de résoudre plus efficacement l'équation en pression. Les équations de transport sont, quant à elles, toujours résolues sur la grille fine. La première partie de cette thèse présente la mise en œuvre et les résultats de tests réalisés sur la méthode des éléments finis mixtes multi-échelles dans un environnement de calcul parallèle. Le second et principal objectif de ce travail est de proposer une nouvelle méthode multi-échelle pour la simulation du transport d'un traceur dans un milieu poreux. Cette méthode est construite à partir de résultats rigoureux d'homogénéisation en milieux périodiques. Cette thèse fournit une vitesse de convergence explicite pour ce procédé d'homogénéisation ainsi qu'une estimation d'erreur pour la méthode numérique multi-échelle associée. Des exemples numériques sont ensuite présentés.

Dans l'objectif de maîtriser le démarrage des moteurs fusées en vol balistique, il est proposé de caractériser le transitoire de remplissage des cavités d'injection des organes de combustion. L'étude s'effectue principalement sur la cavité tampon, appelée dôme. Le comburant est maintenu sous pression en amont d'une vanne dont l'ouverture contrôle son passage vers une chambre tampon qui est liée à la chambre de combustion à travers un réseau d'injecteurs. Afin d'empêcher la remontée du carburant vers la chambre tampon, un gaz balaye la chambre de l'entrée vers les injecteurs. Cette étude expérimentale consiste à décrire la structure spatio-temporelle de l'écoulement diphasique dans la cavité tampon suite à l'ouverture de la vanne en présence de l'écoulement du gaz de balayage. Il s'agit de suivre l'évolution de l'écoulement sur des temps courts (quelques centaine ms) par un ensemble de mesures (débits, pressions, distribution spatiale des phases, suivi de l'interface) sans ou avec transfert thermique aux parois de la cavité. Des fluides de substitutions sont utilisés. Pour la partie expérimentale sans transfert thermique, de l'eau et de l'air sont utilisés à la place du comburant et du gaz inerte et pour la partie non isotherme du "x" (fluorocarbone) et un gaz "y" ont été choisis. Dans le premier cas, les expériences isothermes ont mis en évidence le comportement typique des grandeurs comme la pression dans la cavité et le débit de liquide entrant ainsi que la distribution des phases en sortie d'injection, tandis qu'une analyse des résultats a montré l'importance des différentes échelles de temps qui interviennent pendant le remplissage: temps d'ouverture de la vanne, temps de recouvrement des injecteurs par le liquide et les temps de remplissage et de vidange de la cavité tampon. Dans le deuxième cas, les parois sont chauffées au-dessus de la valeur d'ébullition du liquide, pour la gamme de pression envisagée dans le dôme. Le but est de quantifier les effets d'un possible changement de phase aux parois et d'évaluer leur importance sur l'écoulement. De plus, les conséquences dues à l'échauffement du gaz de balayage constituent une partie significative du programme expérimental. Le gaz est chauffé indépendamment des parois. Un modèle théorique traitant des différents régimes d'écoulements pendant le transitoire de remplissage permet de reproduire le comportement des pressions, débit liquide et fraction volumique de gaz dans la cavité. L'ensemble de ces mesures permettent de comprendre le transitoire de remplissage de la cavité d'injection dôme et l'analyse théorique qui accompagne ces expériences doit permettre l'extrapolation des résultats obtenus en laboratoire aux conditions réelles (fluides cryogéniques sous microgravité). Elle doit aussi fournir les conditions aux limites requises pour les approches numériques développées par ailleurs ainsi que les bases de données permettant de tester ces simulations. / In order to control the ignition of rocket motors during ballistic flight, the transient flow of comburant into a reservoir or buffer cavity (dôme) and then through a grid of injectors must be carefully characterised. The liquid oxygen is held under pressure upstream of a valve which opens into the dome. The valve opening is a control parameter. To avoid any possible flow of carburant from combustion chamber back into reservoir, the latter is swept with an inert gas, thus ensuring that the pressure in the reservoir remains higher than in the combustion chamber. This experimental study has the aim of characterising the spatio-temporal structure of two-phase flow into the dome following opening of principal liquid valve. Filling the dome and forcing the liquid through the injectors has an overall time scale of some hundred milliseconds. High resolution measurements of liquid and gas flow rates, pressure, phase distribution, interface velocity and temperatures are recorded for different values of the key parameters as well as visualisations. For the experimental program with heat transfer, the comburant was replaced with "x". Much work was carried out on the effects of heat transfer from either the gas or the walls or both to the liquid entering the dome once these were above the liquid boiling point. Previous studies in the LEGI using water and air, and without heat transfer brought to light the important variations in dome pressure and liquid flow rate during the transient, while analysis of results indicated the importance of a number of time scales : value opening time, time for the liquid to cover the injectors, time to fill the dome, time to empty it. For the heat transfer experiments, the walls are heated for the pressure range chosen. The sweep gas is heated too. The aim of these experiments is to seek evidence of a phase-change at the walls or during interaction with the gas and to evaluate its importance. To carry out these experiments, specific instrumentation was used. The whole of these mesures enable us to understand the transient filling of the injection cavity. Thus, theoretical analysis have to allow extrapolations of results obteined in laboratory to real cases (cryogenic liquid under microgravity). Then, we have to give a data base to developp and validate numerical simulation.

Cette thèse s’inscrit dans le cadre du projet ANR « Uncertain flow optimization » (UFO) consacré au développement et à l’application de méthodes efficaces de quantification d’incertitudes pour l’analyse et l’optimisation d’écoulements. Dans ce contexte, ces méthodes sont appliquées à des gazosiphons de grande échelle utilisés comme pompe. Plus particulièrement, on s’intéresse à de tels gazosiphons choisis pour constituer l’organe central d’un système géothermique innovant,  compatible avec un environnement urbain. On souhaite en quantifier le potentiel énergétique par voie numérique avec la recherche d’un compromis entre justesse des résultats et efficacité optimale. La simulation de l’écoulement diphasique produit dans le gazosiphon est fondée sur un modèle quasi-1D à flux de dérive et s’appuie sur une démarche de résolution implicite. Les résultats sont validés sur les études expérimentales les plus pertinents de la littérature, dont aucune toutefois n’atteint les longueurs requises de l'ordre du kilomètre. Le code de simulation du gazosiphon fait ensuite l’objet d’une démarche de prise en compte d’incertitudes physiques et de modélisation, précédée par une analyse de deux méthodes de quantification d’incertitude : une méthode non-intrusive de type chaos polynomial, et une méthode plus récente dite semi-intrusive qui fut développée en amont du projet UFO. Cet outil est intégré dans une modélisation simplifiée du système géothermique urbain dans son ensemble impliquant les composants en surface, notamment le compresseur d'air. Il en résulte une optimisation énergétique robuste préliminaire de deux variantes du système géothermique urbain proposé, respectivement de récupération de chaleur et de production d’électricité. / This PhD thesis is part of the ANR project « Uncertain Flow Optimization » (UFO). The project is devoted to the development and application of efficient uncertainty quantification methods for flow analysis and optimization. In this framework, these methods are applied to the study of a large-scale airlift pump. The airlift pump is selected to be part of an innovative geothermal system, which can be exploited within an urban environment. We wish to quantify and optimize the energy potential of this new system with numerical tools. They provide both good accuracy and efficiency properties. The airlift two-phase flow simulation is based on a quasi one-dimensional drift flux model, which is implicitly solved. The solver is validated by comparison with relevant experimental airlift studies from the literature. However, these studies remain below the kilometric-targeted pipe length. Thanks to the analysis of two uncertainty quantification methods, a non-intrusive approach relying on polynomial chaos expansion and a new semi-intrusive method developed ahead of the UFO project, we perform airlift pump simulations taking into account physical and modelling uncertainties. This numerical tool is inserted into a simplified model of the complete urban geothermal system that involves surface devices, such as an air compressor. Finally, a robust preliminary optimization process is performed for two versions of the proposed geothermal urban system. They are designed respectively for heat recovery and electricity production.

La production pariétale de bulles de gaz et son impact sur la dynamique de la phase liquide en canal vertical est étudiée numériquement et expérimentalement. Dans un premier temps, un modèle de mélange 2D stationnaire est utilisé pour décrire l’évolution moyenne des panaches de gaz. Grâce à cette approche, un modèle de couche limite a pu être développé et a permis l’identification des nombres adimensionnels pertinents (analogues aux nombres de Rayleigh et de Prandtl pour la thermique) afin de caractériser les écoulements à bulles dispersées. Dans un second temps, un modèle Eulérien-Lagrangien 3D instationnaire, prenant en compte le couplage quadrilatéral (interactions bullesliquide et bullesbulles) est résolu par Simulation Numérique Directe (DNS) et permet ainsi une description plus fine de l’écoulement à l’échelle de la phase dispersée. Enfin, ces approches numériques sont complétées par des mesures de Spectroscopie d’Impédance Electrochimique (SIE) lors de la production de dihydrogène et de dioxygène par électrolyse alcaline. Les modèles d’écoulement proposés ici montrent globalement un très bon accord avec les résultats expérimentaux tirés de la littérature. Les approches homogènes et DNS présentent toutefois quelques disparités sur l’évaluation du taux de vide dans certaines conditions. Parallèlement, les mesures et simulations de SIE ont montré être clairement affectées par les évolutions du panache de bulles, les spectres d’impédance ont notamment mis en évidence une contribution basse fréquence fortement dépendante de la nature de la phase dispersée (taille de bulle et lois de dispersion). Les trois approches (modèle homogène, DNS et SIE) menées conjointement sont donc fortement complémentaires. Elles permettent non seulement une meilleure compréhension de la physique de l’écoulement diphasique, mais offrent aussi une capacité d’analyse de la pertinence des modèles existants tout en ouvrant la voie à leurs futures améliorations / The wall production of gas bubbles and its impact on the liquid dynamics in a vertical channel is studied by means of numerical simulations and experimentation. First, a 2D stationary mixture model is used to describe the averaged plumes evolutions. Through this approach, a boundary layer model has been developed and identified dimensionless numbers (Raleigh-like and Prandtl-like) characteristic of bubbly flows. Secondly, a 3D non-stationary four-way coupled (with bubblesliquid and bubblesbubbles interactions) Eulerian-Lagrangian model is solved by Direct Numerical Simulation (DNS) and allows a finer description of the two-phase flows at bubble-scale. Finally, the numerical methods are completed by Electrochemical Impedance Spectroscopy (EIS) measurements during hydrogen and oxygen production by alkaline electrolysis.The two-phase flow models are in good agreement with experimental results from literature. There are still some disparities between the homogeneous model and the DNS about the void fraction calculation under certain conditions, though. In the meantime, both EIS measurements and simulations were clearly affected by bubbles plume evolutions, the impedance spectra highlighted a low frequency contribution highly sensible to the nature of the dispersed phase (bubble size a dispersion laws). The 3 approaches (homogeneous model, DNS and EIS) used collectively are strongly complementary. They allow not only a better comprehension of the physics of the two-phase flow, but also serve the analysis of existing models while leading the way for further improvements

Cette thèse traite principalement des connaissances fondamentales en hydrodynamique et des caractéristiques des réactions gaz-liquide dans des microréacteurs capillaires. Dans une première partie, nous avons effectué des essais dans trois microcanaux circulaires en verre placés horizontalement. Les diamètres étudiés étaient de 302, 496 et 916 µm. Les arrivées de gaz et de liquide se font de manière symétrique et forme un angle de 120° entre elles. Une cartographie des écoulements diphasiques gaz-liquide a été systématiquement faite pour des vitesses du liquide comprises entre 0,1 et. 2 m/s et des vitesses du gaz comprises entre 0,01 et 50 m/s Ces essais mettent en évidence l'influence du diamètre des canaux, de la viscosité du liquide et de leur tension superficielle. Ces mesures ont été comparées avec les cartes décrivant les différents régimes d'écoulement (à bulles, en bouchons de Taylor, annulaires ou sous forme de mousse) et confrontés aux modèles de la littérature qui prédisent les transitions entre les différents régimes. Nous avons mis en évidence que ces derniers n'étaient pas totalement satisfaisant et en conséquence, un nouveau modèle de transition prenant en compte les effets de taille du canal, les propriétés physiques du liquide a été proposé. Les pertes de charge engendrées par ces écoulements gaz- ont été étudiées. Nous avons constaté que la chute de pression est très dépendante du régime d'écoulement. Cependant pour décrire l'évolution de la perte de charge il est commode de la scinder en trois régions: une où les forces de tension superficielle sont le paramètre prépondérant et qui correspond aux faibles vitesses superficielle du gaz, une zone de transition et une dans laquelle les forces d'inertie sont dominantes et qui correspond aux grandes vitesses superficielles du gaz. La prédiction de cette chute de pression dans la troisième zone a été faite à partir d'un modèle de Lockhart-Martinelli. Ce modèle qui prend en compte les flux de chaque phase dépend d'un paramètre semi empirique C. Nous avons proposé de le corréler avec les nombres de Reynolds correspondant à chacune des deux phases en présence. Cette méthode permet de bien rendre compte de nos mesures. Les caractéristiques hydrodynamiques en écoulement de Taylor ont été examinées. Il a été montré que la formation des bulles dans un écoulement de Taylor est dominée par un mécanisme d'étranglement en entrée du capillaire. La taille des bulles dépend fortement de la viscosité du liquide et la tension superficielle. La chute de pression dans cette zone, lorsque le nombre capillaire est relativement faible, peut assez être bien décrite par le modèle de Kreutzer modifiée par Walsh et al… En fin dans une dernière partie, nous avons réalisé une réaction chimique en écoulement de Taylor. L'oxydation du 2-hydrogéne-ethyltetrahydroanthraquinone (THEAQH2) pour former du peroxyde d'hydrogène a été expérimentalement étudiée dans un microcanal circulaire horizontal de 900 µm de diamètre et 30 cm de long. La présence d'une réaction chimique ne modifie que très peu les transitions entre les différents régimes d'écoulement ni l'évolution des pertes de charge. Les cinétiques de conversion du peroxyde d'hydrogène sont environ deux fois plus rapides celles obtenues dans les réacteurs gaz liquide utilisés habituellement. Mots-clés: microcanal, écoulement diphasiques, écoulement de Taylor, pertes de charge, réaction gaz-liquide. / This dissertation mainly deals with the fundamental knowledge of hydrodynamics and reaction characteristics in gas-liquid microreactors. Extensive experimental investigations have been performed in horizontal circular microchannels with diameter from 302 µm to 916 µm. Gas-liquid two-phase flow patterns in the microchannel have been systematic experimental investigated, in which the influence of channel diameters, liquid viscosities and surface tension were considered. Flow pattern regime maps in the present microchannels were developed, and the comparison with existing regime maps and flow pattern transition models in literature implied that transitions in present work could not be well predicted. As a result, a new transition model taking the effects of channel size, liquid physical properties into account was proposed. The gas-liquid two-phase pressure drop characteristics in microchannels were studied. It has been found that the pressure drop was highly flow patterns dependent, and the main trend can be divided into three regions: surface tension-dominated region, transitional region and inertia-dominated region. The pressure drop characteristics in surface tension-dominated and inertia-dominated region were discussed respectively. A modified Lockhart-Martinelli separated flow model in which the effects of channel diameter and liquid properties on the C-value are taken into account was proposed, and it showed a good agreement with respect to our experimental data and others' reported in literature. Hydrodynamics characteristics of Taylor flow have been examined. It was shown that the formation of Taylor flow was dominated by squeezing mechanism, on which the effects of liquid viscosity and surface tension were dramatically. The two-phase pressure drop of Taylor flow could be well predicted with the Kreutzer's model modified by Walsh et al., when capillary number was relatively low. Oxidation of hydrogenated 2-ethyltetrahydroanthraquinone (THEAQH2) in a horizontal circular microchannel have been experimental investigated. Results of visualization study on oxygen-anthraquinone working solution two-phase flow in microchannel showed that the flow pattern transition model and pressure drop model for inertia-dominated region proposed in this dissertation had good predicting accuracy. It was indicated that the gas-liquid interfacial area and space-time yield of hydrogen peroxide in the microchannel are at least one to two orders of magnitude higher than those in the conventional gas-liquid reactors. Keywords: microchannel, two-phase flow pattern, pressure drop, gas-liquid reaction, Taylor flow..

L’objectif est d’améliorer la connaissance des écoulements de films minces se développant au sein des colonnes servant à distiller les gaz de l’air sous conditions cryogéniques. Ce travail porte sur la caractérisation expérimentale de l’épaisseur du film liquide dans des conditions hydrodynamiques proches de celles rencontrées dans l’industrie. Tout d’abord, ce travail a permis de déterminer les conditions opératoires permettant la réalisation d’écoulement en similitude hydrodynamique avec les écoulements cryogéniques mais à température et pression proches de l’ambiant. Ensuite, un dispositif expérimental a ensuite été développé afin de permettre la réalisation d’écoulement de film liquide avec et sans contre-courant gaz sur plan incliné. Une métrologie optique utilisant un procédé de fluorescence a été développée et calibrée afin de mesurer des épaisseurs de films minces inférieures à 2 mm. Pour finir, l’étude expérimentale a permis de caractériser l’écoulement du film mince par l’étude statistique de son épaisseur, l’étude des régimes d’ondes de surface ainsi que son aire interfaciale. / The aim of this study is to get a better knowledge of thin film flow inside column used to distil the gases of air under cryogenic conditions. Experimental characterisation of the liquid film thickness submitted to hydrodynamic conditions similar to those found in real processes is investigated in this work. First, operating conditions were defined that allow thin film flow in hydrodynamic similitude with cryogenic flows but under conditions close to standard temperature and pression. Then, an experimental setup was designed in order to realize liquid film flows down an inclined plate with and without counter-current air flow. An optical technique using fluorescence was built and calibrated in order to measure the thickness of liquid film up to 2 mm. Finally, thin liquid film flow characteristics were experimentally studied with special attention given to the statistical variation of film thickness, interfacial wave patterns and interfacial area.

Cette thèse est motivée par le besoin de compléter les connaissances actuelles des phénomènes ayant lieu lors d’un renvoi d’eau dans un lit de débris radioactifs, opération appelée « renoyage » et qui intervient dans une séquence d’accident grave où un cœur de réacteur nucléaire est dégradé suite à une perte prolongée de refroidissement primaire. Notre étude, de nature expérimentale, vise à mieux caractériser la crise d’ébullition en convection forcée, dans un milieu poreux chauffant. Le cœur du dispositif expérimental est un milieu poreux modèle quasibidimensionnel, composé de 276 cylindres disposés entre deux plaques de céramique distantes de 3 mm, dont l’une, transparente, permet de visualiser les écoulements. Les cylindres, de 2 mm de diamètre, sont des sondes thermo-résistives qui ont une double fonction : elles sont utilisées comme éléments chauffants et comme capteurs de température. Une boucle fluide permet de contrôler le débit d’injection de liquide dans la section test, la température d’injection ainsi que la pression. La section test est placée verticalement, le liquide est injecté par le bas à une température proche de la saturation. Dans une première série d’expériences, la puissance thermique dissipée globalement par un ensemble de cylindres chauffants est augmentée de façon progressive jusqu’à atteindre l’assèchement d’une zone du milieu poreux. Les résultats montrent deux types de phénoménologies dans le déclenchement de la crise d’ébullition. Pour des débits d’injection faibles (densités de flux massique de l’ordre de 4 kg.m^-2.s^-1 maximum), l’atteinte de la puissance d’assèchement se traduit par un lent recul du front diphasique jusqu’à sa stabilisation en haut de la zone chauffée ; en aval de la zone chauffée, l’écoulement est monophasique vapeur. Pour des débits d’injection plus élevés, la crise d’ébullition apparaît autour d’un des éléments chauffants, conduisant à une ébullition en film localisée, tandis qu’un écoulement diphasique liquide-vapeur continue de parcourir l’aval de la section test. Les visualisations de ces expériences permettent d’identifier qualitativement la structure des écoulements. D’autres expériences consistent à mesurer le flux critique local autour d’un cylindre choisi, pour différentes configurations d’écoulements. Le débit d’injection est fixé. Une puissance de chauffe est imposée à une ligne horizontale de cylindres en amont du cylindre choisi. Les résultats montrent que le flux critique sur ce cylindre diminue en fonction de la puissance délivrée à la ligne chauffée. La distance du cylindre étudié à la ligne chauffée semble avoir peu d’influence sur le flux critique. Des visualisations expérimentales sont utilisées pour caractériser l’écoulement diphasique en aval de la ligne chauffée, dans le but de mettre en relation le flux critique local avec des paramètres hydrodynamiques (saturations, vitesses des phases). Les images obtenues sont difficiles à exploiter. Afin de calibrer les paramètres des algorithmes de traitement d’images, nous avons reproduit une cellule d’essai de géométrie identique à l’originale, mais où l’on injecte du gaz par une ligne de cylindres en amont de la section test dans une configuration d’écoulement diphasique isotherme. Dans ce dispositif, le débit d’injection de gaz est contrôlé et mesuré. Les visualisations obtenues servent alors de références auxquelles sont comparées les visualisations d’ébullition convective. / This work is motivated by the need to better understand the phenomena occurring while some water is injected into a heated porous debris bed. This reflooding operation is a part of the planned mitigation procedure used during a Loss Of Coolant Accident (LOCA) that may occur into a nuclear power plant and results into a severe core damage. Our experimental study aims to characterize the boiling crisis that can happen in a boiling flow taking place within a heatgenerating model porous medium. The test section is a two-dimensional model porous medium, composed of an array of 276 cylinders placed between two ceramic plates spaced from one another by 3 mm, one of which is transparent and allows visualizations of the flow. The 2 mm diameter cylinders are Pt100 resistance temperature detectors that perform a dual function: they act as heating elements (heated by Joule effect) and are also used as temperature probes. A fluid loop allows controlling the liquid injection flow rate, its inlet temperature as well as its pressure. The test section is held vertically, the liquid injected from bottom at a temperature close to the saturation temperature. In a first series of experiments, the thermal power applied to a bundle of heating cylinders is progressively increased until a dry zone is detected in the porous medium. Two kinds of phenomenology are observed during these “dryout experiments”. First, at low liquid injection rate (4 kg.m^-2.s^-1 maximum mass flux), reaching the dryout power results into a liquid front receding down to the upper limit of the heated zone, while downstream the heated zone, the porous medium is vapour-saturated. Second, at higher flow rate, the boiling crisis happens at the surface of a single heating element, resulting in a local film boiling, whereas a two-phase flow still go through the whole test section. High-speed visualizations allow characterizing the flow regimes. Other experiments focus on determining the local critical heat flux on a given cylinder, for different upstream flow configurations. The inlet liquid flow rate is fixed. A thermal power is uniformly applied to a line of heating cylinders, upstream the cylinder under study. Results show that the local critical heat flux decreases as the power applied to the heated line increases. The distance from the cylinder under study to the heated line seems not to have a significant effect on the critical heat flux. Visualizations are used to characterize the two-phase flow upstream the heated line, aiming at expressing the critical heat flux as a function of the hydrodynamic parameters (saturations, phase velocities). The image analysis is particularly challenging. In order to calibrate the image processing parameters, we use a second model porous medium with the same geometry as the heat generating one, but where an isothermal two-phase flow is obtained by injecting gas into the liquid flow rather than generated by boiling. The gas injection flow rate is controlled and measured. Isothermal two-phase flow visualizations provide a reference case and are compared to flow boiling visualizations.