Dear Randy,rnrnSorry for the late reply. Dividing by rho0 is just doing scaling. You need to do that also for the results you get in OpenLB for which you can use the class UnitConverter. The term -1 comes from the followed implemetation strategy in OpenLB (cf. answer question 1).rnrnMathias