Equivalent finite difference and partial differential equations for the lattice Boltzmann method

A general method for the derivation of equivalent finite difference equations (EFDEs) and subsequent equivalent partial differential equations (EPDEs) is presented for a general matrix lattice Boltzmann method (LBM). The method can be used for both the advection diffusion equations and Navier-Stokes equations in all dimensions. In principle, the EFDEs are derived using a recurrence formula. A computational algorithm is proposed for generating sequences of matrices and vectors that are used to obtain EFDEs coefficients. For all DdQq velocity models, the algorithm is proven to be finite and all coefficients are obtained after q iterations. The resulting EFDEs and EPDEs are derived for selected velocity models and include the single relaxation time, multiple relaxation times, and cascaded LBM collision operators. The algorithm for the derivation of EFDEs and EPDEs is implemented in C++ using the GiNaC library for symbolic algebraic computations. Its implementation is available under the terms and conditions of the GNU general public license (GPL).