Introduction#
Timflow steady is a submodule for the modeling of steady-state multi-layer groundwater flow with analytic elements.
Timflow steady may be applied to an arbitrary number of layers and an arbitrary sequence of aquifers and leaky layers. The head, flow, and leakage between aquifer layers may be computed analytically at any point in the aquifer system. The Dupuit approximation is applied to flow in aquifer layers (i.e., the resistance to flow in the vertical direction is neglected), while flow in leaky layers is approximated as vertical.
Tutorials and how-to guides for getting started with timflow steady.
Timflow steady fundamental concepts explained.
Timflow steady example notebooks.
Cross-sectional models.
Notebooks testing/benchmarking timflow steady implementations.
Quick Example#
In this example a well is modelled near a river in a single aquifer.
# import python packages
import numpy as np
import timflow.steady as tfs
# create model
ml = tfs.ModelMaq(kaq=10, z=[20, 0]) # single layer model
# add a river with a fixed water level
yls = np.arange(-100, 101, 20) # 20 points, so 19 segments
xls = 50 * np.ones_like(yls)
river = tfs.RiverString(ml, xy=list(zip(xls, yls)), hls=0.0)
# add a well
well = tfs.Well(ml, 0, 0, rw=0.3, Qw=1000)
# solve model
ml.solve()
# plot head contours
ml.plots.contour(win=[-30, 55, -30, 30], ngr=40, labels=True, decimals=1)
In this example a well is modelled near a river in a single aquifer.
