Open-Channel Hydraulics: Specific Energy, the Hydraulic Jump, Uniform and Gradually Varied Flow
E = y + V²/2g measured from the channel bed, is a function of depth with a minimum, and the depth at which it is minimum is the critical depth y_c = (q²/g)^(1/3) for a rectangular channel. Above y_c the flow is deep and slow, subcritical, with Fr < 1; below it the flow is shallow and fast, supercritical, with Fr > 1; and for every specific energy above the minimum there are two depths that carry the same discharge — the alternate depths. Which one a channel is actually at is decided by the slope and by what lies downstream, and that is the second organising idea: a subcritical flow is controlled from downstream and a supercritical flow from upstream, because a disturbance travels at √(gy) and can only move upstream if the flow is slower than that. When a supercritical flow must become subcritical it cannot do so smoothly, because the two depths do not join along the energy curve; it does so through a hydraulic jump, which is violent, short, and destroys energy — and must therefore be analysed by momentum rather than by energy, since the loss is exactly the unknown. Uniform flow is the equilibrium depth a long prismatic channel settles to under Manning’s V = (1/n) R^(2/3) S^(1/2), and comparing that normal depth with the critical depth classifies the slope as mild or steep. Everything between two such states is gradually varied flow, whose profile is fixed by the sign of dy/dx = (S₀ − S_f)/(1 − Fr²) — two signs, and the whole classification of M, S and C curves falls out of them.The Hydraulic Jump: Solve It With Momentum, Not Energy
A supercritical flow that has to become subcritical cannot slide along the specific energy curve from one alternate depth to the other, because doing so would require gaining energy on the way. What happens instead is a hydraulic jump: a short, turbulent, air-entraining transition in which the depth rises abruptly and a large part of the energy is destroyed. The analysis has to be done by momentum, not energy — and the reason is exact rather than stylistic. The energy loss is precisely the quantity you do not know, so an energy equation contains two unknowns and cannot be solved; a momentum equation across the jump contains only the pressure forces and the momentum flux, both computable from the depths, so it closes. Equating the specific force on either side of a rectangular channel gives the sequent depth relation y₂/y₁ = ½(√(1 + 8Fr₁²) − 1).
q = 4 m³/s per m arriving at a depth of y₁ = 0.5 m. The velocity is 8 m/s and Fr₁ = 8/√(9.81 × 0.5) = 3.612, firmly supercritical. The sequent depth is y₂ = 0.5 × 0.5 × (√(1 + 8 × 13.048) − 1) = 2.3164 m, a rise of a factor of 4.63, and downstream Fr₂ = 0.362, subcritical as required. Now count the energy. E₁ = 0.5 + 8²/19.62 = 3.762 m and E₂ = 2.3164 + 1.7268²/19.62 = 2.468 m, so the loss is 1.294 m — 34.4 per cent of the energy that arrived. Check it against the closed form ΔE = (y₂ − y₁)³/(4y₁y₂) = 1.8164³/(4 × 0.5 × 2.3164) = 1.294 m ✓, and the agreement to six figures is the sign that the momentum solution and the energy accounting are consistent. Three things follow. The loss grows very fast with Fr₁, as the cube of the depth difference, so a jump at Fr₁ = 2 dissipates only a few per cent while one at Fr₁ = 9 dissipates over 70 — which is why a stilling basin below a spillway is designed to force the jump to occur inside it, where the energy can be destroyed harmlessly on concrete rather than on the riverbed. The jump is a one-way device: it converts supercritical to subcritical and never the reverse, because the reverse would require creating energy. And the position matters as much as the existence — a jump that forms too far downstream has already scoured what it was meant to protect, which is why the tailwater depth is checked against y₂ rather than assumed adequate.Uniform Flow, and What Makes a Slope Mild
In a long prismatic channel of constant slope the depth settles to the value at which the component of gravity down the slope exactly balances the boundary resistance. That is the normal depth y_n, and Manning’s equation gives it: V = (1/n) R^(2/3) S^(1/2), with R = A/P the hydraulic radius. Work an example. A rectangular channel 5 m wide flowing 1.5 m deep at n = 0.015 on a slope of 0.001 has A = 7.5 m², P = 8 m and R = 0.9375 m, so V = 2.019 m/s and Q = 15.15 m³/s. Note that n is not dimensionless and carries units that make the formula work only in SI, and that the exponent on R means a channel’s efficiency depends on its shape, not only its area — which is why the most hydraulically efficient rectangular section is one with b = 2y, maximising R for a given area.
The useful step is the comparison that follows. For this channel q = 15.15/5 = 3.029 m³/s per m, so y_c = (3.029²/9.81)^(1/3) = 0.978 m. Since y_n = 1.5 m > y_c = 0.978 m, the normal flow is subcritical — Fr = 0.526 confirms it — and the slope is classified mild. Check it a second, independent way by computing the slope that would make this same channel critical at this same discharge: S_c = 0.003453, or 1 in 290, against the actual 1 in 1000. The actual slope is milder, and the two tests agree. That classification is not vocabulary; it is what determines which gradually varied profiles can exist on the reach, and therefore which way the water surface curves as it approaches any control.
Gradually Varied Flow: Two Signs Give Every Profile
Between a control and the reach’s normal depth the flow is neither uniform nor abrupt: the depth changes gradually over a long distance, and the governing equation is dy/dx = (S₀ − S_f)/(1 − Fr²). Everything about the shape of the water surface is in the signs of that fraction, and reading them is faster than remembering a table of profile names. The numerator compares the bed slope with the friction slope, so it is positive when the flow is deeper than normal — a deeper flow is slower, so it needs less slope to overcome friction. The denominator is positive when the flow is subcritical and negative when it is supercritical. Combine them: two positives or two negatives give a rising surface, one of each gives a falling one.
| Profile | Depth range | S₀ − S_f | 1 − Fr² | dy/dx | Where you see it |
|---|---|---|---|---|---|
| M1 | y > y_n > y_c | + | + | + backwater | Upstream of a dam or weir |
| M2 | y_n > y > y_c | − | + | − drawdown | Approaching a free overfall |
| M3 | y_n > y_c > y | − | − | + rising | Below a sluice gate, up to a jump |
Two boundary behaviours complete the picture and are worth reasoning out rather than memorising. As the depth approaches the normal depth, S_f → S₀, the numerator goes to zero and dy/dx → 0 — the profile becomes asymptotic to the normal depth line, approaching it but never crossing. As the depth approaches the critical depth, Fr → 1, the denominator goes to zero and dy/dx → ∞ — the water surface becomes vertical, which is the equation announcing that the gradually varied assumption has broken down and something abrupt is about to happen. That is exactly where a hydraulic jump or a free overfall appears, and it is why the theory predicts its own limit. The M3 profile in the table shows both ideas working together: water issuing fast and shallow from under a sluice gate rises along an M3 curve until it can no longer stay supercritical, and the transition to the subcritical M-regime happens through a jump rather than through the curve.
Unsteady flow and mobile channels
Everything above assumed the flow was steady and the bed was rigid. The GATE 2027 revision names both assumptions as things to relax. Unsteady flow — a gate closing, a sluice opening, a flood passing — propagates as a wave, and the quantity that governs it is the celerity of a small disturbance in still water, c = √(gy). At a depth of 2 m that is 4.429 m/s; at 1 m, 3.132; at 4 m, 6.264 — so a deeper channel carries information faster, which is why a flood wave outruns the water in it.
V = 1.5 m/s at a depth of y = 2 m. A surge forms and travels upstream at V + c relative to the bank — or rather, the disturbance moves at c = 4.429 m/s relative to the water, and the water is moving downstream at 1.5, so the wave moving against the flow travels upstream at c − V = 2.929 m/s and the one moving with it at c + V = 5.929 m/s. Both signs are real, and the sign of c − V is the whole content of the Froude number: here Fr = V/c = 1.5/4.429 = 0.339 < 1, the flow is subcritical, and a disturbance can travel both upstream and downstream. Make the flow supercritical, Fr > 1, and V > c: nothing can travel upstream at all, so a gate closing sends no warning to the water approaching it, which is the physical reason a supercritical reach is controlled from upstream and a subcritical one from downstream — the same rule that decided the direction of the gradually-varied-flow computation earlier in this chapter, arrived at from the wave side rather than the energy side. Mobile channels relax the other assumption. Where the bed is erodible the channel adjusts its own geometry until it is in regime with the sediment load it carries — which is the empirical territory of Lacey and Kennedy in the next chapter — and the governing threshold is the critical shear stress at which grains begin to move. The practical consequence is one every canal engineer meets: a channel designed to a fixed section in erodible material does not stay in that section, so it must be designed either to be non-scouring and non-silting at once, or to be lined and take the question away.Key takeaways
E = y + q²/2gy²has a minimum aty_c = (q²/g)^(1/3), whereFr = 1exactly andE_min = 1.5 y_c. Above the minimum everyEis shared by two alternate depths carrying the same discharge.- A disturbance travels at
√(gy), so subcritical flow is controlled from downstream and supercritical from upstream. Get that wrong and every number can be right while the water surface is drawn backwards. - Solve a jump with momentum, because the energy loss is the unknown.
y₂/y₁ = ½(√(1+8Fr₁²) − 1); at Fr₁ = 3.61 the depth rises 4.63× and 34% of the incoming energy is destroyed. - Manning gives the normal depth; comparing it with
y_cclassifies the slope. The worked channel hasy_n = 1.5 > y_c = 0.978, so mild — confirmed independently byS_c = 1 in 290against an actual 1 in 1000. - Every GVF profile falls out of two signs in
dy/dx = (S₀ − S_f)/(1 − Fr²).dy/dx → 0at normal depth (asymptotic) and→ ∞at critical — where the theory announces its own breakdown and a jump or an overfall appears.
Practice questions (8)
Attempt each one before opening the answer. Every explanation names the tempting wrong option as well as the right one, because that is where marks are lost.
A rectangular channel carries a discharge of 4 m³/s per metre width. The critical depth, in m, is ______ (round to three decimals)
Numerical answer — type the value.
Show answer
Answer: 1.177
1.177 m. For a rectangular channely_c = (q²/g)^(1/3) = (16/9.81)^(1/3) = (1.6310)^(1/3) = 1.17711 m. Three checks confirm it, and each is worth running once because they are the identities the topic is built on. First, the Froude number at that depth:V_c = q/y_c = 3.398 m/sandFr = V/√(gy) = 3.398/√(9.81 × 1.177) = 1.000000— critical depth is defined byFr = 1, and the cube-root formula is just that condition rearranged. Second, the velocity head at critical isV_c²/2g = 0.5886 m, which is exactlyy_c/2— a relation specific to rectangular channels and a fast way to spot an arithmetic slip. Third, the minimum specific energy isE_min = y_c + y_c/2 = 1.5 y_c = 1.766 m, so no rectangular channel can pass 4 m³/s per metre on less than 1.766 m of specific energy however its depth is arranged. What the critical depth is not is the depth the channel actually flows at. That is the normal depth, fixed by the slope and roughness through Manning, and the comparison between the two is what classifies the slope:y_n > y_cis mild and subcritical,y_n < y_cis steep and supercritical. Note also thaty_cdepends on the discharge alone for a given channel shape — not on the slope, not on the roughness, not onn. Change the bed slope andy_nmoves whiley_cstays exactly where it was.A hydraulic jump must be analysed using the momentum equation rather than the energy equation because
Show answer
Answer: B — the energy loss across the jump is itself the unknown, so the energy equation cannot be closed
The energy loss is itself the unknown. Write the energy equation across the jump and it containsy₁,y₂andh_L, three quantities of which two are unknown — it cannot be solved. Write the momentum equation and it contains only the pressure forces on the two faces and the momentum flux through them, all computable from the depths and the discharge, so it closes and yields the sequent depthy₂/y₁ = ½(√(1 + 8Fr₁²) − 1). Oncey₂is known, the energy equation becomes useful again — not to find the depths but to compute the loss as a result: fory₁ = 0.5 mandq = 4,E₁ = 3.762 mandE₂ = 2.468 m, soh_L = 1.294 m, which is 34.4 per cent of the incoming energy and matches the closed form(y₂ − y₁)³/(4y₁y₂)to six figures. Option A is loosely worded and wrong as stated: momentum is not conserved in general either — it is conserved here only because the external forces over the short length of the jump are small enough to neglect, and on a sloping or a stepped floor that assumption fails and the analysis changes. Option C is false, since specific energy is perfectly well defined on both sides. Option D misdescribes the flow, which is steady in the mean; the turbulence is intense but the jump does not move once it has positioned itself. The general lesson is worth carrying beyond this topic: use the conservation law whose unknown terms you can evaluate, and keep the other one to compute what you were not able to assume.Water at 4 m³/s per metre width enters a hydraulic jump at a depth of 0.5 m. The sequent (conjugate) depth, in m, is ______ (round to three decimals)
Numerical answer — type the value.
Show answer
Answer: 2.316
2.316 m. First the incoming Froude number:V₁ = q/y₁ = 4/0.5 = 8 m/sandFr₁ = 8/√(9.81 × 0.5) = 8/2.2147 = 3.612, comfortably supercritical. Then the sequent depth relation, which comes from equating specific force across the jump:y₂ = ½ y₁ (√(1 + 8Fr₁²) − 1) = 0.25 × (√(1 + 104.38) − 1) = 0.25 × 9.266 = 2.3164 m. The depth ratio is4.63, and downstreamFr₂ = 0.362, subcritical as it must be. The energy accounting is the interesting part and confirms the answer independently:E₁ = 0.5 + 8²/19.62 = 3.762 m,E₂ = 2.3164 + 1.727²/19.62 = 2.468 m, so 1.294 m — 34.4 per cent of the incoming energy — has been destroyed, and the closed form(y₂ − y₁)³/(4y₁y₂)returns the same 1.294 m. That agreement between a momentum-derived depth and an energy-derived loss is the check to run. Two cautions. Sequent depths are not alternate depths: alternate depths share the same specific energy and are what a smooth transition would connect, while sequent depths share the same specific force and are what a jump connects — the jump exists precisely because the flow cannot move between alternate depths without gaining energy. And the loss grows as the cube of(y₂ − y₁), so it climbs steeply withFr₁: a jump atFr₁ = 2dissipates only a few per cent while one atFr₁ = 9dissipates over 70, which is why a stilling basin below a spillway is sized to hold the jump where the concrete can take it.A rectangular channel 5 m wide flows 1.5 m deep with n = 0.015 on a bed slope of 0.001. The discharge, in m³/s, is ______ (round to two decimals)
Numerical answer — type the value.
Show answer
Answer: 15.15
15.15 m³/s.A = 5 × 1.5 = 7.5 m², the wetted perimeter isP = 5 + 2(1.5) = 8 m— the free surface is not wetted and is never included — soR = A/P = 0.9375 m. Manning givesV = (1/n)R^(2/3)S^(1/2) = (1/0.015)(0.9375)^(2/3)(0.001)^(1/2) = 66.667 × 0.9579 × 0.03162 = 2.019 m/s, andQ = AV = 15.15 m³/s. The wetted perimeter is where most marks are lost: including the water surface givesP = 11,R = 0.68and an answer 12 per cent low. The step worth taking next is the classification. Hereq = 15.15/5 = 3.029 m³/s per m, soy_c = (3.029²/9.81)^(1/3) = 0.978 m. Sincey_n = 1.5 m > y_c = 0.978 m, the flow is subcritical —Fr = 2.019/√(9.81 × 1.5) = 0.526agrees — and the slope is mild. Confirm it a second way by asking what slope would make this same channel critical at this same discharge:S_c = 0.003453, or 1 in 290, against the actual 1 in 1000, so the actual bed is milder than critical and the two tests agree. That classification is not terminology; it fixes which gradually varied profiles can exist on the reach and therefore which way the surface curves near any control. Two properties of Manning worth carrying:nis not dimensionless, so the formula as written holds only in SI units, and theR^(2/3)means efficiency depends on shape — the most efficient rectangular section hasb = 2y, maximisingRfor a given area.In a subcritical flow in an open channel, a disturbance such as a weir placed downstream will
Show answer
Answer: B — affect the depth upstream, because the flow is slower than the wave celerity
It affects the depth upstream, because the flow is slower than the wave celerity. A small surface disturbance travels relative to the water at the celerityc = √(gy). In a subcritical flowV < c, that isFr < 1, so the disturbance can propagate upstream against the current, and the influence of a downstream weir, gate or reservoir reaches back up the channel — sometimes for kilometres, as the M1 backwater curve above a dam. That is the whole meaning of the statement that subcritical flow is controlled from downstream. In supercritical flowV > c, no signal can travel upstream, and the flow is controlled from upstream instead — which is why the water leaving a sluice gate is oblivious to what lies ahead of it until it is forced to notice through a jump. Option A appeals to an intuition that is simply wrong: what travels upstream is not water but information, a pressure wave on the free surface, and it can travel against a current slower than itself. Option C reverses the situation, since a steep slope means supercritical flow and therefore no upstream influence. Option D describes what happens in the opposite case — a jump forms where a supercritical flow is forced to become subcritical, not where a subcritical flow meets an obstruction. The practical value of getting this right is that it tells you which end to start the calculation from: a gradually varied flow profile is integrated upstream from the control in subcritical flow and downstream from it in supercritical, and starting at the wrong end produces a profile that diverges instead of settling to the normal depth.On a mild slope, the water surface profile upstream of a dam is classified as
Show answer
Answer: A — M1, a backwater curve with the depth rising in the downstream direction
M1, a backwater curve. Work it from the signs rather than from a remembered picture. A dam raises the depth above the normal depth, so the reach upstream of it hasy > y_n > y_c— that is the M1 zone. Now readdy/dx = (S₀ − S_f)/(1 − Fr²). The flow is deeper than normal, so it is slower than normal, so the friction slopeS_fis less than the bed slopeS₀and the numerator is positive. The flow is subcritical, soFr < 1and the denominator is positive too. Hencedy/dx > 0: the depth increases in the downstream direction, which is exactly what a backwater curve looks like — the surface rises as it approaches the dam. Option B is the M2 zone,y_n > y > y_c, which occurs where a mild channel approaches a free overfall and the depth is drawn down towards critical. Option C is the M3 zone, below critical, which is what happens under a sluice gate on a mild slope before a jump. Option D confuses the structure with the bed: a dam does not change the slope of the channel it sits in, and the S-profiles belong to channels whose bed slope is steep enough that the normal depth is supercritical. Two limiting behaviours make the M1 curve easy to sketch. Far upstreamy → y_n, soS_f → S₀, the numerator vanishes and the profile becomes asymptotic to the normal depth line. At the dam the surface approaches the horizontal, since the velocity there is nearly zero. Between the two it is a smooth rising curve, and its length — often several kilometres — is why reservoir backwater is computed rather than guessed.For a rectangular channel, the minimum specific energy expressed as a multiple of the critical depth is ______
Numerical answer — type the value.
Show answer
Answer: 1.5
1.5. At critical depthFr = 1, soV_c = √(g y_c)and the velocity head isV_c²/2g = g y_c/2g = y_c/2. ThenE_min = y_c + y_c/2 = 1.5 y_c. The relation is specific to a rectangular channel — for a general section the critical condition isQ²T/(gA³) = 1withTthe top width, and the multiple is not 1.5. Check it on the worked example:q = 4 m³/s per mgivesy_c = 1.177 mandE_min = 1.766 m, and1.766/1.177 = 1.500✓. Two things this result is good for. It gives an instant lower bound: no rectangular channel can pass a givenqon less specific energy than1.5 y_c, whatever depth it flows at, so if a transition offers less than that, the flow must choke and back up until it has enough. That is exactly what happens at a broad-crested weir or a raised bed — the flow passes through critical at the control and the upstream depth rises to supply1.5 y_cplus the step. And it makes the specific energy curve easy to sketch: the minimum sits at(1.5 y_c, y_c), the upper limb approaches the 45° line asymptotically as the velocity head vanishes, and the lower limb runs to infinity asy → 0because the velocity head blows up as1/y². The one-third/two-thirds split is worth noting too — at critical depth the flow carries two-thirds of its specific energy as depth and one-third as velocity head, which is the most efficient partition available and the reason the condition is called critical.Which of the following are true of critical flow in an open channel? (One or more)
Show answer
Answer: A — The specific energy is a minimum for the given discharge; B — The Froude number equals 1; C — The critical depth depends on the discharge but not on the bed slope or roughness
A, B and C. These are three descriptions of the same condition. A is the definition by whichy_cis usually derived: settingdE/dy = 0onE = y + q²/2gy²givesq² = gy³, hencey_c = (q²/g)^(1/3). B is the same condition rearranged —q² = gy³isV²y² = gy³, that isV² = gyandFr = 1— and the two are equivalent statements, not separate facts. C is the one that catches people, and it matters. The critical depth follows from the discharge and the channel shape alone: nonand noS₀appear anywhere iny_c = (q²/g)^(1/3). So changing the bed slope moves the normal depth while leavingy_cexactly where it was, and that is precisely why comparing the two classifies the slope as mild or steep. It also meansy_cis the right reference against which to judge any change in a channel, since it stays put while everything else moves. D is false, and it is the trap. At critical flowV = c = √(gy), so a disturbance travelling upstream at the celerity is exactly held stationary by the current — it neither advances nor is swept away. Upstream propagation requiresV < c, that is subcritical flow. This is why critical flow is a control: a section forced to be critical, at a weir crest or a raised bed, imposes a unique relation between depth and discharge and separates the upstream reach from the downstream one, which is what makes such a section a reliable place to measure flow.