Browse papers
A

Section A: Long Answer Questions

Attempt all questions.

5 questions
1long10 marks

A trapezoidal channel carries water in uniform flow. The channel has a bed width b=6.0 mb = 6.0\ \text{m}, side slopes 1:21:2 (vertical : horizontal), a longitudinal bed slope S0=0.0009S_0 = 0.0009 and Manning's roughness n=0.018n = 0.018. The normal depth of flow is y=1.5 my = 1.5\ \text{m}.

(a) Compute the discharge using Manning's equation. (b) Determine the average boundary shear stress and the state (sub/supercritical) of the flow. (c) Briefly explain the physical meaning of the conveyance of a channel section.

Given: b=6.0 mb = 6.0\ \text{m}, side slope z=2z = 2 (H:V), S0=0.0009S_0 = 0.0009, n=0.018n = 0.018, y=1.5 my = 1.5\ \text{m}, g=9.81 m/s2g = 9.81\ \text{m/s}^2.

(a) Discharge (Manning's equation)

Flow area:

A=(b+zy)y=(6.0+2×1.5)×1.5=(6.0+3.0)×1.5=13.5 m2A = (b + zy)y = (6.0 + 2\times1.5)\times1.5 = (6.0+3.0)\times1.5 = 13.5\ \text{m}^2

Wetted perimeter:

P=b+2y1+z2=6.0+2×1.51+4=6.0+3.0×2.2360680=6.0+6.7082039=12.7082 mP = b + 2y\sqrt{1+z^2} = 6.0 + 2\times1.5\sqrt{1+4} = 6.0 + 3.0\times2.2360680 = 6.0 + 6.7082039 = 12.7082\ \text{m}

Hydraulic radius:

R=AP=13.512.7082=1.06231 mR = \frac{A}{P} = \frac{13.5}{12.7082} = 1.06231\ \text{m}

Manning's equation:

Q=1nAR2/3S01/2Q = \frac{1}{n}A R^{2/3} S_0^{1/2}

R2/3=1.062312/3R^{2/3} = 1.06231^{2/3}. Since ln1.06231=0.060456\ln 1.06231 = 0.060456, 23×0.060456=0.040304\tfrac{2}{3}\times0.060456 = 0.040304, e0.040304=1.041128e^{0.040304} = 1.041128.

S01/2=0.0009=0.03S_0^{1/2} = \sqrt{0.0009} = 0.03.

Q=10.018×13.5×1.041128×0.03=55.5556×13.5×1.041128×0.03Q = \frac{1}{0.018}\times13.5\times1.041128\times0.03 = 55.5556\times13.5\times1.041128\times0.03 Q=55.5556×13.5=750.000;750.000×1.041128=780.846;×0.03Q = 55.5556\times13.5 = 750.000;\quad 750.000\times1.041128 = 780.846;\quad \times0.03 Q=23.43 m3/s\boxed{Q = 23.43\ \text{m}^3/\text{s}}

(b) Boundary shear stress and flow state

Average boundary shear stress:

τ0=γRS0=9810×1.06231×0.0009=9.379 N/m2\tau_0 = \gamma R S_0 = 9810\times1.06231\times0.0009 = 9.379\ \text{N/m}^2 τ09.38 Pa\boxed{\tau_0 \approx 9.38\ \text{Pa}}

Flow state — mean velocity:

V=QA=23.4313.5=1.7356 m/sV = \frac{Q}{A} = \frac{23.43}{13.5} = 1.7356\ \text{m/s}

Top width T=b+2zy=6.0+2×2×1.5=12.0 mT = b + 2zy = 6.0 + 2\times2\times1.5 = 12.0\ \text{m}; hydraulic (mean) depth D=A/T=13.5/12.0=1.125 mD = A/T = 13.5/12.0 = 1.125\ \text{m}.

Froude number:

Fr=VgD=1.73569.81×1.125=1.735611.0363=1.73563.32209=0.5224Fr = \frac{V}{\sqrt{gD}} = \frac{1.7356}{\sqrt{9.81\times1.125}} = \frac{1.7356}{\sqrt{11.0363}} = \frac{1.7356}{3.32209} = 0.5224

Since Fr=0.52<1Fr = 0.52 < 1, the flow is subcritical.

(c) Conveyance

Writing Q=KS01/2Q = K S_0^{1/2}, the conveyance is K=1nAR2/3K = \dfrac{1}{n}A R^{2/3}. It is a purely geometric–roughness measure of a section's ability to carry flow: it groups all terms that depend only on the channel shape, size and roughness, separating them from the energy slope. For a given KK, discharge varies with S0\sqrt{S_0}. Here K=Q/S01/2=23.43/0.03=781.0 m3/sK = Q/S_0^{1/2} = 23.43/0.03 = 781.0\ \text{m}^3/\text{s}.

uniform-flowmanning-equationopen-channel
2long10 marks

Water flows in a rectangular channel 4.0 m4.0\ \text{m} wide at a discharge of Q=12 m3/sQ = 12\ \text{m}^3/\text{s}.

(a) Derive the condition for critical flow in a rectangular channel from the specific-energy principle, and obtain the relation Ec=32ycE_c = \tfrac{3}{2}y_c. (b) Compute the critical depth, critical velocity and minimum specific energy. (c) If the actual depth of flow is 1.8 m1.8\ \text{m}, compute the specific energy and the alternate depth carrying the same discharge with the same specific energy.

Given: b=4.0 mb = 4.0\ \text{m}, Q=12 m3/sQ = 12\ \text{m}^3/\text{s}, g=9.81 m/s2g = 9.81\ \text{m/s}^2. Unit discharge q=Q/b=12/4.0=3.0 m2/sq = Q/b = 12/4.0 = 3.0\ \text{m}^2/\text{s}.

(a) Critical-flow condition

Specific energy: E=y+V22g=y+q22gy2E = y + \dfrac{V^2}{2g} = y + \dfrac{q^2}{2gy^2} (per unit width, rectangular).

For minimum EE at fixed qq, set dEdy=0\dfrac{dE}{dy}=0:

dEdy=1q2gy3=0    yc=(q2g)1/3\frac{dE}{dy} = 1 - \frac{q^2}{g y^3} = 0 \;\Rightarrow\; y_c = \left(\frac{q^2}{g}\right)^{1/3}

At this depth q2gyc3=1\dfrac{q^2}{g y_c^3} = 1, i.e. Vc2gyc=1Fr=1\dfrac{V_c^2}{g y_c} = 1 \Rightarrow Fr=1 (critical). Then

Ec=yc+q22gyc2=yc+yc2=32yc.E_c = y_c + \frac{q^2}{2g y_c^2} = y_c + \frac{y_c}{2} = \tfrac{3}{2}y_c.

(b) Critical depth, velocity, minimum specific energy

yc=(q2g)1/3=(3.029.81)1/3=(9.09.81)1/3=(0.917431)1/3y_c = \left(\frac{q^2}{g}\right)^{1/3} = \left(\frac{3.0^2}{9.81}\right)^{1/3} = \left(\frac{9.0}{9.81}\right)^{1/3} = (0.917431)^{1/3}

ln0.917431=0.086178\ln 0.917431 = -0.086178; /3=0.028726/3 = -0.028726; e0.028726=0.971682e^{-0.028726}=0.971682.

yc=0.9717 m\boxed{y_c = 0.9717\ \text{m}} Vc=qyc=3.00.9717=3.0874 m/sVc3.09 m/sV_c = \frac{q}{y_c} = \frac{3.0}{0.9717} = 3.0874\ \text{m/s}\quad\Rightarrow\boxed{V_c \approx 3.09\ \text{m/s}} Ec=32yc=1.5×0.9717=1.4575 mE_c = \tfrac{3}{2}y_c = 1.5\times0.9717 = \boxed{1.4575\ \text{m}}

(c) Specific energy at y=1.8 my = 1.8\ \text{m} and alternate depth

V=qy=3.01.8=1.6667 m/s,V22g=1.666722×9.81=2.7777819.62=0.14157 mV = \frac{q}{y} = \frac{3.0}{1.8} = 1.6667\ \text{m/s},\qquad \frac{V^2}{2g} = \frac{1.6667^2}{2\times9.81} = \frac{2.77778}{19.62} = 0.14157\ \text{m} E=1.8+0.14157=1.9416 m(>Ec, so depth is valid; flow is subcritical).E = 1.8 + 0.14157 = \boxed{1.9416\ \text{m}}\quad(> E_c,\ \text{so depth is valid; flow is subcritical}).

Alternate depth yy' satisfies E=y+q22gy2E = y' + \dfrac{q^2}{2g y'^2}, i.e.

1.9416=y+9.019.62y2=y+0.458716y2.1.9416 = y' + \frac{9.0}{19.62\,y'^2} = y' + \frac{0.458716}{y'^2}.

