Chapter 6 Gradually-Varied Flow in Open Channels

6.1. Introduction A steady non-uniform flow in a prismatic channel with gradual changes in its watersurface elevation is named as gradually-varied flow (GVF). The backwater produced by a dam or weir across a river and drawdown produced at a sudden drop in a channel are few typical examples of GVF. In a GVF, the velocity varies along the channel and consequently the bed slope, water surface slope, and energy line slope will all differ from each other. Regions of high curvature are excluded in the analysis of this flow. The two basic assumptions involved in the analysis of GVF are: 1. The pressure distribution at any section is assumed to be hydrostatic. This follows from the definition of the flow to have a gradually varied water surface. As gradual changes in the surface curvature give rise to negligible normal accelerations, the departure from the hydrostatic pressure distribution is negligible. 2. The resistance to flow at any depth is assumed to be given by the corresponding uniform flow equation, such as the Manning equation, with the condition that the slope term to be used in the equation is the energy line slope, Se and not the bed slope, S0. Thus, if in a GVF the depth of flow at any section is y, the energy line slope Se is given by, Se = n 2V 2 R4 3 (6.1)

where R = hydraulic radius of the section at depth y. 6.2. Basic Differential Equation for the Gradually-Varied Flow Water Surface Profile Since, S0 = Channel bed slope for uniform flow depth y0, dy = Water depth variation for the dx canal reach, d (V2/2g) = Velocity head variation for the dx reach. Writing the energy equation between the cross-sections 1 and 2 (Fig. 6.1.),

Prof. Dr. Atl BULU

V 2 V2 V2 = ( y + dy ) + + d 2 g + Se dx 2g 2g 2 V S 0 dx Se dx = dy + d 2g S 0 dx + y + S0 Se = S0 Se = dy d V 2 + dx dx 2 g dy d V 2 1 + dx dy 2 g

dy = dx
Se

S0 Se d V 2 1+ dy 2 g

(6.2)

Se d E.L y V y+d S0 d V+d d S0 V2/2g+d(V2/2g

Figure 6. 1.

V2 Q2 = dA = T ( y )dy 2 g 2 gA( y )2 dA Q2 2 A Q2 d Q 2T ( y ) dy = = 4 3 dy 2 gA( y )2 gA( y ) 2 gA( y ) V2 = = Fr2 A( y ) g T (y) Substituting this to Equ. (6.2),

dy S0 S e = dx 1 Fr2
2

(6.3)

Prof. Dr. Atl BULU

Equ. (6.3) is the general differential equation of the water surface profile for the gradually varied flows. dy/dx gives the variation of water depth along the channel in the flow direction.
6.3. Classification of Flow Surface Profiles

For a given channel with a known Q = Discharge, n = Manning coefficient, and S0 = Channel bed slope, yc = critical water depth and y0 = Uniform flow depth can be computed. There are three possible relations between y0 and yc as 1) y0 > yc , 2) y0 < yc , 3) y0 = yc . For horizontal (S0 = 0), and adverse slope ( S0 < 0) channels,
1 1 Q = A R 2 3 S0 2 n Horizontal channel, S0 = 0 Q = 0, Adverse channel , S0 < 0 , Q cannot be computed, For horizontal and adverse slope channels, uniform flow depth y0 does not exist. Based on the information given above, the channels are classified into five categories as indicated in Table (6.1).
Table 6.1. Classification of channels Number Channel category 1 Mild slope 2 Steep slope Symbol Characteristic condition M y0 > yc S yc > y0 Remark

3 4 5

Critical slope Horizontal bed Adverse slope

C H A

yc = y0 S0 = 0 S0 < 0

Subcritical flow at normal depth Supercritical flow at normal depth Critical flow at normal depth Cannot sustain uniform flow Cannot sustain uniform flow

For each of the five categories of channels, lines representing the critical depth (yc ) and normal depth (y0 ) (if it exists) can be drawn in the longitudinal section. These would divide the whole flow space into three regions as: Region 1: Space above the topmost line, Region 2: Space between top line and the next lower line, Region 3: Space between the second line and the bed.

Prof. Dr. Atl BULU

