Permeability and Drainage in Soils

Water flow through soil is a key parameter in civil engineering design. It affects seepage under dams, stability of slopes, settlement of clay layers, and the structural integrity of pavements. Water trapped in pavement layers can lead to base softening, pumping of fines, and premature cracking. On the PE Civil Transportation exam, understanding permeability (hydraulic conductivity) and drainage design is essential.

1. Darcy’s Law

Water flows through soil pores due to differences in total energy (hydraulic head). The rate of water flow is governed by Darcy’s Law, which states that flow velocity is directly proportional to the hydraulic gradient.

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

Where:

  • $q = \text{flow rate (discharge volume per unit time, e.g., ft}^3/\text{sec or m}^3/\text{sec)}$
  • $k = \text{hydraulic conductivity (permeability, e.g., ft/day, cm/sec, m/s)}$
  • $i = \text{hydraulic gradient}$ (dimensionless head loss per unit distance of flow path): $$i = \frac{\Delta h}{L}$$ (where $\Delta h$ is the total head loss across the soil length $L$)
  • $A = \text{cross-sectional area of soil perpendicular to the direction of flow}$

Discharge Velocity vs. Seepage Velocity

  • Discharge Velocity ($v$): The apparent velocity of flow calculated over the gross cross-sectional area of the soil: $$v = \frac{q}{A} = k \cdot i$$
  • Seepage Velocity ($v_s$): The actual velocity of water as it flows through the pore channels. Since water can only flow through the void spaces, the seepage velocity is always greater than the discharge velocity: $$v_s = \frac{v}{n} = v \cdot \left(\frac{1+e}{e}\right)$$ Where:
    • $n = \text{porosity of the soil (expressed as a decimal)}$
    • $e = \text{void ratio of the soil}$

2. Equivalent Permeability in Layered Soils

In nature and engineered subgrades, soils are often stratified. To calculate the overall seepage through layered profiles, we determine the equivalent hydraulic conductivity.

       FLOW PARALLEL (Horizontal)                 FLOW PERPENDICULAR (Vertical)
      ============================                 ============================
      | Layer 1: H1, k1   -----> |                 | Layer 1: H1, k1   |
      ----------------------------                 --------------------|  | Flow
      | Layer 2: H2, k2   -----> |                 | Layer 2: H2, k2   |  v Direction
      ----------------------------                 --------------------|
      | Layer 3: H3, k3   -----> |                 | Layer 3: H3, k3   |
      ============================                 ============================

Flow Parallel to Soil Layers (e.g., Horizontal Seepage)

When flow is parallel to the stratification, the hydraulic gradient is the same in all layers. The equivalent permeability is a weighted average based on layer thicknesses:

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

Flow Perpendicular to Soil Layers (e.g., Vertical Seepage)

When flow is perpendicular to the stratification, the flow rate per unit area ($v$) is the same through all layers, but the hydraulic gradient varies. The equivalent permeability is a harmonic weighted average:

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

Key Note: The equivalent permeability parallel to layering ($k_{eq,p}$) is always greater than the equivalent permeability perpendicular to layering ($k_{eq,v}$).

3. Seepage and Flow Nets

A flow net is a graphical representation of the Laplace equation for two-dimensional flow through an isotropic medium. It consists of:

  1. Flow Lines: The pathways that water particles follow from upstream to downstream.
  2. Equipotential Lines: Lines connecting points of equal total head. Flow lines and equipotential lines must intersect at right angles ($90^\circ$) to form curvilinear squares.

Seepage Discharge from a Flow Net

For an isotropic soil, the total seepage flow rate per unit width of the structure is calculated as:

$$q = k \cdot H \cdot \frac{N_f}{N_d}$$

Where:

  • $k = \text{hydraulic conductivity}$
  • $H = \text{total head loss between upstream and downstream water levels}$
  • $N_f = \text{number of flow channels}$ (spaces between adjacent flow lines)
  • $N_d = \text{number of equipotential drops}$ (number of increments between equipotential lines)

4. Pavement Drainage and Subdrains

Good drainage is essential for extending pavement service life. Water entering the pavement structure from surface cracks or groundwater must be drained away quickly.

Pavement Drainage Components

  • Permeable Base Course: A layer of open-graded aggregate with high permeability ($k \approx 1,000$ to $3,000\text{ ft/day}$) located directly beneath the pavement slab to collect water.
  • Longitudinal Edge Drains: Perforated collector pipes placed in trenches along the edge of the pavement to intercept water flowing horizontally out of the base course.
  • Filter Fabric (Geotextile): Placed around the drain trench to prevent fine soil particles from clogging the drainage aggregate while allowing water to pass.

Hydraulic Transmissivity ($T$)

The capacity of a drainage layer to transmit water per unit width is defined as transmissivity:

$$T = k \cdot H_{\text{drain}}$$

Where $H_{\text{drain}}$ is the thickness of the permeable drainage layer.

5. Worked Examples

Worked Example 1: Darcy’s Law and Seepage Velocity

A constant-head permeability test is performed on a cylindrical sand specimen with a diameter of $4.0\text{ inches}$ and a length of $6.0\text{ inches}$. The head difference between the top and bottom of the specimen is maintained at $12.0\text{ inches}$. During a $5\text{-minute}$ test period, a total of $0.082\text{ ft}^3$ of water is collected. The dry weight of the sand is $5.20\text{ lb}$, and the specific gravity of solids is $2.65$.

Goal: Calculate the hydraulic conductivity ($k$) in $\text{ft/day}$ and the seepage velocity ($v_s$) in $\text{ft/day}$.

