Open Channel Flow: Non-Uniform Flow

Learning Objectives

  • Distinguish rapidly varied, gradually varied, spatially varied, steady, and uniform open-channel flow.
  • Use specific energy and Froude number to classify flow regime and locate critical controls.
  • Derive the general critical-flow condition and rectangular critical-depth relation.
  • Distinguish alternate depths from conjugate depths and choose energy or momentum appropriately.
  • Analyze rectangular hydraulic jumps, energy dissipation, and tailwater control.
  • Derive and interpret the gradually varied flow equation from the energy equation.
  • Classify mild, steep, critical, horizontal, and adverse water-surface profiles and identify their control direction.
  • Apply direct-step and standard-step methods with consistent sign conventions, step sizes, and physical checks.
  • Recognize when lateral inflow, abrupt transitions, or nonhydrostatic effects invalidate constant-discharge GVF assumptions.

Why Non-Uniform Flow Requires Different Tools

Whenever depth changes along a channel, the local balance among gravity, pressure, inertia, and boundary resistance changes. The correct analysis depends on how fast depth changes and whether discharge remains constant. Energy methods are especially useful for controls and gradual profiles; momentum methods are especially useful across short, highly dissipative transitions such as hydraulic jumps.

Rapidly Varied Flow

Open-channel flow in which depth changes appreciably over a longitudinal distance comparable with the flow depth.

Gradually Varied Flow

Open-channel flow in which depth changes slowly enough along the channel that hydrostatic pressure and local uniform-flow resistance approximations remain useful.

Spatially Varied Flow

Open-channel flow in which discharge changes with distance because of distributed lateral inflow or outflow.

Flow-Classification Map

  • RVF: hydraulic jumps, flow under gates, abrupt transitions, spillway toes, and short control regions.
  • GVF: backwater and drawdown curves extending over many flow depths.
  • SVF: side weirs, gutters, irrigation laterals, and channels receiving distributed runoff.
  • Steady and uniform remain separate descriptors: a time-invariant backwater curve is steady but non-uniform.

Specific Energy

Mechanical energy per unit weight measured relative to the local channel bottom at a section.

Specific Energy

Flow depth plus kinetic-energy-corrected velocity head at a channel section.

E=y+αV22gE=y+\alpha\frac{V^2}{2g}E=y+αQ22gA2E=y+\alpha\frac{Q^2}{2gA^2}

Variables

SymbolDescriptionUnit
EESpecific energy relative to the local channel bottomm
yyFlow depthm
α\alphaKinetic-energy correction coefficient-
VVSection-average velocitym/s
QQDischargem3/sm^3/s

Specific-Energy Curve and Alternate Depths

At fixed discharge and geometry, the specific-energy curve has a minimum at critical flow. For an energy level above the minimum, two positive depths may exist: a shallow, high-velocity supercritical depth and a deeper, low-velocity subcritical depth. These alternate depths have equal discharge and specific energy; they do not generally have equal momentum function.

Specific-energy curveCritical depth is the minimum-energy state for fixed discharge.depthspecific energycritical point

Specific-energy curve

Critical depth is the minimum-energy state for fixed discharge.

Hydraulic Depth

Flow area divided by free-surface top width.

Hydraulic Depth

Characteristic depth for general-section Froude-number analysis.

Dh=ATD_h=\frac{A}{T}

Froude Number

A dimensionless measure comparing mean-flow speed with the gravity-wave speed associated with the section hydraulic depth.

Froude Number

Flow-regime parameter for a general open-channel section.

Fr=VgDh=VgA/TFr=\frac{V}{\sqrt{gD_h}} =\frac{V}{\sqrt{gA/T}}

Small-Disturbance Gravity-Wave Speed

Approximate shallow-water wave celerity for a prismatic section.

c=gDhc=\sqrt{gD_h}

Flow Regimes and Information Direction

  • Subcritical, Fr<1Fr<1: V<cV<c, so small gravity disturbances can propagate both upstream and downstream. Downstream controls can influence upstream depth.
  • Critical, Fr=1Fr=1: V=cV=c and specific energy is minimum for the specified discharge under the stated assumptions.
  • Supercritical, Fr>1Fr>1: V>cV>c, so small disturbances cannot travel upstream against the flow. Upstream controls govern the downstream supercritical profile.

