Hi to all,
I am currently implementing Zou/He boundary conditions. Everything seems to work fine, but I notice mass increase over time, for example when imposing a Poisseuille profile on in- and outlet of a pipe. The rate depends on the given velocity, but is of the order of 0.1% per 1000 timesteps. It happens in my D3Q19 code as well as in the D2Q9 code presented on this page. I got around it by setting rho=1 instead of calculating the density from the incoming populations on the inlet, but this can’t be the best way to do it I assume?