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.
Variables
| Symbol | Description | Unit |
|---|---|---|
| Specific energy relative to the local channel bottom | m | |
| Flow depth | m | |
| Kinetic-energy correction coefficient | - | |
| Section-average velocity | m/s | |
| Discharge |
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 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.
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.
Small-Disturbance Gravity-Wave Speed
Approximate shallow-water wave celerity for a prismatic section.
Flow Regimes and Information Direction
- Subcritical, : , so small gravity disturbances can propagate both upstream and downstream. Downstream controls can influence upstream depth.
- Critical, : and specific energy is minimum for the specified discharge under the stated assumptions.
- Supercritical, : , 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
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.
Role of the Kinetic-Energy Correction Coefficient
If is retained as a constant different from one, the critical condition becomes . Introductory calculations commonly take , 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.
Variables
| Symbol | Description | Unit |
|---|---|---|
| Discharge per unit channel width | ||
| Critical depth | m |
How the Critical Condition Is Obtained
For fixed in a prismatic section, differentiate with respect to . Because , the condition 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 , total energy may be conserved while specific energy relative to the new bed decreases by approximately 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.
Variables
| Symbol | Description | Unit |
|---|---|---|
| Momentum function or specific force | ||
| Momentum correction coefficient | - | |
| Depth of the flow-area centroid below the free surface | m |
Rectangular Momentum Function per Unit Width
Specific force per unit width for beta equal to one.
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.
Rectangular Hydraulic-Jump Energy Loss
Specific-energy loss between conjugate depths.
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 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 , and verify conjugate depth ratios in real time.
Open Channel Flume & Hydraulic Jump 3D
Interactive 3D laboratory flume with sluice gate, supercritical jet, surface roller vortex, and energy dissipation.
Hydraulic Parameters
Approximate Hydraulic-Jump Classification
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 is the downstream depth required for a free jump at the assumed approach state. Tailwater lower than 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.
A hydraulic jump requires supercritical approach flow (). The rectangular-channel conjugate-depth equation assumes a horizontal prismatic channel and negligible friction through the short jump region.
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.
Variables
| Symbol | Description | Unit |
|---|---|---|
| Longitudinal coordinate positive downstream | m | |
| Bed slope positive downward in the flow direction | - | |
| Friction 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 , differentiating specific energy with depth gives . 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 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: .
- Steep, S: .
- Critical, C: .
- Horizontal, H: ; 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 and .
- Zone 2: actual depth lies between and .
- Zone 3: actual depth lies below both and .
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 profiles
Normal and critical depths organize profile zones.
Not Every Letter-Number Profile Exists
Profile classification is constrained by the ordering of and . 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.
Variables
| Symbol | Description | Unit |
|---|---|---|
| Signed longitudinal distance from section 1 to section 2 | m | |
| Representative average friction slope between sections | - |
Direct-Step Computation
- Choose a sequence of physically admissible depths starting from the known control depth.
- At each depth compute , , , , , and .
- Average between adjacent sections using a stated method.
- Evaluate with one consistent downstream-positive coordinate and section order.
- Accumulate distances in the physically correct stepping direction.
- Reduce the depth increment near critical flow or where section properties change rapidly.
- Check that the reconstructed profile approaches the expected control or asymptote rather than relying only on arithmetic convergence.
Direct-step reach progression
Energy change and average friction slope estimate reach length.
Direct-Step Sign Convention Must Stay Fixed
A negative 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.
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.
Variables
| Symbol | Description | Unit |
|---|---|---|
| Lateral discharge per unit channel length; positive for inflow under the chosen convention |
Constant-Q GVF Is Not an SVF Equation
When changes with , 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
- Determine whether the dominant depth change is rapid, gradual, or caused by distributed discharge change.
- Compute and identify subcritical, critical, or supercritical regime.
- Compute and, for a positive-slope prismatic reach, .
- Identify physical controls and the direction from which they influence the profile.
- Use specific energy for alternate-depth and control problems.
- Use momentum for hydraulic jumps and other short highly dissipative transitions.
- Use direct step or standard step only where GVF assumptions remain credible.
- Compare calculated depths with tailwater, freeboard, erosion, cavitation, aeration, and structure constraints.
- Verify limiting behavior and conservation relationships before accepting the numerical profile.
- 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 when .
- 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 , , and , not memorized blindly.
- Direct-step and standard-step methods require consistent signs, physically admissible roots, appropriate step size, and independent hydraulic verification.