Browse papers
A

Section A: Long Answer Questions

Attempt all questions.

5 questions
1long10 marks

A trapezoidal channel carries a discharge of Q=18 m3/sQ = 18\ \text{m}^3/\text{s}. The channel has a bed width b=4.0 mb = 4.0\ \text{m}, side slopes of 1.5 H:1 V1.5\ \text{H}:1\ \text{V} (i.e. z=1.5z = 1.5), a longitudinal bed slope S0=0.0009S_0 = 0.0009 and Manning's roughness coefficient n=0.018n = 0.018.

(a) Derive Manning's equation for uniform flow starting from the Chezy equation, clearly stating the assumptions.

(b) Determine the normal depth of flow yny_n in the channel by trial. Carry the trial to a depth accurate to about 0.01 m0.01\ \text{m} and state the mean velocity at normal depth.

(a) Derivation of Manning's equation

Uniform flow means the depth, area, velocity and discharge are constant along the channel; the water surface, energy line and bed are all parallel, so the slope of the energy line equals the bed slope, Sf=S0=SS_f = S_0 = S.

Chezy's equation for uniform flow is:

V=CRSV = C\sqrt{R S}

where VV = mean velocity, CC = Chezy coefficient, R=A/PR = A/P = hydraulic radius, SS = slope.

Manning (and Strickler) found experimentally that the Chezy coefficient varies with the roughness and the hydraulic radius as:

C=1nR1/6C = \frac{1}{n} R^{1/6}

where nn is Manning's roughness coefficient. Substituting into Chezy's equation:

V=1nR1/6RS=1nR2/3S1/2V = \frac{1}{n} R^{1/6}\sqrt{R S} = \frac{1}{n} R^{2/3} S^{1/2}

Multiplying by area AA gives the discharge form:

Q=1nAR2/3S1/2\boxed{Q = \frac{1}{n} A R^{2/3} S^{1/2}}

Assumptions: flow is steady and uniform; channel is prismatic with constant slope and roughness; pressure distribution is hydrostatic; Sf=S0S_f = S_0.

(b) Normal depth by trial

For a trapezoidal section with bed width b=4.0 mb = 4.0\ \text{m} and side slope z=1.5z = 1.5:

A=(b+zy)y=(4+1.5y)yA = (b + z y)\,y = (4 + 1.5y)\,y P=b+2y1+z2=4+2y1+1.52=4+2y(1.8028)=4+3.6056yP = b + 2y\sqrt{1+z^2} = 4 + 2y\sqrt{1+1.5^2} = 4 + 2y\,(1.8028) = 4 + 3.6056\,y R=A/PR = A/P

Manning's equation must satisfy:

Qn=AR2/3S1/2Q n = A R^{2/3} S^{1/2} Qn=18×0.018=0.324Q n = 18 \times 0.018 = 0.324 S1/2=0.0009=0.03S^{1/2} = \sqrt{0.0009} = 0.03

so the target is:

AR2/3=QnS1/2=0.3240.03=10.80A R^{2/3} = \frac{Q n}{S^{1/2}} = \frac{0.324}{0.03} = 10.80

Trial table (compute AR2/3A R^{2/3} and compare with 10.80):

yy (m)A=(4+1.5y)yA=(4+1.5y)yP=4+3.6056yP=4+3.6056yR=A/PR=A/PR2/3R^{2/3}AR2/3A R^{2/3}
1.509.3759.4080.99650.99779.353
1.7011.13510.1301.09921.065611.866
1.6210.4189.8411.05861.038810.822
1.61910.4099.8381.05811.038510.810

At y=1.619 my = 1.619\ \text{m}, AR2/3=10.8110.80A R^{2/3} = 10.81 \approx 10.80. Hence:

yn1.62 m\boxed{y_n \approx 1.62\ \text{m}}

Mean velocity at normal depth (using yn=1.62y_n = 1.62, A=10.418 m2A = 10.418\ \text{m}^2):

V=QA=1810.418=1.728 m/sV = \frac{Q}{A} = \frac{18}{10.418} = 1.728\ \text{m/s} V1.73 m/s\boxed{V \approx 1.73\ \text{m/s}}
uniform-flowmanning-equationtrapezoidal-channel
2long8 marks

A rectangular channel 5.0 m5.0\ \text{m} wide carries a discharge of 15 m3/s15\ \text{m}^3/\text{s} at a depth of 0.60 m0.60\ \text{m}.

(a) Define specific energy and sketch the specific-energy curve, indicating the critical depth and the alternate depths.

(b) For the given flow, compute the specific energy, the critical depth, the Froude number and state whether the flow is sub- or supercritical.

(c) Find the alternate depth corresponding to the same specific energy and discharge.

(a) Specific energy

Specific energy EE is the energy per unit weight of water measured with respect to the channel bed:

E=y+V22g=y+Q22gA2E = y + \frac{V^2}{2g} = y + \frac{Q^2}{2 g A^2}

The specific-energy curve (E on x-axis, y on y-axis for fixed Q) is a curve with two limbs that approaches the line y=Ey = E (asymptote) for large depths and the EE-axis for small depths. The two limbs meet at the critical depth ycy_c, where EE is minimum. For any E>EminE > E_{min} there are two depths — the alternate depths y1y_1 (supercritical, lower) and y2y_2 (subcritical, upper):

  y |        / (y = E asymptote)
    |       /  .  subcritical limb
    |      /.
 y_c|----- *  (critical point, E = E_min)
    |    .
    |  .   supercritical limb
    |._______________ E
        E_min