Figure (6.2) shows these regions in the various categories of channels.

Figure 6.2. Regions of flow profiles

Depending upon the channel category and region of flow, the water surface profiles will have characteristics shapes. Whether a given GVF profile will have an increasing or decreasing water depth in the direction of flow will depend upon the term dy/dx in Equ. (6.3) being positive or negative. dy S0 S e = dx 1 Fr2

(6.3)

Prof. Dr. Atl BULU

For a given Q, n, and S0 at a channel, y0 = Uniform flow depth, yc = Critical flow depth, y = Non-uniform flow depth. The depth y is measured vertically from the channel bottom, the slope of the water surface dy/dx is relative to this channel bottom. Fig. (6.3) is basic to the prediction of surface profiles from analysis of Equ. (6.3).

Figure 6.3

To assist in the determination of flow profiles in various regions, the behavior of dy/dx at certain key depths is noted by studying Equ. (6.3) as follows: y > yc Fr < 1 y = yc Fr = 1 y < yc Fr > 1 And also, y > y0 S e < S 0 y = y0 S e = S 0 y < y0 S e > S 0 1. As y y0 , V V0 , S e = S0 0 dy S 0 S 0 = = =0 2 y y0 dx 1 Fr cons lim The water surface approaches the normal depth asymptotically.

Prof. Dr. Atl BULU

2. As y yc , Fr2 = 1 , 1 Fr2 = 0 , lim dy S 0 S e S 0 S e = = = 0 dx 1 Fr2

y yc

The water surface meets the critical depth line vertically. 3. As y , V = 0 Fr = 0 Se 0

The water surface meets a very large depth as a horizontal asymptote. Based on this information, the various possible gradually varied flow profiles are grouped into twelve types (Table 6.2).
Table 6.2. Gradually Varied Flow profiles Channel Region Condition Type

Mild slope Steep slope Critical slope Horizontal bed Adverse slope

1 2 3 1 2 3 1 3 2 3 2 3

y > y0 > yc y0 > y > yc y0 > yc > y y> yc > y0 yc > y > y0 yc > y0 > y y > y0 = y2 y < y0 = yc y > yc y < yc y > yc y < yc

M1 M2 M3 S1 S2 S3 C1 C3 H2 H3 A2 A3

Prof. Dr. Atl BULU

6.4. Water Surface Profiles 6.4.1. M Curves

Figure 6.4

General shapes of M curves are given in Fig. (6.4). Asymptotic behaviors of each curve will be examined mathematically.
a) M1 Curve

Water surface will be in Region 1 for a mild slope channel and the flow is obviously subcritical. Se < S0 Mild slope channel

y0 > yc Subcritical flow dy S 0 S e = dx 1 Fr2 Fr < 1 Subcritical flow (1 Fr2) > 0 y > y0 Se < S0 dy + = > 0 (Water depth will increase in the flow direction) dx + Asymptotic behavior of the water surface is; Water depth can be between ( > y > y0) for Region 1. The asymptotic behaviors of the water surface for the limit values ( , y0) are;

Prof. Dr. Atl BULU

a)

y , V 0, Fr 0, (1 Fr2) = 1 y , V 0, Se 0 lim dy S 0 0 = = S0 dx 1

The water surface meets a very large depth as a horizontal asymptote. b) y y0 , V V0 , Se S0 lim dy S 0 S 0 = =0 dx 1 Fr2

y y0

The water surface approaches the normal depth asymptotically. The most common of all GVF profiles is the M1 type, which is a subcritical flow condition. Obstructions to flow, such as weirs, dams, control structures and natural features, such as bends, produce M1 backwater curves (Fig. 6.5). These extend to several kilometers upstream before merging with the normal depth.

Figure 6.5. M1 Profile b) M2 Curve

Water surface will be in Region 2 for a mild slope channel and the flow is obviously subcritical. (Fig. 6.4). y0 > y > yc y < y0 V > V0 Se > S0 (S0 Se) < 0 y > yc Fr < 1 (1 Fr2) > 0 dy S 0 S e + = = = (Water depth decrease in the flow direction) dx 1 Fr2

Prof. Dr. Atl BULU

