Manning Open Channel Flow & Water Depth Calculator

Principles of Uniform Open Channel Flow

Open channel hydraulics governs the behavior of fluids flowing with a free, unconfined surface exposed to the ambient atmosphere, where movement is driven purely by the longitudinal component of gravity acting down an inclined channel bed. Formulated by Irish civil engineer Robert Manning in 1889, the Manning Equation constitutes the worldwide engineering standard for designing stormwater drainage culverts, municipal flood conveyance channels, agricultural irrigation canals, sanitary sewer interceptors, highway roadside swales, and river engineering works.

A flow is classified as uniform steady flow when water depth, cross-sectional area, velocity, and volumetric discharge remain completely constant along the longitudinal axis of the conduit ($partial y / partial x = 0$ and $partial v / partial x = 0$). Under uniform conditions, the energy slope $S_f$ (friction slope), the hydraulic grade line (water surface slope $S_w$), and the physical bottom bed slope $S_0$ are parallel ($S_f = S_w = S_0$). The depth associated with uniform flow is termed the Normal Depth ($y_n$).

Mathematical Formulations and Derivations

1. The Universal Manning Equation

In metric SI units, the uniform volumetric discharge $Q$ and mean cross-sectional flow velocity $v$ are given by:

where:

• $n$ is the empirical Manning roughness coefficient ($text{s/m}^{1/3}$).
• $A$ is the flow cross-sectional wetted area ($text{m}^2$).
• $P_w$ is the wetted perimeter ($text{m}$), the linear contact length between fluid and solid boundary.
• $R_h = A / P_w$ is the Hydraulic Radius ($text{m}$), representing the ratio of flow volume to boundary frictional resistance.
• $S_0 = tan theta approx sin theta$ is the dimensionless longitudinal bottom bed slope ($text{m/m}$).

2. Geometric Formulations Across Channel Profiles

• Trapezoidal Section: Bottom width $b$, side slopes $z$ horizontal to $1$ vertical:
  A(y) = (b + z y) y, quad P_w(y) = b + 2 y sqrt{1 + z^2}, quad T(y) = b + 2 z y
• Rectangular Section ($z = 0$):
  A(y) = b y, quad P_w(y) = b + 2 y, quad T(y) = b
• Triangular V-Ditch ($b = 0$):
  A(y) = z y^2, quad P_w(y) = 2 y sqrt{1 + z^2}, quad T(y) = 2 z y
• Circular Pipe / Culvert (Diameter $D$): Subtended wetted angle $theta = 2 arccos(1 - 2 y / D)$ (in radians):
  A(y) = frac{D^2}{8} (theta - sin theta), quad P_w(y) = frac{D theta}{2}, quad T(y) = D sin(theta / 2)

3. Solving for Normal Depth ($y_n$)

Because the wetted area $A(y)$ and hydraulic radius $R_h(y)$ depend non-linearly on depth $y$, calculating normal depth $y_n$ for a specified discharge $Q$ requires finding the root of the governing equation:

This calculator executes a high-precision Newton-Raphson numerical algorithm:

4. Critical Depth ($y_c$) and Froude Number ($Fr$)

The Specific Energy ($E$) of an open channel flow is the total mechanical energy head measured relative to the channel bottom invert:

Differentiating with respect to depth $y$ and setting $dE/dy = 0$ reveals that for any specified discharge $Q$, minimum specific energy occurs at the Critical Depth ($y_c$) where:

where $D_h = A / T$ is the hydraulic depth. The dimensionless Froude Number ($Fr$) compares inertial velocity forces to gravitational wave celerity:

• $Fr < 1.0$ (Subcritical Flow): Water is deep and slow ($y_n > y_c$). Gravity wave speed exceeds flow velocity ($c > v$); surface ripples propagate upstream. Control is situated downstream.
• $Fr = 1.0$ (Critical Flow): Flow is at minimum specific energy ($y_n = y_c$). Wave speed matches velocity ($c = v$).
• $Fr > 1.0$ (Supercritical Flow): Water is shallow and rapid ($y_n < y_c$). Flow velocity exceeds wave speed ($v > c$); disturbances cannot propagate upstream. Control is situated upstream.

Step-by-Step Worked Hydraulic Engineering Example

Design a concrete-lined trapezoidal drainage canal conveying a design flood discharge of $Q = 15.0 text{ m}^3/text{s}$ under the following parameters:

• Channel Bottom Width: $b = 3.0 text{ m}$
• Side Slope: $z = 1.5$ (1.5 Horizontal to 1 Vertical, $theta approx 33.7^circ$)
• Bed Slope: $S_0 = 0.0015 text{ m/m}$ ($0.15%$)
• Lining: Smooth finished concrete, $n = 0.013$

Step 1: Compute Manning Channel Section Factor

Step 2: Solve Iteratively for Normal Depth ($y_n$)

We seek $y_n$ such that $A^{5/3} / P_w^{2/3} = 5.0349$:

• Try $y = 1.50 text{ m}$:
  $A = (3.0 + 1.5 times 1.5)(1.5) = (5.25)(1.5) = 7.875 text{ m}^2$
  $P_w = 3.0 + 2(1.5)sqrt{1 + 1.5^2} = 3.0 + 3.0(1.8028) = 8.408 text{ m}$
  $R_h = 7.875 / 8.408 = 0.9366 text{ m}$
  $A R_h^{2/3} = 7.875 times (0.9366)^{2/3} = 7.875 times 0.9572 = 7.538$
  $Q = frac{1}{0.013} times 7.538 times sqrt{0.0015} = 76.923 times 7.538 times 0.03873 = 22.45 text{ m}^3/text{s}$ (Too high!)
• Try $y = 1.20 text{ m}$:
  $A = (3.0 + 1.5 times 1.2)(1.2) = 4.80 times 1.2 = 5.76 text{ m}^2$
  $P_w = 3.0 + 2.4(1.8028) = 7.327 text{ m} implies R_h = 0.7861 text{ m}$
  $Q = 76.923 times 5.76 times (0.7861)^{2/3} times 0.03873 = 14.61 text{ m}^3/text{s}$ (Slightly low)
• Exact Newton-Raphson solution converges in 4 iterations to:
  y_n = 1.221 text{ meters}

Step 3: Hydraulic Properties at Normal Depth

Step 4: Critical Depth ($y_c$) and Froude Number Calculation

Solving $Q^2 T / (g A^3) = 1$ for critical depth yields:

Evaluating the Froude number at normal depth:

Because $y_n = 1.221 text{ m} > y_c = 1.182 text{ m}$ and $Fr = 0.863 < 1.0$, the canal operates in the subcritical tranquil flow regime on a mild bed slope ($S_0 < S_c$), ensuring hydraulic stability without premature hydraulic jumps.

Manning Roughness ($n$) Engineering Reference Table