A posteriori error estimates and stopping criteria for space-time domain decomposition for two-phase flow between different rock types

Elyes Ahmed (1, 2, 3, 4) , Sarah Ali Hassan (2, 3) , Caroline Japhet (1) , Michel Kern (2, 3) , Martin Vohralík (2, 3)

1 LAGA - Laboratoire Analyse, Géométrie et Applications
2 SERENA - Simulation for the Environment: Reliable and Efficient Numerical Algorithms
3 CERMICS - Centre d'Enseignement et de Recherche en Mathématiques et Calcul Scientifique
4 UiB - Department of Mathematics [Bergen]

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.

Mots clés: Linearization; Finite volume scheme; Two-phase Darcy flow; A posteriori error estimate; Stopping criteria; Optimized Schwarz waveform relaxation; Domain decomposition method; Discontinuous capillary pressure; Robin and Ventcell transmission conditions

Domaines: Analyse numérique [math.NA]

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. SMAI Journal of Computational Mathematics, 2019, 5, pp.195-227. ⟨10.5802/smai-jcm.47⟩. ⟨hal-01540956v3⟩