Hydraulic Analysis and Modelling

Introduction and Context

Hydraulic structures—including spillways, stilling basins, intake works, canals, culverts, energy-dissipation systems, and pipelines—must safely convey, control, measure, or dissipate water. Their performance depends on the interaction between flow resistance, geometry, gravity, pressure, inertia, turbulence, and transient effects. Hydraulic analysis converts these physical processes into equations, models, and design decisions.

This chapter develops the principal analytical and modelling tools used for hydraulic structures. It considers open-channel and pressurized flow, gradually and rapidly varied flow, hydraulic jumps, cavitation, dimensional analysis, similitude, physical modelling, and numerical modelling. Emphasis is placed on assumptions, validity, interpretation, and professional judgement rather than on formula substitution alone.

Learning Objectives

On completing this chapter, the reader should be able to:

1. distinguish open-channel flow from pressurized flow and select appropriate governing equations;

2. analyse uniform, gradually varied, and rapidly varied flow;

3. determine hydraulic-jump characteristics and energy losses;

4. assess cavitation risk using pressure and cavitation criteria;

5. apply dimensional analysis and similitude to hydraulic problems;

6. select, design, and interpret physical and numerical hydraulic models; and

7. evaluate uncertainty, limitations, and ethical responsibilities in hydraulic modelling.

1. Governing Principles and Hydraulic Variables

The fundamental conservation laws are conservation of mass, momentum, and energy. For a one-dimensional steady control volume, continuity is

Q = AV

where Q is discharge (m³/s), A is flow area (m²), and V is mean velocity (m/s). This relation assumes that the velocity distribution has been represented by a suitable mean value.

The total head at a section is commonly written as

H = z + y + αV²/(2g)

for an open channel, where z is bed elevation, y is flow depth, g is gravitational acceleration, and α is the kinetic-energy correction coefficient. In a pressurized conduit, y is replaced by pressure head p/(ρg), so that

H = z + p/(ρg) + αV²/(2g).

Here p is gauge pressure and ρ is fluid density. Between two sections, the energy equation is

H₁ + h_p = H₂ + h_L + h_t,

where h_p is head added by a pump, h_t is head extracted by a turbine, and h_L is total head loss. For a system without machines, h_p = h_t = 0. Energy loss is non-negative in ordinary passive flow.

The momentum equation is especially important where energy is dissipated irreversibly, as in hydraulic jumps. In vector form for a control volume,

ΣF = ρQ(V₂ − V₁),

for steady one-dimensional flow, with forces and velocities evaluated consistently in the selected direction. Pressure forces, weight components, boundary shear, and external reactions must be included.

Important Dimensionless Numbers

The Froude number compares inertia with gravity:

Fr = V/√(gD),

where D = A/T is hydraulic depth and T is top width. Flow is subcritical when Fr < 1, critical when Fr = 1, and supercritical when Fr > 1. In a rectangular channel, D equals the flow depth y.

The Reynolds number compares inertia with viscosity:

Re = VL/ν,

where L is a characteristic length and ν is kinematic viscosity. Low Reynolds numbers indicate stronger viscous influence; most full-scale rivers and hydraulic structures are turbulent. The Weber number, We = ρV²L/σ, compares inertia with surface tension and becomes important in small-scale models, nappe formation, and air–water interfaces.

2. Open-Channel Flow

Open-channel flow has a free surface exposed to atmospheric pressure. Its behaviour is governed by channel geometry, bed slope, roughness, discharge, and gravity. The hydraulic radius is

R = A/P,

where P is wetted perimeter. For uniform turbulent flow, the Manning equation is

Q = (1/n)AR^(2/3)S_f^(1/2),

where n is Manning’s roughness coefficient and S_f is the friction slope. In uniform flow, S_f is equal to the bed slope S₀. The equation is empirical and should not be applied uncritically outside its calibration range or to strongly non-uniform regions.

The specific energy relative to the channel bed is

E = y + αV²/(2g) = y + αQ²/(2gA²).

For a specified discharge and channel geometry, specific energy has a minimum at critical flow. The general critical-flow condition is

Q²T/(gA³) = 1,

when α is treated as approximately constant. Equivalently, Fr = 1. In a rectangular channel of width b, q = Q/b and

y_c = (q²/g)^(1/3).

