Open Channel Flow Basics
Open Channel Flow Basics
CHAPTER I
Introduction
Open channel flow occurs where ever the flow proceeds with the
liquid surface exposed to constant pressure.
This is the most complex flow type, requiring the solution of energy, momentum
and friction equations through time.
Types of Open Channel
  Flume: a channel built (or supported) above the ground to convey fluid
   from one point to another. In the field flumes are made of concrete,
   wood, sheet metal or masonry. Laboratory flumes are usually made of
   wood, metal, glass or a composite of these materials.
  WHY?????
velocity distribution in open channels
Determination of energy and
momentum coefficients
 To determine the values of α and β the
  velocity distribution must have        been
  measured (or be known in some way).
           dA
             V13 A1  V23 A2  V33 A3  Vi Ai
             3
           u                              3
       3                           3
        V A       3
               V ( A1  A2  A3 )     V  Ai
         Q V1 A1  V2 A2  V3 A3
      V                        
                                      V Ai       i
A A1  A2  A3  A i
                                     2
          2
         V A        2
                 V ( A1  A2  A3 )     V  Ai
Uniform flow and the Development of
Friction formulae
 When uniform flow occurs gravitational
 forces exactly balance the frictional
 resistance forces which apply as a shear force
 along the boundary (channel bed and walls).
                          0.00281 1.811
                 41.65            
             C               S        n
                             0.00281 n
                1   41.65         
                                S    R
 where n is known as Kutter’s n. The above formula gives C in British units,
 which could be converted into metric units (exercise)
The Manning equation
     R   2/3
               So              1 A5 / 3
V                          Q
                               n P 2 / 3 S o1 / 2
          n
a few typical values of Manning's n
      Q  AC RS o                          A5 / 3
                      So, K  ACR1/ 2   
      Q  KS o1 / 2                       nP 2 / 3
    Use of conveyance may be made when
    calculating discharge and stage in compound
    channels and also calculating the energy and
    momentum coefficients in this situation.
Computations in uniform flow
 We can use Manning's formula for discharge to calculate steady
  uniform flow.
 Two calculations are usually performed to solve uniform flow
  problems.
1.   Discharge from a given depth
2.   Depth for a given discharge
 In steady uniform flow the flow depth is know as normal depth.
 As we have already mentioned, and by definition, uniform flow
  can only occur in channels of constant cross-section (prismatic
  channels) so natural channel can be excluded.
 However we will need to use Manning's equation for gradually
  varied flow in natural channels - so application to natural/irregular
  channels will often be required.
Optimal Shape of Cross-Section
 The most hydraulically-efficient shape of channel is the one which
  can pass the greatest quantity of flow for any given area or,
  equivalently, the smallest area for a given quantity of flow.
 Rapid change can also occur when there is a change from super-
  critical to sub-critical flow (see later) in a channel reach at a
  hydraulic jump.
The application of Energy Equation
                                      p V 2
 Recalling the Bernoulli equation,     
                                     g 2 g
                                              z  cons tan t
Bernoulli
 And assuming a hydrostatic pressure distribution we can write
  the pressure at a point on a streamline, A say, in terms of the
  depth d (the depth measured from the water surface in a
  direction normal to the bed) and the channel slope.
                         pA =ρgy2
                                            y2
    In terms of the vertical distance d          y1 cos 
                                           cos               So, p A  gy1 cos 2 
                                      y 2  y1 cos 2 
 Thus we have a cubic function with the only unknown being the
 downstream depth, y2. There are three solutions to this but only one is
 correct for this situation. We must find out more about the flow before we
 can decide which it is.
Specific Energy
 Specific   energy, Es, is
  defined as the energy of the
  flow with reference to the                 Fig. 1.10 Specific energy
  channel bed as the datum.                     (definition sketch)
 The concept of specific
  energy was first introduced
  by Bakmenteff (1918).
 With reference to the figure
  the total energy of flow
  with respect to the channel
  bottom is given by
                            p1 V12
      E s1    ( y  y1 )               V12
                            g 2 g    y
                                          2g
