Modeling Groundwater Flow in Unconfined Conditions: Numerical Model and Solvers’ Efficiency
- Authors: Anuprienko D.V.1, Kapyrin I.V.1,2
-
Affiliations:
- Nuclear Safety Institute of the Russian Academy of Sciences
- Marchuk Institute of Numerical Mathematics of the Russian Academy of Sciences
- Issue: Vol 39, No 7 (2018)
- Pages: 867-873
- Section: Article
- URL: https://journals.rcsi.science/1995-0802/article/view/202598
- DOI: https://doi.org/10.1134/S1995080218070053
- ID: 202598
Cite item
Abstract
A mathematical model for variably saturated flow in unconfined conditions is presented. The model is based on pseudo-unsaturated approach using Richards equation with piecewise linear dependencies between hydraulic head, water content and relative permeability. It is implemented in GeRa (Geomigration of Radionuclides) software package, which is designed for modeling groundwater flow and contaminant transport in porous media and uses finite volume methods on unstructured grids. We consider two nonlinear solvers for nonlinear equations arising from discretization of the Richards equation, namely Newton and Picard methods. A special method for correction of hydraulic head values within the iterations of nonlinear solvers is proposed. The developed numerical techniques are applied to two test cases: dam seepage and real-world groundwater flow problems.
About the authors
D. V. Anuprienko
Nuclear Safety Institute of the Russian Academy of Sciences
Author for correspondence.
Email: denis-anuprienko@yandex.ru
Russian Federation, ul. Bol’shaya Tul’skaya 52, Moscow, 115191
I. V. Kapyrin
Nuclear Safety Institute of the Russian Academy of Sciences; Marchuk Institute of Numerical Mathematics of the Russian Academy of Sciences
Email: denis-anuprienko@yandex.ru
Russian Federation, ul. Bol’shaya Tul’skaya 52, Moscow, 115191; ul. Gubkina 8, Moscow, 119991