Integral analysis of debris material and heat transport in reactor vessel lower plenum

Abstract

An integral, fast-running, two-region model has been developed to characterize the debris material and heat transport in the reactor lower plenum under severe accident conditions. The debris bed is segregated into the oxidic pool and an overlying metallic layer. Debris crusts can develop on three surfaces: the top of the molten pool, the RPV wall, and the internal structures. To account for the decay heat generation, the crust temperature profile is assumed to be parabolic. The oxidic debris pool is homogeneously mixed and has the same material composition, and hence the same thermophysical properties, as the crusts, while the metallic constituents are assumed to rise to the top of the debris pool. Steady-state relationships are used to describe the heat transfer rates, with the assessment of solid or liquid state, and the liquid superheat in the pool being based on the average debris temperature. Natural convection heat transfer from the molten debris pool to the upper, lower and embedded crusts is calculated based on the pool Rayleigh number with the conduction heat transfer from the crusts being determined by the crust temperature profile. The downward heat flux is transferred to the lowest part of the RPV lower head through a crust-to-RPV contact resistance. The sideward heat flux is transferred to the upper regions of the RPV lower head as well as to the internal structures. The upward heat flux goes to the metal layer, water, or available heat sink structures above. Quenching due to water ingression is modeled separately from the energy transfer through the crust. The RPV wall temperature distribution and the primary system pressure are utilized to estimate challenges to the RPV integrity. Should the RPV be submerged, the heat removal is enhanced by the ex-vessel cooling due to nucleate boiling. The convection heat transfer correlations for the molten debris pool were validated against available experimental data and theoretical predictions. Testing of the model for a range of conditions in a PWR lower plenum produced consistent results. In addition, a comparison of the integral approach to a more detailed, special purpose model showed good agreement.