Lattice-boltzmann simulations of the convection-diffusion equation with different reactive boundary conditions

We have tested the accuracy and stability of lattice-Boltzmann (LB) simulations of the convection-diffusion equation in a two-dimensional channel flow with reactive-flux boundary conditions. We compared several different implementations of a zero-concentration boundary condition using the Two-Relaxation-Time (TRT) LB model. We found that simulations using an interpolation of the equilibrium distribution were more stable than those based on Multi-Reflection (MR) boundary conditions. We have extended the interpolation method to include mixed boundary conditions, and tested the accuracy and stability of the simulations over a range of Damköhler and Péclet numbers.