LeakyWall#

class timflow.steady.linedoublet.LeakyWall(model, x1=-1, y1=0, x2=1, y2=0, res=0, order=0, layers=0, label=None, addtomodel=True)#

Bases: LineDoubletHoBase, timflow.steady.equation.LeakyWallEquation

Create a segment of a leaky wall, which is simulated with a line-doublet.

The specific discharge through the wall is equal to the head difference across the wall divided by the resistance of the wall.

Parameters:

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

  • x1 (scalar) – x-coordinate of fist point of line-doublet

  • y1 (scalar) – y-coordinate of fist point of line-doublet

  • x2 (scalar) – x-coordinate of second point of line-doublet

  • y2 (scalar) – y-coordinate of second point of line-doublet

  • res (scalar) – resistance of leaky wall

  • order (int (default is 0)) – polynomial order of potential jump along line-doublet (head jump if transmissivity is equal on each side of wall)

  • layers (scalar, list or array) – layer(s) in which element is placed if scalar: element is placed in this layer if list or array: element is placed in all these layers

  • label (str or None) – label of element

See also

LeakyWallString

potinf(x, y, aq=None)#

Can be called with only one x,y value.

Returns:

(nparam, self.aq.naq) with order order 0, layer[0] order 0, layer[1] … order 1, layer[0] order 1, layer[1] etc.

Return type:

array

disvecinf(x, y, aq=None)#

Can be called with only one x,y value.

Returns:

(nparam, self.aq.naq) with order order 0, layer[0] order 0, layer[1] … order 1, layer[0] order 1, layer[1] etc.

Return type:

array

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.

equation()#

Mix-in class that returns matrix rows for leaky wall condition.

Qnormal = resfac * (headin - headout)

Returns:

  • matrix – (nunknowns,neq)

  • rhs – (nunknowns)