L'objectif de ce travail de thèse concerne l'étude et la mise en œuvre de deux modèles de cavitation dans le solveur NSMB (Navier-Stokes-Multi-Blocks): les modèles HEM (Homogeneous Equilibrium Model) et une équation pour le taux de vide: le modèle à transport de taux de vide (TTV). Le phénomène de cavitation est modélisé par différentes équations d'état de mélange liquide-vapeur (EOS). Des simulations numériques sont réalisées sur des écoulements diphasiques compressibles unidimensionnels et bidimensionnels avec des conditions d'interface et comparées à des solutions de référence. De plus, la méthode TTV basée sur le taux de vide incluant les termes source pour la vaporisation et la condensation dans le logiciel libre open source OpenFOAM est également présentée sur la géométrie Venturi pour capturer le phénomène du jet réentrant. La modélisation de la turbulence joue un rôle majeur dans la capture des comportements instationnaires et un limiteur est introduit pour réduire la viscosité turbulente afin de mieux prédire la structure à deux phases. Une comparaison de divers modèles de cavitation couplés avec des modèles de turbulence est étudiée. Les résultats computationnels sont comparés aux données expérimentales existantes. / The objective of this thesis work concerns the study and implement of two cavitation models in the NSMB (Navier-Stokes-Multi-Blocks) flow solver: the Homogeneous Equilibrium Models (HEM) and a void ratio Transport-based Equation Model (TEM). The cavitation phenomenon is modeled by different liquid-vapor mixture equation of state (EOS). Numerical simulation are performed on some one- and two-dimensional compressible two-phase flows with interface conditions and compared with reference solutions. Moreover, The TEM based method for the void ratio including the source terms for vaporization and condensation in the free, open source software OpenFOAM is also presented on the Venturi geometry to capture the re-entrant jet phenomenon. The turbulence modeling plays a major role in the capture of unsteady behaviors and a limiter is introduced to reduce the eddy-viscosity to better predict the two-phase structure. A comparison of various cavitation models coupled with turbulence models are investigated. Computational results are compared with existing experimental data.

La récupération d'hydrocarbures non conventionnels fait partie des enjeux énergétiques majeurs. Ils ne peuvent être extraits par simple forage car la roche qui les contient, constituée essentiellement de nanopores, présente une très faible perméabilité. A l'échelle macroscopique, c'est-à-dire à l'échelle du bassin, les écoulements de fluides sont décrits par la loi de Darcy qui relie le flux à la perméabilité, au gradient de pression et à la viscosité. La perméabilité d'un matériau peut être mesurée expérimentalement ou théoriquement par homogénéisation à partir de l’hydrodynamique continue. Cependant, lorsque la taille des pores devient comparable à celle des molécules de fluide, une telle description n'est pas satisfaisante. D’une part l’hydrodynamique continue, où la nature du fluide n’intervient qu’à travers la viscosité, ne suffit pas forcément pour décrire l’écoulement. D’autre part les interactions au niveau moléculaire entre le fluide et le solide jouent un rôle important. Cette thèse porte sur le transport de fluides à l'échelle moléculaire et revisite la description traditionnelle qui sert de point de départ pour des écoulements à l'échelle macroscopique, en particulier dans le cas des écoulements multiphasiques. Par des simulations de dynamique moléculaire classique, nous avons étudié l'écoulement de systèmes monophasiques et diphasiques, précisant l’influence de la nature des surfaces, ainsi que de la nature et de la concentration des espèces dissoutes. Nous avons également apporté une contribution méthodologique originale pour le calcul des coefficients de diffusion d'espèces. / Unconventional hydrocarbons recovery is one of the major energy challenges. They cannot be extracted by simple drilling because the rock which contains them, consisting essentially of nanopores, has a very low permeability. On the macroscopic scale of the geological basin, the flows of fluids are described by Darcy's law which connects the flux to the permeability, the pressure gradient and the viscosity. The permeability of a material can be measured experimentally or determined theoretically by homogenization from continuous hydrodynamics. However, when the pore size becomes comparable to that of the fluid molecules, such a description is unsatisfactory. On the one hand continuous hydrodynamics, where the nature of the fluid only enters via the viscosity, is not necessarily sufficient to describe the flow. On the other hand, the interactions at the molecular level between the fluid and the solid play an important role. This thesis deals with the transport of fluids on a molecular scale and revisits the traditional description which serves as a starting point for macroscopic flows, in particular in the case of multiphase flows. Using classical molecular dynamics simulations, we study the flow of one- and two-phase systems, specifying the influence of the nature of the surfaces, as well as the nature and concentration of the dissolved species. We also develop an original methodological contribution to the calculation of the diffusion coefficients of confined species, specifying the effects of the system finite size.

