◀ Back to Note

Tony's Intertwingled Memex • Reference Manual

Hydraulic Reservoir Surge Control (Civil Engineering) Solution

Formulas & Definitions

Problem Overview

We model the water accumulation rate inside a flood detention basin modeled as a horizontal triangular prism of fixed length L=100.0 mL = 100.0\text{ m}L=100.0 m and cross-sectional width-to-depth constraint w=hw = hw=h. Inflow during a storm is dVdt=25.0 m3/min\frac{dV}{dt} = 25.0\text{ m}^3/\text{min}dtdV​=25.0 m3/min. We calculate the vertical rise rate dhdt\frac{dh}{dt}dtdh​ at h=4.0 mh = 4.0\text{ m}h=4.0 m and evaluate emergency spillway bypass criteria.


Step-by-Step Mathematical Solution

Step 1: Establish the Volume Formula via Geometric Constraint

The volume of a triangular prism is the cross-sectional triangular area multiplied by the prism length:

V=(12⋅w⋅h)⋅LV = \left(\frac{1}{2} \cdot w \cdot h\right) \cdot LV=(21​⋅w⋅h)⋅L
(This formula and its geometric reduction principles are defined in Geometric Similarity in Related Rates and Related Rates.)

The basin length L=100.0 mL = 100.0\text{ m}L=100.0 m is a static constant (dLdt=0\frac{dL}{dt} = 0dtdL​=0). The cross-sectional geometry enforces the proportional constraint w=hw = hw=h. Substituting both L=100L = 100L=100 and w=hw = hw=h into the volume formula before differentiating:

V=12(h)(h)(100.0)=50.0h2V = \frac{1}{2}(h)(h)(100.0) = 50.0 h^2V=21​(h)(h)(100.0)=50.0h2

Step 2: Differentiate with Respect to Time

Differentiate both sides implicitly with respect to time ttt using the Power Rule and Chain Rule:

dVdt=ddt[50.0h2]=100.0hdhdt\frac{dV}{dt} = \frac{d}{dt}[50.0 h^2] = 100.0 h \frac{dh}{dt}dtdV​=dtd​[50.0h2]=100.0hdtdh​

Isolating the water level rise rate dhdt\frac{dh}{dt}dtdh​:

dhdt=1100.0hdVdt\frac{dh}{dt} = \frac{1}{100.0 h} \frac{dV}{dt}dtdh​=100.0h1​dtdV​

Step 3: Compute Rise Rate at Depth h=4.0 mh = 4.0\text{ m}h=4.0 m

Substitute dVdt=25.0 m3/min\frac{dV}{dt} = 25.0\text{ m}^3/\text{min}dtdV​=25.0 m3/min and h=4.0 mh = 4.0\text{ m}h=4.0 m:

dhdt=25.0100.0(4.0)=25.0400.0=116 m/min=0.0625 m/min\frac{dh}{dt} = \frac{25.0}{100.0 (4.0)} = \frac{25.0}{400.0} = \frac{1}{16}\text{ m/min} = 0.0625\text{ m/min}dtdh​=100.0(4.0)25.0​=400.025.0​=161​ m/min=0.0625 m/min

Converting to centimeters per minute:

dhdt=0.0625 m/min×100 cm/m=6.25 cm/min\frac{dh}{dt} = 0.0625\text{ m/min} \times 100\text{ cm/m} = 6.25\text{ cm/min}dtdh​=0.0625 m/min×100 cm/m=6.25 cm/min

Step 4: Emergency Bypass Protocol Evaluation

  1. Remaining Freeboard Depth to Spillway:
    Δh=Hspillway−h=6.0 m−4.0 m=2.0 m\Delta h = H_{\text{spillway}} - h = 6.0\text{ m} - 4.0\text{ m} = 2.0\text{ m}Δh=Hspillway​−h=6.0 m−4.0 m=2.0 m
  2. Additional Volume Required to Fill Basin from h=4.0 mh = 4.0\text{ m}h=4.0 m to h=6.0 mh = 6.0\text{ m}h=6.0 m:
    V(6.0)=50.0(6.0)2=50.0(36.0)=1,800.0 m3V(6.0) = 50.0(6.0)^2 = 50.0(36.0) = 1,800.0\text{ m}^3V(6.0)=50.0(6.0)2=50.0(36.0)=1,800.0 m3
    V(4.0)=50.0(4.0)2=50.0(16.0)=800.0 m3V(4.0) = 50.0(4.0)^2 = 50.0(16.0) = 800.0\text{ m}^3V(4.0)=50.0(4.0)2=50.0(16.0)=800.0 m3
    ΔV=1,800.0−800.0=1,000.0 m3\Delta V = 1,800.0 - 800.0 = 1,000.0\text{ m}^3ΔV=1,800.0−800.0=1,000.0 m3
  3. Exact Time to Spillway Overtopping at dVdt=25.0 m3/min\frac{dV}{dt} = 25.0\text{ m}^3/\text{min}dtdV​=25.0 m3/min:
    tovertop=ΔVdVdt=1,000.0 m325.0 m3/min=40.0 minutest_{\text{overtop}} = \frac{\Delta V}{\frac{dV}{dt}} = \frac{1,000.0\text{ m}^3}{25.0\text{ m}^3/\text{min}} = 40.0\text{ minutes}tovertop​=dtdV​ΔV​=25.0 m3/min1,000.0 m3​=40.0 minutes
  4. Safety Verdict: Because the true time to overtop is 40.0 minutes40.0\text{ minutes}40.0 minutes (which exceeds the 25.0 minute25.0\text{ minute}25.0 minute mandatory bypass threshold), emergency discharge gates do not need to be opened immediately. However, supervisory control systems must remain on high alert because if inflow surges above 40.0 m3/min40.0\text{ m}^3/\text{min}40.0 m3/min, the safety margin collapses.

Engineering Conclusion

In civil hydraulic engineering, retention basin cross-sections expand with depth: the horizontal water surface area at depth hhh is Asurface(h)=w⋅L=h⋅100=100hA_{\text{surface}}(h) = w \cdot L = h \cdot 100 = 100hAsurface​(h)=w⋅L=h⋅100=100h. Because dVdt=Asurface(h)⋅dhdt\frac{dV}{dt} = A_{\text{surface}}(h) \cdot \frac{dh}{dt}dtdV​=Asurface​(h)⋅dtdh​, the rate of vertical rise dhdt=1100hdVdt\frac{dh}{dt} = \frac{1}{100h}\frac{dV}{dt}dtdh​=100h1​dtdV​ decreases as the reservoir fills.

When water is shallow (h=1 mh = 1\text{ m}h=1 m), the level rises at 25 cm/min25\text{ cm/min}25 cm/min. At h=4 mh = 4\text{ m}h=4 m, the rise rate decelerates to 6.25 cm/min6.25\text{ cm/min}6.25 cm/min. Civil engineers use this related rates relationship in automated supervisory control and data acquisition (SCADA) networks. Acoustic depth sensors continuously measure hhh and compute dhdt\frac{dh}{dt}dtdh​, allowing algorithms to dynamically estimate time-to-crest and modulate automated sluice gates to prevent earthen embankment breach without unnecessarily flooding downstream wash channels.