L’étude de l’ébullition transitoire est un enjeu important pour la sureté nucléaire. Un tel phénomène peut se produire lors d’un accident de type RIA (Reactivity Initiated Accident)dans un réacteur nucléaire où le pic de puissance au niveau d’un crayon de combustible peut déclencher une ébullition transitoire conduisant à une forte augmentation de la température de la gaine et à un risque de rupture. Plusieurs études en conditions réacteurs ont permis d’obtenir des courbes d’ébullition transitoires mais la modélisation qui en découle manque encore de fiabilité. Dans le cadre d’une collaboration avec l’Institut de Radioprotection et de Sûreté Nucléaire (IRSN), une expérience modèle a été construite à l’Institut de Mécanique des Fluides de Toulouse (IMFT). Elle génère un écoulement de réfrigérant HFE7000 dans un canal de section semi-annulaire, simulant l’écoulement autour d’un crayon de combustible, dont la partie intérieure, composée d’une feuille de métal, est chauffée rapidement par effet Joule, simulant l’échauffement de la gaine du crayon. La thermographie infra-rouge permet de mesurer la température de la paroi externe du métal. L’application d’une peinture noire sur le métal augmente son émissivité mais aussi la résistance thermique de la paroi. La précision de la mesure de la température d’intérêt a été optimisée en fonction de l’épaisseur de peinture et une correction sur le bilan d’énergie prend en compte ce paramètre. Ces mesures sont couplées avec une caméra rapide qui permet de visualiser les régimes d’ébullition et d’obtenir des tailles de bulles à l’aide de la mise en place d’algorithmes de traitement d’image. On représente sur un diagramme flux-température les transferts thermiques lors des différents régimes en stationnaire et en transitoire. Chaque régime d’ébullition, en conditions stationnaire ou transitoire, est alors passé en revue : la convection, le déclenchement de l’ébullition, l’ébullition nucléée, la crise d’ébullition, l’ébullition en film et le remouillage. Les régimes stationnaires sont correctement modélisés par des corrélations usuelles. La convection transitoire est caractérisée sur toute la paroi et son évolution se rapproche de la solution quasistationnaire. Il est montré que les transferts thermiques lors du passage vers l’ébullition nucléée sont dépendants de la formation d’une importante poche de vapeur qui se propage sur la paroi. Une étude locale de cette propagation est alors nécessaire. Afin de simuler des transitoires de température durant l’ébullition nucléée, un système d’asservissement de type P.I.D. permet d’imposer des créneaux ou des rampes de températures (de 5 à 500 K.s 1 ). Les résultats en ébullition nucléée sont conformes avec ceux de la littérature, tant en conditions stationnaire que transitoire. L’expérience permet d’étudier le transfert de chaleur lorsqu’un film de vapeur se forme et isole la paroi. Ce régime d’ébullition en film, pendant la chauffe ou le refroidissement de la paroi peut ainsi être stabilisée pendant plusieurs secondes avec ce système. On caractérise ainsi les conditions de déclenchement de l’ébullition en film, la dynamique de sa propagation et les transferts une fois établi. Enfin, l’implémentation des caractéristiques physiques de notre expérience dans le code SCANAIR de l’IRSN, permet de commencer à calculer et comparer nos résultats expérimentaux avec les simulations numériques. Des calculs de conduction instationnaire sont notamment considérés en imposant la température mesurée pour analyser nos résultats lors du régime de convection et après le déclenchement de l’ébullition. / The study of rapid transient boiling is an important issue in the nuclear safety. Such a phenomenon may occur in the case of a RIA (Reactivity Initiated Accident) in the core of a nuclear reactor powerplant, where a power excursion can trigger the formation of a vapour film around the fuel rod, leading to an important rise of the rod temperature and a risk of failure. Some studies in reactor conditions provided transient boiling curves but the modeling lacks of reliability. In collaboration with the IRSN (Institut de Radioprotection et de Sûreté Nucléaire), an experiment model was built at the Institute of Fluid Mechanics of Toulouse. It generates the flow of a refrigerant, HFE7000, in a semi-annular section channel, whose inner wall is made of a metal foil rapidly heated by Joule effect, simulating the heating of a fuel rod. Infrared thermography is used to measure the temperature of the metal foil, painted with a black paint to increase its emissivity, causing also an increase of the wall thermal resistance. The measurement accuracy of the interest temperature has been optimized according to the paint thickness and a correction on the energy balance takes account this parameter. These measurements are coupled with a high-speed camera that allows visualizing the boiling regimes and get bubble sizes using image processing algorithms. On a flux-temperature diagram, the heat transfers are represented both for steady and transient regimes. Each boiling regime is then reviewed : convection, onset of nucleate boiling, nucleate boiling, boiling crisis, film boiling and rewetting. Steady regimes are correctly modeled by usual correlations. Transient convection is characterized over the whole wall and its evolution is closed to the quasi-steady solution. It is shown that heat transfer during the transition to nucleate boiling are strongly related to the formation of a large vapor phase that spreads on the wall. A local study of this propagation is then necessary. In order to simulate and control transient temperature during nucleate boiling, a P.I.D. is implemented to impose a steady or ramps temperature (from 5 to 500 K.s 1 ). The results in nucleate boiling make it possible to recover the results of the literature in both steady and transient conditions. The experiment allows to study the heat transfer when a vapor film is formed and insulates the wall. The film boiling regime during heating or the cooling of the wall can thus be stabilized for several seconds with this system. The conditions for triggering of film boiling are thus characterized, as its spread dynamic and its transfers once established. Finally, the implementation of the physical characteristics of our experience in IRSN’s SCANAIR code allows us to begin to calculate and compare our experimental results with numerical simulations. Unsteady conduction calculations are applied to the measured temperature to analyze our results during the convection regime and after the onset of boiling.

