Direct numerical simulation of turbulent mixing in regular and fractal grid turbulence

Abstract

Turbulent mixing in regular and fractal grid turbulence is investigated in this work by using direct numerical simulation (DNS). Two types of turbulence-generating grids are used: a biplane square grid (regular grid) and a square fractal grid. The thickness ratios tr of the fractal grids are set at 5.0 and 8.5. The grid solidity is maintained at σ=0.36 for all the grids. The mesh Reynolds number, ReM=U0Meff/ν, is set at 2500 for all cases, where U0 is the cross-sectionally averaged mean velocity; Meff, the effective mesh size; and ν, the kinematic viscosity. The grids are numerically generated using the immersed boundary method at 4Meff downstream of the entrance to the computational domain. The computational domain size normalized by Meff is 64×8×8 in the streamwise, vertical and spanwise directions for the regular grid and 64×16×16 for the fractal grids. Scalar mixing layers that initially have a step profile develop downstream of the grids. The Prandtl number is set at Pr=0.71 considering the heat transfer in air flow. Instantaneous temperature fields, instantaneous fluctuating temperature fields and fundamental turbulent statistics are presented. The results show that turbulent mixing is more strongly enhanced in fractal grid turbulence than in regular grid turbulence for the same ReM. In fractal grid turbulence, turbulent mixing is more strongly enhanced at tr=8.5 than at tr=5.0.