3.1 ENERGY PRINCIPLE
Under steady flow, mass and energy are conserved in open-channel flow.
On the other hand, under unsteady flow, mass and momentum are conserved.
Energy is expressed in FL units
With reference to the channel of large slope of Fig. 3-1, the total head H at a section 0 containing point A on a streamline may be written as follows:
in which zA = elevation of point A above a datum plane, dA = depth of point A below the water surface, θ = slope angle of the channel bottom, and VA2/(2g) = velocity head of the flow in the streamline passing through point A.
Due to cross-sectional nonuniformity, the velocity head is likely to vary along the flow depth and width. In practice, at a given cross section, the velocity head is based on the mean velocity V for that cross section, and the Coriolis coefficient (Section 2.2) is used to account for the nonuniformity in the velocity distribution. This leads to the total energy for the channel section:
For a channel of small slope, cos θ ≅ 1. Thus, the total energy reduces to:
The statement of conservation of energy between cross sections 1 and 2 leads to (Fig. 3-1):
in which hf = energy loss.
The line representing the total head is the energy grade line (EGL) or energy line. The slope of the line is the energy gradient, or Sf. The slope of the water surface is denoted by Sw and the slope of the channel bottom is So, where So = tan θ. Under uniform flow conditions, all three slopes are the same: Sf = Sw = So.
For a channel of small slope, Eq. 3-4 reduces to:
When α1 = α2 ≅ 1 (a hydraulically wide channel), and hf = 0, Eq. 3-5 reduces to the statement of conservation of energy, or Bernoulli equation:
3.2 SPECIFIC ENERGY
Equation 3-6 is applicable to a channel of small slope. In the limit, when the channel is horizontal, or else, as an approximation, when the channel is sufficiently short, Eq. 3-6 simplifies to a statement of conservation of specific energy. In effect, for z1 ≅ z2, Eq. 3-6 reduces to:
Equation 3-7 states the constancy of specific energy in a hydraulically wide horizontal channel. The specific energy is:
where E is the energy per (unit of) weight, measured with respect to the channel bottom, defined in terms of the flow depth y and the velocity head V 2/(2g).
Since V = Q /A, the specific energy in terms of discharge is:
Since Q is a constant under steady flow, and flow area A is always a unique function of flow depth y, the specific energy for a given Q is only a function of y.
For a given cross section and discharge Q, a plot of flow depth y versus specific energy E leads to the specific energy curve. A typical curve is shown in Fig. 3-2. For a given Q, the specific energy curve is a type of hyperbolic function, featuring a minimum value of E and two limbs, a lower limb AC and an upper limb BC.
The minimum value of E (at point C) characterizes the critical state of flow, with critical depth. The lower limb AC, with smaller depths, describes supercritical flow. This limb is asymptotic to the horizontal axis. The upper limb BC, with larger depths, describes subcritical flow. This limb is asymptotic to a 45° line OD that originates at the origin (O). At any point P on the curve, the abscissa represents the specific energy E and the ordinate the flow depth y.
For a given specific energy, there are two points on the curve, one corresponding to low stage (small depth) and another corresponding to high stage (large depth). These depths are referred to as alternate depths. At point C, the alternate depths coalesce into one depth, termed the critical depth; at that point, the specific energy is a minimum. When the depth of flow is greater than the critical depth, the flow is termed subcritical. When the depth of flow is smaller than the critical depth, the flow is termed supercritical. For a discharge Q, the specific energy curve is curve AB in Fig. 3-2. For a smaller discharge, it is A'B'; for a larger discharge, it is A"B".
Critical flow criterion
From Fig. 3-2, critical depth (critical flow) corresponds to minimum specific energy. To prove this assertion mathematically, differentiate Eq. 3-9 with respect to y and equate to zero, to yield:
By definition, the differential change in flow area dA is equal to the top width T times the differential change in flow depth dy (Fig. 3-2):
Thus, Eq. 3-10 reduces to:
Since Dc = Ac /Tc , and Vc = Q /Ac , Eq. 3-12 reduces to:
The left-hand side of Eq. 3-13 is effectively the square of the Froude number (Section 1-3). Thus, the condition F 2 = 1, or better yet, F = 1, properly describes the condition of critical flow, for which the specific energy is a minimum. Alternatively, Eq. 3-13 may be expressed as follows:
which states that, under critical flow, the velocity head is equal to one-half of the hydraulic depth. Figure 3-3 shows a channel operating at near critical flow.
Equation 3-14 describes the critical flow condition for a channel of small slope that is hydraulically wide (α = 1). In general, for a channel of large slope and arbitrary cross-sectional shape, the critical flow condition is:
in which Dc is the hydraulic depth corresponding to critical flow, measured normal to the channel bottom. Thus, the general critical flow condition is:
and the general definition for Froude number, applicable to any channel, regardless of slope and cross-sectional shape, is:
3.3 LOCAL PHENOMENA
Changes from supercritical to subcritical flow, or from subcritical to supercritical flow, occur frequently in open-channel flow, depending on the prevailing channel slope and cross section. If the changes occur within a relatively short distance, they constitutes local phenomena. The hydraulic drop and the hydraulic jump are examples of local phenomena.
The hydraulic drop
The hydraulic drop is triggered by a sharp depression in the water surface, usually caused by an abrupt change in bottom elevation. The free overfall shown in Fig. 3-4 is an example of a hydraulic drop.
Flow in the proximity of a free overfall is usually rapidly varied; therefore,
the brink depth is somewhat smaller than the critical depth
computed by the theory of parallel flow. Figure 3-5 shows a schematic of the flow near the brink.
The actual drawdown curve is shown with a solid line,
while the theoretical water-surface curve, assuming parallel flow,
is shown with a dashed line.
For channels of small slope, the computed critical depth is about 1.4 times the brink depth;
The hydraulic jump
The hydraulic jump is triggered by an abrupt rise in the water surface as the flow progresses downstream. In a hydraulic jump the flow changes from supercritical upstream to subcritical downstream, accompanied by an appreciable loss of energy. The amount of energy loss depends on the upstream and downstream flow conditions.
The jump occurs frequently under one of the following situations:
The flow depth before the jump is called the initial depth y1, while the flow after the jump is called the sequent depth y2. The initial and sequent depths y1 and y2 are shown on the specific energy curve of Fig. 3-7. They should be distinguished from the alternate depths y1 and y2', which are two possible depths for the same specific energy.
The specific energy E1 at the initial depth y1 is greater than the specific energy E2 at the sequent depth y2. The energy loss due to the hydraulic jump is the difference between specific energies for initial and sequent depths:
Also shown, to the right of Fig. 3-7, is the specific force curve (Section 3.5), where the sequent depths y1 and y2 have the same specific force: F1 = F2.
3.4 MOMENTUM PRINCIPLE
Momentum M is equal to a force integrated over a period of time, or mass times velocity, Eq. 2-27, repeated here for convenience:
The momentum flux, or force F, of a flow with velocity V through a cross section of area A, Eq. 2-30, repeated here as Eq. 3-20, is:
Since Q = VA, the momentum flux, or force F, exerted by a flow of discharge Q and velocity V is:
According to Newton's second law of motion, the change ΔF in momentum flux through a control volume is equal to the resultant of all the external forces acting on the control volume. The external forces are: (1) body force, and (2) surface forces. The body force is the gravitational force, resolved along the direction of motion (the force labeled W sinθ in Fig. 3-8). In general, there is a nonzero channel bottom slope θ; otherwise, the channel bottom would be horizontal and the gravitational component along the direction of motion would vanish.
The surface forces on the control volume are of three kinds: (1) on the bottom, (2) on the sides, and (3) on the top. The bottom surface force is due to friction, which is always acting in the direction opposite to the flow (the force labeled Ff in Fig. 3-8). There is no such thing as zero friction; under certain conditions, friction may be neglected, but it is never zero.
The side surface forces are two:
one upstream, the force labeled P1 in Fig. 3-8,
and another downstream, the force labeled P2. These forces are
due to the water pressure, which is
either hydrostatic under parallel flow, or nonhydrostatic under convex or concave
The top surface force is due to wind. Under the scales normally considered in open-channel flow, wind forces are small and are usually neglected. However, wind forces may not be negligible in cases on free-surface flow in reservoirs or flow in a wide open space such as the ocean.
The statement of momentum (flux) conservation in a control volume is (Fig. 3-8):
Or, in terms of unit weight:
Following the usual convention of mechanics, the momentum flux difference is equal to the flux at the downstream section 2 minus the flux at the upstream section 1. The forces acting on the control volume are positive in the flow direction and negative against it. Equation 3-23 is known as the momentum flux balance equation or, for short, the momentum equation.
For parallel flow in a rectangular channel of small slope and width b, the force P1 is:
Similarly, the force P2 is:
The weight W of the control volume (Fig. 3-8) is:
The weight of the control volume, resolved along the direction of motion (Fig. 3-8), is:
Generally, under typical gradually varied flow conditions, the friction force Ff along the channel bottom is about equal and opposite in sign to the gravitational force W sin θ (Eq. 3-27). Thus, the friction force may be expressed as follows:
in which hf' = head loss due to friction.
The discharge Q is:
Substituting Eqs. 3-24 to 3-29 into Eq. 3-23 and simplifying terms, the following equation is obtained:
Equation 3-30 differs from Eq. 3-5 in several important respects:
Thus, while Eqs. 3-5 and 3-30 look similar, they are not equivalent.
Equation 3-5 applies to steady gradually varied flow, while
Eq. 3-30 applies for unsteady gradually varied flow.
In other words,
In practice, the momentum principle applies to problems where forces can be shown to play a significant role. Typically, problems of steady gradually varied flow use conservation of energy, while problems of unsteady gradually varied flow use conservation of momentum. The exception is the hydraulic jump, which is rapidly varied. Both energy and momentum principles are used in the solution of the hydraulic jump.
3.5 SPECIFIC FORCE
In horizontal channels, the gravitational force resolved along the direction of motion is effectively zero. As a convenient approximation, for nearly horizontal channels, the gravitational force may be considered small and be neglected on practical grounds.
The frictional force develops along the channel bottom; the longer the channel, the greater the frictional force. Thus, for a short channel, the frictional force may be taken as sufficiently small and neglected on practical grounds. Note that the frictional force is never zero; its neglect is only justified as an approximation, when compared with the other forces that are present in open-channel flow.
The neglect or elimination of the gravitational and frictional forces in the momentum flux balance reduces it to:
Assuming β1 = β2 = 1, Eq. 3-31 reduces to:
Since V1 = Q / A1, and V2 = Q / A2:
The pressure force P acting on a cross-section of area A and distance z̄ from the centroid of the area to the water surface [Fig. 3-9 (b)] is:
Substituting Eqs. 3-35 and 3-36 into Eq. 3-33, and dividing by unit weight γ :
In general, the specific force is defined as follows:
Equation 3-37 states that specific force is conserved in open-channel flow in a short horizontal channel, i.e., F1 = F2. Note that specific force is a force per unit of γ, the weight per unit of volume; therefore, specific force has the units of volume [L3].
The specific force curve shown in Fig. 3-9 (c) is obtained by plotting F in the abscissas and y in the ordinates. This curve is similar to the specific energy curve [Fig. 3-9 (a)], but with significant differences. The limb AC approaches the horizontal axis asymptotically toward the right. The limb BC rises upward and extends without limit to the right.
For a given value of specific force F1, the curve has two possible depths: y1 and y2. These are the initial and sequent depths of a hydraulic jump, respectively. At point C [Fig. 3-9 (c)], the sequent depths coalesce into one depth, termed the critical depth; at that point, the specific force is a minimum. Note that this is the same critical depth obtained by specific energy considerations; see Fig. 3-9 (a).
Critical flow criterion for specific force
As in the case of specific energy, to prove that minimum specific force corresponds to the critical flow criterion, differentiate Eq. 3-38 with respect to y to yield:
For a change in depth dy, the corresponding change d(z̄A) in the static moment of the water area about the free surface is equal to:
As usual in differential calculus, the second-order term in Eq. 3-40 is neglected, to yield:
Therefore, Eq. 3-39 reduces to:
Simplifying Eq. 3-42:
Since dA /dy = T, Q /A = V, and A /T = D, Eq. 3-43 reduces to:
which is the square of the Froude number:
Equation 3-45 is the criterion for critical flow, applicable to both specific energy and specific force (specific momentum) curves.
Note that the sequent depth y2 is always smaller than the high alternate depth y2' (Fig. 3-7). Furthermore, the energy E2 is always smaller than the energy E1, while the specific force F2 remains equal to the specific force F1 [Fig. 3-7 and Fig. 3-9 (c)]. In order to maintain a constant specific force, the flow depth must increase from y1 to y2 at the cost of losing a certain amount of energy. The energy loss is equal to: ΔE = E1 - E2. This situation occurs in the hydraulic jump, where the specific forces before and after the jump are equal, but with a consequent loss of energy (Fig. 3-10).
Specific force per unit of channel width
In a hydraulically wide channel, the specific force per unit of channel width b is:
where q = Q /b.
In terms of mean velocity V = q/y, the specific force per unit of channel width b is:
Specific force in units of force
In units of force, the specific force is:
In units of force, the specific force per unit of channel width b is:
In units of force, and in terms of mean velocity V, the specific force per unit of channel width b is:
3.6 MOMENTUM PRINCIPLE IN NONPRISMATIC CHANNELS
When the external forces are negligible or known in advance, the momentum principle may be applied to nonprismatic channels. The momentum principle is particularly applicable to the hydraulic jump, where the high internal losses that occur cannot be fully evaluated using the energy principle alone.
The following example illustrates how the momentum principle may be applied to the design of a channel transition (Chow, 1959). In this case, the upstream flow is supercritical and the downstream flow subcritical, opening up the possibility that a hydraulic jump may form somewhere within the transition. The latter may be avoided by modifying the shape of the channel cross section in the transition; for example, raising the bed elevation as shown in part (b) below.
Example 3-2: Design of a channel transition, from supercritical to subcritical flow
Assume that a rectangular channel 8-ft wide carries 100 cfs at a depth of 0.5 ft.
The channel is connected through a 50-ft long straight-line transition
to a 10-ft wide channel flowing at a depth of 4 ft.
It is desired to design the transition, including the calculation of the flow profile. For simplicity, neglect friction and other minor losses
Use Eq. 3-9 to calculate
the total energy at the upstream end of the transition:
The difference between these two values is the energy loss, or energy drop: ΔE = Eu - Ed = 6.117 ft. This energy must be dissipated through the transition by some means. The following two cases are considered:
A. Abrupt energy dissipation.
The upstream and downstream Froude numbers (Eq. 1-9) are
Fu = 6.23 and
Table 3-1 shows the computations of the channel transition, with hydraulic jump. Note the following:
Table 3-1 shows that the specific force for the low stages (Col. 5) varies little
along the transition, having a mean value of about 78.8 ft3.
On the other hand, the specific force for the high stages
By varying the shape (width) of the cross sections in the transition, the location of the intersection of the specific force lines may be changed, i.e., the position of the jump. Changing the flow depth in the downstream channel (i.e., at the end of the transition) will also change the position of the jump. Generally, a small increase in the downstream depth will move the jump upstream; conversely, a small decrease in the downstream depth will move the jump downstream.
B. Gradual energy dissipation. The hydraulic jump may be avoided by providing a way to gradually dissipate the energy. This can be accomplished by varying the width or raising the bottom of the transition in such a way that the actual energy grade line (EGL) is a straight line joining the total heads at the two end sections [Fig. 3-12 (b)]. In an actual design, first assume the flow profile and then proportion the dimensions of the raised bottom following applicable hydraulic laws. The following steps are recommended:
Table 3-2 shows the computations of the channel transition, without hydraulic jump. Note the following:
Critical flow over the hump occurs at a distance, measured from the upstream end,
where the flow depth (Col. 9)
is equal to the critical depth (Col. 10). For this example,
the match occurs at a distance X = 45.26 ft, and the flow depth
at the critical flow section is: y = yc = 1.478 ft. [This result was obtained
using ONLINE TRANSITION DESIGN
SUPERCRITICAL B with a resolution of
Figure 3-12 (b) shows the shape of the hump necessary
to provide a gradual drop in the energy gradeline (EGL),
thus avoiding the hydraulic jump.
Note that the specific force lines intersect in the vicinity of the critical flow section
Fig. 3-12 Energy and momentum principles applied to a channel transition:
Figure 3-12 (b) shows the shape of the hump necessary to provide a gradual drop in the energy gradeline (EGL), thus avoiding the hydraulic jump. Note that the specific force lines intersect in the vicinity of the critical flow section (Table 3-2).
Fig. 3-12 Energy and momentum principles applied to a channel transition:
Chow, V. T. 1959. Open-channel Hydraulics. McGraw Hill, New York.
Henderson, F. M. 1966. Open channel flow. Macmillan, New York.
|Documents in Portable Document Format (PDF) require Adobe Acrobat Reader 5.0 or higher to view; download Adobe Acrobat Reader.|