7.1 Manning's Equation & Open Channel Flow Hydraulics

7.1 Manning's Equation & Open Channel Flow Hydraulics

Quick Reference: Open channel flow is liquid motion driven by gravity with a continuous free surface exposed to atmospheric pressure. Unlike closed conduit pressurized pipe flow, open channel flow depth and velocity vary with cross-sectional geometry, bed slope, and boundary friction. Under steady, uniform flow, the channel bed ($S_0$), the water surface or Hydraulic Grade Line (HGL), and the total Energy Grade Line (EGL) are strictly parallel ($S_0 = S_w = S_f = S$). In US Customary engineering, Manning's Equation ($V = \frac{1.486}{n} R^{2/3} S^{1/2}$) is the standard predictive tool for velocity and discharge capacity. Designing temporary swales and ditches requires evaluating flow regimes via the Froude Number ($Fr = \frac{V}{\sqrt{g D}}$) and sizing channel linings using the maximum permissible shear stress method ($\tau_0 = \gamma R S$) rather than legacy velocity limits alone.

---

Governing Physics of Open Channel Flow

In civil engineering and environmental hydraulics, fluid conveyance is divided into two broad categories: pipe flow (where water fills a closed conduit under positive hydrostatic pressure) and open channel flow (where the moving liquid maintains a free surface exposed directly to atmospheric pressure). On construction sites and land development projects, runoff conveyance systems—including diversion ditches, roadside swales, perimeter interceptors, temporary bypass channels, and natural streams—function almost exclusively as open channels.

Forces Governing Open Channel Flow

Open channel flow is driven by the tangential component of gravity acting along the longitudinal slope of the channel bed. For a discrete control volume of water of length $L$, cross-sectional area $A$, and unit weight $\gamma$, the gravitational driving force ($F_g$) parallel to a channel bed inclined at angle $\theta$ is:

$F_g = \gamma \times A \times L \times \sin \theta$

Because channel slopes in civil engineering are typically gentle ($\theta < 10^\circ$), the approximation $\sin \theta \approx \tan \theta = S_0$ (where $S_0$ is the bed slope in $\text{ft/ft}$) is universally applied:

$F_g \approx \gamma \times A \times L \times S_0$

Resisting this gravitational force is the boundary shear resistance ($F_r$) exerted by the channel bed and banks over the wetted perimeter ($P$):

$F_r = \tau_0 \times P \times L$

Where:

• $\tau_0$ = Mean boundary shear stress (tractive force per unit wetted area, in $\text{lbs/ft}^2$ or psf)
• $P$ = Wetted perimeter of the cross-section (ft)
• $L$ = Length of the channel reach (ft)

Steady Uniform Flow vs. Varied Flow

• Steady Flow: Hydraulic parameters at any fixed cross-section (depth $y$, velocity $V$, discharge $Q$) do not change over time ($\partial y / \partial t = 0$, $\partial V / \partial t = 0$).
• Uniform Flow: Hydraulic parameters do not change along the longitudinal length of the channel reach ($\partial y / \partial x = 0$, $\partial V / \partial x = 0$). The depth of water in uniform flow is termed the normal depth ($y_n$).

When flow is both steady and uniform, the gravitational driving force precisely balances the boundary shear resistance ($F_g = F_r$):

$\gamma \times A \times L \times S_0 = \tau_0 \times P \times L$

$\tau_0 = \gamma \times \left( \frac{A}{P} \right) \times S_0 = \gamma \times R \times S_0$

This fundamental equilibrium equation confirms that boundary shear stress is directly proportional to fluid density, channel slope, and the ratio of flow area to wetted perimeter ($R = A / P$).

The Energy Grade Line (EGL) & Hydraulic Grade Line (HGL)

To analyze energy distribution along a channel, engineers apply the Bernoulli Energy Equation. The total hydraulic energy head ($H$) at any channel cross-section, measured relative to an arbitrary datum, is expressed in feet of fluid:

$H = z + y + \frac{V^2}{2g}$

Where:

• $z$ = Channel invert bed elevation above datum (ft)

