A new approach to lattice Boltzmann simulation is introduced which avoids the direct numerical determination of the single-particle distribution from the Boltzmann collision operator. Instead, one starts from the Chapman-Enskog expansion of the single