(b) Computations

Unit discharge: q=Q/b=15/5=3.0 m2/sq = Q/b = 15/5 = 3.0\ \text{m}^2/\text{s}.

Velocity at y=0.60 my = 0.60\ \text{m}:

V=qy=3.00.60=5.0 m/sV = \frac{q}{y} = \frac{3.0}{0.60} = 5.0\ \text{m/s}

Specific energy:

E=y+V22g=0.60+5.022(9.81)=0.60+1.274=1.874 mE = y + \frac{V^2}{2g} = 0.60 + \frac{5.0^2}{2(9.81)} = 0.60 + 1.274 = 1.874\ \text{m} E=1.87 m\boxed{E = 1.87\ \text{m}}

Critical depth (rectangular):

yc=(q2g)1/3=(3.029.81)1/3=(0.9174)1/3=0.972 my_c = \left(\frac{q^2}{g}\right)^{1/3} = \left(\frac{3.0^2}{9.81}\right)^{1/3} = (0.9174)^{1/3} = 0.972\ \text{m} yc=0.972 m\boxed{y_c = 0.972\ \text{m}}

Froude number at y=0.60y = 0.60:

Fr=Vgy=5.09.81×0.60=5.02.426=2.06Fr = \frac{V}{\sqrt{g y}} = \frac{5.0}{\sqrt{9.81 \times 0.60}} = \frac{5.0}{2.426} = 2.06

Since Fr=2.06>1Fr = 2.06 > 1 (and y<ycy < y_c), the flow is supercritical.

(c) Alternate depth

The alternate depth y2y_2 satisfies the same EE and qq:

E=y2+q22gy22=1.874y2+919.62y22=1.874E = y_2 + \frac{q^2}{2 g y_2^2} = 1.874 \Rightarrow y_2 + \frac{9}{19.62\, y_2^2} = 1.874 y2+0.4587y22=1.874y_2 + \frac{0.4587}{y_2^2} = 1.874

Solve by trial for the subcritical root (y2>ycy_2 > y_c):

y2y_2 (m)y2+0.4587/y22y_2 + 0.4587/y_2^2
1.701.859
1.721.875
1.7191.875

At y2=1.72 my_2 = 1.72\ \text{m} the expression equals 1.8751.8741.875 \approx 1.874.

y21.72 m (alternate subcritical depth)\boxed{y_2 \approx 1.72\ \text{m (alternate subcritical depth)}}
specific-energycritical-flowalternate-depths
3long8 marks

A hydraulic jump forms in a horizontal rectangular channel 6.0 m6.0\ \text{m} wide. The supercritical depth before the jump is 0.45 m0.45\ \text{m} and the discharge is 24 m3/s24\ \text{m}^3/\text{s}.

(a) Derive the sequent-depth (conjugate-depth) relation for a hydraulic jump in a rectangular channel using the momentum principle.

(b) Compute the sequent depth, the energy loss in the jump, the efficiency, and the length of the jump (use Lj6y2L_j \approx 6 y_2).

(a) Sequent-depth relation

Applying the momentum equation between sections 1 (before) and 2 (after) a jump on a horizontal frictionless bed, the specific force (momentum function) is conserved:

Q2gA1+zˉ1A1=Q2gA2+zˉ2A2\frac{Q^2}{g A_1} + \bar{z}_1 A_1 = \frac{Q^2}{g A_2} + \bar{z}_2 A_2

For a rectangular channel of width bb, using q=Q/bq = Q/b, A=byA = by, centroid depth zˉ=y/2\bar z = y/2, dividing by bb:

q2gy1+y122=q2gy2+y222\frac{q^2}{g y_1} + \frac{y_1^2}{2} = \frac{q^2}{g y_2} + \frac{y_2^2}{2}

Rearranging and factoring leads to:

y2y1=12(1+8Fr121)\frac{y_2}{y_1} = \frac{1}{2}\left(\sqrt{1 + 8 Fr_1^2} - 1\right)

where Fr1=V1/gy1Fr_1 = V_1/\sqrt{g y_1} is the upstream Froude number. This is the sequent-depth relation.

(b) Computations

Unit discharge: q=Q/b=24/6=4.0 m2/sq = Q/b = 24/6 = 4.0\ \text{m}^2/\text{s}.

Upstream velocity:

V1=qy1=4.00.45=8.889 m/sV_1 = \frac{q}{y_1} = \frac{4.0}{0.45} = 8.889\ \text{m/s}

Upstream Froude number:

Fr1=V1gy1=8.8899.81×0.45=8.8892.101=4.231Fr_1 = \frac{V_1}{\sqrt{g y_1}} = \frac{8.889}{\sqrt{9.81 \times 0.45}} = \frac{8.889}{2.101} = 4.231

Sequent depth:

y2=y12(1+8Fr121)=0.452(1+8(4.231)21)y_2 = \frac{y_1}{2}\left(\sqrt{1 + 8 Fr_1^2} - 1\right) = \frac{0.45}{2}\left(\sqrt{1 + 8(4.231)^2} - 1\right) 8Fr12=8×17.90=143.2144.2=12.0088 Fr_1^2 = 8 \times 17.90 = 143.2 \Rightarrow \sqrt{144.2} = 12.008 y2=0.225×(12.0081)=0.225×11.008=2.477 my_2 = 0.225 \times (12.008 - 1) = 0.225 \times 11.008 = 2.477\ \text{m} y22.48 m\boxed{y_2 \approx 2.48\ \text{m}}

Energy loss in a rectangular jump:

ΔE=(y2y1)34y1y2=(2.4770.45)34(0.45)(2.477)=(2.027)34.459=8.3284.459=1.868 m\Delta E = \frac{(y_2 - y_1)^3}{4 y_1 y_2} = \frac{(2.477 - 0.45)^3}{4 (0.45)(2.477)} = \frac{(2.027)^3}{4.459} = \frac{8.328}{4.459} = 1.868\ \text{m} ΔE1.87 m\boxed{\Delta E \approx 1.87\ \text{m}}

Energy before the jump:

E1=y1+V122g=0.45+8.889219.62=0.45+4.027=4.477 mE_1 = y_1 + \frac{V_1^2}{2g} = 0.45 + \frac{8.889^2}{19.62} = 0.45 + 4.027 = 4.477\ \text{m}

Efficiency (ratio of energy after to energy before):

E2=E1ΔE=4.4771.868=2.609 mE_2 = E_1 - \Delta E = 4.477 - 1.868 = 2.609\ \text{m} η=E2E1=2.6094.477=0.583=58.3%\eta = \frac{E_2}{E_1} = \frac{2.609}{4.477} = 0.583 = 58.3\% η58%\boxed{\eta \approx 58\%}

Fraction of energy dissipated =ΔE/E1=1.868/4.477=41.7%= \Delta E / E_1 = 1.868/4.477 = 41.7\%.

Length of jump:

Lj6y2=6×2.477=14.86 mL_j \approx 6 y_2 = 6 \times 2.477 = 14.86\ \text{m} Lj14.9 m\boxed{L_j \approx 14.9\ \text{m}}

With Fr1=4.23Fr_1 = 4.23 this is a steady jump (well-formed, 4.5>Fr1>2.54.5 > Fr_1 > 2.5 region near steady).

hydraulic-jumpenergy-dissipationsequent-depth
4long8 marks

(a) Derive the dynamic equation of gradually varied flow (GVF) dydx=S0Sf1Fr2\dfrac{dy}{dx} = \dfrac{S_0 - S_f}{1 - Fr^2}, stating the assumptions involved.

(b) Classify and sketch the possible GVF water-surface profiles on a mild slope (M1, M2, M3). For each, indicate the depth ranges relative to normal depth yny_n and critical depth ycy_c and give one physical example.

(a) Dynamic equation of GVF

Assumptions: steady flow; pressure distribution hydrostatic (gradual variation, streamlines nearly straight and parallel); the head loss at a section equals that of uniform flow at the same depth and velocity (so SfS_f from Manning/Chezy applies); slope small (cosθ1\cos\theta \approx 1); channel prismatic; α=1\alpha = 1.

Total energy referred to a horizontal datum:

H=z+y+V22gH = z + y + \frac{V^2}{2g}

Differentiate with respect to xx (channel distance):

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 line slope) and dzdx=S0\dfrac{dz}{dx} = -S_0 (bed slope). Also for a given QQ:

ddx(V22g)=ddy(Q22gA2)dydx=Q2TgA3dydx=Fr2dydx\frac{d}{dx}\left(\frac{V^2}{2g}\right) = \frac{d}{dy}\left(\frac{Q^2}{2gA^2}\right)\frac{dy}{dx} = -\frac{Q^2 T}{g A^3}\frac{dy}{dx} = -Fr^2\frac{dy}{dx}

since Fr2=Q2TgA3Fr^2 = \dfrac{Q^2 T}{g A^3}.

Substituting:

Sf=S0+dydxFr2dydx-S_f = -S_0 + \frac{dy}{dx} - Fr^2 \frac{dy}{dx} S0Sf=dydx(1Fr2)S_0 - S_f = \frac{dy}{dx}\left(1 - Fr^2\right) dydx=S0Sf1Fr2\boxed{\frac{dy}{dx} = \frac{S_0 - S_f}{1 - Fr^2}}

(b) Mild-slope profiles (M-profiles)

On a mild slope yn>ycy_n > y_c. Three zones are defined by yny_n and ycy_c:

ProfileDepth rangedy/dxdy/dxBehaviourExample
M1y>yn>ycy > y_n > y_c++ (backwater)depth increases d/s, surface rises asymptotic to horizontalBackwater behind a dam / weir on a mild river
M2yn>y>ycy_n > y > y_c- (drawdown)depth falls toward ycy_cUpstream of a free overfall / sudden drop
M3yn>yc>yy_n > y_c > y++supercritical, depth rises toward ycy_c then jumpDownstream of a sluice gate on a mild bed
          M1 (y>yn) ---- rising backwater
  ----------------------------------  y = yn (NDL)
          M2 (yc<y<yn) -- drawdown
  - - - - - - - - - - - - - - - - -  y = yc (CDL)
          M3 (y<yc) ---- rising to jump
  ////////// mild bed //////////////

