Multi-scale simulation of microstructural evolution for molten U–Zr–O ternary system

The morphological evolution of the molten corium is critical to the success of the “in-vessel retention” (IVR) strategy, which is of great significance for the safety of the nuclear industry under severe accidents. However, the direct observation of microstructural evolution by experiments is impossible in extreme conditions. Herein, a multi-scale simulation approach was employed to investigate the microstructural evolution and related mechanisms for the molten U–Zr–O ternary system. First, based on the reported calculated phase diagram, a simplified free energy formula describing a phase diagram of U–Zr–O was determined, where only liquid phases are focused. Then, the atomic diffusion coefficients were calculated by ab-initio molecular dynamics simulations. Based on such information, a homemade phase-field model was developed, which was coupled with Cahn–Hillard and Navier–Stokes equations. Accordingly, the microstructural evolution of molten U–Zr–O materials was systematically investigated. The results reveal that the evolution of the molten core was driven by both the decreasing free energy and gravity potential, accessing the final morphology of the system as a two-layer structure, including a metallic layer and an oxidized layer. Moreover, diverse intermediate morphologies could be observed from different initial component distributions, although their average compositions were identical. Finally, some intermediate simulated microstructural patterns agreed well with experimental observation, demonstrating the unstable structures captured by experiments and the validity of the simulations. The present study provides an effective approach to investigate the microstructural evolution of molten core and a guideline for optimizing the IVR strategy based on different morphologies in the corium.