In recent years, the lattice Boltzmann method (LBM) has be come an alternative and robust method for the simulation of complexfluid flows. Unlike the traditional computational fluid dynamics(CFD) methods, which are based on the discretization of macroscopic continuum equations, the LBM is based on the microscopicmodels and mesoscopic kinetic equations. In this method, the macroscopic dynamics of a fluid is obtained from the collective behavior of many microscopic particles in the system. It has been proved that the LBM can recover the N–S equations by using the Chapman Enskog expansion [14]. The key advantages of the LBM include the explicit feature of the governing equation, easy for parallel computation, and simple implementation of boundary conditions on curved boundaries. The LBM is derived from the kinetic theory of gases, which introduces some significant features distinguishing it from the traditional numerical methods. The convection operator of the LBM is linear, which is combined with are laxation process to recover the nonlinear macroscopic advection. In addition, the in compressible N–S equations can be obtained in the nearly in compressible limit of the LBM. The pressure of the LBM is calculated using an equation of state rather than by solving the Poisson equation. The LBM comprises two main processes,namely, streaming and collision. In general, the collision operator is linearized by Bhatnagar–Gross–Krook (BGK) approximation.