Asymptotic behavior of the water surface is; y0 > y > yc y y0, Se S0 , (S0 Se) = 0 dy S 0 S 0 = =0 y y0 dx 1 Fr2 lim The water surface approaches the normal depth asymptotically. y yc , Fr = 1 , (1 Fr2 ) = 0

y yc

lim

dy S 0 S e = = dx 0

The water surface meets the critical depth line vertically. The M2 profiles occur at a sudden drop of the channel, at constriction type of transitions and at the canal outlet into pools (Fig. 6.6).

Figure 6.6 c) M3 Curve

Water surface will be in Region 3 for a mild slope channel and the flow is obviously subcritical. (Fig. 6.4). y0 > yc >0 y < y0 V > V0 Se > S0 (S0 Se) < 0 y < yc Supercritical flow Fr > 1 (1 Fr2) < 0 dy S0 S e = = = + (Water depth will increase in the flow direction) dx 1 Fr2 Asymptotic behavior of the water surface is;

Prof. Dr. Atl BULU

y yc Fr = 1 (1 Fr2) = 0 lim dy S0 Se = = dx 0

y yc

The water surface meets the critical depth line vertically. y 0, Se , (S0 Se) = y 0, Fr =
V = gy

lim
y 0

dy = (Unknown) dx
3

y The angle of the water surface with the channel bed may be taken as S0 0 . y c

Figure 6.7. M3 Profile

Where a supercritical stream enters a mild slope channel, M3 type of profile occurs. The flow leading from a spillway or a sluice gate to a mild slope forms a typical example (Fig. 6.7). The beginning of the M3 curve is usually followed by a small stretch of rapidly varied flow and the downstream is generally terminated by a hydraulic jump. Compared to M1 and M2 profiles, M3 curves are of relatively short length.

10

Prof. Dr. Atl BULU

Example 6.1 : A rectangular channel with a bottom width of 4.0 m and a bottom slope of 0.0008 has a discharge of 1.50 m3/sec. In a gradually varied flow in this channel, the depth at a certain location is found to be 0.30 m. assuming n = 0.016, determine the type of GVF profile. Solution:

a) To find the normal depth y0,

1 A 1 Q = VA = A S 0 2 nP
53

1.50 =

(4 y0 ) 0.00080.5 1 0.016 (4 + 2 y0 )2 3
0.0842 =
5 y0 3 (4 + 2 y0 )2 3

By trial and error, y0 = 0.43 m. b) Critical depth yc, q= Q 1.50 = = 0.375 m 2 sec B 4.0
13

q2 yc = g c) Type of profile,

0.3752 = 9.81

13

= 0.24m

y0 = 0.43 m > yc = 0.24 m (Mild slope channel, M profile) y = 0.30 m y0 > y > yc (Region 2) Water surface profile is of the M2 type.

11

Prof. Dr. Atl BULU

6.4.2. S Curves

Figure 6.8

General shapes of S curves are given in Fig. (6.8). Asymptotic behaviors of each curve will be examined mathematically.
a) S1 Curve

Water surface will be in Region 1 for a steep slope channel and the flow is obviously supercritical. y > yc > y0 y > y0 V < V0 Se < S0 (S0 Se) > 0 y > yc Fr < 1 (1 Fr2) > 0 dy S0 S e + = = =+ + dx 1 Fr2 (Water depth will increase in the flow direction) Asymptotic behavior of the water surface is; > y > yc y , V 0 , Fr = 0 , (1- Fr2) = 1 y , V 0 , Se = 0 lim
y

dy S 0 Se S0 0 = = = S0 dx 1 Fr2 1 0

12

Prof. Dr. Atl BULU

The water surface meets a very large depth as a horizontal asymptote. y yc , Fr = 1 , (1 Fr2) = 0 dy S0 Se S0 S e = = = y yc dx 1 Fr2 0 lim The water surface meets the critical depth line vertically.

Figure 6.9. S1 Profile

The S1 profile is produced when the flow from a steep channel is terminated by a deep pool created by an obstruction, such as a weir or dam (Fig. 6.9). At the beginning of the curve, the flow changes from the normal depth (supercritical flow) to subcritical flow through a hydraulic jump. The profiles extend downstream with a positive water slope to reach a horizontal asymptote at the pool elevation.
b) S2 Curve

