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 Γ. This leads to two problems defined on the subdomains which are coupled through conditions expressing flux and pressure continuity at Γ. 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 convergence of the scheme under mild restrictions on the time step size. 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. After presenting a parametric study that can be used to optimise the proposed scheme, we briefly discuss the effect of parallelisation and give an example of a four-domain implementation.
%0 Journal Article
%1 seus2018linear
%A Seus, David
%A Mitra, Koondanibha
%A Pop, Iuliu Sorin
%A Radu, Florin Adrian
%A Rohde, Christian
%D 2018
%I Elsevier
%J Computer methods in applied mechanics and engineering
%K ubs_10008 ubs_20013 ubs_30123 ubs_40190 unibibliografie
%P 331-355
%R 10.1016/j.cma.2018.01.029
%T A linear domain decomposition method for partially saturated flow in porous media
%V 333
%X 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 Γ. This leads to two problems defined on the subdomains which are coupled through conditions expressing flux and pressure continuity at Γ. 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 convergence of the scheme under mild restrictions on the time step size. 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. After presenting a parametric study that can be used to optimise the proposed scheme, we briefly discuss the effect of parallelisation and give an example of a four-domain implementation.
@article{seus2018linear,
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 Γ. This leads to two problems defined on the subdomains which are coupled through conditions expressing flux and pressure continuity at Γ. 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 convergence of the scheme under mild restrictions on the time step size. 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. After presenting a parametric study that can be used to optimise the proposed scheme, we briefly discuss the effect of parallelisation and give an example of a four-domain implementation.},
added-at = {2020-03-27T14:54:59.000+0100},
author = {Seus, David and Mitra, Koondanibha and Pop, Iuliu Sorin and Radu, Florin Adrian and Rohde, Christian},
biburl = {https://puma.ub.uni-stuttgart.de/bibtex/2eef8eeacad284b9fb6a57b8729ad19eb/unibiblio},
doi = {10.1016/j.cma.2018.01.029},
interhash = {d47c491f240e9476bc8745d110eec997},
intrahash = {eef8eeacad284b9fb6a57b8729ad19eb},
issn = {0045-7825},
journal = {Computer methods in applied mechanics and engineering},
keywords = {ubs_10008 ubs_20013 ubs_30123 ubs_40190 unibibliografie},
language = {eng},
pages = {331-355},
publisher = {Elsevier},
timestamp = {2020-03-27T13:54:59.000+0100},
title = {A linear domain decomposition method for partially saturated flow in porous media},
volume = 333,
year = 2018
}