M1 is asymptotic to NDL upstream and horizontal downstream; M2 is asymptotic to NDL upstream and meets CDL nearly vertically; M3 starts near the bed (gate) and rises, usually terminated by a hydraulic jump.

gradually-varied-flowflow-profilesclassification
5long8 marks

A rectangular channel 8.0 m8.0\ \text{m} wide carries 20 m3/s20\ \text{m}^3/\text{s} with n=0.025n = 0.025 on a bed slope S0=0.0004S_0 = 0.0004. A weir raises the depth at a control section to 2.50 m2.50\ \text{m}.

(a) Verify that the channel slope is mild and identify the resulting profile type.

(b) Using the direct-step method, compute the distance (in one step) from the control section (y=2.50 my = 2.50\ \text{m}) to the section where the depth is 2.30 m2.30\ \text{m}.

(a) Slope classification

Unit discharge q=20/8=2.5 m2/sq = 20/8 = 2.5\ \text{m}^2/\text{s}.

Critical depth:

yc=(q2g)1/3=(2.529.81)1/3=(0.6371)1/3=0.860 my_c = \left(\frac{q^2}{g}\right)^{1/3} = \left(\frac{2.5^2}{9.81}\right)^{1/3} = (0.6371)^{1/3} = 0.860\ \text{m}

Normal depth from Manning (Qn/S1/2=AR2/3Q n / S^{1/2} = A R^{2/3}):

QnS1/2=20×0.0250.0004=0.50.02=25.0\frac{Qn}{S^{1/2}} = \frac{20 \times 0.025}{\sqrt{0.0004}} = \frac{0.5}{0.02} = 25.0

For b=8b=8, A=8yA = 8y, P=8+2yP = 8 + 2y, R=8y/(8+2y)R = 8y/(8+2y). Trial:

yy (m)AAPPRRR2/3R^{2/3}AR2/3A R^{2/3}
1.8014.4011.601.2411.16016.70
2.3018.4012.601.4601.29023.74
2.4019.2012.801.5001.31025.16
2.3919.1212.781.4961.30825.01

So yn2.39 my_n \approx 2.39\ \text{m}. Since yn(2.39)>yc(0.86)y_n (2.39) > y_c (0.86), the slope is mild.

The control depth 2.50 m>yn>yc2.50\ \text{m} > y_n > y_c, so the profile is an M1 (backwater) curve.

(b) Direct-step method (one step, y1=2.50y2=2.30y_1=2.50 \to y_2=2.30)

Direct-step uses:

Δx=E2E1S0Sfˉ\Delta x = \frac{E_2 - E_1}{S_0 - \bar{S_f}}

where E=y+V2/2gE = y + V^2/2g and Sˉf\bar S_f is the mean of SfS_f at the two sections, with Sf=(QnAR2/3)2S_f = \left(\dfrac{Qn}{A R^{2/3}}\right)^2.

Section 1: y1=2.50 my_1 = 2.50\ \text{m}

A1=8(2.50)=20.0 m2,V1=20/20=1.0 m/sA_1 = 8(2.50) = 20.0\ \text{m}^2,\quad V_1 = 20/20 = 1.0\ \text{m/s} E1=2.50+1.0219.62=2.50+0.0510=2.5510 mE_1 = 2.50 + \frac{1.0^2}{19.62} = 2.50 + 0.0510 = 2.5510\ \text{m} P1=8+2(2.5)=13.0, R1=20/13=1.538, R12/3=1.334P_1 = 8 + 2(2.5) = 13.0,\ R_1 = 20/13 = 1.538,\ R_1^{2/3} = 1.334 Sf1=(20×0.02520.0×1.334)2=(0.526.68)2=(0.018740)2=3.512×104S_{f1} = \left(\frac{20 \times 0.025}{20.0 \times 1.334}\right)^2 = \left(\frac{0.5}{26.68}\right)^2 = (0.018740)^2 = 3.512\times10^{-4}

Section 2: y2=2.30 my_2 = 2.30\ \text{m}

A2=8(2.30)=18.40 m2,V2=20/18.40=1.0870 m/sA_2 = 8(2.30) = 18.40\ \text{m}^2,\quad V_2 = 20/18.40 = 1.0870\ \text{m/s} E2=2.30+1.0870219.62=2.30+0.06022=2.3602 mE_2 = 2.30 + \frac{1.0870^2}{19.62} = 2.30 + 0.06022 = 2.3602\ \text{m} P2=8+2(2.30)=12.60, R2=18.40/12.60=1.460, R22/3=1.290P_2 = 8 + 2(2.30) = 12.60,\ R_2 = 18.40/12.60 = 1.460,\ R_2^{2/3} = 1.290 Sf2=(0.518.40×1.290)2=(0.523.74)2=(0.021062)2=4.436×104S_{f2} = \left(\frac{0.5}{18.40 \times 1.290}\right)^2 = \left(\frac{0.5}{23.74}\right)^2 = (0.021062)^2 = 4.436\times10^{-4}

Mean friction slope:

