timflow.steady.inhomogeneity.LeakyBuildingPit3D#

class timflow.steady.inhomogeneity.LeakyBuildingPit3D(model, xy, kaq=1.0, z=None, kzoverkh=1.0, npor=0.3, topboundary='conf', topres=0.0, topthick=0.0, hstar=None, order=3, ndeg=3, layers=None, res=np.inf)#

Bases: LeakyBuildingPit

Element to simulate a building pit with a leaky wall in Model3D.

Parameters:

  • model (Model object) – model to which the element is added

  • xy (array or list) – list or array of (x,y) pairs of coordinates of corners of the inhomogeneity polygonal boundary is automatically closed (so first point is not repeated)

  • kaq (float, array or list) – hydraulic conductivity of each aquifer from the top down if float, hydraulic conductivity is the same in all aquifers

  • kzoverkh (float, array or list) – vertical anisotropy ratio vertical k divided by horizontal k if float, value is the same for all layers, length is number of layers

  • z (array or list) – elevation tops and bottoms of the aquifers from the top down leaky layers may have zero thickness if topboundary=’conf’: length is 2 * number of aquifers if topboundary=’semi’: length is 2 * number of aquifers + 1 as top of leaky layer on top of systems needs to be specified

  • topboundary (string, 'conf' or 'semi' (default is 'conf')) – indicating whether the top is confined (‘conf’) or semi-confined (‘semi’). For a building pit, the ‘conf’ option is generally more applicable.

  • topres (float) – resistance of top leaky layer (read if topboundary=”semi”)

  • topthick (float) – thickness of top semi-confining layer (read if topboundary=’semi’)

  • hstar (float or None (default is None)) – head value above semi-confining top, only read if topboundary=’semi’

  • npor (float, array or list) – porosity of all aquifers and leaky layers from the top down if float, porosity is the same for all layers if topboundary=’conf’: length is 2 * number of aquifers - 1 if topboundary=’semi’: length is 2 * number of aquifers

  • order (int) – polynomial order of flux along each segment

  • ndeg (int) – number of points used between two segments to numerically integrate normal discharge

  • layers (list or np.array) – layers in which leaky wall is present.

  • res (float or np.array, optional) – resistance of leaky wall, if passed as an array must be either shape (n_segments,) or (n_segments, n_layers). Default is np.inf, which simulates an impermeable wall.

Methods#

findlayer(z)

Returns layer-number, layer-type and model-layer-number.

summary()

Get summary of aquifer parameters.