Chapter 4: Non uniform flow in open channels

Transcription

1 Chapter 4: Non uniform flow in open channels

2 Learning outcomes By the end of this lesson, students should be able to: Relate the concept of specific energy and momentum equations in the effect of change in bed level - Broad Crested Weir Relate the concept of specific energy and momentum equations in the effect of lateral contraction of channel ( Venturi Flume) UiTMKS/ FCE/ BCBidaun/ ECW301

3 Introduction Analysis of steady non uniform flow in open channels. Non uniform flow occurs in transitions where there is change in cross section or obstruction in channel. Analysis requires a different approach, requiring the use of the energy equation in a different form. UiTMKS/ FCE/ BCBidaun/ ECW301 3

4 Specific energy & alternative depths of flow Specific energy, E, E D+ (16.1) For a wide rectangular channel, mean velocity is, While the volume rate of flow per unit width, v g Q Q v A BD q Q B q D UiTMKS/ FCE/ BCBidaun/ ECW301 4

5 Substituting v & q into E, q E D+ gd q + g This equation has 3 roots: 1 root is negative & unreal D 3 ED (16.) (16.3) roots are positive & real, which give alternate depths: Larger depth: deep slow flow (subcritical/ tranquil/ streaming flow). 0 Smaller depth: shallow fast flow (supercritical/ shooting flow) UiTMKS/ FCE/ BCBidaun/ ECW301 5

6 UiTMKS/ FCE/ BCBidaun/ ECW301 6

7 UiTMKS/ FCE/ BCBidaun/ ECW301 7

8 Critical depth, D C : Depth at which the roots coincide q q max E E min To find D C : de dd 1 q gd 3 When de/dd 0, Critical depth, D C q g 1 3 Q gb 1 3 (16.4) UiTMKS/ FCE/ BCBidaun/ ECW301 8

9 3 Sub. q gd C from (16.4) to (16.), E Therefore, for rectangular channel, D C gd + gd Differentiating (16.3) & assuming E is constant 3 C 3 C 3 D C D C 3 E q dq dd D [ g( E D) ] g 1 ( E D) 1 1 D ( ) E D 1 (16.5) UiTMKS/ FCE/ BCBidaun/ ECW301 9

10 dq 1 When D DC, 0, ( E DC) DC 0, (16.5) dd becomes, q max 1 E 3 g E E 3 g E Or q max gdc (16.6) 3 UiTMKS/ FCE/ BCBidaun/ ECW301 10

11 Critical velocity: velocity of flow corresponding to critical depth. Sub. E 3 DC, D D C into (16.1), 3 D v C C D C ( gd ) C C v + g UiTMKS/ FCE/ BCBidaun/ ECW301 11

12 D C for non rectangular sections UiTMKS/ FCE/ BCBidaun/ ECW301 1

13 For any shape and cross sectional area the E for any D, v E D+ g Since v Q/A, Q E D+ A g (16.7)

14 For flow at D C and v C, E min, differentiating (16.7), 1 Q A g 3 da dd (16.8) But a change in depth will produce a change in cross-sectional area, therefore da/ddb. Q A 3 B g 1 0 (16.9) UiTMKS/ FCE/ BCBidaun/ ECW301 14

15 For critical flow, v C Q A Ag B 1 From (16.9) v C gd (16.10) where D is the average depth. UiTMKS/ FCE/ BCBidaun/ ECW301 15

16 Froude Number Assume a surface wave of height δz is propagated from left to right of observer. Wave is brought to rest relative to observer by imposing a velocity c equal to wave velocity on the observer, flow will appear steady.

