An adaptive strategy for discontinuous Galerkin simulations of Richards' equation

Jean-Baptiste Clément (1, 2) , Frederic Golay (1) , Mehmet Ersoy (1) , Damien Sous (2, 3)

1 IMATH - Institut de Mathématiques de Toulon - EA 2134
2 MIO - Institut méditerranéen d'océanologie
3 SIAME - Laboratoire des Sciences de l'Ingénieur Appliquées à la Mécanique et au génie Electrique

Numerical solution of Richards' equation remains challenging to get robust, accurate and cost-effective results, particularly for moving sharp wetting fronts. An adaptive strategy for both space and time is proposed to deal with 2D sharp wetting fronts associated with varying and possibly vanishing diffusivity caused by nonlinearity, heterogeneity and anisotropy. Adaptive time stepping makes nonlinear convergence reliable and backward difference formula provides high-order time scheme. Adaptive mesh refinement tracks wetting fronts with an a posteriori error indicator. The novelty of this paper consists in using this technique in combination with a weighted discontinuous Galerkin framework to better approximate steep wetting fronts by a discontinuity. The potential of the overall approach is shown through various examples including analytical and laboratory benchmarks and simulation of full-scale multi-materials dam wetting experiment.

Mots clés: Unsaturated porous media; Weighted Discontinuous Galerkin method; Adaptive Mesh Refinement; A posteriori error estimation; Adaptive time stepping; BDF

Domaines: Mécanique des fluides [physics.class-ph]; Equations aux dérivées partielles [math.AP]

Jean-Baptiste Clément, Frederic Golay, Mehmet Ersoy, Damien Sous. An adaptive strategy for discontinuous Galerkin simulations of Richards' equation. 2021. ⟨hal-03132113v1⟩