  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.
