Mathematical Modeling and Associated Numerical Simulation of Fusion/Solidification Front Evolution in the Context of Severe Accident of Nuclear Power Engineering

Adrien Drouillet, Guillaume Bois, Romain Le Tellier, Raphaël Loubère and Mathieu Peybernes
Additional contact information
Adrien Drouillet: Institut de Mathématiques de Bordeaux, University of Bordeaux, CNRS, UMR 5251, 33400 Talence, France
Guillaume Bois: CEA-DES, Université Paris—Saclay, 91191 Gif-sur-Yvette, France
Romain Le Tellier: CEA-DES-IRESNE-DTN, CEA Cadarache, 13108 Saint Paul-lez-Durance, France
Raphaël Loubère: Institut de Mathématiques de Bordeaux, University of Bordeaux, CNRS, UMR 5251, 33400 Talence, France
Mathieu Peybernes: CEA-DES-IRESNE-DTN, CEA Cadarache, 13108 Saint Paul-lez-Durance, France

Mathematics, 2021, vol. 10, issue 1, 1-23

Abstract: Considering transient processes where liquid/solid phase change occurs, this paper focuses on the associated modeling and numerical treatment in the frame of “Computational Fluid Dynamics” simulations. While being of importance in many industrial applications involving solidification and melting of mixed materials, including power and manufacturing engineering, the first application of this work pertains to the analysis of severe accidents in a nuclear reactor. Indeed, in this context, the molten core materials (a.k.a. corium) can form a high-temperature multiphase liquid pool at the boundary of which fusion and solidification phenomena are of prime importance. In this context, even if materials at play are treated as pure components, it is mandatory to distinguish two different phase change temperatures with a solid fusion temperature and a liquid solidification temperature. Accordingly, in the frame of a sharp interface representation, the paper introduces non-classical heterogeneous conditions at the liquid/solid boundary in such a way that both moving interface (through Stefan conditions associated with fusion or solidification) and static interface (imposing heat flux continuity) are supported at the same time on different spatial locations along this boundary. Within a monolithic resolution of Navier–Stokes and heat conduction equations, this interface is explicitly tracked with combined Front-Tracking and VOF methods. In order to ensure zero velocity in the solid phase, an Immersed Boundary Method and a direct forcing penalization are also introduced. The main relevant features of this combination of numerical methods are discussed along with their implementation in the TrioCFD code taking advantage of the pre-existing code capabilities. Numerical simulations including both verification tests and a case of interest for our industrial application are reported and demonstrate the applicability of the proposed triptych model+methods+code to treat such problems. The numerical tools and the simulation code developed in this work could be used not only in the several accident context but also to simulate melting, solidification and fusion processes occurring in aerodynamics, hypersonic reentry vehicles and laser applications to cite but a few.

Keywords: simulation of front; fusion; solidification; corium crust; 2D heat conduction; phase changes; Navier–Stokes; modeling (search for similar items in EconPapers)
JEL-codes: C (search for similar items in EconPapers)
Date: 2021
References: View complete reference list from CitEc
Citations:

Downloads: (external link)
https://www.mdpi.com/2227-7390/10/1/116/pdf (application/pdf)
https://www.mdpi.com/2227-7390/10/1/116/ (text/html)

Related works:
This item may be available elsewhere in EconPapers: Search for items with the same title.

Export reference: BibTeX RIS (EndNote, ProCite, RefMan) HTML/Text

Persistent link: https://EconPapers.repec.org/RePEc:gam:jmathe:v:10:y:2021:i:1:p:116-:d:715384

Access Statistics for this article

Mathematics is currently edited by Ms. Emma He

More articles in Mathematics from MDPI
Bibliographic data for series maintained by MDPI Indexing Manager ().