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...
Saved in:
| Published in: | Environmental modelling & software : with environment data news Vol. 101; pp. 64 - 72 |
|---|---|
| Main Authors: | , , |
| Format: | Journal Article |
| Language: | English |
| Published: |
Oxford
Elsevier Ltd
01.03.2018
Elsevier Science Ltd |
| Subjects: | |
| ISSN: | 1364-8152, 1873-6726 |
| Online Access: | Get full text |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| Summary: | 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/.
•The model solves groundwater flow in unconfined aquifers with moving.boundary.•The model is easily tailored to be coupled to an overland flow model.•The GPGPU-accelerated implementation provides over 50-fold speedup.•Transient solutions in grids with 5 million points obtained in about 5 min.•Our high-performance implementation is open-source. |
|---|---|
| Bibliography: | SourceType-Scholarly Journals-1 ObjectType-Feature-1 content type line 14 ObjectType-Article-1 ObjectType-Feature-2 content type line 23 |
| ISSN: | 1364-8152 1873-6726 |
| DOI: | 10.1016/j.envsoft.2017.12.004 |