Les accidents graves de centrales nucléaires ont une probabilité d’occurrence très faible, mais compte tenu des risques encourus, il est nécessaire de savoir prédire l’évolution de l’accident. Dans le scénario le plus critique, le dégagement de chaleur induit par la désintégration des produits de fission entraînerait la fusion du cœur et la formation d’un magma (« corium ») qui tomberait sur le radier en béton du bâtiment réacteur, provoquant sa fusion. L’objectif des études est d’évaluer la vitesse de fusion du béton. Dans ce contexte, le travail effectué dans cette thèse se situe dans la continuité du modèle de ségrégation de phases développé par Seiler et Froment, et s’appuie sur les résultats expérimentaux des essais ARTEMIS. D’une part, nous avons développé un nouveau modèle de transferts à travers le milieu interfacial. Ce modèle fait intervenir trois mécanismes de transfert : la conduction, la convection et un dégagement de chaleur latente. D’autre part, nous avons revu la modélisation couplée du bain et du milieu interfacial, ce qui a conduit au développement de deux nouveaux modèles : « le modèle liquidus », pour lequel on suppose qu’il n’y a pas de résistance au transfert de soluté, et le « modèle à épaisseur de milieu interfacial constante », pour lequel on suppose qu’il n’y a pas de dissolution du milieu interfacial. Le modèle à épaisseur de milieu interfacial constante permet de prédire correctement les valeurs expérimentales de la vitesse de fusion du béton et de la température du bain, dans les essais 3 et 4 tandis que le modèle liquidus, appliqué aux essais 2 et 6, prédit correctement l’évolution de la vitesse de fusion et de la température du bain / Severe accidents of nuclear power plants are very unlikely to occur, yet it is necessary to be able to predict the evolution of the accident. In some situations, heat generation due to the disintegration of fission products could lead to the melting of the core. If the molten core falls on the floor of the building, it would provoke the melting of the concrete floor. The objective of the studies is to calculate the melting rate of the concrete floor. The work presented in this report is in the continuity of the segregation phase model of Seiler and Froment. It is based on the results of the ARTEMIS experiments. Firstly, we have developed a new model to simulate the transfers within the interfacial area. The new model explains how heat is transmitted to concrete: by conduction, convection and latent heat generation. Secondly, we have modified the coupled modelling of the pool and the interfacial area. We have developed two new models: the first one is the “liquidus model”, whose main hypothesis is that there is no resistance to solute transfer between the pool and the interfacial area. The second one is “the thermal resistance model”, whose main hypothesis is that there is no solute transfer and no dissolution of the interfacial area. The second model is able to predict the evolution of the pool temperature and the melting rate in the tests 3 and 4, with the condition that the obstruction time of the interfacial area is about 105 s. The model is not able to explain precisely the origin of this value. The liquidus model is able to predict correctly the evolution of the pool temperature and the melting rate in the tests 2 and 6