Sˉf=3.512+4.4362×104=3.974×104\bar S_f = \frac{3.512 + 4.436}{2}\times10^{-4} = 3.974\times10^{-4}

Distance:

Δx=E2E1S0Sˉf=2.36022.55100.00040.0003974=0.19082.6×106\Delta x = \frac{E_2 - E_1}{S_0 - \bar S_f} = \frac{2.3602 - 2.5510}{0.0004 - 0.0003974} = \frac{-0.1908}{2.6\times10^{-6}} Δx=73,400 m (approx)\Delta x = -73{,}400\ \text{m (approx)}

The large magnitude arises because at yy near yn=2.39y_n=2.39 the denominator S0SˉfS_0-\bar S_f is nearly zero (the M1 profile is asymptotic to the NDL). The negative sign means section 2 lies upstream of the control. A practical computation would use small steps; the magnitude shows the backwater extends a very long distance upstream.

Δx73 km upstream (asymptotic M1 behaviour near yn)\boxed{|\Delta x| \approx 73\ \text{km upstream (asymptotic M1 behaviour near } y_n)}

Note: taking the step only to 2.45 m2.45\ \text{m} (well above yny_n) would give a finite, modest distance; this part illustrates why direct-step requires depths kept away from yny_n.

gvf-computationdirect-step-methodbackwater
B

Section B: Short Answer Questions

Attempt all questions.

6 questions
6short6 marks

(a) State the conditions for the most economical (best hydraulic) trapezoidal section.

(b) Design the most economical trapezoidal channel section to carry Q=10 m3/sQ = 10\ \text{m}^3/\text{s} with side slope z=1z = 1 (45°), n=0.020n = 0.020 and S0=0.0006S_0 = 0.0006. Find the bed width, depth and check the velocity.

(a) Conditions for the most economical trapezoidal section

For a given area, the section is most economical when the wetted perimeter (and hence cost/lining) is minimum, which for Manning also maximises discharge. The conditions are:

  1. Half the top width equals the sloping side length: b+2zy2=y1+z2\dfrac{b + 2zy}{2} = y\sqrt{1+z^2}.
  2. The hydraulic radius equals half the depth: R=y/2R = y/2.
  3. A semicircle drawn with centre at the water surface and radius equal to the depth is tangent to the three sides (the section can inscribe a semicircle of radius yy).

(b) Design

Condition 1 gives the bed width for z=1z = 1:

b=2y(1+z2z)=2y(21)=2y(0.41421)=0.8284yb = 2y\left(\sqrt{1+z^2} - z\right) = 2y\left(\sqrt{2} - 1\right) = 2y(0.41421) = 0.8284\,y

Geometry:

A=(b+zy)y=(0.8284y+y)y=1.8284y2A = (b + zy)y = (0.8284y + y)y = 1.8284\,y^2

Using condition 2, R=y/2R = y/2.

Manning:

Q=1nAR2/3S1/2Q = \frac{1}{n} A R^{2/3} S^{1/2} 10=10.020(1.8284y2)(y2)2/3(0.0006)1/210 = \frac{1}{0.020}\,(1.8284 y^2)\left(\frac{y}{2}\right)^{2/3}(0.0006)^{1/2}

Compute constants: (1/2)2/3=0.6300(1/2)^{2/3} = 0.6300; 0.0006=0.024495\sqrt{0.0006} = 0.024495; 1/0.020=501/0.020 = 50.

10=50×1.8284×0.6300×0.024495×y2y2/310 = 50 \times 1.8284 \times 0.6300 \times 0.024495 \times y^{2}\,y^{2/3} 10=50×1.8284×0.6300×0.024495×y8/310 = 50 \times 1.8284 \times 0.6300 \times 0.024495 \times y^{8/3}

Coefficient =50×1.8284=91.42; ×0.6300=57.59; ×0.024495=1.4107= 50 \times 1.8284 = 91.42;\ \times 0.6300 = 57.59;\ \times 0.024495 = 1.4107.

10=1.4107y8/3y8/3=7.08910 = 1.4107\,y^{8/3} \Rightarrow y^{8/3} = 7.089 y=7.0893/8=7.0890.375y = 7.089^{3/8} = 7.089^{0.375} ln7.089=1.9585, ×0.375=0.7344, e0.7344=2.084\ln 7.089 = 1.9585,\ \times 0.375 = 0.7344,\ e^{0.7344} = 2.084 y=2.08 m\boxed{y = 2.08\ \text{m}}

Bed width:

b=0.8284×2.084=1.726 m1.73 mb = 0.8284 \times 2.084 = 1.726\ \text{m} \approx \boxed{1.73\ \text{m}}

Velocity check:

A=1.8284×2.0842=1.8284×4.343=7.941 m2A = 1.8284 \times 2.084^2 = 1.8284 \times 4.343 = 7.941\ \text{m}^2 V=QA=107.941=1.259 m/sV = \frac{Q}{A} = \frac{10}{7.941} = 1.259\ \text{m/s} V1.26 m/s\boxed{V \approx 1.26\ \text{m/s}}

This lies in the non-silting, non-scouring range for a lined/firm channel, so the design is acceptable.

channel-designbest-hydraulic-sectioneconomical-section
7short6 marks

