We present a systematic derivation of a model based on the central moment lattice Boltzmann equation that rigorously maintains Galilean invariance of forces to simulate inertial frame independent flow fields. In this regard, the central moments, i.e. moments shifted by the local fluid velocity, of the discrete source terms of the lattice Boltzmann equation are obtained by matching those of the continuous full Boltzmann equation of various orders. This results in an exact hierarchical identity between the central moments of the source terms of a given order and the components of the central moments of the distribution functions and sources of lower orders. The corresponding source terms in velocity space are then obtained from an exact inverse transformation due to a suitable choice of orthogonal basis for moments. Furthermore, such a central moment based kinetic model is further extended by incorporating reduced compressibility effects to represent incompressible flow. Moreover, the description and simulation of fluid turbulence for full or any subset of scales or their averaged behavior should remain independent of any inertial frame of reference. Thus, based on the above formulation, a new approach in lattice Boltzmann framework to incorporate turbulence models for simulation of Galilean invariant statistical averaged or filtered turbulent fluid motion is discussed.

We analyze the Lattice Boltzmann method for the simulation of fluctuating hydrodynamics by Adhikari et al. [Europhys. Lett., 71 (2005), 473-479] and find that it shows excellent agreement with theory even for small wavelengths as long as a stationary system is considered. This is in contrast to other finite difference and older lattice Boltzmann implementations that show convergence only in the limit of large wavelengths. In particular cross correlators vanish to less than 0.5%. For larger mean velocities, however, Galilean invariance violations manifest themselves.