Open Channel Flow — Manning’s & Chezy’s Equation

1. Introduction — Open Channel vs Pipe Flow

An open channel is any conduit in which a liquid flows with a free surface — a surface in contact with the atmosphere at constant (atmospheric) pressure. Rivers, irrigation canals, drainage ditches, roadside gutters, partially filled sewers, and spillway chutes are all open channels. The defining feature is that the driving force for flow is gravity acting on the sloping liquid surface (or, equivalently, the component of gravity along the channel bed), not an externally applied pressure gradient.

1.1 Types of Open Channel Flow

2. Geometric Properties of Channel Sections

Before applying Manning’s or Chezy’s equation, the geometric properties of the channel cross-section must be computed.

2.1 Formulas for Common Cross-Sections

In the circle formulas, θ = central angle (radians) subtended by the water surface at the centre; θ = 2 cos⁻¹(1 – 2y/D). For a half-full circle: θ = π, A = πD²/8, P = πD/2, R = D/4.

3. Chezy’s Formula

Chezy’s formula (1769) is the earliest empirical formula for uniform open channel flow. It was developed by Antoine Chezy for designing the Paris water supply canals.

3.1 Empirical Formulas for Chezy’s C

4. Manning’s Equation

Manning’s equation (1889), proposed by Irish engineer Robert Manning, is the standard formula for uniform open channel flow in civil engineering practice worldwide — and the primary formula for GATE CE.

4.1 Derivation Link to Chezy’s Formula

4.2 Conditions for Validity of Manning’s Equation

• Steady, uniform flow (depth and velocity constant along the channel)
• Turbulent flow — essentially always satisfied for natural channels and designed canals
• Fully rough turbulent regime — Manning’s n is independent of velocity (unlike laminar flow where f = 64/Re)
• Channel slope not too steep (S < 0.1) — for steep channels, corrections may be needed

5. Manning’s n — Values & Selection

GATE CE standard values (memorise): Concrete pipe/channel: n = 0.013; Earthen canal: n = 0.025. These are the values given or assumed in most GATE numerical problems.

6. Discharge Formulas for Common Sections

6.1 Rectangular Channel

6.2 Trapezoidal Channel (most common in practice)

6.3 Circular Pipe (Partial Flow)

Circular pipes running partially full are common in sewer and drainage design. Key result:

7. Most Economical Channel Section

The most economical (best hydraulic or most efficient) channel section is the one that conveys maximum discharge for a given cross-sectional area A, Manning’s n, and slope S. From Manning’s equation Q = (1/n)A R^2/3 S^1/2, maximising Q for fixed A means maximising R = A/P — which means minimising the wetted perimeter P. This is equivalent to minimising excavation and lining cost for a given conveyance capacity.

7.1 Most Economical Rectangular Section

7.2 Most Economical Trapezoidal Section

7.3 Best Hydraulic Section — Summary

Practical note: The semi-circle is the theoretical optimum, but trapezoidal channels are used in practice because they are easier to excavate, more stable, and more suitable for earthen construction. The half-hexagon trapezoidal section is the standard design for irrigation canals.

8. Froude Number & Flow Classification

8.1 Critical Depth for Rectangular Channel

8.2 Critical Depth for Trapezoidal Channel

9. Gradually Varied Flow (GVF) Profiles

Gradually varied flow (GVF) occurs when the depth changes slowly along the channel due to a change in bed slope, cross-section, or a downstream control (dam, weir, change in slope). The governing equation is:

9.1 GVF Profile Classification

Profiles are named by the channel slope type (M, S, C, H, A) and zone (1, 2, 3). The slope type compares the actual bed slope S₀ to the critical slope S_c (the slope at which normal depth = critical depth):

9.2 Most Important GVF Profiles for GATE CE

10. Worked Examples (GATE CE Level)

Example 1 — Discharge in a Trapezoidal Canal (GATE CE 2021 type)

Problem: A trapezoidal irrigation canal has a base width of 4 m, side slopes 1.5:1 (H:V), depth of flow 1.2 m, bed slope 1/2500, and Manning’s n = 0.025. Find the discharge.

Answer: Q ≈ 4.94 m³/s

Example 2 — Most Economical Rectangular Section (GATE CE type)

