Direct numerical simulation of flow and heat transport in a closely-spaced bare rod bundle

Abstract

Direct Numerical Simulation (DNS) of fully-developed flow and heat transfer is performed for a bare rod bundle with a low pitch-to-diameter ratio using the highly-scalabale spectral element code Nek5000. The simulation is performed at a Reynolds number, Reh, based on the bulk velocity and hydraulic diameter, of 9800, with iso-temperature and iso-flux thermal boundary conditions and Prandtl numbers 0.025, 1.0 and 2.0. The mean velocity and temperature statistics are first validated against the experimental measurements of Hooper et al. (1983) and complementary DNS results of Lai et al. (2019) which are reported in the literature at a higher Reh of 22,600. The mean flow and temperature are demonstrated a similar behaviour at the present low Renolds number and the higher Reynolds number of Lai et al. (2019), indicating the Reynolds number scaling of the mean flow. The anisotropy of the flow is illustrated using distribution of normal and Reynolds stresses, and the Lumley invariant map. It is shown that the turbulence in the subchannel region is characteristic of wall-bounded turbulent flows such as channel flows. In the narrow gap, however, turbulence is suppressed, only enhanced in the streamwise direction due to flow pulsations. Further, temperature statistics with both thermal boundary conditions are presented, demonstrating the effect of the Prandtl number on the mean and fluctuating quantities. The temperature statistics with iso-temperature wall boundary conditions, not reported in literature before, are shown to behave differently in the interstitial subchannel region and the narrow gap between two rods. The distribution of turbulent heat flux is also presented. The streamwise component of turbulent heat flux is again shown to be enhanced, in comparison to the same in the subchannel. The pulsating motion of flow across the narrow gap is discussed further. A frequency analysis of the velocity in the narrow gap elucidates the frequency of flow pulsations. Unlike mean flow statistics, the flow pulsations are shown to exhibit a low Reynolds number effect. The reference DNS data presented in the present work, complementary to that of Lai et al. (2019), serves as a benchmark against which lower order turbulence models may be further developed or improved. In particular, the newly presented anisotropy data and the turbulence heat flux data are essential for further development, improvement and validation of turbulent heat transfer models at different Prandtl numbers.