Froude Number Open Channel Flow: Derivation, Critical Regimes & Solved Examples (2026)

1. Physical Significance of Gravity Waves and Inertial Forces in Open Channels

In open channel hydraulics, flow behavior is governed by the dynamic interplay between fluid inertia and gravitational restoring forces. Analyzing the froude number open channel flow provides hydraulic engineers with the fundamental dimensionless parameter required to categorize free-surface flow regimes, predict surface wave propagation directions, and design stable conveyance structures. Named after the pioneering British naval architect William Froude, this dimensionless index quantifies the ratio of inertial forces to gravitational forces operating within moving water.

                      GRAVITY WAVE DISTURBANCE PROPAGATION
         Subcritical Flow (Fr < 1):       Supercritical Flow (Fr > 1):
           V < c (Waves travel upstream)    V > c (Waves swept downstream)
               <--- c - V                       +-----> V - c
             (Wave spreads both ways)         (Disturbance swept into wedge)

In closed pipe networks under pressure, acoustic pressure wave celerity dominates fluid dynamics. In contrast, free-surface water systems are governed by shallow gravity wave celerity. When analyzing the Froude number in open channels, the magnitude of mean flow velocity relative to surface wave propagation speed determines whether upstream water levels can sense downstream channel controls or if disturbances are swept downstream.

Fundamental Mechanics of Hydraulic Jumps and Stilling Basins

Understanding this dimensionless index is essential for designing spillway chutes, irrigation canals, bridge piers, and energy dissipation basins.

2. Mathematical Derivation of the Froude Number from Wave Celerity

The mathematical origin of the dimensionless ratio stems from the propagation speed of an infinitesimal gravity wave traveling across a shallow water layer of depth $y$.

                   SHALLOW WATER WAVE DERIVATION CONTROL VOLUME
                   ---------------------( Wave Front: c )----
                        Depth y + dy  |   Depth y
                     =======================================> V
                                   Bed (Datum)

Consider a small surface wave of height $dy$ propagating with absolute celerity $c$ across a horizontal, frictionless rectangular channel. By applying the continuity equation across the wave front in a relative coordinate frame moving with the wave at speed $c$:

$$(c – V) y = (c – V – dV)(y + dy)$$

Expanding and neglecting second-order differential products ($dV \cdot dy \approx 0$):

$$c \, y – V \, y = c \, y + c \, dy – V \, y – V \, dy – y \, dV \implies c \, dy – y \, dV = 0 \implies dV = c \frac{dy}{y}$$

Now apply the one-dimensional momentum conservation equation across the wave front:

$$\Delta F_p = \rho Q \Delta V \implies \frac{1}{2} \rho g (y + dy)^2 – \frac{1}{2} \rho g y^2 = \rho (c \, y) (dV)$$

Simplifying hydrostatic pressure terms and substituting $dV = c \frac{dy}{y}$:

$$\rho g y \, dy = \rho c y \left( c \frac{dy}{y} \right) \implies g y \, dy = c^2 \, dy \implies c = \sqrt{g y}$$

The ratio of bulk flow velocity $V$ to small gravity wave celerity $c$ defines the dimensionless parameter:

$$Fr = \frac{V}{c} = \frac{V}{\sqrt{g y}}$$

For non-rectangular channels, water depth $y$ is generalized to hydraulic depth $D_h = A / T$, where $A$ is flow cross-sectional area and $T$ is top water surface width:

$$Fr = \frac{V}{\sqrt{g D_h}} = \frac{Q / A}{\sqrt{g (A / T)}} = \sqrt{\frac{Q^2 T}{g A^3}}$$

3. Flow Regime Classification: Subcritical, Critical, and Supercritical States

