Fluid Mechanics & Hydraulic Machines

Chapter 8
Open-Channel Flow and Hydraulic Structures

Bachelors Program in Mechnical Engineering
Fluid Mechanics & Hydraulic Machines   •   Bachelor

Open-Channel Flow and Hydraulic Structures

Introduction and Context

Open-channel flow is the movement of liquid with a free surface exposed to the atmosphere. Rivers, canals, drains, spillways, irrigation channels, and partially filled conduits are common examples. Unlike pressurized pipe flow, the driving force is primarily the component of gravity along the channel, while the free surface adjusts to satisfy the hydraulic and geometric conditions.

Analysis of open channels is central to irrigation, flood protection, hydropower, urban drainage, river engineering, and environmental restoration. The engineer must determine discharge, depth, velocity, energy losses, flow regime, and the effect of structures such as weirs, gates, transitions, and stilling basins.

Learning Objectives

By the end of this chapter, the reader should be able to:

  • distinguish uniform, gradually varied, rapidly varied, subcritical, and supercritical flow;
  • calculate hydraulic radius using channel geometry;
  • apply the Chezy and Manning equations to uniform flow;
  • use specific energy to determine critical depth and assess channel transitions;
  • analyze hydraulic jumps and estimate energy dissipation;
  • calculate discharge over common weirs and notches;
  • explain the hydraulic function of spillways, stilling basins, and control structures;
  • recognize assumptions, limitations, safety requirements, and field uncertainties in open-channel design.

1. Fundamental Geometry and Terminology

For a channel cross-section, the principal geometric quantities are:

SymbolMeaningTypical unit
AFlow aream²
PWetted perimeter, excluding the free surfacem
RHydraulic radius, A ÷ Pm
TTop width of the free surfacem
yFlow depth measured verticallym
QDischargem³/s
VMean velocity, Q ÷ Am/s
S₀Bed slopedimensionless
S_fFriction or energy slopedimensionless

The hydraulic radius is

R = A ÷ P

It is a measure of the efficiency with which a channel conveys flow. For a rectangular channel of width b and depth y,

A = by, P = b + 2y, R = by ÷ (b + 2y)

For a very wide rectangular channel, where b is much greater than y, P ≈ b and therefore R ≈ y. This approximation should not be used for narrow or deep channels.

The hydraulic depth is

D = A ÷ T

Hydraulic depth is important in wave propagation and critical-flow analysis. In a rectangular channel, D = y.

Common Channel Sections

For a trapezoidal channel with bottom width b, side slope z horizontal to 1 vertical, and depth y,

A = y(b + zy)
P = b + 2y√(1 + z²)

For a triangular channel, b = 0, so A = zy² and P = 2y√(1 + z²).

2. Classification of Open-Channel Flow

Uniform and Nonuniform Flow

Uniform flow has constant depth, area, velocity, and flow resistance along the channel. It is an idealization generally applicable to long, prismatic channels away from entrances, controls, bends, and transitions. In uniform flow,

S₀ = Sf

Nonuniform flow occurs when depth, velocity, or area changes along the channel. It may be gradually varied or rapidly varied.

Gradually Varied Flow

Gradually varied flow changes depth over a relatively long distance, so hydrostatic pressure distribution remains a reasonable approximation. Examples include backwater upstream of a dam and drawdown approaching a control section. The governing one-dimensional equation is

dy ÷ dx = (S₀ − Sf) ÷ (1 − Fr²)

where x is the longitudinal coordinate and Fr is the Froude number. This equation assumes steady flow, a prismatic or known channel geometry, hydrostatic pressure, and negligible local acceleration effects.

Rapidly Varied Flow

Rapidly varied flow changes substantially over a short distance. Pressure is not necessarily hydrostatic, and local geometry strongly affects the flow. Hydraulic jumps, drops, gates, and abrupt transitions are examples. Momentum analysis is usually more reliable than an energy equation because substantial turbulence dissipates mechanical energy.