Explore Energy and Alternate Depths

Use the interactive solver to move across the specific-energy curve, identify the critical minimum, and compare shallow and deep alternate depths at equal discharge and energy.

Specific Energy and Alternate Depths

Learning objective: Relate discharge and flow depth to the specific-energy curve, critical depth, Froude number, and alternate depths.

Rectangular channel with a numerically bracketed alternate-depth solution.

Unit discharge
5.000 m²/s
Critical depth
1.3659 m
Minimum energy
2.0489 m
Selected energy
2.0663 m
Selected regime
Subcritical (Fr=0.869)
Alternate depth
1.24734 m
Alternate regime
Supercritical (Fr=1.146)
Energy residual
-9.69e-11 m
Hydraulics interactive visualizationRelate discharge and flow depth to the specific-energy curve, critical depth, Froude number, and alternate depths.depth y (log scale)energy E (log scale)criticalselectedalternate

For any energy above the minimum, the shallow and deep roots lie on opposite sides of critical depth. The logarithmic axes keep both roots visible when their magnitudes differ greatly.

Critical Depth

The flow depth at which specific energy is minimum for the specified discharge and channel geometry.

General Critical-Flow Condition

Critical condition for a prismatic section when alpha is approximately constant and taken as one.

Q2TgA3=1\frac{Q^2T}{gA^3}=1

Role of the Kinetic-Energy Correction Coefficient

If α\alpha is retained as a constant different from one, the critical condition becomes αQ2T/(gA3)=1\alpha Q^2T/(gA^3)=1. Introductory calculations commonly take α≈1\alpha\approx1, but measured non-uniform velocity distributions can justify a correction.

Rectangular Critical Depth and Minimum Specific Energy

Critical relations for unit discharge q=Q/b with alpha equal to one.

yc=(q2g)1/3y_c=\left(\frac{q^2}{g}\right)^{1/3}Emin⁡=32ycE_{\min}=\frac{3}{2}y_c

Variables

SymbolDescriptionUnit
qqDischarge per unit channel widthm2/sm^2/s
ycy_cCritical depthm

How the Critical Condition Is Obtained

For fixed QQ in a prismatic section, differentiate E(y)E(y) with respect to yy. Because dA/dy=TdA/dy=T, the condition dE/dy=0dE/dy=0 reduces to the general critical relation. The result is geometric as well as dynamic: a nonrectangular channel reaches critical flow at the depth where its changing area and top width satisfy the same minimum-energy condition.

Control Section

A section where the hydraulic state establishes a unique depth-discharge relation that governs an adjacent water-surface profile.

Critical Controls and Choking

Broad-crested weirs, flumes, free overfalls, and suitably shaped contractions may create critical control. If a bed rise or contraction requires a minimum specific energy greater than the energy available to the approaching state, the original state cannot pass unchanged. Upstream depth then adjusts, or discharge changes, until the control can be satisfied. This condition is called choking.

Specific Energy Is Referenced to the Local Bed

Across a bed rise of height Δz\Delta z, total energy may be conserved while specific energy relative to the new bed decreases by approximately Δz\Delta z when losses are neglected. Do not compare specific-energy values from different bed elevations without accounting for the datum change.

Momentum Function

A section function combining momentum flux and the hydrostatic pressure moment, used to relate rapidly varied states when a short-reach momentum balance is appropriate.

General Momentum Function

Specific-force function for a hydrostatic open-channel section.

M=βQ2gA+AyˉM=\beta\frac{Q^2}{gA}+A\bar y

Variables

SymbolDescriptionUnit
MMMomentum function or specific forcem3m^3
β\betaMomentum correction coefficient-
yˉ\bar yDepth of the flow-area centroid below the free surfacem

Rectangular Momentum Function per Unit Width

Specific force per unit width for beta equal to one.

M′=q2gy+y22M'=\frac{q^2}{gy}+\frac{y^2}{2}

Conjugate Depths