Un caloduc pulsé diphasique est un lien thermique composé d'un tube capillaire lisse sous forme de serpentin reliant un évaporateur à un condenseur, séparés par une partie adiabatique. Les conditions de température et de pression du fluide à l'intérieur du caloduc sont proches des conditions de changement de phase. De ce fait, et grâce aux dimensions capillaires du tube, le fluide se distribue en différentes parties liquide et vapeur distribuées de manière alternée. Les instabilités thermo-hydrauliques permanentes sont à l'origine d'un écoulement oscillant qui permet le transfert de chaleur de l'évaporateur jusqu'au condenseur.L'objectif du présent projet de recherche consiste à étudier le comportement thermo-hydraulique de trois caloducs cryogéniques pulsés diphasiques testés avec différents fluides cryogéniques (azote, néon et argon) pour le refroidissement d'aimants à haute température critique. De plus, un code numérique a été développé pour les futures simulations 2D des caloducs pulsés diphasiques.Au cours de ce projet de recherche, de nombreux tests expérimentaux ont été réalisés avec trois fluides cryogéniques différents: azote, néon et argon. Les résultats expérimentaux des tests avec une augmentation de puissance progressive dans l'évaporateur ont révélé des capacités de transfert thermiques très différentes en fonction du fluide, chaque fluide présentant un comportement thermo-hydraulique différent. L'état thermodynamique du fluide lors du fonctionnement stable du PHP et la phase d'assèchement (dry-out) ont été étudiés. Les différences dans le comportement des différents fluides ont été expliquées après l'analyse de leurs propriétés physiques. De plus, les taux de remplissage de fluide dans le PHP donnant les meilleures performances thermiques ont été définis. Ajouté à cela, de nombreux tests réalisés en configuration ouverte (avec le PHP connecté au volume tampon) et en configuration fermé (avec le PHP isolé du volume tampon) ont permis de conclure sur la capacité de régulation du volume tampon en cas de surpression dans le PHP. Aussi, les résultats expérimentaux des longs tests de stabilité ont permis de vérifier la stabilité du système PHP pendant des longues périodes de fonctionnement. Par ailleurs, des tests spécifiques ont été réalisés pour déterminer des conditions optimales de démarrage, l'influence de la température du condenseur dans les performances thermiques du système et l'influence du nombre de tubes en parallèle dans la capacité de transfert thermique du système. Finalement, une série de tests avec une forte puissance thermique imposée au niveau de l'évaporateur imitant une situation de quench dans un aimant supraconducteur ont données des précieuses informations sur les limites thermiques du système. Concernant les simulations numériques, un modèle a été développé avec le solveur Fluent pour des simulations dans une géométrie 2D axisymétrique en utilisant la méthode VOF. La dynamique du fluide dans un tube capillaire a été modélisée et les simulations thermiques ont permis de conclure que les instabilités thermodynamiques restent insuffisantes pour maintenir les oscillations du fluide. Ce modèle est présenté comme une nouvelle plateforme pour de futures modélisations 2D des caloducs pulsés diphasiques. / A pulsating (or oscillating) heat pipe (PHP or OHP) is a heat transfer device composed of a single capillary tube bent in many U-turns, connecting an evaporator to a condenser, separated by an adiabatic part. In the PHP, temperature and pressure conditions of the working fluid are close to phase-change conditions. Due to this and to the capillary dimensions of the tube, the fluid is distributed in alternating liquid slugs and vapor plugs. Permanent thermal instabilities in the PHP create the oscillating flow which allows the transfer of heat from one end (the evaporator) to the other (the condenser).The objective of the present work consists in characterizing the thermo-hydraulic behavior of the meter-scale horizontal cryogenic pulsating heat pipes as a cooling solution for space superconducting magnets. To this, several experiments have been conducted in a cryogenic facility containing three different horizontal pulsating heat pipes. In addition, a numerical 2D model has been proposed for future horizontal pulsating heat pipes simulations.During the research project, numerous tests have been performed using three different working fluids: nitrogen, neon and argon. From experimental results of progressive heat load tests it has been possible to compare the maximum heat load transfer capacity of the PHP with each fluid and the corresponding thermal performance. It has also been noticed that each fluid presents a specific behavior concerning the fluid oscillations. In addition, the thermodynamic state of the fluid in operating conditions and the dry-out process have been characterized. Differences between fluid's behaviors have been partly explained by analyzing the evolution of the fluid physical properties related to the movement and the heat transfer capacity. Furthermore, it has been possible to conclude about the relation between the liquid filling ratio in the PHP and its thermal performance, determining the filling ratios giving the highest thermal performances. Moreover, similar tests have been performed in open configuration (with the PHP connected to the buffer volume) and closed configuration (with the PHP isolated from the buffer volume). From this, it has been possible to conclude about the regulation made by the buffer volume in case of overpressure in the PHP. Also, experimental results from long stability tests have confirmed that these pulsating heat pipe are able to work in stable conditions during long periods as a reliable cooling system. In addition to that, specific tests have been done to determine the optimum start-tup conditions, the influence of the temperature of the condenser in the thermal performance and the influence of the number of turns in the global heat transfer capacity. A final series of tests have been achieved with a sudden extra heat load at the surface of the evaporator while the PHP is operating in stable conditions, simulating a quench event of a superconducting magnet. Experimental results gave us precious information about the transient thermal behavior and operating limits of this kind of device during transient heat loads like quench situations. Concerning the numerical part, a numerical model has been proposed for transient simulations with a pressure-based Fluent solver using the Volume of Fluid (VOF) method in a 2D axisymmetric geometry. Certain characteristics of fluid dynamics in capillary tubes have been confirmed. It has also been noticed that thermodynamic instabilities are not enough to generate the fluid oscillations in capillary tubes. Even if the 2D axisymmetric simulation is still at its early stages, several aspects of the models have been validated after analyzing the evolution of different parameters, suggesting that this kind of model can be considered as a new platform for future 2D pulsating heat pipes simulations.

