Modélisation de l'interaction entre le cœur fondu d'un réacteur à eau pressurisée et le radier en béton du bâtiment réacteur

Modelling of the Molten Core Concrete Interaction (MCCI)

• Auteur / Autrice : Mathieu Guillaumé
• Direction : Hervé Combeau, Jean-Marie Seiler
• Type : Thèse de doctorat
• Discipline(s) : Mécanique et énergétique
• Date : Soutenance le 12/12/2008
• Etablissement(s) : Vandoeuvre-les-Nancy, INPL
• Ecole(s) doctorale(s) : EMMA - Ecole Doctorale Energie - Mécanique - Matériaux
• Partenaire(s) de recherche : Laboratoire : Laboratoire de science et génie des matériaux et de métallurgie (Nancy)
• Jury : Président / Présidente : Lounès Tadrist
• Jury : Examinateurs / Examinatrices : Hervé Combeau, Jean-Marie Seiler, Lounès Tadrist, Florian Fichot, Dominique Gobin, Kresna Atkhen
• Jury : Rapporteur / Rapporteuse : Florian Fichot, Dominique Gobin

Mots clés

Mots clés contrôlés

Fission nucléaire

Centrales nucléaires -- Mesures de sécurité

Mots clés libres

Sûreté nucléaire

Modélisation

Interaction corium- béton

Nuclear safety

Modelling

Molten core concrete interaction

Résumé

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