Problem: Design the most economical rectangular channel to carry 6 m³/s with n = 0.013 and S = 1/1000. Find the width b and depth y.

Answer: y = 1.40 m; b = 2.80 m

Example 3 — Critical Depth in a Rectangular Channel (GATE CE 2020 type)

Problem: A rectangular channel 3 m wide carries a discharge of 9 m³/s. Find (a) the critical depth, (b) the critical velocity, and (c) the minimum specific energy.

Answer: y_c = 0.972 m; V_c = 3.09 m/s; E_min = 1.46 m

Example 4 — Froude Number and Flow Classification

Problem: Water flows at 2.4 m/s in a rectangular channel with a depth of 0.5 m. (a) Find the Froude number. (b) Classify the flow. (c) Find the critical depth for the same discharge.

Answer: Fr = 1.083 (Supercritical); y_c = 0.528 m

Example 5 — Manning’s n from Chezy’s C (GATE concept)

Problem: A concrete-lined trapezoidal canal has a hydraulic radius R = 1.5 m. Chezy’s C = 68 m^1/2/s. Find Manning’s n and the velocity if the bed slope is 1/1600.

Answer: n = 0.016; V = 2.08 m/s

Example 6 — GVF Profile Identification (GATE MCQ type)

Problem: A mild-slope channel has normal depth y_n = 1.8 m and critical depth y_c = 1.0 m. A dam is built downstream, raising the water surface to y = 2.5 m at the dam face. Identify the GVF profile and describe how depth varies moving upstream from the dam.

Answer: M1 (backwater) profile — depth decreases from 2.5 m at dam toward normal depth 1.8 m far upstream.

11. Common Mistakes

Mistake 1 — Confusing Hydraulic Radius R with Pipe Radius r

Error: Using R = D/2 (pipe radius) in Manning’s equation instead of R = A/P (hydraulic radius).
Root Cause: Both are called “R” in different contexts. Hydraulic radius R = A/P is completely different from the geometric radius of a circular cross-section (which is D/2). For a full circular pipe, R_hydraulic = (πD²/4)/(πD) = D/4 — not D/2.
Fix: Always compute R = A/P. For a full pipe: R = D/4. For a rectangular channel of width b and depth y: R = by/(b+2y). Never use D/2 as hydraulic radius.

Mistake 2 — Using Manning’s n = 0.013 for an Earthen Canal

Error: Applying n = 0.013 (correct for smooth concrete) to an unlined earthen irrigation canal.
Root Cause: Memorising n = 0.013 as a single value without noting it applies only to smooth concrete pipes and lined channels.
Fix: Earthen canal: n = 0.020–0.025. Concrete lined: n = 0.013–0.015. Natural river: n = 0.025–0.050. The problem statement always specifies the channel type — read it carefully before selecting n.

Mistake 3 — Applying y_c = (q²/g)^1/3 to Non-Rectangular Channels

Error: Using the rectangular critical depth formula y_c = (q²/g)^1/3 for trapezoidal or circular channels.
Root Cause: The formula y_c = (q²/g)^1/3 is derived specifically for a rectangular channel where T = b = constant. For a trapezoidal channel, T = b + 2zy changes with depth, so the critical condition A³/T = Q²/g must be solved implicitly.
Fix: For rectangular: use y_c = (q²/g)^1/3 directly. For any other shape: use A³/T = Q²/g and solve iteratively (trial and error with a table of A and T values at different y).

Mistake 4 — Misidentifying Normal Depth vs Critical Depth in GVF Problems

Error: Confusing normal depth y_n (depth for uniform flow at given Q, n, S₀) with critical depth y_c (depth at which Fr = 1 for given Q).
Root Cause: Both are special depths — y_n is found from Manning’s equation; y_c is found from the Froude number condition. They are equal only on a critical slope channel.
Fix: y_n → solve Manning’s Q = (1/n)AR^2/3S₀^1/2 for y. y_c → solve A³/T = Q²/g for y. Compare: y_n > y_c means mild slope (M-type); y_n < y_c means steep slope (S-type).

Mistake 5 — Wrong Sign Convention in GVF Profile Identification