Solve iteratively for the supercritical root (y<ycy' < y_c):

  • y=0.6y'=0.6: 0.6+0.458716/0.36=0.6+1.27421=1.874210.6 + 0.458716/0.36 = 0.6+1.27421=1.87421 (low)
  • y=0.58y'=0.58: 0.58+0.458716/0.3364=0.58+1.36361=1.943610.58 + 0.458716/0.3364 = 0.58+1.36361=1.94361 (≈)
  • y=0.581y'=0.581: 0.581+0.458716/0.337561=0.581+1.358918=1.9399180.581 + 0.458716/0.337561 = 0.581+1.358918=1.939918

Interpolating, y0.5804 my' \approx 0.5804\ \text{m}. Check: 0.5804+0.458716/0.336864=0.5804+1.361729=1.942130.5804 + 0.458716/0.336864 = 0.5804+1.361729 = 1.94213.

y0.580 m (supercritical alternate depth)\boxed{y' \approx 0.580\ \text{m}\ (\text{supercritical alternate depth})}
specific-energycritical-flowrectangular-channel
3long8 marks

A hydraulic jump forms in a horizontal rectangular channel 5.0 m5.0\ \text{m} wide carrying Q=30 m3/sQ = 30\ \text{m}^3/\text{s}. The depth before the jump is y1=0.6 my_1 = 0.6\ \text{m}.

(a) Compute the upstream Froude number and classify the jump. (b) Determine the sequent (conjugate) depth y2y_2. (c) Compute the energy loss in the jump and the efficiency of the jump (ratio E2/E1E_2/E_1).

Given: b=5.0 mb = 5.0\ \text{m}, Q=30 m3/sQ = 30\ \text{m}^3/\text{s}, y1=0.6 my_1 = 0.6\ \text{m}, g=9.81 m/s2g = 9.81\ \text{m/s}^2.

Unit discharge q=30/5.0=6.0 m2/sq = 30/5.0 = 6.0\ \text{m}^2/\text{s}.

(a) Upstream Froude number

V1=qy1=6.00.6=10.0 m/sV_1 = \frac{q}{y_1} = \frac{6.0}{0.6} = 10.0\ \text{m/s} Fr1=V1gy1=10.09.81×0.6=10.05.886=10.02.42611=4.1218Fr_1 = \frac{V_1}{\sqrt{g y_1}} = \frac{10.0}{\sqrt{9.81\times0.6}} = \frac{10.0}{\sqrt{5.886}} = \frac{10.0}{2.42611} = 4.1218 Fr1=4.12\boxed{Fr_1 = 4.12}

Since 4.5>Fr1>2.54.5 > Fr_1 > 2.5, this is a steady (oscillating-to-steady boundary) jump — strictly 2.5<Fr1<4.52.5<Fr_1<4.5 is the oscillating range; with Fr1=4.12Fr_1=4.12 it is a well-developed oscillating jump approaching the steady type.

(b) Sequent depth (Belanger equation)

y2y1=12(1+8Fr121)\frac{y_2}{y_1} = \frac{1}{2}\left(\sqrt{1+8Fr_1^2}-1\right) 8Fr12=8×4.12182=8×16.9893=135.914;1+135.914=136.9148Fr_1^2 = 8\times4.1218^2 = 8\times16.9893 = 135.914;\quad 1+135.914 = 136.914 136.914=11.7011\sqrt{136.914} = 11.7011 y2y1=12(11.70111)=10.70112=5.35057\frac{y_2}{y_1} = \frac{1}{2}(11.7011-1) = \frac{10.7011}{2} = 5.35057 y2=5.35057×0.6=3.210 my_2 = 5.35057\times0.6 = \boxed{3.210\ \text{m}}

(c) Energy loss and efficiency

Energy loss in a rectangular jump:

ΔE=(y2y1)34y1y2=(3.2100.6)34×0.6×3.210=(2.610)37.704=17.77967.704=2.3079 m\Delta E = \frac{(y_2-y_1)^3}{4 y_1 y_2} = \frac{(3.210-0.6)^3}{4\times0.6\times3.210} = \frac{(2.610)^3}{7.704} = \frac{17.7796}{7.704} = 2.3079\ \text{m} ΔE2.31 m\boxed{\Delta E \approx 2.31\ \text{m}}

Specific energies:

E1=y1+V122g=0.6+10.0219.62=0.6+5.0968=5.6968 mE_1 = y_1 + \frac{V_1^2}{2g} = 0.6 + \frac{10.0^2}{19.62} = 0.6 + 5.0968 = 5.6968\ \text{m} V2=qy2=6.03.210=1.8692 m/s;E2=3.210+1.8692219.62=3.210+3.4939119.62=3.210+0.17808=3.3881 mV_2 = \frac{q}{y_2} = \frac{6.0}{3.210} = 1.8692\ \text{m/s};\quad E_2 = 3.210 + \frac{1.8692^2}{19.62} = 3.210 + \frac{3.49391}{19.62} = 3.210 + 0.17808 = 3.3881\ \text{m}

Check: E1E2=5.69683.3881=2.3087 mΔEE_1 - E_2 = 5.6968 - 3.3881 = 2.3087\ \text{m} \approx \Delta E (consistent).

Efficiency:

E2E1=3.38815.6968=0.594759.5%\frac{E_2}{E_1} = \frac{3.3881}{5.6968} = 0.5947 \Rightarrow \boxed{59.5\%}

The jump dissipates about 40.5%40.5\% of the upstream specific energy.

hydraulic-jumpenergy-lossrectangular-channel
4long8 marks

(a) Derive the dynamic equation of gradually varied flow (GVF) in the form dydx=S0Sf1Fr2\dfrac{dy}{dx} = \dfrac{S_0 - S_f}{1 - Fr^2}, stating all assumptions. (b) A wide rectangular channel (q=4 m2/sq = 4\ \text{m}^2/\text{s}) has S0=0.004S_0 = 0.004 and n=0.020n = 0.020. Compute the normal depth yny_n and critical depth ycy_c, classify the channel slope, and sketch/name the possible GVF profiles. (c) State whether, for a depth y=1.2 my = 1.2\ \text{m} on this slope, the surface rises or falls in the flow direction (give reasoning).

(a) Derivation of the GVF equation

Total head above a datum: H=z+y+V22gH = z + y + \dfrac{V^2}{2g}.

Assumptions: (i) steady flow; (ii) gradual depth change so pressure is hydrostatic and streamlines nearly parallel; (iii) the head loss at a section equals that of uniform flow at the same depth/velocity, so SfS_f from Manning/Chezy applies; (iv) channel slope small (cosθ1\cos\theta\approx1); (v) prismatic channel, α=1\alpha=1; (vi) constant nn.

Differentiate w.r.t. xx (taking xx positive downstream):

dHdx=dzdx+dydx+ddx ⁣(V22g)\frac{dH}{dx} = \frac{dz}{dx} + \frac{dy}{dx} + \frac{d}{dx}\!\left(\frac{V^2}{2g}\right)

Now dHdx=Sf\dfrac{dH}{dx} = -S_f (energy slope), dzdx=S0\dfrac{dz}{dx} = -S_0 (bed slope). Also ddx ⁣(V22g)=Fr2dydx\dfrac{d}{dx}\!\left(\dfrac{V^2}{2g}\right) = -Fr^2\dfrac{dy}{dx} for a given QQ. Substituting:

Sf=S0+dydx(1Fr2)    dydx=S0Sf1Fr2-S_f = -S_0 + \frac{dy}{dx}(1 - Fr^2)\;\Rightarrow\; \boxed{\frac{dy}{dx} = \frac{S_0 - S_f}{1 - Fr^2}}

(b) Normal and critical depths; slope classification

Wide channel ⇒ RyR \approx y. Manning (per unit width): q=1ny5/3S01/2q = \dfrac{1}{n} y^{5/3} S_0^{1/2}.

yn=(qnS01/2)3/5=(4×0.0200.004)3/5=(0.080.0632456)3/5=(1.264911)3/5y_n = \left(\frac{q\,n}{S_0^{1/2}}\right)^{3/5} = \left(\frac{4\times0.020}{\sqrt{0.004}}\right)^{3/5} = \left(\frac{0.08}{0.0632456}\right)^{3/5} = (1.264911)^{3/5}

ln1.264911=0.234923\ln 1.264911 = 0.234923; ×0.6=0.140954\times0.6 = 0.140954; e0.140954=1.151367e^{0.140954} = 1.151367.

yn=1.151 m\boxed{y_n = 1.151\ \text{m}}

Critical depth (wide rectangular): yc=(q2/g)1/3=(16/9.81)1/3=(1.630989)1/3y_c = (q^2/g)^{1/3} = (16/9.81)^{1/3} = (1.630989)^{1/3}. ln1.630989=0.489149\ln 1.630989 = 0.489149; /3=0.163050/3 = 0.163050; e0.163050=1.177099e^{0.163050} = 1.177099.

yc=1.177 m\boxed{y_c = 1.177\ \text{m}}

Since yn(1.151)<yc(1.177)y_n (1.151) < y_c (1.177), the normal flow is supercritical ⇒ the slope is steep (S).

Possible profiles on a steep slope: S1 (for y>ycy > y_c), S2 (for yc>y>yny_c > y > y_n), S3 (for y<yny < y_n).

ASCII sketch (steep slope zones):

         ~~~~~~ S1  (y > y_c)
----------------------------- critical depth line  y_c=1.177
   \___ S2 (y_n<y<y_c)
----------------------------- normal depth line    y_n=1.151
   ___/  S3 (y < y_n)
============================= bed (steep)

(c) Behaviour at y=1.2 my = 1.2\ \text{m}

Here y=1.2 m>yc=1.177>yn=1.151y = 1.2\ \text{m} > y_c = 1.177 > y_n = 1.151, so this is the S1 zone. Numerator S0SfS_0 - S_f: since y>yny > y_n, the actual depth exceeds normal ⇒ Sf<S0S_f < S_0 ⇒ numerator >0>0. Denominator 1Fr21 - Fr^2: since y>ycy > y_c the flow is subcritical ⇒ Fr<1Fr<1 ⇒ denominator >0>0. Thus dydx>0\dfrac{dy}{dx} > 0: the water surface rises in the flow direction (S1 backwater profile, asymptotic to the horizontal far upstream-of-an-obstruction and tending away from ycy_c).

gradually-varied-flowwater-surface-profilesclassification
5long9 marks

An unlined trapezoidal channel of most-efficient (best hydraulic) section is to be designed to carry Q=20 m3/sQ = 20\ \text{m}^3/\text{s}. The side slopes are fixed at 1:1.51:1.5 (V:H), Manning's n=0.022n = 0.022, and the bed slope is S0=0.0006S_0 = 0.0006.

(a) Derive the conditions for the most efficient trapezoidal section and state the relation between bed width, depth and side slope. (b) Design the channel: find the depth yy, bed width bb, and check the mean velocity against a permissible non-scouring velocity of 1.0 m/s1.0\ \text{m/s}.

Given: Q=20 m3/sQ = 20\ \text{m}^3/\text{s}, side slope z=1.5z = 1.5 (H:V), n=0.022n = 0.022, S0=0.0006S_0 = 0.0006, g=9.81g=9.81.

(a) Most-efficient trapezoidal section

For fixed AA and zz, the section is most efficient when PP is minimum (maximises RR, hence QQ). A=(b+zy)yA = (b+zy)y, P=b+2y1+z2P = b + 2y\sqrt{1+z^2}. From AA: b=A/yzyb = A/y - zy. So

P=Ayzy+2y1+z2.P = \frac{A}{y} - zy + 2y\sqrt{1+z^2}.

Set dP/dy=0dP/dy = 0 at constant AA:

dPdy=Ay2z+21+z2=0A=y2(21+z2z).\frac{dP}{dy} = -\frac{A}{y^2} - z + 2\sqrt{1+z^2} = 0 \Rightarrow A = y^2\left(2\sqrt{1+z^2} - z\right).

Equating to A=(b+zy)yA=(b+zy)y gives the best-section relation:

b+zy=y(21+z2z)    b=y(21+z22z).\boxed{b + zy \,=\, y\left(2\sqrt{1+z^2}-z\right)\;\Rightarrow\; b = y\left(2\sqrt{1+z^2}-2z\right).}

A consequence is R=y/2R = y/2 (hydraulic radius equals half the depth) for the most efficient trapezoid.

(b) Design

With z=1.5z = 1.5: 1+z2=1+2.25=3.25=1.802776\sqrt{1+z^2} = \sqrt{1+2.25} = \sqrt{3.25} = 1.802776.

Bed width: b=y(2×1.8027762×1.5)=y(3.6055513.0)=0.605551yb = y(2\times1.802776 - 2\times1.5) = y(3.605551 - 3.0) = 0.605551\,y.

Area: A=y2(21+z2z)=y2(3.6055511.5)=2.105551y2A = y^2(2\sqrt{1+z^2} - z) = y^2(3.605551 - 1.5) = 2.105551\,y^2.

Hydraulic radius R=y/2=0.5yR = y/2 = 0.5y.

Manning: Q=1nAR2/3S01/2Q = \dfrac{1}{n} A R^{2/3} S_0^{1/2}.

20=10.022(2.105551y2)(0.5y)2/3(0.0006)1/2.20 = \frac{1}{0.022}\,(2.105551\,y^2)(0.5y)^{2/3}(0.0006)^{1/2}.

(0.5)2/3(0.5)^{2/3}: ln0.5=0.693147\ln0.5=-0.693147, ×23=0.462098\times\tfrac23=-0.462098, e0.462098=0.629961e^{-0.462098}=0.629961. So (0.5y)2/3=0.629961y2/3(0.5y)^{2/3}=0.629961\,y^{2/3}. S01/2=0.0006=0.0244949S_0^{1/2}=\sqrt{0.0006}=0.0244949.

20=10.022×2.105551×0.629961×0.0244949  y2+2/320 = \frac{1}{0.022}\times2.105551\times0.629961\times0.0244949\;y^{2+2/3}

Coefficient: 10.022=45.4545\frac{1}{0.022}=45.4545; 45.4545×2.105551=95.706945.4545\times2.105551 = 95.7069; ×0.629961=60.2898\times0.629961 = 60.2898; ×0.0244949=1.476788\times0.0244949 = 1.476788.

20=1.476788y8/3y8/3=201.476788=13.54295.20 = 1.476788\,y^{8/3} \Rightarrow y^{8/3} = \frac{20}{1.476788} = 13.54295.

y=13.542953/8y = 13.54295^{3/8}: ln13.54295=2.605910\ln13.54295 = 2.605910; ×0.375=0.977216\times0.375 = 0.977216; e0.977216=2.656832e^{0.977216} = 2.656832.

y=2.657 m\boxed{y = 2.657\ \text{m}}

Bed width:

b=0.605551×2.657=1.609 mb = 0.605551\times2.657 = \boxed{1.609\ \text{m}}

Velocity check:

A=2.105551×2.6572=2.105551×7.059649=14.86489 m2A = 2.105551\times2.657^2 = 2.105551\times7.059649 = 14.86489\ \text{m}^2 V=QA=2014.86489=1.3455 m/sV = \frac{Q}{A} = \frac{20}{14.86489} = 1.3455\ \text{m/s}

Since V=1.35 m/s>1.0 m/sV = 1.35\ \text{m/s} > 1.0\ \text{m/s} (permissible non-scouring velocity), the design velocity exceeds the permissible value. The section would scour; the bed slope must be flattened (reduce S0S_0) or the channel lined / widened from the strict best section. (If a higher permissible velocity were allowed, e.g. for a lined channel, the section is acceptable.)

channel-designbest-hydraulic-sectiontrapezoidal
B

Section B: Short Answer Questions

Attempt all questions.

6 questions
6short5 marks

Water flows at 2.5 m/s2.5\ \text{m/s} in a rigid steel pipe of length 1200 m1200\ \text{m}. The pressure-wave (celerity) speed is a=1200 m/sa = 1200\ \text{m/s}. A valve at the downstream end is closed.

(a) State the difference between sudden and gradual valve closure (in terms of the pipe period T=2L/aT = 2L/a). (b) Compute the maximum (Joukowsky) pressure rise for instantaneous closure. (c) If the valve is closed in 3.0 s3.0\ \text{s}, is the closure rapid or slow, and estimate the pressure rise (head) for this case.

Given: V0=2.5 m/sV_0 = 2.5\ \text{m/s}, L=1200 mL = 1200\ \text{m}, a=1200 m/sa = 1200\ \text{m/s}, ρ=1000 kg/m3\rho = 1000\ \text{kg/m}^3, g=9.81 m/s2g = 9.81\ \text{m/s}^2.

Pipe period: T=2La=2×12001200=2.0 sT = \dfrac{2L}{a} = \dfrac{2\times1200}{1200} = 2.0\ \text{s}.

(a) Sudden vs gradual closure

Closure is sudden (rapid) if the closure time tcT=2L/at_c \le T = 2L/a — the reflected relief wave has not yet returned to the valve, so the full Joukowsky head develops. It is gradual (slow) if tc>2L/at_c > 2L/a — the returning negative wave partially relieves the pressure, so the peak head is smaller than the Joukowsky value.

(b) Joukowsky (instantaneous) pressure rise

Δp=ρaV0=1000×1200×2.5=3.0×106 Pa=3.0 MPa\Delta p = \rho\,a\,V_0 = 1000\times1200\times2.5 = 3.0\times10^{6}\ \text{Pa} = \boxed{3.0\ \text{MPa}}

As head:

ΔH=aV0g=1200×2.59.81=30009.81=305.81 m  ΔH305.8 m\Delta H = \frac{a V_0}{g} = \frac{1200\times2.5}{9.81} = \frac{3000}{9.81} = 305.81\ \text{m}\;\Rightarrow\boxed{\Delta H \approx 305.8\ \text{m}}

(c) Closure in tc=3.0 st_c = 3.0\ \text{s}

Since tc=3.0 s>T=2.0 st_c = 3.0\ \text{s} > T = 2.0\ \text{s}, the closure is slow (gradual).

Using the (Michaud/Allievi) approximation for slow uniform closure:

ΔH=2LV0gtc=2×1200×2.59.81×3.0=600029.43=203.87 m\Delta H = \frac{2 L V_0}{g\,t_c} = \frac{2\times1200\times2.5}{9.81\times3.0} = \frac{6000}{29.43} = 203.87\ \text{m} ΔH203.9 m\boxed{\Delta H \approx 203.9\ \text{m}}

This is well below the Joukowsky value, confirming the benefit of slower closure. (Equivalently ΔH=ΔHmaxT/tc=305.8×2.0/3.0=203.9 m\Delta H = \Delta H_{max}\cdot T/t_c = 305.8\times2.0/3.0 = 203.9\ \text{m}.)

water-hammersurgespressure-rise
7short5 marks

Using the direct step method, compute one step of a backwater profile in a wide rectangular channel with q=3.0 m2/sq = 3.0\ \text{m}^2/\text{s}, S0=0.0005S_0 = 0.0005, n=0.025n = 0.025. The flow depth increases from y1=2.0 my_1 = 2.0\ \text{m} to y2=2.2 my_2 = 2.2\ \text{m}. Find the distance Δx\Delta x between these two sections.

Given (wide channel, RyR\approx y): q=3.0 m2/sq = 3.0\ \text{m}^2/\text{s}, S0=0.0005S_0 = 0.0005, n=0.025n = 0.025, g=9.81g=9.81.

Direct step: Δx=E2E1S0Sˉf\Delta x = \dfrac{E_2 - E_1}{S_0 - \bar S_f}, where E=y+q22gy2E = y + \dfrac{q^2}{2g y^2} and Sˉf\bar S_f is the mean friction slope.

Specific energy at each section (q22g=9.019.62=0.458716\dfrac{q^2}{2g} = \dfrac{9.0}{19.62} = 0.458716):

Sectionyy (m)V=q/yV=q/y (m/s)V2/2gV^2/2g (m)EE (m)
12.01.50000.458716/4.0=0.1146790.458716/4.0=0.1146792.114679
22.21.36360.458716/4.84=0.0947760.458716/4.84=0.0947762.294776

E2E1=2.2947762.114679=0.180097 mE_2 - E_1 = 2.294776 - 2.114679 = 0.180097\ \text{m}.

Friction slope Sf=n2V2R4/3=n2q2y10/3S_f = \dfrac{n^2 V^2}{R^{4/3}} = \dfrac{n^2 q^2}{y^{10/3}} (wide channel, R=yR=y).

n2q2=0.0252×3.02=0.000625×9.0=0.005625n^2 q^2 = 0.025^2\times3.0^2 = 0.000625\times9.0 = 0.005625.

Section 1: y10/3=2.010/3y^{10/3} = 2.0^{10/3}. ln2=0.693147\ln2=0.693147, ×10/3=2.310490\times10/3=2.310490, e2.310490=10.0794e^{2.310490}=10.0794. Sf1=0.005625/10.0794=5.5807×104S_{f1} = 0.005625/10.0794 = 5.5807\times10^{-4}.

Section 2: y10/3=2.210/3y^{10/3} = 2.2^{10/3}. ln2.2=0.788457\ln2.2=0.788457, ×10/3=2.628190\times10/3=2.628190, e2.628190=13.8497e^{2.628190}=13.8497. Sf2=0.005625/13.8497=4.0614×104S_{f2} = 0.005625/13.8497 = 4.0614\times10^{-4}.

Mean: Sˉf=12(5.5807+4.0614)×104=4.82105×104\bar S_f = \tfrac12(5.5807+4.0614)\times10^{-4} = 4.82105\times10^{-4}.

Step length:

Δx=E2E1S0Sˉf=0.1800970.00050.000482105=0.1800970.0000178950\Delta x = \frac{E_2 - E_1}{S_0 - \bar S_f} = \frac{0.180097}{0.0005 - 0.000482105} = \frac{0.180097}{0.0000178950} Δx=10063.8 m  Δx1.01×104 m (10.06 km)\Delta x = 10063.8\ \text{m}\;\Rightarrow\boxed{\Delta x \approx 1.01\times10^{4}\ \text{m}\ (\approx 10.06\ \text{km})}

The positive Δx\Delta x (with depth rising and Sf<S0S_f<S_0, i.e. y>yny>y_n) indicates an M1 backwater computed in the upstream direction; the large length reflects the gentle slope and the depth being only slightly above normal.

gvf-computationdirect-step-methodbackwater
8short5 marks

A suppressed (full-width) rectangular sharp-crested weir of crest length L=3.0 mL = 3.0\ \text{m} operates under a head H=0.45 mH = 0.45\ \text{m}. Take the coefficient of discharge Cd=0.62C_d = 0.62.

(a) Compute the discharge over the weir (neglecting approach velocity). (b) Including the velocity of approach, if the upstream channel cross-section is 3.0 m×1.2 m3.0\ \text{m}\times1.2\ \text{m}, estimate the corrected discharge (one iteration). (c) State two practical uses of weirs in open-channel systems.

Given: L=3.0 mL = 3.0\ \text{m}, H=0.45 mH = 0.45\ \text{m}, Cd=0.62C_d = 0.62, g=9.81 m/s2g = 9.81\ \text{m/s}^2.

(a) Discharge (no approach velocity)

Q=23CdL2gH3/2Q = \frac{2}{3} C_d L \sqrt{2g}\, H^{3/2}

2g=19.62=4.429446\sqrt{2g} = \sqrt{19.62} = 4.429446. H3/2=0.451.5H^{3/2} = 0.45^{1.5}: 0.45=0.670820\sqrt{0.45}=0.670820, ×0.45=0.301869\times0.45 = 0.301869.

Q=23×0.62×3.0×4.429446×0.301869Q = \tfrac{2}{3}\times0.62\times3.0\times4.429446\times0.301869

23×0.62=0.413333\tfrac23\times0.62 = 0.413333; ×3.0=1.24\times3.0 = 1.24; ×4.429446=5.492513\times4.429446 = 5.492513; ×0.301869=1.65805\times0.301869 = 1.65805.

Q1.658 m3/s\boxed{Q \approx 1.658\ \text{m}^3/\text{s}}

(b) With velocity of approach (one iteration)

Approach area A=3.0×1.2=3.6 m2A = 3.0\times1.2 = 3.6\ \text{m}^2. Approach velocity Va=Q/A=1.65805/3.6=0.46057 m/sV_a = Q/A = 1.65805/3.6 = 0.46057\ \text{m/s}. Approach head ha=Va22g=0.46057219.62=0.21212519.62=0.010812 mh_a = \dfrac{V_a^2}{2g} = \dfrac{0.46057^2}{19.62} = \dfrac{0.212125}{19.62} = 0.010812\ \text{m}.

Corrected discharge:

Q=23CdL2g[(H+ha)3/2ha3/2]Q = \frac{2}{3}C_d L\sqrt{2g}\left[(H+h_a)^{3/2} - h_a^{3/2}\right]

(H+ha)=0.460812(H+h_a) = 0.460812; (0.460812)3/2(0.460812)^{3/2}: 0.460812=0.678832\sqrt{0.460812}=0.678832, ×0.460812=0.312819\times0.460812 = 0.312819. ha3/2=(0.010812)1.5h_a^{3/2} = (0.010812)^{1.5}: 0.010812=0.103981\sqrt{0.010812}=0.103981, ×0.010812=0.0011242\times0.010812 = 0.0011242. Bracket =0.3128190.0011242=0.311695= 0.312819 - 0.0011242 = 0.311695.

Q=5.492513×0.311695=1.71200 m3/sQ = 5.492513\times0.311695 = 1.71200\ \text{m}^3/\text{s} Q1.712 m3/s\boxed{Q \approx 1.712\ \text{m}^3/\text{s}}

(about 3.3%3.3\% higher than the no-approach value).

(c) Two uses of weirs

  1. Flow measurement — a calibrated weir gives discharge from a single head reading HH.
  2. Water-level / flow regulation — weirs raise and maintain upstream water levels for diversion, head regulation at canal intakes, and act as overflow/escape structures (e.g., spillways, drop structures) to pass surplus flow safely.
hydraulic-structuresweirsdischarge-measurement
9short5 marks

A rectangular canal 4.0 m4.0\ \text{m} wide carries a steady flow of depth 1.0 m1.0\ \text{m} at a velocity of 1.5 m/s1.5\ \text{m/s}. Due to sudden partial gate closure downstream, the depth just upstream of the gate rises to 1.6 m1.6\ \text{m}, producing a positive surge that travels upstream.

(a) Compute the absolute speed (celerity) of the surge. (b) Compute the new (reduced) flow velocity behind the surge. (c) Comment on whether the surge is a moving hydraulic jump.

Given: b=4.0 mb = 4.0\ \text{m}, y1=1.0 my_1 = 1.0\ \text{m} (initial), V1=1.5 m/sV_1 = 1.5\ \text{m/s}, y2=1.6 my_2 = 1.6\ \text{m} (after surge), g=9.81 m/s2g = 9.81\ \text{m/s}^2. Positive surge moving upstream with absolute speed VwV_w (take upstream as the surge direction; flow is downstream).

Let the surge move upstream with absolute speed CC (magnitude). Use control-volume relations on the moving surge. For a surge advancing upstream against a flow, with Vw=CV_w = C (speed of surge), the relations are:

Continuity (relative to surge):

(V1+Vw)y1=(V2+Vw)y2(V_1 + V_w)y_1 = (V_2 + V_w)y_2

Momentum / surge celerity (positive surge): relative approach velocity equals

(V1+Vw)=gy22y1(y1+y2)(V_1 + V_w) = \sqrt{\frac{g\,y_2}{2 y_1}(y_1 + y_2)}

(a) Surge speed

V1+Vw=9.81×1.62×1.0(1.0+1.6)=15.6962.0×2.6=7.848×2.6=20.4048=4.51717 m/sV_1 + V_w = \sqrt{\frac{9.81\times1.6}{2\times1.0}(1.0+1.6)} = \sqrt{\frac{15.696}{2.0}\times2.6} = \sqrt{7.848\times2.6} = \sqrt{20.4048} = 4.51717\ \text{m/s} Vw=4.51717V1=4.517171.5=3.0172 m/sV_w = 4.51717 - V_1 = 4.51717 - 1.5 = 3.0172\ \text{m/s} Vw3.02 m/s (upstream)\boxed{V_w \approx 3.02\ \text{m/s (upstream)}}

(b) Velocity behind the surge From continuity:

V2+Vw=(V1+Vw)y1y2=4.51717×1.01.6=2.823231 m/sV_2 + V_w = \frac{(V_1+V_w)\,y_1}{y_2} = \frac{4.51717\times1.0}{1.6} = 2.823231\ \text{m/s} V2=2.823231Vw=2.8232313.0172=0.1940 m/sV_2 = 2.823231 - V_w = 2.823231 - 3.0172 = -0.1940\ \text{m/s} V20.19 m/s\boxed{V_2 \approx -0.19\ \text{m/s}}

The small negative value means the water behind the surge is nearly stagnant / very slightly reversed — consistent with a strong downstream gate restriction.

Check (continuity, absolute): discharge change is absorbed by the rising level; relative discharge (V1+Vw)y1=4.51717 m2/s(V_1+V_w)y_1 = 4.51717\ \text{m}^2/\text{s} and (V2+Vw)y2=2.823231×1.6=4.51717 m2/s(V_2+V_w)y_2 = 2.823231\times1.6 = 4.51717\ \text{m}^2/\text{s} ✓.

(c) Moving hydraulic jump?

Yes. A positive surge is essentially a hydraulic jump translating along the channel. In a reference frame moving with the surge at speed VwV_w, the flow is steady and satisfies the same momentum (Belanger-type) relation as a stationary jump, with the abrupt rise from y1y_1 to y2y_2. Here the relative upstream Froude number Frr=(V1+Vw)/gy1=4.51717/9.81=4.51717/3.13209=1.442>1Fr_r = (V_1+V_w)/\sqrt{g y_1} = 4.51717/\sqrt{9.81} = 4.51717/3.13209 = 1.442 > 1, confirming supercritical relative inflow — the hallmark of a (moving) hydraulic jump.

surgesmoving-hydraulic-jumpunsteady-flow
10short5 marks

Water flows at Q=8 m3/sQ = 8\ \text{m}^3/\text{s} in a rectangular channel 2.5 m2.5\ \text{m} wide at a depth of 1.4 m1.4\ \text{m}. A smooth upward step (hump) of height Δz=0.15 m\Delta z = 0.15\ \text{m} is built across the bed.

(a) Determine the depth of flow over the hump (neglect losses). (b) Determine the maximum hump height for which the upstream depth is unaffected (the choking height).

Given: Q=8 m3/sQ = 8\ \text{m}^3/\text{s}, b=2.5 mb = 2.5\ \text{m}, y1=1.4 my_1 = 1.4\ \text{m}, Δz=0.15 m\Delta z = 0.15\ \text{m}, g=9.81g=9.81. Unit discharge q=8/2.5=3.2 m2/sq = 8/2.5 = 3.2\ \text{m}^2/\text{s}.

Upstream specific energy:

V1=q/y1=3.2/1.4=2.285714 m/s;V122g=2.285714219.62=5.22449019.62=0.266284 mV_1 = q/y_1 = 3.2/1.4 = 2.285714\ \text{m/s};\quad \frac{V_1^2}{2g} = \frac{2.285714^2}{19.62} = \frac{5.224490}{19.62} = 0.266284\ \text{m} E1=1.4+0.266284=1.666284 mE_1 = 1.4 + 0.266284 = 1.666284\ \text{m}

Flow is subcritical (y1>ycy_1 > y_c, see below).

Critical depth: yc=(q2/g)1/3=(3.22/9.81)1/3=(10.24/9.81)1/3=(1.043833)1/3y_c = (q^2/g)^{1/3} = (3.2^2/9.81)^{1/3} = (10.24/9.81)^{1/3} = (1.043833)^{1/3}. ln1.043833=0.042896\ln1.043833 = 0.042896; /3=0.014299/3 = 0.014299; e0.014299=1.014402e^{0.014299}=1.014402. So yc=1.0144 my_c = 1.0144\ \text{m}; Ec=1.5yc=1.52160 mE_c = 1.5 y_c = 1.52160\ \text{m}.

(a) Depth over hump

Energy over the hump: E2=E1Δz=1.6662840.15=1.516284 mE_2 = E_1 - \Delta z = 1.666284 - 0.15 = 1.516284\ \text{m}.

Note E2=1.51628 m>Ec=1.52160E_2 = 1.51628\ \text{m} > E_{c} = 1.52160? Compare: 1.51628<1.521601.51628 < 1.52160. This is marginally below EcE_c, which would mean choking. Re-check EcE_c: Ec=1.5×1.014402=1.521603E_c = 1.5\times1.014402 = 1.521603. Indeed E2(1.51628)<Ec(1.52160)E_2 (1.51628) < E_c (1.52160) by 0.0053 m0.0053\ \text{m} — the proposed hump just barely chokes the flow.

Because E2E_2 falls a hair below the minimum energy, strictly the flow chokes and y1y_1 must rise. For the value to be physically passable we take the limiting case: the depth over the crest is critical, y2=yc=1.014 my_2 = y_c = 1.014\ \text{m}. (If losses/idealisation place E2E_2 exactly at EcE_c, the surface passes through critical on the hump.)

y2yc=1.014 m (flow at/near critical on the hump)\boxed{y_2 \approx y_c = 1.014\ \text{m (flow at/near critical on the hump)}}

(b) Maximum (choking) hump height

Choking occurs when E2=EcE_2 = E_c, i.e.

Δzmax=E1Ec=1.6662841.521603=0.144681 m\Delta z_{max} = E_1 - E_c = 1.666284 - 1.521603 = 0.144681\ \text{m} Δzmax0.145 m\boxed{\Delta z_{max} \approx 0.145\ \text{m}}

Since the proposed Δz=0.15 m>0.145 m\Delta z = 0.15\ \text{m} > 0.145\ \text{m}, the hump exceeds the choking height: the upstream depth will increase (back up) and the flow over the crest is critical. This confirms part (a).

non-uniform-flowchannel-transitionspecific-energy
11short5 marks

(a) Derive the relation between Chezy's coefficient CC and Manning's nn. (b) A rectangular channel 3.0 m3.0\ \text{m} wide flows at depth 1.0 m1.0\ \text{m} with n=0.015n = 0.015; compute Chezy's CC. (c) Explain briefly why, in open channels, the maximum velocity occurs slightly below the free surface rather than at it.

(a) Relation between CC and nn

Chezy: V=CRSfV = C\sqrt{R S_f}. Manning: V=1nR2/3Sf1/2V = \dfrac{1}{n} R^{2/3} S_f^{1/2}. Equating the two velocity expressions:

CRSf=1nR2/3Sf1/2    CR1/2=1nR2/3C\sqrt{R S_f} = \frac{1}{n} R^{2/3} S_f^{1/2} \;\Rightarrow\; C R^{1/2} = \frac{1}{n} R^{2/3} C=1nR1/6\boxed{C = \frac{1}{n} R^{1/6}}

Thus Chezy's CC depends weakly (one-sixth power) on hydraulic radius, whereas nn is treated as constant for a given surface.

(b) Numerical value of CC

Rectangular: A=by=3.0×1.0=3.0 m2A = b y = 3.0\times1.0 = 3.0\ \text{m}^2; P=b+2y=3.0+2.0=5.0 mP = b + 2y = 3.0 + 2.0 = 5.0\ \text{m}; R=A/P=3.0/5.0=0.6 mR = A/P = 3.0/5.0 = 0.6\ \text{m}. R1/6=0.61/6R^{1/6} = 0.6^{1/6}: ln0.6=0.510826\ln0.6 = -0.510826; /6=0.085138/6 = -0.085138; e0.085138=0.918386e^{-0.085138} = 0.918386.

C=10.015×0.918386=66.6667×0.918386=61.226 m1/2/sC = \frac{1}{0.015}\times0.918386 = 66.6667\times0.918386 = 61.226\ \text{m}^{1/2}/\text{s} C61.2 m1/2/s\boxed{C \approx 61.2\ \text{m}^{1/2}/\text{s}}

(c) Location of maximum velocity

In an open channel the free surface is subject to secondary currents and air–water interaction (surface resistance, wind/surface tension and the downward-directed secondary flow cells). These secondary (cross-circulation) currents carry low-momentum fluid from the corners and sides toward the centre of the surface and push high-momentum fluid downward, so the point of maximum velocity is depressed to roughly 0.050.050.250.25 of the depth below the surface (commonly about 0.2y0.2 y below it for wide channels). Consequently the mean velocity occurs near 0.60.6 of the depth from the surface, and the average of velocities at 0.2y0.2 y and 0.8y0.8 y approximates the mean — the basis of the standard velocity-measurement rules.

chezy-manningvelocity-distributionconcepts

Frequently asked questions

Where can I find the BE Civil Engineering (IOE, TU) Hydraulics (IOE, CE 604) question paper 2077?
The full BE Civil Engineering (IOE, TU) Hydraulics (IOE, CE 604) 2077 (regular) question paper is available free on Kekkei. You can read every question online and attempt the paper under timed exam conditions.
Does the Hydraulics (IOE, CE 604) 2077 paper come with solutions?
Yes. Every question on this Hydraulics (IOE, CE 604) past paper includes a step-by-step solution, plus instant AI feedback when you attempt it on Kekkei.
How many marks is the BE Civil Engineering (IOE, TU) Hydraulics (IOE, CE 604) 2077 paper?
The BE Civil Engineering (IOE, TU) Hydraulics (IOE, CE 604) 2077 paper carries 80 full marks and is meant to be completed in 180 minutes, across 11 questions.
Is practising this Hydraulics (IOE, CE 604) past paper free?
Yes — reading and attempting this Hydraulics (IOE, CE 604) past paper on Kekkei is completely free.