Mathematical model of groundwater head variation processes in heterogeneous porous media
Received: 2026-10-14
Published: 2026-09-30
Abstract
ter head in heterogeneous porous media, and illustrates them with a transient multi-well simulation. The relevance of the study stems from
the need for efficient management of groundwater resources, for ecological sustainability, and for forecasting hydrological processes under
a changing climate. Starting from Darcy’s law and the law of mass conservation, a second-order parabolic partial differential equation with
variable coefficients is derived that describes filtration in anisotropic heterogeneous media; anisotropy is represented through a diagonal
hydraulic-conductivity tensor and heterogeneity through spatially variable parameters. Boundary conditions are formulated for all six faces
of the domain, including river–aquifer interaction, infiltration, evapotranspiration, and irrigation return flow. The differential problem is
discretized by an explicit finite-difference scheme in which the interface conductivities in the heterogeneous medium are evaluated by the
harmonic mean, which correctly reproduces series-connected hydraulic resistances. The stability of the scheme is analysed by the von
Neumann method and a Courant–Friedrichs–Lewy stability condition is obtained; convergence follows from the Lax equivalence theorem.
A transient numerical experiment is carried out for a 2000 × 2000 m aquifer containing two extraction wells and one injection well: the
computations reproduce the progressive formation and interaction of the depression cones and the recharge mound, the transient head
response at the wells, and the three-dimensional structure of the head surface. The computational algorithm parallelizes naturally and can
be applied to groundwater-resource management, to forecasting the hydrodynamic state of aquifers, and to the design of wells and drainage systems.
Keywords
List of references
-
Lu C, Xu Z, Yang J, et al. An improved numerical model for groundwater flow simulation with MPFA method on arbitrary polygon grids. J Hydrol. 2022;606:127425. https://doi.org/10.1016/j.jhydrol.2021.127425.
-
Mohtashami A, Akbarpour A, Mollazadeh M. Development of a hybrid meshless local Petrov–Galerkin and particle filter model for unconfined aquifer simulation. Groundwater. 2022;60(3):375–388. https://doi.org/10.1111/gwat.13161.
-
Swathi B, Eldho TI. Generalized finite difference method for free surface flow modeling in unconfined aquifers. Eng Anal Bound Elem. 2024;158:214–229. https://doi.org/10.1016/j.enganabound.2023.10.024.
-
Bianchi M, Scheidegger J, Hughes A, et al. Simulation of national-scale groundwater dynamics in geologically complex aquifer systems: an example from Great Britain. Hydrol Sci J. 2024;69(5):572–591. https://doi.org/10.1080/02626667.2024.2320847.
-
Secci D, Godoy VA, Gómez-Hernández JJ. Physics-informed neural networks for solving transient unconfined groundwater flow. Comput Geosci. 2024;182:105494. https://doi.org/10.1016/j.cageo.2023.105494.
-
Alakayleh Z, Fang X, Clement TP. A comprehensive performance evaluation of physics-informed neural networks for groundwater flow modeling. Groundwater. 2024;62(2):177–192.
-
Dai H, Chen X, Ye M, et al. Deep learning in groundwater modeling: methods, applications and challenges. J Hydrol. 2025;650:132456.
-
Wang N, Chang H, Zhang D. Data-driven methods for flow and transport in porous media: a review. Int J Heat Mass Transf. 2024;235:126149.
-
Gotovac H, Malenica L, Gotovac B. Control volume isogeometric analysis for groundwater flow modeling in heterogeneous porous media. Adv Water Resour. 2021;148:103838. https://doi.org/10.1016/j.advwatres.2020.103838.
-
Chakraborty A, Rabinovich A, Chaudhuri A. Reduced-order modeling of groundwater flow in heterogeneous porous media. J Comput Phys. 2020;412:109426.
-
Ravshanov N, Aminov S, Shafiev T. Mathematical modeling of fluid filtration in porous media: numerical algorithms and software. Int J Appl Math. 2024;37(2):145–162.
-
Khuzhayorov B, Usmonov A, Nik Long NMA, et al. Anomalous solute transport in a two-zone heterogeneous porous medium. Appl Math Model. 2023;124:1–18.
-
Kresic N, Mikszewski A. Recent developments and trends in groundwater flow and transport modeling. Hydrogeol J. 2023;31:1249–1271.
-
Darcy H. Les Fontaines Publiques de la Ville de Dijon. Victor Dalmont; 1856.
-
Theis CV. The relation between the lowering of the piezometric surface and the rate and duration of discharge of a well using ground-water storage. Trans Am Geophys Union. 1935;16(2):519–524.
-
Bear J, Cheng AHD. Modeling Groundwater Flow and Contaminant Transport. Springer; 2010.
-
de Marsily G. Quantitative Hydrogeology: Groundwater Hydrology for Engineers. Academic Press; 1986.
-
Patankar SV. Numerical Heat Transfer and Fluid Flow. Hemisphere Publishing; 1980.
-
Peaceman DW. Interpretation of well-block pressures in numerical reservoir simulation. Soc Pet Eng J. 1978;18(3):183–194.
-
Hantush MS, Jacob CE. Non-steady radial flow in an infinite leaky aquifer. Trans Am Geophys Union. 1955;36(1):95–100.
-
von Neumann J, Richtmyer RD. A method for the numerical calculation of hydrodynamic shocks. J Appl Phys. 1950;21(3):232–237.
-
Lax PD, Richtmyer RD. Survey of the stability of linear finite difference equations. Commun Pure Appl Math. 1956;9(2):267–293.
About the Authors
License
Copyright (c) 2026 Uzbek Journal “Problems of Mechanics”
This work is licensed under a Creative Commons Attribution 4.0 International License.