Coupled neutronic-thermal-hydraulics analysis in a coolant subchannel of a PWR using CFD techniques

The high capacity of Computational Fluid Dynamics code to predict multi-dimensional thermal-hydraulics behaviour and the increased availability of capable computer systems are making that method a good tool to simulate phenomena of thermal-hydraulics nature in nuclear reactors. However, since there are no neutron kinetics models available in commercial CFD codes to the present day, the application of CFD in the nuclear reactor safety analyses is still limited. The present work proposes the implementation of the point kinetics model (PKM) in ANSYS - Fluent to predict the neutronic behaviour in a Westinghouse Sequoyah nuclear reactor, coupling with the phenomena of heat conduction in the rod and thermal-hydraulics in the cooling fluid, via the reactivity feedback. Firstly, a mesh convergence and turbulence model study was performed, using the Reynolds-Average Navier-Stokes method, with square arrayed rod bundle featuring pitch to diameter ratio of 1:32. Secondly, simulations using the k-! SST turbulence model were performed with an axial distribution of the power generation in the fuel to analyse the heat transfer through the gap and cladding, and its in fluence on the thermal-hydraulics behaviour of the cooling fluid. The wall shear stress distribution for the centre-line rods and the dimensionless velocity were evaluated to validate the model, as well as the in fluence of the mass flow rate variation on the friction factor. The coupled model enabled to perform a dynamic analysis of the nuclear reactor during events of insertion of reactivity and shutdown of primary coolant pumps. (author)