Numerical study of the stratification erosion benchmark for NPPs containment using CFD code GASFLOW-MPI

During a severe accident involving a Nuclear Power Plant (NPP), a hydrogen-enriched stratification may occur in the dome of the containment. This is due to the gaseous density difference of a steam/hydrogen turbulent jet released from the break. This accident scenario is one of the high-ranking phenomena identified in the Phenomena Identification and Ranking Table (PIRT) for containment safety thermal hydraulics. In this work, the temporal evolution of the stratification erosion is analysed using the parallel CFD code GASFLOW-MPI with comparison to the experimental data from the THAI TH-20 benchmark, where helium is used as a hydrogen simulate. Identifying a suitable turbulence model is an important issue for a successful stratification erosion simulation due to turbulent mixing between the stratification and turbulent jet featured with the sharp velocity and density gradients. Three well-known turbulence models: 1. The classic k-e model, 2. Large Eddy Simulation (LES) model and 3. Detached Eddy Simulation (DES) model are employed to evaluate their performances. Due to the negative buoyancy effect resulting from the light gas layer, the jet gradually loses its initial momentum and cannot penetrate the stratified layer in depth. The interaction Froude (Fr) number which is used to categorize the mixing mechanism is close to 1 in this work, indicating that the erosion process of the stratified layer is dominated by both the molecular diffusion and the momentum transport. The experimental data measured at different heights in THAI TH-20 benchmark, from the pure air region at the lower containment region up into the large helium concentration gradient region, and further up to high helium concentration region at the top of the containment. The calculated results show that the LES and DES turbulence models can capture time histories of helium concentration at all of the sensor points accurately, while the k-e turbulence model predictions are always delayed due to its overestimating the turbulence mixing in the stratified layer. Moreover, compared with the k-e model, more detailed turbulent eddy structure is resolved by the LES model and DES model, whose frequency satisfies the negative −5/3 energy decay rate for a wide region of frequency. Therefore, the LES and DES models are recommended for the stratification erosion simulation. This can provide guidance for nuclear engineers analysing hydrogen distribution phenomena during severe accidents.