7.3 Governing Equations for Unconfined Groundwater Flow

Unconfined flow equations are non-linear in that the transmissivity of the aquifer depends on the saturated thickness and the saturated thickness varies in the direction of flow because the water table slopes. By definition, a confined aquifer is fully saturated, so transmissivity values are constant for a confined aquifer of constant thickness. Under water table conditions the water table slopes, flow is parallel to the water table and the saturated thickness decreases in the direction of flow (Figure 55). If the slope of the water table is small, Darcy’s Law can be applied to develop governing equations by using the Dupuit simplification or Dupuit assumptions. Dupuit’s simplification uses the approximate gradient (dh over the distance x, –dh/dx) rather than the true gradient (dh over the flow path length, –dh/dL) by assuming the flow is horizontal (no vertical components of flow) as shown by the solid arrows in Figure 55 and by Equation 72.

\displaystyle Q=-K\frac{dh}{dx}A     rather than     \displaystyle Q=-K\frac{dh}{dL}A (72)

Expressing the flow area as the product of the height of the water table and the unit width of the system in the y direction, the flow is shown in Equation 73.

\displaystyle Q= -K\frac{dh}{dx}h\Delta y (73)
Figure showing the approximation used by Dupuit's assumption
Figure 55 – Dupuit’s simplification mathematically approximates unconfined flow as horizontal by using the gradient –dh/dx (blue solid arrows) instead of the gradient along the flow path –dh/dL (orange dashed lines). The calculated heads and flow rates are sufficiently accurate if the slope of the water table is small.

If the bottom of the unconfined aquifer is used as the datum, then the head defines the saturated thickness. To include this dependency in the flow equations, Equation 70 is adjusted so that aquifer thickness, b, is replaced with h, and the varying value of h has to be inside the derivative. In addition, to represent unconfined flow, specific yield, Sy, is used as the aquifer storativity. Thus, unconfined, two-dimensional (plan view), transient, anisotropic, heterogeneous conditions of groundwater flow are represented by Equation 74.

\displaystyle S_y\frac{\partial h}{\partial t}\ =\frac{\partial}{\partial x}K_x\left(h\frac{\partial h}{\partial x}\right)+\frac{\partial}{\partial y}K_y\left(h\frac{\partial h}{\partial y}\right) (74)

For unconfined, two-dimensional (plan view), transient, anisotropic, homogeneous conditions of groundwater flow, hydraulic conductivities do not need to be within the derivative, resulting in Equation 75.

\displaystyle S_y\frac{\partial h}{\partial t}\ =K_x\frac{\partial}{\partial x}\left(h\frac{\partial h}{\partial x}\right)+K_y\frac{\partial}{\partial y}\left(h\frac{\partial h}{\partial y}\right) (75)

Unconfined, two-dimensional, plan view, transient, isotropic, homogeneous flow is represented using only one value of K as shown in Equation 76.

\displaystyle S_y\frac{\partial h}{\partial t}\ =\frac{K}{2}\left(\frac{\partial^{2} h^2}{\partial x^2}+\frac{\partial^{2} h^2}{\partial y^2}\right) (76)

License

Hydrogeologic Properties of Earth Materials and Principles of Groundwater Flow Copyright © 2020 by The Authors. All Rights Reserved.