L'objectif de cette étude est d'examiner les écoulements diphasiques liquide/gaz dans les canaux des plaques bipolaires des piles à combustible afin de comprendre et de trouver des solutions au problème d'engorgement. L'influence de la section du canal et du matériau utilisé sur les pertes de charge (PDC) et sur la structure de l'écoulement d'eau liquide est plus particulièrement étudiée dans une expérience hors pile. Les mesures des PDC ont permis de mettre en évidence des séquences de bouchages et de débouchages, la fréquence de ces dernières augmentant avec le débit. Par ailleurs il est montré que le rapport des PDC diphasiques moyennées par les PDC en air sec décroit avec le débit et ceci indépendamment de la dimension du canal et qu'il est d'autant plus faible que le revêtement est hydrophile. Enfin différents régimes d'écoulements diphasiques (stratifié et de gouttes) ont pu être mis en évidence dans les différentes zones du canal et un modèle pour chacun d'eux a été établi, confirmant les résultats expérimentaux / This work focuses on the gas-liquid flows in the cathode plate, with the objective to observe their patterns, to understand their behavior, to estimate the pressure drops (PD) and eventually, to reduce clogging and its possible consequences in term of oxygen starvation downstream. A special emphasis is put on the effect of the channel section (typically between 0.5 and 1 mm²) and on the surface properties of the flow field plate materials. The experiments are performed ex-situ. The PD is measured locally along the channel as well as globally between the inlet and outlet, which put forward the existence of clogging/unclogging sequences. The characteristic frequency of these sequences increases with the air flow rate. The results show that the ratio of PD in two-phase flow to PD in dry flow decreases with the air flow rate while it does not seem to depend on the channel size (within the tested range). Moreover this ratio is lower with hydrophilic coating. Finally two main flow patterns (slug and annular flow) have been observed depending on the distance from the inlet and they have been simulated

