when simulating multiphase flow in 3d pipes I get perfect volume and mass balance
when I add obstacles or I simulate flow in porous media the balance between the injected and the discharged phases does not hold
his this a general problem with lattice Boltzmann methods due to compressibility ?
A programming error in my code is possible but unlikely given all the specific check I have done
I have a given velocity profile ad the inlet and I set rho=1 at the discharge using Zou boundary conditions
thank you for your attention, giovanna