Critical Depth Open Channel Flow: Specific Energy Derivation & Solved Examples (2026)

1. Fundamentals of Specific Energy and Critical State Hydraulics

In free-surface hydraulic engineering, accurate water surface profile calculations require identifying hydraulic control sections. Calculating the critical depth open channel flow provides hydraulic engineers with the foundational baseline for analyzing specific energy minimums, channel transitions, and water surface controls. Critical depth represents the unique water depth at which a given discharge passes through a channel cross-section with minimum total specific energy head.

                    SPECIFIC ENERGY RELATIONSHIP
         Depth (y)
             ^
             |       / Subcritical Regime (Fr < 1, y > yc)
             |      /
         yc -+-----o  <-- Critical Depth yc (Minimum Specific Energy E_min)
             |      \
             |       \ Supercritical Regime (Fr > 1, y < yc)
             +--------+----------------------------> Specific Energy (E)
             0       E_min

Whenever open channel flow transitions between tranquil subcritical flow and rapid supercritical flow, the water surface must pass through this singular depth. Evaluating critical depth in open channels allows engineers to size flow metering weirs, calculate backwater profiles, and predict hydraulic jump locations across civil infrastructure works.

Comprehensive Guide to Froude Number and Critical Regimes

The critical state represents an unstable water surface condition where small energy changes produce large depth fluctuations, making it ideal for discharge measurement.

2. Mathematical Derivation of Critical Depth from Minimum Specific Energy

Specific energy $E$ is defined as the total mechanical energy head of flowing water per unit weight, measured relative to the channel bottom datum:

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

where $y$ is flow depth, $V$ is mean flow velocity, $Q$ is total discharge, $A$ is flow cross-sectional area, and $g$ is gravitational acceleration.

To establish the mathematical condition for minimum specific energy at constant discharge $Q$, differentiate $E$ with respect to depth $y$ and set the derivative to zero:

$$\frac{dE}{dy} = \frac{d}{dy} \left( y + \frac{Q^2}{2g A^2} \right) = 1 – \frac{2 Q^2}{2g A^3} \frac{dA}{dy} = 1 – \frac{Q^2}{g A^3} \frac{dA}{dy} = 0$$

From channel geometry, the differential area element $dA$ corresponding to depth change $dy$ is $dA = T \, dy$, where $T$ is the top water surface width:

$$\frac{dA}{dy} = T$$

Substituting $\frac{dA}{dy} = T$ into the derivative yields the universal governing criterion for critical flow:

$$1 – \frac{Q^2 T}{g A^3} = 0 \implies \frac{Q^2 T_c}{g A_c^3} = 1$$

Rearranging terms in terms of mean velocity $V_c = Q / A_c$ and hydraulic depth $D_{h,c} = A_c / T_c$:

$$\frac{V_c^2}{g (A_c / T_c)} = \frac{V_c^2}{g D_{h,c}} = Fr^2 = 1 \implies Fr = \frac{V_c}{\sqrt{g D_{h,c}}} = 1.0$$

This thermodynamic proof demonstrates that critical depth occurs precisely when the Froude number equals unity.

3. Governing Formulations for Rectangular, Trapezoidal, and Triangular Channels

Evaluating the governing criterion $\frac{Q^2 T_c}{g A_c^3} = 1$ across standard cross-sectional geometries yields distinct analytical and numerical solutions:

                    CROSS SECTIONAL CRITICAL DEPTH FORMULAS
   1. Rectangular Channel:      2. Triangular Channel:       3. Trapezoidal Channel:
      yc = (q^2 / g)^(1/3)         yc = (2 Q^2 / g z^2)^(1/5)   Q^2 (b + 2 z yc) / [g (b + z yc)^3 yc^3] = 1

3.1 Rectangular Cross-Section

For a rectangular channel of constant bottom width $b$, flow area is $A = b y$ and top width is $T = b$. Defining unit discharge $q = Q / b$:

$$\frac{Q^2 b}{g (b y_c)^3} = 1 \implies \frac{Q^2}{g b^2 y_c^3} = 1 \implies \frac{q^2}{g y_c^3} = 1 \implies y_c = \sqrt[3]{\frac{q^2}{g}}$$

The minimum specific energy in a rectangular channel is: $E_{min} = 1.5 y_c$.

3.2 Symmetrical Triangular Cross-Section

For a triangular flume with side slopes $z$ ($1\text{ V}:z\text{ H}$), area is $A = z y^2$ and top width is $T = 2 z y$:

$$\frac{Q^2 (2 z y_c)}{g (z y_c^2)^3} = 1 \implies \frac{2 z Q^2 y_c}{g z^3 y_c^6} = 1 \implies \frac{2 Q^2}{g z^2 y_c^5} = 1 \implies y_c = \left( \frac{2 Q^2}{g z^2} \right)^{1/5}$$

The minimum specific energy in a triangular channel is: $E_{min} = 1.25 y_c$.

3.3 Symmetrical Trapezoidal Cross-Section

For a trapezoidal channel with bottom width $b$ and side slopes $z$, area is $A = (b + z y)y$ and top width is $T = b + 2 z y$. The governing relationship is implicit:

$$\frac{Q^2 (b + 2 z y_c)}{g [(b + z y_c) y_c]^3} = 1$$

This non-linear algebraic equation is solved iteratively using the Newton-Raphson method.

Bélanger Conjugate Depth Equation Derivation

4. Section Factor for Critical Flow and Dimensionless Relationships

Hydraulic engineers define the section factor for critical flow computation $Z$ as:

$$Z = A \sqrt{D_h} = A \sqrt{\frac{A}{T}} = \frac{A^{1.5}}{\sqrt{T}}$$

At the critical flow state, substituting the governing equation gives: $Z_c = \frac{Q}{\sqrt{g}}$. This formulation isolates hydraulic discharge and gravity on the right-hand side, allowing engineers to determine critical depth by matching geometric section factor $Z(y)$ directly to flow demand $Q / \sqrt{g}$.

5. Hydraulic Control Sections: Weirs, Sluice Gates, and Overfalls

A control section is defined as a specific channel location where a unique, stable relationship exists between water depth and discharge. Because critical flow exhibits minimum specific energy, it acts as a physical barrier preventing downstream disturbances from propagating upstream.

                    BROAD-CRESTED WEIR FLOW CONTROL
            Head H
          +---------\                                     Critical Depth yc
          |          \-------------------\              - - - - - - - - -
          |          |   Weir Crest       \------------> Supercritical Flow
          +----------+                     \
                                            \==========> Stilling Basin

1. Broad-Crested Weirs: As water accelerates over a raised weir crest, specific energy constraints force flow through critical depth $y_c$, establishing a dependable stage-discharge rating curve: $Q = C_d b \sqrt{g} \left( \frac{2}{3} H \right)^{1.5}$.
2. Free Overfalls (Brink Depth): At an abrupt drop at the end of a mild channel, flow passes through critical depth shortly upstream of the brink ($x \approx 3.5 y_c$ to $4 y_c$), while brink depth drops to $y_b \approx 0.715 y_c$.
3. Channel Chokes: Channel width constrictions or bottom humps reduce available specific energy. If energy drops below $E_{min}$, flow chokes, creating upstream backwater to increase upstream specific energy.

6. Critical Slope and Channel Classification Criteria

Channel bed slope $S_0$ relative to critical slope $S_c$ governs normal depth $y_n$ and flow regimes:

$$S_c = \frac{g n^2 D_{h,c}}{R_{h,c}^{4/3}}$$

where $n$ is Manning roughness coefficient and $R_{h,c}$ is hydraulic radius at critical depth.

Classification and Stilling Basin Design for Hydraulic Jumps

7. Solved Numerical Example 1: Rectangular Spillway Chute

A concrete rectangular spillway chute has bottom width $b = 6.0\text{ m}$ and carries design discharge $Q = 45.0\text{ m}^3/\text{ s}$. Manning roughness is $n = 0.013$. Calculate critical depth, critical velocity, minimum specific energy, and critical slope using critical flow hydraulics equations.

Compute unit discharge $q$: $q = \frac{Q}{b} = \frac{45.0}{6.0} = 7.50\text{ m}^2/\text{ s}$.

Calculate critical depth $y_c$:

$$y_c = \sqrt[3]{\frac{q^2}{g}} = \sqrt[3]{\frac{(7.50)^2}{9.81}} = \sqrt[3]{\frac{56.25}{9.81}} = \sqrt[3]{5.7339} = 1.7900\text{ m}$$

Calculate critical velocity $V_c$: $V_c = \frac{q}{y_c} = \frac{7.50}{1.7900} = 4.1899\text{ m/s}$.
Calculate minimum specific energy $E_{min}$: $E_{min} = 1.5 y_c = 1.5 \times 1.7900 = 2.6850\text{ m}$.
Critical slope: $A_c = 10.74\text{ m}^2$, $P_c = 9.58\text{ m}$, $R_{h,c} = 1.1211\text{ m}$. $S_c = \left( \frac{0.013 \times 4.1899}{(1.1211)^{2/3}} \right)^2 = 0.002549 \approx 0.255\%$.