Les turbines aéronautiques doivent satisfaire à des normes d'émissions polluantes toujours en baisse. La qualité du mélange du carburant et de l'air dans la chambre de combustion est responsable de la formation de polluants nocifs pour l'environnement. La simulation aux grandes échelles (LES) permet d'étudier les mécanismes de mélanges turbulents de l'air et du carburant. La prise en compte de l'aspect liquide du carburant injecté devient nécessaire pour prédire correctement l'apparition de vapeur de carburant au sein du foyer. Le but de cette thèse est évaluer la fiabilité des simulations LES Euler-Euler dans une configuration complexe. Les processus d'injection, et d'évaporation du carburant liquide sont analysés et modélisés dans les simulations LES car ils pilotent la formation de vapeur de carburant. Les méthodes numériques pour résoudre les équations continues de la phase dispersée doivent permettre des simulations précises et robustes dans une configuration représentative d'une chambre de combustion. Les simulations présentées dans ces travaux reproduisent l'écoulement diphasique évaporant non-réactif du banc d'essai Mercato. Ce banc est équipé d'un système d'injection d'air vrillé et d'un atomiseur pressurisé-swirlé de kérosène typiques des foyers aéronautiques réels. Dans ces travaux, le modèle pour l'injection de liquide FIM-UR a été développé pour définir les conditions limites conduisant à un spray issu d'un atomiseur préssurisé-swirlé. Le kérosène employé dans les campagnes expérimentales est modélisé dans les simulations par un composé permettant d'obtenir des temps d'évaporation réalistes. Trois stratégies numériques ont été mises en place sur la configuration MERCATO. Les comparaisons des résultats numériques aux mesures expérimentales ont permis d'évaluer la stratégie numérique conduisant à la meilleure précision. L'utilisation du schéma centré TTGC associé à un opérateur de viscosité artificielle localisée par un senseur adapté est optimale lorsque l'équation sur l'énergie décorrélée des gouttes est résolue. Cette stratégie permet de contrôler la localisation et les niveaux de viscosité par rapport à un schéma décentré. Les termes sources liés au mouvement mésoscopique permettent de redistribuer l'énergie dans les zones de compression ou de détente de la phase dispersée, et d'obtenir les bonnes répartitions des fluctuations dans la chambre de combustion. La stratégie retenue est comparée aux statistiques de la dynamique du spray résolu par une approche Lagrangienne employant la même injection monodispersse. Le méthode Euler-Euler conduit à la même précision de la dynamique de la phase dispersée que la méthode Euler-Lagrange. L'accès à l'évolution instationnaire de l'écoulement permet d'identifier les mêmes mécanismes de dispersion et de mélange dans les deux simulations. Des différences sur la répartition de diamètre moyen et de carburant dans la chambre ont été mis en évidence et reliés à la polydispersion locale qui n'est pas résolue dans l'approche Euler-Euler monodisperse et qui apparaît naturellement dans l'approche Euler-Lagrange malgré l'injection monodisperse.

Ce travail est consacré à l'analyse et au développement de concepts et d'outils mathématiques en vue de leur application à des problématiques propres aux sites de stockage géologique profond de déchets radioactifs. La première partie porte sur l'estimation en champ lointain de la concentration de radionucléides issus du relâchement des colis de confinement, lorsque les incertitudes sur le relâchement sont prises en compte. En s'appuyant sur les travaux de A. Bourgeat et A. Piatniski sur l'homogénéisation d'une équation de convection-diffusion avec second membre aléatoire, on développe des outils numériques permettant d'approcher le comportement probabiliste du champ de concentration dans une configuration du type site de stockage. Dans une seconde partie, on s'intéresse à la migration de gaz dans et autour d'un site de stockage. Après une revue sur la modélisation physique des écoulements diphasiques de type eau/hydrogène en milieu poreux, on propose une nouvelle formulation mathématique du problème qui décrit, dans un même jeu d'équations, les écoulements à une (liquide) et deux (liquide/gaz) phases. Une étude de l'existence de solutions de cette formulation est menée à l'aide de la théorie générale des équations différentielles quasilinéaires elliptiques-paraboliques introduite par H.W. Alt et S. Luckhaus. Une méthode de résolution numérique du problème est mise en oeuvre pour la simulation de différents cas test, des plus simples au plus représentatif d'un site de stockage géologique. Enfin, l'homogénéisation périodique du modèle est effectuée et appliquée à la simulation de l'exercice Couplex-Gaz proposé par l'ANDRA

Nous étudions expérimentalement le transport des bouchons liquides dans les réseaux branchés de microcanaux. Les bouchons liquides sont poussés sous différentes conditions dans un réseau arborescent de microcanaux de largeur croissante ou décroissante avec le numéro de génération. Le motif global d'écoulement peut être symétrique ou asymétrique, avec des bouchons se divisant de manière synchrone ou asynchrone, selon la nature et l'intensité de la force de poussée et la géométrie du réseau. Ce comportement s'explique par des interactions entre bouchons à l'oeuvre au sein d'un élément fondamental du réseau composé de trois bifurcations adjacentes. Quand un unique bouchon est poussé à pression constante, ses descendants atteignent les sorties du réseau uniquement si la pression est supérieure aux seuils du réseau. Les descendants d'un bouchon poussé à débit constant, atteignent toujours les sorties d'un réseau de largeur décroissante, même lorsque le débit est faible. Inversement, dans un réseau à largeur croissante, seuls certains descendants peuvent atteindre les sorties alors que les autres restent bloqués aux bifurcations intermédiaires du réseau. De plus, nous montrons qu'une relation linéaire existe entre la pression et le débit total dans le réseau, dont la validité est vérifiée en présence de bouchons successifs. Dans ce cas, le débit total varie avec la distance initiale entre les bouchons, à pression et longueurs de bouchons fixées. Finalement, des résultats préliminaires sont présentés, concernant la réouverture sous haute pression de réseaux par ruptures des bouchons.

