Transient heat transfer and crust evolution during debris bed melting process in the hypothetical severe accident of HPR1000

In the late in-vessel phase of a nuclear reactor severe accident, the internal heat transfer and crust evolution during the debris bed melting process have important effects on the thermal load distribution along the vessel wall, and further affect the reactor pressure vessel (RPV) failure mode and the state of melt during leakage. This study coupled the phase change model and large eddy simulation to investigate the variations of the temperature, melt liquid fraction, crust and heat flux distributions during the debris bed melting process in the hypothetical severe accident of HPR1000. The results indicated that the heat flow towards the vessel wall and upper surface were similar at the beginning stage of debris melting, but the upward heat flow increased significantly as the development of the molten pool. The maximum heat flux towards the vessel wall reached 0.4 MW/m2. The thickness of lower crust decreased as the debris melting. It was much thicker at the bottom region with the azimuthal angle below 20° and decreased rapidly at the azimuthal angle around 20–50°. The maximum and minimum thicknesses were 2 and 90 mm, respectively. By contrast, the distribution of upper crust was uniform and reached stable state much earlier than the lower crust, with the thickness of about 10 mm. Moreover, the sensitivity analysis of initial condition indicated that as the decrease of time interval from reactor scram to debris bed dried-out, the maximum debris temperature and melt fraction became larger, the lower crust thickness became thinner, but the upper crust had no significant change. The sensitivity analysis of in-vessel retention (IVR) strategies indicated that the passive and active external reactor vessel cooling (ERVC) had little effect on the internal heat transfer and crust evolution. In the case not considering the internal reactor vessel cooling (IRVC), the upper crust was not obvious.