  Understanding the coolant thermal hydraulics in rod bundles is essential to
the design of nuclear reactors. However, flows with low Reynolds numbers
present serious modeling challenges, especially in heat transfer and natural
convection. They are difficult to analyze through standard Computational Fluid
Dynamics (CFD) tools. High-fidelity simulations, such as Direct Numerical
Simulations (DNS), can provide invaluable insight into flow physics, supporting
experiments in developing a deeper understanding and eventually enabling the
accurate simulation of this class of flows. Data generated from these
high-fidelity methods can then be used to benchmark available turbulence models
and deliver cheap, faster running methods. In the present work, the convective
heat transfer in hexagonal subchannels was studied through a DNS approach,
using the high-order spectral element method code Nek5000, developed at Argonne
National Laboratory. First, the geometric model composed of two hexagonal
arrayed rod bundle subchannels with a pitch-to-diameter ratio of 1.5 is built,
and then, the mesh is generated. These unusually high P/D and low Reynolds
numbers represent conditions of interest for gas-fast reactors (GFRs). To our
knowledge, there is no available dataset in these conditions. In this work, we
detailed the development of the numerical benchmark and a series of preliminary
LES simulations. Periodic boundary conditions are applied in the streamwise and
spanwise directions and non-slip boundary conditions at the wall. Four cases
are studied, with Reynolds of Re=2500, 5000, 7500, and 10000. All calculations
have been performed with the Prandtl number of 0.61, corresponding to helium
conditions of interest for Gas Fast Reactor applications. The results are
analyzed with a polynomial-order convergence study, and the Reynolds stresses,
and the turbulent kinetic budgets are presented and discussed.
