A wall function approach in lattice Boltzmann method: algorithm and validation using turbulent channel flow

Abstract

In the lattice Boltzmann method (LBM), the widely utilized wall boundary is the bounce-back (BB) boundary, which corresponds to the no-slip boundary. The BB boundary prevents the LBM from capturing the accurate shear drag on the wall when addressing high Reynolds number flows using coarse-grid systems. In this study, we proposed the "wall-function bounce (WFB)" boundary, a general framework to incorporate wall functions into the LBM's boundary condition, independent of specific information of discrete velocity schemes and collision functions. The WFB boundary calculates the appropriate shear drag on the wall using a wall function model, and thereafter just modifies partial diagonal distribution functions to reflect the shear drag. The Spalding's law was utilized as the wall function in WFB. Simulations of turbulent channel flow at Re_τ=640 and 2003 using the LBM-based large-eddy simulation (LBM-LES) were conducted to validate the effectiveness of the proposed boundary condition. The results indicate that the BB boundary underestimated the time-averaged velocity in the buffer layer at Re_τ=640, and the averaged velocity in the entire domain at Re_τ=2003, when using coarse-grid systems. However, WFB obtained the proper shear drag on the wall and thus, compensated for the underestimation and agreed better with the experimental or direct numerical simulation data, especially at the first-layer grid. In addition, WFB improved the Reynolds normal stress in the near-wall region to some extent. The distributions of shear stress on the wall by WFB were analogous to those by the wall model function in the finite volume method.