Direct numerical simulation and large eddy simulation on a turbulent wall-bounded flow using lattice Boltzmann method and multiple GPUs

Direct numerical simulation (DNS) and large eddy simulation (LES) were performed on the wall-bounded flow at Re τ = 180 using lattice Boltzmann method (LBM) and multiple GPUs (Graphic Processing Units). In the DNS, 8 K20M GPUs were adopted. The maximum number of meshes is 6.7 × 10 7, which results in the nondimensional mesh size of Δ + = 1.41 for the whole solution domain. It took 24 hours for GPU-LBM solver to simulate 3 × 10 6 LBM steps. The aspect ratio of resolution domain was tested to obtain accurate results for DNS. As a result, both the mean velocity and turbulent variables, such as Reynolds stress and velocity fluctuations, perfectly agree with the results of Kim et al. (1987) when the aspect ratios in streamwise and spanwise directions are 8 and 2, respectively. As for the LES, the local grid refinement technique was tested and then used. Using 1.76 × 10 6 grids and Smagorinsky constant (C s) = 0.13, good results were obtained. The ability and validity of LBM on simulating turbulent flow were verified.