CollectorWell#

class timflow.steady.linesink.CollectorWell(model, xy, Qw=100.0, rw=0.1, res=0.0, layers=0, order=0, label=None)#

Bases: DitchString

Collector well: collection of line sinks with a specified total discharge.

Collection of (discontinuous) line sinks with specified total discharge and unknown but uniform head.

Parameters:

model (Model object) – model to which the element is added
xy (np.array) – array of shape (N, 4) with start and end coordinates of the line sinks on each row: [(x1, y1, x2, y2), …]
Qw (float) – total discharge of the collector well
rw (float) – radius of the collector well arms
res (float) – resistance of the well screen
layers (int, array or list) – layer (int) or layers (list or array) where well is screened
order (int) – order of the line sink elements
label (string, optional) – label of the collector well

Examples

Create a collector well with two arms:

ml = timflow.steady.Model3D(kaq=10, z=np.arange(20, -1, -2), kzoverkh=0.1)
xy = [(1, 0, 10, 0), (0, 1, 0, 10)]
w = timflow.steady.CollectorWell(ml, xy=xy, Qw=1000, layers=np.arange(5, 10))
ml.solve()

potinf(x, y, aq=None)#

Compute unit potential influence of element.

Returns:

linesink 0, order 0, layer[0]: order 0, layer[1] … order 1, layer[0] order 1, layer[1] …
linesink 1, order 0, layer[0]: order 0, layer[1] … order 1, layer[0] order 1, layer[1] …

Return type:

array

disvecinf(x, y, aq=None)#: Returns array of size (2, nparam, naq).

discharge_per_linesink()#

Discharge of the linesinks in each layer.

Returns:: rv – array of shape (nlay, nlinesinks)
Return type:: np.array

discharge()#: Discharge of the element in each layer.

potinflayers(x, y, layers, aq=None)#

Returns array of size (len(layers),nparam).

Only used in building equations.

potentiallayers(x, y, layers, aq=None)#: Returns array of size len(layers) only used in building equations.

disvec(x, y, aq=None)#: Returns array of size (2, nparam, naq).

disvecinflayers(x, y, layers, aq=None)#

Returns two arrays of size (len(layers),nparam).

Only used in building equations.

disveclayers(x, y, layers, aq=None)#: Returns two arrays of size len(layers) only used in building equations.