|
|||||
|
|
Solving the nonstationary Richards equation with adaptive hp-FEM |
|
| Authors: | Pavel Solin Michal Kuraz |
| Institution: | a Department of Mathematics and Statistics, University of Nevada, Reno, USA b Institute of Themomechanics, Prague, Czech Republic c Department of Water Resources and Environmental Modeling, Czech University of Life Sciences Prague, Czech Republic |
| Abstract: | This paper examines the potential of the adaptive hp-FEM method for the numerical solution of time-dependent variably saturated Darcian flow problems described by the Richards equation. The method is illustrated on three model problems: a benchmark with known exact solution, groundwater seepage into a dry lysimeter box with time-dependent boundary conditions, and capillary barrier behavior under an intense infiltration. In the second part of the paper we present the weak formulation of the Richards equation for the Newton’s and Picard’s methods, give a brief overview of adaptive hp-FEM with emphasis on aspects that are usually not discussed in the literature, and we briefly introduce the open source adaptive hp-FEM library HERMES that was used to generate numerical results for this paper. All computations that we present are easily reproducible. |
| Keywords: | Richards equation Adaptive hp-FEM Dynamically-changing meshes Reproducible research |
| 本文献已被 ScienceDirect 等数据库收录！ | |