Froude Number and Flow Regime

The Froude number compares inertial effects with gravity-wave effects:

Fr = V ÷ √(gD)

where g is gravitational acceleration, approximately 9.81 m/s², and D is hydraulic depth.

  • Fr < 1: subcritical flow; gravity effects dominate, and disturbances can travel upstream.
  • Fr = 1: critical flow.
  • Fr > 1: supercritical flow; inertia dominates, and disturbances cannot propagate upstream.

The classification is essential when deciding whether an upstream or downstream boundary condition controls the water surface.

3. Uniform Flow: Chezy and Manning Equations

For steady uniform flow in a prismatic channel, the downslope component of the water weight is balanced by boundary resistance. The Chezy equation is

V = C√(R Sf)

and therefore

Q = AC√(R Sf)

where C is the Chezy coefficient, with units m^(1/2)/s. The coefficient depends on roughness, channel shape, Reynolds number, vegetation, sediment, and degree of irregularity.

The Manning equation is more widely used in engineering practice:

V = (1 ÷ n)R^(2/3)Sf^(1/2)
Q = (1 ÷ n)AR^(2/3)Sf^(1/2)

where n is Manning’s roughness coefficient, with units s/m^(1/3). For uniform flow in a prismatic channel, Sf = S₀.

The equations are empirical and presume steady flow, approximately uniform conditions, one-dimensional mean velocity, and a known roughness. They should not be blindly applied near hydraulic structures, sharp bends, rapidly changing sections, or heavily nonuniform vegetation.

Worked Example 1: Uniform Flow in a Rectangular Channel

A concrete rectangular channel is 3.0 m wide and carries water at a depth of 1.0 m. The bed slope is 0.0016 and Manning’s n is 0.015. Determine the discharge and mean velocity, assuming uniform flow.

The area and wetted perimeter are

A = 3.0 × 1.0 = 3.0 m²

P = 3.0 + 2(1.0) = 5.0 m

Thus,

R = 3.0 ÷ 5.0 = 0.60 m

Using Manning’s equation,

V = (1 ÷ 0.015)(0.60)^(2/3)(0.0016)^(1/2)

V ≈ 2.01 m/s

Therefore,

Q = AV = 3.0 × 2.01 ≈ 6.03 m³/s

The result is appropriate only if the channel is sufficiently long and unobstructed for uniform-flow conditions to develop.

4. Specific Energy and Critical Flow

Specific energy is the energy per unit weight measured relative to the channel bed. For a channel with a uniform velocity distribution approximation,

E = y + αV² ÷ (2g)

where α is the kinetic-energy correction coefficient. In many elementary engineering calculations α is taken as 1. Since V = Q ÷ A,

E = y + αQ² ÷ (2gA²)

For a fixed discharge, specific energy varies with depth. The minimum specific energy occurs at critical flow. The general critical-flow condition is

Q²T ÷ (gA³) = 1

when α = 1. Equivalently,

Fr = 1

For a rectangular channel of width b, A = by and T = b. The critical depth is

yc = (q² ÷ g)^(1/3)

where q = Q ÷ b is discharge per unit width, in m²/s. The minimum specific energy is

Ec = 3yc ÷ 2

for a rectangular channel with α = 1.

A channel transition may convert subcritical flow to supercritical flow or vice versa. A smooth contraction can reduce depth and increase velocity. If the available specific energy is insufficient to pass the discharge, upstream water rises, producing a control or backwater effect. A transition should be designed to avoid separation, excessive losses, cavitation, and undesirable flow regimes.

Worked Example 2: Critical Depth

A rectangular channel 4 m wide carries 12 m³/s. Find the critical depth and minimum specific energy.

Discharge per unit width is

q = 12 ÷ 4 = 3.0 m²/s

Hence,

yc = (3.0² ÷ 9.81)^(1/3) ≈ 0.972 m

and

Ec = 3(0.972) ÷ 2 ≈ 1.458 m

