Analysis of hemispherical vessel ablation failure involving natural convection by MPS method with corrective matrix

In a severe accident of a light water reactor, the reactor pressure vessel (RPV) lower head may fail due to ablation at the vessel wall boundary involving natural convection of molten core materials. Accurate prediction of RPV lower head failure is essential for assessing severe accident progression and improving accident management because it greatly influences the subsequent ex-vessel accident progressions. However, there have been still large uncertainties about RPV lower head failure mode in the Fukushima Daiichi Nuclear Accident in 2011. The Lagrangian based MPS (moving particle semi-implicit) method has advantage of analyzing such phenomena involving complex interfaces and liquid-solid phase changes over other Eulerian mesh-based method. In the preceding study, small-scale Pb–Bi hemisphere vessel ablation experiment, with silicone oil as simulated molten core, was reproduced qualitatively by original MPS method. However, ablation mechanism associated with natural convection of the high temperature liquid could not be discussed because of significant influence of numerical discretizing error. In this study, the improved MPS method coupling corrective matrix in the particle interaction model which largely suppress the numerical fluctuation was adopted to analyze the experiment. The results show that the ablated metal relocation may enhance convective heat transfer in the downstream. As a result, ablation of the vessel wall extends from the level, close to the silicone oil surface down to the bottom of the vessel rather than previously simulated localized ablation near the silicone oil surface.