The specific force, or momentum function, is useful in rapidly varied flow:

M = Q²/(gA) + A\bar{y},

for a prismatic channel, where \bar{y} is the depth of the centroid of the flow area below the free surface. The complete expression must be adapted when pressure distribution, channel slope, or non-prismatic geometry is significant.

Uniform and Non-Uniform Flow

Uniform flow has constant depth and velocity along the channel. It is an idealization approached in sufficiently long, prismatic reaches away from controls. Non-uniform flow occurs when depth changes because of changes in geometry, slope, roughness, discharge, or downstream control.

The normal depth y_n is the depth corresponding to uniform flow at a specified Q, geometry, roughness, and slope. The critical depth y_c is determined by the minimum-specific-energy condition. Their relative values classify channel slopes. For example, a mild slope generally has y_n > y_c, whereas a steep slope generally has y_n < y_c. This classification supports, but does not replace, profile analysis.

3. Gradually Varied Flow

Gradually varied flow (GVF) changes depth over a distance sufficiently long that the pressure distribution is approximately hydrostatic and the flow remains nearly one-dimensional. The standard differential equation for a prismatic channel is

dy/dx = (S₀ − S_f)/(1 − Fr²).

The sign convention depends on the direction selected for x; the equation above assumes x increases in the direction of flow. S₀ is bed slope, S_f is friction slope, and Fr is the local Froude number.

The equation becomes singular near critical flow, indicating that the assumptions or numerical treatment require special care. Profiles are conventionally classified into zones according to the relative positions of y, y_n, and y_c. Backwater profiles occur where depth is greater than normal depth and are often caused by downstream controls such as weirs, gates, or reservoirs. Drawdown profiles occur where depth decreases toward a control.

The direct-step method uses the energy equation between two sections. For a prismatic channel with no external work,

Δx = (E₂ − E₁)/(S₀ − \bar{S_f}),

where \bar{S_f} is the average friction slope over the step. The method is useful for hand checks and preliminary design. Standard-step methods are more flexible for irregular cross-sections and are widely implemented in one-dimensional software. Numerical solvers must manage boundary conditions, hydraulic jumps, wetting and drying, and transitions between subcritical and supercritical regimes.

Boundary Conditions and Controls

Subcritical flow is influenced primarily by downstream conditions, while supercritical flow is influenced primarily by upstream conditions. A critical section can act as a control, allowing discharge or depth to be related through critical-flow equations. Incorrect boundary-condition selection is a common source of apparently plausible but physically incorrect model results.

4. Rapidly Varied Flow and Hydraulic Jumps

Rapidly varied flow changes over a short distance, so hydrostatic pressure and gradual-variation assumptions may fail. Examples include flow beneath gates, over sharp transitions, at spillway toes, and through hydraulic jumps. Momentum conservation is generally more reliable than energy conservation across a jump because substantial energy is dissipated.

For a hydraulic jump in a horizontal rectangular channel, the sequent-depth relation is

y₂/y₁ = \frac{1}{2}[\sqrt{1 + 8Fr₁²} − 1],

where y₁ is the incoming depth, y₂ is the conjugate downstream depth, and Fr₁ = V₁/√(gy₁). The relation assumes a rectangular channel, steady flow, approximately hydrostatic pressure at the sections, negligible bed-slope effect across the jump, and a fully developed jump away from the roller boundaries.

The energy loss across the jump is

ΔE = (y₂ − y₁)³/(4y₁y₂)

for a rectangular channel. The downstream depth is substantially greater than the upstream depth when the incoming Froude number is large. Jump classification—undular, weak, oscillating, steady, or strong—depends on the incoming Froude number and affects basin performance, turbulence, spray, and residual energy.

Stilling-Basin Design Considerations

A stilling basin must accommodate the expected range of discharge and tailwater level. Design considerations include:

• matching the required conjugate depth with available tailwater;
• providing sufficient basin length for roller development;
• using chute blocks, baffle blocks, or an end sill where appropriate;
• checking uplift, scour, abrasion, vibration, and structural loads; and
• considering air entrainment and operation outside the design discharge.

A jump that is forced too far downstream may cause scour, whereas excessive tailwater may drown the jump and produce unstable recirculation. Physical modelling is often justified where three-dimensional effects or a wide operating range is important.

