Fulltext

A multidomain discretization of the Richards equation in layered soil

Berninger, Heiko ; Kornhuber, Ralf ; Sander, Oliver

In: Computational Geosciences, 2015, vol. 19, no. 1, p. 213-232

Ajouter à la liste personnelle
    Titre
    • A multidomain discretization of the Richards equation in layered soil
    Auteur
    • Berninger, Heiko. Section de Mathématiques, Université de Genève, 2-4 rue du Lièvre, CP 64, 1211, Genève 4, Switzerland
    • Kornhuber, Ralf. Freie Universität Berlin, Institut für Mathematik, Arnimallee 6, 14195, Berlin, Germany
    • Sander, Oliver. Institut für Geometrie und Praktische Mathematik, RWTH Aachen, Templergraben 55, 52062, Aachen, Germany
    Type de document
    • Postprint
    Langue
    • Anglais
    Publié dans
    • Computational Geosciences, 2015, vol. 19, no. 1, p. 213-232. Springer International Publishing
    Autre version électronique
    Identifiant OAI-PMH
    • oai:doc.rero.ch:332750
    Summary
    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.