Solution:

  1. Calculate Flow Rate ($q$): $$q = \frac{\text{Volume}}{\text{Time}} = \frac{0.082\text{ ft}^3}{5\text{ min}} = 0.0164\text{ ft}^3/\text{min}$$ Convert to $\text{ft}^3/\text{day}$: $$q = 0.0164\text{ ft}^3/\text{min} \cdot 1440\text{ min/day} = 23.616\text{ ft}^3/\text{day}$$
  2. Calculate Cross-Sectional Area ($A$): Diameter $D = 4\text{ inches} = 0.333\text{ ft}$. $$A = \frac{\pi D^2}{4} = \frac{\pi (0.333)^2}{4} = 0.0873\text{ ft}^2$$
  3. Calculate Hydraulic Gradient ($i$): Head loss $\Delta h = 12\text{ inches}$, Length $L = 6\text{ inches}$. $$i = \frac{\Delta h}{L} = \frac{12.0}{6.0} = 2.0$$
  4. Calculate Hydraulic Conductivity ($k$): From Darcy’s Law: $q = k \cdot i \cdot A \implies k = \frac{q}{i \cdot A}$ $$k = \frac{23.616\text{ ft}^3/\text{day}}{2.0 \cdot 0.0873\text{ ft}^2} = 135.26\text{ ft/day}$$
  5. Calculate Porosity ($n$) to find Seepage Velocity:
    • Volume of solid particles ($V_s$): $$V_s = \frac{W_s}{G_s \cdot \gamma_w} = \frac{5.20\text{ lb}}{2.65 \cdot 62.4\text{ lb/ft}^3} = 0.0314\text{ ft}^3$$
    • Total volume of sample ($V = A \cdot L$): $$L = 6\text{ inches} = 0.50\text{ ft}$$ $$V = 0.0873\text{ ft}^2 \cdot 0.50\text{ ft} = 0.0437\text{ ft}^3$$
    • Porosity ($n$): $$n = \frac{V_v}{V} = \frac{V - V_s}{V} = \frac{0.0437 - 0.0314}{0.0437} = \frac{0.0123}{0.0437} = 0.281$$
  6. Calculate Seepage Velocity ($v_s$):
    • Discharge velocity ($v = k \cdot i$): $$v = 135.26\text{ ft/day} \cdot 2.0 = 270.52\text{ ft/day}$$
    • Seepage velocity ($v_s$): $$v_s = \frac{v}{n} = \frac{270.52\text{ ft/day}}{0.281} = 962.7\text{ ft/day}$$

Result: Hydraulic conductivity is $135.3\text{ ft/day}$ and seepage velocity is $962.7\text{ ft/day}$.

Worked Example 2: Equivalent Permeability in Layered Soils

A highway embankment foundation consists of two soil layers over an impermeable bedrock.

  • Layer 1 (Top): Silt, thickness $H_1 = 6\text{ ft}$, permeability $k_1 = 0.25\text{ ft/day}$.
  • Layer 2 (Bottom): Sand, thickness $H_2 = 4\text{ ft}$, permeability $k_2 = 85.0\text{ ft/day}$.

Goal: Determine the equivalent horizontal permeability ($k_{eq,h}$) and vertical permeability ($k_{eq,v}$) of the system.

Solution:

  1. Equivalent Horizontal Permeability ($k_{eq,h}$): $$k_{eq,h} = \frac{k_1 H_1 + k_2 H_2}{H_1 + H_2}$$ $$k_{eq,h} = \frac{(0.25 \cdot 6) + (85.0 \cdot 4)}{6 + 4} = \frac{1.5 + 340.0}{10} = \frac{341.5}{10} = 34.15\text{ ft/day}$$
  2. Equivalent Vertical Permeability ($k_{eq,v}$): $$k_{eq,v} = \frac{H_1 + H_2}{\frac{H_1}{k_1} + \frac{H_2}{k_2}}$$ $$\frac{H_1}{k_1} = \frac{6}{0.25} = 24.0\text{ days}$$ $$\frac{H_2}{k_2} = \frac{4}{85.0} = 0.047\text{ days}$$ $$k_{eq,v} = \frac{10}{24.0 + 0.047} = \frac{10}{24.047} = 0.416\text{ ft/day}$$

Result: Equivalent horizontal permeability is $34.15\text{ ft/day}$ (dominated by the sand layer) and vertical permeability is $0.416\text{ ft/day}$ (restricted by the silt layer).

6. Exam Pitfalls and Tips

  • Pore Velocity vs. Bulk Velocity: Make sure to distinguish between $v$ (discharge velocity) and $v_s$ (seepage velocity). PE questions will often ask for the “velocity of flow through the voids” or “water travel time,” which requires using $v_s$, not $v$.
  • Porosity vs. Void Ratio: In the seepage velocity equation $v_s = v/n$, remember that $n$ is porosity. If the problem gives you void ratio ($e$), you must convert it to porosity first using $n = e / (1+e)$.
  • Unit Conversions: Permeability is often given in laboratory units like $\text{cm/s}$, while the dimensions of the structure are in feet or meters. A standard conversion to memorize is: $$1.0\text{ cm/sec} \approx 2835\text{ ft/day}$$
  • Flow Net Counting: When counting flow channels ($N_f$) and equipotential drops ($N_d$) in a flow net, remember:
    • $N_f$ is the number of channels (lanes of flow), which is usually the number of flow lines minus one.
    • $N_d$ is the number of drops (steps in potential), which is the number of equipotential lines minus one (do not count the upstream and downstream boundary contours twice).