5. Pressurized Flow

Pressurized flow occurs when a conduit is full and the pressure is not everywhere atmospheric. For steady flow in a pipe, the Darcy–Weisbach head loss is

h_f = f(L/D)(V²/2g),

where f is the Darcy friction factor, L is pipe length, and D is internal diameter. Local losses are represented by

h_m = K(V²/2g),

where K is a dimensionless loss coefficient for an entrance, bend, valve, contraction, expansion, or other fitting. Total loss is the sum of distributed and local losses, subject to the applicability of the selected coefficients.

The friction factor depends mainly on Reynolds number and relative roughness ε/D. For laminar flow, f = 64/Re. For turbulent flow, it is obtained from the Colebrook–White relation or an explicit approximation such as the Haaland equation. Because conventions differ regarding Darcy and Fanning friction factors, the selected definition must be stated.

A hydraulic grade line is z + p/(ρg), while the energy grade line adds αV²/(2g). Their separation is the velocity head. These lines are valuable for identifying negative pressures, inadequate clearances, air-valve requirements, and possible cavitation zones.

Transients and Water Hammer

Rapid valve closure or pump failure can create pressure waves. A simplified Joukowsky relation is

Δp = ρaΔV,

where a is wave speed and ΔV is the change in velocity. The relation assumes a rapid change relative to the pressure-wave travel time and a suitably defined effective wave speed. Real systems require consideration of pipe elasticity, fluid compressibility, air pockets, friction, valve characteristics, surge tanks, and wave reflections. Transient analysis is therefore normally performed using the method of characteristics or a validated specialized model.

6. Cavitation and Aeration

Cavitation occurs when local absolute pressure falls to, or below, the vapour pressure of the liquid, allowing vapour cavities to form. Subsequent collapse can cause pitting, vibration, noise, loss of performance, and structural damage. Cavitation is distinct from ordinary air entrainment, although air can modify cavity behaviour.

A commonly used cavitation index is

σ = (p_abs − p_v)/(ρV²/2),

where p_abs is local absolute pressure, p_v is vapour pressure, and V is a representative velocity. The index is dimensionless. Cavitation risk increases as σ decreases, but the critical value depends on surface roughness, geometry, turbulence, dissolved air, scale, and duration. Consequently, a universal critical σ should not be assumed.

For spillways and high-speed conduits, designers should check pressure distributions, surface finish, offsets, joints, curvature, aeration devices, and operating envelopes. Aerators introduce air into low-pressure regions and can reduce damage, but they also create spray, vibration, noise, and possible downstream operational problems. Absolute pressure must be used; substituting gauge pressure in the cavitation criterion is a serious analytical error.

7. Dimensional Analysis and Similitude

Dimensional analysis reduces the number of independent variables and helps identify governing similarity requirements. The Buckingham π theorem states that a physically meaningful relationship involving n dimensional variables expressed in terms of k independent fundamental dimensions can be rewritten using n − k independent dimensionless groups.

For a spillway discharge, one possible functional relation is

Q = f(H, g, ρ, μ, L, σ),

where H is head, μ is dynamic viscosity, L is characteristic length, and σ is surface tension. A dimensionless form may include discharge coefficient, Froude number, Reynolds number, and Weber number. The precise groups depend on the chosen variables and physical regime.

Similarity Types

Geometric similarity requires corresponding dimensions to have a constant scale ratio. Kinematic similarity requires corresponding velocity and time ratios to be consistent. Dynamic similarity requires ratios of corresponding forces to be equal.

For free-surface flows dominated by gravity, Froude similarity is usually most important:

Fr_m = Fr_p.

If the length scale ratio is λ_L = L_m/L_p, then under Froude scaling, velocity ratio is λ_V = √λ_L and discharge ratio is

λ_Q = λ_L^(5/2)

for geometrically similar models. Time scales as λ_t = √λ_L. Subscripts m and p denote model and prototype.

Exact simultaneous equality of Froude and Reynolds numbers is generally impossible with the same fluid and a reduced geometric scale. This is a scale-effect problem. For a spillway model, maintaining adequate model Reynolds number may be more important than exact Reynolds similarity, while surface-tension effects must be checked near nappe edges and small openings. Model results should therefore be interpreted using distortion, sensitivity, and prototype validation—not copied mechanically.

