Dear All.

It seems from the book of Dieter A. Wolf-Gladrow:

“Lattice-Gas Cellular Automata and Lattice Boltzmann Models - An Introduction”, pag 173-174, that the collision operator and the body forces don’t contribute to the mass conservation equation. By running the matlab code “Rayleigh-Bénard convection” i tried to consider their contribution. By summing the 9 terms for the distribution equation:

fOut(i,:, = fIn(i,:, - omegaNS .* (fIn(i,:,:)-fEq(i,:,:)) + force(i,:,:);

it seems that, while the body force sum is zero, the streaming and the collision operator contribution have the same value. I get confused about that so maybe someone may help me. Thanks,

Stefano