Numerics of the lattice Boltzmann method: Effects of collision models on the lattice Boltzmann simulations

We conduct a comparative study to evaluate several lattice Boltzmann (LB) models for solving the near incompressible Navier-Stokes equations, including the lattice Boltzmann equation with the multiple-relaxation-time (MRT), the two-relaxation-time (TRT), the single-relaxation-time (SRT) collision models, and the entropic lattice Boltzmann equation (ELBE). The lid-driven square cavity flow in two dimensions is used as a benchmark test. Our results demonstrate that the ELBE does not improve the numerical stability of the SRT or the lattice Bhatnagar-Gross-Krook (LBGK) model. Our results also show that the MRT and TRT LB models are superior to the ELBE and LBGK models in terms of accuracy, stability, and computational efficiency and that the ELBE scheme is the most inferior among the LB models tested in this study, thus is unfit for carrying out numerical simulations in practice. Our study suggests that, to optimize the accuracy, stability, and efficiency in the MRT model, it requires at least three independently adjustable relaxation rates: one for the shear viscosity nu (or the Reynolds number Re), one for the bulk viscosity zeta, and one to satisfy the criterion imposed by the Dirichlet boundary conditions which are realized by the bounce-back-type boundary conditions.