A multidomain discretization of the Richards equation in layered soil

Research output: Contribution to journalResearch articleContributedpeer-review

Contributors

  • Heiko Berninger - , Free University of Berlin (Author)
  • Ralf Kornhuber - , Free University of Berlin (Author)
  • Oliver Sander - , RWTH Aachen University (Author)

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 semilinear 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.

Details

Original languageEnglish
Pages (from-to)213-232
Journal Computational geosciences : modeling, simulation and data analysis
Volume19
Publication statusPublished - 2014
Peer-reviewedYes
Externally publishedYes

External IDs

Scopus 84940001754
ORCID /0000-0003-1093-6374/work/142250577

Keywords