In the mixed formulation developed in [20], both the pressure p and the relative density ϱ are kept as two independent variables and the complementarity conditions (10) are enforced using a non-smooth constraint function. Actually, for convenience, the relative density ϱ is replaced by the void fraction λ ¼ 1 ϱ, for which the complementarity conditions take the standard form