• $y$ = Flow depth normal to the bed (ft), representing potential pressure head

• $\frac{V^2}{2g}$ = Velocity head (dynamic kinetic energy head, ft), where $g = 32.2\text{ ft/s}^2$ and $V$ is mean velocity in ft/s

• Hydraulic Grade Line (HGL): Represents the piezometric head ($z + y$). In open channel flow, because the water surface is open to the atmosphere (gauge pressure $p = 0$), the HGL is identical to the free water surface profile.

• Energy Grade Line (EGL): Positioned at a vertical distance of $\frac{V^2}{2g}$ above the HGL, representing the total mechanical energy available per unit weight of fluid.

In steady, uniform flow, the velocity head remains constant along the reach. Consequently, the channel bed slope ($S_0$), the water surface slope ($S_w$ or HGL slope), and the energy friction slope ($S_f$ or EGL slope) are perfectly parallel:

$S_0 = S_w = S_f = S$

When flow is non-uniform (such as backwater curves approaching a constricted culvert or drawdown curves approaching an abrupt drop), $S_0 \neq S_w \neq S_f$, requiring gradually varied flow (GVF) standard step routing.

---

Manning's Equation in US Customary Units

First formulated by the Irish engineer Robert Manning in 1889 and subsequently adapted by civil engineers worldwide, Manning's Equation is the premier empirical relationship for evaluating uniform, steady open channel flow.

The Governing Formulations

In US Customary units (English Engineering system), Manning's Equation for mean cross-sectional velocity is:

$V = \frac{1.486}{n} \times R^{2/3} \times S^{1/2}$

By multiplying mean velocity ($V$) by the cross-sectional flow area ($A$), the volumetric discharge capacity ($Q$) is obtained:

$Q = V \times A = \frac{1.486}{n} \times A \times R^{2/3} \times S^{1/2}$

Where:

• $V$ = Mean cross-sectional flow velocity ($\text{ft/s}$)
• $Q$ = Volumetric flow rate or discharge capacity ($\text{ft}^3\text{/s}$ or $\text{cfs}$)
• $n$ = Manning's roughness coefficient (dimensionless empirical parameter)
• $A$ = Cross-sectional area of flow ($\text{ft}^2$)
• $P$ = Wetted perimeter of the cross-section ($\text{ft}$)
• $R$ = Hydraulic radius ($\text{ft}$), defined as $R = A / P$
• $S$ = Longitudinal channel energy gradient or bed slope ($\text{ft/ft}$, dimensionless decimal format)

Derivation of the 1.486 Conversion Factor

In the International System of Units (SI), Manning's equation is written without a numerical coefficient:

$V = \frac{1.0}{n} \times R^{2/3} \times S^{1/2} \quad \text{[SI Metric, with } V \text{ in m/s, } R \text{ in m]}$

To preserve the exact same numerical values for the roughness coefficient $n$ across both metric and English systems, an empirical conversion factor is required for US Customary units. Because $1\text{ meter} = 3.28084\text{ feet}$, converting the length term $R^{2/3}$ from meters to feet requires multiplying by:

$(3.28084)^{1/3} \approx 1.48599 \approx 1.486$

Thus, the constant $1.486$ is purely a dimensional conversion multiplier allowing engineers to use standardized $n$-values while working in feet and seconds.

---

Geometric Elements of Open Channels

Accurately computing flow capacity requires breaking down the cross-sectional geometry into five fundamental hydraulic parameters:

1. Cross-Sectional Flow Area ($A$, $\text{ft}^2$): The plane area occupied by water, measured perpendicular to the direction of flow.
2. Wetted Perimeter ($P$, $\text{ft}$): The total linear distance of the channel boundary (bed and side slopes) in direct contact with the flowing water. Critically, the top width of the air-water free surface is excluded because air provides negligible shear drag compared to the solid boundary.
3. Hydraulic Radius ($R$, $\text{ft}$): The ratio of flow area to wetted perimeter ($R = A / P$). The hydraulic radius serves as an index of hydraulic efficiency. A channel shape that minimizes wetted perimeter relative to area reduces frictional drag, thereby maximizing mean velocity and discharge.
4. Top Width ($T$, $\text{ft}$): The linear horizontal width of the water surface across the channel.
5. Hydraulic Depth ($D$, $\text{ft}$): The ratio of flow area to top width ($D = A / T$). The hydraulic depth represents the average depth of the water column and is the governing characteristic length used in wave speed and Froude number calculations.