Water flows at 2.5 m/s2.5\ \text{m/s} in a rigid steel pipeline 1200 m1200\ \text{m} long. The pressure wave (celerity) in the pipe is 1100 m/s1100\ \text{m/s}. A valve at the downstream end is operated.

(a) Distinguish between rapid (sudden) and gradual valve closure for water hammer.

(b) Compute the maximum pressure rise for instantaneous closure. (c) If the valve is closed in 3.0 s3.0\ \text{s}, determine whether the closure is rapid or gradual and estimate the pressure rise for that case. Take ρ=1000 kg/m3\rho = 1000\ \text{kg/m}^3.

(a) Rapid vs gradual closure

The critical (pipe-period) time is:

Tc=2LaT_c = \frac{2L}{a}

where LL = pipe length and aa = wave celerity. It is the time for the pressure wave to travel to the reservoir and back.

  • Rapid / sudden closure: valve closing time TTcT \le T_c. The full Joukowsky pressure rise develops because the reflected wave has not yet returned: Δp=ρaV0\Delta p = \rho a \,V_0.
  • Gradual closure: T>TcT > T_c. The reflected (relief) wave returns before closure is complete, so the pressure rise is smaller than the Joukowsky value and is estimated from rigid-column or interpolation formulae.

(b) Maximum (instantaneous) pressure rise

Joukowsky equation:

Δp=ρaV0=1000×1100×2.5=2.75×106 Pa\Delta p = \rho\, a\, V_0 = 1000 \times 1100 \times 2.5 = 2.75 \times 10^6\ \text{Pa} Δp=2.75 MPa\boxed{\Delta p = 2.75\ \text{MPa}}

As head:

ΔH=Δpρg=2.75×1061000×9.81=280.3 m\Delta H = \frac{\Delta p}{\rho g} = \frac{2.75\times10^6}{1000 \times 9.81} = 280.3\ \text{m} ΔH280 m of water\boxed{\Delta H \approx 280\ \text{m of water}}

(c) Closure in T=3.0 sT = 3.0\ \text{s}

Critical time:

Tc=2La=2×12001100=24001100=2.182 sT_c = \frac{2L}{a} = \frac{2 \times 1200}{1100} = \frac{2400}{1100} = 2.182\ \text{s}

Since T=3.0 s>Tc=2.18 sT = 3.0\ \text{s} > T_c = 2.18\ \text{s}, the closure is gradual.

For gradual closure, estimate the pressure rise from the rigid-column / Michaud (Allievi linear) formula:

ΔH=2LV0gT=2×1200×2.59.81×3.0=600029.43=203.9 m\Delta H = \frac{2 L V_0}{g\, T} = \frac{2 \times 1200 \times 2.5}{9.81 \times 3.0} = \frac{6000}{29.43} = 203.9\ \text{m} ΔH204 m of water\boxed{\Delta H \approx 204\ \text{m of water}}

Corresponding pressure rise:

Δp=ρgΔH=1000×9.81×203.9=2.00×106 Pa=2.00 MPa\Delta p = \rho g \Delta H = 1000 \times 9.81 \times 203.9 = 2.00\times10^6\ \text{Pa} = \boxed{2.00\ \text{MPa}}

The gradual-closure pressure (2.00 MPa) is lower than the instantaneous value (2.75 MPa), confirming that slower valve operation reduces water-hammer over-pressure.

water-hammersurgepressure-rise
8short6 marks

A rectangular channel 5.0 m5.0\ \text{m} wide carries water at a depth of 1.2 m1.2\ \text{m} with a velocity of 1.5 m/s1.5\ \text{m/s}. Due to sudden partial closure of a downstream gate, a positive surge (moving hydraulic jump) travels upstream, raising the depth to 1.8 m1.8\ \text{m}.

(a) Briefly classify surges in open channels.

(b) Determine the absolute velocity (celerity) of the surge and the new flow velocity behind the surge.

(a) Classification of surges

A surge is a moving wave front (rapidly varied unsteady flow) caused by a sudden change in discharge or depth. Types:

  • Positive surge: depth increases as the wave passes (a moving hydraulic jump). It can move upstream (e.g. gate closing) or downstream (e.g. gate opening downstream).
  • Negative surge: depth decreases as the wave passes; it is a flattening wave (e.g. gate opening upstream, dam-break receding wave).

Each can travel either upstream or downstream, giving four sub-cases.

(b) Upstream positive surge

Let the surge move upstream with absolute speed VwV_w (positive in the downstream direction, so an upstream-moving surge has VwV_w negative; we solve for its magnitude). Use a reference frame moving with the surge to make the flow steady — then continuity and momentum of a stationary hydraulic jump apply.

Let downstream-positive be the flow direction. Initial: y1=1.2 my_1 = 1.2\ \text{m}, V1=1.5 m/sV_1 = 1.5\ \text{m/s} (toward the gate). After surge: y2=1.8 my_2 = 1.8\ \text{m}, velocity V2V_2 (unknown). Surge absolute velocity =Vw= V_w (upstream, so we take its magnitude CC upstream).

Continuity in the moving frame (relative velocities V1+CV_1 + C and V2+CV_2 + C, since the surge moves upstream against the flow):

y1(V1+C)=y2(V2+C)(1)y_1 (V_1 + C) = y_2 (V_2 + C) \quad (1)

