A Linear Domain Decomposition Method for Non-equilibrium Two-Phase Flow Models

Abstract

We consider a model for two-phase flow in a porous medium posed in a domain consisting of two adjacent regions. The model includes dynamic capillarity and hysteresis. At the interface between adjacent subdomains, the continuity of the normal fluxes and pressures is assumed. For finding the semi-discrete solutions after temporal discretization by the θ-scheme, we proposed an iterative scheme. It combines a (fixed-point) linearization scheme and a non-overlapping domain decomposition method. This article describes the scheme, its convergence and a numerical study confirming this result. The convergence of the iteration towards the solution of the semi-discrete equations is proved independently of the initial guesses and of the spatial discretization, and under some mild constraints on the time step. Hence, this scheme is robust and can be easily implemented for realistic applications.