timflow.steady.linesink.RadialCollectorWell#

class timflow.steady.linesink.RadialCollectorWell(model, x=0, y=0, narms=5, nls=10, L=10.0, angle=0, rcaisson=1.0, rw=0.1, Qw=100.0, res=0.0, layers=0, label=None)#

Bases: CollectorWell

Radial collector well.

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

Parameters:

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

  • x (float) – x-coordinate of the center of the collector well

  • y (float) – y-coordinate of the center of the collector well

  • L (float) – length of each arm

  • narms (int) – number of arms

  • rcaisson (float) – radius of the caisson

  • rw (float) – radius of the arms

  • nls (int) – number of line sinks per arm

  • Qw (float) – total discharge of the collector well

  • res (float) – resistance of the arms

  • layers (int, array or list) – layer(s) in which the well is screened

  • label (string, optional) – label of the collector well

Examples

Create a radial collector well with 5 arms:

ml = timflow.steady.Model3D(kaq=10, z=np.arange(20, -1, -2), kzoverkh=0.1)
w = timflow.steady.RadialCollectorWell(ml, x=0, y=0, narms=5, nls=10, angle=0,
                                       rcaisson=2.0, rw=0.1, Qw=1000, layers=5)
ml.solve()

Methods#

compute_xy(x, y, rcaisson, narms, nls, L, angle, ...)

Compute the x,y-coordinates array for the radial collector well.

discharge_per_arm()

Discharge of each arm.

potinf(x, y[, aq])

Compute unit potential influence of element.

disvecinf(x, y[, aq])

Returns array of size (2, nparam, naq).

discharge_per_linesink()

Discharge of the linesinks in each layer.

discharge()

Discharge of the element in each layer.

potinflayers(x, y, layers[, aq])

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

potentiallayers(x, y, layers[, aq])

Returns array of size len(layers) only used in building equations.

disvec(x, y[, aq])

Returns array of size (2, nparam, naq).

disvecinflayers(x, y, layers[, aq])

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

disveclayers(x, y, layers[, aq])

Returns two arrays of size len(layers) only used in building equations.