17 UiTMKS/ FCE/ BCBidaun/ ECW ( ) ( ) ( ) ( ) ( ) (5.3) 0 Z B Z B Z wave) per unit time (right of Mass wave) per unit time(left of Mass u Z Z Z u c Z c u Z Z u u Z c u u c u u δ δ δ δ δ δ δ δ δ ρ δ δ ρ ( ) ( ) ( ) u c Z g u u δ δ δ δ ρ δ c BZ u - ZBZ g velocity Change of Mass per unit time Z Hydrostatic force due to

18 Sub. δu to (5.3), ( c u) δz gδz ( Z + δz) ( c u) ( c u) ( Z+ δz) gz If wave height δz is small, Velocity of the propation of the wave relative to the fluid c- u gz g UiTMKS/ FCE/ BCBidaun/ ECW301 18

19 Ratio of the stream velocity u to the propagation velocity c-u is known as Froude Number Fr. u Fr c u (5.34) Fr can be used to determine the type of flow for open channel. u gz UiTMKS/ FCE/ BCBidaun/ ECW301 19

20 For critical flow conditions, the Froude Number is, Fr If v < v c subcritical flow If v > v c supercritical flow Important difference, v gd Subcritical : disturbances can travel upstream and downstream thus enabling downstream conditions to determine the behavior of the flow. Supercritical: disturbances cannot travel upstream and thus downstream cannot control the behavior of the flow. v C gd C 1 UiTMKS/ FCE/ BCBidaun/ ECW301 0

21 Tranquil Critical Shooting Subcritical Critical Supercritical Depth D > D C D D C D < D C Velocity v < v C v v C v > v C Fr Fr < 1 Fr 1 Fr > 1 Channel slope Mild Critical Steep Control Downstream - Upstream Disturbance Wave can travel upstream Standing waves Waves cannot travel upstream UiTMKS/ FCE/ BCBidaun/ ECW301 1

22 Figure 4.1 shows that when flow is in the region of D c, small changes of E and q results in relatively large changes in D. Small surface waves are therefore easily formed but since velocity of propagation v p v c, these waves will be stationary or standing waves. Their presence therefore, an indication of critical flow conditions UiTMKS/ FCE/ BCBidaun/ ECW301

23 Figure 4.1a: Plot of D vs q for a constant E UiTMKS/ FCE/ BCBidaun/ ECW301 3

24 Line OA from figure below can be drawn at 45 through the origin. If scales for E and D are the same, horizontal distances from vertical axis to line OA will be equal to D, and the distance from OA to the specific energy curve will be v /g. If q is constant, Tranquil flow: D increases, v increases, E curve is asymptotic to OA Shooting flow: D decreases, v increases, E curve will be asymptotic to the E axis. UiTMKS/ FCE/ BCBidaun/ ECW301 4

25 Which of the two alternate depths for a given E will occur at a cross section depends on the slope of the channel. Critical slope, s c, is defined as the slope of the channel which will maintain flow at critical depth, DC. Uniform tranquil: s < s c mild slope Uniform shooting flow: s > s c steep slope UiTMKS/ FCE/ BCBidaun/ ECW301 5

26 Example 4.1 A rectangular channel 8 m wide conveys water at a rate of 15 m 3 /s. If the velocity in the channel is 1.5 m/s, determine; a) E b) D C c) v c d) E min e) Type of flow UiTMKS/ FCE/ BCBidaun/ ECW301 6

27 Example 4. Determine the critical depth in the trapezoidal channel shown below if the discharge in the channel is 0.34 m 3 /s. The channel has side slopes with a vertical to horizontal ratio of 1:1. UiTMKS/ FCE/ BCBidaun/ ECW301 7

28 Example 4.3 Determine the critical depth in a channel of triangular cross section conveying water at a velocity of.75 m/s and at a depth of 1.5 m. The channel has side slopes of 1:. UiTMKS/ FCE/ BCBidaun/ ECW301 8

29 Exercise A channel has a trapezoidal cross-section with a base width of 0.6 m and sides sloping at When the flow along the channel is 0 m 3 min -1, determine the critical depth. (0.7 m) UiTMKS/ FCE/ BCBidaun/ ECW301 9

30 Control sections Control sections cross sections at which the flow passes through the critical depth. Such sections are limiting factor in the design of channel. and some of the cases in which they occur are: Transition from tranquil to shooting flow Entrance to a channel of steep slope from a reservoir Free outfall from a channel with a mild slope Change in bed level or change in width of channel UiTMKS/ FCE/ BCBidaun/ ECW301 30

31 Transition from tranquil to shooting flow UiTMKS/ FCE/ BCBidaun/ ECW301 31

32 May occur where there is a change of bed slope s. Upstream slope is mild and s is less than the critical slope s c. Over s considerable distance the depth will change smoothly from D 1 to D. At the break in the slope, the depth will pass through D C forming a control section which regulates the upstream depth. At the tail end, the reverse transition from shooting to tranquil flow occurs suddenly by means of a hydraulic jump. UiTMKS/ FCE/ BCBidaun/ ECW301 3

33 Entrance to a channel of steep slope from a reservoir If depth of flow in the channel is less than DC for the channel, water surface must pass through DC in the vicinity of the entrance, since conditions in the reservoir correspond to tranquil flow.

34 Free Outfall from a Channel with a Mild Slope

35 If slope s of the channel is less than sc the upstream flow will be tranquil. At the outfall, theoretically, the depth will be critical, D C. In practice, gravitational acceleration will cause an increase of velocity at the brink so that D < D C. While experiments indicate that depending on the slope upstream: D C occurs at distance of between 3D C to 10D C from the brink. D at the brink is approximately 0.7D C. UiTMKS/ FCE/ BCBidaun/ ECW301 35

36 Change in bed level or channel width Flow over a broad-crested weir Effect of lateral contraction of a channel UiTMKS/ FCE/ BCBidaun/ ECW301 36

37 Flow over a broad-crested weir Broad-crested weir is an obstruction in the form of a raised portion of the bed extending across the full width of the channel with a flat upper surface or crest sufficiently broad in the direction of flow for the surface of the liquid to become parallel to the crest. Upstream edge is rounded to avoid separation losses that occur at a sharp edged. The flow upstream of the weir is tranquil and the conditions downstream of the weir allow a free fall over the weir. UiTMKS/ FCE/ BCBidaun/ ECW301 37

