Over the last decade, lattice Boltzmann methods have proven to be reliable and efficient tools for the numerical simulation of complex flows. The specifics of such methods as turbulence solvers, however, are not yet completely documented. This paper provides results of direct numerical simulations (DNS), by a lattice Boltzmann scheme, of fully developed, incompressible, pressure-driven turbulence between two parallel plates. These are validated against results from simulations using a standard Chebyshev pseudo-spectral method. Detailed comparisons, in terms of classical one-point turbulence statistics at moderate Reynolds number, with both numerical and experimental data show remarkable agreement.
Consequently, the choice of numerical method has, in sufficiently resolved DNS computations, no dominant effect at least on simple statistical quantities such as mean flow and Reynolds stresses. Since only the method-independent statistics can be credible, the choice of numerical method for DNS should be determined mainly through considerations of computational efficiency. The expected practical advantages of the lattice Boltzmann method, for instance against pseudo-spectral methods, are found to be significant even for the simple geometry and the moderate Reynolds number considered here. This permits the conclusion that the lattice Boltzmann approach is a promising DNS tool for incompressible turbulence.