Seepage & Flow Nets

1. Introduction to Seepage

Seepage is the flow of water through soil under a hydraulic gradient. It occurs in virtually every geotechnical situation where water is present: beneath concrete dams, through earthen embankments, around sheet pile walls, under canal linings, and around basement walls. Seepage analysis answers three critical engineering questions:

• How much water is seeping through (quantity)?
• What is the pore water pressure and uplift force on structures?
• Is the exit gradient high enough to cause piping or erosion failure?

For simple one-dimensional flow (water flowing straight down through a layer), Darcy’s Law (q = kiA) is sufficient. For two-dimensional problems — flow under a dam, around a sheet pile, through a zoned embankment — the flow lines are curved and the head varies in two dimensions. The flow net is the graphical tool that solves these 2D problems.

2. Laplace’s Equation for Seepage

Combining Darcy’s Law with the continuity equation (conservation of mass) for steady, incompressible flow through a homogeneous isotropic porous medium gives Laplace’s equation:

Laplace’s equation has two orthogonal families of solution curves: equipotential lines (lines of constant head h) and flow lines (stream lines along which water flows). These two families of curves are always mutually perpendicular and together form the flow net.

For anisotropic soil (k_x ≠ k_z), the governing equation is: k_x ∂²h/∂x² + k_z ∂²h/∂z² = 0. This can be reduced to the Laplace form by a coordinate transformation (see Section 8).

3. Flow Net — Properties and Rules

3.1 Components

3.2 Rules for a Valid Flow Net

1. Flow lines and equipotential lines must intersect at right angles (90°) everywhere.
2. The mesh formed must be composed of curvilinear squares — each “square” has equal dimensions in both directions (Δl = Δb) so that the width-to-length ratio is 1:1 in each element.
3. Boundary conditions must be satisfied:
   • Upstream and downstream boundaries (free water surfaces in contact with soil) are equipotential lines.
   • Impermeable boundaries (bedrock, sheet pile, concrete base) are flow lines.
   • The phreatic surface (free surface in an earth dam) is both a flow line and an equipotential at atmospheric pressure.
4. The ratio N_f/N_d must be consistent — adding more flow lines requires proportionally more potential drops to maintain squares.

3.3 Typical N_f/N_d Values

N_f is rarely an integer — a fractional flow channel is acceptable (e.g., N_f = 3.7). This simply means the outermost channel carries a fraction of the flow of a full channel.

4. Flow Net Construction

Flow nets are traditionally drawn by hand using trial-and-error sketching (Casagrande’s method). Modern practice uses finite element or finite difference software, but hand-drawn flow nets remain important for understanding the physics and for quick estimates.

4.1 Steps for Hand-Drawn Flow Net

1. Draw the geometry to scale — soil body, structures, boundaries.
2. Identify and mark boundary conditions (upstream and downstream equipotentials; impermeable boundaries as flow lines).
3. Sketch 3–4 flow lines that smoothly curve from the upstream to the downstream boundary, respecting the impermeable boundaries.
4. Draw equipotential lines perpendicular to the flow lines, spacing them to form approximate squares.
5. Adjust iteratively — squeeze flow lines closer together where gradients are high (e.g., near corners of sheet piles), spread them apart where gradients are low.
6. Verify: all intersections are 90°; all elements are approximately square; boundary conditions are satisfied.

A correctly drawn flow net has N_f and N_d values that are consistent regardless of the level of refinement — halving the element size doubles both N_f and N_d, keeping N_f/N_d constant.

5. Seepage Quantity

5.1 Derivation

Consider one flow channel and one potential drop. The head drop per potential interval = H/N_d. For a curvilinear square of dimension b × b, the hydraulic gradient i = (H/N_d)/b and the flow area = b × 1 (per unit length). Flow in one square = k × (H/N_d)/b × b = kH/N_d. Total flow through all N_f channels: q = N_f × kH/N_d = kH(N_f/N_d). ✓

Note that q depends only on the ratio N_f/N_d, not on the absolute size of the flow net squares. This is why a coarser or finer flow net gives the same answer as long as the ratio is maintained.

6. Uplift Pressure and Seepage Force

6.1 Pressure Head at Any Point

At any point in the flow net, the total head h is known from the position of the equipotential line passing through it. The pore water pressure at that point is:

6.2 Head Loss Distribution

Each potential drop represents a head loss of Δh = H/N_d. If a point lies on the n-th equipotential line (counting from downstream), the residual head at that point = n × (H/N_d) above the downstream head.

6.3 Total Uplift Force on Base of Dam

6.4 Seepage Force

7. Exit Gradient and Piping

7.1 Exit Gradient

