Manning's Open Channel & Pipe Flow Velocity Calculator

1. Gauckler-Manning-Strickler Formula for Open Channel Flow

In civil engineering hydraulics, drainage design, and irrigation engineering, open channel flow occurs whenever liquid travels with a free surface exposed to atmospheric pressure. The Gauckler-Manning formula (developed by Philippe Gauckler in 1867 and popularized by Robert Manning in 1889) is the globally accepted empirical standard for calculating uniform steady gravity flow:

Where $V$ is mean cross-sectional flow velocity, $n$ is the dimensionless Manning roughness coefficient, $R_h$ is the hydraulic radius ($m$), $S_0$ is the longitudinal channel bottom slope ($m/m$), and $k = 1.0$ for SI metric units ($k = 1.4859$ for US customary units).

Total volumetric flow rate (discharge $Q$) is obtained from the continuity equation:

2. Geometric Formulations by Channel Cross-Section

1. Trapezoidal & Rectangular Channels:

For a symmetrical trapezoidal canal with bottom width $B$, water depth $y$, and side slopes $z:1$ (horizontal to vertical run):

• Cross-Sectional Area: $A = (B + z y) y$
• Wetted Perimeter: $P = B + 2 y sqrt{1 + z^2}$
• Top Surface Width: $T = B + 2 z y$
• Hydraulic Radius: $R_h = frac{A}{P} = frac{(B + z y) y}{B + 2 y sqrt{1 + z^2}}$

For a rectangular channel, simply set $z = 0$, giving $A = B y$ and $P = B + 2 y$.

2. Partially Filled Circular Pipes:

For a circular sewer or culvert barrel of internal diameter $D$ filled to depth $y le D$, defining the central wetted angle $theta$ in radians:

• Cross-Sectional Area: $A = frac{D^2}{8} (theta - sintheta)$
• Wetted Perimeter: $P = frac{D theta}{2}$
• Top Surface Width: $T = D sinleft(frac{theta}{2}right) = 2 sqrt{y (D - y)}$
• Hydraulic Radius: $R_h = frac{D}{4} left( 1 - frac{sintheta}{theta} right)$

3. The Circular Pipe Hydraulic Anomaly ($y/D = 0.938$)

A unique mathematical feature of circular gravity conduits is that maximum discharge does not occur when the pipe is 100% full:

• At $y/D = 0.81$, velocity reaches its maximum: $V_{max} = 1.14 times V_{text{full}}$.
• At $y/D = 0.938$, discharge reaches its maximum: $Q_{max} = 1.076 times Q_{text{full}}$.

Between $y/D = 0.938$ and $1.0$, as water touches the upper pipe crown, the wetted perimeter rapidly increases from $sim 85%$ to $100%$ without adding meaningful cross-sectional area. The added surface boundary drag causes flow to decelerate.

4. Froude Number and Critical Flow Dynamics

The hydraulic depth is defined as $D_h = A / T$. The Froude number of the open channel flow is:

• $Fr < 1.0$ (Subcritical Flow): Tranquil, deep flow. Waves travel upstream; water surface profiles are controlled from downstream boundaries. Desired for municipal canals and navigation.
• $Fr = 1.0$ (Critical Flow): The flow possesses minimum specific energy. Critical slope $S_c$ is the bed slope that produces $yn = yc$.
• $Fr > 1.0$ (Supercritical Flow): Rapid, shooting flow. Highly susceptible to standing waves, roll waves, and sudden hydraulic jumps. Designers limit highway chutes to $Fr < 0.8$ or $Fr > 1.2$ to avoid unstable pulsating transitions.

5. Step-by-Step Worked Numerical Example: Concrete Irrigation Canal

Consider a trapezoidal concrete canal carrying irrigation water:

• Bottom width $B = 3.0text{ m}$, side slope $z = 1.5$ ($1.5text{H}:1text{V}$)
• Normal water depth $y = 1.20text{ m}$
• Longitudinal bed slope $S_0 = 0.0015$ ($0.15%$)
• Manning's roughness $n = 0.014$ (smooth finished concrete)

