  This article is the first of two in which we develop a relaxation finite
volume scheme for the convective part of the multiphase flow models introduced
in the series of papers [10, 9, 4]. In the present article we focus on
barotropic flows where in each phase the pressure is a given function of the
density. The case of general equations of state will be the purpose of the
second article. We show how it is possible to extend the relaxation scheme
designed in [8] for the barotropic Baer-Nunziato two-phase flow model to the
multiphase flow model with N-where N is arbitrarily large-phases. The obtained
scheme inherits the main properties of the relaxation scheme designed for the
Baer-Nunziato two phase flow model. The approximated phase fractions and phase
densities are proven to remain positive and a discrete energy inequality is
also proven under a classical CFL condition. For the same level of refinement,
the relaxation scheme is shown to be much more accurate than Rusanov's scheme,
and for a given level of approximation error, the relaxation scheme is shown to
perform much better in terms of computational cost than Rusanov's scheme.
Moreover, contrary to Rusanov's scheme which develops strong oscillations when
approximating vanishing phase solutions, the numerical results show that the
relaxation scheme remains stable in such regimes.
