9.3 Soil Permeability, Seepage, Flow Nets, and Effective Stress

9.3 Soil Permeability, Seepage, Flow Nets, and Effective Stress

Water flow through soil media dictates critical geotechnical design problems—including seepage losses beneath concrete gravity dams, hydraulic uplift forces on basements and dry docks, piping erosion at weir toes, slope stability under transient seepage, and 1D consolidation settlement under building footings. In the CELE Hydraulics and Geotechnical Engineering examination cluster, seepage and effective stress principles represent the most analytically rigorous problem sets.

---

Subsurface Water & Darcy's Law

Subsurface soil water is classified into three physical categories:

1. Gravitational Water (Free Water): Water that moves through interconnected void channels under the driving influence of gravity and hydraulic gradients, governed by fluid flow principles.
2. Capillary Water: Water held in void pores above the phreatic surface (groundwater table) by surface tension forces acting along soil grain menisci. Capillary rise induces negative pore water pressure (matric suction): $u = -\gamma_w h_c$, which mathematically increases effective stress and apparent cohesion.
3. Adsorbed (Hygroscopic) Water: Dipolar water molecules held tenaciously to the electrically charged surfaces of clay minerals by hydrogen bonding and van der Waals forces. Adsorbed water does not circulate freely and cannot transmit hydrostatic pressure.

Total Hydraulic Head & Darcy's Formulations

In geotechnical hydraulics, the total head ($h$) at any subsurface point is expressed through Bernoulli's equation. Because groundwater seepage occurs at very low velocities, the kinetic velocity head term ($\frac{v^2}{2g} \approx 0$) is negligible:

$h = \frac{u}{\gamma_w} + z$

where $\frac{u}{\gamma_w}$ is the pressure head ($h_p$), $z$ is the elevation head relative to a chosen datum, and $h$ is the piezometric (total) head.

In 1856, Henry Darcy formulated the empirical governing law for laminar flow through saturated sand filters:

$v = k \cdot i$

$q = v \cdot A = k \cdot i \cdot A$

where:

• $v$ = discharge velocity (or superficial Darcy velocity), calculated across the gross cross-sectional area $A$ of the soil specimen.
• $k$ = hydraulic conductivity (or coefficient of permeability) with units of velocity ($cm/s$ or $m/s$).
• $i$ = hydraulic gradient, defined as head loss per unit flow length: $i = \frac{\Delta h}{L}$.

Discharge Velocity vs. Actual Seepage Velocity

Water does not flow through the solid soil skeleton; it is constrained strictly to the tortuous void pathways. By equating total discharge across gross area $A$ and void area $A_v$:

$q = v \cdot A = v_s \cdot A_v$

$v_s = v \cdot \frac{A}{A_v} \approx v \cdot \frac{V}{V_v} = \frac{v}{n} = v \left(\frac{1+e}{e}\right)$

where $v_s$ is the seepage velocity (actual interstitial pore velocity) and $n$ is porosity. Because $n < 1.0$, the true seepage velocity is strictly greater than Darcy's discharge velocity ($v_s > v$).

---

Laboratory Permeability Testing Methods

Laboratory hydraulic conductivity is measured using two standardized testing apparatuses depending on soil texture:

1. Constant Head Permeability Test (ASTM D2434)

Utilized for coarse-grained soils with relatively high hydraulic conductivity ($k > 10^{-4}\text{ cm/s}$, such as clean sands and gravels). A constant head difference ($h$) is maintained across a soil specimen of length $L$ and area $A$. The water volume $Q$ collected over elapsed time $t$ is measured:

$Q = q \cdot t = (k \cdot i \cdot A) t = k \left(\frac{h}{L}\right) A \cdot t$

$k = \frac{Q \cdot L}{A \cdot h \cdot t}$

2. Falling (Variable) Head Permeability Test (ASTM D5856)

Utilized for fine-grained soils with low permeability ($k < 10^{-4}\text{ cm/s}$, such as silts and clays). Water flows from a narrow standpipe of area $a$ through a soil specimen of area $A$ and length $L$. As water seeps into the soil, the hydraulic head in the standpipe drops from $h_1$ at time $t = 0$ to $h_2$ at time $t$:

$\text{Inflow Rate: } q_{in} = -a \frac{dh}{dt}$

$\text{Outflow Rate (Darcy): } q_{out} = k \left(\frac{h}{L}\right) A$

Equating inflow and outflow gives the governing separable differential equation:

$-a \frac{dh}{dt} = k \left(\frac{h}{L}\right) A \implies -\frac{dh}{h} = \left(\frac{k A}{a L}\right) dt$

Integrating between limits $[h_1, h_2]$ and $[0, t]$:

$-\int_{h_1}^{h_2} \frac{dh}{h} = \frac{k A}{a L} \int_0^t dt \implies \ln\left(\frac{h_1}{h_2}\right) = \frac{k A t}{a L}$

$k = \frac{a \cdot L}{A \cdot t} \ln\left(\frac{h_1}{h_2}\right) = 2.303 \frac{a \cdot L}{A \cdot t} \log_{10}\left(\frac{h_1}{h_2}\right)$

---

Equivalent Permeability in Stratified / Layered Soils

Natural sedimentary and alluvial soil deposits are horizontally stratified. Because permeability varies drastically across consecutive strata, equivalent hydraulic conductivities must be derived for directional seepage.

1. Flow Parallel to Stratification ($k_{h,eq}$)

When groundwater flows horizontally parallel to horizontal soil layers of individual thicknesses $H_1, H_2, \dots, H_n$ and permeabilities $k_1, k_2, \dots, k_n$:

• The hydraulic gradient is uniform across all layers: $i_1 = i_2 = \dots = i$.
• Total discharge is the sum of individual layer discharges: $q = \sum q_i = \sum (v_i A_i) = \sum (k_i i H_i \cdot 1)$.

$k_{h,eq} = \frac{\sum k_i H_i}{\sum H_i} = \frac{k_1 H_1 + k_2 H_2 + \dots + k_n H_n}{H_1 + H_2 + \dots + H_n}$

2. Flow Perpendicular to Stratification ($k_{v,eq}$)

When groundwater seeps vertically through horizontal soil layers:

• The discharge velocity is identical across all layers to preserve fluid mass continuity: $v_1 = v_2 = \dots = v$.
• Total head loss is the sum of individual layer head losses: $\Delta h = \sum \Delta h_i = \sum (i_i H_i) = \sum \left(\frac{v H_i}{k_i}\right)$.

$k_{v,eq} = \frac{\sum H_i}{\sum \left(\frac{H_i}{k_i}\right)} = \frac{H}{\frac{H_1}{k_1} + \frac{H_2}{k_2} + \dots + \frac{H_n}{k_n}}$

---

Two-Dimensional Seepage & Flow Nets

Combining Darcy's law with the 2D fluid continuity equation for an incompressible fluid moving through isotropic, homogeneous soil yields the classical Laplace Equation of steady-state seepage:

$\frac{\partial^2 h}{\partial x^2} + \frac{\partial^2 h}{\partial z^2} = 0$

A flow net is a graphical solution to Laplace's equation comprising two mutually orthogonal families of curves:

1. Flow Lines ($\Psi$): Streamlines tracing the trajectory of seeping water particles. The region bounded between two adjacent flow lines is a flow channel. Total number of flow channels = $N_f$.
2. Equipotential Lines ($\Phi$): Contours connecting points of equal total hydraulic head ($h$). The head drop between any two adjacent equipotential lines is the potential drop $\Delta h = \frac{H}{N_d}$, where $H$ is the total driving head (upstream minus downstream) and $N_d$ is the total number of potential drops.

Upstream Water (H1)                     Downstream Tailwater (H2)
   |                                                |
===+========= Concrete Dam / Impervious Base =======+===
   |                                                |
   +~~~~~~~~~~~~~~~~ Soil Foundation ~~~~~~~~~~~~~~~+
   |  -> Flow Line 1 --------------------------->   |  (Nf = Flow Channels)
   |  -> Flow Line 2 --------------------------->   |
   |   |        |        |        |        |        |  (Nd = Potential Drops)
   |  Φ1       Φ2       Φ3       Φ4       Φ5       Φ6

1. Seepage Quantity Under Hydraulic Structures