At critical flow, the velocity is approximately V = √(gyc) ≈ 3.09 m/s, and Fr ≈ 1.

5. Gradually Varied Flow and Water-Surface Profiles

The gradually varied flow equation describes the change in depth along a channel. Its interpretation depends on bed slope, normal depth yₙ, and critical depth y_c. Typical profiles include backwater curves upstream of a dam, drawdown curves near a steep drop, and profiles behind sluice gates.

A practical solution may use the standard-step method. Between two cross-sections, the energy equation is written as

z₁ + y₁ + αV₁² ÷ (2g) = z₂ + y₂ + αV₂² ÷ (2g) + hf

where z is bed elevation and h_f is friction loss. For a reach of length L, a common approximation is

hf ≈ L(Sf,1 + Sf,2) ÷ 2

The method is iterative because area, velocity, hydraulic radius, and friction slope depend on the unknown depth. Computational packages often automate this calculation, but the engineer must still select valid boundary conditions and inspect the resulting profile.

6. Hydraulic Jumps and Energy Dissipation

A hydraulic jump is a rapidly varied transition from supercritical to subcritical flow. It is accompanied by a sudden increase in depth, intense turbulence, air entrainment, and energy dissipation. Hydraulic jumps are deliberately created in stilling basins downstream of spillways and outlets to protect the channel bed.

For a rectangular channel, conservation of momentum gives the conjugate-depth relationship:

y₂ ÷ y₁ = 1/2[√(1 + 8Fr₁²) − 1]

where y₁ is the supercritical incoming depth, y₂ is the downstream conjugate depth, and Fr₁ is the upstream Froude number based on y₁. The equation assumes a horizontal rectangular channel, hydrostatic pressure at sections away from the roller, negligible external forces, and approximately uniform velocity distribution.

The energy loss across the jump is

∆E = (y₂ − y₁)³ ÷ (4y₁y₂)

for a rectangular channel. This loss is mechanical energy converted mainly into turbulence and heat; it is not a violation of conservation of energy.

Jump classification is often based on upstream Froude number. Weak jumps occur for values only slightly above 1, while stronger jumps occur at larger values. A stable stilling basin requires compatibility between the incoming flow, tailwater depth, basin geometry, and energy-dissipation devices.

Worked Example 3: Hydraulic Jump

Water enters a rectangular stilling basin at y₁ = 0.20 m and V₁ = 8.0 m/s. Determine the upstream Froude number, conjugate depth, and energy loss.

For a rectangular channel, D = y₁. Therefore,

Fr₁ = 8.0 ÷ √(9.81 × 0.20) ≈ 5.72

The conjugate-depth ratio is

y₂ ÷ y₁ = 1/2[√(1 + 8(5.72)²) − 1] ≈ 7.59

Thus,

y₂ ≈ 7.59 × 0.20 = 1.52 m

The energy loss is

∆E = (1.52 − 0.20)³ ÷ [4(0.20)(1.52)] ≈ 1.88 m

The large loss confirms that a powerful jump is present and that substantial structural protection is required.

7. Weirs and Notches

A weir is a barrier over which water flows, whereas a notch is an opening cut into a thin plate or wall. These structures measure discharge, regulate water level, or provide overflow protection. Their equations derive from the energy equation and require discharge coefficients determined experimentally.

For a sharp-crested rectangular weir with width b and head H above the crest, a common model is

Q = (2 ÷ 3)Cd b√(2g)H^(3/2)

where C_d is the discharge coefficient. The formula assumes a sharp crest, free aerated nappe, adequate upstream approach, and head measured sufficiently upstream to avoid local drawdown effects.

For a 90° V-notch, the idealized expression is

Q = (8 ÷ 15)Cd√(2g)tan(θ ÷ 2)H^(5/2)

where θ is the notch angle. V-notches are sensitive to small discharges and are often useful in laboratory and small-stream measurements.

