The Lattice Boltzmann Method (LBM), a weakly compressible approach, has attracted considerable attention for turbulent flow simulations but remains limited in high-Reynolds-number applications due to its reliance on uniform meshes. This study investigates the interpolation-based Lattice Boltzmann Method (IBLBM) for direct numerical simulation (DNS) of turbulent channel and duct flows using nonuniform grids. Simulations are conducted for turbulent Poiseuille channel flows at \(\boldsymbol{Re_{\tau} = 180, 395, 640},\) and 1000, as well as for turbulent duct flow at \(\boldsymbol{Re_{\tau} = 360}\) . The predicted mean velocity profiles, Reynolds stresses, vorticity distributions, and turbulent kinetic energy budgets show excellent agreement with benchmark DNS data. The IBLBM scheme accurately resolves near-wall turbulence and secondary flow structures in ducts while maintaining similar computational overhead across different interpolation orders. For the 7- and 200-million-grid cases at \(\boldsymbol{Re_{\tau} = 180}\) and 640, LBM with \(\boldsymbol{\Delta^+ \sim 3}\) is approximately 1.9 and 3.7 times faster, respectively, than sixth-order IBLBM; however, to match IBLBMs wall resolution ( \(\boldsymbol{\Delta^+ \sim 0.3}\) ) using \(\boldsymbol{\Delta^+ \sim 1}\) , LBM becomes four to six times slower. These results demonstrate the robustness and accuracy of IBLBM for high-Reynolds-number wall-bounded turbulence, offering an optimal balance between computational efficiency and physical fidelity through the use of nonuniform grids.