The upstream and downstream depths of a hydraulic jump that satisfy the applicable momentum balance for the same discharge.

Hydraulic Jump

A rapid transition from supercritical to subcritical flow accompanied by strong turbulence and irreversible mechanical-energy dissipation.

Rectangular Sequent-Depth Relation

Conjugate-depth relation for a horizontal rectangular jump with beta approximately one and negligible external streamwise force over the short jump.

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

Rectangular Hydraulic-Jump Energy Loss

Specific-energy loss between conjugate depths.

ΔE=E1−E2=(y2−y1)34y1y2\Delta E=E_1-E_2=\frac{(y_2-y_1)^3}{4y_1y_2}

Why Momentum Is Used across a Jump

A hydraulic jump converts organized mechanical energy into turbulence, heat, air entrainment, surface waves, and sound, so energy is not conserved between the two uniform sections. Over the short jump reach, however, momentum can be balanced accurately when pressure is hydrostatic at the end sections and bed friction and other external streamwise forces are small relative to momentum and pressure forces.

Hydraulic jump transitionShallow supercritical flow transitions to deeper subcritical flow.supercriticalrollersubcritical

Hydraulic jump transition

Shallow supercritical flow transitions to deeper subcritical flow.

Interactive 3D Flume & Hydraulic Jump Simulation

Explore rapidly varied flow in a 3D glass flume. Adjust the upstream sluice gate and discharge to induce hydraulic jumps across undular, weak, oscillating, and steady regimes. Trace the spatial EGL and HGL ribbons, energy dissipation ΔE\Delta E, and verify conjugate depth ratios y2/y1y_2/y_1 in real time.

Open Channel Flume & Hydraulic Jump 3D

Interactive 3D laboratory flume with sluice gate, supercritical jet, surface roller vortex, and energy dissipation.

Loading 3D Hydraulic Flume Simulation…

Hydraulic Parameters

Discharge (QQ)0.25 m³/s
Gate Opening (aa)12 cm
Channel Width (bb)1.00 m
Bed Slope (S0S_0)0.20%
Initial Depth y1y_10.074 m
Initial Froude Fr1Fr_13.93 (Supercritical)
Conjugate Depth y2y_20.378 m
Sequent Froude Fr2Fr_20.34 (Subcritical)
Head Loss ΔHL\Delta H_L0.249 m
Dissipation Efficiency38.4%

Approximate Hydraulic-Jump Classification

Upstream Froude numberCommon descriptionTypical behavior
1.01.0–1.71.7UndularSurface waves; limited dissipation
1.71.7–2.52.5WeakSmall roller; low energy loss
2.52.5–4.54.5OscillatingUnstable jet and large surface oscillations
4.54.5–99StableWell-developed roller and effective dissipation
>9>9StrongHighly energetic jump with large loads and aeration

The boundaries are empirical rather than exact and can shift with geometry, approach turbulence, tailwater, and basin configuration.

Tailwater Determines Where the Jump Can Stand

The calculated sequent depth y2y_2 is the downstream depth required for a free jump at the assumed approach state. Tailwater lower than y2y_2 tends to move the jump downstream; sufficiently high tailwater can move it upstream or submerge it. Stilling-basin design therefore compares the sequent-depth requirement with the downstream rating curve over the full operating range.

Explore Hydraulic-Jump Behavior

Use the interactive simulation to vary the approach state, observe the conjugate depth and jump classification, and compare energy loss with the momentum-function balance.

Hydraulic Jump Simulator

Learning objective: Connect upstream Froude number to sequent depth, momentum balance, jump classification, and mechanical-energy dissipation.

Hydraulics interactive visualizationConnect upstream Froude number to sequent depth, momentum balance, jump classification, and mechanical-energy dissipation.y₁y₂V₁V₂

A hydraulic jump requires supercritical approach flow (Fr1>1Fr_1>1). The rectangular-channel conjugate-depth equation assumes a horizontal prismatic channel and negligible friction through the short jump region.

Steady jump
Upstream
V₁ = 10.00 m/s
Fr₁ = 4.52
Downstream
y₂ = 2.953 m
Fr₂ = 0.31
Specific-energy loss2.498 m
Energy retained55.4%
Momentum-function check8.9e-16