Sous la direction de SNECMA, un nouveau moteur cryotechnique pour Ariane 5 est en phase de développement. Ce moteur, VINCI, fonctionne à l'oxygène et à l'hydrogène liquides et devra être rallumable en cours de vol. Le transitoire de remplissage d'une cavité intermédiaire par laquelle transite l'oxygène liquide avant d'arriver dans la chambre de combustion s'avère être une étape critique qu'il faut s'efforcer de maîtriser. Cette cavité, appelée dôme LOX, est directement reliée à la chambre de combustion par les injecteurs à oxygène. Des outils numériques sont actuellement en cours de conception : ils permettront à terme de simuler le remplissage de cette cavité dans l'espace. Afin de valider ces outils numériques, un programme expérimental a vu le jour, impliquant SNECMA, le CNES (Centre National d'Etudes Spatiales) et le LEGI (Laboratoire des Ecoulements Géophysiques et Industriels). Il s'agit de mener à bien des expériences " simples " et reproductibles afin de disposer d'une base de données expérimentales qui servira de cas tests pour les simulations. Un banc d'essais expérimental a été progressivement mis en place et instrumenté au LEGI. Le fluide utilisé en substitution de l'oxygène liquide est de l'eau. On a choisi de respecter une similitude du nombre de Weber entre le cas " réel " et les conditions expérimentales. Ce dernier équivaut au rapport des forces d'inertie sur les forces de tension de surface. Deux campagnes d'essais ont été réalisées, dans lesquelles on s'est centré sur l'étude des aspects hydrodynamique du transitoire de remplissage : variation du débit total et des pressions au cours d'un essai, évaluation du taux de vide dans la cavité, visualisation de l'écoulement dans la cavité et en sortie des injecteurs. L'instrumentation à disposition est constituée d'un débitmètre à effet Coriolis, de capteurs de pression, d'une sonde optique, de caméras rapides et d'un laser pour l'imagerie. La première campagne d'essais a visé le remplissage de la cavité en eau seule. Le paramètre de contrôle principal était la pression génératrice de l'écoulement liquide. Dans la deuxième campagne d'essais on injecte simultanément dans la cavité un écoulement d'eau et un écoulement d'air. Cela se rapproche plus des conditions réelles du transitoire de remplissage, au cours duquel la cavité est balayée par un écoulement d'hélium. Le paramètre de contrôle supplémentaire est le débit de gaz initial. Ces campagnes ont également souligné l'importance du profil d'ouverture de la vanne de l'écoulement liquide sur le transitoire de remplissage de la cavité. Ces campagnes d'essais constituent une première étape dans la compréhension du transitoire de remplissage du dôme LOX. Elles ont permis de visualiser la forme de l'écoulement dans la cavité et en sortie des injecteurs et d'identifier certains phénomènes intervenant dans le remplissage de la cavité. En particulier, nous avons mis en évidence l'existence d'un délai de mise en place de l'écoulement par les injecteurs, qui peut être responsable d'un pic de pression dans la cavité au cours du transitoire. L'influence de la fraction gazeuse sur l'écoulement dans les injecteurs a été soulignée mais reste à quantifier de façon précise. La prochaine étape de l'étude concerne les aspects énergétiques du transitoire de remplissage, notamment les transferts thermiques ayant lieu, dans la réalité, entre l'hélium, l'oxygène et les parois du dôme LOX.

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.