Hydraulically Optimal Sections

A cross-section that conveys the maximum possible discharge ($Q$) for a given cross-sectional area ($A$), slope ($S$), and roughness ($n$) is termed the hydraulically optimal (or most efficient) cross-section. Mathematically, this corresponds to the shape with the minimum wetted perimeter ($P$), which maximizes the hydraulic radius ($R$):

• Theoretical Optimum: A semicircular channel cross-section has the lowest possible perimeter-to-area ratio. However, semicircular shapes are practically impossible to excavate and stabilize in unreinforced earthen soils.
• Trapezoidal Optimum: A half-hexagon profile with side slopes angled at $60^\circ$ to the horizontal ($z = 1/\sqrt{3} \approx 0.577:1$) and bottom width $b = 2 y / \sqrt{3}$. For this shape, $R = y / 2$.
• Rectangular Optimum: A rectangle whose bottom width is exactly twice its flow depth ($b = 2y$), yielding $R = y / 2$.
• Triangular Optimum: A $90^\circ$ V-notch with $1:1$ side slopes ($z = 1.0$), yielding $R = y / (2 \sqrt{2}) \approx 0.354 y$.

Channel Geometry Reference Formulas

The following table summarizes the standard analytical expressions used in CPESC hydraulic modeling across common geometric cross-sections:

Note on Trapezoidal Notation: In civil engineering convention, side slopes are expressed as horizontal-to-vertical ratio ($z:1$, meaning $z$ units horizontal to 1 unit vertical). For example, a $2:1$ side slope has $z = 2.0$.

---

Manning's Roughness Coefficient ($n$)

The roughness coefficient $n$ represents the total integrated resistance offered by the channel bed and banks to fluid motion. Manning's $n$ is not purely a material friction constant; it is influenced by multiple hydraulic and physical mechanisms.

Physical Determinants of Manning's $n$

1. Surface Roughness & Grain Size: The physical size and angularity of bed material grains ($d_{50}$). Frictional resistance increases as grain diameter increases.
2. Channel Cross-Section Irregularity: Variations in channel shape, bottom dips, and gouges increase turbulence, raising the effective $n$-value.
3. Vegetation Density & Retardance: Flexible turfgrass, emergent sedges, and rigid woody vegetation exert substantial drag. In vegetated swales, the roughness coefficient varies dynamically with flow depth: when water is shallow and grass stands erect, $n$ is exceptionally high ($n > 0.08\text{–}0.12$); as flow depth and velocity increase, grass bends and submerges, flattening against the bed and reducing effective roughness ($n \approx 0.035\text{–}0.045$).
4. Channel Meandering (Sinuosity): Curves and sharp bends create secondary spiral flow cells that consume hydraulic energy, increasing effective $n$ by $5%\text{ to }30%$.
5. Obstructions & Sediment Deposition: Boulders, fallen debris, roots, and localized sediment bars obstruct streamlines and increase turbulent losses.

Cowan's Formulative Method for Composite Roughness

Where complex natural channels must be evaluated, Cowan (1956) formulated a systematic additive procedure adopted by the USGS and NRCS:

$n = (n_0 + n_1 + n_2 + n_3 + n_4) \times m_5$

Where:

• $n_0$ = Base material roughness (bare soil = $0.020$, rock = $0.025$, fine gravel = $0.024$)
• $n_1$ = Degree of surface irregularity ($0.000\text{ to }0.020$)
• $n_2$ = Variations in channel cross-section ($0.000\text{ to }0.015$)
• $n_3$ = Effect of obstructions ($0.000\text{ to }0.060$)
• $n_4$ = Vegetative retardance class ($0.005\text{ to }0.100$)
• $m_5$ = Meandering correction factor ($1.00\text{ for straight channels up to }1.30\text{ for severe meandering}$)

