Lattice Boltzmann method with two relaxation times for advection–diffusion equation: Third order analysis and stability analysis

The objectives of this study are to investigate the third order accuracy and linear stability of the lattice Boltzmann method (LBM) with the two-relaxation-time collision operator (LTRT) for the advection–diffusion equation (ADE) and compare the LTRT model with the single-relaxation-time (LBGK) model. While the LBGK has been used extensively, the LTRT appears to be a more flexible model because it uses two relaxation times. The extra relaxation time can be used to improve solution accuracy and/or stability. This study conducts a third order Chapman–Enskog expansion on the LTRT to recover the macroscopic differential equations up to the third order. The dependency of third order terms on the relaxation times is obtained for different types of equilibrium distribution functions (EDFs) and lattices. By selecting proper relaxation times, the numerical dispersion can be significantly reduced. Furthermore, to improve solution accuracy, this study introduces pseudo-velocities to develop new EDFs to reduce the second order numerical diffusion. This study also derives stability domains based on the lattice Peclet number and Courant number for different types of lattices, EDFs and different values of relaxation times, while conducting linear stability analysis on the LTRT. Numerical examples demonstrate the improvement of the LTRT solution accuracy and stability by selecting proper relaxation times, lattice Peclet number and Courant number.