Gradually Varied Flow Assumptions

The classical GVF equation assumes steady, one-dimensional, slowly varying flow in a prismatic channel; approximately hydrostatic pressure; known roughness; small bed slope; and no significant lateral inflow or outflow. Abrupt controls, strong curvature, hydraulic jumps, and rapidly changing geometry require other formulations.

Gradually Varied Flow Equation

Depth gradient for a prismatic reach with downstream-positive x coordinate.

dydx=S0−Sf1−Fr2\frac{dy}{dx}=\frac{S_0-S_f}{1-Fr^2}

Variables

SymbolDescriptionUnit
xxLongitudinal coordinate positive downstreamm
S0S_0Bed slope positive downward in the flow direction-
SfS_fFriction or energy slope-

Derivation Logic for the GVF Equation

Write total head as channel elevation plus specific energy. Along the flow direction, energy decreases at the friction slope, while bed elevation changes at the bed slope. For a prismatic section with α≈1\alpha\approx1, differentiating specific energy with depth gives dE/dy=1−Fr2dE/dy=1-Fr^2. Combining these relations yields the GVF equation and shows why the denominator becomes singular at critical flow.

Reading the Sign of dy/dx

The numerator S0−SfS_0-S_f indicates whether gravity supplies more or less slope than local resistance demands. The denominator changes sign between subcritical and supercritical flow. Evaluating both signs predicts whether depth rises or falls downstream before any numerical integration is attempted.

Normal Depth

The uniform-flow depth corresponding to the specified discharge, channel geometry, roughness, and positive bed slope.

Slope Classification from Normal and Critical Depths

  • Mild, M: yn>ycy_n>y_c.
  • Steep, S: yn<ycy_n<y_c.
  • Critical, C: yn=ycy_n=y_c.
  • Horizontal, H: S0=0S_0=0; ordinary positive-slope normal depth does not exist.
  • Adverse, A: the bed rises in the flow direction; ordinary positive-slope normal depth does not exist.

Profile Zones and Common Shapes

  • Zone 1: actual depth lies above both yny_n and ycy_c.
  • Zone 2: actual depth lies between yny_n and ycy_c.
  • Zone 3: actual depth lies below both yny_n and ycy_c.

Typical examples are M1 backwater upstream of a downstream control, M2 drawdown toward a critical control, S2 supercritical flow downstream of an upstream critical control, and S3 flow rising toward steep-slope normal depth. Horizontal and adverse slopes do not have finite normal depth, so only the physically meaningful profile labels are used.

Gradually varied flow profilesNormal and critical depths organize profile zones.normal depthcritical depthprofile

Gradually varied flow profiles

Normal and critical depths organize profile zones.

Not Every Letter-Number Profile Exists

Profile classification is constrained by the ordering of yny_n and ycy_c. Determine slope type first, then locate actual depth relative to the available reference depths. Memorizing all letter-number combinations without this logic is a common source of errors.

Control Direction for Profile Computation

Subcritical GVF is normally computed upstream from a known downstream control because gravity-wave information can propagate upstream. Supercritical GVF is normally computed downstream from an upstream control because information cannot travel upstream against the flow. A hydraulic jump may connect the two regimes, but the GVF equation cannot be stepped through the jump itself.

Direct-Step Method

Reach length between two known depths in a prismatic constant-discharge channel.

Δx=E2−E1S0−Sˉf\Delta x=\frac{E_2-E_1}{S_0-\bar S_f}

Variables

SymbolDescriptionUnit
Δx\Delta xSigned longitudinal distance from section 1 to section 2m
Sˉf\bar S_fRepresentative average friction slope between sections-

Direct-Step Computation

  1. Choose a sequence of physically admissible depths starting from the known control depth.
  2. At each depth compute AA, RR, VV, FrFr, EE, and SfS_f.
  3. Average SfS_f between adjacent sections using a stated method.
  4. Evaluate Δx=(E2−E1)/(S0−Sˉf)\Delta x=(E_2-E_1)/(S_0-\bar S_f) with one consistent downstream-positive coordinate and section order.
  5. Accumulate distances in the physically correct stepping direction.
  6. Reduce the depth increment near critical flow or where section properties change rapidly.
  7. Check that the reconstructed profile approaches the expected control or asymptote rather than relying only on arithmetic convergence.