Design Values of Manning's $n$ for CPESC Applications

---

Flow Regimes & The Froude Number ($Fr$)

The Froude Number ($Fr$) is a dimensionless parameter that represents the ratio of inertial forces to gravitational forces within an open channel flow field. It governs the hydraulic state of the flow, dictates how surface disturbances propagate, and indicates whether catastrophic hydraulic jumps may occur.

Mathematical Definition

$Fr = \frac{V}{\sqrt{g \times D}}$

Where:

• $V$ = Mean cross-sectional flow velocity ($\text{ft/s}$)
• $g$ = Gravitational acceleration constant ($32.2\text{ ft/s}^2$)
• $D$ = Hydraulic depth ($\text{ft}$), defined as cross-sectional area divided by top width ($D = A / T$)
• $\sqrt{g D} = c$ = Wave celerity ($\text{ft/s}$), which is the speed at which a shallow gravity wave or ripple travels relative to the fluid

                    ┌─── Fr < 1.0 : Subcritical (Tranquil, deep, slow; V < c)
                    │
Froude Number (Fr) ─┼─── Fr = 1.0 : Critical (Minimum energy head, unstable surface)
                    │
                    └─── Fr > 1.0 : Supercritical (Rapid, shallow, shooting; V > c)

The Three Hydraulic Flow Regimes

1. Subcritical Flow ($Fr < 1.0$):

   • Physics: Inertial forces are subordinate to gravitational forces ($V < c$). Because fluid velocity is slower than wave celerity, surface waves and disturbances can travel upstream against the current.
   • Operational Behavior: Flow is deep, tranquil, and slow. The depth of flow at any point is governed by downstream controls (such as a culvert inlet, weir, drop structure, or receiving reservoir backwater).
   • Design Objective: Civil engineers design the vast majority of temporary diversion swales, roadside ditches, and grassed waterways to operate in the subcritical regime ($Fr \le 0.80$) to avoid standing waves, high-velocity liner detachment, and uncontrolled bed erosion.

2. Critical Flow ($Fr = 1.0$):

   • Physics: Flow velocity exactly matches the wave celerity ($V = c = \sqrt{g D}$). Surface disturbances remain stationary, forming standing waves.
   • Energy Relationship: For a specified volumetric discharge ($Q$), critical flow corresponds to the absolute minimum specific energy head ($E = y + V^2 / 2g$) capable of conveying that flow. Conversely, for a given specific energy head, critical flow conveys the maximum possible discharge.
   • Critical Depth ($y_c$): The flow depth at $Fr = 1.0$. In a rectangular channel, critical depth simplifies to $y_c = \sqrt[3]{q^2 / g}$, where $q = Q / b$ is unit discharge. Critical flow is inherently unstable; minor bed irregularities cause severe water surface undulations.

3. Supercritical Flow ($Fr > 1.0$):

   • Physics: Inertial forces dominate gravitational forces ($V > c$). Fluid moves faster than surface gravity waves, meaning disturbances cannot propagate upstream. Flow depth is governed entirely by upstream controls (such as a steep chute crest or sluice opening).
   • Operational Behavior: Flow is shallow, rapid, and shooting. Common on steep cut slopes, riprap chutes, and paved flumes.
   • Construction Hazards: Supercritical flow exerts immense tractive shear stress on channel liners. Furthermore, whenever supercritical flow encounters a downstream obstacle, a flatter slope, or an expanding cross-section, it cannot transition smoothly back to subcritical flow. Instead, it undergoes an abrupt, violent hydraulic jump.

The Hydraulic Jump Hazard

A hydraulic jump is a rapid, turbulent open channel phenomenon where supercritical flow ($Fr_1 > 1.0$) spontaneously transitions into subcritical flow ($Fr_2 < 1.0$). This transition involves massive internal turbulence, intense roller waves, air entrainment, and significant energy dissipation. On construction sites, an uncontrolled hydraulic jump inside an unlined or turf-lined swale creates severe standing waves that can overtop diversion dikes, gouge out the channel banks, and instantly tear rolled erosion control blankets from their anchoring staples. Downstream of steep chutes, an engineered energy dissipator (such as a riprap plunge pool or baffled stilling basin) must be provided to force and contain the hydraulic jump within an armored boundary.