The exit gradient is the hydraulic gradient at the downstream exit face of the seepage flow (e.g., the toe of an embankment or the downstream side of a sheet pile). It is the most critical location because upward seepage here directly lifts soil particles against gravity.

7.2 Critical Hydraulic Gradient

7.3 Factor of Safety Against Piping

7.4 Piping Failure Mechanism

When i_e ≥ i_cr, the upward seepage force equals or exceeds the submerged weight of soil. Individual soil particles are lifted and carried away by the seeping water — a process called piping or internal erosion. Once a small pipe or channel forms, it grows progressively backward (called regressive or backward erosion) toward the upstream face, eventually creating a continuous channel through which water rushes uncontrolled. Piping is the most common cause of earthen dam failure worldwide and develops rapidly — often without visible warning.

8. Seepage Through Anisotropic Soil

Natural soils typically have higher horizontal permeability than vertical (k_x > k_z). The Laplace equation is not directly applicable, but the problem can be transformed into an isotropic one.

8.1 Coordinate Transformation

If k_x > k_z: the transformed domain is compressed in the x-direction (shorter, taller geometry). The flow net is drawn in this compressed space as if the soil were isotropic, then the seepage quantity is computed using k_eff.

9. Seepage Through an Earth Dam — Phreatic Line

In an earth dam, the top flow line (phreatic surface or line of saturation) is a free surface — it is both a flow line (no flow crosses it) and a boundary where pore pressure = atmospheric (u = 0, so h = z at all points on this line).

9.1 Casagrande’s Correction for the Phreatic Line

For a homogeneous earth dam with a horizontal downstream filter, Casagrande approximated the phreatic line as a parabola. The basic parabola has its focus at the lower end of the downstream slope toe.

The basic parabola entry is corrected at the upstream face (where the phreatic line must be tangent to the upstream slope) and at the downstream exit (where the line meets the slope or drain at a specific angle depending on the downstream slope angle α).

9.2 Seepage Through Earth Dam (Dupuit-Kozeny Formula)

10. Worked Examples

Example 1 — Seepage Quantity Under a Sheet Pile

Problem: A flow net for seepage under a sheet pile has N_f = 4 flow channels and N_d = 12 equipotential drops. The upstream head is 6 m and downstream head is 0 m. Soil permeability k = 4 × 10⁻⁴ m/s. Find the seepage quantity per metre length of sheet pile.

Example 2 — Uplift Pressure Under a Concrete Dam

Problem: A concrete dam sits on a 20 m wide soil foundation. Flow net: N_f = 4, N_d = 10. Upstream head = 8 m (above base), downstream head = 0. Find: (a) head drop per potential interval, (b) pore pressure at the mid-point of the dam base (5th equipotential from downstream), (c) total uplift force per unit length assuming linear pressure variation from heel to toe.

(a) Head Drop Per Interval

(b) Pore Pressure at Mid-Point (5th equipotential from downstream)

(c) Total Uplift Force

Note: Real uplift diagrams from flow nets are not linear but follow the equipotential distribution. Linear variation is a conservative approximation.

Example 3 — Exit Gradient and FOS Against Piping

Problem: A flow net for seepage under a dam shows the last (exit) potential square has a head drop of 0.9 m and a length of 0.6 m at the downstream exit. The soil has G_s = 2.66 and e = 0.70. Find: (a) exit gradient, (b) critical hydraulic gradient, (c) FOS against piping.

(a) Exit Gradient

(b) Critical Hydraulic Gradient

(c) FOS Against Piping

FOS < 1.0 → Piping failure is occurring. Immediate remedial measures needed: downstream filter blanket, relief wells, or upstream impermeable blanket to increase the seepage path.

Example 4 — GATE-Style: Seepage Through Stratified Section

Problem (GATE CE type): Seepage occurs vertically downward through two horizontal soil layers. Layer 1: H₁ = 2 m, k₁ = 6 × 10⁻⁴ m/s. Layer 2: H₂ = 3 m, k₂ = 2 × 10⁻⁵ m/s. Total head loss H = 1.0 m across both layers. Find: (a) flow rate per unit area, (b) head loss in each layer.

(a) Equivalent Vertical Permeability

(b) Head Loss in Each Layer

Nearly all the head loss (97.8 %) occurs in the low-permeability Layer 2, even though it is only 60 % of the total thickness — confirming that the less permeable layer dominates vertical seepage resistance.

11. Common Mistakes

Mistake 1: Treating N_f/N_d as N_d/N_f in the Seepage Formula

What happens: q = kH(N_d/N_f) is written instead of q = kH(N_f/N_d). Since N_d > N_f typically, this gives a seepage quantity that is much larger than the correct value.

