Abstract:Richards equation is often used to represent two-phase fluid flow in an unsaturated porous medium when one phase is much heavier and more viscous than the other. However, it cannot describe the fully saturated flow for some capillary functions without specialized treatment due to degeneracy in the capillary pressure term. Mathematically, gravity-dominated variably saturated flows are interesting because their governing partial differential equation switches from hyperbolic in the unsaturated region to elliptic in the saturated region. Moreover, the presence of wetting fronts introduces strong spatial gradients often leading to numerical instability. In this work, we develop a robust, multidimensional mathematical model and implement a well-known efficient and conservative numerical method for such variably saturated flow in the limit of negligible capillary forces. The elliptic problem in saturated regions is integrated efficiently into our framework by solving a reduced system corresponding only to the saturated cells using fixed head boundary conditions in the unsaturated cells. In summary, this coupled hyperbolic-elliptic PDE framework provides an efficient, physics-based extension of the hyperbolic Richards equation to simulate fully saturated regions. Finally, we provide a suite of easy-to-implement yet challenging benchmark test problems involving saturated flows in one and two dimensions. These simple problems, accompanied by their corresponding analytical solutions, can prove to be pivotal for the code verification, model validation (V&V) and performance comparison of simulators for variably saturated flow. Our numerical solutions show an excellent comparison with the analytical results for the proposed problems. The last test problem on two-dimensional infiltration in a stratified, heterogeneous soil shows the formation and evolution of multiple disconnected saturated regions.

Modeling an Aquifer: Numerical Solution to the Groundwater Flow Equation

A Mathematical Model and its Semi-analytical Solution for the Flow of Anisotropic Confined Aquifer with Two Radial Layers Considering Heterogeneity of Well Skin

Coupled Finite-Volume Model for 2D Surface and 3D Subsurface Flows

A Physically-Based Integrated Numerical Model for Flow, Upland Erosion, and Contaminant Transport in Surface-Subsurface Systems

Modelling of shallow aquifers in interaction with overland water

Analytical modeling of 2D groundwater flow in a semi-infinite heterogeneous domain with variable lateral sources

2d Well-Balanced Positivity Preserving Shallow Flow Model With Wetting And Drying On Triangular Grids

Mathematical Analysis of a Subsurface Flow Model

Numerical Approach of a Coupled Pressure-Saturation Model Describing Oil-Water Flow in Porous Media

A well-balanced positivity preserving two-dimensional shallow flow model with wetting and drying fronts over irregular topography

Comprehensive two-dimensional analytical modeling of groundwater levels in bi-directional sloping heterogeneous aquifers under variable recharge conditions

Modeling of groundwater flow and transport in coastal karst aquifers

Analytical 1-D dual-porosity equivalent solutions to 3-D discrete single-continuum models. Application to karstic spring hydrograph modelling

A hyperbolic-elliptic PDE model and conservative numerical method for gravity-dominated variably-saturated groundwater flow

An analytical and meshless model for 3D transient flow in a confined aquifer with nonuniform thickness: Application to stream depletion due to groundwater extraction

Méthode Des Volumes Finis Et Précision Des Modèles Numériques Des Écoulements Souterrains

Mathematical Development and Verification of a Non-Orthogonal Finite Volume Model for Groundwater Flow Applications

Modeling of unsaturated flow through porous media using meshless methods

A three-dimensional numerical model for variably saturated groundwater flow using meshless weak-strong form method

Geodesic Distance Integration in Analytical Frameworks for Aquifer Hydraulic Modeling