Water surface will be in Region 2 for a steep slope channel and the flow is supercritical. y0 < y < yc y > y0 Se < S0 (S0 Se) > 0 y < yc Supercritical flow Fr > 1 (1 Fr2) < 0 dy S0 S e + = = = dx 1 Fr2 (Water depth will decrease in the flow direction)

13

Prof. Dr. Atl BULU

Asymptotic behavior of the water surface is; y0 < y < yc y yc , Fr 1 , (1 Fr2) = 0 dy S0 Se S0 S e = = = y yc dx 1 Fr2 0 lim The water surface meets the critical depth line vertically. y y0 , V V0 , Se S0 , (S0 Se) = 0 lim 0 dy = =0 dx 1 Fr2

y y0

The water surface approaches the normal depth asymptotically.

Figure 6.10. S2 Profile

Profiles of the S2 type occur at the entrance region of a steep channel leading from a reservoir and at a brake of grade from mild slopes to steep slope (Fig. 6.10). Generally S2 profiles are short of length.
c) S3 Curve

Water surface will be in Region 3 for a steep slope channel and the flow is supercritical. 0 < y < y0 y < y0 Se > S0 (S0 Se) < 0 y < yc Supercritical flow Fr > 1 (1 Fr2) < 0 dy S0 S e = = =+ dx 1 Fr2 (Water depth will increase in the flow direction)

14

Prof. Dr. Atl BULU

Asymptotic behavior of the water surface is; 0 < y < y0 y 0, Se , (S0 Se) = - y 0 , Fr = lim
y 0

V = gy

dy = dx

(Unknown)
3

y The angle of the water surface with the channel bed may be taken as S0 0 . y c y y0 , Se S0 , (S0 Se) = 0 0 dy = =0 y y0 dx 1 Fr2 lim The water surface approaches the normal depth y0 asymptotically.

Figure 6.11. S3 Profile

Figure 6.12. S3 Profile

Free flow from a sluice gate with a steep slope on its downstream is of the S3 type (Fig. 6.11). The S3 curve also results when a flow exists from a steeper slope to a less steep slope (Fig. 6.12).

15

Prof. Dr. Atl BULU

6.4.3. C Curves

Figure 6.13

General shapes of C curves are given in Fig. (6.13). Asymptotic behaviors of each curve will be examined mathematically. Since the flow is at critical stage, y0 = yc, there is no Region 2.
a) C1 Curve

Water surface will be in Region 1 for a critical slope channel. y0 = yc < y < y > yc Se < S0 (S0 Se) > 0 y > yc Subcritical flow Fr < 1 (1 Fr2) > 0 dy S 0 S e + = = =+ + dx 1 Fr2 (Water depth will increase in the flow direction) Asymptotic behavior of the water surface is; y , Se 0 , ( S0 Se) = S0 y , V = 0 Fr = 0 (1 Fr2) = 1 lim dy S 0 0 = = S0 y dx 1 0

The water surface meets a very large depth as a horizontal asymptote. y yc , Fr = 1 (1 Fr2) = 0 y yc = y0 , Se = S0 = Sc (S0 Se ) =0 dy 0 = y yc = y0 dx 0 lim (Unknown)

C1 and C3 profiles are very rare and highly unstable 16

6.4.4.

H Curves

Figure 6.14

General shapes of H curves are given in Fig. (6.14). For horizontal slope channels, uniform flow depth y0 does not exist. Critical water depth can be computed for a given discharge Q and therefore critical water depth line can be drawn. Since there is no uniform water depth y0, Region 1 does not exist.
a) H2 Curve

Water surface will be in Region 2 for a horizontal slope channel. > y > yc y > yc Se < S0 = 0 (S0 Se) < 0 y > yc subcritical flow Fr < 1 (1 Fr2) > 0 dy = = (Water depth will decrease in the flow direction) dx + Asymptotic behavior of the water surface is; y , Se 0 , ( S0 Se) = S0 y , V = 0 Fr = 0 (1 Fr2) = 1 lim dy S 0 0 = = S0 y dx 1 0

