Three-dimensional cascaded lattice Boltzmann method: Improved implementation and consistent forcing scheme

The cascaded or central-moment-based lattice Boltzmann method (CLBM) proposed in Phys. Rev. E 73, 066705 (2006)PLEEE81539-375510.1103/PhysRevE.73.066705 possesses very good numerical stability. However, two constraints exist in three-dimensional (3D) CLBM simulations. First, the conventional implementation for 3D CLBM involves cumbersome operations and requires much higher computational cost compared to the single-relaxation-time (SRT) LBM. Second, it is a challenge to accurately incorporate a general force field into the 3D CLBM. In this paper, we present an improved method to implement CLBM in 3D. The main strategy is to adopt a simplified central moment set and carry out the central-moment-based collision operator based on a general multi-relaxation-time (GMRT) framework. Next, the recently proposed consistent forcing scheme for CLBM Fei and Luo, Phys. Rev. E 96, 053307 (2017)2470-004510.1103/PhysRevE.96.053307 is extended to incorporate a general force field into 3D CLBM. Compared with the recently developed nonorthogonal CLBM Rosis, Phys. Rev. E 95, 013310 (2017)2470-004510.1103/PhysRevE.95.013310, our implementation is proved to reduce the computational cost significantly. The inconsistency of adopting the discrete equilibrium distribution functions in the nonorthogonal CLBM is analyzed and validated. The 3D CLBM developed here in conjunction with the consistent forcing scheme is verified through numerical simulations of several canonical force-driven flows, highlighting very good properties in terms of accuracy, convergence, and consistency with the nonslip rule. Finally, the techniques developed here for 3D CLBM can be applied to make the implementation and execution of 3D MRT-LBM more efficient.