A GPGPU-accelerated implementation of groundwater flow model in unconfined aquifers for heterogeneous and anisotropic media

The application of computer simulation models plays a significant role in the understanding of water dynamics in basins. The recent and explosive growth of the processing capabilities of General-Purpose Graphics Processing Units (GPGPUs) has resulted in widespread interest in parallel computing from the modelling community. In this paper, we present a GPGPU implementation of finite-differences solution of the equations of the 2D groundwater flow in unconfined aquifers for heterogeneous and anisotropic media. We show that the GPGPU-accelerated solution implemented using CUDA1 C/C++ largely outperforms the corresponding serial solution in C/C++. The results show that the GPGPU-accelerated implementation is capable of providing up to a 56-fold speedup in the solution using an ordinary office computer equipped with an inexpensive GPU2 card. The code developed for this research is available for download and use at http://modelagemambientaluffs.blogspot.com.br/.

[1]  Keith Beven,et al.  Using interactive recession curve analysis to specify a general catchment storage model , 1997 .

[2]  Luca Formaggia,et al.  Domain Decomposition Techniques , 2006 .

[3]  Xiaohui Ji,et al.  CUDA-based solver for large-scale groundwater flow simulation , 2011, Engineering with Computers.

[4]  S. Kanae,et al.  Estimating monthly total nitrogen concentration in streams by using artificial neural network. , 2011, Journal of environmental management.

[5]  J. Crank,et al.  A practical method for numerical evaluation of solutions of partial differential equations of the heat-conduction type , 1947 .

[6]  Kazuaki Hiramatsu,et al.  A comparison of techniques for hydrograph recession analysis , 2004 .

[7]  J. Crank,et al.  A practical method for numerical evaluation of solutions of partial differential equations of the heat-conduction type , 1947, Mathematical Proceedings of the Cambridge Philosophical Society.

[8]  D. Allen,et al.  Modeling coupled surface water – Groundwater processes in a small mountainous headwater catchment , 2014 .

[9]  E. Newman.  Water Movement Through Root Systems , 1976 .

[10]  Peter A. Forsyth,et al.  A parallel computational framework to solve flow and transport in integrated surface-subsurface hydrologic systems , 2012, Environ. Model. Softw..

[11]  Raghavan Srinivasan,et al.  A parallelization framework for calibration of hydrological models , 2012, Environ. Model. Softw..

[12]  R. Allan Freeze,et al.  Theoretical analysis of regional groundwater flow: 1. Analytical and numerical solutions to the mathematical model , 1966 .

[13]  Luca Ridolfi,et al.  Effect of vegetation–water table feedbacks on the stability and resilience of plant ecosystems , 2006 .

[14]  P. L. Libardi,et al.  Análise da determinação da condutividade hidráulica do solo pelo método do perfil instantâneo , 2013 .

[15]  D. E. Prudic,et al.  GSFLOW - Coupled Ground-Water and Surface-Water Flow Model Based on the Integration of the Precipitation-Runoff Modeling System (PRMS) and the Modular Ground-Water Flow Model (MODFLOW-2005) , 2008 .

[16]  P. Painlevé,et al.  Mémoire sur la transformation des équations de la Dynamique , 1894 .

[17]  Luca Ridolfi,et al.  Riparian vegetation distribution induced by river flow variability: A stochastic approach , 2006 .

[18]  A. W. Harbaugh,et al.  The History of MODFLOW , 2003, Ground water.

[19]  Ji Chen,et al.  Parallelization of a hydrological model using the message passing interface , 2013, Environ. Model. Softw..

[20]  R. Allan Freeze,et al.  Three-Dimensional, Transient, Saturated-Unsaturated Flow in a Groundwater Basin , 1971 .

[21]  Jun Zhou,et al.  Multi-GPU Implementation of a 3D Finite Difference Time Domain Earthquake Code on Heterogeneous Supercomputers , 2013, ICCS.

[22]  Marc F. P. Bierkens,et al.  Dynamic attribution of global water demand to surface water and groundwater resources: Effects of abstractions and return flows on river discharges , 2013 .

[23]  A. W. Harbaugh.  MODFLOW-2005 : the U.S. Geological Survey modular ground-water model--the ground-water flow process , 2005 .

[24]  Tom G. Chapman,et al.  Modeling groundwater flow over sloping beds , 1980 .

[25]  Carl Tim Kelley,et al.  Numerical simulation of water resources problems: Models, methods, and trends , 2013 .

[26]  Weixing Guo.  Transient groundwater flow between reservoirs and water-table aquifers , 1997 .

[27]  C. Colvin,et al.  A review of information on interactions between vegetation and groundwater , 1999 .

[28]  Ralph C. Heath,et al.  Basic ground-water hydrology , 1983 .

[29]  Kevin Skadron,et al.  Scalable parallel programming , 2008, 2008 IEEE Hot Chips 20 Symposium (HCS).

[30]  J. Wainwright,et al.  Environmental Modelling: Finding Simplicity in Complexity , 2013 .

[31]  Richard Barrett,et al.  Templates for the Solution of Linear Systems: Building Blocks for Iterative Methods , 1994, Other Titles in Applied Mathematics.

[32]  Kengo Nakajima.  Large-scale Simulations of 3D Groundwater Flow using Parallel Geometric Multigrid Method , 2013 .

[33]  J. Delleur.  The Handbook of Groundwater Engineering, Second Edition , 2006 .

[34]  Willem J. de Lange,et al.  An operational, multi-scale, multi-model system for consensus-based, integrated water management and policy analysis: The Netherlands Hydrological Instrument , 2014, Environ. Model. Softw..

[35]  P. Willemsen,et al.  Accelerating urban fast response Lagrangian dispersion simulations using inexpensive graphics processor parallelism , 2011, Environ. Model. Softw..

[36]  Praveen Kumar,et al.  GPU-based high-performance computing for integrated surface-sub-surface flow modeling , 2015, Environ. Model. Softw..

[37]  J. C. Ker.  LATOSSOLOS DO BRASIL: UMA REVISÃO , 1997 .

[38]  Raghavan Srinivasan,et al.  Efficient multi-objective calibration of a computationally intensive hydrologic model with parallel computing software in Python , 2013, Environ. Model. Softw..

[39]  Keith Beven,et al.  Including spatially variable effective soil depths in TOPMODEL , 1997 .

[40]  Xiaohui Ji,et al.  Parallelization of MODFLOW Using a GPU Library , 2014, Ground water.

[41]  M. Parlange,et al.  Recession flow analysis for aquifer parameter determination , 1998 .

[42]  Guillermo A. Mendoza,et al.  Multi-criteria decision analysis in natural resource management: A critical review of methods and new modelling paradigms , 2006 .

[43]  H. V. D. Vorst,et al.  Conjugate gradient type methods and preconditioning , 1988 .

[44]  Timothy R. Ginn,et al.  Richards Equation–Based Modeling to Estimate Flow and Nitrate Transport in a Deep Alluvial Vadose Zone , 2012 .