The water surface meets a very large depth as a horizontal asymptote. y yc , Fr = 1 (1 Fr2) = 0 y yc , Se = Sc S0 = 0, (S0 Se ) = -Se dy S 0 S e 0 S e = = = y yc dx 0 1 Fr2 The water surface meets the critical depth line vertically. lim

17

Prof. Dr. Atl BULU

b) H3 Curve

0 > y > yc y < yc Se > S0 = 0 (S0 Se) = -Se

y yc , supercritical flow , Fr > 1 (1 Fr2 ) < 0 dy S 0 S e = = =+ dx 1 Fr2 (Water depth will increase in the flow direction) Asymptotic behavior of the water surface is; y 0, Se , (S0 Se) = - y 0 , Fr = lim dy = y 0 dx
V = gy

(Unknown)

y yc Fr = 1 (1 Fr2) = 0 y yc , Se = Sc S0 = 0, (S0 Se ) = -Se dy S 0 S e 0 S e = = = y yc dx 0 1 Fr2 lim The water surface meets the critical depth line vertically.

Figure 6.15

A horizontal channel can be considered as the lower limit reached by a mild slope as its bed slope becomes flatter. The H2 and H3 profiles are similar to M2 and M3 profiles respectively (Fig. 6.15). However, the H2 curve has a horizontal asymptote. 18
Prof. Dr. Atl BULU

6.4.5.

A Curves

Figure 6.16

General shapes of A curves are given in Fig. (6.16). For adverse slope channels, uniform flow depth y0 does not exist. Critical water depth can be computed for a given discharge Q and therefore critical water depth line can be drawn. Since there is no uniform water depth y0, Region 1 does not exist as well as in A curves. A2 and A3 curves are similar to H2 and H3 curves respectively.

Figure 6.17

Adverse slopes are rather rare (Fig. 6.17). These profiles are of very short length.
6.5. Control Sections

A control section is defined as a section in which a fixed relationship exists between the discharge and depth of flow. Weirs, spillways, sluice gates are some typical examples of structures which give rise to control sections. The critical depth is also a control point. However, it is effective in a flow profile which changes from subcritical to supercritical flow. In the reverse case of transition from supercritical flow to subcritical flow, a hydraulic jump is usually formed bypassing the critical depth as a control point. Any GVF profile will have at least one control section.

19

Prof. Dr. Atl BULU

(d)
Figure 6.18

In the synthesis of GVF profiles occurring in a serially connected channel elements, the control sections provide a key to the identification of proper profile shapes. A few typical control sections are shown in Fig. (6.18, a-e). It may be noted that subcritical flows have controls in the downstream end, while supercritical flows are governed by control sections existing at the upstream end of the channel section. In Figs. (6.18 a and b) for the M1 profile, the control section (indicated by a dark dot in the figures) is just at the upstream of the spillway and sluice gate respectively. In Figs. (6.18 b and d) for M3 and S3 profiles respectively, the control point is at the vena contracta of the sluice gate flow. In subcritical flow reservoir offtakes (Fig. 6.18 c), even though the discharge is governed

20

Prof. Dr. Atl BULU

by the reservoir elevation, the channel entry section is not strictly a control section. The water surface elevation in the channel, will be lower than the reservoir elevation by a headloss amount equivalent to (1 + )V2/2g where is the entrance loss coefficient. The true control section will be at a downstream location in the channel. For the situation shown in Fig. (6.18 c) the critical depth at the free overflow at the channel end acts as the downstream control. For a sudden drop (free overflow) due to curvature of the streamlines the critical depth usually occurs at distance of about 4yc upstream of the drop. This distance, being small compared to GVF lengths, is neglected and it is usual to perform calculations by assuming yc to occur at the drop.

(e) For a supercritical canal intake (Fig. 6.18 e), the reservoir water surface falls to the critical depth at the head of the canal and then onwards the water surface follows the S2 curve. The critical depth occurring at the upstream end of the channel is the control for this flow.

Figure 6.19

21

Prof. Dr. Atl BULU

6.5.1. Feeding a Pool for Subcritical Flow