Direct-step reach progressionEnergy change and average friction slope estimate reach length.section 1energy changesection 2

Direct-step reach progression

Energy change and average friction slope estimate reach length.

Direct-Step Sign Convention Must Stay Fixed

A negative Δx\Delta x does not automatically indicate an algebra error; it may simply mean section 2 lies upstream of section 1 under a downstream-positive coordinate. Label section order and coordinate direction before substitution and interpret the sign physically.

Standard-Step Energy Equation

Energy balance between two surveyed stations with losses and kinetic-energy correction factors.

z1+y1+α1V122g=z2+y2+α2V222g+hLz_1+y_1+\alpha_1\frac{V_1^2}{2g} = z_2+y_2+\alpha_2\frac{V_2^2}{2g}+h_L

When Standard Step Is Preferred

Standard-step analysis is suited to surveyed natural channels, nonprismatic reaches, changing roughness, bridges, culverts, and structures because the station spacing is known and the unknown depth is solved iteratively. Modern one-dimensional water-surface-profile software extends this same energy-balance framework with contraction, expansion, ineffective-flow, and structure-loss models.

Step Size Is a Hydraulic Modeling Choice

Large steps can skip controls, hide multiple roots, or violate the slowly varying assumption. Near critical flow, abrupt geometry changes, bridges, or rapidly changing friction, shorten the step and use a solver that brackets physically admissible depths rather than relying only on unconstrained Newton iteration.

Spatially Varied Discharge

Continuity relation for distributed lateral inflow or outflow.

dQdx=ql\frac{dQ}{dx}=q_l

Variables

SymbolDescriptionUnit
qlq_lLateral discharge per unit channel length; positive for inflow under the chosen conventionm2/sm^2/s

Constant-Q GVF Is Not an SVF Equation

When QQ changes with xx, entering or leaving water carries momentum and may have a different velocity direction from the main flow. Continuity, energy, and momentum relations require additional terms; ordinary constant-discharge GVF formulas are not sufficient by themselves.

Non-Uniform Flow Decision Workflow

  1. Determine whether the dominant depth change is rapid, gradual, or caused by distributed discharge change.
  2. Compute FrFr and identify subcritical, critical, or supercritical regime.
  3. Compute ycy_c and, for a positive-slope prismatic reach, yny_n.
  4. Identify physical controls and the direction from which they influence the profile.
  5. Use specific energy for alternate-depth and control problems.
  6. Use momentum for hydraulic jumps and other short highly dissipative transitions.
  7. Use direct step or standard step only where GVF assumptions remain credible.
  8. Compare calculated depths with tailwater, freeboard, erosion, cavitation, aeration, and structure constraints.
  9. Verify limiting behavior and conservation relationships before accepting the numerical profile.
Key Takeaways
  • RVF, GVF, and SVF are distinguished by the spatial rate of depth change and whether discharge changes along the reach.
  • Specific energy is referenced to the local bed; alternate depths share energy, while hydraulic-jump conjugate depths share momentum function under the jump assumptions.
  • Critical flow is the minimum-specific-energy state and satisfies Q2T/(gA3)=1Q^2T/(gA^3)=1 when α≈1\alpha\approx1.
  • Froude number also encodes information direction: subcritical flow can respond to downstream controls, while supercritical flow is controlled from upstream.
  • Hydraulic jumps dissipate energy but can be related through momentum; their location depends on tailwater as well as the approach state.
  • The GVF equation follows directly from the energy equation and becomes singular near critical flow, where the gradual-flow model loses validity.
  • Profile labels follow from slope type and depth zone; they should be reasoned from yny_n, ycy_c, and FrFr, not memorized blindly.
  • Direct-step and standard-step methods require consistent signs, physically admissible roots, appropriate step size, and independent hydraulic verification.