Lors d'un hypothétique accident grave de réacteur à eau sous pression, un mélange de matériaux fondus, appelé corium, issu de la fusion du cœur peut se relocaliser dans le puits de cuve constitué par un radier en béton. Les codes d'évaluation réacteur pour simuler la phénoménologie de l'interaction corium-béton sont basés sur une description à grande échelle des échanges qui soulève de nombreuses questions, tant sur la prise en compte des phénomènes multi-échelles mis en jeu que sur la structure adoptée de la couche limite au voisinage du front d'ablation. Dans ce contexte, l'objectif principal de ce travail consiste à aborder le problème de la structure de la couche limite par simulation numérique directe. Ce travail s'inscrit dans le cadre plus général d'une description et d'une modélisation multi-échelle des échanges, c'est-à-dire de l'échelle locale associée au voisinage du front d'ablation jusqu'à l'échelle du code d'évaluation réacteur. Une telle description multi-échelle des échanges soulève le problème de la description locale de l'écoulement multiphasique multiconstituant mais aussi le problème du changement d'échelle et en particulier le passage de l'échelle locale à l'échelle de description supérieure dite macroscopique associée aux mouvements convectifs dans le bain de corium. Parmi les difficultés associées au changement d'échelle, nous nous intéressons à la problématique de la construction de conditions aux limites effectives ou lois de parois pour les modèles macroscopiques. Devant la complexité du problème multiphasique multiconstituant posé au voisinage du front, cette contribution a été abordée sur un problème modèle. Des conditions aux limites dites effectives ont été construites dans le cadre d'une méthode de décomposition de domaine puis testées pour un problème d'écoulement laminaire de convection naturelle sur parois rugueuses. Mˆeme si le problème traité reste encore éloigné des applications visées, cette contribution offre de nombreuses perspectives et constitue une première étape d'une modélisation multiéchelle des échanges pour la problématique de l'interaction corium-béton. Dans le cas plus complexe des écoulements multiphasiques multiconstituants et devant les difficultés expérimentales associées, le développement de lois de parois pour les outils existants aux échelles de description supérieures nécessite, au préalable, de disposer d'un outil de simulation numérique directe de l'écoulement au voisinage du front d'ablation. L'outil développé dans ce travail correspond à un modèle de Cahn-Hilliard/Navier-Stokes pour un mélange diphasique (liquide-gaz) compositionnel (corium-béton fondu) s'appuyant sur une description du système selon trois paramètres d'ordre associés respectivement aux fractions volumiques du gaz et aux deux espèces miscibles de la phase liquide ainsi que sur une décomposition de l'énergie libre selon une contribution diphasique et compositionnelle. Les équations de transport sont dérivées dans le cadre de la thermodynamique des processus irréversibles et résolues sur la base d'une application éléments finis de la plate-forme PELICANS. Plusieurs expériences numériques illustrent la validité et les potentialités d'application de cet outil sur des problèmes diphasiques et/ou compositionnels. Enfin, à partir de l'outil développé, nous abordons par simulation numérique directe une étude de la structure de la couche limite au voisinage du front d'ablation pour des bétons siliceux et silico-calcaire. / In the late phases of some scenario of hypothetical severe accident in Pressurized Water Reactors, a molten mixture of core and vessel structures, called corium, comes to interact with the concrete basemat. The safety numerical tools are lumped parameter codes. They are based on a large averaged description of heat and mass transfers which raises some uncertainties about the multi-scale description of the exchanges but also about the adopted boundary layer structure in the vicinity of the ablation front. In this context, the aim of this work is to tackle the problem of the boundary layer structure by means of direct numerical simulation. This work joins within the more general framework of a multi-scale description and a multi-scale modeling, namely from the local scale associated with the vicinity of the ablation front to the scale associated with the lumped parameter codes. Such a multi-scale description raises not only the problem of the local description of the multiphase multicomponent flow but also the problem of the upscaling between the local- and the macro-scale which is associated with the convective structures within the pool of corium. Here, we are particularly interested in the building of effective boundary conditions or wall laws for macro-scale models. The difficulty of the multiphase multicomponent problem at the local scale leads us to consider a relatively simplified problem. Effective boundary conditions are built in the frame of a domain decomposition method and numerical experiments are performed for a natural convection problem in a stamp shaped cavity to assess the validity of the proposed wall laws. Even if the treated problem is still far from the target applications, this contribution can be viewed as a first step of a multi-scale modeling of the exchanges for the molten core concrete issue. In the more complicated case of multiphase multicomponent flows, it is necessary to have a direct numerical simulation tool of the flow at the local scale to build wall laws for macro-scale models. Here, the developed tool corresponds to a Cahn-Hilliard/Navier-Stokes model for a two-phase compositional system. It relies on a description of the system by three volume fractions and on a free energy composed by a two-phase part and a compositional part. The governing equations are derived in the frame of the thermodynamic of irreversible processes. They are solved on the basis of a finite element application of the object-oriented software component library PELICANS. Several numerical experiments illustrate the validity and the potentialities of application of this tool on two-phase compositional problems. Finally, using the developed tool, we tackle by means of direct numerical simulation the problem boundary layer structure in the vicinity of the ablation front for limestone-sand and siliceous concretes.