A broad-crested weir has a crest long enough for nearly parallel flow over the crest. Its discharge relation is often written in the form

Q = C b H^(3/2)

where C incorporates gravity and the relevant coefficient. Because coefficients depend strongly on crest geometry and approach conditions, calibrated or standard-design values should be used.

Important practical corrections include approach-velocity head, end contractions, submergence, crest condition, sediment deposition, and nappe aeration. A weir that becomes submerged no longer follows the free-discharge equation without a submergence correction.

8. Spillways, Gates, and Channel Transitions

A spillway safely conveys excess reservoir water downstream. Major types include ogee, broad-crested, chute, side-channel, shaft, and labyrinth spillways. Selection depends on discharge, reservoir geometry, available elevation difference, foundation conditions, sediment, debris, and downstream energy dissipation.

An ogee spillway is shaped to approximate the lower nappe of a free overfall at its design head. At heads different from the design value, separation or negative pressure may develop. Chute spillways convey high-velocity flow down a steep slope and normally terminate in a hydraulic jump basin, flip bucket, or plunge pool.

Sluice gates regulate discharge and commonly produce supercritical jets downstream. Gate openings, contraction, ventilation, cavitation resistance, and downstream tailwater must be considered. A gate should never be designed solely from a local velocity calculation; structural vibration and transient loads may govern.

Channel transitions connect sections of different width, bed elevation, or alignment. A gradual transition reduces separation and local losses. In subcritical flow, downstream conditions commonly influence the upstream water surface; in supercritical flow, upstream conditions dominate. If a transition drives the flow toward critical conditions, it acts as a control section.

9. Surges and Unsteady Effects

A surge is a moving change in water depth caused by rapid gate movement, sudden discharge variation, obstruction, or flood-wave propagation. A positive surge raises the water surface; a negative surge lowers it. Rapid operational changes can create dangerous waves in canals and conduits.

Surge analysis requires the continuity and momentum equations for unsteady open-channel flow. In simplified form, the Saint-Venant equations are

∂A ÷ ∂t + ∂Q ÷ ∂x = 0

and

∂Q ÷ ∂t + ∂(Q² ÷ A) ÷ ∂x + gA∂y ÷ ∂x = gA(S₀ − Sf)

where t is time. These equations assume one-dimensional flow and require suitable initial and boundary conditions. Numerical models are generally used for real systems. Gate operation should be gradual enough to limit surge amplitude, while emergency controls and freeboard must accommodate credible transients.

10. Applied Case: Illustrative Irrigation Canal Upgrade

Background

An illustrative irrigation scheme uses a 4 m wide concrete-lined canal to deliver water to farms. The design discharge is 10 m³/s, but operators report overtopping during gate adjustments and erosion downstream of a control gate.

Problem

Field observations show that the canal operates near critical depth during peak delivery. A sudden gate opening creates a supercritical jet, followed by an unstable hydraulic jump whose toe moves through the unprotected reach.

Analysis

The discharge per unit width is q = 10 ÷ 4 = 2.5 m²/s. The rectangular critical depth is

yc = (2.5² ÷ 9.81)^(1/3) ≈ 0.86 m

Operating depth near this value means that small changes in discharge or depth can cause large changes in flow regime. The existing downstream tailwater is insufficient to hold the jump inside the short basin. The gate movement also produces a positive surge that temporarily raises the upstream water level.

Approach and Solution

The proposed upgrade includes a longer stilling basin, a designed end sill, riprap and concrete cutoff protection, a water-level sensor, and a slower gate-opening protocol. A one-dimensional unsteady model is used to test several operation scenarios. The model is checked against field depth and discharge observations rather than accepted without calibration.

Outcome

The redesigned basin places the expected hydraulic jump within the protected area for the design discharge and a specified range of tailwater levels. The operating rule limits gate movement rate, reducing surge amplitude. Freeboard is increased to accommodate uncertainty in roughness, sediment, and emergency inflow.

Learning Points