A mild slope channel discharging into a pool of variable surface elevation is indicated in Fig. (6.19). Four cases are shown. 1) The pool elevation is higher than the elevation of the normal (uniform) water depth line at B. This gives rise to a drowning of the channel end. A profile of the M1 type is produced with the pool level at B as control. 2) The pool elevation is lower than the elevation of the normal depth line but higher than the critical depth line at B. The pool elevation acts as a control for the M2 curve. 3) The pool elevation has dropped down to that of the critical depth line at B and the control is still at the pool elevation. 4) The pool elevation has dropped lower than the elevation of the critical depth line at B. The water surface cannot pass through a critical depth at any location other than B and hence a sudden drop in the water surface at B is observed. The critical depth at B is the control of this flow.
6.5.2. Feeding a Pool for Supercritical Flow

When the slope of the channel is steep, that is greater than the critical slope, the flow in the channel becomes supercritical. In practical applications, steep channels are usually short, such as raft and log chutes that are used as spillways.

Figure 6.20

As the control section in a channel of supercritical flow is at the upstream end, the delivery of the channel is fully governed by the critical discharge at section 1.

22

Prof. Dr. Atl BULU

1) When the tailwater level B is less than the outlet depth (normal depth) at section 2, the flow in the canal is unaffected by the tailwater. The flow profile passes through the critical water depth near C and approaches the normal depth by means of a smooth drawdown curve of the S2 type. 2) When the tailwater level B is greater than the outlet depth, the tailwater will raise the water level in the upstream portion of the canal to form an S1 profile between j and b`, producing a hydraulic jump at the end j of the profile. However, the flow upstream from the jump will not be affected by the tailwater. Such a hydraulic jump is objectionable and dangerous, particularly when the canal is raft chute or some other structure intended to transport a floating raft from the upstream reservoir to a downstream pool. A neutralizing reach may be suggested for the solution of this problem. (Chow, 1959). (Fig. 6.21)

Figure 6.21. Elimination of hydraulic jump by neutralizing reach.

In this reach the bottom slope of the channel is made equal to the critical slope. According to a corresponding case of C1 profile in Fig. (6.21), the tailwater levels will be approximately horizontal lines which intersect the surface of flow in the canal without causing any disturbance. At the point of intersection, theoretically, there is a jump of zero height. 3) As the tailwater rises further, the jump will move upstream, maintaining its height and form in the uniform flow zone nb, until it reaches point n. The height of the jump becomes zero when it reaches the critical depth at c. In the meantime, the flow profile reaches its theoretical limit cb`` of the S1 profile. Beyond this limit the incoming flow will be directly affected by the tailwater. In practical applications, the horizontal line cb``` may be taken as the practical limit of the tailwater stage.

23

Prof. Dr. Atl BULU

6.6.

Analysis of Flow Profile

A channel carrying a gradually varied flow can in general contain different prismoidal channel cross-sections of varying hydraulic properties. There can be a number of control sections at various locations. To determine the resulting water surface profile in a given case, one should be in a position to analyze the effects of various channel sections and controls connected in series. Simple cases are illustrated to provide information and experience to handle more complex cases.
6.6.1. Break in Grades

Simple situations of a series combination of two channel sections with differing bed slopes are considered. In Fig. (6.22.a), a break in grade from a mild channel to a milder channel is shown. It is necessary to first draw the critical-depth line (CDL) and the normal-depth line (NDL) for both slopes. Since yc does not depend upon the slope for a taken Q = discharge, the CDL is at a constant height above the channel bed in both slopes. The normal depth y01 for the mild slope is lower than that of the of the milder slope (y02). In this case, y02 acts as a control, similar to the weir or spillway case and an M1 backwater curve is produced in the mild slope channel. Various combinations of slopes and the resulting GVF profiles are presented in Fig. ( 6.22, a-h). It may be noted that in some situations there can be more than one possible profiles. For example, in Fig. (6.22 e), a jump and S1 profile or an M3 profile and a jump possible. The particular curve in this case depends on the channel and its flow properties. In the examples indicated in Fig. (6.22), the section where the grade changes acts a control section and this can be classified as a natural control. It should be noted that even though the bed slope is considered as the only variable in the above examples, the same type of analysis would hold good for channel sections in which there is a marked change in the roughness characteristics with or without change in the bed slope. A long reach of unlined canal followed by a line reach serves as a typical example for the same. A change in the channel geometry (the bed width or side slope) beyond a section while retaining the prismoidal nature in each reach also leads to a natural control section.

