Open Channel Flow: Non-Uniform Flow

Learning Objectives

  • Understand non-uniform flow and specific energy concepts.
  • Analyze critical flow and the Froude number.
  • Evaluate hydraulic jumps and their energy dissipation.
  • Understand gradually varied flow profiles.

Specific energy, critical depth, hydraulic jumps, and gradually varied flow profiles.

Non-Uniform Flow Overview

Non-uniform flow occurs when the depth of flow changes along the length of the channel. This happens due to changes in slope, cross-section, or obstructions.

Specific Energy (EE)

Specific energy is the energy per unit weight of fluid relative to the channel bottom.

Specific Energy

Calculates the specific energy as the sum of flow depth and velocity head.

E=y+V22g=y+Q22gA2E = y + \frac{V^2}{2g} = y + \frac{Q^2}{2g \cdot A^2}

Variables

SymbolDescriptionUnit
EESpecific energym
yyFlow depthm
VVMean velocitym/s
ggAcceleration due to gravitym/s2m/s^2
QQDischarge (volumetric flow rate)m3/sm^3/s
AACross-sectional area of flowm2m^2

Specific Energy Curve

For a given discharge QQ, the specific energy curve (EE vs yy) is an asymptotic curve showing two possible depths for any energy greater than the minimum (E>EminE > E_{min}). The upper leg of the curve represents subcritical flow, while the lower leg represents supercritical flow.

Flow Depth Types

  1. Subcritical Depth (y2>ycy_2 > y_c): Slow, deep flow. (Froude Number Fr<1Fr < 1)
  2. Supercritical Depth (y1<ycy_1 < y_c): Fast, shallow flow. (Froude Number Fr>1Fr > 1)

Interactive Simulation

Use the simulation below to explore the relationship between specific energy and flow depth. You can find the critical depth (ycy_c) where specific energy is minimized.

Specific Energy Curve (EE vs yy)

Unit Discharge (q=Q/bq=Q/b):2.50 m²/s
Critical Depth (ycy_c):0.860 m
Min Specific Energy (EminE_{min}):1.291 m
Hover over graph to inspect

The Specific Energy curve shows two possible depths for a given energy E>EminE > E_{min}: a subcritical depth (slow, deep) and a supercritical depth (fast, shallow). ycy_c represents the transition point.

Critical Flow

Critical flow is the flow state where specific energy is minimum for a given discharge. At critical flow, the Froude number is 11.

Froude Number (FrFr)

The ratio of inertial forces to gravity forces. In open channels, it determines whether the flow is subcritical, critical, or supercritical.

Froude Number

Calculates the Froude number using flow velocity and hydraulic depth.

Fr=VgDhFr = \frac{V}{\sqrt{g \cdot D_h}}

Variables

SymbolDescriptionUnit
FrFrFroude numberdimensionless
VVMean velocitym/s
ggAcceleration due to gravitym/s2m/s^2
DhD_hHydraulic depth (A/TA/T where TT is top width)m

Froude Number Simplification for Rectangular Channels

For a rectangular channel, the hydraulic depth DhD_h is exactly equal to the flow depth yy, which simplifies the Froude number calculation to:

Fr=VgyFr = \frac{V}{\sqrt{g \cdot y}}

Critical Depth Overview

The critical depth (ycy_c) corresponds to the depth of minimum specific energy for a given discharge. The formula depends on the channel's cross-sectional geometry.

  • Rectangular Channel: Defined in terms of unit discharge qq.
  • Triangular Channel: Defined in terms of total discharge QQ and side slope mm.

Critical Depth in Rectangular Channels

Calculates critical depth for a rectangular channel using unit discharge.

yc=q2g3y_c = \sqrt[3]{\frac{q^2}{g}}

Variables

SymbolDescriptionUnit
ycy_cCritical depthm
qqDischarge per unit width (Q/bQ/b)m2/sm^2/s
ggAcceleration due to gravitym/s2m/s^2

Critical Depth in Triangular Channels

Calculates critical depth for a symmetric triangular channel.

yc=(2Q2gm2)1/5y_c = \left( \frac{2 \cdot Q^2}{g \cdot m^2} \right)^{1/5}

Variables

SymbolDescriptionUnit
ycy_cCritical depthm
QQTotal dischargem3/sm^3/s
ggAcceleration due to gravitym/s2m/s^2
mmSide slope ratio (Horizontal to Vertical, H:1)dimensionless

Hydraulic Jump

A hydraulic jump is a phenomenon where flow transitions abruptly from supercritical (shallow, fast) to subcritical (deep, slow), resulting in significant energy dissipation, high turbulence, and a rapid rise in the water surface.

  • Applications: Hydraulic jumps are widely used as energy dissipators downstream of spillways, sluice gates, and drop structures to prevent downstream erosion and scouring. They also promote mixing in water treatment facilities.

Sequent Depths (y1,y2y_1, y_2)

The depths immediately before and after the hydraulic jump, also referred to as conjugate depths.

Belanger Equation for Sequent Depths

Relates the downstream sequent depth to the upstream depth and Froude number in a rectangular channel.

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

Variables

SymbolDescriptionUnit
y1y_1Upstream flow depth (supercritical)m
y2y_2Downstream flow depth (subcritical)m
Fr1Fr_1Upstream Froude number (Fr1>1Fr_1 > 1)dimensionless

Energy Dissipation in Hydraulic Jumps

