Analysis of dynamical flow structure in a square arrayed rod bundle

Large eddy simulation (LES) of turbulent flow in a bare rod bundle was performed, and a new concept about the flow structure that enhances heat transport between subchannels was proposed. To investigate the geometrical effect, the LES was performed for three different values of rod diameter over pitch ratio (D/P = 0.7, 0.8, 0.9). The computational domain containing 4 subchannels was large enough to capture large-scale structures wide across subchannels. Lateral flow obtained was unconfined in a subchannel, and some flows indicated a pulsation through the rod gap between subchannels. The gap flow became strong as D/P increased, as existing experimental studies had reported. Turbulence intensity profile in the rod gap suggested that the pulsation was caused by the turbulence energy transferred from the main flow to the wall-tangential direction. This implied that the flow pulsation was an unsteady mode of the secondary flow and arose from the same geometrical effect of turbulence. This implication was supported by the analysis results: two-points correlation functions of fluctuating velocities indicated two length-scales, P-D and D, respectively of cross-sectional and longitudinal motions; turbulence stress in the cross-sectional mean flow contained a non-potential component, which represented energy injection through the unsteady longitudinal fluid motion.