@dpiponi Hmm, looking at the code, his rho_0 for that k * (rho - rho_0) is the average fluid density in the sim grid. So it's going to have similar non-physical effects that you get trying to enforce global conservation of energy, mass, etc by reinjecting it evenly. But probably okay for small closed systems or if you tune it by hand.