1. Calculate Cross-Sectional Flow Area ($A$): $$A = (B + z y) y = (3.0 + 1.5 times 1.20) times 1.20 = (3.0 + 1.80) times 1.20 = 4.80 times 1.20 = 5.76text{ m}^2$$
2. Calculate Wetted Perimeter ($P$): $$P = B + 2 y sqrt{1 + z^2} = 3.0 + 2 times 1.20 times sqrt{1 + (1.5)^2} = 3.0 + 2.40 times sqrt{1 + 2.25} = 3.0 + 2.40 times 1.8028 approx 7.327text{ m}$$
3. Calculate Hydraulic Radius ($R_h$): $$R_h = frac{A}{P} = frac{5.76text{ m}^2}{7.327text{ m}} approx 0.7861text{ m}$$ $$R_h^{2/3} = (0.7861)^{0.6667} approx 0.8517text{ m}^{2/3}$$
4. Calculate Flow Velocity ($V$): $$V = frac{1.0}{0.014} times 0.8517 times sqrt{0.0015} = 71.4286 times 0.8517 times 0.03873 approx 2.356text{ m/s} quad (7.73text{ ft/s})$$
5. Calculate Volumetric Discharge ($Q$): $$Q = A cdot V = 5.76text{ m}^2 times 2.356text{ m/s} approx 13.57text{ m}^3/text{s} quad (13{,}570text{ L/s} approx 479text{ cfs})$$
6. Calculate Top Width and Froude Number ($Fr$): $$T = B + 2 z y = 3.0 + 2 times 1.5 times 1.20 = 3.0 + 3.60 = 6.60text{ m}$$ $$D_h = frac{A}{T} = frac{5.76}{6.60} approx 0.8727text{ m}$$ $$Fr = frac{V}{sqrt{g D_h}} = frac{2.356}{sqrt{9.80665 times 0.8727}} = frac{2.356}{sqrt{8.558}} = frac{2.356}{2.925} approx 0.805 quad (< 1.0text{ Subcritical})$$

6. Municipal Storm Sewer & Channel Design Guidelines

• Self-Cleansing Velocity: Gravity sewers must achieve a minimum velocity of $0.6–0.9text{ m/s}$ ($2.0–3.0text{ ft/s}$) at least once daily to prevent heavy grit and organic solids from settling on the pipe invert.
• Maximum Non-Erosive Velocity: For unlined earthen canals, limit velocities to $< 0.8 - 1.2text{ m/s}$ to prevent bed scour. For reinforced concrete linings, velocities up to $3.5–4.5text{ m/s}$ are permissible.
• Freeboard Allowance: Always provide adequate freeboard height above maximum design water stage ($F_b approx 0.5 + 0.25 y$ for earthen canals, or at least $0.3text{ m}$ for concrete channels) to accommodate wind-generated wave run-up and sediment accumulation.

7. The Best Hydraulic Section (Most Economical Geometry)

In open channel design, excavation and concrete lining costs are proportional to the wetted perimeter $P$. A channel cross-section that conveys maximum discharge $Q$ for a given cross-sectional area $A$ (or equivalently, minimizes wetted perimeter $P$ for a given area $A$) is known as the Best Hydraulic Section:

• Semi-Circular Channel: Mathematically the most efficient of all possible open cross-sections. Hydraulic radius equals half the radius: $R_h = r / 2 = y / 2$. However, constructing circular curves in earth or concrete is often labor-intensive.
• Trapezoidal Channel (Semi-Hexagon): The most efficient trapezoidal section forms half of a regular hexagon, where side slopes are $z = 1/sqrt{3} approx 0.5774$ ($60^circ$ angle with horizontal), bottom width is $B = frac{2}{sqrt{3}} y approx 1.155 y$, and hydraulic radius is exactly $R_h = y / 2$.
• Rectangular Channel: The best rectangular section is a half-square where width is exactly twice the depth: $B = 2 y$, and $R_h = y / 2$.

8. Composite Roughness in Compound Channels (Horton & Pavlovskii Methods)

Real channels often have composite boundaries (e.g., concrete bed with riprap banks, or a main river channel with a heavily vegetated overbank floodplain). To apply Manning's equation, an equivalent composite roughness $n_e$ must be determined:

• Horton-Einstein Formulation: Assumes mean velocity in all sub-areas is identical: $$n_e = left( frac{sum_{i=1}^N P_i n_i^{3/2}}{P} right)^{2/3}$$
• Pavlovskii Formulation: Assumes shear forces are additive: $$n_e = left( frac{sum_{i=1}^N P_i n_i^2}{P} right)^{1/2}$$

For compound river floodplains where depths on the floodplain are much shallower than the main incised channel, dividing the channel into separate parallel sub-sections and summing their discharges ($Q_{text{total}} = Q_{text{main}} + Q_{text{left}} + Q_{text{right}}$) prevents severe underestimation of flood carrying capacity.