Toward Direct Numerical Simulation of Turbulent and Transitional Flow in Hexagonal Subchannels for Helium Conditions

TL;DR

This paper presents high-fidelity DNS simulations of turbulent and transitional helium flows in hexagonal subchannels, providing new benchmark data for low Reynolds number conditions relevant to gas-fast reactors.

Contribution

The study develops the first DNS benchmark dataset for low Reynolds number helium flows in hexagonal subchannels, aiding turbulence model validation and reactor coolant analysis.

Findings

Reynolds stresses and turbulent kinetic budgets characterized

Polynomial-order convergence confirmed accuracy of simulations

Benchmark data for Re=2500 to 10000 provided

Abstract

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.