Hydraulics & Water Resources

A structured learning path from fluid fundamentals through advanced hydraulic design. Master the analysis of open channels, pipe systems, pumps, and hydraulic structures.

Start Learning Full Handbook

Level 1

Beginner — Fluid Properties and Hydrostatics

Start here if you are new to fluid mechanics.

Fluid Properties and Fundamental Concepts

Water, the primary fluid in hydraulic engineering, has well-defined physical properties that govern its behavior. Density (ρ = 1000 kg/m³ at 4°C), specific weight (γ = ρg = 9.81 kN/m³), dynamic viscosity (μ ≈ 1.0 × 10⁻³ Pa·s at 20°C), and kinematic viscosity (ν = μ/ρ ≈ 1.0 × 10⁻⁶ m²/s) are the fundamental fluid properties used in every hydraulic calculation. Viscosity varies significantly with temperature — a crucial factor in cold-climate hydraulics.

Flow regimes are classified as laminar or turbulent based on the Reynolds number (Re = ρVD/μ). Laminar flow (Re < 2000) is smooth and orderly, governed entirely by viscous forces. Turbulent flow (Re > 4000) is chaotic with eddies and mixing, governed primarily by inertial forces. The transition region (2000 < Re < 4000) is unstable and should be avoided in design. Most water engineering applications involve turbulent flow, with the exception of slow seepage through fine-grained soils.

Hydrostatic Pressure and Forces

Hydrostatic pressure increases linearly with depth according to Pascal's law: p = γh, where h is the depth below the free surface. The pressure distribution on a submerged vertical surface is triangular, with the resultant force acting at the center of pressure (located h/3 from the base for a rectangular surface). For inclined surfaces, the center of pressure is given by y_p = y_c + I_c/(y_c A), where y_c is the centroid depth and I_c is the second moment of area.

Hydrostatic forces are critical for designing retaining walls, dams, gates, and tanks. The total force on a submerged gate equals the product of pressure at the centroid and the area. Stability analysis of gravity dams requires calculating the resultant of hydrostatic forces (horizontal and vertical components) and checking that it falls within the middle third of the base. Buoyancy forces, governed by Archimedes' principle, affect submerged structures and pipelines.

Continuity, Energy, and Momentum Principles

Three fundamental conservation principles govern all hydraulic analysis. The continuity equation (Q = A₁V₁ = A₂V₂) expresses conservation of mass for incompressible flow. The Bernoulli equation (P₁/γ + V₁²/2g + z₁ = P₂/γ + V₂²/2g + z₂ + h_L) expresses conservation of energy, with head losses (h_L) accounting for friction and minor losses. The momentum equation (ΣF = ρQ(V₂ - V₁)) is essential for analyzing forces on hydraulic structures.

The energy grade line (EGL) and hydraulic grade line (HGL) are graphical representations of the Bernoulli equation along a flow system. The EGL plots total head (P/γ + V²/2g + z), while the HGL plots piezometric head (P/γ + z). The vertical distance between them represents velocity head (V²/2g). Drawing EGL and HGL profiles is an essential skill for identifying flow problems such as cavitation (where HGL drops below the pipe) and inadequate pressure heads.

Level 2

Intermediate — Open Channel Flow and Pipe Networks

Build on fundamentals with practical flow analysis.

Open Channel Flow and Manning's Equation

Open channel flow occurs when water flows with a free surface exposed to atmospheric pressure. Manning's equation (V = (1/n) R^(2/3) S^(1/2)) is the most widely used empirical formula for uniform flow in open channels. The Manning roughness coefficient (n) depends on the channel lining material: concrete (0.012-0.015), earth (0.020-0.030), gravel (0.025-0.035), and natural streams (0.030-0.050). Reliable Manning's n values are critical for accurate flow capacity estimation.

The most efficient cross-section for a given area is the one that maximizes the hydraulic radius (R = A/P). For a given area, a semicircular section has the maximum hydraulic radius, but trapezoidal sections with side slopes of 2:1 (H:V) to 3:1 are more practical. The most efficient trapezoidal section occurs when the hydraulic depth equals half the bottom width. Normal depth (y_n) is the depth at which flow is uniform for a given discharge, slope, and channel geometry. Use the Manning's Equation Calculator for channel design.

Specific Energy and Critical Flow

Specific energy (E = y + V²/2g) is the energy per unit weight of water measured relative to the channel bottom. For a given discharge, the specific energy curve shows two possible depths — subcritical (deep, slow) and supercritical (shallow, fast) — for the same specific energy. The critical depth (y_c) occurs at the minimum specific energy and corresponds to Froude number (Fr) = 1. The Froude number (Fr = V/√(gy)) characterizes flow regime: subcritical (Fr < 1), critical (Fr = 1), or supercritical (Fr > 1).

Critical depth is essential for flow measurement and channel design. For rectangular channels, y_c = (q²/g)^(1/3), where q = Q/b is the discharge per unit width. Controls (sections where critical depth occurs) govern the upstream water surface profile. The critical slope (S_c) is the slope at which normal depth equals critical depth for a given discharge. Mild slopes (S < S_c) produce subcritical flow; steep slopes (S > S_c) produce supercritical flow.

Pipe Flow and Network Analysis

Pipe flow analysis uses the Darcy-Weisbach equation (h_f = f (L/D) (V²/2g)) for head loss due to friction. The Darcy friction factor (f) depends on the Reynolds number and relative pipe roughness (ε/D) through the Moody chart or the Colebrook equation (1/√f = -2 log[ε/(3.7D) + 2.51/(Re√f)]). The Hazen-Williams equation (V = 0.849 C R^0.63 S^0.54) provides an alternative empirical approach widely used in water supply design, where C is the Hazen-Williams roughness coefficient.