Fix: Recall the logic: more flow channels (N_f) → more flow; more potential drops (N_d) → smaller gradient per drop → less flow. So q ∝ N_f/N_d. Think of it as: seepage quantity increases with flow channels and decreases with seepage path length (more N_d = longer effective path).

Mistake 2: Computing Exit Gradient Using Average Gradient Instead of Local Last-Square Gradient

What happens: i_e = H/L_total (average gradient over entire seepage path) is used instead of i_e = Δh/l (head drop in the last square / length of the last square). The local exit gradient is typically much higher than the average.

Fix: The exit gradient must always be the local gradient at the downstream exit face, read from the last potential square of the flow net. The average gradient bears no relation to piping risk.

Mistake 3: Forgetting that Each Potential Drop = H/N_d, Not H/N_f

What happens: Head loss per interval is computed as H/N_f (using number of flow channels) instead of H/N_d (using number of potential drops). This gives incorrect pore pressures and uplift forces.

Fix: Equipotential lines divide the total head loss into N_d equal drops. Flow lines divide the total flow into N_f equal channels. Head loss per drop = H/N_d. Flow per channel = q/N_f.

Mistake 4: Using k_H Instead of k_eff = √(k_xk_z) for Anisotropic Seepage

What happens: The horizontal permeability k_H is used in the seepage formula for an anisotropic soil, overestimating seepage quantity for cases where k_x > k_z.

Fix: For anisotropic soils, transform the domain (compress x by √(k_z/k_x)), draw the flow net in the transformed space, and use k_eff = √(k_xk_z) in the seepage formula.

Mistake 5: Treating the Phreatic Line as a Straight Line in Earth Dam Seepage

What happens: The phreatic surface in a homogeneous earth dam is assumed to be a straight inclined line from the upstream water level to the downstream face. This overestimates the seepage quantity and gives incorrect pressure distributions.

Fix: The phreatic line is a parabola (Casagrande’s approximation), not a straight line. In GATE problems involving earth dam seepage, always use the parabolic phreatic line or the Dupuit formula.

12. Frequently Asked Questions

Q1. Why must flow lines and equipotential lines intersect at exactly 90° in a valid flow net?

This is a mathematical consequence of the Laplace equation. Equipotential lines (h = constant) represent surfaces of equal energy. Water flows in the direction of maximum head decrease — i.e., perpendicular to the equipotential lines. This is directly analogous to heat flow perpendicular to isotherms, or electric current flow perpendicular to equipotential surfaces. If the intersection is not 90°, it implies that a component of flow is running parallel to an equipotential, which is physically impossible — there is no gradient driving flow along a constant-head surface. A flow net with non-perpendicular intersections either violates boundary conditions or contains drawing errors.

Q2. What are “curvilinear squares” and how do you check for them in a flow net?

A curvilinear square is a four-sided figure whose sides are curves (flow lines and equipotential lines) rather than straight lines, but which is approximately square in the sense that the average length equals the average width. Practically, you check by trying to inscribe a circle inside each element — if the circle fits snugly touching all four sides, the element is a curvilinear square. In practice, elements near corners (e.g., at the base of a sheet pile) may be quite distorted; the usual tolerance is that no element should have a length-to-width ratio exceeding 1.5:1. The overall accuracy of the flow net depends more on correctly satisfying boundary conditions than on achieving perfect squares.

Q3. How does a downstream filter blanket prevent piping?

A downstream filter blanket (or toe drain) placed at the exit face of a dam embankment prevents piping in two ways. First, it increases the seepage path length by directing seepage downward through a filter rather than allowing it to exit directly at the embankment face — this reduces the exit gradient i_e. Second, a properly designed filter (satisfying Terzaghi’s filter criteria: D₁₅(filter)/D₈₅(base) < 5 and D₁₅(filter)/D₁₅(base) > 4) retains fine soil particles while freely draining water, preventing the internal erosion that leads to piping. The filter acts as a protective blanket — water passes through it freely but soil particles cannot migrate into or through it.

Q4. What is the relationship between seepage, effective stress, and slope stability?

Seepage has a direct and often critical effect on slope stability. When water seeps through a slope toward the downstream face, it generates two effects: first, pore water pressure reduces effective stress (σ′ = σ − u), reducing shear strength. Second, the seepage force j = iγ_w acts in the direction of flow (often parallel to and toward the slope face), adding a disturbing force to the potential failure mass. Both effects reduce the factor of safety against sliding. Slopes that are stable when dry can fail when saturated — a common trigger for landslides after heavy rainfall. This is why seepage analysis (usually via a flow net or finite element seepage model) is always performed as part of slope stability analysis for embankments, natural slopes near water bodies, and cut slopes in areas of high groundwater.