8. Physical Hydraulic Models

Physical models are valuable when analytical solutions are inadequate, geometry is complex, or consequences of failure are severe. They can investigate flow distribution, free-surface profiles, hydraulic jumps, scour tendencies, pressure fields, air entrainment, sediment movement, and gate operation.

A sound modelling programme begins with a clearly defined decision question. The designer then selects the relevant dimensionless criteria, scale, fluid, measurement system, boundary conditions, and acceptance criteria. Instrumentation may include point gauges, pressure transducers, electromagnetic or ultrasonic velocity meters, acoustic Doppler systems, imaging, and discharge measurement devices.

Model verification confirms that the apparatus and instruments reproduce known or benchmark conditions. Validation assesses whether the model represents the prototype sufficiently for the intended decision. Repeatability, uncertainty, calibration drift, blockage, boundary roughness, model distortion, and scale effects should be documented.

Physical models do not automatically provide truth. A model can produce precise measurements of the wrong phenomenon if its similarity requirements are inappropriate. In particular, small models may exaggerate surface tension, viscous resistance, air entrainment differences, and sediment-scale effects.

9. Numerical Modelling

Numerical models solve governing equations using discretized spatial and temporal domains. One-dimensional models simplify cross-sections and are efficient for long river systems and network studies. Two-dimensional depth-averaged models resolve horizontal variation and are useful for floodplains, bends, bifurcations, and complex planform flow. Three-dimensional models resolve vertical structure and are justified where secondary currents, buoyancy, detailed turbulence, intakes, sediment suspension, or strongly three-dimensional recirculation are central.

The shallow-water equations are commonly used for depth-averaged free-surface flow. In conservative form, the continuity equation is

∂h/∂t + ∂(hu)/∂x + ∂(hv)/∂y = 0,

where h is water depth and u and v are depth-averaged velocity components. Momentum equations add flux gradients, bed slope, friction, and other source terms. Their validity requires hydrostatic pressure and a horizontal length scale sufficiently large relative to flow depth. They are not automatically appropriate for highly vertical jets, strongly dispersive waves, or detailed air–water turbulence.

The Navier–Stokes equations provide a more general basis, but turbulence closure, free-surface tracking, computational cost, and uncertain boundary conditions remain significant challenges. Reynolds-averaged, large-eddy, and direct numerical simulation represent increasing levels of turbulence resolution and computational demand. No turbulence model is universally superior.

Numerical Quality Assurance

A credible numerical study includes:

1. conceptual-model definition and stated assumptions;

2. mesh or grid-convergence testing;

3. time-step sensitivity assessment where unsteady flow is modelled;

4. conservation and stability checks;

5. calibration using independent field or laboratory data;

6. validation against observations not used for calibration; and

7. sensitivity and uncertainty analysis.

Calibration should not compensate blindly for incorrect geometry or missing physics. Over-calibration can produce excellent agreement at one location while reducing predictive skill elsewhere. Results should be reported with ranges and limitations, not excessive precision.

10. Integrated Applied Case: Spillway and Stilling-Basin Assessment

The following is an illustrative case based on a hypothetical 1,200 m³/s spillway rehabilitation project.

Background and Problem

A concrete gravity dam requires increased flood-release capacity. Preliminary calculations indicate a high-velocity chute discharging into a stilling basin. The owner is concerned about downstream scour, low pressures near a chute transition, and performance during both design and partial-gate releases.

Analysis

A one-dimensional energy analysis establishes approximate water-surface elevations and velocity heads. The chute velocity is then used to estimate Reynolds and Froude numbers. Because the flow is supercritical and the transition is geometrically complex, the one-dimensional results are treated as screening calculations rather than final evidence.

A hydraulic-jump calculation is made for each relevant discharge using the incoming depth and Froude number. Tailwater-rating uncertainty is included because the basin jump location depends strongly on tailwater. A cavitation assessment uses absolute pressure measurements from a physical model and computes the cavitation index at the lowest-pressure locations.

Approach and Solution

A geometrically similar physical model is constructed using Froude similarity. The model is sufficiently large to reduce viscous and surface-tension scale effects, and pressure taps are installed at the chute transition and basin floor. Tests cover the design flood, lower discharges, asymmetric gate operation, and high tailwater.

