CFD simulation of corium flow through an end fitting of a pressurised heavy water reactor

Abstract

This study investigates the numerical simulation of corium relocation during a hypothetical accident, focussing on the modelling of core melt, its solidification and redistribution. This paper presents results of a high-fidelity computational fluid dynamics (CFD) simulation of an ex-vessel corium flow through a horizontal end fitting geometry, a pressurized heavy water reactor (PHWR) relevant configuration. The objective of this work is to simulate the corium progression within the end fitting of a PHWR to assess the suitability of two viscosity models using a CFD code; Siemens STAR-CCM+. This first-of-a-kind, exploratory study demonstrates the potential of the existing CFD models for modeling corium spreading within the end fitting geometry. Mixed convection turbulent corium flow in the end-fitting has been modeled in the three-dimensional CFD simulations. Transient simulations using the volume-of-fluid approach and melting-solidification models were performed to predict the corium melt spread within the end fitting. Two suitable temperature-dependant models for corium viscosity were tested, one being empirical viscosity model by Ramacciotti and the other mechanistic viscosity model available in STAR-CCM+. The unique flow features such as the Rayleigh-Bénard instability could be observed using the Ramacciotti empirical viscosity modelling approach. A qualitative assessment of the results was undertaken, showing that the CFD model has predicted the corium spreading behaviour that is consistent with the expectations based on engineering judgement and experience. This kind of simulations are expected to serve as a guidance in planning lab-scale experiments and aid in CFD-informed modeling improvements for integral (or system) codes going forward.