The case demonstrates that open-channel design is not only a matter of inserting values into Manning’s equation. Critical-flow sensitivity, rapidly varied flow, transient operation, structural protection, monitoring, and uncertainty must be considered together.

11. Design Practice, Limitations, and Safety

Open-channel equations are simplified representations of complex three-dimensional flows. Major sources of uncertainty include roughness variation, weeds, sediment, debris, nonuniform cross-sections, secondary currents, air entrainment, changing tailwater, and inaccurate survey data.

Good practice includes:

  1. survey the channel and identify all hydraulic controls;
  2. define design discharges, flood frequencies, operating levels, and allowable velocities;
  3. calculate normal and critical depths before selecting a computational method;
  4. distinguish energy-based analysis from momentum-based analysis;
  5. check freeboard, erosion, uplift, seepage, cavitation, debris blockage, and structural loads;
  6. calibrate discharge coefficients and roughness values where field measurements are available;
  7. model unsteady operation when gates or rapidly changing inflows are important;
  8. inspect structures after floods and revise operating procedures when evidence changes.

Safety and regulatory compliance are essential. Spillways must safely pass the governing design flood, and failure consequences should be assessed. Water-level controls should include safe access, emergency shutdown provisions, and protection against unanticipated gate movement. Environmental design should consider fish passage, sediment continuity, downstream habitat, and ecological flow requirements. Automated sensors and AI-based forecasting can improve warning and operation, but they do not replace hydraulic verification, sensor maintenance, cybersecurity, or professional judgment.

Chapter Summary

Open-channel flow is governed by gravity, geometry, resistance, and free-surface conditions. Hydraulic radius connects flow geometry to resistance equations, while the Chezy and Manning equations estimate uniform-flow velocity and discharge under stated assumptions. Specific energy provides a powerful framework for critical depth, channel controls, and transitions. The Froude number distinguishes subcritical, critical, and supercritical regimes and indicates the direction of disturbance propagation.

Gradually varied flow is analyzed through the water-surface-profile equation or step methods, whereas rapidly varied flow is commonly analyzed with momentum principles. Hydraulic jumps convert supercritical flow to subcritical flow and dissipate energy, making them valuable in stilling basins but potentially destructive in unprotected channels. Weirs and notches measure or regulate discharge, although their coefficients and approach conditions must be carefully controlled. Spillways, gates, transitions, and surge-control systems must be designed as integrated hydraulic and structural systems. Reliable practice combines equations, field data, numerical modeling, safety checks, environmental responsibility, and professional judgment.

Key Terms

Hydraulic radius: Flow area divided by wetted perimeter.

Hydraulic depth: Flow area divided by free-surface width.

Uniform flow: Flow with constant depth and velocity along a channel.

Nonuniform flow: Flow whose depth or velocity varies spatially.

Gradually varied flow: Flow whose depth changes over a relatively long distance.

Rapidly varied flow: Flow with substantial depth change over a short distance.

Manning coefficient: Empirical roughness parameter n used in the Manning equation.

Chezy coefficient: Empirical resistance coefficient C in the Chezy equation.

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

Critical depth: Depth corresponding to minimum specific energy for a given discharge.

Froude number: Ratio of inertial effects to gravity-wave effects.

Subcritical flow: Flow with Froude number less than one.

Supercritical flow: Flow with Froude number greater than one.

Hydraulic jump: Rapid transition from supercritical to subcritical flow.

Conjugate depths: Two depths satisfying the momentum relationship across a hydraulic jump.

Weir: A barrier over which water flows.

Notch: An opening in a thin plate or wall used for measurement or control.

Spillway: Structure that safely conveys excess water from a reservoir or channel.

Surge: A moving water-level disturbance caused by unsteady flow.

Stilling basin: Protected structure designed to dissipate flow energy, often through a hydraulic jump.

Review Questions