38 UiTMKS/ FCE/ BCBidaun/ ECW301 38

39 The discharge over the weir will be, therefore, be the maximum possible and flow over the weir will take place at D C. For a rectangular channel, D D C Q C Q Q gb ( ) BgD E 3 C B g E BE 3 (16.11) UiTMKS/ FCE/ BCBidaun/ ECW301 39

40 The specific energy, E, measured above the crest of the weir will be (assuming no losses), H is the height of the upstream water level above the crest and v is the mean velocity at a point upstream where flow is uniform. If the upstream depth is large compared with the depth over the weir, (v /g) is negligible, therefore, Rewriting (16.11), E H E v + g H Q1.705BH 3 (16.1) UiTMKS/ FCE/ BCBidaun/ ECW301 40

41 A single measurement of the head, H above the crest of the weir would then be sufficient to determine Q. Since D (, the depth over the crest of /gb ) 1 3 C Q the weir is fixed, irrespective of its height. Any increase in the weir height will not change D C but will cause an increase in the depth of the flow upstream. Therefore, maximum height of the weir, z D 1 D UiTMKS/ FCE/ BCBidaun/ ECW301 41

42 If the level of the flow downstream is raised, the surface level will be drawn down over the hump, but the depth may not fall to the critical depth. The rate of flow can be calculated by applying Bernoulli s s Equation and continuity equation and depends on the difference in surface level upstream and over the weir. UiTMKS/ FCE/ BCBidaun/ ECW301 4

43 Example 4.4 A broad crested weir 500 mm high is used to measure the discharge in a rectangular channel. The width of the channel is 0 m and the height of the channel upstream of the weir is 1.5 m. What is the discharge in the channel if water falls freely over the weir? Assume that the velocity upstream is very small. Determine the difference in water level between upstream and over the top of the weir. UiTMKS/ FCE/ BCBidaun/ ECW301 43

44 Effect of lateral contraction of a channel When width of a channel is reduced while bed remains flat, q increases. As channel narrows - neglecting losses, E remains constant for tranquil flow, depth will decrease while for shooting flow depth will increases.

45 Free surface does not pass through D C Lateral contraction Free surface passes through D C With hump UiTMKS/ FCE/ BCBidaun/ ECW301 45

46 Free surface does not pass through D C UiTMKS/ FCE/ BCBidaun/ ECW301 46

47 Arrangement forms a venturi flume (venturi meter). Applying energy equation between upstream & throat and ignoring losses, UiTMKS/ FCE/ BCBidaun/ ECW ) (applying BD D B gh D B Q Dv B Q BD D B gh v Q Q h D D D B D B g v g v D g v D

48 Actual discharge, Q CdB D 1 gh BD BD 1 1 C d is a coefficient of discharge 0.95 to 0.99 UiTMKS/ FCE/ BCBidaun/ ECW301 48

49 Free surface passes through D C UiTMKS/ FCE/ BCBidaun/ ECW301 49

50 The flowrate is given by, Q B B B Dv D ( E D ) Assuming that the upstream velocity head is negligible, g E BE 3 1 g E 3 Q 1.705BH (16.15) where H is the of the upstream free surface above bed level at the throat. C 3 UiTMKS/ FCE/ BCBidaun/ ECW301 50

51 Lateral contraction with hump UiTMKS/ FCE/ BCBidaun/ ECW301 51

52 Height of upstream water level above the hump, H D 1 Z When upstream conditions are tranquil and the bed slope is the same downstream as upstream, impossible for shooting flow to be maintained for any great distance from the throat. Revert to tranquil flow downstream by means of a hydraulic jump or standing wave. Venturi flume operating in this mode is known as standing wave flume. UiTMKS/ FCE/ BCBidaun/ ECW301 5

53 Example 4.5 A venturi flume is constructed in a channel which is 3.5 m wide. If the throat width in the flume is 1. m and the depth upstream from the constriction is 1.5 m, calculate the discharge in the channel when the depth at the throat is 1. m. If the conditions are such that a standing wave is formed, what is the discharge? UiTMKS/ FCE/ BCBidaun/ ECW301 53

54 Example 4.6 A Venturi flume is.5m wide and 1.4m deep upstream with a throat width of 1.3m. Assuming that a standing wave form downstream, calculate the rate of flow of water if the discharge coefficient is Do not ignore the velocity of approach. UiTMKS/ FCE/ BCBidaun/ ECW301 54

55 Review of past semesters questions

56 APR 010 A 10 m wide channel conveys 5 m3/s of water at a depth of 1.6 m. Determine : i) specific energy of the flowing water ii) critical depth, critical velocity and minimum specific energy UiTMKS/ FCE/ BCBidaun/ ECW301 56

57 APR 010 A venturi flume is 1.40 m wide at the entrance and 0.7 m wide at the throat. Determine the flow if the depths at the entrance and at the throat is 0.8 m and 0.6 m respectively. Neglect all losses. UiTMKS/ FCE/ BCBidaun/ ECW301 57