A second order turbulence model based on a Reynolds stress approach for two-phase boiling flow and application to fuel assembly analysis

High-thermal performance PWR (pressurized water reactor) spacer grids require both low-pressure loss and high critical heat flux (CHF) properties. Numerical investigations on the effect of angles and position of mixing vanes and to understand in more details the main physical phenomena (wall boiling, entrainment of bubbles in the wakes, recondensation) are required. In the field of fuel assembly analysis or design by means of CFD codes, the overwhelming majority of the studies are carried out using two-equation Eddy Viscosity Models (EVM), especially the standard K– ɛ model, while the use of Reynolds Stress Transport Models (RSTM) remains exceptional. The simulation of swirling flow generated by the mixing vanes plays an important role for the prediction of the CHF for the fuel assemblies. For this reason, according to Mimouni et al. (2008b, 2009b), rotation effects and RSTM model are more specifically addressed in the paper. Before comparing performance of EVM and RSTM models on fuel assembly geometry, we performed calculations with simpler geometries, the DEBORA case and the ASU-annular channel case. ASU-annular channel case has already been addressed in Mimouni et al. (2008b, 2009b). Then, a geometry closer to actual fuel assemblies is considered. It consists of a rectangular test section in which a 2 × 2 rod bundle equipped with a simple spacer grid with mixing vanes is inserted. The influence of the turbulence model on target variables linked to CHF limitation will be discussed. Moreover, the sensitivity to the mesh refinement will be particularly examined. The study of this case is a further step towards the modelling of the two-phase boiling flow in real-life grids and rod bundles.