Open Channel Flow

An in-depth introduction to the principles governing the flow of water in open channels, a critical aspect of civil and environmental engineering.

Overview

Open channel flow is the flow of a liquid within a conduit with a free surface subjected to atmospheric pressure. This applies to natural rivers, engineered canals, partially full storm sewers, and spillways. The driving force is gravity, and the resisting force is boundary friction. Key concepts include Manning's Equation, Tractive Force Method, Specific Energy (and channel transitions), the Froude Number, and rapidly varied flows like the Hydraulic Jump.

  1. Uniform Flow and Manning's Equation

Uniform flow occurs when the depth of flow, water area, and velocity remain constant along the length of the channel. This implies that the gravitational forces pushing the water downhill are exactly balanced by the frictional forces from the channel bed and banks.
The universally accepted empirical formula to analyze uniform open channel flow is Manning's Equation (historically predated by the Chezy equation, v=CRSv = C \sqrt{RS}):

Formula

Mathematical expression.

v=1nR2/3S1/2(SI Units) v = \frac{1}{n} R^{2/3} S^{1/2} \quad \text{(SI Units)}

Variables

SymbolDescriptionUnit
vvAverage Velocitym/s
nnManning's Roughness Coefficientdimensionless
RRHydraulic Radius ($A/P$)m
SSChannel Bed Slopem/m

Formula

Mathematical expression.

Q=1nAR2/3S1/2(SI Units) Q = \frac{1}{n} A R^{2/3} S^{1/2} \quad \text{(SI Units)}

Variables

SymbolDescriptionUnit
QQFlow Rate (Discharge)m³/s
AACross-sectional Area of Flow
nnManning's Roughness Coefficientdimensionless
RRHydraulic Radiusm
SSChannel Bed Slopem/m
Where:
  • v: Mean velocity (m/s)
  • Q: Discharge or flow rate (m³/s)
  • n: Manning's roughness coefficient (dimensionless)
  • A: Cross-sectional area of flow (m²)
  • P: Wetted perimeter (m)
  • R: Hydraulic radius (m), defined as R=A/PR = A / P
  • S: Longitudinal slope of the channel bed (m/m)

Most Efficient Channel Section

For a given flow area (AA) and slope (SS), the most efficient section provides the maximum discharge (QQ). This occurs when the hydraulic radius (RR) is maximized, meaning the wetted perimeter (PP) must be minimized. For a rectangular channel, the most efficient section is when the width is exactly twice the depth (b=2yb = 2y).

  1. Tractive Force Method for Unlined Channels

When designing unlined, erodible earth canals, simply using Manning's equation is insufficient because high velocities will scour the channel.

Tractive Force (τ0\tau_0)

The average shear stress exerted by the flowing water on the channel boundary. The Tractive Force Method ensures that the actual shear stress acting on the bed and banks does not exceed the permissible shear stress of the specific soil material.

Formula

Mathematical expression.

τ0=γwRS \tau_0 = \gamma_w \cdot R \cdot S

Variables

SymbolDescriptionUnit
τ0\tau_0Average Boundary Shear StressN/m²
γw\gamma_wSpecific Weight of WaterN/m³
RRHydraulic Radiusm
SSChannel Bed Slopem/m
Where γw\gamma_w is the specific weight of water. The channel dimensions must be sized so that τ0τpermissible\tau_0 \le \tau_{permissible}.

  1. Specific Energy and Channel Transitions

Specific energy (EE) in an open channel is defined as the total energy head of the water measured relative to the channel bottom. It is the sum of the flow depth (yy) and the velocity head.

Formula

Mathematical expression.

E=y+v22g=y+Q22gA2 E = y + \frac{v^2}{2g} = y + \frac{Q^2}{2gA^2}

Variables

SymbolDescriptionUnit
EESpecific Energym
yyFlow Depthm
vvAverage Velocitym/s
QQDischargem³/s
AACross-sectional Area
ggAcceleration due to Gravitym/s²
If you plot Specific Energy (EE) against depth (yy), you obtain a curve showing that for any given specific energy greater than a minimum value, there are two possible flow depths called alternate depths.

Critical Depth (ycy_c)

The specific depth at which the specific energy is at an absolute minimum for a given discharge. For a rectangular channel, yc=(q2/g)1/3y_c = (q^2 / g)^{1/3}, where qq is the discharge per unit width.

Channel Transitions and Choking

Specific energy is critical for analyzing localized transitions like a raised channel bottom (a bump) or a narrowed width (a contraction). As flow goes over a bump, specific energy decreases. If the bump is too high, the specific energy drops to the absolute minimum (critical depth). If the bump is raised any further, a choke condition occurs: the flow must back up upstream to gain enough energy to pass over the obstruction.

The Froude Number (FrF_r)

The flow state relative to critical depth is classified using the Froude Number (ratio of inertial to gravity forces).

Formula

Mathematical expression.

