Numerical simulation of subcooled flow boiling for nuclear engineering applications using OpenFOAM

Abstract This work is focused on the development and validation of models and methods for the simulation of wall boiling in nuclear engineering applications with the computational fluid dynamics (CFD) code OpenFOAM. The new chtMultiRegionReactingTwoPhaseEulerFoam solver was developed based on the reactingTwoPhaseEulerFoam solver of OpenFOAM Foundation version 7. The solver is used for the simulation of two-phase flow under consideration of wall boiling and conjugate heat transfer (CHT) between solid structure and two-phase fluid regions. The Euler–Euler approach for two-phase flows was used. The heat flux during wall boiling was calculated with the help of the extended Rensselaer Polytechnic Institute wall heat flux partitioning model, in which the convective heat flux between solid wall and two-phase flow with high void fractions was also considered. The solver was validated against experimental data from the OECD/NEA PWR Subchannel and Bundle Tests benchmark. This Nuclear Power Energy Corporation (NUPEC) database provides data for different fuel assembly subchannel geometries at different thermal-hydraulic conditions. 10 experimental runs with different boundary conditions of the benchmark exercise I-1 were simulated with the chtMultiRegionReactingTwoPhaseEulerFoam solver. The solver showed good numerical stability in all examined cases, which captured different boiling regimes with up to cross-section averaged void fractions of 0.6. The results were compared with measured data for the averaged over the cross-section of the investigated geometry void fractions. Good agreement with experimental data was observed.