The magnitude of the dimensionless parameter establishes three distinct hydraulic regimes:

                    HYDRAULIC FLOW REGIME SPECTRUM
      Fr < 1.0                Fr = 1.0                Fr > 1.0
   [Subcritical] ----------> [Critical] ----------> [Supercritical]
   - Deep, tranquil flow     - Minimum energy head   - Shallow, rapid shooting flow
   - Upstream wave travel    - Control section       - Hydraulic jump potential

In subcritical flow ($Fr 1$), fluid moves faster than surface wave speed, making upstream feedback impossible until a hydraulic jump forms.

Derivation and Applications of Critical Depth in Open Channels

4. Geometric Generalization: Hydraulic Depth across Non-Rectangular Channels

In natural rivers and engineered channels, cross-sectional geometry varies significantly. Calculating the parameter requires computing hydraulic depth $D_h = A / T$:

1. Rectangular Channel: $A = b y$, $T = b \implies D_h = \frac{b y}{b} = y \implies Fr = \frac{V}{\sqrt{g y}}$.
2. Trapezoidal Channel: With bottom width $b$ and side slope $z$ ($1\text{ V}:z\text{ H}$):

$$A = (b + z y) y, \quad T = b + 2 z y \implies D_h = \frac{(b + z y) y}{b + 2 z y}$$

3. Triangular Channel: Symmetrical side slopes $z$:

$$A = z y^2, \quad T = 2 z y \implies D_h = \frac{z y^2}{2 z y} = \frac{y}{2} \implies Fr = \frac{V}{\sqrt{g (y/2)}} = \frac{V \sqrt{2}}{\sqrt{g y}}$$

4. Circular Conduit: Diameter $D$, flow angle $\theta$ (radians):

$$A = \frac{D^2}{8}(\theta – \sin \theta), \quad T = D \sin\left(\frac{\theta}{2}\right) \implies D_h = \frac{D(\theta – \sin\theta)}{8 \sin(\theta/2)}$$

Comprehensive Guide to Hydraulic Jump Regimes and Energy Dissipation

5. Specific Energy Relationship and the Critical State Singularity

Specific energy $E$ defines total energy per unit weight of fluid measured relative to the channel bed:

$$E = y + \frac{V^2}{2g} = y + \frac{Q^2}{2g A^2}$$

Differentiating specific energy with respect to depth $y$ at constant discharge $Q$:

$$\frac{dE}{dy} = 1 – \frac{Q^2}{g A^3} \frac{dA}{dy}$$

Because incremental area $dA = T \, dy \implies \frac{dA}{dy} = T$:

$$\frac{dE}{dy} = 1 – \frac{Q^2 T}{g A^3} = 1 – Fr^2$$

Setting $\frac{dE}{dy} = 0$ to identify minimum specific energy:

$$1 – Fr^2 = 0 \implies Fr = 1$$

This fundamental mathematical result demonstrates that minimum specific energy occurs precisely at critical flow ($Fr = 1.0$).

                    SPECIFIC ENERGY CURVE (E vs y)
         Depth (y)
             ^
             |       / Subcritical Branch (Fr < 1)
             |      /
         yc -+-----o  <-- Critical State (Fr = 1.0, Minimum E)
             |      \
             |       \ Supercritical Branch (Fr > 1)
             +--------+----------------------------> Specific Energy (E)
             0       E_min

6. Solved Numerical Example 1: Rectangular Concrete Irrigation Flume

A rectangular concrete channel has bottom width $b = 4.0\text{ m}$ and carries discharge $Q = 18.0\text{ m}^3/\text{ s}$ at flow depth $y = 1.2\text{ m}$. Compute the dimensionless flow parameter, wave celerity, and classify the regime using free-surface hydraulic flow principles.

Compute flow area and mean velocity:

• Cross-sectional area: $A = b y = 4.0 \times 1.2 = 4.80\text{ m}^2$.

• Mean flow velocity: $V = \frac{Q}{A} = \frac{18.0}{4.80} = 3.750\text{ m/s}$.

• Top surface width: $T = b = 4.0\text{ m}$.