Fr=vgD F_r = \frac{v}{\sqrt{gD}}

Variables

SymbolDescriptionUnit
FrF_rFroude Numberdimensionless
vvAverage Velocitym/s
ggAcceleration due to Gravitym/s²
DDHydraulic Depth (Area / Top Width)m

Flow Regimes

  • Subcritical Flow (Fr<1F_r < 1): Flow depth is greater than critical depth (y>ycy > y_c). Velocity is slow, flow is tranquil. Waves can propagate upstream.
  • Critical Flow (Fr=1F_r = 1): Flow depth equals critical depth.
  • Supercritical Flow (Fr>1F_r > 1): Flow depth is less than critical depth (y<ycy < y_c). Velocity is high, flow is rapid. Waves cannot propagate upstream.

Specific Energy and Flow Regime Simulator

Adjust the unit discharge (qq) and the flow depth (yy) in a rectangular channel to see how the Specific Energy (EE), Froude Number (FrF_r), and flow regime change.

Calculated Parameters

  • Velocity (vv):3.33 m/s
  • Specific Energy (EE):2.07 m
  • Critical Depth (ycy_c):1.37 m
  • Froude Number (FrF_r):0.87
  • Flow Regime:Subcritical
3.3 m/s
3.3 m/s
3.3 m/s
3.3 m/s
3.3 m/s
y_c = 1.37m

  1. Gradually Varied Flow (GVF)

Analyzing steady flows where the depth changes gradually over long distances.
Unlike uniform flow where depth is constant, Gradually Varied Flow (GVF) occurs when depth changes slowly due to backwater effects from dams or changes in slope. Engineers calculate the continuous water surface profile. Profiles are classified based on the channel slope (Mild, Steep) and the depth relative to normal (yny_n) and critical depth (ycy_c). Common profiles like the M1 curve (backwater behind a dam) determine levee heights.

  1. The Hydraulic Jump and Specific Force

A hydraulic jump is a rapidly varied flow phenomenon that occurs when a high-velocity, supercritical flow abruptly transitions into a slow-moving, subcritical flow, resulting in a churning surface and massive energy loss.
Engineers intentionally induce hydraulic jumps at the base of spillways to safely dissipate kinetic energy before water enters fragile natural riverbeds.
Unlike specific energy (which decreases across a jump), Specific Force (or Momentum Flux) is conserved across the jump. The depths before (y1y_1) and after (y2y_2) the jump are called conjugate depths. They are related by the momentum equation:

Formula

Mathematical expression.

y2y1=12(1+8Fr121) \frac{y_2}{y_1} = \frac{1}{2} \left( \sqrt{1 + 8 F_{r1}^2} - 1 \right)

Variables

SymbolDescriptionUnit
y1y_1Upstream Depthm
y2y_2Downstream Depth (Conjugate Depth)m
Fr1F_{r1}Upstream Froude Numberdimensionless
The energy lost (ΔE\Delta E) during the jump in a rectangular channel is:

Formula

Mathematical expression.

ΔE=E1E2=(y2y1)34y1y2 \Delta E = E_1 - E_2 = \frac{(y_2 - y_1)^3}{4y_1y_2}

Variables

SymbolDescriptionUnit
ΔE\Delta EEnergy Loss (Head Loss)m
E1,E2E_1, E_2Specific Energy Upstream and Downstreamm
y1y_1Upstream Depthm
y2y_2Downstream Depthm

Hydraulic Jump Energy Dissipation

Observe how supercritical flow abruptly transitions to subcritical flow, dissipating immense kinetic energy. Adjust the upstream conditions to see the resulting jump height and energy loss.

Jump Characteristics

Upstream (Supercritical)

v₁: 8.00 m/s

Fr₁: 3.61

E₁: 3.76 m

Downstream (Subcritical)

v₂: 1.73 m/s

Fr₂: 0.36

y₂ (Conjugate): 2.32 m

Energy Loss (ΔE): 1.29 m
y₁
y₂

Engineering Insight

In Water Resources Engineering, the practical application of theoretical formulas often requires careful consideration of real-world variables, such as varying friction coefficients, unpredictable environmental conditions, and changing climate patterns. A rigorous approach to empirical validation and an understanding of the safety margins involved are paramount for resilient infrastructure design.
Key Takeaways
  • Uniform Flow: Governed by Manning's equation, where gravity balances friction. Erodible channels require Tractive Force design to prevent scour.
  • Specific Energy: Energy relative to the bed. It dictates how flow reacts to bumps or contractions, potentially causing a "choke".
  • Flow Regimes: The Froude number classifies flow as subcritical (tranquil) or supercritical (rapid).
  • Gradually Varied Flow: Analyzes long-distance backwater profiles created by downstream obstructions like dams.
  • Hydraulic Jump: A rapid transition from supercritical to subcritical flow, conserving Momentum (Specific Force) while dissipating massive Energy.