anaflow.flow.thiem

thiem(rad, r_ref, transmissivity, rate=- 0.0001, h_ref=0.0)[source]

The Thiem solution.

The Thiem solution for steady-state flow under a pumping condition in a confined and homogeneous aquifer. This solution was presented in [Thiem1906].

Parameters
  • rad (numpy.ndarray) – Array with all radii where the function should be evaluated.

  • r_ref (float) – Reference radius with known head (see h_ref).

  • transmissivity (float) – Transmissivity of the aquifer.

  • rate (float, optional) – Pumpingrate at the well. Default: -1e-4

  • h_ref (float, optional) – Reference head at the reference-radius r_ref. Default: 0.0

Returns

head – Array with all heads at the given radii.

Return type

numpy.ndarray

References

Thiem1906

Thiem, G., ‘’Hydrologische Methoden, J.M. Gebhardt’’, Leipzig, 1906.

Notes

The parameters rad, r_ref and transmissivity will be checked for positivity. If you want to use cartesian coordiantes, just use the formula r = sqrt(x**2 + y**2)

Examples

>>> thiem([1,2,3], 10, 0.001, -0.001)
array([-0.3664678 , -0.25615   , -0.19161822])