For a flow net constructed with curvilinear squares (element aspect ratio $\frac{\Delta s}{\Delta l} \approx 1.0$), the seepage discharge ($q$) per unit length ($1.0\text{ m}$ perpendicular to the cross-section) is:

$q = k \cdot H \left(\frac{N_f}{N_d}\right)$

2. Uplift Pressure Under Concrete Structures

At any point $i$ along the base of a concrete weir or gravity dam, the pore water pressure ($u_i$) is computed from total head and elevation:

$h_i = h_{upstream} - n_i \cdot \left(\frac{H}{N_d}\right)$

$u_i = (h_i - z_i)\gamma_w$

where $n_i$ is the number of potential drops experienced from the upstream reservoir to point $i$, and $z_i$ is the elevation of the base point above the chosen datum.

3. Exit Gradient & Quicksand / Piping Phenomenon

At the downstream toe of a hydraulic structure, water emerges vertically. The exit hydraulic gradient ($i_{exit}$) across the final flow field of length $l_{exit}$ is:

$i_{exit} = \frac{\Delta h}{l_{exit}} = \frac{H / N_d}{l_{exit}}$

As upward seepage force increases, it counteracts the gravitational submerged weight of the soil particles. The critical hydraulic gradient ($i_{cr}$) occurs when the upward seepage force per unit volume ($j = i \gamma_w$) exactly equals the submerged unit weight ($\gamma'$):

$i_{cr} \cdot \gamma_w = \gamma' \implies i_{cr} = \frac{\gamma'}{\gamma_w} = \frac{\gamma_{sat} - \gamma_w}{\gamma_w} = \frac{G_s - 1}{1+e}$

