Numerical schemes obtained from lattice Boltzmann equations for advection diffusion equations

Stability and accuracy of the numerical schemes obtained from the lattice Boltzmann equation (LBE) used for numerical solutions of two-dimensional advection-diffusion equations are presented. Three kinds of velocity models are used to determine the moving velocity of particles on a squre lattice. A system of explicit finite difference equations are derived from the LBE based on the Bhatnagar, Gross and Krook (BGK) model for individual velocity model. In order to approximate the advecting velocity field, a linear equilibrium distribution function is used for each of the moving directions. The stability regions of the schemes in the special case of the relaxation parameter omega in the LBE being set to omega=1 are found by analytically solving the eigenvalue problems of the amplification matrices corresponding to each scheme. As for the cases of general relaxation parameters, the eigenvalue problems are solved numerically. A benchmark problem is solved in order to investigate the relationship between the accuracy of the numerical schemes and the order of the Peclet number. The numerical experiments result in indicating that for the scheme based on a 9-velocity model we can find the parameters depending on the order of the given Peclet number, which generate accurate solutions in the stability region.