Implements pressure and stress correction for multiphase model with density contrast.
Implements pressure and stress correction for a multiphase model with density contrast. Pressure is regularized and hence different from the form used by Safari and Fakahri.