When $i_{exit} \ge i_{cr}$, the effective stress drops to zero ($\sigma' = 0$). Cohesionless soil completely loses its internal shear strength, transforming into a boiling suspension known as quicksand. This triggers subsurface retrogressive erosion called piping. The design Factor of Safety against piping must satisfy:

$FS_{piping} = \frac{i_{cr}}{i_{exit}} \ge 3.0 \text{ to } 5.0$

---

Terzaghi's Principle of Effective Stress

Formulated by Karl Terzaghi in 1925, the Principle of Effective Stress is the single most important concept in geotechnical engineering. Soil deformations (compression, consolidation) and shearing resistance are governed not by total stress alone, but by intergranular stresses transmitted across solid grain contact points:

$\sigma = \sigma' + u \implies \sigma' = \sigma - u$

where:

• $\sigma$ = total vertical stress, representing the total weight of soil and water above a given depth per unit area: $\sigma = \sum \gamma_i h_i$.
• $u$ = pore water pressure, the hydrostatic fluid pressure in void spaces: $u = \gamma_w h_w$.
• $\sigma'$ = effective vertical stress, representing the net skeletal stress resisting deformation.

Seepage Influence on Effective Stress Profiles

1. Hydrostatic (No Flow) Condition: $\sigma = \gamma_{sat} \cdot z$ $u = \gamma_w \cdot z$ $\sigma' = \sigma - u = (\gamma_{sat} - \gamma_w) z = \gamma' \cdot z$
2. Upward Seepage Condition (Hydraulic gradient $i$): Upward flow exerts an upward drag on soil particles, increasing pore pressure by $\Delta u = +i z \gamma_w$: $\sigma' = \gamma' z - i z \gamma_w = z(\gamma' - i \gamma_w)$
   • At $i = i_{cr} = \frac{\gamma'}{\gamma_w}$, effective stress vanishes ($\sigma' = 0$), producing boiling/liquefaction.
3. Downward Seepage Condition (Hydraulic gradient $i$): Downward flow exerts a downward drag on soil particles, decreasing pore pressure by $\Delta u = -i z \gamma_w$: $\sigma' = \gamma' z + i z \gamma_w = z(\gamma' + i \gamma_w)$
   • Downward seepage enhances soil stability and effective stress.

---

CELE Board-Exam Worked Situational Problem

Problem Statement

A concrete gravity dam retains an upstream head of water $H_1 = 12.0\text{ m}$ and a downstream tailwater head $H_2 = 1.5\text{ m}$. The dam foundation rests on a homogeneous sand layer with thickness $H = 9.0\text{ m}$, underlain by impermeable basalt bedrock. The sand has hydraulic conductivity $k = 4.0 \times 10^{-5}\text{ m/s}$, void ratio $e = 0.65$, and solid specific gravity $G_s = 2.66$. Take $\gamma_w = 9.81\text{ kN/m}^3$.

A flow net drawn for the foundation geometry reveals $N_f = 4$ flow channels and $N_d = 12$ equipotential drops. The final flow element at the downstream exit has a length of $l_{exit} = 1.60\text{ m}$.

Calculate:

1. The total under-seepage rate ($q$) per meter length of dam in cubic meters per day ($m^3/day/m$).
2. The critical hydraulic gradient ($i_{cr}$) and the exit hydraulic gradient ($i_{exit}$) at the downstream toe.
3. The Factor of Safety ($FS$) against piping, and the effective vertical stress ($\sigma'$) at a depth of $z = 3.0\text{ m}$ directly beneath the downstream toe under this upward seepage condition.

Step-by-Step Solution

Part 1: Under-Seepage Rate Calculation

1. Net Driving Head ($H$): $H = H_1 - H_2 = 12.0\text{ m} - 1.5\text{ m} = 10.5\text{ m}$
2. Seepage Discharge per Unit Width ($q$): $q = k \cdot H \left(\frac{N_f}{N_d}\right) = (4.0 \times 10^{-5}\text{ m/s}) \times 10.5\text{ m} \times \left(\frac{4}{12}\right)$ $q = 4.0 \times 10^{-5} \times 10.5 \times 0.3333 = 1.400 \times 10^{-4}\text{ m}^3\text{/s/m}$
3. Conversion to Daily Discharge: $q_{day} = 1.400 \times 10^{-4}\text{ m}^3\text{/s} \times 86,400\text{ s/day} = 12.096\text{ m}^3\text{/day/m} \approx 12.10\text{ m}^3\text{/day/m}$

Part 2: Critical and Exit Hydraulic Gradients

1. Critical Hydraulic Gradient ($i_{cr}$): $i_{cr} = \frac{G_s - 1}{1+e} = \frac{2.66 - 1}{1 + 0.65} = \frac{1.66}{1.65} = 1.0061 \approx 1.01$
2. Head Drop per Equipotential Step ($\Delta h$): $\Delta h = \frac{H}{N_d} = \frac{10.5\text{ m}}{12} = 0.875\text{ m}$
3. Exit Gradient ($i_{exit}$): $i_{exit} = \frac{\Delta h}{l_{exit}} = \frac{0.875\text{ m}}{1.60\text{ m}} = 0.5469 \approx 0.547$

Part 3: Factor of Safety and Effective Stress

1. Factor of Safety against Piping ($FS$): $FS = \frac{i_{cr}}{i_{exit}} = \frac{1.0061}{0.5469} = 1.84$ (Note: Because $FS = 1.84 < 3.0$, the dam toe is vulnerable to piping; downstream protective filters or relief wells are mandatory).
2. Submerged Unit Weight of Sand ($\gamma'$): $\gamma_{sat} = \frac{(G_s + e)\gamma_w}{1+e} = \frac{(2.66 + 0.65) \times 9.81}{1.65} = \frac{3.31 \times 9.81}{1.65} = 19.68\text{ kN/m}^3$ $\gamma' = \gamma_{sat} - \gamma_w = 19.68 - 9.81 = 9.87\text{ kN/m}^3$
3. Effective Vertical Stress with Upward Seepage at $z = 3.0\text{ m}$: At the downstream toe, seepage is directed upward with gradient $i \approx i_{exit} = 0.5469$: $\sigma' = z(\gamma' - i \gamma_w) = 3.0\text{ m} \times [9.87\text{ kN/m}^3 - (0.5469 \times 9.81\text{ kN/m}^3)]$ $\sigma' = 3.0 \times [9.87 - 5.365] = 3.0 \times 4.505 = 13.52\text{ kPa}$ Hydrostatic check (without seepage): $\sigma'_{hydro} = 3.0 \times 9.87 = 29.61\text{ kPa}$. Notice that upward seepage has reduced the in-situ effective stress by over 54%.

---

CELE Board Examination Traps & Critical Pitfalls