  A three-dimensional (3D) lattice Boltzmann method based on central moments is
derived. Two main elements are the local attractors in the collision term and
the source terms representing the effect of external and/or self-consistent
internal forces. For suitable choices of the orthogonal moment basis for the
three-dimensional, twenty seven velocity (D3Q27), and, its subset, fifteen
velocity (D3Q15) lattice models, attractors are expressed in terms of
factorization of lower order moments as suggested in an earlier work; the
corresponding source terms are specified to correctly influence lower order
hydrodynamic fields, while avoiding aliasing effects for higher order moments.
These are achieved by successively matching the corresponding continuous and
discrete central moments at various orders, with the final expressions written
in terms of raw moments via a transformation based on the binomial theorem.
Furthermore, to alleviate the discrete effects with the source terms, they are
treated to be temporally semi-implicit and second-order, with the implicitness
subsequently removed by means of a transformation. As a result, the approach is
frame-invariant by construction and its emergent dynamics describing fully 3D
fluid motion in the presence of force fields is Galilean invariant. Numerical
experiments for a set of benchmark problems demonstrate its accuracy.
