Modeling Groundwater Flow in Heterogeneous Porous Media with YAGMod

Computation

Laura Cattaneo¹  Alessandro Comunian¹  Giovanna de Filippis¹  Mauro Giudici¹  Chiara Vassena¹  Qinjun Kang² 

[1] Università degli Studi di Milano, Dipartimento di Scienze della Terra “A. Desio”, via Cicognara 7, I-20129 Milano, Italy;;Università degli Studi di Milano, Dipartimento di Scienze della Terra “A. Desio”, via Cicognara 7, I-20129 Milano, Italy

关键词: groundwater;    forward problem;    inverse problem;   

DOI  :  10.3390/computation4010002

来源: mdpi

【 摘 要 】

Modeling flow and transport in porous media requires the management of complexities related both to physical processes and to subsurface heterogeneity. A thorough approach needs a great number of spatially-distributed phenomenological parameters, which are seldom measured in the field. For instance, modeling a phreatic aquifer under high water extraction rates is very challenging, because it requires the simulation of variably-saturated flow. 3D steady groundwater flow is modeled with YAGMod (yet another groundwater flow model), a model based on a finite-difference conservative scheme and implemented in a computer code developed in Fortran90. YAGMod simulates also the presence of partially-saturated or dry cells. The proposed algorithm and other alternative methods developed to manage dry cells in the case of depleted aquifers are analyzed and compared to a simple test. Different approaches yield different solutions, among which, it is not possible to select the best one on the basis of physical arguments. A possible advantage of YAGMod is that no additional non-physical parameter is needed to overcome the numerical difficulties arising to handle drained cells. YAGMod also includes a module that allows one to identify the conductivity field for a phreatic aquifer by solving an inverse problem with the comparison model method.

【 授权许可】

CC BY   
© 2015 by the authors; licensee MDPI, Basel, Switzerland.