Dans le cadre de la recherche sur les accidents graves pour la deuxième et la troisième génération de réacteurs nucléaires, certains aspects de l'ablation de béton dans le puits de cuve au cours de l'interaction corium-béton (ICB) restent encore inexpliquées. La détermination d'échange de chaleur le long de la région interfaciale entre un bain de corium et un béton est importante pour l'évaluation de la progression d'ablation du béton et, éventuellement, la percée de fondation. Le projet CLARA s'inscrit une recherche expérimentale sur la thermohydraulique au sein d'un bain de liquide agitée par des bulles de gaz. Les essais CLARA sont réalisés avec des matériaux simulants. Ils permettent de mettre en évidence l'influence de la vitesse superficielle du gaz, de la viscosité du liquide et de la géométrie sur le coefficient d'échange de chaleur entre le bain de liquide chauffé et les parois verticales et horizontales de la piscine qui sont maintenues à une température uniforme. La première campagne d'essais a été réalisée avec la configuration du bain de petite taille (50 cm × 25 cm × 25 cm). Les essais ont été réalisés avec des liquides couvrant un large éventail de viscosité dynamique, d'environ 1 mPa s à 10000 mPa s. La vitesse superficielle du gaz est modifiée jusqu'à 8 cm/s. Cette thèse comporte une brève description de la phénoménologie de l'ICB, une synthèse bibliographique sur les corrélations d'échange de chaleur existantes pour l'écoulement diphasique et le taux de vide, une description de l'installation CLARA, les résultats des essais et leur interprétation. Les résultats expérimentaux sont comparés avec les modèles existants et certains nouveaux modèles pour l'évaluation du coefficient d'échange de chaleur dans un écoulement diphasique. / In the context of severe accident research for the second and the third generation of nuclear power plants, there are still open issues concerning some aspects of the concrete cavity ablation during the molten corium - concrete interaction (MCCI). The determination of heat transfer along the interfacial region between the molten corium pool and the ablating basemat concrete is crucial for the assessment of concrete ablation progression and eventually the basemat meltthrough. For the purpose of experimental investigation of thermalhydraulics inside a liquid pool agitated by gas bubbles, the CLARA project has been launched. The CLARA experiments are performed using simulant materials and they reveal the influence of superficial gas velocity, liquid viscosity and pool geometry on the heat transfer coefficient between the internally heated liquid pool and vertical and horizontal pool walls maintained at uniform temperature. The first test campaign has been conducted with the small pool configuration (50 cm × 25 cm × 25 cm). The tests have been performed with liquids covering a wide range of dynamic viscosity from approximately 1 mPa s to 10000 mPa s and the superficial gas velocity is varied up to 8 cm/s. This thesis comprises a brief description of MCCI phenomenology, literature reviews on the existing heat transfer correlations for twophase flow and the void fraction, a description of CLARA setup, experimental results and their interpretation. The experimental results are compared with existing models and some new models for the assessment of heat transfer coefficient in two-phase flow. / U kontekstu istraživanja teških nesre´ca u nuklearnim elektranamadruge i tre´ce generacije, neka pitanja vezana za ablaciju temelja kontejnmentatijekom interakcije rastaljenog korijuma i betona i dalje ostajuotvorena. Odred¯ivanje prijenosa topline u površinskom podrucˇjuizmed¯u bazena rastaljenog korijuma i betona kljucˇno je za odred¯ivanjenapredovanja ablacije i u konaˇcnici procjene vremena rastapanjatemelja kontejnmenta. U svrhu eksperimentalnog istraživanja prijenosatopline u tek´cinama miješanima ubrizgavanjem zraka, pokrenutje projekt nazvan CLARA.CLARA eksperimenti izvode se koriste´ci imitacijske materijale i otkrivajuutjecaj fiktivne brzine plina, viskoznosti teku´cine i geometrijebazena na koeficijent prijenosa topline izmed¯u grijanog bazena te njegovihvetrikalnih i horizontalnih stijenki ˇcija se temperatura održavana konstantnoj temperaturi. Prva serija eksperimenata provedena je sbazenom male konfiguracije (50 cm × 25 cm × 25 cm). Eksperimentisu izvedeni s teku´cinama dinamiˇcke viskoznosti od približno 1 mPas do 10000 mPa s, dok je maksimalna fiktivna brzina plina 8 cm/s.Ova disertacija sadrži kratak opis fenomenologije procesa interakcijerastaljenog korijuma i betona, pregled postoje´cih korelacija zaviprijenos topline u dvofaznom toku i korelacija za poroznost, opisCLARA eksperimentalne postave, rezultate eksperimenta i njihovuinterpretaciju. Rezultati eksperimenta su uspored¯eni s predvid¯anjimaprema postojec´im modelima. Predloženi su takod¯er i neke nove korelacijeza odred¯ivanje koeficijenta prijenosa topline u dvofaznom toku.