That sounds good – a common way to deal with the mass is to project the two avergage densities in the considered domain to a certain values. It is like the pressure in the Navier-Stokes equations which is definded up to a constant. So you basically fix that constant. There is a technical report on how it works for single component flows on our website. Best Mathias