Short Answer Questions

  1. Define hydraulic radius and explain why it is useful in open-channel-flow calculations.
  2. Distinguish between gradually varied and rapidly varied flow.
  3. What physical interpretation is given to a Froude number less than one?
  4. State the assumptions underlying the Manning equation for uniform flow.
  5. Why is a hydraulic jump useful downstream of a spillway?

Long Answer Questions

  1. Derive the rectangular-channel critical-depth equation from the minimum-specific-energy condition.
  2. Explain the differences between the Chezy and Manning equations, including coefficients, units, and limitations.
  3. Describe the standard-step approach for calculating a gradually varied water-surface profile.
  4. Explain the momentum analysis of a hydraulic jump and discuss the assumptions behind the conjugate-depth equation.
  5. Compare sharp-crested weirs, V-notches, and broad-crested weirs as discharge-measuring devices.

Analytical / Application Questions

  1. A rectangular channel 5 m wide carries 15 m³/s. Calculate its critical depth and minimum specific energy. State whether a measured depth of 0.60 m would represent subcritical or supercritical flow.
  2. A concrete channel has width 3 m, depth 0.9 m, slope 0.002, and Manning coefficient 0.014. Calculate its normal-flow discharge and assess whether the result is consistent with a design discharge of 8 m³/s.
  3. A spillway discharges a high-velocity jet into a downstream channel where erosion is occurring. Develop a hydraulic investigation plan covering critical depth, Froude number, hydraulic-jump location, tailwater, energy dissipation, transient operation, and field validation.

Multiple Choice Questions with Answers

Question 1

For a rectangular channel, the hydraulic radius is defined by which expression?

A. A ÷ T

B. A ÷ P

C. P ÷ A

D. Q ÷ A

Correct Answer: B. A ÷ P

Question 2

A flow with Froude number greater than one is classified as:

A. Subcritical

B. Critical

C. Supercritical

D. Uniform only

Correct Answer: C. Supercritical

Question 3

Which equation is commonly used to estimate uniform open-channel discharge using Manning roughness?

A. Q = (1 ÷ n)AR^(2/3)S_f^(1/2)

B. Q = AV²

C. Q = gA y

D. Q = C_d A ÷ g

Correct Answer: A. Q = (1 ÷ n)AR^(2/3)S_f^(1/2)

Question 4

For a rectangular channel, critical flow occurs when:

A. Fr = 0

B. Fr < 1

C. Fr = 1

D. Fr > 2 only

Correct Answer: C. Fr = 1

Question 5

A hydraulic jump changes the flow from:

A. Subcritical to supercritical

B. Supercritical to subcritical

C. Uniform to pressurized only

D. Laminar to stagnant

Correct Answer: B. Supercritical to subcritical

Question 6

A rectangular channel 4 m wide carries 12 m³/s. Its discharge per unit width is:

A. 0.33 m²/s

B. 3.0 m²/s

C. 8.0 m²/s

D. 48 m²/s

Correct Answer: B. 3.0 m²/s

Question 7

The primary purpose of a stilling basin downstream of a spillway is to:

A. Increase reservoir storage

B. Eliminate all friction

C. Dissipate excess flow energy safely

D. Convert water into vapor

Correct Answer: C. Dissipate excess flow energy safely

Question 8

The idealized discharge relation for a sharp-crested rectangular weir varies with head approximately as:

A. H^(1/2)

B. H

C. H^(3/2)

D. H^3

Correct Answer: C. H^(3/2)

Question 9

In gradually varied flow, the pressure distribution is generally approximated as:

A. Hydrostatic

B. Vacuum everywhere

C. Independent of depth

D. Atmospheric at every interior point

Correct Answer: A. Hydrostatic

Question 10

For uniform flow in a prismatic channel, the usual relationship between bed slope and friction slope is:

A. S₀ = 0

B. S₀ = S_f

C. S₀ = 1 ÷ S_f

D. S₀ = −S_f

Correct Answer: B. S₀ = S_f

Fluid Mechanics & Hydraulic Machines
Chapter 8 of 14