Minor losses from fittings, valves, and bends are expressed as h_m = K V²/(2g), where K is the loss coefficient. Pipe network analysis for water distribution systems uses the Hardy Cross method (iterative balancing of head losses around loops) or the more efficient gradient method (Newton-Raphson solution of the system equations). Pump selection requires matching the system head-discharge curve to the pump performance curve, with attention to net positive suction head (NPSH) requirements to prevent cavitation. Use the Hazen-Williams Calculator and Pump Power Calculator for pipe system design.

Level 3

Advanced — Hydraulic Structures and Transient Flow

For senior students and practicing engineers.

Hydraulic Jump and Energy Dissipation

A hydraulic jump occurs when supercritical flow transitions to subcritical flow, dissipating kinetic energy through turbulence. The jump characteristics are governed by the upstream Froude number: undular jump (Fr = 1-1.7), weak jump (Fr = 1.7-2.5), oscillating jump (Fr = 2.5-4.5), steady jump (Fr = 4.5-9.0), and strong jump (Fr > 9.0). The sequent depth ratio (y₂/y₁) and energy loss (ΔE/E₁) are functions of the upstream Froude number.

Energy dissipation structures — stilling basins, baffle blocks, and riprap aprons — are designed to contain and stabilize the hydraulic jump. The USBR stilling basin types (I through IV) correspond to different Froude number ranges. Stilling basin length, tailwater depth requirements, and scour protection are critical design parameters. Jump location is controlled by the relationship between the sequent depth and the tailwater depth. Use the Hydraulic Jump Calculator for detailed analysis.

Weirs, Flumes, and Culverts

Weirs are overflow structures used for flow measurement and water level control. Sharp-crested weirs (rectangular, V-notch, or Cipolletti) provide accurate discharge measurement using standard equations: Q = C_d (2/3) √(2g) L H^(3/2) for rectangular weirs, and Q = C_d (8/15) √(2g) tan(θ/2) H^(5/2) for V-notch weirs. Broad-crested weirs are more robust and commonly used in irrigation systems. The discharge coefficient (C_d) accounts for approach velocity and contraction effects.

Parshall flumes measure flow in open channels by creating a critical flow section. The flume dimensions are standardized, and discharge is related to the upstream head through empirical equations. Culverts are closed conduits conveying flow under embankments. Culvert hydraulics depends on inlet control (the inlet limits capacity) versus outlet control (friction and backwater limit capacity). Inlet control analysis uses nomographs based on the US Federal Highway Administration's HDS-5 methods. Use the Weir Flow Calculator and Stormwater Runoff Calculator for design.

Water Hammer and Transient Flow

Water hammer (hydraulic transient) occurs when a sudden change in flow velocity — from valve closure, pump startup, or power failure — generates pressure waves that travel through the pipe system. The Joukowsky equation (ΔH = a ΔV/g) gives the maximum pressure rise, where a is the wave celerity. Wave celerity depends on pipe elasticity and fluid compressibility: a = √(K/ρ) / √(1 + (K/E)(D/e)), where K is water bulk modulus, E is pipe modulus, D is diameter, and e is wall thickness.

Transient control measures include surge tanks (provide a free surface to absorb pressure waves), air chambers (compressible air cushions), pressure relief valves, and controlled valve closure times. Surge analysis requires solving the water hammer equations (momentum and continuity partial differential equations) using the method of characteristics (MOC). For simple systems, the Allievi charts provide graphical solutions for valve closure transients. Use the Water Hammer Calculator to estimate surge pressures in pipeline systems.

Practice Exercises

Exercise 1: Hydrostatic Force on a Gate

A vertical rectangular gate 2 m wide and 1.5 m high is submerged in water with its top edge 3 m below the free surface. Calculate the total hydrostatic force on the gate and the depth to the center of pressure. Sketch the pressure distribution diagram.

Exercise 2: Open Channel Design

A trapezoidal channel with bottom width 3 m and side slopes 2:1 (H:V) carries a discharge of 15 m³/s on a longitudinal slope of 0.001. The Manning roughness coefficient is 0.025. Calculate the normal depth. Determine whether the flow is subcritical or supercritical by computing the Froude number. Verify with the Manning's Equation Calculator.

Exercise 3: Pipe Network Head Loss

Water flows through a 200 mm diameter cast iron pipe (ε = 0.26 mm) at a rate of 0.12 m³/s over a length of 500 m. Using the Darcy-Weisbach equation, calculate the friction head loss. Determine the friction factor using the Colebrook equation. Compare with the result from the Hazen-Williams Calculator.

Exercise 4: Hydraulic Jump

A rectangular channel 4 m wide carries 20 m³/s at a supercritical depth of 0.5 m. Calculate the upstream Froude number, the sequent depth after the hydraulic jump, and the energy dissipated in the jump. Use the Hydraulic Jump Calculator to verify your results.

References

Chow, V.T. Open-Channel Hydraulics. McGraw-Hill, 1959 (Classic reference).
Finnemore, E.J. and Franzini, J.B. Fluid Mechanics with Engineering Applications. 10th ed., McGraw-Hill, 2002.
Streeter, V.L., Wylie, E.B., and Bedford, K.W. Fluid Mechanics. 9th ed., McGraw-Hill, 1998.
Chaudhry, M.H. Applied Hydraulic Transients. 3rd ed., Springer, 2014.
USBR. Design of Small Dams. U.S. Bureau of Reclamation, 1987.
FHWA HDS-5. Hydraulic Design of Highway Culverts. Federal Highway Administration, 2012.
Civil Engineering Handbook — Hydraulics chapter with detailed design guidance.
Engineering Formula Library — Manning's, Hazen-Williams, and continuity formulas.
Engineering Glossary — Definitions of hydraulic engineering terms.