• Hydraulic depth: $D_h = \frac{A}{T} = 1.20\text{ m}$.

Compute gravity wave celerity:

$$c = \sqrt{g D_h} = \sqrt{9.81 \times 1.20} = \sqrt{11.772} = 3.431\text{ m/s}$$

Calculate the dimensionless flow parameter:

$$Fr = \frac{V}{c} = \frac{3.750}{3.431} = 1.093$$

Because $Fr = 1.093 > 1.0$, the flow is supercritical (rapid). Flow velocity exceeds wave celerity ($3.75\text{ m/s} > 3.43\text{ m/s}$), meaning surface waves cannot propagate upstream.

7. Solved Numerical Example 2: Symmetrical Trapezoidal Drainage Canal

A trapezoidal earthen canal has bottom width $b = 5.0\text{ m}$, side slopes $z = 1.5$ ($1\text{ V}:1.5\text{ H}$), and flow depth $y = 1.8\text{ m}$. The canal conveys discharge $Q = 35.0\text{ m}^3/\text{ s}$. Determine the flow regime.

Calculate geometric parameters:

• Cross-sectional area: $A = (b + z y) y = [5.0 + 1.5(1.8)](1.8) = (5.0 + 2.7)(1.8) = 7.7 \times 1.8 = 13.860\text{ m}^2$.

• Top surface width: $T = b + 2 z y = 5.0 + 2(1.5)(1.8) = 5.0 + 5.4 = 10.400\text{ m}$.

• Hydraulic depth: $D_h = \frac{A}{T} = \frac{13.860}{10.400} = 1.3327\text{ m}$.

• Mean flow velocity: $V = \frac{Q}{A} = \frac{35.0}{13.860} = 2.5253\text{ m/s}$.

Compute wave celerity:

$$c = \sqrt{g D_h} = \sqrt{9.81 \times 1.3327} = \sqrt{13.0737} = 3.6158\text{ m/s}$$

Calculate the dimensionless flow index:

$$Fr = \frac{V}{\sqrt{g D_h}} = \frac{2.5253}{3.6158} = 0.6984$$

Because $Fr = 0.698 < 1.0$, the flow is subcritical (tranquil). Any surface disturbance travels upstream at speed $c – V = 3.616 – 2.525 = 1.091\text{ m/s}$.

8. Solved Numerical Example 3: Partially Filled Circular Storm Culvert

A precast concrete culvert with internal diameter $D = 1.50\text{ m}$ operates at flow depth $y = 0.90\text{ m}$ with discharge $Q = 2.50\text{ m}^3/\text{ s}$. Evaluate the flow state using dimensionless open channel analysis.

Calculate central angle $\theta$ subtended by water surface:

$$\cos\left(\frac{\theta}{2}\right) = 1 – \frac{2y}{D} = 1 – \frac{2(0.90)}{1.50} = 1 – 1.20 = -0.20 \implies \frac{\theta}{2} = \arccos(-0.20) = 1.7722\text{ rad} \implies \theta = 3.5443\text{ rad}$$

Calculate geometric section properties:

• Flow area: $A = \frac{D^2}{8}(\theta – \sin\theta) = \frac{(1.5)^2}{8}[3.5443 – \sin(3.5443)] = 0.28125[3.5443 – (-0.3887)] = 0.28125(3.9330) = 1.1062\text{ m}^2$.

• Top surface width: $T = D \sin(\theta/2) = 1.50 \sin(1.7722) = 1.50(0.9798) = 1.4697\text{ m}$.

• Hydraulic depth: $D_h = \frac{A}{T} = \frac{1.1062}{1.4697} = 0.7527\text{ m}$.

• Mean velocity: $V = \frac{Q}{A} = \frac{2.50}{1.1062} = 2.2599\text{ m/s}$.

Calculate the dimensionless flow parameter:

$$Fr = \frac{V}{\sqrt{g D_h}} = \frac{2.2599}{\sqrt{9.81 \times 0.7527}} = \frac{2.2599}{\sqrt{7.3837}} = \frac{2.2599}{2.7173} = 0.8317$$

Because $Fr = 0.832 < 1.0$, flow in the circular culvert is subcritical.

9. Practical Hydraulic Engineering Applications: Jumps, Weirs, and Chutes

The dimensionless flow index serves as the primary governing parameter across major hydraulic infrastructure designs:

9.1 Hydraulic Jump Formation

When rapid supercritical flow ($Fr_1 > 1.0$) transitions into tranquil subcritical flow ($Fr_2 < 1.0$), a hydraulic jump occurs. The upstream value $Fr_1$ dictates jump classification, energy dissipation efficiency, and stilling basin geometry according to USBR Engineering Monograph No. 25.

9.2 Critical Flow Measurement Structures

Broad-crested weirs, flumes, and channel drop structures force flow through the critical transition ($Fr = 1.0$). Because a unique relationship connects discharge to critical depth, measuring upstream depth enables accurate flow metering without backwater interference.

9.3 Spillway Chutes and Cavitation Protection

Steep dam spillway chutes accelerate flow to high values ($Fr > 4.0$). High velocities demand careful evaluation of surface aeration and boundary layer development per ASCE Environmental & Water Resources Guidelines to prevent cavitation damage on concrete surfaces. Standards from NIST Technical Resources provide calibrated flow measurement tolerances.

10. Froude Number Synthesis

Mastering the froude number open channel flow provides civil and environmental engineers with the fundamental hydrodynamic foundation required to analyze free-surface systems. By linking bulk flow velocity to shallow water wave celerity, this dimensionless parameter establishes clear boundaries between tranquil subcritical states and rapid supercritical regimes. Whether evaluating flood conveyance in natural river channels, optimizing weir control sections, or sizing spillway stilling basins, this parameter remains the primary governing metric in open channel hydraulic design.

References & Standards Cited

• USBR: Design of Small Dams, 3rd ed., United States Bureau of Reclamation, Denver, CO.

• ASCE: Hydraulic Design of Stilling Basins and Energy Dissipators, American Society of Civil Engineers, Reston, VA.

• NIST: Fluid Mechanics Standards and Flow Measurement Guidelines, National Institute of Standards and Technology, Gaithersburg, MD.

• Chow, V. T.: Open-Channel Hydraulics, McGraw-Hill, New York.

• Henderson, F. M.: Open Channel Flow, Macmillan Publishing, New York.

• Sturm, T. W.: Open Channel Hydraulics, 2nd ed., McGraw-Hill, New York.

Frequently Asked Questions (FAQ)

What does a Froude number equal to 1.0 physically indicate?

A value of 1.0 indicates critical flow, where flow velocity exactly matches surface gravity wave celerity ($V = sqrt{g D_h}$). At this state, specific energy is minimized for a given discharge.

Why can surface waves not travel upstream in supercritical flow?

In supercritical flow ($Fr > 1$), mean water velocity exceeds wave propagation speed ($V > c$). The downstream bulk movement sweeps surface ripples downstream faster than they can travel upstream.

What is the difference between hydraulic depth and hydraulic radius?

Hydraulic depth ($D_h = A / T$) is the ratio of flow area to top surface width and governs wave celerity and gravity forces. Hydraulic radius ($R_h = A / P$) is the ratio of flow area to wetted perimeter and governs boundary friction in Manning equations.

How does channel slope influence the flow regime?

Mild slopes sustain subcritical flow, steep slopes generate supercritical flow, and critical slopes maintain critical flow.

Why is the Froude parameter used instead of the Reynolds parameter for open channel transitions?

Free-surface flows are dominated by gravitational restoring forces and surface wave dynamics rather than viscous shear forces. The dimensionless ratio of inertial to gravity forces governs free-surface profiles and hydraulic transitions.