---

Permissible Velocity vs. Maximum Permissible Shear Stress Method

Historically, civil engineers designed drainage channels using the permissible velocity method. In modern CPESC practice, this empirical approach has been superseded by the maximum permissible shear stress (tractive force) method established in FHWA Hydraulic Engineering Circular No. 15 (HEC-15).

Limitations of the Permissible Velocity Method

The permissible velocity method specifies a single maximum allowable mean velocity ($V_{allow}$) for a given soil type or liner material (for example, stating that sandy loam channels should not exceed $2.0\text{ ft/s}$). While conceptually straightforward, this method suffers from a major physical deficiency: it fails to account for flow depth.

Hydrodynamic boundary drag is not dictated by mean velocity alone. A deep channel flowing at $3.0\text{ ft/s}$ exerts far greater erosive shear force on its bed than a very shallow sheet flow moving at $4.0\text{ ft/s}$. Because the permissible velocity method treats all flow depths identically, it dangerously underestimates erosive forces in deep flows and over-designs shallow swales.

The Maximum Permissible Shear Stress (Tractive Force) Method

Under tractive force theory, channel stability is governed by the hydrodynamic drag force exerted by the moving fluid on the channel perimeter per unit surface area. The channel liner is deemed stable if the maximum shear stress generated by the design storm ($\tau_{max}$) is less than or equal to the permissible shear stress ($\tau_{allow}$) of the soil or lining:

$\tau_{max} \le \tau_{allow}$

Boundary Shear Stress Formulations

For a wide channel under steady uniform flow, the mean boundary shear stress ($\tau_0$) is:

$\tau_0 = \gamma \times R \times S$

Where:

• $\tau_0$ = Mean boundary shear stress ($\text{lbs/ft}^2$ or psf)
• $\gamma$ = Unit weight of water ($62.4\text{ lbs/ft}^3$)
• $R$ = Hydraulic radius ($\text{ft}$)
• $S$ = Longitudinal energy slope or bed slope ($\text{ft/ft}$)

In trapezoidal and rectangular channels, boundary shear stress is not distributed uniformly across the wetted perimeter. Shear stress peaks at the center of the bed and decreases along the side slopes:

• Maximum Bed Shear Stress ($\tau_{bed}$): $\tau_{bed} = \gamma \times y \times S$ (Where $y$ is maximum flow depth in feet)

• Maximum Side Slope Shear Stress ($\tau_{side}$): $\tau_{side} = K_{side} \times \gamma \times y \times S$ (Where $K_{side}$ is a side shear factor typically ranging from $0.75\text{ to }0.78$ for standard trapezoidal channels with $z \ge 1.5$)

Because soil particles on side slopes are also subjected to gravitational downhill forces, the allowable shear stress on side slopes is lower than on the channel bed. The side slope safety factor must be evaluated separately from the channel bed.

---

Comprehensive Worked Engineering Calculation

Problem Statement

A CPESC professional is designing a temporary trapezoidal drainage swale to convey stormwater runoff from an active 18-acre mass grading site. The channel must be sized for the 10-year design storm discharge of $Q_{10} = 28.5\text{ cfs}$.

Preliminary site geometry establishes:

• Bottom width ($b$): $3.0\text{ ft}$
• Side slopes ($z$): $2:1$ ($z = 2.0\text{ ft horizontal per 1.0 ft vertical}$)
• Longitudinal bed slope ($S$): $0.015\text{ ft/ft}$ ($1.5%\text{ grade}$)
• Specified channel liner: Advanced Turf Reinforcement Mat (TRM) with an unvegetated allowable shear stress of $\tau_{allow} = 2.50\text{ psf}$ and a design Manning's roughness of $n = 0.035$

The engineer must evaluate a trial flow depth of $y = 1.25\text{ ft}$ to verify:

1. Cross-sectional Area ($A$), Wetted Perimeter ($P$), and Hydraulic Radius ($R$);
2. Mean flow velocity ($V$) and volumetric flow capacity ($Q$);
3. Top width ($T$), hydraulic depth ($D$), and Froude Number ($Fr$);
4. Maximum bed shear stress ($\tau_{bed}$) and maximum side shear stress ($\tau_{side}$) compared to $\tau_{allow}$.

Step 1: Geometric Elements Calculation

$\text{Area: } A = (b + z y) y = (3.0\text{ ft} + (2.0 \times 1.25\text{ ft})) \times 1.25\text{ ft} = (3.0 + 2.50) \times 1.25 = 5.50 \times 1.25 = 6.875\text{ ft}^2$

$\text{Wetted Perimeter: } P = b + 2 y \sqrt{1 + z^2} = 3.0 + (2 \times 1.25) \sqrt{1 + 2.0^2} = 3.0 + 2.50 \sqrt{5} = 3.0 + 2.50(2.23607) = 3.0 + 5.5902 = 8.590\text{ ft}$

$\text{Hydraulic Radius: } R = \frac{A}{P} = \frac{6.875\text{ ft}^2}{8.590\text{ ft}} = 0.8003\text{ ft}$

Step 2: Velocity & Discharge Capacity via Manning's Equation

$R^{2/3} = (0.8003)^{0.6667} = 0.8619$

$S^{1/2} = (0.015)^{0.5} = 0.12247$

$V = \frac{1.486}{n} \times R^{2/3} \times S^{1/2} = \frac{1.486}{0.035} \times 0.8619 \times 0.12247$

$V = 42.4571 \times 0.8619 \times 0.12247 = 4.482\text{ ft/s} \approx 4.48\text{ ft/s}$

$Q = V \times A = 4.482\text{ ft/s} \times 6.875\text{ ft}^2 = 30.81\text{ cfs} \approx 30.8\text{ cfs}$

Hydraulic Capacity Evaluation: Calculated capacity ($30.8\text{ cfs}$) exceeds the design requirement ($Q_{10} = 28.5\text{ cfs}$), providing adequate conveyance capacity.

Step 3: Flow Regime & Froude Number Evaluation

$\text{Top Width: } T = b + 2 z y = 3.0 + (2 \times 2.0 \times 1.25) = 3.0 + 5.0 = 8.0\text{ ft}$

$\text{Hydraulic Depth: } D = \frac{A}{T} = \frac{6.875\text{ ft}^2}{8.0\text{ ft}} = 0.8594\text{ ft}$

$\text{Wave Celerity: } c = \sqrt{g \times D} = \sqrt{32.2\text{ ft/s}^2 \times 0.8594\text{ ft}} = \sqrt{27.673} = 5.260\text{ ft/s}$

$Fr = \frac{V}{c} = \frac{4.482\text{ ft/s}}{5.260\text{ ft/s}} = 0.852$

Flow Regime Evaluation: Because $Fr = 0.852 < 1.0$, flow is in the subcritical regime. Disturbances can travel upstream, and the channel will not generate spontaneous destructive hydraulic jumps during uniform flow conditions.

Step 4: Tractive Force & Shear Stress Stability Analysis

$\tau_{bed} = \gamma \times y \times S = 62.4\text{ lbs/ft}^3 \times 1.25\text{ ft} \times 0.015\text{ ft/ft} = 1.170\text{ lbs/ft}^2 \text{ (psf)}$

$\tau_{side} \approx 0.76 \times \tau_{bed} = 0.76 \times 1.170\text{ psf} = 0.889\text{ lbs/ft}^2 \text{ (psf)}$

$\text{Factor of Safety (Bed): } FS_{bed} = \frac{\tau_{allow}}{\tau_{bed}} = \frac{2.50\text{ psf}}{1.170\text{ psf}} = 2.14$

Design Verification: The maximum bed shear stress of $1.17\text{ psf}$ is well below the TRM allowable limit of $2.50\text{ psf}$ ($FS = 2.14 > 1.30$). The channel configuration and lining are hydraulically and structurally stable.