Momentum (moving frame, treat as stationary jump):

(V1+C)2=gy22y1(y1+y2)(2)(V_1 + C)^2 = \frac{g\, y_2}{2 y_1}(y_1 + y_2) \quad (2)

From (2):

(V1+C)2=9.81×1.82×1.2(1.2+1.8)=17.6582.4×3.0=7.3575×3.0=22.073(V_1 + C)^2 = \frac{9.81 \times 1.8}{2 \times 1.2}(1.2 + 1.8) = \frac{17.658}{2.4}\times 3.0 = 7.3575 \times 3.0 = 22.073 V1+C=22.073=4.698 m/sV_1 + C = \sqrt{22.073} = 4.698\ \text{m/s} C=4.6981.5=3.198 m/sC = 4.698 - 1.5 = 3.198\ \text{m/s} Surge celerity (absolute, upstream)3.20 m/s\boxed{\text{Surge celerity (absolute, upstream)} \approx 3.20\ \text{m/s}}

New velocity behind the surge from (1):

V2+C=y1(V1+C)y2=1.2×4.6981.8=5.63761.8=3.132 m/sV_2 + C = \frac{y_1 (V_1 + C)}{y_2} = \frac{1.2 \times 4.698}{1.8} = \frac{5.6376}{1.8} = 3.132\ \text{m/s} V2=3.1323.198=0.066 m/sV_2 = 3.132 - 3.198 = -0.066\ \text{m/s} V20.07 m/s (nearly stagnant; slight reverse)\boxed{V_2 \approx -0.07\ \text{m/s (nearly stagnant; slight reverse)}}

The flow behind the upstream-moving surge is almost brought to rest, consistent with a near-complete downstream gate closure. The surge advances upstream at about 3.2 m/s3.2\ \text{m/s}.

surgesmoving-hydraulic-jumpunsteady-flow
9short6 marks

(a) Explain, with a sketch, the flow over a broad-crested weir and the concept of critical flow at the crest.

(b) A broad-crested weir 4.0 m4.0\ \text{m} long (across the flow) has a measured upstream head over the crest of 0.50 m0.50\ \text{m}. Taking the coefficient of discharge Cd=0.85C_d = 0.85 and neglecting the approach velocity, compute the discharge using the standard broad-crested-weir formula.

(a) Flow over a broad-crested weir

A broad-crested weir has a crest long enough (in the flow direction) that the streamlines become horizontal and parallel over the crest, so the pressure is hydrostatic and the flow passes through critical depth somewhere on the crest. Because critical flow gives maximum discharge for the available specific energy, the discharge is controlled by the head over the crest.

   ~~~~~~~~  H (upstream head)
        ____________  <- water surface drops over crest
  ----->|  yc (critical depth on crest)
  flow  |============| crest (broad)
        |   weir     |
  ______|____________|_____

At the crest the specific energy H\approx H and depth =yc=23H= y_c = \tfrac{2}{3}H (rectangular). Substituting Vc=gycV_c = \sqrt{g y_c} into q=Vcycq = V_c y_c gives the weir formula.

(b) Discharge computation

Standard broad-crested-weir formula:

Q=Cd2323g  LH3/2Q = C_d\,\frac{2}{3}\sqrt{\frac{2}{3}g}\; L\, H^{3/2}

Evaluate the constant:

2323×9.81=236.54=23(2.5573)=1.7049\frac{2}{3}\sqrt{\frac{2}{3}\times 9.81} = \frac{2}{3}\sqrt{6.54} = \frac{2}{3}(2.5573) = 1.7049

With Cd=0.85C_d = 0.85, L=4.0 mL = 4.0\ \text{m}, H=0.50 mH = 0.50\ \text{m}:

H3/2=0.501.5=0.35355H^{3/2} = 0.50^{1.5} = 0.35355 Q=0.85×1.7049×4.0×0.35355Q = 0.85 \times 1.7049 \times 4.0 \times 0.35355 Q=0.85×1.7049=1.4492;×4.0=5.7967;×0.35355=2.0494Q = 0.85 \times 1.7049 = 1.4492;\quad \times 4.0 = 5.7967;\quad \times 0.35355 = 2.0494 Q2.05 m3/s\boxed{Q \approx 2.05\ \text{m}^3/\text{s}}

Check via critical-flow basis: yc=23(0.50)=0.3333 my_c = \tfrac{2}{3}(0.50) = 0.3333\ \text{m}, Vc=9.81×0.3333=1.808 m/sV_c = \sqrt{9.81 \times 0.3333} = 1.808\ \text{m/s}, ideal q=1.808×0.3333=0.6027 m2/sq = 1.808 \times 0.3333 = 0.6027\ \text{m}^2/\text{s}, ideal Q=0.6027×4.0=2.411Q = 0.6027 \times 4.0 = 2.411; with Cd=0.85C_d=0.85, Q=2.05 m3/sQ = 2.05\ \text{m}^3/\text{s} — consistent.

hydraulic-structuresbroad-crested-weirdischarge-measurement
10short4 marks

Water flows in a rectangular channel 3.0 m3.0\ \text{m} wide at a depth of 1.5 m1.5\ \text{m} and discharge 9 m3/s9\ \text{m}^3/\text{s}. A smooth upward step (hump) is built in the bed. Determine the maximum height of the hump that can be placed without causing the upstream depth to change (i.e. the height that just brings the flow to critical over the hump).

