Wrong momentum density calculation for compressible methods with forces.
In conservedquantitycomputation.py (line 210f) the calculation of the momentum density uses the macroscopic velocity shift
mom_density_eq_coll = apply_force_model_shift('macroscopic_velocity_shift', dim, mom_density_eq_coll, self._forceModel, self._compressible)
which defaults to . When using this momentum density to calculcate the velocity (as it is done in all generated lattice models in waLBerla), we divide by the density again, resulting in the overall velocity equation , which is obviously wrong.
Affected are all lbmpy simulations that work directly on the momentum denisty, and all waLBerla simulations that use generated lattice models and calculate the velocity at some point (in particular, the VTK output will be slightly off).