English

Constructing relaxation systems for lattice Boltzmann methods

Numerical Analysis 2022-09-14 v3 Numerical Analysis Analysis of PDEs

Abstract

We present the first top-down ansatz for constructing lattice Boltzmann methods (LBM) in d dimensions. In particular, we construct a relaxation system (RS) for a given scalar, linear, d-dimensional advection-diffusion equation. Subsequently, the RS is linked to a d-dimensional discrete velocity Boltzmann model (DVBM) on the zeroth and first energy shell. Algebraic characterizations of the equilibrium, the moment space, and the collision operator are carried out. Further, a closed equation form of the RS expresses the added relaxation terms as prefactored higher order derivatives of the conserved quantity. Here, a generalized (2d+1)x(2d+1) RS is linked to a DdQ(2d+1) DVBM which, upon complete discretization, yields an LBM with second order accuracy in space and time. A rigorous convergence result for arbitrary scaling of the RS, the DVBM and conclusively also for the final LBM is proven. The top-down constructed LBM is numerically tested on multiple GPUs with smooth and non-smooth initial data in d=3 dimensions for several grid-normalized non-dimensional numbers.

Keywords

Cite

@article{arxiv.2208.14976,
  title  = {Constructing relaxation systems for lattice Boltzmann methods},
  author = {Stephan Simonis and Martin Frank and Mathias J. Krause},
  journal= {arXiv preprint arXiv:2208.14976},
  year   = {2022}
}
R2 v1 2026-06-28T00:30:14.201Z