Published online by Cambridge University Press: 31 October 2018
The Wasserstein gradient flow structure of the partial differential equation system governing multiphase flows in porous media was recently highlighted in Cancès et al. [Anal. PDE10(8), 1845–1876]. The model can thus be approximated by means of the minimising movement (or JKO after Jordan, Kinderlehrer and Otto [SIAM J. Math. Anal.29(1), 1–17]) scheme that we solve thanks to the ALG2-JKO scheme proposed in Benamou et al. [ESAIM Proc. Surv.57, 1–17]. The numerical results are compared to a classical upstream mobility finite volume scheme, for which strong stability properties can be established.
C. C. was supported by the French National Research Agency (ANR) through grant ANR-13-JS01-0007-01 (project GEOPOR) and ANR-11-LABX0007-01 (Labex CEMPI). L. M. was partially supported by the Portuguese Science Foundation through FCT grant PTDC/MAT-STA/0975/2014. T. O. G. was partially supported by the Fonds de la Recherche Scientifique – FNRS under grant MIS F.4539.16.