GeoPython2019

Modeling of Subsurface Flow and Transport with Dynamic Boundary Conditions
2019-06-25, 14:00–14:30, Room 2

A newly developed model makes groundwater modeling customizable by applying dynamic boundary conditions.
The user of the model can implement the behavior of such boundary conditions by writing Python plugins.


Boundary conditions are essential for groundwater models.
Typically, there are three major types of boundary conditions
(1) specified head (Dirichlet),
(2) specified flow (Neumann),
and (3) the combination of both former ones into a boundary condition of flow with resistance (Cauchy).
The user can specify values for these boundary conditions such as a well
at a certain location with a given pumping rate for a specified duration.
For some special applications, however, the specified values may further
depend on internal model conditions.
For example, the flow rate of an infiltration well that re-infiltrates water is equal to
the pumping rate of the extraction well.
This can be useful for geothermal applications within groundwater bodies.
The newly developed model, ueflow, allows the user to implement such a scheme by writing a plugin.
In addition to just using the pumping rate as infiltration rate, the user can incorporate
other constrains such as energy costs for pumping, capacities of water treatment facilities,
maintenance schedules for pumps based on pumping regimes, or other technical constrains.

The talk gives a short overview of ueflow that is based on the finite volume model framework
FiPy (Guyer et al. 2009).
FiPy is implemented in Python and offers multiple, high-performance solvers as well as
several tools for generating grids and other input data.
The next steps in developing ueflow involve adding a graphical user interface based on Python GIS libraries.
GeoPython is a good place to connect to scientists having similar tasks, leading to a fruitful exchange and collaboration for this open-source project.

Guyer, J. E., Wheeler, D., Warren, J. A. (2009). FiPy: Partial Differential Equations with Python. Computing in Science & Engineering 11(3) pp. 6—15 (2009), doi:10.1109/MCSE.2009.52, http://www.ctcms.nist.gov/fipy