24

Prof. Dr. Atl BULU

25

Prof. Dr. Atl BULU

26

Prof. Dr. Atl BULU

Figure 6.22

6.6.2. Serial Combination of Channel Sections

To analyze a general problem of many channel sections and controls, the following steps are to be applied. 1. Draw the longitudinal section of the system. 2. Calculate the critical depth and normal depths of various reaches and draw the CDL and NDL in all reaches. 3. Mark all the controls, both the imposed as well as natural controls. 4. Identify the possible profiles.
Example 6. 2: Identify and sketch the GVF profiles in three mild slopes which could be described as mild, steeper mild and milder. The three slopes are in series. The last slope has a sluice gate in the middle of the reach and the downstream end of the channel has a free overfall. Solution: The longitudinal profile of the channel, critical depth line and normal depth lines for the various reaches are shown in the Fig. (6.23). The free overfall at E is obviously a control. The vena contracta downstream of the sluice gate at D is another control. Since for subcritical flow the control is at the downstream end of the channel, the higher of the two normal depths at C acts as a control for the reach CB, giving rise to an M1 profile over CB. At B, the normal depth of the channel CB acts as a control giving rise to an M2 profile over AB. The controls are marked distinctly in Fig. (6.23). With these controls the possible flow profiles are: an M2 profile on channel AB, M1 profile on channel BC, M3 profile and M2 profile through a jump on the stretch DE.

27

Prof. Dr. Atl BULU

Figure 6.23 Example 6.3: A trapezoidal channel has three reaches A, B, and C connected in series with the following physical characteristics. Reach Bed width (B) 4.0 m 4.0 m 4.0 m Side slope (m) 1.0 1.0 1.0 Bed slope S0 0.0004 0.009 0.004 n

1 2 3

0.015 0.012 0.015

For a discharge Q = 22.5 m3/sec through this channel, sketch the resulting water surface profiles. The length of the reaches can be assumed to be sufficiently long for the GVF profiles to develop fully.
Solution: The normal depths and critical water depths in the various reaches are calculated as:

y 1 m

B=4m

A = (B + my ) y P = B + 2 y 1 + m2 T = B + 2my

28

Prof. Dr. Atl BULU

A = (4 + y ) y P = 4 + 2 y 1 + 12 = 4 + 23 2 y T = 4 + 2y

A A 0 Q = AV = S 0 .5 n P

23

[(4 + y0 ) y0 ]5 3 22.5 =
n

(4 + 2

0 S 0 .5 32

y0 )

23

Uniform flow depths for the given data for every reach are calculated by trial and error method; Reach A: Reach B: Reach C: S0A = 0.0004 , nA=0.015 y0A = 2.26 m S0B = 0.009 , nB =0.012 y0B = 0.81 m S0C = 0.004 , nC = 0.015 y0C =1.17 m

Since the channel is prismoidal (the geometry does not change in the reaches), there will be only one critical water depth. Q 2Tc =1 gAc3

22.52 (4 + 2 yc ) =1 3 9.81 [(4 + yc ) yc ] yc =1.32 m

Reach A is a mild slope channel as y0A = 2.26 m > y0c = 1.32 m and the flow is subcritical. Reach B and C are steep slope channels as y0B, y0C < yc and the flow is supercritical on both reaches. Reach B is steeper than reach C. The various reaches are schematically shown in Fig. (6.24). The CDL is drawn at a height of 1.32 m above the bed level and NDLs are drawn at the calculated y0 values. The controls are marked in the figure. Reach A will have an M2 drawdown curve, reach B an S2 drawdown curve and reach C a rising curve as shown in the figure. It may be noted that the resulting profile as above is a serial combination of Fig. (6.22 d and f).

29

Prof. Dr. Atl BULU

Figure 6.23

30

Prof. Dr. Atl BULU