  Cascaded lattice-Boltzmann method (Cascaded-LBM) employs a new class of
collision operators aiming to improve numerical stability. It achieves this and
distinguishes from other collision operators, such as in the standard single or
multiple relaxation time approaches, by performing relaxation process due to
collisions in terms of moments shifted by the local hydrodynamic fluid
velocity, i.e. central moments, in an ascending order-by-order at different
relaxation rates. In this paper, we propose and derive source terms in the
Cascaded-LBM to represent the effect of external or internal forces on the
dynamics of fluid motion. This is essentially achieved by matching the
continuous form of the central moments of the source or forcing terms with its
discrete version. Different forms of continuous central moments of sources,
including one that is obtained from a local Maxwellian, are considered in this
regard. As a result, the forcing terms obtained in this new formulation are
Galilean invariant by construction. The method of central moments along with
the associated orthogonal properties of the moment basis completely determines
the expressions for the source terms as a function of the force and macroscopic
velocity fields. In contrast to the existing forcing schemes, it is found that
they involve higher order terms in velocity space. It is shown that the
proposed approach implies "generalization" of both local equilibrium and source
terms in the usual lattice frame of reference, which depend on the ratio of the
relaxation times of moments of different orders. An analysis by means of the
Chapman-Enskog multiscale expansion shows that the Cascaded-LBM with forcing
terms is consistent with the Navier-Stokes equations. Computational experiments
with canonical problems involving different types of forces demonstrate its
accuracy.
