In this paper, we consider a class of models describing multiphase fluids in the framework of mixture theory. The considered systems, in their more general form, contain both the gradient of a hydrostatic pressure, generated by an incompressibility constraint, and a compressible pressure depending on the volume fractions of some of the different phases. To approach these systems, we propose an approximation based on the Leray projection, which involves the use of a symbolic symmetrizer for quasi-linear hyperbolic systems and related paradifferential techniques. In two space dimensions, we prove the well-posedness of this approximation and its convergence to the unique classical solution to the original system. In the last part, we shortly discuss the three dimensional case.