A two-dimensional numerical model is independently developed for downstream flow distribution and scour-risk screening. It is calibrated against measured water levels and velocities, then validated at a discharge not used in calibration. The model is not used to replace the physical model in the chute and aeration region, where three-dimensional free-surface effects are dominant.

Results show that the initial basin configuration produces an unstable jump at partial discharge and local low pressures near a geometric offset. The design is modified by smoothing the offset, adding an aeration device, and revising the end-sill geometry. The physical model demonstrates improved pressure margins and a more stable jump over the operating range. The numerical model indicates reduced concentration of downstream velocity, but residual scour risk remains and is addressed through bed protection and monitoring.

Outcome and Learning Points

The final decision is based on agreement among analytical checks, model observations, numerical simulations, and structural and operational constraints. The case illustrates model hierarchy: simple equations provide screening, physical models resolve complex local phenomena, and numerical models extend understanding over larger domains and operating conditions. It also demonstrates that uncertainty in tailwater, gate operation, roughness, and scale effects must be carried into the decision rather than hidden in a single design value.

11. Limitations, Governance, and Contemporary Practice

Hydraulic modelling is subject to epistemic uncertainty from incomplete knowledge and aleatory variability from changing flow conditions. Important sources include uncertain roughness, sediment transport, debris, air entrainment, climate-adjusted floods, bathymetric change, and operational decisions.

Professional practice requires traceable data, documented assumptions, independent checking, version control, and clear communication of uncertainty. Safety-critical hydraulic decisions should not rely on an opaque software output or an unverified artificial-intelligence system. AI and machine-learning methods can assist with surrogate modelling, image-based flow measurement, anomaly detection, parameter estimation, and optimization. However, training-data bias, extrapolation failure, conservation violations, and lack of interpretability require strict physical validation. AI should support—not replace—governing equations, engineering judgement, and accountable review.

Best practice is to use a staged workflow: define the decision, identify dominant physics, perform hand calculations, select the simplest adequate model, verify implementation, validate against independent evidence, quantify uncertainty, and monitor prototype performance after construction. Model complexity should be justified by the decision risk and the physics that must be resolved.

Chapter Summary

Hydraulic analysis combines conservation of mass, energy, and momentum with empirical resistance laws and dimensionless reasoning. Open-channel flow is characterized by free-surface effects and the Froude number, while pressurized systems require pressure-head, friction-loss, and transient analysis. Uniform flow provides normal depth; critical flow corresponds to minimum specific energy; gradually varied flow is described by the gradually varied flow equation; and rapidly varied flow, especially hydraulic jumps, is analysed primarily through momentum conservation.

Cavitation assessment requires absolute pressure and vapour pressure, with the cavitation index interpreted in relation to geometry and scale. Dimensional analysis identifies controlling dimensionless groups, while similitude determines which force ratios should be preserved in a physical model. Froude similarity is normally dominant in free-surface hydraulic structures, but Reynolds, Weber, and other scale effects must be assessed.

Physical and numerical models are complementary. Physical models are powerful for complex local free-surface and air–water processes; numerical models efficiently investigate spatial alternatives and operating conditions. Both require verification, validation, uncertainty assessment, and transparent reporting. Sound hydraulic engineering uses a hierarchy of models and never treats computational precision as a substitute for physical validity or professional responsibility.

Key Terms

Open-channel flow: Flow with a free surface exposed to atmospheric pressure.

Pressurized flow: Flow in a closed or full conduit where pressure is not everywhere atmospheric.

Hydraulic radius: Flow area divided by wetted perimeter.

Froude number: Ratio expressing the relative importance of inertia and gravity.

Reynolds number: Ratio expressing the relative importance of inertia and viscosity.

Specific energy: Energy per unit weight measured relative to the channel bed.

Critical flow: Flow at minimum specific energy for a given discharge.

Normal depth: Depth associated with uniform flow under specified conditions.

Gradually varied flow: Flow whose depth changes slowly enough for hydrostatic pressure to remain a reasonable approximation.

Rapidly varied flow: Flow whose depth changes over a short distance, often with non-hydrostatic effects.

Hydraulic jump: Abrupt transition from supercritical to subcritical flow with energy dissipation.

Sequent depths: Upstream and downstream conjugate depths of a hydraulic jump.

