headwell: Create a analytic element of a well with a constant head

In raem: Analytic Element Modeling of Steady Single-Layer Groundwater Flow

Create a analytic element of a well with a constant head

Description

headwell() creates an analytic element of a well with a constant, specified head. The discharge into the well is computed by solving the corresponding aem model. The head can be specified at the well or at any other location.

Usage

headwell(xw, yw, hc, rw = 0.3, xc = xw, yc = yw, rc = rw, resistance = 0, ...)

Arguments

xw	numeric, x location of the well.
yw	numeric, y location of the well.
hc	numeric, specified hydraulic head at the collocation point.
rw	numeric, radius of the well. Defaults to 0.3 (meter).
xc	numeric, x location of the collocation point. See details. Defaults to xw.
yc	numeric, y location of the collocation point. See details. Defaults to yw.
rc	numeric, radius of the collocation point. See details. Defaults to rw.
resistance	numeric, hydraulic resistance at the collocation point. Defaults to 0 (no resistance).
...	ignored

Details

The discharge from the well at location xw - yw is computed by solving the aem model given the specified head hc. This head is specified at xc + rc - yc, called the collocation point. This can be used to compute the discharge of the well by specifying the head at some other location. By default, the location of the well and the collocation point are the same.

The hydraulic resistance of the well screen at the collocation point can be increased for a well in poor connection with the aquifer. If the aquifer is unconfined (i.e. has a variable saturated thickness), the system of equations becomes non-linear with respect to the hydraulic head and iteration is required to solve the model.

Value

Analytic element of a well with constant head which is an object of class headwell and inherits from well.

See Also

well()

Examples

hw <- headwell(xw = 400, yw = 300, hc = 20, rw = 0.3)
hw <- headwell(xw = 400, yw = 300, hc = 20, rw = 0.3, resistance = 10)
hw <- headwell(xw = 400, yw = 300, hc = 20, rw = 0.3, xc = 500, yc = 500, rc = 0)

raem documentation built on Sept. 11, 2024, 6:50 p.m.