Skip to main content

Seepage with impervious foundation: derivation

Short answer

Under Dupuit’s assumptions and Darcy’s law, q dx = –k y dy integrates from headwater h₁ to tailwater h₂ over horizontal distance L to give specific discharge q = k(h₁² – h₂²) / (2L) and the phreatic surface y = √[h₁² – (h₁² – h₂²)x/L].

Derive the expression for the specific discharge through a homogeneous earthen dam section with an impervious foundation where the phreatic line emerges at the downstream tailwater.

Dupuit's assumptions: (i) the hydraulic gradient at any point on the phreatic line equals the slope of the line, i=dydxi=\dfrac{dy}{dx}; (ii) the gradient is constant over the whole vertical section; (iii) flow is horizontal and Darcy's law holds.

   h1 ___
      \  \ phreatic line         y
       \  \____                  |
        \       \____ h2 ........|  tail water
  ______ \  dam       \_____     |
 /////////// impervious base ////|///
  |<--------- L ---------------->|

Derivation

Take xx along the flow from the u/s end, and yy as the height of the phreatic line above the impervious base at distance xx. The line falls in the flow direction, so dy/dxdy/dx is negative. Per unit length of dam, the discharge through the vertical section of height yy is

q=k i A=−k dydx (y×1)q = k\, i\, A = -k\,\frac{dy}{dx}\,(y\times 1)

Hence

q dx=−k y dyq\,dx = -k\,y\,dy

Integrate from the u/s end (x=0x=0, y=h1y=h_1) to the exit at tail water (x=Lx=L, y=h2y=h_2), where h1h_1 and h2h_2 are the water depths above the base on the u/s and d/s side, and LL is the horizontal seepage length:

q∫0Ldx=−k∫h1h2y dyq L=k2(h12−h22)\begin{aligned} q\int_0^L dx &= -k\int_{h_1}^{h_2} y\,dy \\ q\,L &= \frac{k}{2}\left(h_1^2-h_2^2\right) \end{aligned}  q=k(h12−h22)2L \boxed{\,q=\frac{k\left(h_1^{2}-h_2^{2}\right)}{2L}\,}

The phreatic line is a parabola: integrating only up to the section xx gives q x=k2(h12−y2)q\,x = \frac{k}{2}(h_1^2-y^2), so

y=h12−(h12−h22)L xy=\sqrt{h_1^{2}-\frac{(h_1^{2}-h_2^{2})}{L}\,x}

Here LL is measured from the corrected entry point (0.3 times the horizontal projection of the wetted u/s slope inside the water line) to the point where the line meets the tail water. If the line emerges above tail water on the d/s face, the exit zone is treated by Schaffernak's or Casagrande's correction. For h2=0h_2=0 this becomes q=kh122Lq=\dfrac{k h_1^2}{2L}.

✓ Answer checked by Claude Opus 5.5 (high effort)

All Hydropower Engineering chapter 4 questions →