A posteriori error estimates and stopping criteria for space-time domain decomposition for two-phase flow between different rock types
The SMAI journal of computational mathematics, Volume 5 (2019) , pp. 195-227.

We consider two-phase flow in a porous medium composed of two different rock types, so that the capillary pressure field is discontinuous at the interface between the rocks. This is a nonlinear and degenerate parabolic problem with nonlinear and discontinuous transmission conditions on the interface. We first describe a space-time domain decomposition method based on the optimized Schwarz waveform relaxation algorithm (OSWR) with Robin or Ventcell transmission conditions. Complete numerical approximation is achieved by a finite volume scheme in space and the backward Euler scheme in time. We then derive a guaranteed and fully computable a posteriori error estimate that in particular takes into account the domain decomposition error. Precisely, at each iteration of the OSWR algorithm and at each linearization step, the estimate delivers a guaranteed upper bound on the error between the exact and the approximate solution. Furthermore, to make the algorithm efficient, the different error components given by the spatial discretization, the temporal discretization, the linearization, and the domain decomposition are distinguished. These ingredients are then used to design a stopping criterion for the OSWR algorithm as well as for the linearization iterations, which together lead to important computational savings. Numerical experiments illustrate the efficiency of our estimates and the performance of the OSWR algorithm with adaptive stopping criteria on a model problem in three space dimensions. Additionally, the results show how a posteriori error estimates can help determine the free Robin or Ventcell parameters.