Error: Saying depth increases upstream for an M1 profile (it decreases upstream — the profile asymptotes to normal depth upstream).
Root Cause: The dy/dx in the GVF equation is with respect to the downstream direction (x increasing downstream). dy/dx > 0 means depth increases downstream. An M1 profile has dy/dx > 0 → depth increases in the downstream direction (toward the dam) and decreases moving upstream.
Fix: Always state the direction clearly. “Depth increases downstream” and “depth decreases upstream” describe the same M1 profile. Draw a sketch with the dam on the right, flow from left to right, and show depth profile approaching y_n asymptotically on the left.

12. Frequently Asked Questions

Q1. Why does maximum discharge in a circular pipe occur at about 94% full, not completely full?

When a circular pipe flows full, the wetted perimeter P equals the full circumference πD. As the water level drops slightly below full (say to 94% of D), the water surface becomes a chord cutting across the top of the circle. At this point, the wetted perimeter decreases faster than the flow area A decreases — the top arc of the pipe is removed from P while a relatively small area is lost. Since R = A/P and Manning’s Q = (1/n)AR^2/3S^1/2, the hydraulic radius R actually increases as depth drops from 100% to about 94% full. This larger R more than compensates for the smaller A, giving a net increase in Q. At depths below 81% full, both A and R decrease rapidly, and Q falls below the full-pipe value. This counterintuitive result is critical for sewer design — sewers are typically designed to run at 80–90% full under peak flow to maintain a margin against surcharge while benefiting from the near-maximum discharge capacity.

Q2. What is the physical basis of Manning’s roughness coefficient n, and why does it not have SI units?

Manning’s n is an empirically determined coefficient that captures the combined effect of channel boundary roughness, irregular cross-section, channel alignment (bends and meanders), obstructions (vegetation, boulders), and scale effects on flow resistance. Dimensionally, n has units of s/m^1/3 (seconds per metre to the one-third power), but engineers treat it as dimensionless by convention — the unit is absorbed into the formula. When the Manning equation is used in imperial units (feet and seconds), a conversion factor of 1.486 appears: V = (1.486/n) R^2/3 S^1/2. The n values in the standard tables (0.013 for concrete, 0.025 for earthen canals) are in SI (metric) units. Using these n values with the metric formula gives correct results without any conversion factor. The “dimensionless” treatment works because engineers always use the same unit system consistently within a calculation.

Q3. How are GVF profiles used in the design of irrigation canals and flood control channels?

GVF profiles are essential for several practical design tasks. In irrigation canals, M1 backwater profiles determine how far upstream a regulating structure (cross regulator or head regulator) raises the water level — this is critical for ensuring adequate head for offtake structures serving distributary canals. If the backwater extends too far, it can cause overtopping of canal banks upstream. M2 drawdown profiles (just upstream of a fall or free overfall) determine where a canal bed lining must be extended to protect against scour in the transition zone. In flood control, M1 profiles behind embankments or weirs define the flood inundation extent and the required embankment height upstream. S2 profiles (just downstream of a slope break from mild to steep) and M3 profiles (below a sluice gate on a mild slope) both require a hydraulic jump to return flow to subcritical conditions — the location and energy loss in the jump are computed from the GVF profile. Standard Backwater Computation Methods (Step Method, Direct Integration) are specified in IS 1205 and are used for flood plain mapping in India.

Q4. When should Chezy’s formula be used instead of Manning’s equation in practice?

Manning’s equation has largely replaced Chezy’s formula in modern civil engineering practice because Manning’s n is more stable (nearly independent of hydraulic radius for a given surface material) and more extensively tabulated. Chezy’s C, by contrast, varies with hydraulic radius (C = R^1/6/n) — a larger channel has a higher C even with the same lining material. However, Chezy’s formula is still encountered in: (1) older IS codes and Irrigation Department records from pre-independence era India that specified design discharges using Kutter’s formula (which gives C); (2) situations where a Chezy C value is directly measured or specified in a contract document; (3) academic derivations where the C form is more compact; and (4) some coastal and tidal hydraulics applications. When given a Chezy C in a GATE problem, the simplest approach is to convert to Manning’s n = R^1/6/C and proceed with the more familiar Manning’s equation, or directly use V = C√(RS).