The primary purpose of engineered hydraulic jumps is to safely dissipate destructive kinetic energy.

Head Loss Concept

The specific energy loss across a hydraulic jump in a horizontal rectangular channel is a function only of the upstream (y1y_1) and downstream (y2y_2) sequent depths.

  • The power dissipated per unit width (P/bP/b) is γqΔE\gamma \cdot q \cdot \Delta E.
  • A larger jump (higher y2/y1y_2/y_1 ratio) dissipates exponentially more energy, making it an extremely effective stilling basin mechanism.

Energy Loss in a Hydraulic Jump

Calculates the specific energy head loss across a hydraulic jump in a rectangular channel.

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

Variables

SymbolDescriptionUnit
ΔE\Delta ESpecific energy head lossm
E1E_1Specific energy before the jumpm
E2E_2Specific energy after the jumpm
y1y_1Upstream sequent depthm
y2y_2Downstream sequent depthm

Interactive Simulation

Interact with the simulation below to visualize a hydraulic jump and see how different upstream flow conditions (Froude number, initial depth) influence the downstream sequent depth and the total energy loss.

Hydraulic Jump Simulator

Hydraulic jump is occurring. Upstream Froude number is 0.00 which is supercritical. Sequent depth is 0.00 meters. Downstream Froude number is 0.00 which is subcritical. Energy loss is 0.00 meters.

Input Parameters

5.0 m²/s
0.50 m

Note: A jump only forms if the upstream flow is supercritical (Fr1>1Fr_1 > 1).

Flow Characteristics

Upstream Froude (Fr₁)
0.00
Supercritical
Sequent Depth (y2y_2)
0.00 m
Downstream Froude (Fr₂)
0.00
Subcritical
Energy Loss (ΔE)
0.00 m

ℹ️ Understanding the Hydraulic Jump

A hydraulic jump occurs when a high-velocity, supercritical flow (Fr>1Fr > 1) transitions into a low-velocity, subcritical flow (Fr<1Fr < 1). This abrupt transition causes a sudden rise in water surface elevation and creates significant turbulence.

Engineers use this phenomenon to intentionally dissipate kinetic energy at the base of spillways and dams, preventing downstream erosion. The conjugate depth equation is used to predict the downstream depth (y2y_2) given the upstream conditions.

Water Surface Profiles Computation

Water surface profile computation methods calculate the change in depth over a specific channel length.

Standard Step Method Overview

The Standard Step Method is the most widely used numerical technique for computing gradually varied flow (GVF) profiles. It involves solving the energy equation between two adjacent cross-sections (stations). This method is iterative for unknown depths (y2y_2) given a known starting point (y1y_1) and distance (Δx\Delta x).

Standard Step Method Energy Equation

The fundamental energy balance equation between two adjacent channel cross-sections.

y1+z1+V122g=y2+z2+V222g+hf+hey_1 + z_1 + \frac{V_1^2}{2g} = y_2 + z_2 + \frac{V_2^2}{2g} + h_f + h_e

Variables

SymbolDescriptionUnit
y1,y2y_1, y_2Flow depth at stations 1 and 2m
z1,z2z_1, z_2Elevation of the channel bottom above a datum at stations 1 and 2m
V1,V2V_1, V_2Mean velocity at stations 1 and 2m/s
ggAcceleration due to gravitym/s2m/s^2
hfh_fFriction head loss between sections (calculated as Sf,avgΔxS_{f,avg} \cdot \Delta x)m
heh_eEddy loss (minor losses due to channel expansion or contraction)m

Gradually Varied Flow (GVF)

Gradually varied flow occurs where the depth changes gradually over a long distance, meaning the streamlines are practically parallel and hydrostatic pressure distribution holds. The water surface profile is classified based on the bed slope (S0S_0) and the actual depth relative to critical (ycy_c) and normal (yny_n) depths.

  • M Profiles (Mild Slope, yn>ycy_n > y_c):
    • M1 (Backwater curve): Depth is greater than normal depth (y>yn>ycy > y_n > y_c). Commonly occurs upstream of a dam or a weir where water pools. Depth increases in the direction of flow.
    • M2 (Drawdown curve): Depth is between normal and critical depth (yn>y>ycy_n > y > y_c). Flow accelerates as it approaches a free overfall or sudden channel expansion. Depth decreases in the direction of flow.
    • M3: Flow is supercritical on a mild slope (yn>yc>yy_n > y_c > y), typically occurring immediately downstream of a sluice gate before a hydraulic jump forms.
  • S Profiles (Steep Slope, yc>yny_c > y_n): Profiles (S1, S2, S3) describe flow adjusting on steep inclines. S1 curves follow a hydraulic jump, while S2 and S3 occur near channel slope transitions.
Key Takeaways
  • Specific Energy (EE): The total energy head relative to the channel bottom, composed of depth (yy) and velocity head (V2/2gV^2/2g).
  • Froude Number (FrFr): Measures the flow state (Fr=1Fr = 1 is critical, Fr<1Fr < 1 is subcritical, Fr>1Fr > 1 is supercritical).
  • Hydraulic Jump: A highly turbulent, rapid transition from supercritical to subcritical flow, engineered to dissipate energy.
  • Belanger Equation: Quantifies the relationship between upstream and downstream sequent depths of a hydraulic jump.
  • Gradually Varied Flow (GVF): Describes slow water surface depth variations along a channel, categorized into profiles (e.g., M1, M2, S2) based on bed slope and relative depths.