Hydraulic grade line: Elevation plus pressure head in a pressurized system.

Cavitation: Formation and collapse of vapour cavities caused by low absolute pressure.

Cavitation index: Dimensionless measure comparing pressure margin above vapour pressure with velocity head.

Dimensional analysis: Method of expressing physical relations using dimensionless variables.

Similitude: Correspondence between model and prototype geometry, motion, and forces.

Scale effect: Difference caused when a model cannot reproduce all relevant similarity conditions.

Verification: Demonstration that a model or implementation solves its intended equations correctly.

Validation: Assessment of whether a model adequately represents the real system for its intended use.

Review Questions

Short Answer Questions

1. Distinguish between normal depth and critical depth.

2. Why is momentum conservation preferred to energy conservation across a hydraulic jump?

3. Define the Froude number and explain its significance in open-channel flow.

4. What is the difference between the hydraulic grade line and the energy grade line?

5. Why must absolute pressure be used in a cavitation assessment?

Long Answer Questions

1. Derive and explain the critical-flow condition for a general open channel.

2. Discuss the assumptions, interpretation, and limitations of the gradually varied flow equation.

3. Explain hydraulic-jump theory and the implications for stilling-basin design.

4. Compare Froude and Reynolds similarity in reduced-scale hydraulic models.

5. Explain a rigorous verification, calibration, and validation process for a two-dimensional numerical hydraulic model.

Analytical / Application Questions

1. A rectangular channel is 6 m wide and carries 30 m³/s. Determine the critical depth using the rectangular-channel critical-depth equation and comment on the flow regime if the actual depth is 1.8 m.

2. A hydraulic jump enters a rectangular channel at depth 0.40 m with Froude number 5. Determine the sequent depth and energy loss using the hydraulic-jump relations. State the assumptions.

3. A spillway model is built at a geometric length scale of 1:25 under Froude similarity. Determine the velocity, discharge, and time scale ratios and discuss two likely scale effects that require checking.

Multiple Choice Questions with Answers

Question 1

Which dimensionless number primarily compares inertia with gravity in open-channel flow?

A. Reynolds number

B. Froude number

C. Weber number

D. Mach number

Correct Answer: B. Froude number

Question 2

For a rectangular channel, the critical depth is given by which expression?

A. (Q²/g)^(1/2)

B. (Q/g)^(1/3)

C. (q²/g)^(1/3)

D. gq²

Correct Answer: C. (q²/g)^(1/3)

Question 3

Which conservation principle is generally preferred for analysing a hydraulic jump?

A. Mass only

B. Momentum

C. Temperature

D. Viscosity

Correct Answer: B. Momentum

Question 4

In a hydraulic jump, the flow normally changes from:

A. Subcritical to supercritical

B. Uniform to pressurized

C. Supercritical to subcritical

D. Laminar to inviscid

Correct Answer: C. Supercritical to subcritical

Question 5

The Darcy–Weisbach equation calculates:

A. Critical depth

B. Pipe friction head loss

C. Hydraulic-jump length

D. Cavitation pressure

Correct Answer: B. Pipe friction head loss

Question 6

Cavitation is most directly associated with local pressure falling toward:

A. Atmospheric pressure only

B. The pipe roughness pressure

C. Vapour pressure

D. The stagnation pressure

Correct Answer: C. Vapour pressure

Question 7

Under Froude similarity, if the model length scale is 1:25, the model-to-prototype velocity ratio is:

A. 1/5

B. 1/25

C. 1/125

D. 5

Correct Answer: A. 1/5

Question 8

Which model type generally resolves horizontal variation while using depth-averaged equations?

A. One-dimensional model

B. Two-dimensional shallow-water model

C. Zero-dimensional model

D. Static structural model

Correct Answer: B. Two-dimensional shallow-water model

Question 9

The hydraulic grade line in a pressurized pipe is:

A. z + p/(ρg)

B. z + V²/(2g)

C. pV/ρ

D. z − p/(ρg)

Correct Answer: A. z + p/(ρg)

Question 10

A numerical model is validated when:

A. Its mesh contains many cells regardless of results

B. It reproduces independent observations adequately for its intended use

C. It uses the most complex turbulence model available

D. It has no calibration parameters

Correct Answer: B. It reproduces independent observations adequately for its intended use