Specific Energy…
 Thus the specific energy at an open channel section is equal to the sum of
  the flow depth and the velocity head.
 In the above equation V1 denotes the velocity of flow at the point of
  interest, in the figure above at point 1. In practice it is easier to use the
  average velocity of flow at the section and speak about the specific energy
  of the flow at a section.
 How ever the velocity of flow changes from point to point with in the flow
  and as a result the specific energy changes from stream line to stream line.
  It is common to use the average velocity of flow with a correction factor.
 The specific energy computed using the average velocity is taken to apply
  for all points in the section, i.e. is taken as the specific energy of the
  section. For steady flow this can be written in terms of discharge Q
                (Q / A) 2                                        q 2
      Es  y                                          Es  y 
                   2g                                             2 gy 2
 All point in the channel between point 1 and 2 must lie on the
  specific energy curve between point A and B or B'.
 To reach point B' then this implies that Es1 - Es2 > ∆z which is not
  physically possible. So point B on the curve corresponds to the
  specific energy and the flow depth at section 2.
Critical, Sub-critical and super critical flow
 The specific energy change with depth was plotted above for
  a constant discharge Q, it is also possible to plot a graph with
  the specific energy fixed and see how Q changes with depth.
For a fixed discharge:
                              Vc2       y
                 E sc  y c       yc  c
                              2g         2
                       2
                 y c  E sc
                       3
Variation of the Discharge with depth for a given specific
                       energy value
The Froude number
                                                     V
 The Froude number is defined          FN 
                                                    gDm
  for channels as:
 Its physical significance is the             Inertial Force
                                     FN2 
  ratio of inertial forces to                Gravitational Force
                    Q                        Q
               Q                  Q
                   y1b                      y2b
                                y 2  y1 3
                        E 
                                 4 y1 y 2
These are useful results and which can be used in gradually varied flow
calculations to determine water surface profiles.
The higher the upstream Froude number the higher the jump and the
greater the loss of energy in the jump.
         CHAPTER FIVE
         5.0 RAPIDLY VARIED FLOW
 Characteristics of RVF
 STEADY NON-UNIFORM FLOW
 Pronounced curvature of streamlines.
 Abrupt change of flow profile (virtually broken)
 Resulting the state of high turbulance
 Example: Hydraulic Jump
   RVF vs GVF
   In view of contrast with UF & GVF the following characteristics should be noted.
   Pronounced curvature  hydrostatic pressure distribution can not be assumed
   Rapid variation in flow regime takes place in a very short distance
   Effect of boundary friction is comparatively small, which would play a primary role in a
    GVF
   In RVF the velocity-distribution coefficients  and  are much greater than unity and can
    not be accurately determined.
   Flow is actually confined by separation zones as well as solid boundaries. (Because profiles
    could be broken).
   Approach to the problem
   The theory that assumes:
   Parallel flow,
   Hydrostatic distribution of pressure- Does not apply in RVF computation.
   For RVF of continuous flow profile a mathematical equation can be established,
   Approach to the solution of such equation
   Graphical method (e.g. flow-net analysis)
   Numerical method (e.g. method of relaxation)
   No satisfactory general solution has yet been obtained
   Practical approach
   No general solution yet been found
   Various RVF phenomena are treated as isolate cases
   Each with own semi-empirical/empirical treatment
     Experimental results are used empirically
     Flow interpreted qualitatively using energy principle, momentum principle, geometry
       plus sometimes dimensional analysis
   Three isolated cases
     Flow over spillway
     Hydraulic jump
     Flow under gate
 Flow over spillways
   Definition:
   Spillway: is a structure over or through a dam for discharging flood flows; overflow
    channel; opening built into a dam or the side of a reservoir to release (to spill) excess
    floodwater.
 SHARP-CRESTED WEIR (SCW) VS BROAD CRESTED WEIR (BCW)
   BCW
            Overflow structure with horizontal crest
            above which the deviation from a hydrostatic pressure distribution
            because of centripetal acceleration may be neglected.
    stream-lines are parallel and straight
   Criteria 0.5  H1/L0.07
     If 0.07 H1/L the energy loss above the crest can not be neglected
     0.5  H1/L so that the hydrostatic pressure distribution can be assumed
   L = length of the weir crest in the direction of flow,
   H1 total energy head over the weir crest
   SCW
     Overflow structure (H1/L > 15)
     The crest length in the direction of the flow is short enough not to influence the H-Q
      relationship of a weir
     In practice, 0.002m L so that even at a minimum head of 0.03m the nappe is
      completely free from the weir body after passing the weir  no adhered nappe can
      occur
     An air pocket beneath the nappe form from which a quantity of air is removed
      continuously by the over falling jet.
     Therefore, Precaution is required not to ensure that the pressure in the air pocket is not
      reduced. Otherwise resulting undesirable effects:
        Owing to the increase of the under pressure the curvature of the over falling jet will
         increase, causing increase of the discharge coefficient
        Irregular supply of air to the pocket will cause vibration of the jet resulting an
         unsteady flow
 SCW is the simplest form of overflow spillway
 Motto:
   Spillways must discharge the peak flow under smallest possible head.
   Negative pressure on the crest must be limited to avoid danger of cavitation on the
    crest or vibration of the structure.
   Theoretically, there should be atmospheric pressure on the crest
   In same time t, the particle will travel1 a vertical distance y (taking y is positive downward)
                     y   Vo t Sin   gt 2
                                        (2)2
   Move the origin up so that it coincides with the peak.
                                         1
                y   Vo ts in           gt 2  C '
                                         2
                                                                         x    
                                                                t           
                                                                      V Cos 
        Eliminating t from 1 & 2, from 12                            o        
                     x                1  x 
          y   Vo           Sin   g           C '
                     V0 Cos          2  Vo Cos  
