A linear domain decomposition method for partially saturated flow in porous media

Abstract

The Richards equation is a nonlinear parabolic equation that is commonly used for modelling saturated/unsaturated flow in porous media. We assume that the medium occupies a bounded Lipschitz domain partitioned into two disjoint subdomains separated by a fixed interface $\Gamma$. This leads to two problems defined on the subdomains which are coupled through conditions expressing flux and pressure continuity at $\Gamma$. After an Euler implicit discretisation of the resulting nonlinear subproblems a linear iterative ($L$-type) domain decomposition scheme is proposed. The convergence of the scheme is proved rigorously. In the last part we present numerical results that are in line with the theoretical finding, in particular the unconditional convergence of the scheme. We further compare the scheme to other approaches not making use of a domain decomposition. Namely, we compare to a Newton and a Picard scheme. We show that the proposed scheme is more stable than the Newton scheme while remaining comparable in computational time, even if no parallelisation is being adopted. Finally we present a parametric study that can be used to optimize the proposed scheme.

34 pages, 13 figures, 7 tables

Cited by in corpus (6)

• Anderson accelerated fixed-stress splitting schemes for consolidation of unsaturated porous media
• Upscaling of unsaturated flow in fractured porous media
• Two-phase hyperbolic model for porous media saturated with a viscous fluid and its application to wavefields simulation
• The Multiscale Perturbation Method for Two-Phase Reservoir Flow Problems
• Splitting-based domain decomposition methods for two-phase flow with different rock types
• Monolithic and splitting based solution schemes for fully coupled quasi-static thermo-poroelasticity with nonlinear convective transport