Práce komplexně popisuje vodní a vodovzdušné chlazení pomocí metod CFD (Computational Fluid Dynamics) konkrétně s využitím software ANSYS FLUENT. Skládá se ze dvou hlavních částí, z nichž první se zabývá numerickým popisem jediné vodní kapky a druhá popisem směsí kapek představující paprsek válcové a ploché trysky. Je založena převážně na vícefázových modelech proudění a vlastních uživatelsky definovaných funkcí (User Defined Functions, UDF) představujících stěžejní část práce. Uvedené výpočtové modely jsou ve většině případů verifikovány pomocí experimentálních dat nebo jiných numerických modelů. V první části práce jsou teoreticky postupně rozebrány všechny tři použité vícefázové modely proudění. První z nich, Volume Of Fluid model (VOF), byl použit pro modelování jediné kapky (mikromodel). Zatímco zbývající dva, Euler-Euler model a Euler-Lagrange model, byly aplikovány v modelu celého paprsku trysky (makromodel). Mikromodel popisuje dynamiku volného pádu vodní kapky. Pro malé průměry kapek (~100µm) standardní model povrchového napětí (Continuum Surface Force, CSF) způsoboval tzv. parazitní proudy. Z toho důvodu je v práci rozebrána problematika výpočtu normál, křivostí volných povrchů a povrchového napětí jako zdroje objemových sil v pohybových rovnicích. Makromodel se zabývá studiem dynamiky celého paprsku tj. oblastí od ústí trysky po dopad na horký povrch, bere v úvahu kompletní geometrii, tzn. např. podpůrné válečky, bramu, spodní část krystalizátoru apod. V práci je rozebrána 2D simulace dopadu paprsku válcové trysky pomocí VOF modelu Euler-Lagrange modelu na horký povrch. Pro případ s VOF modelem byl navržen model blánového varu. Euler-Euler model a Euler-Lagrange model byly využity pro simulaci paprsku ploché trysky horizontálně ostřikující horkou bramu přímo pod krystalizátorem nad první řadou válečků. Pro Euler-Euler model byl navržen model sekundárního rozpadu paprsku založený na teorii nejstabilnější vlnové délky (Blob jet model). Jelikož diskrétní Lagrangeovy částice tvořily v určitých místech spíše kontinuální fázi, byl navržen a otestován model pro konverzi těchto částic do VOF.

<p>In this study, universal databases and semi-empirical correlations are developed for cryogenic two-phase heat transfer and pressure drop in heated tubes undergoing flow boiling.</p>