hv   Dividing each term by the total head H above the crest
H
         Y           X       gH                             C'
             tan                           HX 2   
         H           H   2Vo2 Cos 2                        H
 
                       gH                                           C'
     Let
              A    , 2 CosB=
                            2 - tan, &                          C 
                   2V 0                                             H
                             2
              Y    X     X
                 A   B   C
     Then,   H     H
                         H
      = General equation for the lower surface of the nappe in dimension less term
     Since the horizontal velocity component is constant, the vertical thickness of the nappe T may be
      assumed constant and          T
                         D 
                                   H
      Adding a term               to the above equation the
                                                         2   general equation for the upper surface of
      the nappe is                           Y      X        X
                                                 A   B        CD
                                             H       H
                                                             H
     These equations are quadratic hence, the nappe Surfaces are theoretically parabolic.
     Several experimental studies on the nappe over a sharp-crested weir have been made.
       Reputed works has been done by US Bureau of Reclamation, they developed the following
         equations for the constants in
                                     hv the general nappe equations.hv          hv 
                                                                                      2
                                                                                         hv
                                      H                             H    1.568   0.892  0.127
            A = - 0.425 + 0.25       B = 0.411 - 1.603 -                      H       H
                           hv
                          H
     C = 0.150 – 0.45           D = 0.57 – 0.02 (10m) 2 exp (10m)
   Where:
         hv = the velocity
                          hv head of the approach flow
         m= - 0.208 H
   For high weirs, the velocity of approach is relatively small and can be ignored (hv  0)
          A = 0.425
          B = 0.055
          C = 0.150
          D = 0.559
   Experimental data have indicated that these equations are not valid When,
        X                             hv
      H  < 0.5  and that          >  0.2
                                      H
   For < 0.5, The pressure with in the nappe in the Vicinity of the weir crest is > Patm
   because of the convergence of the streamlines.
   Consequently, forces other than gravity are acting on the nappe,
   which makes the principle of the projectile invalid.
   N.B: The above theory and equations apply only if the approach flow is sub critical.
   For Supercritical flow, or Fr >1, the nappe profile becomes essentially a function of the
    Froude number rather than a function of the boundary geometry as described above.
 Aeration of the Nappe
   In the preceding discussion the over falling nappe is considered a crated; i.e., The upper
    and lower nappe surfaces are subject to full atmospheric pressure.
   In practice,
        - Usually insufficient aeration below the nappe occurs due to removal of air by over
    falling set.
    Effects of reduction of pressure
      Increase in pressure difference on the spillway itself
      Change in the shape of the nappe for which the spillway crest is designed
      Increase in discharge, sometimes accompanied by fluctuation or pulsation of the
         nappe, which may be very objectionable if the weir or spillway is used for measuring
         purposes.
      Unstable performance of the hydraulic model
 Crest Shape of Overflow Spillways
   Earliest shapes were based on a simple parabola designed to fit the trajectory of the falling nappe (the
    equation for the lower surface of the nappe).
   Bazin’s made comprehensive laboratory investigation for nappe shapes.
   the used of Bazin’s data in design will produce a crest shape that conincides with lower surface of as
    aerated nappe over a sharp-crested weir.
   Such a profile is known as Bazin profile Advantage
      Should cause no negative pressure on the crest (the presence of negative pressure will lead to
        danger of cavitation damage).
      In selecting a suitable profile avoidance of negative pressure should be considered an objective,
      along with such other factors as maximum hydraulic efficiency, practicability, stability &
        economy.
   Extensive experiments on the shape of the nappe over-sharp crested weir were conducted
    by U. S Bureau of Reclamation; including Bazin’s,
   The Bureau has developed coordinates of the nappe surface for various slope faced weirs
   On the basis of the Bureau data, The U.S Army. Corps of Engineers has developed several
    standard shapes at its Waterways Experimental Station.
   Such shapes designed as the WES standard spillway shapes, can be expressed by the
    following equation:-
                                                          n 1
                                    X n  K Hd                   Y
     X and Y are Coordinates of the crest profile with the origin at the highest point of the
      crest.
     Hd is the design head excluding the velocity head of the approach flow
   K & n are parameters depending on the slope of the upstream face. values of k & N are
    given as flows:
   Slope of upstream face              k                n
   Vertical                    2.000             1.850
   3 :1 (V = H)                1.936             1.836
   3:2                        1.939             1.810
   3:3                       1.873              1.776
   For intermediate slopes: approximate value of k and n may be obtained by plotting the
    above values against the corresponding slopes and interpolating
   from the plot the required values for any given slope within the plotted range.
   The upstream face of the spillway crest may some times be designed to set back, as shown
    by the dashed lines
   Discharge of WES Spillway
   The discharge over a spillway can be computed by an equation in the form of SCW/BCW
   Q = CLHe1.5
   He- the total energy head on the crest, including the velocity head in the approach canal.
   The effect of the approach velocity is negligible when height h of the spillway is greater
    than 1.33Hd (h > 1.33 Hd), where the design head exclude the approach velocity head.
   Under this condition, i.e. h/Hd > 1.33, He = Hd can be taken (the approach velocity head is
    negligible) and the coefficient of discharge C has been found to be C = 2.21 (if is in ft C He =
    4.03)
                              y1  2              
   This has been verified with experiments
   Types of Jump
   Hydraulic Jumps on horizontal floor are of several distinct types. They can be conveniently classified
    according to Froud Number Fr1 of the incoming flow as follows.
   Fr1 =1          critical flow no jump can form
   1< Fr1 < 1.7       the water surface shows undulation (undular jump)
   1.7 < Fr1 < 2.5 a series of small rollers develop on the surface of the jump,
                  but the d/s water surface remains smooth.
                  The velocity throughout is fairly uniform, and the energy loss low.
                  The jump is called weak jump.
    2.5 < Fr1 < 4.5 there is an oscillating jet entering the jump bottom to surface and back again with no
                                                     periodicity.
                 Each oscillation produces a large wave of irregular period which, very common in canals,
                     can travel for miles doing unlimited damage to earth banks and ripraps.
              This jump is called Oscillating Jump.
   4.5 < Fr1 < 9.0 steady Jump:-
   The down stream extremity of the surface roller and the point at which the high-velocity jet tends to
    leave the flow occur at practically the same vertical section.       The action and position of this jump
    are least sensitive to variation in tail-water depth.
    The jump is well balanced and the performance is at its best.
      The energy dissipation ranges from 45 to 70%.
    Fr > 9.0 Strong jump:-
    The high-velocity jet grabs intermittent slugs of water rolling down the front face of the
     jump,
    generating waves down-stream and a rough surface can prevail.
    The jump action is rough but effective since the energy dissipation may reach 85%.
    N.B. It should be noted that the ranges of the Froude Number given above for the various
     types of jump are not clear-cut but overlap to a certain extent depending on local
     conditions.
                                                                        4 y1 y 2
                                                                                 E
                                             Re lative loss : the ratio
                                                                                 E1
    Efficiency: the ratio of the specific energy after the jump to that before the jump is defined
     as the efficiency of the jump.
                                              E2
                                                 
                                                         2
                                                               3
                                                   8F1  1 2  4 F1  1
                                                                             2
                                              E1         2
                                                      8F1 2  F1
                                                                  2
                                                                            
   The relative loss is equal to ,                                                   E2
                                                                               1
                                                                                      E1
   this also is a dimensionless function. of Fr1.
 Height of Jump:- the difference between the depths after and before the jump.
   Hj = y2 – y1
   Expressing each term as a ratio with respect to initial specific energy.                                y2
                                                                                                hj
                                                                                                                       y1
                                                                 h1 y2   y
                                                                         1                                 E1
                                                                                                E1
                                                                                                                       E1
                                                                  E1 E1 E1
   Where          is the relative height,
     is the relative initial depth, and
                                                                                                                  2
                                                                                         hj           1  8 F1        3
                                                                                                            2
     is the relative sequence depth.                                                    E1               F1  2
 Length of Jump:
   The length of a jump (also length of stilling basin) is empirically given as,
              L= K(y2-y1)
   Where, k – is a coefficient derived from laboratory and filed experiment. 4.5 < k < 5.5 where the lower k = 4.5 applies
    of Fr2 > 10 and the higher k=5.5 for Fr2 < 3.
   Gates
     Flow inunder
              canals are mainly used as water level regulators.
                     Gates
   Sometimes, gates are used as discharge regulator (measuring device).
   They are under-shot or underflow structures. Example slice gate, radial gate roller gate
   The design of underflow gate focuses on head-discharge relationship (Q-H).
   The objective is to minimize head loss; this means that the gate has to be lifted out off the
    water for design discharge.
   The other concern of the design is the pressure distribution over the gate as a function of
    opening and gate form.
   The H-Q relationship for gate depends on the shape and dimension of the control section
    and the resulting curvature of the streamlines.
   For gated structures the control section is defined by the vena contract, being the smallest
    cross section just down steam of the gate.
    In the vena contract, streamlines are straight and parallel.
   In gate flow 3 flow types can be distinguished.
 1) Free flow: the opening is relatively small (h1/a >2), and the contraction of the steam-
  lines in vertical direction is strong.
 The down stream water level (h2) won’t affect the flow underneath the gate and a
  hydraulics jump will occur down stream of the vena contra.
 The discharge depends up on the gate opening the upstream water level and the
  contraction coefficient.
 Submerged flow: the d/s water level influences the flow underneath the gate. The
    hydraulic jump is drowned and the jet underneath the gate is submerged. The discharge
    depends upon the upstream and downstream water level and the gate opening.
   The boundary between free and submerged flow is a sharp one, which can be cleanly found
    from the gate opening and the two water levels.
 Weir flow: on off gate
   The equation for a free flow underneath a sharp edged gate is:
                Q  Cd Ba 2gh1
   Cd   = discharge coefficient
   B    = Width of gate opening
   a    = height of gate opening
   h1   = upstream water depth
                                                Cc
                                     Cd 
   The discharge coefficient Cd                      a
                                             1  Cc
                                                      h1
   Where,
   CC   = Contraction coefficient of the jet depending on the shape of the gate and
   d    = diameter of the rounded bottom edge
   For       d
                 4.7  Cc  0.99 Rounded edged gates 
              a
   For submerged flow, some equations include the difference between the upstream and
    downstream depths and others use the upstream water level only. The general equation is
    offer given as.                       Q  C Ba 2gh     2           1
   Where,
   a     = vertical opening of the gate( a< 0.67h1)
   h1    = Upstream water depth
   B     = Effective width of the opening
   C2 = discharge    coefficient. h h
                  a  0.67 * h1 or 1  0.67 1  1  1.5
   For values of                           a   athe discharge follows from the equation for a broad-
.
SHARP CRESTED WEIR
BROAD CRESTED WEIR
Sketch of
a Venturi
flume
Irrigation Canal
Grand Canal of China
Beijing and Hangzhou
Hydroelectric Power Plant
SPILLWAY