Solution

Unit discharge: q=9/3=3.0 m2/sq = 9/3 = 3.0\ \text{m}^2/\text{s}.

Upstream conditions:

V1=qy1=3.01.5=2.0 m/sV_1 = \frac{q}{y_1} = \frac{3.0}{1.5} = 2.0\ \text{m/s} E1=y1+V122g=1.5+2.0219.62=1.5+0.2039=1.7039 mE_1 = y_1 + \frac{V_1^2}{2g} = 1.5 + \frac{2.0^2}{19.62} = 1.5 + 0.2039 = 1.7039\ \text{m}

Froude number: Fr1=2.0/9.81×1.5=2.0/3.836=0.521Fr_1 = 2.0/\sqrt{9.81 \times 1.5} = 2.0/3.836 = 0.521 (subcritical).

Critical conditions over the hump (the flow reaches minimum specific energy):

yc=(q2g)1/3=(3.029.81)1/3=(0.9174)1/3=0.9717 my_c = \left(\frac{q^2}{g}\right)^{1/3} = \left(\frac{3.0^2}{9.81}\right)^{1/3} = (0.9174)^{1/3} = 0.9717\ \text{m} Ec=32yc=32(0.9717)=1.4576 mE_c = \frac{3}{2} y_c = \frac{3}{2}(0.9717) = 1.4576\ \text{m}

Maximum hump height (critical flow corresponds to minimum specific energy on the hump; energy on the hump E2=E1ΔzE_2 = E_1 - \Delta z, set E2=EcE_2 = E_c):

Δzmax=E1Ec=1.70391.4576=0.2463 m\Delta z_{max} = E_1 - E_c = 1.7039 - 1.4576 = 0.2463\ \text{m} Δzmax0.246 m\boxed{\Delta z_{max} \approx 0.246\ \text{m}}

For any hump higher than this, the upstream depth must rise (choking occurs); a hump of 0.246 m0.246\ \text{m} just brings the flow to critical over the crest without disturbing the upstream depth.

non-uniform-flowcritical-flowchannel-transition
11short4 marks

(a) Define conveyance KK and the section factor for uniform flow.

(b) A circular concrete sewer of diameter 1.0 m1.0\ \text{m} (n=0.013n = 0.013) flows half full on a slope S0=0.001S_0 = 0.001. Compute the discharge.

(a) Conveyance and section factor

Manning's equation can be written Q=KS1/2Q = K\,S^{1/2}, where the conveyance:

K=1nAR2/3K = \frac{1}{n} A R^{2/3}

represents the discharge-carrying capacity of a section per unit S\sqrt{S}; it depends only on geometry and roughness. The section factor for uniform flow is AR2/3A R^{2/3} (=nK= nK), while the section factor for critical flow is Z=AA/TZ = A\sqrt{A/T}.

(b) Half-full circular sewer

Diameter D=1.0 mD = 1.0\ \text{m}, radius r=0.5 mr = 0.5\ \text{m}. At half full the water surface is along the diameter:

Area (half circle):

A=12πD24=πD28=π(1.0)28=0.3927 m2A = \frac{1}{2}\,\frac{\pi D^2}{4} = \frac{\pi D^2}{8} = \frac{\pi (1.0)^2}{8} = 0.3927\ \text{m}^2

Wetted perimeter (half the circumference):

P=12πD=π(1.0)2=1.5708 mP = \frac{1}{2}\,\pi D = \frac{\pi (1.0)}{2} = 1.5708\ \text{m}

Hydraulic radius:

R=AP=0.39271.5708=0.2500 m(=D/4, as expected for half-full)R = \frac{A}{P} = \frac{0.3927}{1.5708} = 0.2500\ \text{m}\quad (= D/4,\ \text{as expected for half-full})

Discharge (Manning):

Q=1nAR2/3S1/2Q = \frac{1}{n} A R^{2/3} S^{1/2} R2/3=0.252/3=0.3969R^{2/3} = 0.25^{2/3} = 0.3969 S1/2=0.001=0.031623S^{1/2} = \sqrt{0.001} = 0.031623 Q=10.013×0.3927×0.3969×0.031623Q = \frac{1}{0.013} \times 0.3927 \times 0.3969 \times 0.031623 Q=76.923×0.3927×0.3969×0.031623Q = 76.923 \times 0.3927 \times 0.3969 \times 0.031623 =76.923×0.3927=30.21; ×0.3969=11.99; ×0.031623=0.3792= 76.923 \times 0.3927 = 30.21;\ \times 0.3969 = 11.99;\ \times 0.031623 = 0.3792 Q0.38 m3/s\boxed{Q \approx 0.38\ \text{m}^3/\text{s}}
uniform-flowconveyancesection-factor

Frequently asked questions

Where can I find the BE Civil Engineering (IOE, TU) Hydraulics (IOE, CE 604) question paper 2076?
The full BE Civil Engineering (IOE, TU) Hydraulics (IOE, CE 604) 2076 (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) 2076 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) 2076 paper?
The BE Civil Engineering (IOE, TU) Hydraulics (IOE, CE 604) 2076 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.