The paper is devoted to the computation of two-phase flows in a porous mediumwhen applying the two-fluid approach. The basic formulation is presented first, together with the main properties of the model. A few basic analytic solutions are then provided, some of them correspondingto solutions of the one-dimensional Riemann problem. Three distinct Finite-Volume schemes are then introduced. The first two schemes, which rely on the Rusanov scheme,are shown to give wrong approximations in some cases involving sharp porous profiles.The third one, which is an extension of a scheme proposed by Kröner and Thanh [SIAM J. Numer. Anal. 43 (2006) 796–824]for the computation of single phase flows in varying cross section ducts,provides fair results in all situations. Properties of schemes and numerical results are presented. Analytic tests enable to compute the L 1 norm of the error.

We consider the system of partial differential equations governingthe one-dimensional flow of two superposed immiscible layers of shallow water. The difficulty in this system comes from the coupling terms involving some derivatives of the unknowns that make the system nonconservative, and eventually nonhyperbolic. Due to these terms, a numerical scheme obtained by performing an arbitrary scheme to each layer, and using time-splitting or other similar techniques leads to instabilities in general. Here we use entropy inequalities in order to control the stability. We introduce a stable well-balanced time-splitting scheme for the two-layer shallow water system that satisfies a fully discrete entropy inequality. In contrast with Roe type solvers, it does not need the computation of eigenvalues, which is not simple for the two-layer shallow water system. The solver has the property to keep the water heights nonnegative, and to be able to treat vanishing values.