Published online: 2019-12-12
DOI: https://doi.org/10.5802/smai-jcm.47
Classification: 65M08,  65M15,  65M50,  65M55,  76S05
Keywords: two-phase Darcy flow, discontinuous capillary pressure, finite volume scheme, domain decomposition method, optimized Schwarz waveform relaxation, Robin and Ventcell transmission conditions, linearization, a posteriori error estimate, stopping criteria
@article{SMAI-JCM_2019__5__195_0,
     author = {Elyes Ahmed and Sarah Ali Hassan and Caroline Japhet and Michel Kern and Martin Vohral\'\i k},
     title = {A posteriori error estimates and stopping criteria  for space-time domain decomposition  for two-phase flow between different rock types},
     journal = {The SMAI journal of computational mathematics},
     publisher = {Soci\'et\'e de Math\'ematiques Appliqu\'ees et Industrielles},
     volume = {5},
     year = {2019},
     pages = {195-227},
     doi = {10.5802/smai-jcm.47},
     language = {en},
     url = {smai-jcm.centre-mersenne.org/item/SMAI-JCM_2019__5__195_0/}
}
Elyes Ahmed; Sarah Ali Hassan; Caroline Japhet; Michel Kern; Martin Vohralík. A posteriori error estimates and stopping criteria  for space-time domain decomposition  for two-phase flow between different rock types. The SMAI journal of computational mathematics, Volume 5 (2019) , pp. 195-227. doi : 10.5802/smai-jcm.47. https://smai-jcm.centre-mersenne.org/item/SMAI-JCM_2019__5__195_0/

[1] Elyes Ahmed; Sarah Ali Hassan; Caroline Japhet; Michel Kern; Martin Vohralík A posteriori error estimates and stopping criteria for space-time domain decomposition for two-phase flow between different rock types (2017) (https://hal.inria.fr/hal-01540956)

[2] Elyes Ahmed; Jérôme Jaffré; Jean E. Roberts A reduced fracture model for two-phase flow with different rock types, Math. Comput. Simulation, Volume 137 (2017), pp. 49-70 | Article | MR 3624734

[3] Elyes Ahmed; Caroline Japhet; Michel Kern Global-in-time domain decomposition for a nonlinear diffusion problem (2019) (https://hal.inria.fr/hal-02263280, accepted in the Proceedings of the 25th International Conference on Domain Decomposition Methods)

[4] Elyes Ahmed; Caroline Japhet; Michel Kern Space-time domain decomposition for two-phase flow between different rock types (2019) (https://hal.inria.fr/hal-02275690)

[5] Clarisse Alboin; Jérôme Jaffré; Jean E. Roberts; Xuewen Wang; Christophe Serres Domain decomposition for some transmission problems in flow in porous media, Numerical treatment of multiphase flows in porous media (Beijing, 1999) (Lecture Notes in Physics) Volume 552, Springer, 2000, pp. 22-34 | Article | MR 1876007 | Zbl 1010.76050

[6] Sarah Ali Hassan; Caroline Japhet; Michel Kern; Martin Vohralík A posteriori stopping criteria for optimized Schwarz domain decomposition algorithms in mixed formulations, Comput. Methods Appl. Math., Volume 18 (2018) no. 3, pp. 495-519 | Article | MR 3824776 | Zbl 1398.65272

[7] Sarah Ali Hassan; Caroline Japhet; Martin Vohralík A posteriori stopping criteria for space-time domain decomposition for the heat equation in mixed formulations, Electron. Trans. Numer. Anal., Volume 49 (2018), pp. 151-181 | Article | MR 3819380 | Zbl 1397.65167

[8] Boris Andreianov; Konstantin Brenner; Clément Cancès Approximating the vanishing capillarity limit of two-phase flow in multi-dimensional heterogeneous porous medium, ZAMM, Z. Angew. Math. Mech., Volume 94 (2014) no. 7-8, pp. 655-667 | Article | MR 3230275 | Zbl 1297.76152

[9] Mario Arioli; Daniel Loghin; Andrew J. Wathen Stopping criteria for iterations in finite element methods, Numer. Math., Volume 99 (2005) no. 3, pp. 381-410 | Article | MR MR2117733 | Zbl 1069.65124

[10] Khalid Aziz; Antonín Settari Petroleum Reservoir Simulation, Applied Science Publishers, 1979

[11] Roland Becker; Claes Johnson; Rolf Rannacher Adaptive error control for multigrid finite element methods, Computing, Volume 55 (1995) no. 4, pp. 271-288 | Article | MR 1370103

[12] Daniel Bennequin; Martin J. Gander; Laurence Halpern A homographic best approximation problem with application to optimized Schwarz waveform relaxation, Math. Comput., Volume 78 (2009) no. 265, pp. 185-223 | Article | MR 2448703 | Zbl 1198.65179

[13] Heiko Berninger; Sébastien Loisel; Oliver Sander The 2-Lagrange multiplier method applied to nonlinear transmission problems for the Richards equation in heterogeneous soil with cross points, SIAM J. Sci. Comput., Volume 36 (2014) no. 5, p. A2166-A2198 | Article | MR 3261611 | Zbl 1355.74056

[14] Michiel Bertsch; Roberta Dal Passo; Cornelis J. Van Duijn Analysis of oil trapping in porous media flow, SIAM J. Math. Anal., Volume 35 (2003) no. 1, pp. 245-267 | Article | MR 2001474 | Zbl 1049.35108

[15] Konstantin Brenner; Clément Cancès; Danielle Hilhorst Finite volume approximation for an immiscible two-phase flow in porous media with discontinuous capillary pressure, Comput. Geosci., Volume 17 (2013) no. 3, pp. 573-597 | Article | MR 3050007 | Zbl 1392.76035

[16] Filipa Caetano; Martin J. Gander; Laurence Halpern; Jérémie Szeftel Schwarz waveform relaxation algorithms for semilinear reaction-diffusion equations, Netw. Heterog. Media, Volume 5 (2010) no. 3, pp. 487-505 | Article | MR 2670651 | Zbl 1262.35009

[17] Clément Cancès Nonlinear parabolic equations with spatial discontinuities, NoDEA, Nonlinear Differ. Equ. Appl., Volume 15 (2008) no. 4-5, pp. 427-456 | Article | MR 2465972 | Zbl 1180.35534

[18] Clément Cancès Finite volume scheme for two-phase flows in heterogeneous porous media involving capillary pressure discontinuities, M2AN Math. Model. Numer. Anal., Volume 43 (2009) no. 5, pp. 973-1001 | Article | Numdam | MR MR2559741 | Zbl 1171.76035

[19] Clément Cancès; Thierry Gallouët; Alessio Porretta Two-phase flows involving capillary barriers in heterogeneous porous media, Interfaces Free Bound., Volume 11 (2009) no. 2, pp. 239-258 | Article | MR MR2511641 | Zbl 1178.35196

[20] Clément Cancès; Iuliu Sorin Pop; Martin Vohralík An a posteriori error estimate for vertex-centered finite volume discretizations of immiscible incompressible two-phase flow, Math. Comput., Volume 83 (2014) no. 285, pp. 153-188 | Article | MR 3120585 | Zbl 06227551

[21] Guy Chavent; Jérôme Jaffré Mathematical models and finite elements for reservoir simulation, Studies in Mathematics and Its Applications, Volume 17, North-Holland, 1986 | Zbl 0603.76101

[22] Daniele A. Di Pietro; Eric Flauraud; Martin Vohralík; Soleiman Yousef A posteriori error estimates, stopping criteria, and adaptivity for multiphase compositional Darcy flows in porous media, J. Comput. Phys., Volume 276 (2014), pp. 163-187 | Article | MR 3252574 | Zbl 1349.76323

[23] Daniele A. Di Pietro; Martin Vohralík; Soleiman Yousef Adaptive regularization, linearization, and discretization and a posteriori error control for the two-phase Stefan problem, Math. Comput., Volume 84 (2015) no. 291, pp. 153-186 | Article | MR 3266956 | Zbl 1307.65123

[24] Guillaume Enchéry; R. Eymard; Anthony Michel Numerical approximation of a two-phase flow problem in a porous medium with discontinuous capillary forces, SIAM J. Numer. Anal., Volume 43 (2006) no. 6, pp. 2402-2422 | Article | MR 2206441 | Zbl 1145.76046

[25] Alexandre Ern; Martin Vohralík A posteriori error estimation based on potential and flux reconstruction for the heat equation, SIAM J. Numer. Anal., Volume 48 (2010) no. 1, pp. 198-223 | Article | MR 2608366 | Zbl 1215.65152

[26] Alexandre Ern; Martin Vohralík Adaptive inexact Newton methods with a posteriori stopping criteria for nonlinear diffusion PDEs, SIAM J. Sci. Comput., Volume 35 (2013) no. 4, p. A1761-A1791 | Article | MR 3072765 | Zbl 1362.65056

[27] Robert Eymard; Thierry Gallouët; Raphaèle Herbin Finite volume methods, Handbook of Numerical Analysis, Vol. VII, North-Holland, 2000, pp. 713-1020 | Article | MR MR1804748 | Zbl 0981.65095

[28] Robert Eymard; Thierry Gallouët; Raphaèle Herbin Finite volume approximation of elliptic problems and convergence of an approximate gradient, Appl. Numer. Math., Volume 37 (2001) no. 1-2, pp. 31-53 | Article | MR 1825115

[29] Martin J. Gander Optimized Schwarz methods, SIAM J. Numer. Anal., Volume 44 (2006) no. 2, pp. 699-731 | Article | MR 2218966 | Zbl 1117.65165

[30] Martin J. Gander; Laurence Halpern; Frédéric Nataf Optimal Schwarz waveform relaxation for the one dimensional wave equation, SIAM J. Numer. Anal., Volume 41 (2003) no. 5, pp. 1643-1681 | Article | MR 2035001 | Zbl 1085.65077

[31] Benjamin Ganis; Kundan Kumar; Gergina Pencheva; Mary F. Wheeler; Ivan Yotov A global Jacobian method for mortar discretizations of a fully implicit two-phase flow model, Multiscale Model. Simul., Volume 12 (2014) no. 4, pp. 1401-1423 | Article | MR 3267144 | Zbl 1312.76027

[32] Florian Haeberlein; Laurence Halpern; Anthony Michel Newton-Schwarz optimised waveform relaxation Krylov accelerators for nonlinear reactive transport, Domain decomposition methods in science and engineering XX (Lect. Notes Comput. Sci. Eng.) Volume 91, Springer, 2013, pp. 387-394 | Article | MR 3243014

[33] Laurence Halpern; Florence Hubert A finite volume Ventcell-Schwarz algorithm for advection-diffusion equations, SIAM J. Numer. Anal., Volume 52 (2014) no. 3, pp. 1269-1291 | Article | MR 3213798 | Zbl 1382.65360

[34] Rainer Helmig Multiphase Flow and Transport Processes in the Subsurface, Springer, 1997

[35] Thi-Thao-Phuong Hoang; Jérôme Jaffré; Caroline Japhet; Michel Kern; Jean E. Roberts Space-time domain decomposition methods for diffusion problems in mixed formulations, SIAM J. Numer. Anal., Volume 51 (2013) no. 6, pp. 3532-3559 | Article | MR 3144798 | Zbl 1295.65095

[36] Thi-Thao-Phuong Hoang; Caroline Japhet; Michel Kern; Jean E. Roberts Ventcell conditions with mixed formulations for flow in porous media, Decomposition Methods in Science and Engineering XXII (Lecture Notes in Computational Science and Engineering) Volume 104 (2016), pp. 531-540 | Article | MR 3643584

[37] Caroline Japhet; Frédéric Nataf The best interface conditions for domain decomposition methods: absorbing boundary conditions, Absorbing Boundaries and Layers, Domain Decomposition Methods, Nova Sci. Publ., 2001, pp. 348-373 | MR 2039948

[38] Pavel Jiránek; Zdeněk Strakoš; Martin Vohralík A posteriori error estimates including algebraic error and stopping criteria for iterative solvers, SIAM J. Sci. Comput., Volume 32 (2010) no. 3, pp. 1567-1590 | Article | MR 2652091 | Zbl 1215.65168

[39] Hanyu Li; Mary F. Wheeler Sequential Refinement Solver using Space-Time Domain Decomposition for Non-linear Multiphase Flow Problems (2019) (https://arxiv.org/abs/1901.09436)

[40] Knut-Andreas Lie; Stein Krogstad; Ingeborg S. Ligaarden; Jostein R. Natvig; Halvor M. Nilsen; Bård Skaflestad Open source MATLAB implementation of consistent discretisations on complex grids, Comput. Geosci., Volume 16 (2012) no. 2, pp. 297-322 | Article | Zbl 1348.86002

[41] Véronique Martin An optimized Schwarz waveform relaxation method for the unsteady convection diffusion equation in two dimensions, Appl. Numer. Math., Volume 52 (2005) no. 4, pp. 401-428 | Article | MR 2112867 | Zbl 1070.65088

[42] Ricardo H. Nochetto; Alfred Schmidt; Claudio Verdi A posteriori error estimation and adaptivity for degenerate parabolic problems, Math. Comput., Volume 69 (2000) no. 229, pp. 1-24 | Article | Numdam | MR MR1648399 | Zbl 0942.65111

[43] Jan Papež; Ulrich Rüde; Martin Vohralík; Barbara Wohlmuth Sharp algebraic and total a posteriori error bounds for h and p finite elements via a multilevel approach (2017) (https://hal.inria.fr/hal-01662944, submitted for publication)

[44] Jan Papež; Zdeněk Strakoš; Martin Vohralík Estimating and localizing the algebraic and total numerical errors using flux reconstructions, Numer. Math., Volume 138 (2018) no. 3, pp. 681-721 | Article | MR 3767698 | Zbl 06857783

[45] Gergina Pencheva; Martin Vohralík; Mary F. Wheeler; Tim Wildey Robust a posteriori error control and adaptivity for multiscale, multinumerics, and mortar coupling, SIAM J. Numer. Anal., Volume 51 (2013) no. 1, pp. 526-554 | Article | MR 3033022 | Zbl 1267.65165

[46] Valentine Rey; Pierre Gosselet; Christian Rey Strict lower bounds with separation of sources of error in non-overlapping domain decomposition methods, Internat. J. Numer. Methods Engrg., Volume 108 (2016) no. 9, pp. 1007-1029 | Article | MR 3573714 | Zbl 1296.65167

[47] Valentine Rey; Christian Rey; Pierre Gosselet A strict error bound with separated contributions of the discretization and of the iterative solver in non-overlapping domain decomposition methods, Comput. Methods Appl. Mech. Engrg., Volume 270 (2014), pp. 293-303 | Article | MR 3154030 | Zbl 1296.65167

[48] David Seus; Koondanibha Mitra; Iuliu Sorin Pop; Florin Adrian Radu; Christian Rohde A linear domain decomposition method for partially saturated flow in porous media, Comput. Methods Appl. Mech. Engrg., Volume 333 (2018), pp. 331-355 | Article | MR 3771899

[49] Gurpreet Singh; Mary F. Wheeler A space-time domain decomposition approach using enhanced velocity mixed finite element method, J. Comput. Phys., Volume 374 (2018), pp. 893-911 | Article | MR 3864481 | Zbl 1416.65333

[50] Jan Ole Skogestad; Eirik Keilegavlen; Jan M. Nordbotten Domain decomposition strategies for nonlinear flow problems in porous media, J. Comput. Phys., Volume 234 (2013), pp. 439-451 | Article | MR 2999786 | Zbl 1286.65179

[51] Jan Ole Skogestad; Eirik Keilegavlen; Jan M. Nordbotten Two-scale preconditioning for two-phase nonlinear flows in porous media, Transp. Porous Media, Volume 114 (2016) no. 2, pp. 485-503 | Article | MR 3540014 | Zbl 1327.76138

[52] Cornelis J. van Duijn; Johannes Molenaar; M. J. de Neef The effect of capillary forces on immiscible two-phase flow in heterogeneous porous media, Transport in Porous Media, Volume 21 (1995) no. 1, pp. 71-93 | Article

[53] Martinus Th. van Genuchten A closed form for predicting the hydraulic conductivity of unsaturated soils, Soil Sci. Soc. Amer. J., Volume 44 (1980), pp. 892-898 | Article

[54] Martin Vohralík; Mary F. Wheeler A posteriori error estimates, stopping criteria, and adaptivity for two-phase flows, Comput. Geosci., Volume 17 (2013) no. 5, pp. 789-812 | Article | MR 3104635 | Zbl 1393.76069

[55] Ivan Yotov A mixed finite element discretization on non-matching multiblock grids for a degenerate parabolic equation arising in porous media flow, East-West J. Numer. Math., Volume 5 (1997) no. 3, pp. 211-230 | MR 1485115 | Zbl 0897.76057

[56] Ivan Yotov Interface solvers and preconditioners of domain decomposition type for multiphase flow in multiblock porous media, Scientific computing and applications (Adv. Comput. Theory Pract.) Volume 7, Nova Sci. Publ., Huntington, NY, 2001, pp. 157-167 | MR 2032407 | Zbl 1153.76432