Go to ⇒ Presentation format 
1.1 OPENCHANNEL FLOW
Openchannel flow has a free surface and it is, therefore, subject to atmospheric pressure
In closedconduit flow, the cross section is fixed by the pipe boundaries. On the other hand, in openchannel flow, the flow cross section is not fixed, varying with the flow. In closedconduit flow, the roughness varies from smooth brass to corroded pipes; in openchannel flow, it varies from acrylic glass or lucite^{®} (a very smooth type of plastic), to that of natural stream channels and their neighboring flood plains.
In closedconduit flow, the
hydraulic pressure at the center of the pipe defines the hydraulic grade line (HGL in Fig. 12).
The hydraulic pressure (head of water) above the centerline of the pipe is referred to as the piezometric head.
The energy grade line includes the velocity head
In openchannel flow, the flow depth measured above the
channel bottom defines the water surface elevation,
which is equivalent to the hydraulic grade line of closedconduit flow; see Fig. 13.
The total energy grade line includes the velocity head
Note the difference between closedconduit and openchannel flow. In closedconduit flow, water will rise in a piezometer tube up the level where it defines the hydraulic grade line associated with the hydraulic pressure in the conduit. On the other hand, in openchannel flow, the water surface is the hydraulic grade line, which is at atmospheric pressure. 1.2 TYPES OF FLOW
There are two general types of openchannel cross sections:
Artificial, or humanmade channels, are usually prismatic, featuring a constant shape and size, at least for a certain length of channel. Conversely, natural channels are typically nonprismatic, i.e., the shape and size of the cross section varies along the channel. Artificial channels are also referred to as canals. Several geometric and hydraulic properties help describe an open channel (Fig. 14). These are:
In prismatic channels, the flow depth d is often referred to as y, particularly when it cannot be confused with stage. Also, the channel side slope is often referred to as z H : 1 V, particularly when it cannot be confused with bed elevation.
Openchannel flow may be classified as follows:
The flow is steady when the hydraulic variables (discharge, flow area, mean velocity, flow depth, and so on) do not vary in time. Conversely, the flow is unsteady when the hydraulic variables vary in time and space. Steady flow is relatively simple to calculate, compared to unsteady flow. The flow is uniform when the channel is prismatic and the hydraulic variables (Q, A, V, d, and so on) are constant in time and space. The flow is in equilibrium when the channel is nonprismatic and the hydraulic variables are approximately constant in time and space. The calculation of uniform flow is relatively straight forward when compared to that of other states of flow. The flow is gradually varied when the discharge Q is constant but the other hydraulic variables (A, V, d, and so on) vary gradually in space. Under gradually varied flow, the pressure distribution in the vertical direction, normal to the flow, is very close to hydrostatic, i.e., proportional to the flow depth. The flow is rapidly varied when the discharge is constant but the other hydraulic variables (A, V, d, and so on) vary rapidly in space, in such a way that a hydrostatic pressure distribution cannot be assumed in the vertical direction normal to the flow. While the calculation of gradually varied flow is somewhat involved but doable, the calculation of rapidly varied flow is generally more complex, in practice being based on empirical formulas, for lack of a theoretical solution. The flow is spatially varied when the discharge Q varies in space only, i.e., along the channel, typically due to lateral inflow or outflow.
Steady uniform flow occurs in a prismatic channel (Fig. 15); steady equilibrium flow occurs in a nonprismatic channel. Unsteady uniform flow does not exist in nature, because the flow cannot be uniform and unsteady at the same time. The word "unsteady" implies nonequilibrium; thus, unsteady equilibrium flow does not exist.
Steady gradually varied flow is represented by the water surface profiles, also referred to as backwater (or drawdown) computations (Chapter 7). Unsteady gradually varied flow is the calculation of flood flows, or flood routing (Chapter 10). Steady rapidly varied flow is represented by the flow over spillways or the hydraulic jump. Unsteady rapidly varied flow is represented by the moving hydraulic jump, surges, roll waves, kinematic shocks, and tidal bores. Figure 16 shows a train of roll waves in a steep irrigation canal. Spatially varied flow occurs in an artificial canal when the discharge is varying along the channel, due to lateral water extractions or channel overflow.
1.3 STATE OF FLOW
The state of openchannel flow may be described in terms of certain characteristic velocities and
viscosities.
Velocity is the ratio of length (distance) over time, with units
There are three characteristic velocities in openchannel flow:
The mean velocity of the steady uniform flow using the Manning equation (SI units) is:
in which n = Manning's friction coefficient, R = hydraulic radius, and S = friction slope. The mean velocity of the steady uniform flow using the Chezy equation is:
in which C = Chezy friction coefficient. In general, four forces are active in a control volume in openchannel flow. These forces are due to friction, gravity, the pressure (flow depth) gradient, and inertia. Kinematic waves are those where the momentum balance is expressed in terms of the frictional and gravitational forces only (Lighthill and Whitham, 1955). The celerity of kinematic waves, or Seddon celerity, is (Seddon, 1900; Chow, 1959; Ponce, 1989):
in which β = exponent of the dischargeflow area rating, defined as follows:
Dynamic waves are those where the momentum balance is expressed in terms of the pressuregradient and inertial forces only. The celerity of dynamic waves is:
in which g = gravitational acceleration, and D = hydraulic depth, D = A /T. From Eq. 13, the relative celerity of kinematic waves is:
From Eq. 15, the (absolute value of the) relative celerity of dynamic waves is:
For rectangular channels, for which D = d, or for hydraulically wide channels, for which D ≅ d, the relative celerity of dynamic waves is:
Equation 18 is known as the Lagrange (relative) celerity equation, after Lagrange (1788), who first derived it.
The Froude number is defined as follows (Chow, 1959):
The Froude number characterizes the condition of:
Under subcritical flow, surface waves (perturbations) can travel upstream, because their upstream celerity w is greater than the mean flow velocity u. Under critical flow, surface waves (perturbations) remain stationary, because their (absolute) celerity w is equal to the mean flow velocity u. Under supercritical flow, surface waves (perturbations) can travel downstream only, because their upstream celerity w is smaller than the mean flow velocity u.
The Vedernikov number is defined as follows (Vedernikov, 1945; 1946; Powell, 1948; Craya, 1952):
The Vedernikov number characterizes the following states of flow:
Under stable flow, the relative kinematic wave celerity v is smaller than the relative dynamic wave celerity w and, therefore, surface waves (perturbations) are able to attenuate (dissipate). Under neutrally stable flow, the relative kinematic wave celerity v is equal to the relative dynamic wave celerity w and, therefore, surface waves (perturbations) neither attenuate nor amplify. Amplification amounts to negative dissipation. Under unstable flow, the relative kinematic wave celerity v is greater than the relative dynamic wave celerity w. Therefore, surface waves (perturbations) are subject to amplification. In practice, the condition V ≥ 1 leads to the development of roll waves, a train of waves that travel downstream, typically in artificial channels of steep slope (Cornish, 1907) (Fig. 17).
The three velocities u, v, and w lead to only two independent velocity ratios, the Froude (Eq. 18) and Vedernikov (Eq. 19) numbers. The third ratio:
is the dimensionless relative kinematic wave celerity, equal to the exponent of the dischargearea rating minus 1. Thus, it is seen that the exponent β in Eq. 14 is a function of both the Froude and Vedernikov numbers. The value of β varies with the type of friction regime (laminar, transitional, or turbulent; and turbulent Manning or Chezy) and crosssectional shape. For laminar flow, β = 3. For turbulent flow, under Manning friction: 1 ≤ β ≤ 5/3, depending on the shape of the cross section. For turbulent flow, under Chezy friction: 1 ≤ β ≤ 3/2, depending on the shape of the cross section.
There are three asymptotic crosssectional shapes in open channels:
Neutrally stable flow For neutral stability: V = 1. Therefore, from Eq. 111, the Froude number corresponding to neutrally stable flow is:
Table 11 shows values of F_{ns} for selected values of β. It is seen that as β varies from β = 3 (laminar flow) to β = 1 (inherently stable channel), the value of F_{ns} varies from F_{ns} = 1/2 to F_{ns} = ∞. In other words, as β ⇒ 1, F_{ns} ⇒ ∞. In practice, since friction has effectively a lower bound, the Froude number is restricted to an upper bound, which is seldom likely to exceed F = 25. Therefore, in most cases, a value of β = 1.04 would be already stable for practical purposes.
As shown in Table 11, values of β in openchannel
and overland flow are limited in the range
There are three characteristic viscosities in openchannel flow:
The kinematic viscosity ν of the fluid varies as a function of temperature (Appendix A). The concept of hydraulic diffusivity ν_{h} is due to Hayami (1951). Hayami combined the governing equations of openchannel flow (Chapter 10) to develop a single convectiondiffusion equation, i.e., an equation describing the convection (firstorder process) and diffusion (secondorder process) of a flood wave. The hydraulic diffusivity is defined as follows:
in which q_{o} = equilibrium unitwidth discharge, and S_{o} = friction (energy) slope. It is seen that flood wave diffusion is directly proportional to unitwidth discharge and inversely proportional to friction (energy) slope. Equation 113 can be expressed in terms of velocity and flow depth as follows:
A related value of diffusivity, which is independent of slope, is:
In general, for an arbitrary crosssectional shape:
in which R_{o} = hydraulic radius. In kinematic wave theory, the characteristic reach length is defined as follows (Lighthill and Whitham, 1955):
in which L_{o} is the length of channel in which the equilibrium flow drops a head equal to its depth. Thus, in terms of the characteristic reach length, the hydraulic diffusivity is:
In a manner resembling the hydraulic diffusivity, the wave diffusivity is conveniently defined as follows:
in which L = wavelength of the disturbance.
The Reynolds number R is (Chow, 1959):
The Reynolds number R describes the flow regime as either:
Under steady flow conditions in openchannel flow, laminar flow occurs for R ≤ 500 and turbulent flow for R > 2000. Transitional flow occurs in the intermediate range: 500 < R ≤ 2000. Under unsteady flow, the mixed laminarturbulent flow described in Table 11 is akin to transitional flow, featuring a comparable range of Reynolds numbers. In practice, most openchannel flow cases are in the turbulent regime. Conversely, most overland flow cases (i.e., freesurface plane flow) are in the laminar or mixed laminarturbulent regime.
The dimensionless wavenumber σ is defined as follows (Ponce and Simons, 1977):
The wavenumber σ describes the dimensionless length scale of the wave, as shown in Fig. 111, in terms of: (a) kinematic waves, (b) dyamic waves, and (c) mixed kinematicdynamic waves. Figure 11 is applicable for the case of Chezy friction in hydraulically wide channels. Under kinematic flow, depicted on the left side of Fig. 111:
Under dynamic flow, depicted on the right side of Fig. 111:
Under mixed kinematicdynamic flow, depicted by the middle section of Fig. 111:
Dynamic hydraulic diffusivity
The dynamic hydraulic diffusivity, which unlike the Hayami hydraulic diffusivity of
For low Vedernikov numbers, V ⇒ 0, the dynamic hydraulic diffusivity reduces to the expression for kinematic hydraulic diffusivity, i.e., Eq. 118. Conversely, for high Vedernikov numbers, V ⇒ 1, and the dynamic hydraulic diffusivity vanishes. Under this flow condition, the total absence of wave attenuation is conducive to the development of roll waves (Figs. 16 and 17).
1.4 FLOW REGIMES
The flow regimes in openchannel flow are:
The flow regimes are characterized by the Reynolds number R, Eq. 120. In openchannel flow, the laminar regime prevails for R ≤ 500, the transitional regime for 500 < R ≤ 2000, and the turbulent regime for R > 2000. The flow regimes vary with roughness of the channel surface. Figure 112 shows the relation between Reynolds number R and DarcyWeisbach friction factor f for flow in smooth channels. Figure 113 shows the relation between Reynolds number R and DarcyWeisbach friction factor f for flow in rough channels. The DarcyWeisbach friction formula, developed in connection with flow in pipes, is:
in which h_{f} = frictional head loss; f = DarcyWeisbach friction factor; L = length of the pipe; d_{o} = pipe diameter; V = mean flow velocity in the pipe; and g = gravitational acceleration.
The examination of Figs. 112 and 113 enables the following conclusions:
Note that the Prandtlvon Karman equation can be expressed in explicit form as follows:
The DarcyWeisbach formula, Eq. 123, is strictly applicable to closedconduit (pipe) flow. In pipe flow, the characteristic frictional length is the pipe diameter d_{o}. On the other hand, in openchannel flow, the characteristic frictional length is the hydraulic radius R, i.e., the ratio of flow area to wetted perimeter:
Since the flow area of a circular pipe (flowing full) is A = π d_{o}^{2}/4, and the wetted perimeter is P = π d_{o}, it follows that the hydraulic radius is equal to 1/4 of the pipe diameter or, conversely, that the diameter is equal to 4 times the hydraulic radius. Therefore, the DarcyWeisbach formula applicable to openchannel flow is:
in which V = mean flow velocity in the channel. In openchannel flow, the energy slope, which under steady flow is the same as friction, bed, or bottom slope, is:
For any arbitrary crosssectional shape, the Froude number is:
in which D = hydraulic depth: D = A /T. Equation 130 can be expressed in terms of the Froude number as follows:
Equation 132 states the proportionality between energy slope and Froude number, with the proportionality factor being a function of the DarcyWeisbach friction factor and the shape factor D /R. For a hydraulically wide channel, for which D ≅ R, Eq. 132 reduces to:
Thus, for a hydraulically wide channel, the proportionality factor between energy slope and Froude number is only a function of the DarcyWeisbach friction factor. In essence, for application to openchannel flow, a modified DarcyWeisbach friction factor f, equal to 1/8 of the conventional DarcyWeisbach friction factor f is applicable. The modified DarcyWeisbach equation for openchannel flow is:
Table 12 shows approximate values of DarcyWeisbach friction factor f and corresponding modified friction factor f for selected values of R in the turbulent range.
QUESTIONS
PROBLEMS
REFERENCES
Chow, V. T. 1959. Openchannel Hydraulics. McGraw Hill, New York. Cornish, V. 1907. Progressive waves in rivers. The Geographical Journal. Vol. 29, No. 1, January, 2331. Craya, A. 1952. The criterion for the possibility of roll wave formation. Gravity Waves, Circular 521, 141151, National Institute of Standards and Technology, Gaithersburg, Md. Dooge, J. C. I., W. B. Strupczewski, and J. J. Napiorkoswki. 1982. Hydrodynamic derivation of storage parameters in the Muskingum model. Journal of Hydrology, 54, 371387. Fread, D. 1993. "Flow Routing," Chapter 10 in Handbook of Hydrology, D. R. Maidment, editor, McGrawHill, New York. Hayami, I. 1951. On the propagation of flood waves. Bulletin, Disaster Prevention Research Institute, No. 1, December. Lagrange, J. L. de. 1788. Mécanique analytique, Paris, part 2, section II, article 2, p 192. Lighthill, M. J., and G. B. Whitham. 1955. On kinematic waves: I. Flood movement in long rivers. Proceedings, Royal Society of London, Series A, 229, 281316. Ponce, V, M., and D. B. Simons. 1977. Shallow wave propagation in openchannel flow. Journal of the Hydraulics Division, ASCE, Vol. 103, No. HY12, December, 14611476. Ponce, V. M. 1989. Engineering Hydrology: Principles and Practices. PrenticeHall, Englewood Cliffs, New Jersey. Ponce, V. M. 1991a. The kinematic wave controversy. Journal of Hydraulic Engineering, ASCE, Vol. 117, No. 4, April, 511525. Ponce, V. M. 1991b. New perspective on the Vedernikov number. Water Resources Research, Vol. 27, No. 7, 17771779, July. Ponce, V. M., and P. J. Porras. 1995. Effect of crosssectional shape on freesurface instability. Journal of Hydraulic Engineering, ASCE, Vol. 121, No. 4, April, 376380. Powell, R. W. 1948. Vedernikov's criterion for ultrarapid flow. Transactions, American Geophysical Union, Vol. 29, No. 6, 882886. Seddon, J. A. 1900. River hydraulics. Transactions, ASCE, Vol. XLIII, 179243, June. Vedernikov, V. V. 1945. Conditions at the front of a translation wave disturbing a steady motion of a real fluid, Dokl. Akad. Nauk SSSR, 48(4), 239242. Vedernikov, V. V. 1946. Characteristic features of a liquid flow in an open channel, Dokl. Akad. Nauk SSSR, 52(3), 207210.

Documents in Portable Document Format (PDF) require Adobe Acrobat Reader 5.0 or higher to view; download Adobe Acrobat Reader. 