358	RWTH Publication No: 230329        2013        IGPM358.pdf
TITLE	A Multidomain Discretization of the Richards Equation in Layered Soil
AUTHORS	Heiko Berninger, Ralf Kornhuber, Oliver Sander
ABSTRACT	We consider the Richards equation on a domain that is decomposed into nonoverlapping layers, i.e., the decomposition has no cross points. We assume that the saturation and permeability functions are space-independent on each subdomain. Kirchhoff transformation of each subdomain problem separately then leads to a set of semi-linear equations, which can each be solved efficiently using monotone multigrid. The transformed subdomain problems are coupled by nonlinear continuity and flux conditions. This nonlinear coupled problem can be solved using substructuring methods like the Dirichlet-Neumann or Robin iteration. We give several numerical examples showing the discretization error, the solver robustness under variations of the soil parameters and a hydrological example with four soil layers and surface water.
KEYWORDS	Richards equation, domain decomposition, layered soil, nonlinear discretization, robustness, surface water
DOI	10.1007/s10596-014-9461-8
PUBLICATION	Computational Geosciences (2015) Volume 19, Issue 1, pp 213–232