On recovering the second-order convergence of the lattice Boltzmann method with reaction-type source terms

This study analyses a method to consistently recover the second-order convergence of the lattice Boltzmann method (LBM), which is frequently degraded by the improper discretisation of required source terms. The work focuses on advection-diffusion models in which the source terms are dependent on the intensity of transported fields. The main findings are applicable to a wide range of formulations within the LBM framework. All considered source terms are interpreted as contributions to the zeroth-moment of the distribution function. These account for sources in a scalar field, such as density, concentration or temperature. In addition to this, certain immersed boundary methods can be interpreted as a source term in their formulation, highlighting a further application for this work. This paper is dedicated to three aspects regarding proper inclusion of the source term in LBM schemes. Firstly, it identifies the differences observed between the ways source terms are included in the LBM schemes present in the literature. The algebraic manipulations are explicitly presented in this paper to clarify the differences observed, and identify their origin. Secondly, it analyses in full detail, the implicit relation between the value of the transported macroscopic field, and the sum of the LBM densities. Moreover, three equivalent forms of the collision operator are presented. Finally, closed-form solutions of this implicit relation are shown for a variety of common models. The second-order convergence of the proposed LBM schemes is verified on both linear and non-linear source terms. The pitfalls of the commonly used acoustic and diffusive scaling are identified and discussed. Moreover, for a simplified case, the competing errors are shown visually with isolines of error in the space of spatial and temporal resolutions.