Accuracy of Discrete-Velocity BGK Models for the Simulation of the Incompressible Navier-Stokes Equations

Abstract

Two discretizations of a 9-velocity Boltzmann equation with a BGK collision operator are studied. A Chapman-Enskog expansion of the PDE system predicts that the macroscopic behavior corresponds to the incompressible Navier-Stokes equations with additional terms of order Mach number squared. We introduce a fourth-order scheme and compare results with those of the commonly used lattice Boltzmann discretization and with finite-difference schemes applied to the incompressible Navier-Stokes equations in primitive-variable form. We numerically demonstrate convergence of the BGK schemes to the incompressible Navier-Stokes equations and quantify the errors associated with compressibility and discretization effects. When compressibility error is smaller than discretization error, convergence in both grid spacing and time step is shown to be second-order for the LB method and is confirmed to be fourth-order for the fourth-order BGK solver. However, when the compressibility error is simultaneously reduced as the grid is refined, the LB method behaves as a first-order scheme in time.