8. Solved Numerical Example 2: Trapezoidal Canal Iteration

A trapezoidal channel has bottom width $b = 4.0\text{ m}$ and side slopes $z = 2.0$ ($1\text{ V}:2\text{ H}$). It conveys discharge $Q = 28.0\text{ m}^3/\text{ s}$. Compute critical depth $y_c$ to four decimal places using open channel critical state iteration.

Compute target section factor: $Z_{target} = \frac{Q}{\sqrt{g}} = \frac{28.0}{\sqrt{9.81}} = 8.9397\text{ m}^{2.5}$.

Define residual function: $F(y) = [(4 + 2y)y]^{1.5} – 8.9397 \sqrt{4 + 4y} = 0$.

Evaluating roots gives $y_c = 1.3562\text{ m}$.
Check governing equation: $A_c = [4 + 2(1.3562)](1.3562) = 9.1032\text{ m}^2$, $T_c = 4 + 4(1.3562) = 9.4248\text{ m}$.
$\frac{Q^2 T_c}{g A_c^3} = \frac{(28.0)^2 \times 9.4248}{9.81 \times (9.1032)^3} = \frac{7389.04}{7400.31} \approx 1.0000$.

9. Solved Numerical Example 3: Symmetrical Triangular Storm Drainage Swale

A roadside triangular swale with side slopes $z = 2.5$ ($1\text{ V}:2.5\text{ H}$) carries storm discharge $Q = 6.0\text{ m}^3/\text{ s}$. Compute critical depth and minimum specific energy.

Calculate critical depth using the triangular formula:

$$y_c = \left( \frac{2 Q^2}{g z^2} \right)^{1/5} = \left( \frac{2 \times (6.0)^2}{9.81 \times (2.5)^2} \right)^{1/5} = \left( \frac{72.0}{9.81 \times 6.25} \right)^{1/5} = (1.17431)^{0.2} = 1.0326\text{ m}$$

Calculate parameters: $A_c = 2.5(1.0326)^2 = 2.6657\text{ m}^2$, $T_c = 2(2.5)(1.0326) = 5.1630\text{ m}$. $D_{h,c} = 0.5163\text{ m}$. $V_c = \frac{6.0}{2.6657} = 2.2508\text{ m/s}$. $Fr = \frac{2.2508}{\sqrt{9.81 \times 0.5163}} = 1.0001 \approx 1.0$. Minimum specific energy: $E_{min} = 1.25 y_c = 1.25(1.0326) = 1.2908\text{ m}$.

10. Critical Depth Synthesis

Mastering critical depth open channel flow provides hydraulic engineers with the core analytical foundation needed to design open channel conveyance networks, flow measurement weirs, and energy dissipation structures. By defining the thermodynamic state of minimum specific energy, critical depth establishes the singular transition threshold between subcritical and supercritical regimes. Whether evaluating natural river hydraulics or designing spillway chutes per USBR Design Standards and ASCE Hydraulic Guidelines, critical depth analysis guarantees robust, stable, and predictable fluid management. Calibration standards from NIST Flow Standards ensure metering accuracy.

References & Standards Cited

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

• ASCE: Energy Dissipation in Hydraulic Structures, American Society of Civil Engineers, Reston, VA.

• NIST: Fluid Flow Measurement Standards in Open Channels, 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.

• French, R. H.: Open-Channel Hydraulics, McGraw-Hill, New York.

Frequently Asked Questions (FAQ)

Why is specific energy minimized at critical depth?

Specific energy combines potential head and kinetic velocity head. As depth decreases, potential head drops while velocity head increases exponentially. The derivative $dE/dy = 1 - Fr^2$ equals zero when $Fr = 1$, creating a global minimum energy point.

Can critical depth occur in a uniform channel with constant slope?

Yes. Critical flow occurs uniformly when the channel bed slope equals the critical slope ($S_0 = S_c$). Under this condition, normal depth equals critical depth ($y_n = y_c$).

How does channel width constriction induce critical depth?

When channel width narrows, discharge per unit width ($q = Q/b$) increases, elevating required minimum specific energy. If available energy equals $E_{min}$, flow reaches critical depth at the throat.

What is the physical significance of brink depth at a free overfall?

At a free outfall, water accelerates rapidly toward the drop, causing severe streamline curvature. Non-hydrostatic pressure drops brink depth to approximately $71.5%$ of true critical depth ($y_b approx 0.715 y_c$).

Why is critical flow inherently unstable in physical channels?

Because the $E$-$y$ curve has a horizontal tangent at minimum energy ($dE/dy = 0$), minor changes in channel roughness or bed elevation create large surface waves, resulting in characteristic water surface undulations.