-
Notifications
You must be signed in to change notification settings - Fork 0
Home
GwMove is a physically-distributed process oriented subsurface groundwater flow simulator. It has been designed following modern / generic object oriented principles to allow for components reuse, community collaboration, and to trigger extensions as required by other application types in groundwater hydrology. The core components (backend) are written in C++ computer language. Currently, the code has the ability to simulate groundwater flow in three-dimensions including confined and unconfined systems. The upper aquifer unit, when unconfined, is modelled with a vertically moving mesh technique to iteratively adjust the water table position. The positions of the upper free nodes of the mesh are adjusted to water table dynamics and the positions of the underlying nodes are found by simple linear interpolation. The method has the advantage of keeping fixed the initial mesh topology during the iterative process, solving only a sequence of linear systems, completely avoiding nonlinear iterations and thus increasing computational robustness. Though an additional cost is required because the global conductance and storativity matrices are recomputed at every mesh iteration. However, this cost is not quite substantial and remains small by comparison to the time spent in the linear solver. GwMove uses a Galerkin finite element method for spatial discretization and a backward-Euler finite difference method for temporal discretization. It supports logically-oriented meshes denoted hereafter by IJK-oriented meshes. Although not as powerful as fully unstructured meshes, our implementation allows for relativeley complex geometries and irregular boundaries to be dealt with. Supported types of boundary conditions include fixed head, fixed flux, and head-dependent flux types. These are divided into 12 individual categories: (1) fixed head (2) fixed flow rate (3) fixed flux (4) river (5) leakage (6) infiltration (7) drainage (8) abstraction (9) seepage (10) recharge (11) evapotranspiration and (12) sea seepage face. Explanations, guidance and intuitive examples demonstrating the use of these boundary conditions for typical situations will be given in this wiki.