Computer Applications in Hydraulic Engineering Ebook
Computer Applications in Hydraulic Engineering Ebook
)ZESBVMJD&OHJOFFSJOH
9th Edition
Brian E. Whitman
Thomas M. Walski
Thomas E. Barnard
S.Rocky Durrans
Steve Lowry
Michael E. Meadows
$PNQVUFS"QQMJDBUJPOTJO)ZESBVMJD&OHJOFFSJOH
9th Edition
Bentley, the “B” Bentley logo, and Bentley Institute Press are either registered or
unregistered trademarks or service marks of Bentley Systems, Incorporated or one of its
direct or indirect wholly owned subsidiaries. Other brands and product names are trademarks
of their respective owners.
The publisher does not warrant or guarantee any of the products described herein or perform
any independent analysis in connection with any of the product information contained herein.
The publisher also does not assume, and expressly disclaims, any obligation to obtain and
include information other than that provided to it by the manufacturer.
The reader is expressly warned to consider and adopt all safety precautions that might
be indicated by the activities herein and to avoid all potential hazards. By following the
instructions contained herein, the reader willingly assumes all risks in connection with such
instructions.
The publisher makes no representation or warranties of any kind, including but not limited
to, the warranties of fitness for particular purpose of merchantability, nor are any such
representations implied with respect to the material set forth herein, and the publisher takes
no responsibility with respect to such material. The publisher shall not be liable for any
special, consequential, or exemplary damages resulting, in whole or part, from the readers’
use of, or reliance upon, this material.
ISBN - 978-1-934493-70-0
ISBN (ePub) - 978-1-934493-71-7
LCCN – 2021930961
Published by:
Bentley Institute Press
Bentley Systems, Incorporated
685 Stockton Drive
Exton, PA 19341, United States
www.bentley.com
Printed in the United States
Authors and Editors
instructions on how to download the software licenses for each program con-
tained on the academic CD.
In the Eighth Edition, we added chapter in transient analysis with tutori-
als using the HAMMER software (Chapter 7). Tutorials throughout the text
have been updated and expanded per the latest software versions, and current
updates have been incorporated into the software on the CD-ROM. CD-ROM
has been replaced by a procedure to download software from Bentley website.
In this Ninth Edition, we added material about Low Impact Development
(LID) for stormwater management to chapter 9, and the textbook CD-ROM
has been replaced by a procedure to download the latest software from Bentley
website. We took all the tutorials out of the print book and placed them on-line
so they can be easily updated as the software changes and is improved with new
features. We will also placed additional problems to go with all the chapters
on this same on-line location. We would like to thank Brandon Whitman who
completed, reviewed, and provided improvements for all the tutorials that were
updated to the latest Bentley CONNECT editions.
The following table lists the software that accompanies each chapter.
Chapter Accompanying Software
Chapter 1 FlowMaster
Chapter 2 not required
Chapter 3 StormCAD or SewerGEMS
Chapter 4 CulvertMaster
Chapter 5 PondPack
Chapter 6 WaterGEMS
Chapter 7 HAMMER
Chapter 8 SewerCAD or SewerGEMS
Chapter 9 SewerGEMS
We welcome feedback regarding our products and services; we feel this is the
best way for us to continue providing software, training, publications, and
support that are tailored to the needs of professionals in the field of hydraulics
and hydrology today.
Included with this book is access to download and install academic versions
of our award-winning software, WaterGEMS, HAMMER, SewerGEMS,
SewerCAD, StormCAD, PondPack, CulvertMaster, and FlowMaster. They
have been included to assist you with working the problems contained in
this book and to introduce you to the latest developments in computer-aided
hydraulic and hydrologic modeling. Details on accessing the software is
included in Appendix A.
We hope that once you try the companion software packaged with CAiHE,
you’ll discover many ways in which it can enhance your work outside of this
book. You can get more products and upgrade information on any of our
software by calling our sales department or by visiting the Bentley’s OpenFlows
website:
Sales
https://www.bentley.com/en/about-us/contact-us/
sales-contact-requestorcalling
OpenFlows
https://www.bentley.com/en/products/brands/openflows
North America
1-800-BENTLEY (1 800 236 8539)
+1 610 458 5000
Asia
400 842 7516
+86 10 5929 7000
Bentley also works to bring the global civil engineering community together
with our online forums. Get in touch with fellow professionals worldwide and
discuss the latest tips and tricks in hydraulic and hydrologic modeling. Post
your engineering questions and share your unique modeling experiences with
an audience of thousands of professionals. Participants from around the world
include many industry-recognized experts and users like you. Get connected
today by visiting:
http://communities.bentley.com
https://communities.bentley.com/products/hydraulics___hydrology/
If you have any comments regarding this publication and/or any of Bentley
products and services, please contact us at:
Bentley Systems, Inc. has developed hydrology and hydraulics software for civil
engineers and provides textbooks, training, and technical support to tens of
thousands of professional civil engineers, modelers, and universities.
Along the way, we have learned a lot about our clients and their professional
and educational backgrounds. Our continuing education program offers
courses for professionals who need to quickly get up-to-speed with various
numerical methods and practices.
Why is this important?
Our experience has shown us a great deal about the areas where engineers are
being trained, and it has also revealed a significant gap in this training—the
link between hydraulic theory and practical computer applications.
Hundreds of textbooks exist that offer enormous detail in the areas of engi-
neering history, equation derivations, and hand calculation methods. There
are also hundreds of published theses and articles that deal with computer
applications—unfortunately, most of these are highly research oriented and
are usually tied to a specific case study or an unusual set of circumstances. Both
of these publication types are very important to the civil engineering industry,
but in the majority of cases, they fail to address the issues most commonly
encountered by professionals in day-to-day practice.
This publication is intended as an introduction to the more common applica-
tions of water resources engineering software, and it demonstrates the types of
situations that an engineer will most likely come across on a daily basis in the
real world. It shows the true benefits of computer software: increased efficiency,
better flexibility, and—most important—an increased ability to try different
and better designs.
It is our hope that engineers, technicians, and students will find this book to
be challenging, but also easy to understand and very practical. Combined with
standard hydraulic references, we believe that this text provides many of the
tools needed to successfully proceed with a career in the fields of hydrology
and hydraulics.
Bibliography 261
Bentley Systems Software 267
Scenario Management 275
Index 289
System of Units
The basis for the problems in this book is real-world systems designed and con-
structed using U.S. customary units. Those units are retained in the problems.
In response to the goal to convert to a globally consistent system of units, many
agencies in the United States changed their regulations to require engineering
calculations and measurements in the International System of Units (SI). Even
though a new design may require SI units, likely the existing engineering sys-
tem, of which the new design is to become a part, was designed and constructed
using U.S. customary units. Consequently, proper system evaluation will require
the engineer to use both systems of units. This requirement will persist for many
years. Therefore it is essential that students, engineers, regulators, and planners
recognize and use contemporary units and, often, a dual system of units.
Conduit wall
Depth
Velocity distribution
Flow area
Wetted perimeter
The variation of flow velocity within a cross section complicates the hydraulic
analysis, so the engineer usually simplifies the situation by looking at the aver-
age (mean) velocity of the section for analysis purposes. This average velocity is
defined as the total flow rate divided by the cross-sectional area and is in units
of length per time.
Q
V =
A
where
V = average velocity (m/s, ft/s)
Q = flow rate (m3/s, ft3/s)
A = area (m2, ft2).
The hydraulic radius can often be related directly to the geometric properties
of the channel. For example, the hydraulic radius of a full circular pipe (such as
a pressure pipe) can be directly computed as:
A
R=
Pw
or
π ⋅ D2 /4 D
R circular = =
π⋅D 4
where
R = hydraulic radius (m, ft)
A = cross-sectional area (m2, ft2)
Pw = wetted perimeter (m, ft)
D = inside pipe diameter (m, ft).
Flow Conveyance
Water travels downhill from points of higher energy to points of lower energy
(unless forced to do otherwise) until it reaches a point of equilibrium, such as
an ocean. This tendency is facilitated by the presence of natural conveyance
channels such as brooks, streams, and rivers. The water’s journey may also be
aided by man-made structures such as drainage swales, pipes, culverts, and
canals. Hydraulic concepts can be applied equally to both man-made struc-
tures and natural features.
If the flow in a conveyance section is open to the atmosphere, such as in a
culvert flowing partially full or in a river, it is said to be open-channel flow or
free-surface flow. If a channel is flowing completely full, as with a water dis-
tribution pipe, it is said to be operating under full-flow conditions. Pressure
flow is a special type of full flow in which forces on the fluid cause it to push
against the top of the channel as well as the bottom and sides. These forces may
result from, for example, the weight of a column of water in a backed up sewer
manhole or elevated storage tank.
Steady Flow
Speaking in terms of “flow,” the word steady indicates that a constant flow
rate is assumed throughout an analysis. In other words, the flow velocity does
not change with respect to time at a given location. For most hydraulic calcu-
lations, this assumption is reasonable. A minimal increase in model accuracy
does not warrant the time and effort that would be required to perform an
analysis with changing (unsteady) flows over time.
When analyzing tributary and river networks, storm sewers, and other collec-
tion systems in which it is desirable to vary the flow rate at different locations
throughout the system, the network can often be broken into segments that
can be analyzed separately under steady flow conditions.
Uniform Flow
Uniform flow refers to the hydraulic condition in which the discharge and
cross-sectional area (and therefore velocity) are constant throughout the length
of the channel. For a pipe flowing full, the only required assumptions are that
the pipe be straight and have no contractions or expansions. For an open chan-
nel, additional assumptions include the following:
• The depth of flow must be constant [i.e., the hydraulic grade line (HGL)
must be parallel to the channel bed]. This depth of flow is called normal
depth.
• Because the velocity is constant, the velocity head does not change through
the length of the section; therefore the energy grade line (EGL) is parallel
to both the HGL and the channel bed.
In channels that are prismatic, the flow conditions will typically approach uni-
form flow if the channel is sufficiently long. When this occurs, the net force
on the fluid approaches zero because the gravitational force is equal to the
opposing friction forces from the channel bottom and walls. Nonuniform flow
is the hydraulic condition in which the discharge and cross-sectional area (and
therefore velocity) are not constant throughout the length of the channel. In
open channels, the velocity profile and cross-sectional area of flow will gradually
change. The analysis of gradually varied flow is described in Chapters 3.6 and 9.4
The velocity at any given point within the turbulent section will be closer to
the mean velocity of the entire section than with laminar flow conditions.
Turbulent flow velocities are closer to the mean velocity because of the con-
tinuous mixing of flow, particularly the mixing of low-velocity flow near the
channel walls with the higher velocity flow toward the center.
4VR
Re =
ν
where
Re = Reynolds number (unitless)
V = average velocity (m/s, ft/s)
R = hydraulic radius (m, ft)
ν = kinematic viscosity (m2/s, ft2/s).
If the Reynolds number is below 2,000, the flow is generally laminar. For flow
in closed conduits, if the Reynolds number is above 4,000, the flow is gener-
ally turbulent. Between 2,000 and 4,000, the flow may be either laminar or
turbulent, depending on how insulated the flow is from outside disturbances.
In open channels, laminar flow occurs when the Reynolds number is less than
500 and turbulent flow occurs when it is above 2,000. Between 500 and 2,000,
the flow is transitional.
Solution
From the section’s shape (rectangular), we can easily calculate the area as the
rectangle’s width multiplied by its depth. Note that the depth used should be
the actual depth of flow, not the total height of the cross section. The wetted
perimeter can also be found easily through simple geometry.
From fluids reference tables, we find that the kinematic viscosity for water at
20°C is 1.00 × 106 m2/s. Substituting these values into the formula to compute
the Reynolds number results in
1.2 Energy
The Energy Principle
The first law of thermodynamics states that for any given system, the change
in energy (ΔE) is equal to the difference between the heat transferred to the
system (Q) and the work done by the system on its surroundings (W ) during
a given time interval.
The energy referred to in this principle represents the total energy of the sys-
tem, which is the sum of the potential energy, kinetic energy, and internal
(molecular) forms of energy such as electrical and chemical energy. Although
internal energy may be significant for thermodynamic analyses, it is commonly
neglected in hydraulic analyses because of its relatively small magnitude.
In hydraulic applications, energy values are often converted into units of
energy per unit weight, resulting in units of length. Using these length equiv-
alents gives engineers a better “feel” for the resulting behavior of the system.
When using these length equivalents, the engineer is expressing the energy
of the system in terms of “head.” The energy at any point within a hydraulic
system is often expressed in three parts, as shown in Figure 1.4:
• pressure head p/g
• elevation head z
• velocity head V 2 2 g
where
p = pressure (N/m2, lb/ft2)
γ = specific weight (N/m3, lb/ft3)
V12 HL
2g Energy g
rade line V22
p1 Hydraulic grade line 2g
g p2
g
Z1 Z2
Datum
p1 V2 p V2
+ z1 + 1 + HG = 2 + z2 + 2 + HL
g 2g γ 2g
where
p = pressure (N/m2, lb/ft2)
γ = specific weight of the fluid (N/m3, lb/ft3)
z = elevation above a datum (m, ft)
V = fluid velocity (m/s, ft/s)
g = gravitational acceleration (m/s2, ft/s2)
HG = head gain, such as from a pump (m, ft)
HL = combined head loss (m, ft).
Hydraulic Grade
The hydraulic grade is the sum of the pressure head (p/γ) and elevation head (z).
For open-channel flow (in which the pressure head is zero), the hydraulic grade
elevation is the same as the water surface elevation. For a pressure pipe, the
hydraulic grade represents the height to which a water column would rise in
a piezometer (a tube open to the atmosphere rising from the pipe). When the
hydraulic grade is plotted as a profile along the length of the conveyance sec-
tion, it is referred to as the HGL.
Energy Grade
The energy grade is the sum of the hydraulic grade and the velocity head
(V 2/2g). This grade is the height to which a column of water would rise in a
pitot tube (an apparatus similar to a piezometer, but also accounting for fluid
velocity). When plotted in profile, this parameter is often referred to as the
EGL. For a lake or reservoir in which the velocity is essentially zero, the EGL
is equal to the HGL.
Neglecting the velocity simplifies the energy equation even further, and we can
now solve for head loss as
HL = 540 m − 434 m − (586,000 N/m2)/9,810 N/m3 = 46.27 m
The rate of head loss (or friction slope) can now be computed as
Friction slope = 46.27 m/(2 × 1,000 m) = 0.023 m/m, or 23 m/km
V = kCR x S y
where
V = mean velocity
C = flow resistance factor
R = hydraulic radius
S = friction slope
x, y = exponents
k = factor to account for empirical constants, unit conversion, etc.
The lining material of the flow channel usually determines the flow resistance
or roughness factor, C. However, the ultimate value of C may be a function of
the channel shape, depth, and fluid velocity.
Manning’s Equation
Manning’s equation is the most commonly used open-channel flow equation.
The roughness component, C, is typically assumed to be constant over the full
range of flows and is represented by a Manning’s roughness value, n. These
n-values have been experimentally determined for various materials and should
not be used with fluids other than water. Manning’s equation is:
k 2 / 3 1/ 2
V = R S
n
where
V = mean velocity (m/s, ft/s)
k = 1.49 for U.S. customary units, 1.00 for SI units
n = Manning’s roughness value
V = C RS
where
V = mean velocity (m/s, ft/s)
C = roughness coefficient (see following calculation)
R = hydraulic radius (m, ft)
S = friction slope (m/m, ft/ft).
The roughness coefficient, C, is related to Kutter’s n through Kutter’s equa-
tion. Note that the n-values used in Kutter’s equation are actually the same as
Manning’s n coefficients.
where
C = roughness coefficient
n = Manning’s roughness value
R = hydraulic radius (m, ft)
S = friction slope (m/m, ft/ft).
Hazen–Williams Equation
The Hazen–Williams equation is most frequently used in the design and analysis
of pressure pipe systems. The equation was developed experimentally and there-
fore should not be used for fluids other than water (and only within temperatures
normally experienced in potable water systems). The Hazen–Williams equation is:
where
V = mean velocity (m/s, ft/s)
k = 1.32 for U.S. customary units, or 0.85 for SI units
8g
V = RS
f
where
V = flow velocity (m/s, ft/s)
g = gravitational acceleration (m/s2, ft/s2)
f = Darcy–Weisbach friction factor (unitless)
R = hydraulic radius (m, ft)
S = friction slope (m/m, ft/ft).
The Darcy–Weisbach friction factor, f, can be found using the Colebrook–
White equation for fully developed turbulent flow, as follows:
where
k = roughness height (m, ft)
R = hydraulic radius (m, ft)
Re = Reynolds number (unitless).
This iterative search for the correct value of “f ” can become quite time-consuming
for hand computations and computerized solutions of many pipes. Another
method, developed by Swamee and Jain, solves directly for f in full-flowing cir-
cular pipes. This equation is:
1.325
f = 2
⎡ ⎛ k 5.74 ⎞⎟⎤
⎢ ln ⎜⎜ ⎥
⎢⎣ ⎜⎝ 3.7 D + Re0.9 ⎟⎟⎠⎥
⎦
where
f = friction factor (unitless)
k = roughness height (m, ft)
D = pipe diameter (m, ft)
Re = Reynolds number (unitless).
Discussion
If the same system is analyzed with 2,000–3,000 gal/min of flow, the difference
in head loss between the two equations becomes almost 10 ft.
Why such a big difference? For starters, the two methodologies are completely
unrelated, and the estimated roughness coefficients were taken from a list of
approximate values. If the Hazen–Williams equation is used with a roughness
value of 125, the results are much closer. This difference should emphasize the
fact that models are only as good as the data that is input into them, and the
engineer needs to fully understand all of the assumptions that are being made
before accepting the results.
V2
E= y+
2g
where
E = specific energy (m, ft)
y = depth of flow (m, ft)
V = mean velocity (m/s, ft/s)
g = gravitational acceleration (m/s2, ft/s2).
If we assume the special case of an infinitely short section of open channel
(with essentially no friction losses and no change in elevation), we see that the
general energy equation can be reduced to an equality of specific energies. In
other words,
V12 V2
E1 = y1 + = y2 + 2 = E2
2g 2g
Channel depth, y
45°
Subcritical flow, low velocity, F < 1
Critical depth, F = 1
Supercritical flow, high velocity, F > 1
Specific energy, E
Recall that the velocity of the section is directly related to the area of flow,
and that the area of flow is a function of channel depth. This means that, for a
given discharge, the specific energy at each point is solely a function of channel
depth and more than one depth may exist with the same specific energy. If the
channel depth is plotted against specific energy for a given flow rate, the result
is similar to the graph shown in Figure 1.5.
As this figure shows, a depth exists for which the specific energy is at a mini-
mum. This depth is called the critical depth. If the velocity is higher than crit-
ical velocity (that is, the depth is less than critical depth), the flow is considered
supercritical. If the velocity is lower than critical velocity (the depth is greater
than critical depth), the flow is subcritical.
The velocity at critical depth is equal to the wave celerity—the speed at which
waves will ripple outward from a pebble tossed into the water. A unitless value
called the Froude number, F, represents the ratio of actual fluid velocity to
wave celerity. The Froude number is computed as follows:
V
F =
gD
where
F = Froude number (unitless)
D = hydraulic depth of the channel, defined as A/T
A = flow area (m2, ft2)
T = top width of flow (m, ft)
V = fluid velocity (m/s, ft/s)
g = gravitational acceleration (m/s2, ft/s2).
By definition, when the flow is at critical depth (that is, the velocity is equal to
the wave celerity), the Froude number must be equal to 1. The equation can
therefore be rewritten and refactored to form the following equality:
A3 Q 2
=
T g
where
A = flow area (m2, ft2)
T = top width of flow (m, ft)
Q = channel flow rate (m3/s, ft3/s)
g = gravitational acceleration (m/s2, ft/s2).
This equation can now be used to determine the depth for which this equality
holds true, which is critical depth. For simple geometric shapes, the solution is
relatively easy to determine. However, quite a few iterations may be required to
find the solution for an irregularly shaped channel such as a natural streambed.
In fact, several valid critical depths may exist for irregular channels.
Example 1.4: Uniform Open-Channel Flow
A concrete trapezoidal channel has a bottom width of 4 m and 45-degree side
slopes. If the channel is on a 1% slope and is flowing at a depth of 1 m throughout
its length, how much flow is being carried (use Manning’s equation)? How much
flow would the same channel carry if it were a rectangular channel 4 m wide?
Solution
Because the channel is flowing at the same depth throughout, we can assume
that normal depth has been achieved (that is, the friction slope is equal to the
channel slope). We will assume a Manning’s n of 0.013 for concrete.
From the trapezoidal geometry, we can easily calculate the area and wetted
perimeter, and then the hydraulic radius, as follows:
A = (4 m × 1 m) + 2 × (0.5 × 1 m × 1 m) = 5.00 m2
Pw = 4 m + 2 × (1 m × 20.5) = 6.83 m
R = A/Pw = 5.00 m2/6.83 m = 0.73 m
Manning’s equation for velocity can then be solved. The discharge can be com-
puted as
V = (1.00/0.013) × 0.732/3 × 0.011/2 = 6.25 m/s
Q = V × A = 6.25 m/s × 5.00 m2 = 31.2 m3/s
To answer the second part of the question, we simply repeat the steps for a
rectangular section shape.
A = (4 m × 1 m) = 4 m2
Pw = 4 m + 2 × (1 m) = 6 m
R = 4 m2/6 m = 0.67 m
V = (1.00/0.013) × 0.672/3 × 0.011/2 = 5.87 m/s
Q = 5.87 m/s × 4 m2 = 23.5 m3/s
As we would expect, this discharge is less than the discharge of the trapezoidal
section.
Example 1.5: Critical Depth
What is the critical depth for a grassy triangular channel with 2H:1V side
slopes and a 0.5% slope when the flow is 3.00 m3/s? If the channel is actually
flowing at a depth of 1.2 m, is the flow critical, subcritical, or supercritical?
Solution
For such a simple geometry, we can quickly create a relationship between the
flow area, top width, and depth of flow:
T = 4y, A = 0.5 × T × y = 0.5 × 4y × y = 2y2 m2
Inserting these values into the previous equation for critical depth, we can
algebraically solve for the channel depth:
(2y2)3/4y = Q2/g
8y6/4y = Q2/g
2y5 = Q2/g
y5 = Q2/2g
y = (Q2/2g)1/5 = [(3.00m3/s)2/(2 × 9.8 m/s2)]0.2 = (0.46 m5)0.2 = 0.86 m
The critical depth for this section is 0.86 m. The actual flow depth of 1.2 m is
greater than critical depth, so the flow is subcritical.
Orifices
Orifices are regularly shaped, submerged openings through which flow is pro-
pelled by the difference in energy between the upstream and downstream sides
of the opening. The stream of flow expelled from the orifice is called the jet.
When the jet exits the orifice, adverse velocity components cause it to contract
to a point after which the flow area remains relatively constant and the flow
lines become parallel (see Figure 1.6). This point is called the vena contracta.
Orifices and the orifice equations have the following applications:
• regulating the flow out of detention ponds,
• regulating the flow through channels in the form of radial and sluice gates,
Vena contracta
Jet
p1 V2 p V2
+ z1 + 1 + HG = 2 + z2 + 2 + HL
g 2g γ 2g
P
g =H
Datum 1 2
Taking these known variables and assumptions into account and solving for
V2, the energy equation becomes:
V2 = 2 g (H − HL )
To find the flow exiting the structure at point 2, multiply both sides of the
equation by the orifice area, A.
AV2 = Q = A 2 g (H − HL )
Q = CA 2 gH
Weirs
Weirs are notches or gaps over which fluid flows. The lowest point of structure
surface or edge over which water flows is called the crest, whereas the stream of
water that exits over the weir is called the nappe (see Figure 1.8). Depending
on the weir design, flow may contract as it exits over the top of the weir, and,
as with orifices, the point of maximum contraction is called the vena contracta.
This contraction can be counteracted or suppressed by designing the weir such
that its shape conforms to the shape of the channel. This type of weir is called
a suppressed weir. With a contracted weir, the crest and nappe vary from the
channel to such a degree that a significant contraction of flow area does occur.
In addition to suppressed and contracted weir types, weirs are also distin-
guished as either sharp-crested or broad-crested. A sharp-crested weir has a
sharp upstream edge formed so that the nappe flows clear of the crest. Broad-
crested weirs have crests that extend horizontally in the direction of flow far
enough to support the nappe and fully develop hydrostatic pressures for at least
a short distance.
Weirs can also be distinguished by their shapes. The most common shapes are
shown in Figure 1.9. The effects of weir shape and other factors previously
mentioned are accounted for with modifications to the weir equation (derived
in Example 1.7), such as adjustments the weir coefficient. Table 1.3 contains
information on coefficients for V-notch weirs.
Weirs have the following applications:
• serving as emergency spillways for regulating high-return event flows over-
topping dams and detention ponds,
• regulating the flow in channels,
• measuring flow,
• approximating the flow over roadways acting as broad-crested weirs when
flow exceeds a culvert’s capacity (see Chapter 4),
• approximating the interception capacity of unsubmerged drainage inlets in
swales (see Chapter 3),
• approximating the flow allowed through an unsubmerged culvert operating
under inlet control (see Chapter 4).
Nappe
Crest
C varies between
q H 3/2 0.611 and 0.570
V-notch Q = C 8 2g tanq H
15 depending on
2
H and Q*
L Metric Metric
Q = CLH3/2 C = 1.86
Cipolletti H 4:1
English English
Q = CLH 3/2 C = 3.367
Cd is a function
Hr
crested
p1 V2 p V2
+ z1 + 1 + HG = 2 + z2 + 2 + HL
γ 2g γ 2g
V2 = 2 g (H − HL )
To find the flow, multiply both sides by the flow area, which in this case is the
length of the weir, L, multiplied by the height of the head, H.
P
=H
g
Datum 1 2
AV2 = Q = LH 2 g ( H − HL )
32
Q = CLH
cross section are frequently negligible, and a more detailed analysis may not be
the most efficient or effective use of time and resources.
There are circumstances when an engineer may be called upon to provide a
more detailed analysis, including unsteady flow computations. For a storm
sewer, the flows may rise and fall over time as a storm builds and subsides. For
water distribution piping, a pressure wave may travel through the system when
a valve is closed abruptly (the same “water-hammer” effect can probably be
heard in your house if you close a faucet quickly).
As an engineer, it is important to understand the purpose of an analysis; other-
wise, appropriate methods and tools to meet that purpose cannot be selected.
1.9 Problems
Solve the following problems using the FlowMaster computer program. Unless
stated otherwise, assume water is at 20°C.
1. The cross section of a rough, rectangular, concrete (k = 0.2 × 10−2 ft)
channel measures 6 × 6 ft. The channel slope is 0.02 ft/ft. Using the
Darcy–Weisbach friction method, determine the maximum allowable
flow rate through the channel to maintain one foot of freeboard (free-
board is the vertical distance from the water surface to the overtopping
level of the channel). For these conditions, find the following characteris-
tics (note that FlowMaster may not directly report all of these):
(a) Flow area
(b) Wetted perimeter
(c) Hydraulic radius
(d) Velocity
(e) Froude number
2. A 450-mm circular concrete (n = 0.013) pipe constructed on a 0.6%
slope carries 0.1 m3/s.
(a) Using Manning’s equation and normal depth assumptions, what are
the depth and velocity of flow?
(b) What would the velocity and depth be if the pipe were constructed
of corrugated metal (n = 0.024) instead of concrete?
3. A trapezoidal channel carries 2.55 m3/s at a depth of 0.52 m. The channel
has a bottom width of 5 m, a slope of 1.00%, and 2H:1V side slopes.
(a) What is the appropriate Manning’s roughness coefficient?
(b) How deep would the water be if the channel carried 5 m3/s?
4. Use Manning’s equation to analyze an existing brick-in-mortar (n = 0.015)
triangular channel with 3H:1V side slopes and a 0.05 longitudinal slope.
The channel is intended to carry 7 cfs during a storm event.
(a) If the maximum depth in the channel is 6 in., is the existing design
acceptable?
7. A paved highway drainage channel has the geometry shown in the follow-
ing figure. The maximum allowable flow depth is 0.75 ft (to prevent the
flow from encroaching on traffic), and the Manning’s n-value is 0.018 for
the type of pavement used.
20 ft 5 ft 10 ft
1 ft
0.5 ft
Not to scale
10. Develop a performance curve for the pipe in Problem 9 that shows the
available flow to the industry with residual pressures ranging from 20 to
80 psi (assume the source can maintain 64 psi regardless of flow rate).
Create similar curves for 10- and 12-in. diameter pipes and compare the
differences in flow.
11. Using the Darcy–Weisbach equation, find the minimum size of circular
corrugated metal storm drain (k = 1.0 × 10−3 ft) that will carry 1.5 cfs
with a maximum depth of 6 in. The drain carries water down a hill 3-ft
high to a pond with a free outfall 75 ft away. What pipe size should be
used? Assume pipes are available in 3-in. increments. What would the
maximum capacity of this pipe be? What would the capacity of the pipe
be when it is flowing full?
12. A channel with the cross section shown in the following figure has a
Manning’s coefficient of 0.040 from stations 0 to 3 and 0.054 from sta-
tions 3 to 8. The flow through the channel is 13 m3/s, and the water
surface is 1.7-m high. Find the following:
(a) Weighted Manning’s coefficient
(b) Slope of the channel
(c) Top width
(d) Wetted perimeter
(e) Flow regime (supercritical or subcritical)
2.00
Elevation, m
1.50
1.00
0.50
0.00
0 1 2 3 4 5 6 7 8
Station, m
(c) How would the capacity of the channel be affected if you were to pave
the center of the channel (n = 0.013) between stations 3 and 5?
14. A rectangular concrete channel with a width of 1 m and a height of 0.5 m
is on a slope of 0.008 m/m. Design a concrete circular channel for which
the depth is half of the diameter and the flow area is the same as that of
the rectangular channel. Which channel is more efficient and by how
much?
15. A weir is placed in a rectangular channel to measure the flow. The dis-
charge from the rectangular channel enters a trapezoidal channel with a
stony bottom. The trapezoidal channel is 0.50-m wide at the base with
2:1 (H:V ) equal side slopes. The weir is a sharp-crested, V-notch weir
with a crest 0.43 m above the channel bottom, a weir coefficient of 0.58,
and a notch angle of 1.57 radians. The height of the water above the weir
is 0.70 m, and the depth of water in the trapezoidal channel is measured
to be 0.40 m.
What is the flow rate? What is the slope of the trapezoidal channel (using
Manning’s formula)? If the discharge is increased until the elevation of
the water surface in the trapezoidal channel reaches 0.61 m, what will the
headwater elevation be at the weir?
16. The outlet structure on a pond is used to regulate the flow out of the
pond for different storm events. An outlet structure must be designed
to discharge 2.20 m3/min when the water surface elevation in the pond
reaches 1.52 m, and 6.29 m3/min when the water surface elevation reaches
2.60 m. The outlet structure will be a circular orifice and a sharp-crested
rectangular weir combination, with the centroid of the orifice at an ele-
vation of 0.90 m and the weir crest at an elevation of 2.50 m. Both will
discharge to free outfall conditions.
Assume an orifice coefficient of 0.6. Find the orifice diameter needed to
supply the correct discharge when the water surface reaches the first spec-
ified elevation. What will the discharge from the orifice be when the water
surface reaches the second specified elevation? Find the width of the weir
needed to supply the extra discharge necessary to meet the requirement.
Use Manning’s formula where necessary.
17. An approximately trapezoidal, clean, natural stream carries the discharge
from a pond down a 0.001 slope. The maximum depth in the channel is
0.5 m. The channel has equal side slopes of 3.0 (H:V) and a bottom width
of 1.0 m. The pond discharges water through a circular orifice into the
channel. The centroid of the orifice is located 1.0 m above the bottom of
the channel.
Assume an orifice coefficient of 0.6. Design the orifice to discharge the
maximum flow rate possible without exceeding the maximum allowed
depth in the channel when the water surface in the pond reaches 4.6 m
above the channel bottom. Use Manning’s formula when necessary.
Hydrology is the study of water on and below the Earth’s surface and in the
atmosphere. The discussion in this chapter is limited to rainfall and the result-
ing surface runoff. Data on surface runoff flow rates and volumes provides
the loads used to design and evaluate storm sewers, culverts, ditches, and
detention ponds. Several of these structures are covered in more detail in other
chapters.
2.1 Rainfall
Rainfall data is fundamental for determining the amount of stormwater gener-
ated by a particular storm event. This section describes the properties of rainfall
used in modeling.
For example, a 5-year return period represents a storm event that is expected
to occur once every 5 years on average. This does not mean that two storm
events of that size will not occur in the same year, nor does it mean that the
next storm event of that size will not occur for another 5 years. It just means
that the average will be once every 5 years (that is, there is a 20% chance of
occurrence in any given year).
Solution
For a storm event with a 20% probability of being equal to or exceeded in a
given year, the recurrence interval is computed as
1/0.20 = 5 years
1/0.02 = 50 years
Intensity–Duration–Frequency Data
For a selected storm duration, a rainfall intensity exists that corresponds to
a given exceedance probability or recurrence interval. A rainfall IDF curve
illustrates the average rainfall intensities corresponding to a particular storm
225
200
10
0
ye
ar
175
s
150
Intensity, mm/h
10
ye
125 ar
s
100 2y
ear
s
75
50
25
0
0 15 30 45 60
Duration, min
recurrence interval for various storm durations (see Figure 2.1). These curves
are the result of the statistical analysis of rainfall data for a particular area.
Given the information on the graph shown in Figure 2.1, you can determine
that the average 1-h rainfall intensity expected to be equaled or exceeded, on
average, once every 100 years is 56.0 mm/h.
Although graphical rainfall curves are acceptable for hand calculations, they are
not well-suited to computer analyses. Data is therefore input into hydrologic
software either as equations or in a tabular format. Creating a rainfall table
from a set of IDF curves is a simple matter of manually picking values from
the curves. For example, Table 2.1 can be created from the data presented in
the IDF curves of Figure 2.1.
Rain depth, in
r
in/h
=5
1.0 in
.i
Fills Average intensity for
Avg
over a given time slice
time 0.2 hr can be computed using
the incemental depth
for that time slice.
0.7
0.6
Incrmental rainfall depth, in
0.5
0.4
0.3
0.2
0.1
0
0.0 0.10.20.3 0.4 0.50.60.7 0.80.9 1.01.11.21.3 1.4 1.5 1.61.71.81.92.02.1 2.22.32.42.5
Time, h
about 5 min. In small urban drainage basins where it is often necessary to use
time increments as small as 1 or 2 min, the data must be extrapolated.
In addition to selecting an appropriate Δt, the engineer must select the total
duration to be used when developing a design storm hyetograph. In many
cases, the storm duration will be specified by the review agency having jurisdic-
tion over the area in which a stormwater conveyance facility will be built; this
approach promotes consistency from one design to another.
Many methods have been proposed for distributing a total rainfall depth
throughout a storm to develop a design storm hyetograph. The NRCS devel-
oped one of the more commonly used methods in the United States (SCS,
1986). With this method, Table 2.2 is used to find fractions of the total
1.0
0.9
II
0.8
0.7
Rain fraction
0.6
III
0.5
0.4
IA I
0.3
0.2
0.1
0.0
0 3 6 9 12 15 18 21 24
Time, h
accumulated rainfall depth for storms with 24-h durations. (Figure 2.4 depicts
Table 2.2 graphically.) The storms are classified into various types, with each
type being recommended for use in a certain U.S. geographical region, as
shown in Figure 2.5. If necessary, interpolation may be used to obtain values
not shown in Table 2.2. Nonlinear interpolation methods are recommended
for this purpose.
Solution
Figure 2.5 illustrates that a Type III storm distribution is a reasonable choice
for Boston. From TP 40, the total depth of a 50-year, 24-h storm in Boston
is estimated to be 6.0 in. Table 2.3 illustrates the calculation of the storm
hyetograph.
The first column of the table is the time, in hours, since the beginning of the
storm, and is tabulated in 1-h increments for the total storm duration of 24 h.
(In actuality, the Δt used in the calculations would be 0.1 h; the 1-h increment
is used here for brevity.) The second column is the fraction of the total storm
depth that has accumulated at each time during the storm. These values are
obtained by interpolation from Table 2.2 for the Type III storm distribution.
The third column contains the cumulative rainfall depths for each time during
the storm and is obtained by multiplying each fraction in the second column
by the total storm depth of 6.0 in. The fourth column contains the incremental
IA
III
III
I
Rainfall
distribution
Type I
Type IA
Type II
Type III
III
III
Figure 2.5: Coverage of NRCS (SCS) rainfall distributions (Adapted from SCS, 1986).
5
Precipition, in
0
0 4 8 12 16 20 24
Time, h
depths of rainfall within each time interval during the storm; these values are
computed as the difference between the current and preceding values in the
third column.
The resulting graph of cumulative precipitation is shown in Figure 2.6, and the
hyetograph is shown in Figure 2.7. The height of each bar on the hyetograph
is the average rainfall intensity during that time interval, and the area of each
bar is the incremental rainfall depth during that time interval. Because the time
increment is 1 h (1-h increment is shown for simplicity; actual Δt is 0.1 h), the
value for the height of the bar (in units of in./h) is equal to the incremental
depth for that time increment (in in).
1.6
1.4
1.2
Intensity, in/h
1.0
0.8
0.6
0.4
0.2
0.0
0 to 1
2 to 3
4 to 5
6 to 7
8 to 9
10 to 11
12 to 13
14 to 15
16 to 17
18 to 19
20 to 21
22 to 23
Time, h
Watershed Area
A watershed is an area that drains to a single point of discharge. The first step
in computing runoff at a point of interest (such as a culvert, inlet, or property
outfall) is to delineate its watershed on a contour map. Because water flows
downhill, delineating a watershed is simply a matter of identifying an outfall
and locating the watershed boundary such that any rain that falls within the
boundary will be directed toward that point of discharge. The delineated area
is then measured using a planimeter, or graphical or computer-aided meth-
ods. Because river and stream systems collect water, a watershed may have any
number of subwatersheds within it. The focus of the analysis and the determi-
nation of how large a watershed must be to be analyzed separately depend on
the scope and purpose of the project at hand. Figure 2.8 shows the collection
channels for two typical natural watersheds.
Discharge
points
(outfalls)
After the watershed is delineated and measured, the modeler can proceed to
the next step of computing rainfall abstractions and runoff depth. To deter-
mine the total runoff volume, runoff depth is multiplied by the watershed area.
Rainfall Abstractions
Rainfall that does not contribute to direct surface runoff may be intercepted
by vegetation, infiltrated into the ground surface, retained as depression stor-
age in puddles and small irregularities in the land surface, or returned to the
atmosphere through transpiration and evaporation. Collectively, these losses of
rainfall are called abstractions. The rainfall that remains after abstractions have
occurred comprises the surface runoff and is called effective precipitation.
The four basic types of abstractions are interception, depression storage, infil-
tration, and evaporation. Evaporation, however, is not typically considered in
modeling stormwater conveyances.
Interception
Interception refers to the capture of rainfall on the leaves and stems of vege-
tation before it reaches the ground surface. Water intercepted by vegetation is
returned to the atmosphere by evaporation during dry-weather periods.
Depression Storage
Excess water begins to pond on the land surface when the rainfall intensity
exceeds the infiltration capacity of a soil during a storm event. The ponded
water fills small depressions and irregularities in the ground surface, and
additional water is held on the surface through the phenomenon of sur-
face tension. The water held in depressions is called depression storage, and
it either evaporates during dry-weather periods or infiltrates into the soil.
Noninfiltrated rainfall that remains after surface depressions and irregularities
have been filled contributes to surface runoff.
Infiltration
When rainfall occurs on a pervious surface, some of the rainwater infiltrates
into the ground. The infiltrated water may contribute to groundwater recharge,
or it may be taken up by the roots of vegetation and subsequently transpired
through the leaves of plants. Infiltrated water also may be evaporated from the
soil during dry-weather periods between storm events, or it may move laterally
through the near-surface soils and reappear as surface water in a stream.
Vr = Dr × A
where
Vr = runoff volume,
A = contributing drainage area, and
Dr = total depth of runoff.
De
Dr
Drainage area
Ds
Dp
Dli
Di
Volume of runnoff, Vr
Solution
Total runoff depth is computed as
Dr = 2.0−0.3−0.7−0.2−0 = 0.8 in
f (t ) = fc + ( f0 − fc )e−k (t −t0 )
where
f (t) = the infiltration rate (in/h or mm/h) at time t (min),
f c = a steady-state infiltration rate that occurs for sufficiently large t,
f 0 = the initial infiltration rate at the time that infiltration begins (that is, at
time t = t0), and
k (min−1) = a decay coefficient.
It can be shown theoretically that the steady-state infiltration rate fc is equal to
the saturated vertical hydraulic conductivity of the soil.
fo − fc
Fmax = fc t +
k
(1 − e−kt )
Time
Figure 2.10: Rainfall intensity can limit infiltration rate during the early
part of the storm.
Solution
From Table 2.4, f0 = 3.5 in/h; fc = 2.6 in/h; and k = 0.02 min−1. Thus f (2 h) =
2.6 in/h + (3.5−2.6 in/h)e−0.02(2−0.5) = 3.47 in/h.
Solution
From Table 2.4, f0 = 3.5 in/h; fc = 2.6 in/h; and k = 0.02 min−1. Thus:
Runoff Coefficient
The concept of the runoff coefficient assumes that the effective rainfall intensity
ie is a fraction of the gross rainfall intensity i, or
ie = Ci
Solution
ie = Ci = 0.70 × 2.0 in/h = 1.4 in/h
Table 2.5: Runoff Coefficients for Use in the Rational Method (Schaake
et al., 1967)
Type of Area or Development C
Types of Development
Urban business 0.70–0.95
Commercial office 0.50–0.70
Residential development
Single-family homes 0.30–0.50
Condominiums 0.40–0.60
Solution
From Table 2.5, the runoff coefficients for the park and residential areas are
estimated to be 0.20 and 0.40, respectively. The composite runoff coefficient
for the entire drainage basin of 45 ac is therefore
F = P−Ia−Pe
where
F = equivalent depth of retention (in., mm),
P = total rainfall depth in storm (in., mm),
Ia = equivalent depth of initial abstractions (in., mm), and
Pe = depth of effective precipitation (in., mm).
An assumption made in the development of the curve number method is
F Pe
=
S P − Ia
( P − I a )2
Pe =
( P − Ia ) + S
for values of P > Ia. Data analyzed by the NRCS indicated that Ia is related
to S and on average supported the use of the relationship Ia = 0.2S. Thus the
equation becomes
( P − 0.2S )2
Pe =
P + 0.8S
where P > 0.2S (Pe = 0 when P ≤ 0.2S). Because the initial abstraction Ia
consists of interception, depression storage, and infiltration prior to the onset
of direct runoff, it may be appropriate in some applications to assume that
Ia = 0.1S or Ia = 0.3S instead of Ia = 0.2S. For example, the relationship
Ia = 0.1S might be appropriate in a heavily urbanized area where there is little
opportunity for initial abstractions to occur. The preceding equation must be
modified when the relationship between Ia and S is assumed to be different
from Ia = 0.2S.
The use of this equation for estimating the depth of effective rainfall during
a storm requires an estimate of the maximum possible retention S. NRCS
conducted research to approximate S for various soil and cover conditions.
In order to provide engineers with tables having a manageable range of coeffi-
cients that varying between 0 and 100, the original values for S were rearranged
using the following simple relationship:
1000
CN =
S + 10
where
CN = runoff curve number and
S = maximum possible retention (in).
Rearranging, S is related to the curve number, CN, as
1000
S= − 10
CN
Practical values of CN range from about 30 to 98, with larger values being
associated with more impervious land surfaces. The NRCS has tabulated curve
numbers as a function of soil type, land use, hydrologic condition of the drain-
age basin, and antecedent moisture condition.
Soil Group
Soils in the United States have been classified by the NRCS into four hydro-
logic groups: Groups A, B, C, and D. Group A soils have high infiltration rates
(low runoff potential), even when they are thoroughly wetted. Typical Group A
soils are well-drained sands and gravels. Group D soils are at the opposite end
of the spectrum, having low infiltration rates (high runoff potential). Typical
Group D soils are clays, shallow soils over nearly impervious material, and soils
with a high water table. Group B and Group C soils are in the midrange of
the spectrum.
In the United States, hydrologic soil group information can be determined
from the NRCS soil survey of the county in which the project is located. Note
that when a drainage basin undergoes urbanization, the hydrologic soil group
may change due to compaction of the soil by heavy construction equipment or
mixing of soils as a consequence of grading operations.
Cover Type
The surface conditions of a drainage area have a significant impact on runoff.
For example, in the exaggerated case of a sandy (Group A) soil completely
paved with asphalt, the soil itself will have no impact on the amount of runoff.
Even for pervious conditions, cover type plays a significant role in the amount
of runoff from a site. For example, a heavily forested area will yield runoff
volumes that differ from those of a lawn or plowed field.
Hydrologic Condition
The hydrologic condition of rangeland, meadow, or pasture is defined to be
“good” if it is lightly grazed and has vegetative cover on more than 75% of the
area. Conversely, a “poor” hydrologic condition corresponds to a heavily grazed
area with vegetation covering less than 50% of the surface.
4.2CN II
CN I =
10 − 0.058CN II
23CN II
CN III =
10 + 0.13CN II
where CNI , CNII , and CNIII = curve numbers for AMC-I, -II, and -III,
respectively.
Table 2.6: Runoff Curve Numbers for Urban Areas (SCS, 1986)1
Example 2.8: Estimating Runoff Depth and Volume Using the NRCS
(SCS) Curve Number Method (Modified from SCS, 1986)
Estimate the curve number, depth of runoff (effective precipitation), and
runoff volume for a 1,000-ac drainage basin if the total depth of precipitation
is 5.0 in. All soils in the basin are in hydrologic soil group C. The proposed land
use is 50% detached houses with 0.25-ac lots; 10% townhouses with 0.125-
ac lots; 25% schools, parking lots, plazas, and streets with curbs and gutters;
and 15% open space, parks, and schoolyards with good grass cover. Use an
antecedent soil moisture condition of AMC-III. The detached housing and
townhouse areas have directly connected impervious area percentages similar
to those assumed in Table 2.6.
Solution
The composite curve number corresponding to AMC-II conditions is com-
puted as a weighted average, as shown in Table 2.7.
Thus the composite CN is computed as
CN = 86,100 ÷ 1,000 = 86
23(86 )
CN III = = 93
10 + 0.13(86 )
S = (1,000/93) − 10 = 0.75 in
Ia = 0.2S = 0.15 in
Time of Concentration
The maximum amount of flow coming from any watershed is related to the
amount of time it takes for the entire watershed to be contributing to the flow.
In other words, it may start raining right now, but it could be several minutes
(or even hours) before the water that falls on some parts of the watershed
actually makes its way to the discharge point.
Some places in a watershed are hydraulically closer to the discharge point
than others, but for peak flow generation only the most hydraulically remote
location is considered crucial. The amount of time that it takes for the first
drop of water falling on this location to work its way to the discharge point is
called the time of concentration.
Many methods for calculating time of concentration are given in various
private, federal, and local publications. Although each of these methods
is different (in some cases only slightly), they are all based on the type of
ground cover, the slope of the land, and the distance along the flow path.
In most locales, there is also a minimum time of concentration (typically
5–10 min) recommended for small watersheds, such as a section of a park-
ing lot draining to a storm sewer. Some methods predict the response time
directly, whereas others predict the velocity of flow. The predicted velocity
coupled with estimates of the flow path length can then be used to estimate
the response time.
With few exceptions, methods for prediction of basin response time are
empirical in nature. Consequently, large errors in response time estimates
can be expected to occur if these methods are not carefully selected and
applied, and these errors can significantly affect peak runoff estimates. The
method selected for estimation of tc should be one that was developed for
basin conditions similar to those in the drainage basin for which an estimate
is desired.
Table 2.8 lists several commonly used methods for estimation of the time of
concentration of a drainage basin. Figure 2.11 illustrates average overland flow
velocities as a function of land use characteristics and surface slope. When flows
Table 2.8: Commonly Used Methods for Estimation of the Time of Concentration, in minutes
Equation Source Remarks
0.4
tc = 60LA /DS
0.2 Williams (1922) L = basin length (mi), A = basin area (mi2),
D = diameter (mi) of a circular basin of
area A, and S = basin slope (%). The basin
area should be smaller than 50 mi2
tc = KL
0.77
/S n Kirpich (1940) Developed for small drainage basins in
Tennessee and Pennsylvania, with basin
areas from 1 to 112 ac. L = basin length
(ft), S = basin slope (ft/ft), K = 0.0078
and n = 0.385 for Tennessee; K = 0.0013
and n = 0.5 for Pennsylvania. The esti-
mated tc should be multiplied by 0.4 if the
overland flow path is concrete or asphalt,
or by 0.2 if the channel is concrete-lined
tc = (2 LN /3S
0.5 0.47
) Hathaway (1945), Kerby (1959) Drainage basins with areas of less than 10 ac
and slopes of less than 0.01. This is an over-
land flow method. L = overland flow length
from basin divide to a defined channel (ft),
S = overland flow path slope (ft/ft), and
N is a flow retardance factor (N = 0.02
for smooth impervious surfaces; 0.10 for
smooth, bare packed soil; 0.20 for poor
grass, row crops, or moderately rough bare
surfaces; 0.40 for pasture or average grass;
0.60 for deciduous timberland; and 0.80 for
conifer timberland, deciduous timberland
with deep ground litter, or dense grass)
tc = 300( L /S )0.5 Johnstone and Cross (1949) Developed for basins in the Scotie and
Sandusky River watersheds with areas
between 25 and 1,624 mi2. L = basin
length (mi), and S = basin slope (ft/mi)
tc = {(0.007 I + c )/S
0.33
}× Izzard (1946) Hydraulically derived formula. I = effec-
tive rainfall intensity (in/h), S = slope of
( IL /43, 200)−0.67 L /60
overland flow path (ft/ft), L = length of
overland flow path (ft), and c is a roughness
coefficient (c = 0.007 for smooth asphalt,
0.012 for concrete pavement, 0.017 for tar
and gravel pavement, and 0.060 for dense
bluegrass turf)
tc = 0.94 I −0.4 ( Ln /S 0.5 )0.6 Henderson and Wooding (1964) Based on kinematic wave theory for flow
on an overland flow plane. I = rainfall
intensity (in/h), L = length of overland
flow (ft), n = Manning’s roughness coeffi-
cient, S = overland flow plane slope (ft/ft)
AV2 = Q = A 2 g (H − HL ) Federal Aviation Agency (1970) Developed based on airfield drainage data. C =
Rational Method runoff coefficient, L = over-
land flow length (ft), and S = slope (%)
50
ow
n
tio
ow
ad
n
law
iva
me
rfl
t
tte
and
cul
and
d
20
hal way
un
gu
ge
re
ver
gro
low
tilla
r
stu
ate
xo
are
pa
ow sed w
m
nd
10
yb
ds
imu
ss,
rou
arl
) an
gr a
min
s
yg
hee Gra
Ne
Slope, %
ort
eav
or
5
tfl
Sh
low
,h
t
Fal
res
a (s
Fo
are
2
ed
Pav
1
0.5
0.1 0.2 0.5 1 2 5 10 20
Velocity, ft/s
Solution
For the reach from A to B, the average flow velocity is V = 0.7 ft/s. The travel
time for that reach is therefore
50
0
B 900
ft
ft C
2,0
00
ft
D
3,
00
0
ft
E
Similarly, for the reach from B to C, the average flow velocity is V = 2.8 ft/s.
The travel time for that reach is therefore
2/3 2/3
1.49 ⎛⎜ D ⎞⎟ 1.49 ⎛⎜ 3 ⎞⎟
V = ⎜ ⎟ S1/ 2 = ⎜⎜⎝ ⎟⎟⎠ (0.015) = 10 ft/s
1/ 2
n ⎜⎝ 4 ⎟⎠ 0.015 4
Travel time in the open channel from D to E is computed using the bank-full
velocity, again found via Manning’s equation:
1.49 2 / 3 1/ 2 1.49
V = R S = (1.78)2 / 3 (0.005)1/ 2 = 8.2 ft/s
n 0.019
The time of concentration is the sum of the four individual travel times and is
tc = 1,601 s = 0.44 h
Q = CiA
where
Q = runoff rate (ac.in/h, ha mm/h),
C = runoff (abstractions) coefficient,
i = rainfall intensity (in/h, mm/h), and
A = drainage area (ac, ha).
Because 1 ac.in/hr = 1.008 cfs ≈ 1 cfs, engineers performing calculations by hand
typically ignore the conversion factor and simply assume that the discharge Q
is in units of cfs. In SI units, a conversion factor of 0.278 will yield Q in units
of m3/s.
The time of concentration is the smallest time for which the entire basin
is contributing runoff to the basin outlet; therefore the storm duration
must be at least as long as the time of concentration if a steady-state con-
dition is to be achieved. Steady-state conditions also dictate that the storm
intensity be spatially and temporally uniform. It is not reasonable to expect
that rainfall will be spatially uniform over a large drainage basin, or that it
will be temporally uniform over a duration at least as long as the time of
concentration when tc (and hence A) is large. Therefore these conditions
limit the applicability of the Rational Method to small drainage basins. An
upper limit of 200 ac has been suggested by some, but the limit should
really depend on the storm characteristics of the particular locale. These
local characteristics may limit the applicability of the Rational Method to
basins smaller than 10 ac.in some cases.
When several drainage basins (or subbasins) discharge to a common facility
such as a storm sewer or culvert, the time of concentration should be taken
as the longest of all the individual times of concentration and should include
pipe travel times when appropriate. Further, the total drainage area served (the
sum of the individual basin areas) should be no larger than the 200-ac limit (or
smaller where applicable) of the Rational Method.
The basic steps for applying the Rational Method are as follows:
Example 2.10: Computing Flows for Multiple Subbasins with the Rational
Method
Figure 2.13 is a plan view of a storm sewer system draining three subbasins.
Use the Rational Method to determine the peak discharge in each pipe and
size each pipe assuming the pipes flow full. Assume also that the pipes will be
concrete with n = 0.013. Perform the calculations for a storm recurrence inter-
val of 25 years. Subbasin and pipe characteristics and IDF data for the 25-year
event are tabulated as follows:
Subbasin A
Pipe 1
Subbasin B Subbasin C
Pipe 2
Pipe 3
To outfall
Solution
Flow into Pipe 1 occurs from Subbasin A only. Using the time of concentra-
tion as the storm duration, the 25-year rainfall intensity is 5.26 in/h. The peak
discharge used in sizing Pipe 1 is therefore
Q = 0.6(5.26)(6.0) = 19 cfs
Assuming that Pipe 1 is flowing full, its required diameter D may be found
using Manning’s equation as follows:
3/8
⎛ Qn ⎞⎟3 / 8 ⎛⎜ 19(0.013) ⎞⎟
D = ⎜⎜ ⎟ = ⎜ ⎟ = 1.87 ft
⎝ 0.464S ⎟⎠
1 / 2
⎝⎜ 0.464(0.01) ⎠⎟
1 / 2
Pipe 2 is treated the same way as Pipe 1, recognizing that runoff from Subbasin
C only enters Pipe 2. The peak discharge from Subbasin C is Q = 22 cfs, and
the required diameter of Pipe 2 is D = 24 in. The travel time in Pipe 2 is
t = 45 s = 0.75 min.
Pipe 3 must be sized to handle the runoff from all three of the subbasins, which
have a total area of A = 14.5 ac. The runoff coefficient for the combined areas
is computed as a composite value and is
The time of concentration is computed as the longest of the travel times to the
upstream end of Pipe 3. These travel times are (1) the time of concentration
of Subbasin B (10 min), (2) the time of concentration of Subbasin A plus the
travel time in Pipe 1 (20 + 1.0 = 21 min), and (3) the time of concentration
of Subbasin C plus the travel time in Pipe 2 (15 + 0.75 = 15.75 min). Thus
the time of concentration for Pipe 3 is 21 min, and the corresponding rainfall
intensity (by interpolation) is 5.17 in/h.
Q = 0.72(5.17)(14.5) = 54 cfs
The required diameter of Pipe 3 (rounded to the nearest standard size) is 36 in.
Qp = qu APe
where
Qp = the peak flow rate (cfs),
qu = the unit peak flow rate (cfs/mi2/in),
A = the drainage basin area (mi2), and
Pe = the depth of effective precipitation (in).
The unit peak flow rate is estimated as
qu = 10K
where
K = C0 + C1 log10 tc + C2 (log10 tc )2
C0, C1, and C2 are coefficients listed in Table 2.9 as a function of the SCS 24-h
design storm distribution type and the ratio Ia /P, and tc is the time of concentra-
tion (h). P is the 24-h rainfall depth, and Ia is the depth of initial abstractions.
For convenience, the ratio Ia /P is tabulated in Table 2.10 as a function of the
Table 2.9: Coefficients for SCS (NRCS) Peak Discharge Method (SCS, 1996)
Table 2.10: Ia/P for Selected Rainfall Depths and Curve Numbers (Adapted from SCS, 1996)
Curve Number, CN
P (in)
40 45 50 55 60 65 70 75 80 85 90 95
0.50 0.50 0.50 0.50 0.50 0.50 0.50 0.50 0.50 0.50 0.50 0.44 0.21
1.00 0.50 0.50 0.50 0.50 0.50 0.50 0.50 0.50 0.50 0.35 0.22 0.11
1.50 0.50 0.50 0.50 0.50 0.50 0.50 0.50 0.44 0.33 0.24 0.15 0.10
2.00 0.50 0.50 0.50 0.50 0.50 0.50 0.43 0.33 0.25 0.18 0.11 0.10
2.50 0.50 0.50 0.50 0.50 0.50 0.43 0.34 0.27 0.20 0.14 0.10 0.10
3.00 0.50 0.50 0.50 0.50 0.44 0.36 0.29 0.22 0.17 0.12 0.10 0.10
3.50 0.50 0.50 0.50 0.47 0.38 0.31 0.24 0.19 0.14 0.10 0.10 0.10
4.00 0.50 0.50 0.50 0.41 0.33 0.27 0.21 0.17 0.13 0.10 0.10 0.10
4.50 0.50 0.50 0.44 0.36 0.30 0.24 0.19 0.15 0.11 0.10 0.10 0.10
5.00 0.50 0.49 0.40 0.33 0.27 0.22 0.17 0.13 0.10 0.10 0.10 0.10
5.50 0.50 0.44 0.36 0.30 0.24 0.20 0.16 0.12 0.10 0.10 0.10 0.10
6.00 0.50 0.41 0.33 0.27 0.22 0.18 0.14 0.11 0.10 0.10 0.10 0.10
6.50 0.46 0.38 0.31 0.25 0.21 0.17 0.13 0.10 0.10 0.10 0.10 0.10
7.00 0.43 0.35 0.29 0.23 0.19 0.15 0.12 0.10 0.10 0.10 0.10 0.10
7.50 0.40 0.33 0.27 0.22 0.18 0.14 0.11 0.10 0.10 0.10 0.10 0.10
8.00 0.38 0.31 0.25 0.20 0.17 0.13 0.11 0.10 0.10 0.10 0.10 0.10
8.50 0.35 0.29 0.24 0.19 0.16 0.13 0.10 0.10 0.10 0.10 0.10 0.10
9.00 0.33 0.27 0.22 0.18 0.15 0.12 0.10 0.10 0.10 0.10 0.10 0.10
9.50 0.32 0.26 0.21 0.17 0.14 0.11 0.10 0.10 0.10 0.10 0.10 0.10
10.00 0.30 0.24 0.20 0.16 0.13 0.11 0.10 0.10 0.10 0.10 0.10 0.10
10.50 0.29 0.23 0.19 0.16 0.13 0.10 0.10 0.10 0.10 0.10 0.10 0.10
11.00 0.27 0.22 0.18 0.15 0.12 0.10 0.10 0.10 0.10 0.10 0.10 0.10
11.50 0.26 0.21 0.17 0.14 0.12 0.10 0.10 0.10 0.10 0.10 0.10 0.10
12.00 0.25 0.20 0.17 0.14 0.11 0.10 0.10 0.10 0.10 0.10 0.10 0.10
12.50 0.24 0.20 0.16 0.13 0.11 0.10 0.10 0.10 0.10 0.10 0.10 0.10
13.00 0.23 0.19 0.15 0.13 0.10 0.10 0.10 0.10 0.10 0.10 0.10 0.10
13.50 0.22 0.18 0.15 0.12 0.10 0.10 0.10 0.10 0.10 0.10 0.10 0.10
14.00 0.21 0.17 0.14 0.12 0.10 0.10 0.10 0.10 0.10 0.10 0.10 0.10
14.50 0.21 0.17 0.14 0.11 0.10 0.10 0.10 0.10 0.10 0.10 0.10 0.10
15.00 0.20 0.16 0.13 0.11 0.10 0.10 0.10 0.10 0.10 0.10 0.10 0.10
curve number. The ratio Ia/P should not be less than 0.1, nor larger than 0.5.
The time of concentration should be in the range of 0.1 to 10 hours.
When ponding or swampy areas exist in a drainage basin, the peak discharge
should be reduced to account for the temporary storage of runoff. In this case,
an adjusted peak flow Qpa should be computed as
Qpa = FpQp
where Fp is an adjustment factor obtained from Table 2.11. Ponds and swampy
areas lying along the flow path used for computation of the time of concen-
tration should not be included in the area percentage in Table 2.11. Neither
channel nor reservoir routing can be accommodated with this method.
Example 2.11: Computing Peak Discharge Using the NRCS Method
Compute the peak discharge from the 1,000-ac (1.56 mi2) drainage basin
described in Example 2.8, but using AMC-II soil moisture conditions. The
rainfall distribution is an SCS Type II distribution, and the time of concen-
tration of the drainage basin is 0.9 h. There are no ponds or swampy areas in
the basin.
Solution
From Example 2.8, P = 5.0 in and CN = 86 for AMC-II. For this condition,
the maximum possible retention S = 1.63 in, and the effective precipitation is
Pe = 3.46 in.
From Table 2.10, the ratio Ia/P is 0.10; and from Table 2.9, the coefficients C0,
C1, and C2 are C0 = 2.55323, C1 = −0.61512, and C2 = −0.16403, respec-
tively. Therefore
K = 2.55323−0.61512 log10(0.9)−0.16403[log10(0.9)]2 = 2.58104
De
Drainage area Dr
Ds
Dp
Dli
Volume runoff, Vr Di
Peak
flow
Flow rate
Runoff volume
(integrated area)
Time
t (min) P (in)
0–10 0.24
10–20 0.46
20–30 1.17
30–40 0.58
40–50 0.35
50–60 0.17
Solution
The interception capacity of 0.3 in is subtracted first. Because 0.24 in of
rainfall occurs during the first 10 min of the storm, all of that rainfall plus
an additional 0.06 in of rainfall occurring in the second 10 min of the storm
t (min) P (in)
0–10 0
10–20 0.4
20–30 1.17
30–40 0.58
40–50 0.35
50–60 0.17
10 0 1.5
20 10 1.1 0.22
30 20 0.84 0.16
40 30 0.66 0.13
50 40 0.54 0.1
60 50 0.46 0.08
Subtraction of the infiltration depth in each time interval from the correspond-
ing rainfall depth remaining after interception leads to the following hyetograph
(any negative values produced should be set equal to zero):
t (min) P (in)
0–10 0
10–20 0.18
20–30 1.01
30–40 0.45
40–50 0.25
50–60 0.09
t (min) Pe (in)
10–20 0
10–20 0
20–30 0.99
30–40 0.45
40–50 0.25
50–60 0.09
The rainfall hyetograph and the effective rainfall hyetograph are illustrated in
Figures 2.15 and 2.16. Note that the effective rainfall intensities are less than the
actual rainfall intensities. Also note that effective rainfall, and hence surface run-
off, does not begin (in this example) until 20 min after the beginning of the storm.
8
7
6
Intensity, in/h
4
3
2
1
0
0 to 10 10 to 20 20 to 30 30 to 40 40 to 50 50 to 60
Time, min
8
7
Unit Hydrographs
A unit hydrograph is an important intermediate step in computing a full
runoff hydrograph for design or analysis purposes. After a unit hydrograph
is established for a particular basin, the data can be used to compute the
complete runoff hydrographs resulting from various rainfall events. The basic
theory rests on the assumption that the runoff response of a drainage basin
to an effective rainfall input is linear. Practically speaking, this assumption
means that the concepts of proportionality and superposition can be applied.
For example, it means that the direct runoff caused by 2 in. of effective rain-
fall in a given time interval is four times as great as that caused by 0.5 in. of
effective rainfall in the same amount of time. It also means that the total
amount of time for the basin to respond to each of these rainfall depths is the
same. Consequently, the base length, in hours, of each of the direct runoff
hydrographs is the same.
tL
1 in of effective
rainfall in Δt time units
0 Δt tp tB
Time
1in
runoff Contributing area Flo
Ru w
no
Total ff Flow
volu
water me
applied
Abstractions
Runoff volume
the sprinklers are turned off. Δt represents the duration of time that the sprin-
klers are on. A stream gauge is set on the downstream end of the contributing
drainage area to record the runoff hydrograph resulting from the 1 in of runoff.
The resulting hydrograph represents the unit hydrograph generated for the
drainage area.
The time tL in Figure 2.17, which is the length of time from the center of Δt to the
peak discharge on the unit hydrograph, is the basin lag time. This response time
can be thought of as an approximate average of all possible travel times for runoff
in a drainage basin. In practice (and as shown in the figure), it is usually assumed
to be the amount of time between the center of mass of a pulse of effective rainfall
and the peak of the resultant hydrograph. The basin lag time is often used instead
of time of concentration when estimating the complete runoff hydrograph.
Two basic categories of unit hydrographs exist: (1) unit hydrographs for gauged
watersheds and (2) synthetic unit hydrographs. Most commonly, rainfall and
runoff records for a drainage basin do not exist, and one must resort to syn-
thesis of a unit hydrograph based on information that can be gathered about
the basin. Extensive literature exists on various ways to calculate synthetic unit
hydrographs, including procedures proposed by Clark, Snyder, and Singh. This
text introduces the development of synthetic unit hydrographs by applying the
NRCS (SCS) method. Further details on these and other synthetic unit hydro-
graphs methods can be found in Stormwater Conveyance Modeling and Design
(Haestad Methods and R. Durrans, 2003)
1.2
1.0
0.8
Q / Qp
0.6
0.4
0.2
0.0
0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0
t / tp
1.6 0.560
Δt
tp = + tL
2
The peak discharge Qp (in cfs) of the synthetic unit hydrograph is calculated as
645.33 KAQ
Qp =
tp
where
645.33 = conversion factor = (1 h/3,600 s) (1 ft/12 in) [(5,280 ft)2/1 mi2];
K = 0.75 is constant based on geometric shape of dimensionless unit hydro-
graph (NEH-4);
Q = 1 in. is the runoff depth for unit hydrograph calculation;
A = the drainage basin area (mi2); and
tp = the time to peak (h).
Simplifying yields
484 A
Qp =
tp
Solution
Employing the guidelines for estimation of Δt, its value should be no larger
than
0.133tc = 0.133(90) = 12 min
tL = 0.6tc = 54 min
Table 2.13: SCS (NRCS) Synthetic Unit Hydrograph Computations for Example 2.12
Ordinates of the synthetic unit hydrograph are given in Table 2.13. The first
column of the table is the time t, in minutes, and is tabulated in Δt = 10 min
intervals. The second column is the dimensionless time ratio t/tp, where tp =
59 min. The third column is the dimensionless discharge ratio and is deter-
mined using the dimensionless time ratio and interpolation from Table 2.12.
The fourth and last column contains the ordinates of the 10-min unit hydro-
graph, which are computed as the products of the dimensionless discharge
ratios and Qp = 740 cfs (Figure 2.20).
800
700
Discharge, cfs
600
500
400
300
200
100
0
0 50 100 150 200 250 300
Time, min
Discrete Convolution
A unit hydrograph represents the runoff hydrograph from a drainage basin sub-
jected to 1 in. of runoff (effective precipitation) applied over a duration of Δt.
However, a runoff (effective rainfall) hyetograph typically contains many
Δt time periods, each of which has its own associated runoff depth. In discrete
convolution, the direct runoff hydrograph resulting from a complete rainfall
hyetograph is computed by applying a unit hydrograph to each discrete time
step within the hyetograph.
To conceptualize discrete convolution, refer to Figures 2.21 and 2.22. The
first figure illustrates a runoff (effective rainfall) hyetograph consisting of
np = 3 rainfall pulses, each of duration Δt = 10 min, with the depth of each
pulse denoted by Pi (i = 1, 2, …, np). The second figure shows a unit hydro-
graph with nu = 11 nonzero ordinates shown at Δt time intervals, with the
nonzero ordinates denoted by Uj (j = 1, 2, …, nu). (Note that the ordinates
of the unit hydrograph are in units of cfs per inch of runoff, and thus the
notation of this subsection is different from the previous subsection in which
unit hydrograph ordinates were denoted by discharge Q.) The time incre-
ment Δt used for development of the runoff (effective rainfall) hyetograph
must be the same as the Δt duration of excess rainfall used to create the unit
hydrograph.
Because unit hydrograph theory assumes that drainage basins behave linearly,
the ordinates Qj of the runoff hydrograph produced by P1 inches of effective
rainfall must be equal to Qj = P1Uj. These ordinates are shown in the third
column of Table 2.14 (the first and second columns of the table reproduce the
information in Figure 2.22). These ordinates were calculated by multiplying
each unit hydrograph ordinate by the runoff depth for pulse P1 = 0.2 in. The
values in column three are shifted down such that the time at which the runoff
hydrograph in column three begins is the same as the time at which the first
pulse of effective rainfall begins.
The fourth and fifth columns in Table 2.14 are computed by multiplying the
unit hydrograph ordinates in the second column by the effective rainfall depths
Runoff hyetograph
1.6
1.4
1.2
Runoff depth, in
1.0
0.8 P2
0.6
0.4
P3
0.2
P1
0.0
0 10 20 30 40 50 60
Time, min
Unit hydrograph
300
250 U4
Discharge, cfs/in
200 U5
U3
150
U6
100 U2 U7
U8
50 U1 U9
U10 U11
0
0 20 40 60 80 100 120
Time, min
600
500
Total runoff
400
Flow, cfs
Pulse P2
300
Pulse P3
200
Pulse P1
100
0
0 20 40 60 80 100 120 140 160
Time, min
P2 = 1.5 in and P3 = 0.7 in, respectively. Again, the time at which each runoff
pulse’s hydrograph begins corresponds to the time at which each runoff (effec-
tive rainfall) pulse begins. The sum of the runoff hydrograph ordinates in each
row of the table (shown in the sixth column) is an ordinate of the direct runoff
hydrograph caused by the complete rainfall event.
Figure 2.23 shows the runoff hydrograph resulting from each rainfall pulse and
the total runoff hydrograph for the entire rainfall event.
Example 2.14: Computing a Runoff Hydrograph
Compute the direct runoff hydrograph for the Memphis, Tennessee drain-
age basin described in Example 2.13. Use the unit hydrograph developed for
the basin in that example, and use the runoff (effective rainfall) hyetograph
developed in Example 2.12.
Solution
There are a total of np = 4 runoff (effective rainfall) pulses, and a total of
nu = 26 nonzero unit hydrograph ordinates. Thus there will be a total of nq = 29
nonzero direct runoff hydrograph ordinates.
Table 2.15 illustrates the tabular calculation (discrete convolution). The first
column is the time since the beginning of rainfall in Δt = 10-min increments
(the duration of the unit hydrograph and the duration of the effective rain-
fall pulses). The second column contains the ordinates of the unit hydrograph
from Example 2.13. From the third to the sixth columns are the unit hydro-
graph ordinates multiplied by the effective rainfall depths. Note that the first
entry (the first zero) in each column corresponds to the time at which the
corresponding effective rainfall pulse begins (t = 20 min for the first pulse,
t = 30 min for the second pulse, and t = 40 min for the third pulse; see Figure
2.16). The seventh column is the sum of the previous four and is the direct
runoff hydrograph.
2.5 Problems
1. The following data is from Huff and Angle (1992). The numbers indicate
the total rainfall (in) expected in Chicago, Illinois for storm recurrence
intervals of 2, 5, 10, 25, 50, and 100 years and durations from 5 min to
24 h. Plot a series of IDF curves showing the rainfall intensities for each
recurrence interval for storm durations between 5 min and 2 h.
2. Use the data in Problem 1 and the SCS Type II distribution to develop a
25-year, 24-h storm for Chicago.
3. During a 40-mm rainfall event, the interception capacity of the water-
shed is estimated to be 11 mm, the depression storage is 9 mm, and the
infiltration is 7 mm. What is the total volume of runoff from the 2.56 ha
watershed?
4. A proposed development consists of the following land use areas:
Table 3.1: Suggested Minimum Design Frequencies and Spreads for Gutter Sections on Grade
Road Classification Design Return Design Spreads
Period (Year)
High volume or <70 km/h (45 mph) 10 Shoulder + 1 m (3 ft)
bidirectional >70 km/h (45 mph) 10 Shoulder
Sag point 50 Shoulder + 1 m (3 ft)
Collector Low volume 10 1/2 Driving lane
High volume 10 Shoulder
Sag point 10 1/2 Driving lane
Local streets Low volume 5 1/2 Driving lane
High volume 10 1/2 Driving lane
Sag point 10 1/2 Driving lane
T
T
Qs Qw
Sx
Sx
a Sw
where
Kc = 0.376 (0.56 in U.S. units),
n = Manning’s coefficient, and
SL = road longitudinal slope (m/m, ft/ft).
Using this equation, the spread, T, in a uniform gutter section can be explicitly
calculated for a given flow rate, Q. However, in the case of composite gutters,
T can no longer be expressed as an explicit function of Q. Rather, an iterative
process is required to calculate the spread.
On a road with a grade, the spread will be at a maximum just upstream of the
inlet. Note that the spread at this location is independent of the inlet’s ability to
capture flow and its efficiency. The spread for a specific roadway is a function of
the discharge in the gutter. To decrease the flow to an inlet, reduce inlet spacing
so that they serve as collection points for smaller watersheds.
For composite sections, the variable E0 is introduced to account for the added
conveyance of a depressed gutter. E0 is the ratio of flow in the depressed sec-
tion, Qw, to the total gutter flow, and is expressed as:
−1
⎧⎪ ⎡⎛ ⎞⎟2.67 ⎤⎫⎪⎪
⎪ S ⎢ ⎜ S S ⎥ ⎪
E0 = ⎪⎨1 + W ⎢⎜⎜1 +
W X ⎟
⎟ − 1⎥⎬
⎪⎪ SX ⎢⎜⎝ (T W ) − 1⎟⎟⎠ ⎥ ⎪⎪
⎪⎩ ⎢⎣ ⎥⎦ ⎪⎭
In the case of a grate inlet in a uniform section, the variable E0 is also useful,
with the frontal flow, Qw, now defined as the gutter flow contained in the
width of the grate. In this case, E0 becomes:
2.67
⎛ Wg ⎞⎟
E0 = 1 − ⎜⎜1 − ⎟⎟
⎜⎝ T ⎠
where
E = inlet efficiency,
Qi = intercepted flow (m3/s, cfs), and
Q = total gutter flow (m3/s, cfs).
Flow that is not intercepted by a drainage inlet on grade is bypassed and carried
over to another inlet downstream or is “lost” to a stream or pond, for example.
W
L
Reticuline
Figure 3.2: Grate inlet in a gutter and some typical grate types.
grate inlets have a higher propensity to clog than other types of inlets. If debris
is prevalent in runoff at the point of design, adequate provisions must be made
to account for inlet clogging, such as utilizing a combination (grate and curb)
inlet at these points.
As shown in Figure 3.1, the total gutter flow is composed of frontal flow, Qw
(the flow in the depressed gutter or over the grate width), and side flow, Qs. The
total efficiency of the grate inlet is determined by calculating the grate’s ability
to capture frontal flow and side flow.
Rf is the ratio of intercepted frontal flow to the total frontal flow and is
expressed as:
R f = 1 − K cf (V − V0 )
where
Kcf = 0.295 (0.090 in U.S. units),
V = average velocity in the gutter at the location of the inlet
(m/s, ft/s), and
V0 = splash-over velocity of the inlet (m/s, ft/s).
The splash-over velocity is the minimum velocity of the gutter flow capable of
inducing enough momentum for some of the flow to skip over the grate open-
ing and be carried over downstream. The splash-over velocity is a function of the
grate type and the grate length and can be found in Appendix A of the HEC-22
manual (Brown et al., 2009) or obtained from the grate manufacturer. If the
gutter velocity is less than the splash-over velocity, all frontal flow is intercepted
and Rf equals 1.0.
The ratio, Rs, of side flow intercepted to total side flow is expressed as:
−1
⎛ K V 1.8 ⎞
Rs = ⎜⎜⎜1 + cs 2.3 ⎟⎟⎟
⎜⎝ Sx L ⎟⎠
where
Kcs = 0.0828 (0.15 in U.S. units) and
L = grate length (m, ft).
The total intercepted flow is expressed as:
Qi = Qw R f + Qs Rs
Sx
a Sw
a’ S w’
Sw’ = Slope of local depression (m/m, ft/ft)
a = Gutter depression (mm, in)
a’ = Total depression (mm, in)
The efficiency of a curb inlet is based on the ratio of the actual inlet length
to the inlet length necessary to capture 100% of the total runoff. Curb inlets
are often inset into the pavement to create a local depression. A local depres-
sion, as shown in Figure 3.4, is a depression of the gutter at the location of
the inlet only, as opposed to a gutter depression, which is continuous along
the curb.
The curb opening length, LT, that would be required to intercept 100% of a
flow, Q, on a roadway section with a uniform cross-slope is defined as:
⎛ 1 ⎞⎟0.6
LT = K C Q 0.42SL0.3 ⎜⎜ ⎟
⎜⎝ nSx ⎟⎟⎠
Se = Sx + Sw ′ E0
where
Sw′ = gutter cross-slope at the inlet location measured from the pavement
cross-slope, Sx (m/m, ft/ft) and
E0 = ratio of flow in the depressed section to the total gutter flow upstream of
the inlet (does not account for local depression).
To calculate LT with a composite gutter, replace the road cross-slope, Sx, with
the equivalent slope, Se, in the equation solving for LT.
The efficiency of a curb opening on grade shorter than the required length for
total interception is:
1.8
⎛ L⎞
E = 1 − ⎜⎜1 − ⎟⎟⎟
⎜⎝ LT ⎟⎠
(A) Lc
W
e Lg
charg g
Dis
Lc
(B)
W Lg
e g
arg
D isch
The curb inlet is often extended upstream of the grate for more efficient
removal of debris than when the curb opening length equals the grate length.
The total flow intercepted by this configuration is calculated as the sum of
• the flow intercepted by the portion of the curb opening located upstream
of the grate and
• the flow that bypassed the upstream curb opening and is intercepted by
the grate alone. (the flow intercepted by the portion of the curb opening
adjacent to the grate is neglected.)
Qiw = Cw Pd 1.5
where
Qiw = flow intercepted by the inlet operating as a weir (m3/s, ft3/s);
Cw = weir coefficient, which varies depending on the flow condition and inlet
structure;
P = perimeter of the inlet (m, ft); and
d = flow depth at the curb (m, ft).
Note that for a grate inlet, the perimeter does not include the length along the
curb. Also, if the gutter is depressed (locally or continuously), the perimeter, P,
of the grate is calculated as:
P = L + 1.8W
where
L = grate length (m, ft) and
W = grate width (m, ft).
The depth, d, for both types of inlets is measured from the projected pavement
cross-slope. For a curb inlet, the perimeter is equivalent to the length of the
inlet.
Qio = Co A (2 gd0 )
0.5
where
Qio = flow intercepted by the inlet operating as an orifice (m3/s, ft3/s),
C0 = orifice coefficient (varies based on the class of inlet and its configuration),
A = area of the opening (m2, ft2),
g = 9.81 m/s2 (32.16 ft/s2 in U.S. customary units), and
d0 = effective head at the orifice (m, ft).
Note that for a grate inlet, the effective head, d0, is simply the water depth
along the curb. For a curb inlet, the effective head is expressed as:
h
d0 = di − (sinθ)
2
where
di = depth at lip of curb opening (m, ft),
h = curb throat opening height (m, ft) (see Figure 3.6), and
θ = inclination of the curb throat measured from the vertical direction
as shown in Figure 3.6.
do
di
h
Horizontal throat
do = di – (h/2)
θ = p/2
do
θ
h
Inclined throat
do = di – h/2(sinθ)
do
Vertical throat
do = di
θ=0
h
Grates alone are not typically recommended for installation in sags because of
their propensity to clog and exacerbate ponding during severe weather. A com-
bination inlet may be a better alternative. At low flow depths, the capacity of a
combination inlet for which the grate inlet length equals the curb inlet length
is equivalent to the capacity of the grate inlet alone. At higher flow depths for
the same type of inlet, both the curb inlet and grate inlet act as orifices. The
total intercepted flow is then calculated as the sum of the flows intercepted by
the grate and curb openings.
Set the variable you want to solve for by selecting it from the pull-down menu
in the Solve For: field. Then, click Solve. The calculated values will appear in
the fields highlighted in yellow.
Normal depth
Critical depth
depth throughout the length. For example, a high tailwater elevation may force
the depth to be above normal depth at the downstream end of a pipe, as shown
in Figure 3.7.
When the difference in water depth from one end of a conduit to the other is sig-
nificant, a gradually varied flow profile analysis is required. Gradually varied flow
profile analysis is the process of splitting a channel length into segments and ana-
lyzing each segment separately. It is based on the following several assumptions:
• The head loss within any given calculation segment is the same as for uni-
form flow conditions.
• The velocity is the same across the entire cross section.
• The slope of the conduit is less than 10%.
• The roughness coefficient is constant throughout the reach under consider-
ation and is independent of the depth of flow.
• The depth of flow changes gradually along the length of the conduit, start-
ing from some controlling boundary condition (usually the tailwater eleva-
tion); there are no sudden increases or decreases in depth.
• If the pipe is sufficiently long, the depth of flow will approach normal depth.
Flow Classification
The first step in performing a gradually varied flow analysis is to identify the
flow classification that is expected to occur in the conduit based on the slope of
the channel, the normal depth, the critical depth, and the controlling bound-
ary condition.
Slope Classification
Once the normal depth and critical depth have been computed for the section,
the conduit slope can be determined. If the normal depth is above critical
depth, the slope is said to be mild. If the normal depth is equal to the critical
depth, the slope is said to be critical. If normal depth is below critical depth,
the slope is said to be steep.
There are two other slope types to consider: a horizontal channel and a chan-
nel on an adverse slope (an “uphill”-sloping channel). For these slope types,
normal depth is undefined.
For a gradually varied flow profile, the first letter of the slope type is used in
the identification of the profile. For example, a channel with a hydraulically
steep slope is a type S, a channel with a hydraulically mild slope is a type M,
and so forth.
Profile Classifications
Once the slope classification and flow zone have been determined, the profile
type can be defined, and the engineer can determine how to proceed with the
hydraulic grade computations. Figure 3.8 shows the basic profile types.
To perform the computations, the engineer must determine from the profile
type whether the flow is subcritical or supercritical (based on the location of
actual depths compared to critical depth). In order to prevent excessive veloc-
ities that could cause pipe scour or channel erosion, most storm sewers are
designed with mild slopes to carry subcritical flows. For this type of flow, the
hydraulic control is at the downstream end of the section, and the profile cal-
culation proceeds toward the upstream end. When the flow depth is above
normal depth (as in an M1 profile), this type of analysis is called a backwater
analysis. When the flow depth is between critical depth and normal depth, it is
called a drawdown analysis.
When supercritical flows are encountered, the controlling section is at the
upstream end of the conduit, and the computations proceed from upstream to
downstream. This calculation type is a frontwater analysis.
M1 M2 M3
yn yn yn
Mild slope
yn > yc
yc yc yc
C1 C2 unstable C3
Critical slope
S1 S2 S3
yc yc yc
Steep slope
y n < yc
yn yn yn
H1 none H2 H3
Horizontal slope
yc yc yc
A1 none A2 A3
Adverse slope
yc yc yc
Energy Balance
Even for gradually varied flow, the solution is still a matter of balancing the
total energy between the two ends of each segment. The energy equation as it
relates to each end of a segment can be written as (note that the pressures for
both ends are zero, because it is free surface flow):
V12 V2
Z1 + = Z 2 + 2 + HL
2g 2g
where
Z1 = hydraulic grade at the upstream end of the segment (m, ft),
V1 = velocity at the upstream end (m/s, ft/s),
Z2 = hydraulic grade at the downstream end of the segment (m, ft),
V2 = velocity at the downstream end (m/s, ft/s),
HL = head loss due to friction (other losses are assumed to be zero) (m, ft), and
g = gravitational acceleration constant (m/s2, ft/s2).
The friction loss is computed based on the average rate of friction loss along the
segment and the segment length. This relationship is:
S1 + S2
HL = SAvg 3Δx = Δx
2
where
HL = loss across the segment (m, ft),
SAvg = average friction loss (m/m, ft/ft),
S1 = friction slope at the upstream end of the segment (m/m, ft/ft),
S2 = friction slope at the downstream end of the segment (m/m, ft/ft), and
Δx = length of the segment being analyzed (m, ft).
The conditions at one end of the segment are known (through assumption or
from a previous calculation). Because the friction slope is a function of velocity,
which is a function of depth, the depth at the other end of the segment can be
found through iteration. There are two primary methods to obtain the solu-
tion: the standard-step method and the direct-step method.
Standard-Step Method
This method involves dividing the channel into segments of equal known
lengths and solving for the unknown depth at one end of the segment (starting
with a known or assumed depth at the other end). The standard-step method
is the most popular method of determining the flow profile because it can be
applied to any channel, not just prismatic channels.
Direct-Step Method
The direct-step method is based on the same basic energy principles as the
standard-step method, but it takes a slightly different approach to the solution.
Instead of assuming a segment length and solving for the depth at the end of
the segment, the direct-step method assumes a depth and then solves for the
segment length.
Hydraulic jump S1
S2
Critical depth
Normal depth
discussed. This section explores those cases in which profile types are mixed within
the same section and the steps that can be taken to analyze these occurrences.
Sealing Conditions
There may be conditions whereby part of a pipe or culvert section is flowing
full, while part of the flow remains open. These conditions are called sealing
conditions and are analyzed in separate parts. The portion of the section flow-
ing full is analyzed as pressure flow, and the remaining portion is analyzed as
gradually varied flow.
Hydraulic Jumps
When flow passes rapidly from supercritical to subcritical flow, a phenomenon
called a hydraulic jump occurs. In storm sewer networks, this often occurs when
a steep pipe discharges into a particularly high tailwater, as shown in Figure 3.9.
There are significant losses associated with hydraulic jumps due to the amount
of mixing and turbulence that occur. These forces are also highly erosive, so
engineers typically try to prevent jumps from occurring in storm sewer sys-
tems, or at least predict the location of jumps in order to provide adequate
channel, pipe, or structure protection.
Hydrology Model
As the runoff from a storm event travels through a storm sewer, it combines
with other flows, and the resulting flows are based on the overall watershed
characteristics. As with a single watershed, the peak flow at any location within
the storm sewer is assumed to occur when all parts of the watershed are con-
tributing to the flow. Therefore the rainfall intensity that produces the largest
peak flow at a given location is based on the controlling system time at that
same location.
The controlling system time is the larger of the local time of concentration
(to a single inlet) and the total upstream system time (including pipe travel
times). The controlling time is used for computing the intensity (and therefore
the flow) in the combined system.
Solution
The total CA can be found by simply summing the CA values from the two
inlets. The storm duration, however, must be found by comparing the local
time of concentration at the second inlet to the total time for flow from the
upstream inlet to reach the downstream inlet.
Total CA = 1.23 ac + 0.84 ac = 2.07 ac
Upstream time = 8 min + 2 min = 10 min
The total upstream flow time of 10 min is greater than the local time of con-
centration at the downstream inlet (9 min). The 10-min value is therefore the
controlling time and should be used as the duration of the storm event.
that include circular pipes, pipe arches, boxes, and more. Flow calculations
handle pressure and varied flow situations, including hydraulic jumps, back-
water, and drawdown curves. StormCAD’s flexible reporting features allow you
to customize and print the design and analysis results in report format or as a
graphical plot.
Profiles
StormCAD also includes an option to automatically generate storm sewer
profiles—longitudinal plots of the storm sewer. Profiles allow the design engineer,
the reviewing agency, the contractor, and others to visualize the storm system. They
are useful for viewing the hydraulic grade line and determining if the proposed
storm sewer is in conflict with other existing or proposed underground utilities.
3.10 Problems
Solve Problems 1–3 using FlowMaster and Problems 4–11 using StormCAD.
1. The spread on a proposed road paved with smooth asphalt and having a
1.0% longitudinal slope is limited by local regulations to 3.0 m for the
sake of automobile safety. The road cross-slope is 2.5%.
(a) What is the maximum allowable discharge on the road if the gutter
section is uniform? What is the depth of the flow measured at the
curb for this rate of discharge?
(b) An alternative design proposes the addition of a gutter depression
with a width of 1.5 m and a gutter cross-slope of 4.0%. What is the
maximum allowable discharge in the gutter under these conditions?
What is the depth of the flow at this rate of discharge?
(c) Does the addition of the gutter depression increase or decrease the
capacity of the gutter? Why?
2. A proposed gutter section has a road cross-slope of 0.015 m/m and a
gutter cross-slope of 0.045 m/m. The gutter is 1.20-m wide and is on a
longitudinal slope of 0.02 m/m. The Manning’s coefficient for the road is
0.013. The predetermined runoff in the gutter is 0.05 m3/s.
(a) What is the efficiency of a P-50 mm×100 mm grate inlet that is
0.70-m long and 0.5-m wide? Assume no clogging. What is the
intercepted flow? What is the bypassed flow?
(b) What is the efficiency of a 0.70-m-long curb inlet on the same gutter
section with no local depression? How does this efficiency compare
to the efficiency of the grate inlet in part (a)?
(c) For the curb inlet in part (b), add a local depression that is 20-mm
deep and 0.7-m wide. What is the efficiency of the inlet? How does
this efficiency compare to the grate inlet in part (a) and the curb
inlet in part (b)?
3. An inlet in sag must collect the 50-year peak runoff of 8.0 cfs. The gut-
ter section has a road cross-slope of 2.0% and a gutter cross-slope of
4.0%. The gutter width is 3.0 ft. The local depression is 2.0-in. deep and
3.0-ft wide.
(a) What length combination inlet is necessary to maintain a spread of
8.0 ft at the section? The grate is a 3.0-ft wide, P-50-mm type, and
the curb inlet has an opening height of 0.7 ft and a vertical throat
configuration. Assume the length of the grate will equal the length
of the curb opening.
(b) How big would the curb inlet from part (a), acting alone, have to be
to maintain the spread width of 8.0 ft under the same conditions?
How long would the grate inlet in part (a) have to be to maintain the
same spread? Which inlet type is the best solution?
Solve the following problems using StormCAD.
4. Lay out the storm sewer system shown below in StormCAD and enter
the data for the network from the tables below. Calculate the results using
Manning’s equation and the 10-year storm event data tables that follow.
Inlet I-1 is on grade with a longitudinal slope of 3%, whereas inlet I-2 is
in sag. The tailwater condition at the outlet is free outfall. Assume that
there is no clogging of the inlets.
Fill out an answer table like the one below for each of the following
situations:
(a) Assume uniform gutters with a slope of 0.02 m/m.
(b) Assume continuously depressed gutters with a road cross-slope of
0.02 m/m, a gutter cross-slope of 0.04 m/m, and a gutter width
of 0.8 m.
(c) Assume continuously depressed gutters as described in (b), as well as
inlet lengths increased from 1.1 to 1.6 m.
(d) Explain the reasons for the differences between the three resulting
tables.
5. Enter and calculate the storm sewer network shown below using
Manning’s equation and the rainfall data used in Problem 4. The tailwa-
ter condition is free outfall. Assume no clogging of the inlets. Answer the
questions that follow.
I-3
P-2
I-2
P-3
P-1
O-1
I-1
(d) The above analyses are run using a default Manning’s n of 0.013.
Many drainage design manuals propose a less conservative design
roughness of 0.012. Reanalyze the improved system under 25-year
flows using n = 0.012. How does this change influence the pre-
dicted performance of the system?
I2
P2
P1 P3 Outlet
I1
I3
7. You have been asked by the lead project engineer for a water supply util-
ity to design the stormwater collection system for the proposed ground
storage tank and pump station facility shown in the layout. Pipe lengths
for P-1, P-2, P-3, and P-4 are 88, 92, 185, and 46 ft, respectively. See the
CAD drawing and data for the system layout. Assume C = 0.3 for pervi-
ous areas and C = 0.9 for impervious areas.
(a) Using the StormCAD program’s Automatic Design feature, size
the system using the following design data. Use concrete pipe
(n = 0.013) and the 25-year intensity–duration–frequency data
provided in Problem 6 (Hint: StormCAD can mix SI and U.S.
customary units). The top of bank elevation at the outfall ditch is
846.1 ft. The outfall pipe invert must be located at or above eleva-
tion 838.0 ft. Assume that the water surface elevation at the outfall
is 842.0 ft, and that the pipes should have matching soffit (crown)
elevations at every structure. Present your design in tabular form and
provide a profile plot of your design.
(b) During agency review, the county engineer requests that the water
utility and the county work cooperatively to accommodate the
planned construction of an elementary school nearby by increas-
ing the size of the proposed storm system so that it can handle the
design runoff from the school. The county engineer performs his
own calculations and asks that you increase the size of pipes P-3
and P-4 to handle an external CA of 9.5 acres with a time of con-
centration of 12 min. Using StormCAD, introduce the additional
flow at inlet I-3 and revise the facility design using the Automatic
Design functionality of the program. Are all the design constraints
met? What can you say about the flow conditions in pipe P-3 and
pipe P-4? In the rest of the pipes?
(c) If necessary, manually fine-tune and revise the design to meet all
design criteria. Document your design as in part (a).
Water
storage
tank
I-1 P-1
I-2
P-2
I-4
I-3 P-3
P-
4
Outlet
9. The system described below drains into a river. Because the water surface
elevation of the river varies, analyze the following system with different
tailwater elevations. Use the Darcy–Weisbach equation because it is more
appropriate for pressure flow than Manning’s equation. All pipes are con-
crete (roughness = 0.122 mm), circular pipes. The pipe and inlet data
follow. The ground elevation at the outlet is 312.0 m.
What is the hydraulic grade at the entrance of inlet I-2 under the follow-
ing tailwater conditions for the 10-year storm event?
(a) Tailwater elevation = 310.0 m
(b) Tailwater elevation = 311.0 m
(c) Tailwater elevation = 313.0 m
I-3
I-1
P-2
P-1
I-2
P-3
I-4
4
P-
Outlet
I-1
P-1
P-2 I-2
I-3
P-3
Outlet
(b) If the existing system were constructed using CMP (n = 0.024) instead
of PVC, how would the system perform in a 25-year storm event?
(c) If the outfall were discharging to a pond that had a water surface elevation
equal to that of the ground (using the original PVC pipes), how would
the system perform in a 25-year storm? A 50-year storm?
(d) Using the tailwater conditions specified in part (a), if the park were paved
(with a C = 0.9 and time of concentration = 8 min), how would the
system perform in a 25-year storm? In a 50-year storm?
11. Use the same network and data provided in part (a) of Problem 10.
During a 25-year storm, it was determined that inlet 1 is not capable of
capturing all of the flow produced by the runoff for its catch basin. A
flow of 5 cfs is bypassed to inlet 2, which is in sag (all gutter flow is cap-
tured). Hint: You can apply negative and positive additional carryovers
to inlets I-1 and I-2 under the Catchment tab.
(a) What is the gutter flow captured by inlet 1?
(b) What is the gutter flow captured by inlet 2?
(c) How is the system performing?
(d) A flow monitoring study was performed and the following data was
collected:
Analyze the system using these known flow values for both the 25- and 50-year
storms. How is the system performing? Report the flows, upstream velocities,
and hydraulic grade lines for each pipe.
Culvert Hydraulics
Obtaining accurate solutions in culvert hydraulics can be a formidable compu-
tational task. Culverts act as a significant constriction to flow and are subject to
a range of flow types, including both gradually varied and rapidly varied flows.
It is this mix of flow conditions and the highly transitional nature of culvert
hydraulics that makes the hydraulic solutions so difficult. For this reason, the
documented approach is to simplify the hydraulics problem and analyze the
culvert system using two different assumptions of flow control:
• Outlet control assumption—Computes the upstream HW depth using
conventional hydraulic methodologies that consider the predominant
losses due to the culvert barrel friction, as well as the minor entrance and
exit losses. The tailwater condition during the design storm has an impor-
tant effect on the culvert system.
• Inlet control assumption—Computes the upstream HW depth resulting
from the constriction at the culvert entrance, while neglecting the culvert
barrel friction, tailwater elevation, and other minor losses.
The controlling HW depth is the larger value of the computed inlet control and
outlet control HW depths. Because the culvert system may operate under inlet
control conditions for a range of flow rates and under outlet control conditions
for another range, calculations must be performed for both control conditions.
V2 V2
HWO + 2ug = TW + 2dg + HL
where
HWO = headwater depth above the outlet invert (m, ft);
Vu = approach velocity (m/s, ft/s);
TW = tailwater depth above the outlet invert (m, ft);
Water (A)
surface
H
W.S.
HW
(B)
H
W.S. W.S.
HW
(C)
W.S.
H
HW
W.S.
(D)
W.S. E.G.
H
HW
W.S.
(E)
E.G. H
W.S.
HW W.S.
HWO = TW + HL
where
HWO = headwater depth above the outlet invert (m, ft),
TW = tailwater depth above the outlet invert (m, ft), and
HL = sum of all losses as listed above (m, ft).
Friction Losses
Culverts are frequently hydraulically short (entrance and exit losses exceed pipe
friction losses), and uniform flow depths are not always achieved. For this rea-
son, gradually varied flow methods are well suited to the analysis. For a more
detailed description of gradually varied flow, see Chapters 3.6 and 9.4
⎛V 2 ⎞
He = ke ⎜⎜⎜ ⎟⎟⎟
⎜⎝ 2 g ⎟⎠
where
He = entrance loss (m, ft),
ke = entrance loss coefficient,
V = velocity just inside the barrel entrance (m/s, ft/s), and
g = gravitational acceleration (m/s2, ft/s2).
The entrance loss coefficient varies depending on the type of inlet that is pres-
ent. The smoother the transition from the channel or pond into the culvert,
the lower the loss of coefficient is. Values for the coefficients are presented in
Table 4.1.
*
Note: “End-section conforming to fill slope,” made of either metal or concrete, are the sections commonly available from
manufacturers. From limited hydraulic tests, they are equivalent in operation to a headwall in both inlet and outlet
control. Some end sections with a closed taper in their design have superior hydraulic performance.
⎡ 2 Vd 2 ⎤
HO = 1.0 ⎢V − 2g ⎥
⎢⎣ 2 g ⎦⎥
where
Vd = velocity of the outfall channel (m/s, ft/s),
V = velocity just inside the end of the culvert barrel (m/s, ft/s), and
g = gravitational acceleration (m/s2, ft/s2).
Unsubmerged Flow
There are two equations for unsubmerged (weir) flow. The first is based on the
specific head at critical depth (with correction factors), whereas the second
equation is more closely related to a weir equation. Although either equation
will produce adequate results, the second equation is more commonly used
during hand calculations because it is easier to apply.
Note that these equations were developed by the U.S. FHWA and as such are
only intended for use in U.S. customary units. SI units should be converted
before applying these equations.
(A)
HW Water surface
W.S.
(B)
Water surface
HW W.S.
Inlet unsubmerged
outlet submerged
(C)
Water surface
HW
W.S.
Inlet submerged
outlet unsubmerged
(D)
Median drain
HW W.S.
Water surface
Unsubmerged, Form 1
HWi H ⎡ Q ⎤M
= c +K ⎢ AD0.5 ⎥ − 0.5S
D D ⎢⎣ ⎥⎦
Unsubmerged, Form 2
HWi ⎡ Q ⎤M
=K⎢ ⎥
D ⎢⎣ AD0.5 ⎥⎦
where
HWi = headwater depth above the control section invert (ft);
D = interior height of the culvert barrel (ft);
Hc = specific head at critical depth, yc + Vc2/2g (ft), where yc = critical depth
(ft) and Vc = velocity at critical depth (ft/s);
Q = culvert discharge (ft3/s);
A = full cross-sectional area of the culvert barrel (ft2);
S = culvert barrel slope (ft/ft); and
K, M = constants from Table 4.2.
These equations are applicable up to approximately Q/AD0.5 = 3.5. When
using the first equation with mitered inlets, use a slope correction factor (last
term of Form 1 equation) of +0.7S instead of −0.5S.
Submerged Flow
The equation for submerged (orifice) flow is:
HWi ⎡ Q ⎤2
= c⎢ ⎥ + Y − 0.5S
D ⎢⎣ AD0.5 ⎥⎦
where
HWi = headwater depth above the control section invert (ft),
D = interior height of the culvert barrel (ft),
Hc = specific head at critical depth, yc + Vc2/2g (ft),
Q = culvert discharge (ft3/s),
A = full cross-sectional area of the culvert barrel (ft2),
S = culvert barrel slope (ft/ft), and
c, Y = constants from Table 4.2.
This equation for submerged flow is applicable above approximately Q/AD0.5 =
4.0. When using this equation with mitered inlets, use a slope correction factor
(last term of equation) of +0.7S instead of −0.5S.
4.5 Problems
Solve the following problems using CulvertMaster’s Quick Culvert Calculator.
1. A culvert is 11-m long and has upstream and downstream inverts of
263.4 and 263.1 m, respectively. The downstream tailwater elevation is
below the downstream pipe invert.
For a Ke of 0.5 and a Manning’s n value of 0.013, what minimum diam-
eter concrete circular culvert (in mm) is required to pass 1.4 m3/s under a
roadway with a maximum allowable headwater elevation of 265.2 m?
What is the headwater elevation for the selected culvert?
6,970
72
68
Outlet
74
structure
64
62
66
0
66
6,96
68
Outflow
Inflow hydrograph
hydrograph
Spillway
6,9
72 70
Storage
facility can be designed. Despite the wide variation in design practices, the
underlying philosophy of most regulations is that peak runoff rates discharged
from the undeveloped site for one or more design storm events should not be
exceeded after the site development.
Several basic categories of detention facilities exist. Ponds may be wet or dry
between storm events; detention may occur on the surface or underground;
and a facility may lie in a main drainage path or be removed from it. The type
of facility may be dictated by regulations or by the needs and limitations of the
individual site or drainage area.
caused by land development. Although the peak flow rate from a site for a
given storm event may be the same as it was prior to development, the duration
of the discharge can be substantially longer.
The criteria described previously must be met for one or more specified design
storm events. For most localities, these design events will be specified in their
regulations. For example, the regulations may specify that postdevelopment
peak flows be less than predevelopment peaks for 24-h, synthetic storms hav-
ing return frequencies of 2, 5, and 10 years. (It is advisable, however, to check
the design for storms of different durations.) The regulations may further
state that an emergency overflow structure capable of handling flows from the
100-year storm be constructed. The magnitude of the largest storm that the
facility must be designed to accommodate sometimes depends on factors such
as site area and land use.
Inflow
hydrograph
Attenuation
Discharge
Change in
storage volume
Outflow
hydrograph
Lag
Time
event. If this volume is greater than the storage capacity of the pond, the pond
will overflow.
The following three scenarios summarize the relationship of a pond’s inflow
rate, outflow rate, storage, and water surface elevation for any given time inter-
val within a stormwater runoff event:
• If the average inflow rate is greater than average outflow rate during a time
interval, the volume of stored water in the pond increases during the time
interval, and the water surface elevation in the pond increases.
• If the average inflow rate is equal to the average outflow rate during a time
interval, the volume of stored water does not change during the interval,
and the water surface elevation remains constant.
• If the average outflow rate is greater than the average inflow rate during a
time interval, the volume of stored water decreases during the time interval,
and the water surface elevation in the pond decreases.
Most detention ponds have uncontrolled outlets (i.e., their outlets have no
valves or gates that can be opened or closed). For a pond with an uncontrolled
outlet into a free outfall, the peak of the outflow hydrograph from the pond
will occur at the point where the outflow hydrograph intersects the receding
limb of the inflow hydrograph (see Figure 5.2).
This concept can be used as a check on the validity of a set of routing calcu-
lations and can be reasoned as follows. Prior to the time of the intersection
point, the inflow rate is larger than the outflow rate and the volume of water
in the pond is increasing. After that intersection time, the volume is decreasing
because the outflow rate is greater than the inflow rate. Thus the maximum
stored volume of water in the pond occurs at the time where the two hydro-
graphs intersect. Because the outflow rate from a pond increases as the head-
water depth in the pond increases, it follows that the maximum outflow rate
must occur at the maximum depth, which also corresponds to the maximum
storage during the event. Finally, at the point of intersection, the outflow is
equal to the inflow rate, resulting in a zero rate of change of stored volume at
the moment in time when the water surface crests in the pond.
Dam Embankments
In many cases, detention facilities are formed by excavating and grading a
depression in the land surface, and the resulting pond has no dam embank-
ment per se. In other cases, a stream or localized depression is dammed to
create a detention pond. Most typically, a combination of these two approaches
is used in which an existing depression is enlarged and dammed. The design
of the dam should conform to accepted geotechnical engineering practices.
Governmental dam safety and permitting requirements may be applicable if
the embankment height or storage exceeds a regulatory limit.
Freeboard
Freeboard for a storm event is defined as the vertical distance from the maxi-
mum water surface elevation (stage) to the top of the dam. The storage availa-
ble within the freeboard elevation range adds a “safety” volume to the required
detention storage volumes that exist below this elevation. Some jurisdictions
and reviewing agencies mandate freeboard requirements for the maximum
design event.
Outlet Structures
The detention pond outlet structure allows flow to discharge from the pond at
a controlled rate. Early detention ponds were usually designed to control only a
single runoff event having a specified recurrence interval; thus these ponds had
very simple outlet structures. This practice provided little to no benefit with
smaller or larger storm events.
Modern detention design typically focuses on control of multiple storm events
and often requires several outlet openings at different pond stages. These mul-
tiple openings can consist of multiple pipe outlets with differing diameters or
invert elevations, orifices at various levels, overflow weirs, or any combination
of these elements. Figure 5.3 depicts a profile view of an outlet structure con-
figuration consisting of several structural components.
Culvert
Pond invert
Overflow Spillways
All detention ponds, regardless of the type, should be fitted with overflow
spillway facilities to convey flows in excess of those for which the detention
facility was designed. The overflow spillway may be an extension of a nor-
mal weir outlet, or it may be completely separate. To prevent undue damages
downstream, consideration must be given to the paths that pond overflows will
take. These outflow pathways are part of the major drainage system referred to
in previous chapters.
Average-End-Area Method
The average-end-area method uses the same concept frequently seen in
computations for estimating earthwork volumes. In this method, the incre-
mental storage volume available between two water surface elevations is
computed as:
70
68
Elevation, ft
66
64
h2
h1 62
Positive change in volume causes rise
in water surface.
60
58
0 1 2 3 4 5 6 7
ΔS
Volume, ac ft
S1 S2
A 1+ A 2
ΔS = (h 2− h 1 )
2
where
h1 = water surface elevation 1
h2 = water surface elevation 2 (h2 > h1)
A1 = available water surface area corresponding to stage h1
A2 = available water surface area corresponding to stage h2
ΔS = ncremental storage volume available between h1 and h2.
Computation of incremental storage amounts for a series of known stages and
surface areas, followed by summing of the incremental storages, yields the
desired relationship between h and S.
One limitation of the average-end-area method is that it applies a linear aver-
aging technique to describe area, whereas the change in area is actually a sec-
ond-order function of height for a pond with sloped sides. However, the error
associated with numerical volume integration is well within the error bounds of
other key design assumptions, such as the selected design storm and watershed
homogeneity.
Solution
The difference in elevation along the length of the floor of the vault is
75(0.02) = 1.50 ft
Because the floor elevation at its lower end is at an elevation (stage) of 123.00 ft,
the elevation at its high end is 124.50 ft. For any stage between these elevations,
the rectangular water surface area is equal to the width of the vault (W = 40 ft)
times the length of the water surface area. That length depends on the depth of
water in the vault and can be expressed in terms of the elevation, h, as
(h − 123.00)
L=
0.02
40(h − 123.00)
As = WL =
0.02
These expressions apply to values of h ranging from 123.00 to 124.50 ft. For
stages above 124.50 ft, the surface area is constant because of the vertical walls
and is equal to
The first two columns of Table 5.1 give the elevation versus surface area rela-
tionship for elevations from 123.00 to 130.00 ft.
To compute the elevation–storage relationship using the average-end-area
method, adjacent area values in Column 2 are added and then divided by the
stage increment for those areas. The steps are described next, and results are
reported in Table 5.1. The relationship is plotted in Figure 5.5.
• Column 3 = [Column 2 (current row)]+[Column 2 (previous row)]
• Column 4 = [Column 1 (current row)]−[Column 1 (previous row)]
• Column 5 = [Column 3 (current row)/2]×[Column 4 (current row)]
• Column 6 = [Column 5 (current row)]+[Column 6 (previous row)]
131
130
129
128
Elevation, ft
127
126
125
124
123
122
0 5,000 10,000 15,000 20,000
Storage, cu ft
Conic Method
This method applies a conic geometry to each stage increment to approximate
the nonlinear relationship between stage and water surface area, as illustrated
in Figure 5.6. The incremental storage volume between two stages h1 and h2
(h2>h1) is defined as
⎛ h 2− h 1 ⎞⎟
ΔS = ⎜⎜ ⎟( A + A2 + A 1 A 2 )
⎜⎝ 3 ⎟⎠ 1
Solution
Table 5.3 is used to determine the cumulative volume at each elevation. Each
column represents a step in the calculations used for the conic method. A plot
of the elevation–storage relationship is given in Figure 5.7.
• Columns 1 and 2 are the elevation versus surface area relationship and are
taken from Table 5.2.
• Column 3 = {[Column 1 (current row)] − [Column 1 (previous row)]}/3.
Area A2
h2
h2–h1
Area A1
h1
108
107
106
105
Elevation ,ft
104
103
102
101
100
99
0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0
Storage, ac ft
70
68
66
Elevation, ft
64
h2
h1 62
Rise in water surface
increases outflow
60
58
0 20 40 60 80 100 120 140 160 180 200
Q1 Q2 Discharge, cfs
Overflow Spillways
An overflow spillway has the purpose of safely conveying discharges when the
design storage capacity of a detention pond is exceeded during large storm
events. Overflow spillways are commonly referred to as emergency spillways
because they are designed to operate under extreme rainfall conditions to pre-
vent dam overtopping.
The storm recurrence interval for which an overflow spillway should be
designed is an issue that requires exercise of professional judgment, as well
as consideration of the adverse effects of downstream flooding or damages.
In areas where the major drainage system is evaluated for a designated storm
recurrence interval, the same recurrence interval might be used for overflow
spillway design. On the other hand, if failure of a dam embankment could
cause considerable damages or loss of life downstream, regulatory requirements
pertaining to dams may be appropriate.
Detailed coverage of spillway hydraulics is beyond the scope of this book. For
more advanced discussion on spillway and outlet design, refer to Stormwater
Conveyance Modeling and Design (Haestad Methods and R. Durrans, 2003).
Structures in Parallel
The most basic type of composite outlet structure is one that has different
outlet devices working independently of one another. In this situation, the
hydraulic rating characteristics for each device are simply summed at each ele-
vation to compute the total outflow. Figure 5.9 displays a composite set of
outlet structures working in parallel, and Figure 5.10 shows the rating curve
resulting from this outlet configuration.
A
Top of dam
Weir spillway
Culvert
A
Front view Section A–A
10
way
spill
rt +
C ulve
Spillway crest
Pond depth, ft
5
rt
lve
Cu
0
0 50 100 150 200
Discharge, cfs
Figure 5.10: Depth versus discharge for a composite outlet with outlet
structures discharging in parallel.
Calculate and graph the elevation versus discharge rating curve from elevations
100.00 to 107.00 ft for each 0.50-ft increment. (Note that a smaller increment
would yield a smoother, more precise rating curve.)
Solution
Table 5.4 shows the stage–storage calculation results described by the following
steps. Figure 5.11 is a plot of the composite structure rating curve.
• Orifice centroid elevation = 100.00 ft + (1/2 × 0.5 ft) = 100.25 ft.
• Orifice area = πD2/4 = 0.196 ft2.
• Column 1 is the range of rating table elevations (100.0 –107.0 ft, in
0.5-ft increments).
• Column 2 is the head on the orifice = (Column 1)−(orifice centroid
elevation).
• Column 3 is the orifice discharge = solution to orifice equation for C =
0.60, A = 0.196, h = Column 2 (the orifice equation is Q = CA(2gH )0.5,
as introduced in Chapter 1.).
• Column 4 is the head on the weir = (Column 1)−(weir crest elevation).
• Column 5 is the weir discharge = solution to rectangular weir equation for
C = 2.6, L = 15.0, h = Column 4 (the weir equation is Q = CLH 3/2, as
introduced in Chapter 1.).
• Column 6 is the composite structure discharge = (Column 3) + (Column 5).
108
107
106
Elevation, ft
105
104
103
102
101
100
99
0 20 40 60 80 100 120
Flow, cfs
Figure 5.11: Rating curve for outlet structure in Example 5.3.
Top of dam
A
Weir spillway
Culvert
Culvert
Structures in Series
When outlet structures are in series, it is possible that backwater effects caused
by hydraulic controls at downstream structures will affect the discharges
through connecting upstream structures. In these cases, hydraulic grade line
calculations are required to determine the interactions between individual
structure components.
Figure 5.12 illustrates a pond outlet with structures in series. In this example,
as the stage rises, the culvert eventually throttles (controls) the combined flow
through the v-notch weir and the inlet box, causing the inlet box riser to be
fully submerged. This phenomenon is shown graphically near the center por-
tion of the stage versus discharge relationship in Figure 5.13.
This text provides only an introduction to basic outlet hydraulics. For more
detailed explanations on analyzing complex outlet structures such as outlets in
series, refer to Stormwater Conveyance Modeling and Design (Haestad Methods
and R. Durrans, 2003).
10
way
spill
ert +
Culv
Spillway crest
rt
V-no
lve
Cu
0
0 50 100 150 200
Flow, cfs
Figure 5.13: Rating curve for a composite outlet with outlet structures
discharging in parallel and in series.
sections) and on the initial water surface elevation in the pond when the inflow
hydrograph begins.
Fundamentally, routing of a pond amounts to a solution of the conservation
of mass equation:
dS
= I (t ) − O(t )
dt
where
S = the volume of water in the pond (ft3, m3);
I(t) = specified inflow at time t (cfs, m3/s); and
O(t) = the outflow at time t (cfs, m3/s).
This equation states that the change in storage during a time period is equal to
the difference between the inflow and the outflow. There are two unknowns in
this equation (S and O); thus information on both storage characteristics and
hydraulic characteristics is required to solve it. Various routing methods are
available, but the most commonly used method in simple stormwater deten-
tion is the storage indication method, which is based on a finite difference
approximation of the conservation of mass equation.
Figure 5.14 shows the inflow and outflow hydrographs for a pond. The graph
magnifies a time slice (step) for which the inflow is greater than the outflow. The
x-axis of the slice has units of time (t), and the y-axis of the slice, which represents
Δ Volume
Inflow Storage
Pond crests
Flow
Outflow
Δ Time Time, h
flow rate, has units of volume (L3) per time, or L3/t. Thus the bounded area of
the slice has units of T × (L3/t) = L3, or volume. Because the inflow is greater
than the outflow for this time slice, the change in volume is positive, which
means the pond water surface (stage) is rising during this time step.
At the point where the inflow equals the outflow, the change in storage is zero,
the water surface elevation (stage) is constant, and the pond is cresting. After the
pond has crested, the outflow rate is greater than the inflow rate, resulting in a net
loss in volume that causes the water surface (stage) to fall during the time step.
Because the shaded area of Figure 5.14 represents the bounded area for which
inflow is greater than outflow before the pond crests, this area represents the
total volume stored while the pond is filling. As stated previously, the shaded
area represents the maximum storage volume.
Figure 5.15 shows a portion of the inflow and outflow hydrographs for a pond.
The graph illustrates a time period in which the inflow rate is greater than
the outflow rate, but the relationships that follow are applicable for any time
period, regardless of the relative magnitudes of the inflow and outflow rates.
Over a time interval Δt = tn+1 − tn, the inflow rate varies from In to In+1, and
In+1 Inflow
In
Discharge, cfs
Outflow
On+1
On
Dt
tn Time tn+1
therefore the average inflow rate over the time interval can be approximated
as (In + In + 1)/2. Over the same time interval, the average outflow rate can be
approximated as (On + On + 1)/2. The change in storage over the time interval
is ΔS = Sn + 1 − Sn, whereas the time rate of change of storage may be approx-
imated as ΔS/Δt = (Sn + 1 − Sn)/Δt. Substitution of the average inflow and
outflow rates and the approximate rate of change of storage into the conserva-
tion of mass equation yields
Sn +1 − Sn I + In +1 On + On +1
= n −
Δt 2 2
⎛ 2Sn +1 ⎞⎟ ⎛ 2Sn ⎞⎟
⎜⎜ ⎜
⎜⎝ Δt + On +1 ⎟⎟⎠ = ⎜⎜⎝ Δt − On ⎟⎟⎠ + ( In + In +1 )
Thus knowledge of the quantities Sn, On, and In at the beginning of a time
interval (at time tn), as well as knowledge of the inflow rate In + 1 at the end
of the time interval, can be used to compute the value of the left side of the
equation. However, there are two unknowns on the left side (On + 1 and Sn + 1).
Additional information must be established to expose the relationship between
On + 1 and Sn + 1 so that these two values can be determined.
A storage indication curve is a graph delineating the relationship between O
and (2S/Δt + O). Such a graph can be constructed by realizing that for any
stage h in a pond, there are corresponding values of O and S; hence, there are
corresponding values of O and (2S/Δt + O).
In performing storage indication calculations for a specified time interval n+1,
the storage indication equation is applied to compute the numerical value of
the quantity (2Sn + 1/Δt + On + 1). Next, the value of (2Sn + 1/Δt − On + 1)
for the time interval n + 1 must be found. Because the value of (2Sn + 1/Δt +
On + 1) is known, and this value can be used with the storage indication curve
to find On+1, (2Sn + 1/Δt − On + 1) can be computed as
⎛ 2Sn +1 ⎞⎟ ⎛ 2Sn +1 ⎞⎟
⎜⎜ ⎜
⎜⎝ Δt − On +1 ⎟⎟⎠ = ⎜⎜⎝ Δt + On +1 ⎟⎟⎠ − 2On +1
The term on the left side of the equation can now be used in the next time
step (n + 2) to compute (2Sn + 2/Δt + On + 2), and the process repeats through
every time step.
Section 5.7 describes the full procedure used in storage indication routing and
presents an example problem.
Solution
The results from Examples 5.2 and 5.3 are combined in Table 5.6. Note the
importance of calculating the same elevations for storage and discharge in
order to directly correlate flow and storage at each elevation.
Follow the five steps outlined for storage indication pond routing:
1. Compute the postdeveloped pond inflow hydrograph using the proce-
dures presented in Chapter 2 or another method. In this example, the
final inflow hydrograph is given.
2. Calculate the stage versus storage curve from the pond grading plan or
storage structure geometry. This information was calculated in Example
5.2 (no calculations are necessary for this example).
3. Calculate the stage versus outflow curve for the pond outlet structure.
This information was calculated in Example 5.3.
4. Calculate and graph the storage indication curve that correlates 2S/Δt + O
and outflow. To calculate this curve, the rating elevations for the stage
versus storage curve must be identical to those used for the stage versus
outflow curve. Use a 10-min time step in this example.
Step 4 results for 2S/Δt + O are given in Table 5.7. Values given in
Column 4 were computed using
Table 5.7: Relationship between Elevation and 2S/Δt + O for Example 5.4
400
350
300
(2S/t) + O
250
200
150
100
50
0
0 5 10 15 20
Flow, cfs
25
20
Inflow
Flow, cfs
15
10
Outflow
5
0
0 50 100 150 200 250 300 350
Time, min
What Is PondPack?
PondPack is a stormwater modeling program capable of analyzing a wide
range of situations, from simple sites to complex networked watersheds. The
program analyzes pre- and postdeveloped watershed conditions and sizes
ponds.
PondPack also computes outlet structure rating curves with tailwater effects,
pond infiltration, pond detention times, and channel routing effects. It
includes the capability to model interconnected ponds with divergent (mul-
tiple) outfalls.
The program can use any rainfall duration or distribution to compute hydro-
graphs for multiple events. Hydrographs can be added at junctions and routed
through multiple reaches and ponds. PondPack helps automate almost every
aspect of drainage design computations.
PondPack graphically displays such items as watershed diagrams, hydrographs,
rainfall curves, I–D–F curves, outlet rating curves, volume curves, time ver-
sus elevation curves, time versus volume curves, cross sections, channel rating
curves, and a wide variety of other output diagrams.
5.9 Problems
The following problems can be solved with the PondPack computer program.
1. A detention pond has the shape of a trapezoidal basin. The specifications are:
Construct an elevation versus volume graph for the basin. What is the volume
at an elevation of 120 ft?
2. Given the elevation-area data for a detention pond, construct an elevation
versus volume graph for the pond. What is the volume of the pond at an
elevation of 205 m?
3. An outlet structure for a pond consists of a riser pipe with two orifices.
The riser is a 3-ft diameter pipe with a rim elevation of 105 ft. A 1-ft
diameter orifice has an invert of 101 ft and a 0.5-ft orifice has an invert
of 100 ft. Use 3.33 for the weir coefficient and 0.614 for the orifice coef-
ficient. Construct an elevation versus flow diagram for a headwater range
of 100–120 ft, with an increment of 1 ft. What is the discharge at a head-
water elevation of 120 ft?
4. A composite outlet has a 4-m broad crested weir at an elevation of
203.5 m and a 50-cm diameter orifice at an elevation of 200 m. The weir
coefficient is 1.84 and the orifice coefficient is 0.614. Construct an ele-
vation versus discharge curve for the structure over a headwater range
of 200–205 m with a step of 1 m. What is the discharge at a headwater
elevation of 205 m?
5. Use the detention pond described in Problem 5.2 and the outlet struc-
ture described in Problem 5.4 to route the following hydrograph. What
is the peak discharge from the pond? What is the maximum elevation in
the pond?
Water Demands
Just as storm sewer analysis is driven by the watershed runoff flow rate, water
distribution system analysis is driven by customer demand. Water usage rates
and patterns vary greatly from system to system and are highly dependent on
climate, culture, and local industry. Every system is different, so the best source
of information for estimating demands is directly recorded system data.
Metered Demand
Metered demands are often a modeler’s best tool and can be used to calculate
average demands, minimum demands, peak demands, and so forth. This data
can also be compiled into daily, weekly, monthly, and annual reports that show
how the demands are influenced by weather, special events, and other factors.
Unfortunately, many systems still do not have complete system metering. For
these systems, the modeler is often forced to use other estimation tools (including
good engineering judgment) to obtain realistic demands.
Demand Patterns
A pattern is a function relating water use to time of day. Patterns allow the
user to apply automatic time-variable changes within the system. Different
categories of users, such as residential or industrial customers, will typically be
assigned different patterns to accurately reflect their particular demand varia-
tions. A diurnal curve is a type of pattern that describes changes in demand
over the course of a daily cycle, reflecting times when people are using more or
less water than average. Most patterns are based on a multiplication factor ver-
sus time relationship, whereby a multiplication factor of 1.0 represents the base
value (often the average value). In equation form, this relationship is written as:
Q t = A t× Qbase
where
Qt = demand at time t,
At = multiplier for time t, and
Qbase = baseline demand.
Using a representative diurnal curve for a residence (Figure 6.1), we see that
there is a peak in the diurnal curve in the morning as people take showers and
prepare breakfast, another slight peak around noon, and a third peak in the
evening as people arrive home from work and prepare dinner. Throughout the
night, the pattern reflects the relative inactivity of the system, with very low
flows compared to the average. (Note that this curve is conceptual and should
not be construed as representative of any particular network.)
There are two basic forms for representing a pattern: stepwise and continuous.
A stepwise pattern is one that assumes a constant level of usage over a period
of time and then jumps instantaneously to another level where it again remains
2.0
Continuous
Stepwise
Multiplication factor
1.5
Average
1.0
0.5
0.0
0 6 12 18 24
Time of day
steady until the next jump. A continuous pattern is one for which several
points in the pattern are known and sections in between are transitional, result-
ing in a smoother pattern. Notice that, for the continuous pattern in Figure 6.1,
the magnitude and slope of the pattern at the start and end times are the
same—a continuity that is recommended for patterns that repeat.
Because of the finite time steps used in the calculations, most computer pro-
grams convert continuous patterns into stepwise patterns for use by the algo-
rithms, with the duration of each step equal to the time step of the analysis.
Qleaks = K ∑ i =1 Pin
N
where
N = total number of demand nodes,
Pi = Pressure head at node i, and
n = exponent.
Proportional methods—There are other various methods that use an emit-
ter coefficient, the node pressure head, and proportional values for the pipe
length, and metered consumption. These methods use weighted factors for
pipe length, and metered consumption in addition to the pressure-dependent
leakage flow rate at each node.
Several hydraulic components that are unique to pressure piping systems, such
as regulating valves and pumps, add complexity to the analysis.
Minor Losses
Localized areas of increased turbulence cause energy losses within a pipe, creat-
ing a drop in the energy and hydraulic grades at that point in the system. These
disruptions are often caused by valves, meters, or fittings (such as the pipe
entrance in Figure 6.2), and are generally called minor losses. These minor
losses are often negligible relative to friction losses and may be ignored during
analysis.
Although the term “minor” is a reasonable generalization for most large-scale
water distribution models, these losses may not always be as minor as the name
implies. In piping systems that contain numerous fittings relative to the total
length of pipe, such as heating or cooling systems, the minor losses may actu-
ally have a significant impact on the energy loss.
The equation most commonly used for determining the loss in a fitting, valve,
meter, or other localized component is:
V2
Hm = K
2g
where
Hm = minor loss (m, ft),
K = minor loss coefficient for the specific fitting,
V = velocity (m/s, ft/s), and
g = gravitational acceleration (m/s2, ft/s2).
Typical values for the fitting loss coefficient are included in Table 6.1. As can be
seen with similar fitting types, the K-value is highly dependent on bend radius,
contraction ratios, and so forth. Gradual transitions create smoother flow lines
and smaller head losses than sharp transitions because of the increased turbu-
lence and eddies that form near a sharp change in the flow pattern. Figure 6.2
shows flow lines for a pipe entrance with and without rounding.
Pum
pc
urv
e
Head
Operating point
urve
em c
Syst HL
Hs
Flow rate
Variable-Speed Pumps
A centrifugal pump’s characteristic curve is fixed for a given motor speed and
impeller diameter but can be determined for any speed and any diameter by
applying the affinity laws. For variable-speed pumps, these affinity laws are
presented as:
2
Q1 n1 H1 ⎛⎜ n1 ⎞⎟
= and = ⎜ ⎟⎟
Q2 n2 H2 ⎜⎝ n2 ⎟⎠
180
160 n = 0.95
n = 0.90
140
n = 0.85
120 n = 0.80
Head, ft 100
80
60
Best
efficiency
40
point
20 Speed
n=
Full speed
0
0 100 200 300 400
Flow, gpm
where
Q = pump flow rate (m3/s, ft3/s),
H = pump head (m, ft), and
n = pump speed (rpm).
Thus pump discharge rate is proportional to pump speed, and the pump dis-
charge head is proportional to the square of the speed. Using this relationship,
once the pump curve is known, the curve at another speed can be predicted.
Figure 6.4 illustrates the affinity laws applied to a variable-speed pump. The
line labeled “Best efficiency point” indicates how the best efficiency point
changes at various speeds.
Path from A to B
Conservation of Energy
Chapter 1 introduced the application of the energy equation to hydraulic
analysis. The principle of conservation of energy dictates that the head losses
through the system must balance at each point (Figure 6.5). For pressure net-
works, this means that the total head loss between any two nodes in the system
must be the same regardless of the path taken between the two points. The
head loss must be “sign consistent” with the assumed flow direction (that is,
head loss occurs in the direction of flow, and head gain occurs in the direction
opposite that of the flow).
Although the equality can become more complicated with minor losses
and controlling valves, the same basic principle can be applied to any path
between two points. As shown in Figure 6.5, the combined head loss around
a loop must equal zero in order to compute the same hydraulic grade for a
given point.
Extended-Period Simulation
An extended-period simulation is used to determine the behavior of the
system over time. This type of analysis allows the user to model tanks fill-
ing and draining, regulating valves opening and closing, and pressures
and flow rates changing throughout the system in response to varying
demand conditions and automatic control strategies formulated by the
modeler.
While a steady-state model may tell the user whether the system has the capa-
bility to meet a specific demand, an extended-period simulation indicates
whether the system has the ability to provide acceptable levels of service over
a period of minutes, hours, or days. Extended-period simulations can also
be used for energy consumption and cost studies, as well as for water quality
modeling.
Data requirements for an extended-period simulation go beyond what is
needed for a steady-state analysis. The user must determine water usage pat-
terns, provide more detailed tank information, and enter operational rules for
pumps and valves.
ΣQiCi
Ca =
ΣQi
where
Ca = average (mixed) constituent concentration,
Q i = inflow rates, and
Ci = constituent concentrations of the inflows.
Age
Water age provides a general indication of the overall water quality at any given
point in the system. Age is typically measured from the time that the water
enters the system from a tank or reservoir until it reaches a junction.
x
Aj = Aj −1 +
V
where
Aj = age of water at jth node,
x = distance from node j − 1 to node j, and
V = velocity from node j − 1 to node j.
If there are several paths for water to travel to the jth node, the water age is
computed as a weighted average using the equation:
AAj =
∑Q i ⎡⎣⎢ AA i + ( x / V )i ⎤⎦⎥
∑Q i
where
AAj = average age at the node immediately upstream of node j and
Qi = flow rate to the jth node from the ith node.
Trace
Identifying the origin of flow at a point in the system is referred to as flow
tracking or trace modeling. In systems that receive water from more than one
source, trace studies can be used to determine the percentage of flow from each
source at each point in the system. These studies can be very useful in delineat-
ing the area influenced by an individual source, observing the degree of mixing
of water from several sources, and viewing changes in origins over time.
Constituents
Reactions can occur within pipes that cause the concentration of substances to
change as the water travels through the system. Based on conservation of mass
for a substance within a link (for extended-period simulations only):
∂c ∂c
=V + θ(c )
∂t ∂x
where
c = substance concentration as a function of distance and time,
t = time increment,
V = velocity,
x = distance along the link, and
θ (c) = substance rate of reaction within the link.
In some applications, there is an additional term for dispersion, but this term
is usually negligible (plug flow is assumed through the system).
Assuming that complete and instantaneous mixing occurs at all junction
nodes, additional equations can be written for each junction node with the
following conservation of mass equation:
Ck x =0 =
∑Q jC j x =L + QeCe
∑Q j + Qe
where
Ck = concentration at node k,
j = pipe flowing into node k,
L = length of pipe j,
Qj = flow in pipe j,
Cj = concentration in pipe j,
Qe = external source flow into node k, and
Ce = external source concentration into node k.
Once the hydraulic model has solved the network, the velocities and the mix-
ing at the nodes are known. Using this information, the water quality behavior
can be derived using a numerical method.
Initial Conditions
Just as a hydraulic simulation starts with some amount of water in each storage
tank, initial conditions must be set for a water age, trace, or constituent con-
centration analysis. These initial water quality conditions are usually unknown,
so the modeler must estimate these values from field data, a previous water
quality model, or some other source of information.
To overcome the problem of unknown initial conditions at the vast majority
of locations within the water distribution model, the duration of the analysis
must be long enough for the system to reach equilibrium conditions. Note that
a constant value does not have to be reached for equilibrium to be achieved;
rather, equilibrium conditions are reached when a repeating pattern in age,
trace, or constituent concentration is established.
Pipes usually reach equilibrium conditions in a short time, but storage tanks are
much slower to show a repeating pattern. For this reason, extra care must be taken
when setting a tank’s initial conditions, in order to ensure the model’s accuracy.
Numerical Methods
There are several theoretical approaches available for solving water qual-
ity models. These methods can generally be grouped as either Eulerian or
Lagrangian in nature, depending on the volumetric control approach that is
taken. Eulerian models divide the system into fixed pipe segments, and then
track the changes that occur as water flows through these segments. Lagrangian
models also break the system into control volumes but then track these water
volumes as they travel through the system. This chapter presents two alterna-
tive approaches for performing water quality constituent analyses.
Time-Driven Method
The Time-Driven Method (TDM) is an example of a Lagrangian approach
(Figure 6.7). This method also breaks the system into segments, but rather
than using fixed control volumes as in Eulerian methods, the concentration
and size of water parcels are tracked as they travel through the pipes. With
each time step, the farthest upstream parcel of each pipe elongates as water
travels into the pipe, and the farthest downstream parcel shortens as water
exits the pipe.
Similar to the DVM, the reactions of a constituent within each parcel are cal-
culated, and the mass and flow entering each node are summed to determine
Original concentrations
Reacted concentrations
Original concentrations
Reacted concentrations
Model Calibration
Model calibration is the process of modifying parameters or values in a model
so it better matches what is happening in the real system. The calibration of
water distribution models is very complicated. There are many unknown val-
ues and parameters that are needed at any one time to reduce the discrepancy
between the model and real system. Many times the pipe roughness value is
adjusted to have the model results match the measured or expected values in
the real system. However, there are many other parameters that could influ-
ence the modeled results. The water demand at junctions and the status of
pipes and valves in the system could also be adjusted when calibrating a model.
Calibrating a model relies on accurate field measurement data. Field measure-
ments of pressures in the system, pipe flow rates, water levels in tanks, valve status,
and pump operating status and speed are all used to calibrate models. Critical to
all these measurements is the time for which the measurements are made. The
time of these measurements must all be synchronized to the time frame of the
model. In addition, because the conditions within a real system changes through-
out the day or year, field data should be collected for many different conditions
and times. The calibration process is used to adjust the model to simulate multiple
demand loading and operational boundary conditions. Only then can the mod-
eler be confident that the model is valid for many different conditions.
WaterGEMS uses Darwin Calibrator to assist in optimizing the model to
match field measurement data. The Darwin Calibrator allows the modeler
to input field data then request the software to determine the optimal solu-
tion of pipe roughness values, junction demands, or status (on/off). Pipes
that have the same hydraulic characteristics where one roughness values is
assigned to all pipes can be grouped together. Junctions can also be grouped
based on the demand pattern and location. Caution should be used group-
ing pipes and junctions since this could greatly affect the model calibration
accuracy.
System Design
The goal of water distribution system design is to maximize the benefits of
the system while minimizing the cost. The optimal solution is a design that
meets all the needs of the system at minimal cost. Some planning is needed to
account for additional future needs of the system, including potential growth
of the system in terms of demand and its location. The modeler must work
with the system owner and planning groups to account for both the current
and future needs.
WaterGEMS uses Darwin Designer to assist engineers to plan and design water
distribution networks. The Darwin Designer can be used to size new pipe and/
or rehabilitate old pipes to minimize cost, maximize benefit, or create a sce-
nario to trade-off cost and benefit. The least cost optimization is used to deter-
mine the pipe material and size to satisfy the needed design requirements. The
maximum benefit optimization is used to determine the most beneficial solu-
tion based on a known budget. The Darwin Designer will generate a number
of solutions that meet the design requirements at minimal cost or maximum
benefit. In either case, the best solution of new pipe or rehabilitation of old
pipe will be based upon the following input hydraulic criteria:
• minimum and maximum allowable pressures;
• minimum and maximum allowable pipe flow velocity;
• additional demand requirements; and
• pipe, pump, tank, valve, etc. status change requirements.
Critical to creating an accurate designed system is time and peak demand require-
ments. The peak demand and fire flow conditions are used to size pipes since the
pipe network must work for all conditions. Using average demand values to size
pipe without accurately accounting for peaking factors can create networks that
are either undersized and will not deliver the required water needs, or oversized
and much more expensive than need be. The daily and seasonal variations can
also greatly affect the final design. Demand variations need to be synchronized in
the model to accurately reflect what could happen in the real system.
6.10 Problems
Solve the following problems using the WaterGEMS computer program.
1. The ductile iron pipe network shown below carries water at 20°C.
Assume that the junctions all have an elevation of 0 m and the reservoir
P-
P-4
2
P-3
J-4 J-3
2. A pressure gage reading of 288 kPa was taken at J-5 in the pipe network
shown below. Assuming a reservoir elevation of 100 m, find the appro-
priate Darcy–Weisbach roughness height (to the hundredths place) to
bring the model into agreement with these field records. Use the same
roughness value for all pipes.
P-3
P-6
J-1 J-2
Pipe Diameter (mm) Length (m) Junction Elevation (m) Demand (l/min)
P-1 250 1,525 J-1 55 950
P-2 150 300 J-2 49 1,060
P-3 150 240 J-3 58 1,440
P-4 150 275 J-4 46 1,175
P-5 150 245 J-5 44 980
P-6 200 230
P-15
PMP-1 P-2
J-1 P-3 J-2 P-4 J-3
Hotel
P-11
P-12
J-11
P-5
10 J-9
P-
P-13 P-14
J-8
P- P-7 P-6
9 J-5
8 J-6 J-4
P-
J-7 J-10
Entrance Administration
and Maintenance
Pipe network: The pipe network consists of the pipes listed in the following
tables. The diameters shown are based on the preliminary design and may not
be adequate for the final design. For all pipes, use ductile iron as the material
and a Hazen–Williams C-factor of 130.
(a) Fill in or reproduce the Results Summary table after each run to get
a feel for some of the key indicators during various scenarios.
(b) For the average day run, what is the pump discharge?
(c) If the pump has a best efficiency point at 300 gpm, what can you say
about its performance on an average day?
(d) For the peak hour run, the velocities are fairly low. Does this mean
you have oversized the pipes? Explain.
(e) For the minimum hour run, what was the highest pressure in the
system? Why would you expect the highest pressure to occur during
the minimum hour demand?
(f) Was the system (as currently designed) acceptable for the fire
flow case with the sprinkled building? On what did you base this
decision?
(g) Was the system (as currently designed) acceptable for the fire
flow case with all the flow provided by hose streams (no sprin-
klers)? If not, how would you modify the system so that it will
work?
4. A ductile iron pipe network (C = 130) is shown below. Use the Hazen–
Williams equation to calculate friction losses in the system. The junctions
and pump are at an elevation of 5 ft and all pipes are 6 in. in diameter.
(Note: Use a standard, three-point pump curve. The data for the pump,
junctions, and pipes are in the tables below.) The water surface of the
reservoir is at an elevation of 30 ft.
PMP-1
P-
P-4
2
P-3
J-4 J-3
(a) What are the resulting flows and velocities in the pipes?
(b) What are the resulting pressures at the junction nodes?
(c) Place a check valve on pipe P-3 such that the valve only allows flow
from J-3 to J-4. What happens to the flow in pipe P-3? Why does
this occur?
(d) When the check valve is placed on pipe P-3, what happens to the
pressures throughout the system?
(e) Remove the check valve on pipe P-3. Place a 6-in. flow control valve
node at an elevation of 5 ft on pipe P-3. The FCV should be set so
that it only allows a flow of 100 gpm from J-4 to J-3 (Hint: a check
valve is a pipe property). What is the resulting difference in flows in
the network? How are the pressures affected?
(f) Why doesn’t the pressure at J-1 change when the FCV is added?
(g) What happens if you increase the FCV’s allowable flow to
2,000 gpm? What happens if you reduce the allowable flow to zero?
5. A local country club has hired you to design a sprinkler system that will
water the greens of their nine-hole golf course. The system must be able to
water all nine holes at once. The water supply has a water surface elevation
of 10 ft. All pipes are PVC (C = 150, use the Hazen–Williams equation
to determine friction losses). Use a standard, three-point pump curve for
the pump, which is at an elevation of 5 ft. The flow at the sprinkler is
modeled using an emitter coefficient. The data for the junctions, pipes,
and pump curve are given in the tables that follow. The initial network
layout is shown below.
P-
5
P-9
J-6
J-3
J-7
P-7
P-4
R-1 P-1 PMP-1
J-10 J-5
P-8
6
P-
P- J-2
P-2
11
P-3
J-1
Emitter Coefficient
Junction Label Elevation (ft) Pipe Label Diameter (in.) Length (ft)
(gpm/psi0.5)
J-1 – 10 P-1 4 10
Hole 1 8 7 P-2 4 1,000
Hole 2 10 7 P-3 4 800
Hole 3 15 40 P-4 3 750
Hole 4 12 5 P-5 3 500
Hole 5 8 5 P-6 3 700
Hole 6 8 15 P-7 2 400
Hole 7 10 20 P-8 4 800
Hole 8 15 10 P-9 3 500
Hole 9 8 12 P-10 2 400
P-11 2 500
Elm St.
J-8
J-2
J-7
P-1
P-8
P-9
J-3 J-1
Street D
Street E
Town
P-6
Street C
P-7
P-5
Street B
h
ighwa
Common land
open
y #64
P-4
P-3 J-4
P-2 J-5
J-9 J-6 Street A
Junction Label Number of Lots Serviced Pipe Label Length (m) Diameter (mm)
J-1 5 P-1 60.0 150
J-2 4 P-2 60.0 150
J-3 4 P-3 110.5 150
J-4 5 P-4 164.0 150
J-5 6 P-5 152.5 150
J-6 6 P-6 204.0 100
J-7 6 P-7 148.0 150
P-8 61.0 100
P-9 194.0 150
Currently, a model of the entire water system does not exist. However, hydrant
tests were conducted using hydrants located on two water mains, one in Town
Highway 64 and the other in Elm Street. The following data were obtained:
Town Highway 64 Hydrant Test
Static pressure 310.3 kPa
Residual pressure 98.5 kPa at 32 l/s
Elevation of pressure gauge 190 m
Elm Street Hydrant Test
Static pressure 413.7 kPa
Residual pressure 319.3 kPa at 40 l/s
Elevation of pressure gauge 191.5 m
The subdivision will connect to existing system mains in these streets at nodes
J-8 and J-9. (Hint: Model the connection to an existing water main with a
reservoir and a pump.)
(a) What are the demands at each of the junction nodes? What is the
total demand?
(b) Does the present water distribution system have enough capacity to
supply the new subdivision?
(c) Which connection to the existing main is supplying more water to
the subdivision? Why?
(d) Are the proposed pipe sizes adequate to maintain velocities between
0.15 and 2.44 m/s, and pressures of at least 140 kPa?
(e) Would the subdivision have enough water if only one connection
were used? If so, which one?
(f) What do you think are some possible pitfalls of modeling two con-
nections to existing mains within the same system, as opposed to
modeling back to the water source?
7) Use the pipe sizes given in the table below for the subdivision in
Problem 6.
Pipe Information for Problem 7
R-1 PMP-1
P-7
P-
3 J-9
P-1
8
P-11
P-10
P-12
J-1
P-2
P-14
J-8
J-2 P-3 J-3
P-9
T-1
(a) Perform an age analysis on the system using a duration of 300 h and
a time step of 2 h. Fill in the results table, indicating the maximum
water age at each junction and tank after the system reaches equi-
librium (a pattern of average water age vs. time becomes evident).
What point in the system generally has the oldest water? Explain
why the water is oldest at this location.
(b) Perform a constituent analysis using the same duration and time step
as in part (a). Fill in the results table, indicating the minimum chlo-
rine concentration for each junction and tank after the system has
reached equilibrium (a pattern of concentration versus time becomes
evident). What point in the system has the lowest chlorine concen-
tration? Explain why the chlorine residual is lowest at this location.
(c) From the above table and graphs of demand, age, and concentration
versus time generated within WaterGEMS, determine the following
correlations:
(1) age and chlorine concentration
(2) demand and chlorine concentration at a junction
(3) demand and water age at a junction
(d) Why is it necessary to run the model for such a long time? Do you
feel that 300 h is too long or too short a time period for testing the
model? Why?
9. A planning commission has indicated a new industry may be connected
to the water system described in Problem 8. You are to determine the pipe
diameters in the network to minimize the installation cost assuming all
the pipes are ductile iron. Use the Darwin Designer to determine the total
cost and size each pipe for each of the following conditions. Use the pipe
cost information from Tutorial 4 for the ductile iron pipe.
(a) Size the pipes using a demand multiplier (peaking factor) of 3.2.
The pressure must remain between 170 and 550 kPa during peak
demand. Exclude pipes P-7 and P-8 in your analysis when deter-
mining the pipe sizes. Hint: you will need to specify an additional
demand of zero (0 l/min) with the default pressure constraints at the
junctions or a fatal error will occur.
(b) It is expected that a new industry with an expected additional
demand of 2,000 l/min with a required minimum pressure of
260 kPa will be added to the system. It could be tapped into the
network at either junction 6, 7, or 8. Size the pipes for the condi-
tions in part (a) above along with the industry added to all proposed
junctions. You will need to analyze the network three times, once
for the industry at J-6, again with the industry at J-7, then finally
with the industry at J-8.
a. Indicate which option(s) would work.
b. Which junction should the industry be tapped into to be the
least costly and what is the expected cost?
c. What is the size of each pipe for the best solution for the least
costly option with the industry added?
d. What is the calculated minimum pressure at the industry for
the best solution?
Transient Flow
Transient flow is the transition from an initial flow condition to the final flow
condition due to a change in the system. A simple illustration of this is when a
valve suddenly closes in a single pipe fed by a tank (Figure 7.1). When the valve
in the pipe instantaneously closes, the fluid velocity suddenly drops behind the
valve. A positive pressure wave is formed and travels up the pipe toward the tank.
100 m
Globe valve
Q (closes in 20 s)
a = = 330
0
1,0 l
00 /s
m/
s
300
t = 24.2 s
Packing
220
Head, m
140
ΔH0
60
–20
–100
0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0
Distance, m (×10,000)
300
t = 30.2 s Packing
220
Head, m
140
ΔH0
60
–20
–100
0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0
Distance, m (×10,000)
As the pressure wave in the pipe moves toward the tank, the velocity behind
the wave is zero and the velocity in front of the wave is approximately equal to
the initial velocity. When the pressure wave reaches the tank, water from the
pipe reenters the tank causing a small increase of the water level in the tank. At
that instant, the pressure wave reaches the tank and the pipe velocity is zero,
but the head in the pipe is higher than in the tank. This head imbalance causes
flow to reverse in the pipe and the pressure in the pipeline to drop, creating a
“pressure reduction wave” that moves back toward the valve. The pressure wave
will oscillate back and forth in the pipe until it is dampened by the system and
reaches the final steady-state condition (Figure 7.2).
p(t)
pmax
p2
p1
pmin
Valve closure
TT x
TM
t
Start of valve closure
On the downstream side of the valve, the water will initially move away from
the closed valve causing a region where the pressure is suddenly reduced. There
is the potential for the pressure to drop to a level where dissolved gasses can
come out of the water, and also below the vapor pressure of the fluid itself to
create vaporous cavitation (column separation). If the flow (liquid columns)
reverses, these pockets of gasses or vapor could suddenly collapse creating
a strong transient with the potential to damage the pipe lining, or rupture
the pipe.
Common Causes
The majority of hydraulic transients are the result of changes at system bound-
aries, typically at the upstream or downstream ends of the systems, or at high
points (locations of high elevation) in the pipe network. Three common causes
of transients in hydraulic systems are described below.
Pumps—A pump delivers energy from the pump impeller to the fluid as the
fluid passes from the suction side to the discharge side of the pump. Pumps
force water through the piping network and components on the discharge side
of the pump. When the pump starts, it can greatly increase the pressure rap-
idly, or when a pump sags or fails, a sudden drop on pressure can propagate
downstream.
Turbines—A turbine is located on the downstream end of a penstock (the
feeding conduit). As water flows through the penstock and enters the turbine,
the water energy rotates a turbine impeller, which creates an electrical current
from the generator. If electricity from the turbine is no longer needed, the
flow through the penstock is stopped, which can create a pressure increase that
propagates upstream into the penstock.
Valves—A valve can change, start, or stop flow suddenly. The rate of closing
or opening of the valve can allow for a gradual increase or decrease of flow
to reduce the magnitude of the transient pressure wave. However, if a valve
suddenly closes (valve slam), there is a sudden high-pressure pulse that prop-
agates upstream of the valve and a sudden low-pressure pulse that propagates
downstream. If the low-pressure pulse is severe enough, gas or vapor pockets
may form, or external water could be pulled into the pipe network. Also, there
is the potential for the pipe to rupture or collapse.
Sudden changes in flow (fluid momentum) create a pressure pulse that prop-
agates away from the disruption in every direction and throughout the entire
pipe network. If no other transient event occurs, than the pressure wave will
continue to propagate as unsteady flow until the transient is dampened by fric-
tion. After the transient is dampened, the system reaches the final steady-state
condition.
The derivation of the complete equations for transient analysis is beyond the
scope of this book but can be found in Almeida and Koelle (1992) and Wylie
and Streeter (1993). Also, for the remainder of this book, only the Elastic
Model will be used to describe hydraulic transient flow.
∂H a2 ∂Q
+ =0
∂t gA ∂x
where
∂H / ∂t = change in hydraulic head (m, ft) with time (s),
∂Q / ∂x = change in flow (m3/s, cfs) with position (m, ft),
a = pressure wave celerity in the fluid (m/s, ft/s),
A = cross-sectional area of the fluid (m2, ft2), and
g = gravity constant (m/s2, ft/s2).
The momentum equation for a fluid can be derived from considering the forces
acting upon a small fluid volume (control volume), including the shear stresses
generated by the fluid motion and viscosity. The three-dimensional momen-
tum equations for a fluid are also known as the Navier–Stokes equations. The
simplified general form of the equation for unsteady one-dimensional flow that
is used in transient flow models is as follows:
∂Q ∂H fQ Q
+ gA + =0
∂t ∂x 2 DA
where
∂Q / ∂t = change in flow (m3/s, cfs) with time (s),
∂H / ∂x = change in hydraulic head (m, ft) with position (m, ft),
Q = fluid flow rate (m3/s, cfs),
D = inside pipe diameter (m, ft), and
f = Darcy–Weisbach friction coefficient.
The last term in the momentum equation represents the friction losses in the
direction of fluid flow.
E f ρf
a=
1 + (DE f eE )ψ
where
a = pressure wave celerity (m/s, fps),
Ef = modulus of elasticity for the fluid (Pa, lbf/ft2),
ρf = fluid density (kg/m3, slugs/ft3),
D = inside pipe diameter (mm, in.),
e = pipe wall thickness (mm, in.),
E = modulus of elasticity for the pipe material (Pa, lbf/ft2), and
Ψ = pipeline support factor.
The pipeline support factor (Ψ) is determined from the Poisson’s ratio of the
pipe material. For thin-walled pipes, if the pipeline has functioning expand-
ing joints, Ψ = 1; if the pipeline anchored against axial (longitudinally) move-
ment, Ψ = 1−μ2; if the pipeline is anchored at one end only, Ψ = 5/4−μ
(see Advanced Water Distribution Modeling and Management. Heastad Press,
2003). Properties of pipe materials and common fluids can be found in
Tables 7.1 and 7.2.
In a single pipeline, the pressure wave will reach the other end of the pipe in a
time interval equal to L/a seconds where L in the pipe length (m, ft). For the
pressure wave to return to the original location that caused the transient, the time
interval will be 2L/a seconds. The quantity 2L/a is known as the characteristic
time for the pipeline. Flow control operations that occur faster than the charac-
teristic time are considered “rapid” and have a high risk to create very dangerous
transients, flow control operations that are slower than the characteristic time
are considered “gradual,” and operations with a time much greater than 2L/a
are considered “slow.” The characteristic time can be used to determine which
hydraulic transient flow model to use. For “slow” operations, it is possible to use
the Rigid-Column Model, and for all other conditions, the Elastic Model is used.
ΔP = aρ f Δv
where
ΔP = pressure change due to hydraulic transient (Pa, psf ) and
Δv = change in fluid velocity (m/s, fps).
For the case where the fluid velocity goes to zero rapidly due to the fluid stop-
ping suddenly, the maximum pressure near the control that stopped the flow
can be estimated:
where
ΔPmax = maximum pressure change (Pa, psf ) and
vi = initial fluid velocity (m/s, fps).
Solution
Determine the pressure wave celerity using the fluid and pipe material prop-
erties. Since the joints are allowed to expand, the pipeline support factor (Ψ)
is equal to 1.
For the fluid velocity going to zero, vi = 6.5 fps the maximum pressure due to
the hydraulic transient can be determined.
The total pressure in the pipeline would be the operating pressure plus
the pressure increase caused by the hydraulic transient.
Design Strategies
When designing pipe network systems, the engineer will need to consider the
technical requirements and demands on the system. Sufficient safety factors,
such as the pipe’s wall thickness, sufficient storage, and emergency operations
need to be evaluated to minimize transient conditions. A long pipe with a
high flow velocity and minimal branching is most vulnerable to hydraulic tran-
sients. Emergency flow control scenarios should be analyzed and tested during
the design phase to ensure that the proper pipe diameter and wall thickness,
and system components are specified.
The topographic conditions of the pipe network need to be considered. It is
common to have local high points along the pipe network. It may be more
desirable to reroute the pipeline to avoid high points where low-pressure con-
ditions exist that are prone to the low-pressure pulses caused by transients. The
presence of high points could lead to pipe collapse or column separation within
the pipe if the pressure pulse is extreme.
Selecting larger diameter pipe will lower the fluid velocity within the pipe.
This could reduce the risk of damage caused by transients; however, the eco-
nomic cost would increase, and water age within the pipe would be higher
which could reduce the water quality and disinfectant residual. For long pipe-
line systems, the diameter should be selected to optimize the construction and
operational cost. In most cases, long pipelines will require transient protection
devices.
Steel, plastic, and thin-walled ductile iron pipes are most prone to damage
caused by hydraulic transients. The pipe will be weakened by repeated exposure
to transients either due to the pressure wave itself or by vibrations caused by the
transient. A pipe weakened by age or corrosion would also be vulnerable. The
engineer may specify a more expensive, thick-walled pipe in areas that are most
vulnerable to transients. Also, the type of joint specified could be influenced
by the possibility of transients.
With the use of a transient model, the engineer can test different valve and
pump operations to minimize the effect of hydraulic transients. The most severe
transients with pumps typically occur with a power outage when the pump
suddenly stops. For normal pump operations, a slow opening/closing valve will
gradually change the fluid velocity, which reduces the pressure pulse created by
transients. Different ramping speeds with soft starts, or using variable speed
drive motors on pumps, can be tested to minimize transients, although this
is typically accomplished much more economically with slow opening/clos-
ing valves near the pump. The control valve after the pump should be opened
slowly after the pump started and closed slowly before the pump shuts down.
Protection Devices
When possible, design strategies should be used to minimize transients.
However, if transients cannot be avoided, then specific protection devises to
control transients may be needed.
Surge tanks—Surge tanks work by allowing water to bleed out of the system
during high pressures and add water to the system during low pressures. A
surge tank is a small tank that has a water height equal to the hydraulic grade
line (HGL) and feeds the protected pipe by gravity (see Figure 7.3). The surge
tank is sized to provide adequate hydraulic capacity during low-pressure con-
ditions, and the connecting valve is sized to prevent the tank from filling or
draining too slowly or quickly. The transient is dissipated as water enters or
leaves the tank. These devices should be equipped to allow for water to circu-
late into the tank to keep water fresh and if necessary, these devices should be
protected from freezing in colder climates.
A one-way tank is a type of surge tank. It will store water under atmospheric
conditions, and through a check valve allow water to enter the protected pipe
during low-pressure transients. This will reduce the low-pressure transient.
Water level
Bypass
Air compressor
x
(A) Air chamber (B) Surge tank
After the transient event, a float valve opens and allows for the gradual refilling
of the one-way tank. The advantage of using the one-way tank is the water
level in the one-way tank can be much lower than the HGL. It is common to
place one-way tanks on the discharge pipe after a pump in locations that are
vulnerable to low-pressure transients after the pump shuts off or experiences a
sudden failure (e.g., power failure).
Hydropneumatic tanks—A hydropneumatic tank is also known as an air
chamber or gas vessel. A hydropneumatic tank is a closed tank that con-
tains compressed air, or nitrogen, above a volume of water (see Figure 7.3).
The gas or air pressure in the chamber can be preset to gauge zero or higher.
During conditions of low pressure in the protected pipe, the compressed air
will expand and force water from the tank into the pipe, which could avoid
or reduce column separation and surges. During conditions of high pressure
in the protected pipe, the air or gas in the tank acts as a “cushion” against the
transient. These devices should be protected from freezing in colder climates.
There are three types of hydropneumatic tanks in HAMMER. A sealed tank
is simply a closed pressure vessel allowing water to move between the tank
and protected pipe based on the pressure within the pipe. A hydropneumatic
tank can also be vented with an air valve to allow air to enter the tank during
conditions of low pressure and effectively become a one-way surge tank. They
can also be vented with an air valve and dipping (ventilation) tube inside. If
the HGL drops and lowers the tank water level below the dipping tube, the air
valve opens and the tank effectively becomes a surge tank open to the atmos-
phere. When the HGL later increases, the air is expelled from the tank through
the air valve and the water level will go above the bottom of the dipping tube,
and the tank will act as a sealed pressure tank again.
Surge relief valves—In the case of a high-pressure transient, a surge relief (pres-
sure relief ) valve can be placed to allow the transient pressure to be relieved.
These valves are set to open and allow water to discharge when the pressure in
the protected pipe exceeds a predetermined value. This controls the maximum
pressure within the protected pipe. After the high-pressure transient is relieved,
the valve will close slowly to avoid causing another transient. Similar to a surge
relief valve is a surge-anticipator relief valve. This valve is controlled by the sys-
tem emergency operations and will start to open immediately in anticipation
of a high transient, usually because of a power failure.
Air inlet valves—Air inlet valves are also called vacuum breakers. These valves
are installed in high points of the pipe network to control vacuum conditions.
During low-pressure transients, these valves will allow air to enter the system.
This will control the minimum pressure within the protected pipe to avoid col-
umn separation or pipe collapse. An adequate period of time should be allowed
for the air that entered the system to be expelled before normal operations are
allowed to restart.
Rupture disk—A rupture disk are also called bursting disk or burst diaphragm.
A rupture disk is a one-time use diaphragm or membrane that will rupture
when a predetermined high or low pressure ruptures the disk. This will create a
“designed” opening in the protected pipe that can be later be easily repaired by
replacing the ruptured disk with a new one. Many times a rupture disk is used
as a backup protection device to a surge relief or air inlet valve.
• reduce the risk of external water entering a piping network that can con-
taminate distribution water; and
• identify locations or pipes that are of high risk to damage caused by a
hydraulic transient.
HAMMER can be used as a stand-alone program, or interface with
MicroStation, AutoCAD or ArcGIS, HAMMER shares the same file format
of WaterCAD and WaterGEMS so models can be easily shared between these
applications. It will also import files from EPANET.
7.6 Problems
Solve the following problems using the HAMMER computer program.
1. A reservoir discharges to atmosphere through 5,000 ft of 12-in. pipeline,
as shown in the following diagram.
Reservoir
elevation = 350 ft
Valve
elevation = 100 ft
Typical pressure drop = 100 psi
Typical flow = 3,000 gpm
Discharge element type = valve
Initial status = open
5,000 ft of 12-in. diameter ductile iron pipe Expansion joints throughout Discharge to
atmosphere
Hazen–Williams C = 120 Wall thickness = 0.28 in.
2. A pump draws water from reservoir R-1 and delivers it to reservoir R-2
through 10,500 ft of 12-in. pipeline. The pipeline passes over a high point
at Junction J-1, as shown in the following diagram.
Reservoir R-2
Junction J-1 Elevation = 1,000 ft
Elevation = 900 ft
Pipe data
ft
3,
00
12-in. ductile iron 00
2,0
Hazen–Williams C = 120 0
ft
Wall thickness = 0.28 in.
Expansion joints throughout
ft
00
5,0
Junction J-2
Elevation = 400 ft
Reservoir R-1
Elevation = 500 ft
500 ft
Pump
Elevation = 490 ft
Assume the fluid is water at 68°F. Set the Transient Friction Method to
Unsteady, Run Duration Type to Time, Run Duration to 240 s, and use a
Vapor Pressure of −14.2 psi.
(a) What are the steady-state pressure and hydraulic grade line (HGL) at
Junction J-1 with the pump running?
(b) What are the maximum pressure and maximum HGL at Junction J-1
following a power failure?
(c) If Junction J-1 is replaced with a combination air/vacuum valve (4-in.
air inlet diameter, 1/4-in. air outlet diameter), what are the maximum
pressure and maximum HGL at the air/vacuum valve following a power
failure?
8.2 Loading
Sanitary sewer loads are divided into two categories: sanitary loads and
wet-weather loads.
• Sanitary loads (dry-weather loads) result from human activity and are not
weather dependent. Common sources of sanitary loads are various residen-
tial, commercial, recreational, and industrial uses.
• Wet-weather loads are related to rainfall activity such as groundwater
infiltration (water leaking into a pipe through cracks, joints, and defects)
and structure inflow (surface water entering a structure through openings
around the cover, or due to a missing cover).
Time-based loads can be classified as either wet-weather or sanitary loads, and
they may have hydrograph- or pattern-based formats.
5.5
5.0
4.5
4.0
3.5
Flow, cfs
3.0
2.5
2.0
1.5
1.0
0.5
0.0
0 5 10 15 20 25 30 35
Time, h
Pattern Loads
A pattern load comprises a base load and an associated loading pattern. The
pattern is a series of multipliers relating sanitary sewer loading to time, and
it allows the user to automatically apply time-variable changes within the
system.
Different categories of users, such as residential or industrial customers, may
be assigned different patterns to accurately reflect their particular load varia-
tions. A diurnal curve is a type of pattern that describes change in loading
over the course of a daily cycle. A dry-weather diurnal curve for sanitary load-
ing often resembles the diurnal curve for water use at the same location (see
Section 6.1), reflecting times when people are using more and less water than
average.
Patterns are typically based on a multiplication factor versus time relationship,
in which a multiplication factor of 1.0 represents the base value (usually the
average value). This relationship can be expressed:
Qt = At × QBase
where
Q t = load at time t,
At = multiplier for time t, and
Q Base = baseline load.
For the residential usage, diurnal curve in Figure 8.2, we see that there is a peak
in the diurnal curve in the morning, another slight peak around noon, and a
third peak in the evening. During the nighttime hours, the pattern reflects the
relative inactivity of the system. (Note that this curve is conceptual and should
not be construed as representative of any particular network.)
As explained in Chapter 6, there are two basic forms for representing a pattern:
stepwise and continuous. A stepwise pattern is one that assumes a constant
level of usage over a period of time and then jumps instantaneously to another
level where it again remains steady until the next jump. A continuous pat-
tern is one for which several points in the pattern are known and sections in
between are transitional, resulting in a smoother pattern. For the continuous
pattern in Figure 8.2, the magnitude and slope of the pattern at the start and
end times are the same—a continuity that is recommended for patterns that
repeat.
2.0
1.8
1.6
1.4
1.2
Multiplier
1.0
0.8
0.6
0.4
0.2
0.0
0 5 10 15 20 25
Time, h
Base Load
Unit Sanitary Load Loading Unit
(l/day per Loading Unit)
Airport Passenger 10
Bar Customer 8
Hotel Employee 40
Restaurant Meal 10
The most common type of EFF is the variable peaking factor (PF). In this
method, a computed EFF is multiplied by the specified base flow.
where
Qpeak = transformed flow (l/day, gpd),
Qbase = base flow (l/day, gpd), and
EFF = extreme flow factor (unitless).
• Babbitt (1958)
5.0
PF =
(P 1, 000)
0.2
• Harmon (1918)
14.0
PF = 1.0 +
4.0 + ( P 1, 000)
0.5
P
18 +
PF = 1000
P
4+
1000
2.69
PF =
Q 0.121
4.0
3.5
Solution
First, to get the average base load, find the amount of load per unit for the
three loading categories mentioned in Table 8.1. Then, multiply this unit flow
by the number of units.
14
PF = 1.0 + = 4.23
4 + (110 1, 000)
0.5
Apply this PF to the base loads calculated from the subdivision. Then, apply
the given PF of 2 to the load derived from the shopping center to compute the
approximate peak flow.
4.23 × (22,400 l/day + 11,400 l/day) + 2 × 400 l/day = 144,000 l/day
Wet-Weather Loading
Wet-weather loading represents the intrusion of ground and surface water into
the sewer system and is divided into two components: infiltration and inflow.
Infiltration represents groundwater entering pipes through cracks, holes,
and other defects. Inflow represents wet-weather flow that enters the sewer
through manhole covers, roof drains, catch basins, sump pumps, and so forth.
Collectively, these flows are referred to as I/I. It is important to minimize intru-
sion of rain water into the sewer system because of the cost of treating the added
flow. Also, wet-weather loads increase the potential for overloading the system,
which causes flow to back up into houses or overflow into the environment.
Infiltration causes noticeable increases in sewer flow, especially after rainstorms
and subsequent rises in the groundwater table. Sewer pipes placed in low areas,
such as near creek beds, are more subject to infiltration than those placed in
higher ground because they are more likely to be located below the water table.
Older sanitary sewers are more susceptible to infiltration than newer sewer
systems because modern techniques and materials used in newer systems can
generally prevent or greatly reduce infiltration.
Infiltration at an individual pipe is generally hard to quantify. Instead, infiltra-
tion is estimated for the whole system or a portion of the system. This global
infiltration is then distributed to individual pipes using one of the following
methods:
• proportional to pipe length
• proportional to pipe diameter/length
• proportional to pipe surface area
• proportional to some user-defined count (e.g., number of defects in
the pipe)
Inflows are unwanted flows entering from the surface and are modeled based
on estimates or requirements from local or state regulatory agencies. Like infil-
tration, inflows can be quite significant and can vary considerably depending
on the depth of runoff, the size of the manhole cover, and the number of
openings in the manhole cover.
Routing Overview
As a hydrograph is routed through a conduit, it undergoes changes in shape
and temporal distribution caused by translation and storage effects.
SewerCAD uses one of two methods to determine the shape and distribution
of a hydrograph routed through a gravity pipe.
• convex routing
• weighted translation routing
Convex Routing
The underlying assumptions of the convex routing method are that the routed
outflow for a time step is based on the inflow and outflow for the previous
time step, and that the flow does not back up in the pipe (i.e., no reverse
flow or reduced flow due to tailwater effects exists). Each outflow ordinate is
calculated as:
where
Ot + Δt = outflow at time t + Δt;
t = current time (s, min),
Δt = hydrologic time step (s, min),
c = convex routing coefficient,
It = inflow at time t(l/s, gpm), and
Ot = outflow at time t(l/s, gpm).
V Δt
c = Δt =
L tt
where
Δt = hydrologic time step (s),
tt = travel time (s),
V = velocity established for representative flow (m/s, ft/s), and
L = length of pipe (m, ft).
In SewerCAD, the velocity used to calculate the coefficient is either the normal
velocity or full velocity generated for a user-specified percentage of the peak
of the inflow hydrograph. In other words, if the percentage of the peak flow
is greater than the capacity of the pipe, the full flow velocity is used. If the
percentage of the peak flow is less than the pipe capacity, the flow velocity for
normal depth is used.
The higher the percentage of flow, the faster the velocity used to calculate the
convex routing coefficient, and the closer the routed hydrograph will be to a
pure translation of the inflow hydrograph.
The user-specified percentage can be modified in the calculation options. A
typical value is around 75%, but this value may be modified for oddly shaped
hydrographs with sharp, uncharacteristic peaks or for calibration purposes.
1 ⎛ 1⎞
Ot = It −Δt + ⎜⎜⎜1 − ⎟⎟⎟ It
c ⎝ c⎠
where
Ot = outflow at current time step (l/s, gpm),
c = convex routing coefficient
It − Δt = inflow at previous time step (l/s, gpm), and
It = inflow at current time step (l/s, gpm).
• perform EPSs that include time-variable loads and hydrologic routing; and
• animate plans and profiles showing sanitary sewer system performance
over time.
The theory and background used in SewerCAD are presented in this chapter
and in the SewerCAD online help system.
SewerCAD can analyze complex sewage collection systems under a variety of
loading conditions. SewerCAD’s powerful loading capabilities can be used to
track system response to an unlimited range of sanitary and wet-weather load-
ing combinations. It can be used for
• pipe sizing
• pump sizing
• master planning
• operational studies
• rehabilitation studies
8.5 Problems
1. In the sanitary sewer represented by the following diagram, a load from
a development of high-rise apartments (at wet well WW-1) with 10,000
residents is pumped to the top of a hill, where it is then transported
via a circular gravity line to a treatment plant (represented by O-1). At
manhole MH-1 is a resort apartment with 300 guests. The load gen-
erated here flows down a circular gravity pipe to the same gravity line
mentioned above. A bar that serves, on average, 50 people per day and 2
large cafeterias with 20 employees each are located near manhole MH-2.
MH-1
P-1
-2 MH-2 P-2
FM
PMP-1
-1
FM
O-1
WW-1
MH-1 MH-2
2
P-
P-
1
MH-3
JC-1
P-
3 4
P-
O-1
(a) List the diameter and slope for each newly designed pipe. Are there
any problems with the designed system?
(b) What percentage of the total flow is wet-weather flow? Which pipe
has the most infiltration?
(c) For a larger magnitude storm, the inflow rate into each manhole
is estimated at 100 l/day. Analyze the model using the previously
designed system and apply the larger wet-weather loading. What
percentage of the total flow is wet-weather loading?
(d) Do you consider the amount of wet-weather flow into the system
significant? What are some methods for alleviating infiltration?
3. The following network is an initial design for a system of force mains in
a sanitary sewer. All pipes are ductile iron (Hazen–Williams’ C = 130).
The pump, PMP-1, is at an elevation of 691 m. Enter the system data
given in the tables below and answer the following questions. The ground
elevation at the outlet is 715 m and the sump elevation is 712 m.
Hints: Make sure that the Fixed Level in Steady State box is checked under
the “Section” tab of the Wet Well dialog. In addition, before running the
model, make the following modification to the calculation options: within the
GO dialog, click the Options button, scroll to the right, and in the Pressure
Hydraulics tab, check the Use Pump Loads box.
FM-5
FM-7
FM-8
J-4
J-5 FM-6
O-1
(d) If the minimum velocity required in the force main to keep particles
from settling is 0.6 m/s, which areas are going to have problems?
(e) What are some possible changes to the design to fix the problem
portions of the system?
(f) Implement a solution suggested in part (e). Describe the fix(es) and
fill in the following chart.
MH-6
MH-10
P-6
P-10
MH-11 MH-7 MH-1
P-11
P-7
P-1
P-8
MH-8
MH-9 MH-2
P-9
P-2
O-1 P-5 P-4 MH-4 P-3 MH-3
MH-5
-2
FM
1
P-12 FM- PMP-1
MH-12 WW-1
All manholes have a sanitary base load of 100 gpm and a continuous diurnal
pattern applied to them as defined below. Each manhole also has wet-weather
loading given by the following hydrograph.
Variable Steady State EPS Sanitary Only EPS Dry and Wet
Maximum flow at outlet (gpm)
Time of max. flow at outlet (h) N/A
Maximum velocity in system at
this time (ft/s)
9.1 Introduction
It was a typical summer, but one accentuated with periods of extreme heat
and intense rainfall. One area in the city was nearing the completion of a
2-year beautification and utility project that had local residents and mer-
chants convinced would solve a history of minor, inconvenience flooding.
Unfortunately, one of the thunderstorms resulted in damaging flooding at
multiple locations. The headlines in the local paper the next morning read
“2.7 Inches of Rain in 70 Minutes: $50 Million to End Flooding.” The State
Climatologist reported the event return period was between 5 and 10 years.
Certainly, a modern urban drainage system should accommodate runoff
from such an event, particularly since local regulations require sewer systems
designed to accommodate peak flows from a 25-year event. What happened?
Why was there so much flooding?
When they designed the local drainage system, engineers used traditional meth-
ods to size the storm sewer pipes, culverts, and open channels to accommodate
estimated peak flows. They did not consider the impact of limited storage in
conduits, and the additive effects and timing of local inflows. Nonuniform
flow profiles were determined to delineate the floodplain. Otherwise, down-
stream boundary conditions and interactions among the different drainage
system components were not considered. All calculations were performed
assuming steady flows.
The crux of the problem was that the steady flow methods did not account for
the dynamic response of the watershed and its internal drainage system. No true
picture, that is, understanding of the drainage system performance was estab-
lished. The results were insufficient and potentially misleading. Consequently,
there was a lack of adequate planning and flood prevention measures.
To overcome these deficiencies, the city undertook a study to identify the
causes of flooding and possible solutions. That study included a hydraulic
analysis of the drainage system using a dynamic model. The results provided
a clear understanding of the performance of individual components and their
interactions at a system wide scale.
The purposes for this chapter are to introduce dynamic modeling and demon-
strate its application to the analysis and design of hydraulic systems using the
SewerGEMS computer software. The option to include low-impact develop-
ment (LID) to reduce the overall surface water runoff due to land planning and
green infrastructure will be discussed. Example applications include stormwa-
ter conveyance systems that include green infrastructure design to size storm-
water sewer networks.
accurate because it improves the hydraulic wave effects observed in the field
during periods when there are significant changes in loading patterns, such as
weekends and holidays when water consumption and wastewater production
patterns change, sometimes dramatically.
This discussion has identified some of the reasons and applications for unsteady
modeling. Although steady flow equations are significantly simpler and easier
to use, their solutions are accurate only for those situations in which the chan-
nel and flow conditions do not vary greatly. When faced with more complex
modeling problems, the unsteady flow equations provide far more appropriate
solutions.
The design objectives for LID in urban areas is to manage the peak flow rates
and maximum water level of stored water, while simultaneously protecting
the urban infrastructure (buildings, roads, and other structures), human lives,
and surrounding environment. The use of LIDs, since they tend to be smaller
and decentralized, requires the hydraulic modeler to carefully examine the pre-
dicted runoff peak flows rates, volumes, and resulting hydrographs throughout
the system design since potentially many LID structures are in use at one time.
As the smaller LID structures fill and overflow into other traditional structures
used to manage the runoff, a dynamic model with detail of how all these struc-
tures are hydraulically connected is needed.
A challenge to add LIDs to existing urban areas is to find locations to place pro-
posed LIDs. The acquisition of land to build these structures is limited in many
urbanized areas and if a site is found, it may not be optimally placed. Also, the
soil for many urban areas tends to have a low infiltration rate or additional
water in the urban subsurface could have undesirable consequences, which
would make LIDs that use infiltration impractical to implement. Regardless of
the challenges, the use of LIDs for managing stormwater is gaining popularity,
and it is common to evaluate many design alternatives when these structures
are considered. A detailed hydraulic model of the stormwater system can be
used to evaluate how LIDs should be placed and sized to manage the runoff
and improve the water quality leaving urban areas.
∂Q ∂( A + Ao )
+ −q =0
∂x ∂t
∂Q ∂(βQ 2 / A ) ⎛ ∂y ⎞⎟
+ + gA ⎜⎜⎜ + So + S f + Si ⎟⎟ + L = 0
∂t ∂x ⎝ ∂x ⎠
where
Q = discharge (cfs, m3/s)
Numerical Solution
There is no known analytical solution to the St. Venant equations; therefore
they are solved numerically. The solution method in SewerGEMS uses a finite
difference scheme that converts the system of partial differential equations to
a system of algebraic equations involving the unknowns at discrete points in
time and space.
Specifically, a weighted four-point implicit scheme is used. An implicit scheme
is preferred over explicit since these schemes have the advantage of maintaining
numerical stability for large computational time steps and exhibit robustness in
modeling systems that integrate the complex hydraulic interactions encountered
in gravity sewer systems. The scheme was chosen because it handles unequal dis-
tance steps, its stability–convergence properties can be conveniently modified,
and the internal and external boundary conditions can be easily applied.
The numerical analog uses the following four-point finite-difference quotients
to approximate the first- and zero-order derivative terms:
∂f f j +1 + fi +j +1 1 − fi j − fi +j 1
= i
∂t Δt
∂f f j +1 + fi +j +1 1 − fi j − fi +j 1
= i
∂t Δt
θ( fi j +1 + fi +j +1 1 ) + (1 − θ )( fi j + fi j+1 )
f =
2
dx
C+ : =V −c
dt
and
dx
C− : =V −c
dt
where x,t are longitudinal distance and time, v is cross-sectional average veloc-
ity, and c is the speed at which the wave effect created by the disturbance trav-
els. This wave speed is called celerity and is evaluated as
c= gD
where g is gravity and D is the hydraulic depth of water (see Chapter 1).
For flow in a rectangular channel, the Froude number is expressed in terms of
wave celerity as
V V
F= =
gD c
the occurrence of the inflow and outflow hydrograph peaks or the centers of
mass. Attenuation is the reduction in peak flow rate resulting from the spread
of the hydrograph shape. Spreading occurs due to diffusion and dispersion
caused by varying conduit conditions that alternately accelerate and decelerate
the flow, for the cases of channel and sewer routing, and hydrograph storage
and release, for the cases of pond and pump station routing. Spreading results
in elongation of the time base, reduction of the peak flow rate, and skewing
of the hydrograph shape. This leads to the development of a “tail” to the
recession limb. Diffusion is spreading that does not involve mass (or volume)
shift, whereas dispersion involves mass transfer from the rising limb to the
recession limb.
Diffusion Wave
Two approximations to the dynamic wave model are the diffusion and kine-
matic wave models. The diffusion wave model that assumes the inertia terms
in the momentum equation are negligible compared with the force terms for
pressure, friction, and gravity. The diffusion wave equations are continuity and
the following reduced form of momentum.
∂y
= So − S f
∂x
For prismatic channels, the diffusion wave equations can be combined into the
single equation
∂Q ∂Q ∂Q
2
+c =D
∂t ∂x ∂x
2
Kinematic Wave
The kinematic wave model neglects all dynamic terms, which includes the
derivative terms for local and convective acceleration, and pressure imbalance.
This reduces the momentum equation to
S f ≈ So
which means it can be approximated with one of the uniform flow formulae.
The kinematic wave model accounts only for hydrograph translation and not
for any diffusion or dispersion. The most common application for kinematic
wave modeling is overland flow during periods of stormwater runoff.
Muskingum Routing
The Muskingum routing method was developed for use in a flood control
study in the Muskingum River basin in Ohio, hence the origin of its name.
This method uses a lumped parameter model that treats the channel reach of
interest as a single spatial element, that is, the reach is not divided into shorter
segments, and the channel cross-sectional geometry, slope, and roughness are
assumed constant. Model input is the inflow hydrograph at the upstream end
and model output is the outflow hydrograph from the downstream end. This
method does not compute flows at intermediate locations along the reach as do
hydraulic models.
The model governing equations are the hydrologic form of continuity
dS
I−Q =
dt
S = K [θ I + (1 − θ )Q ]
The subscripts refer to inflow and outflow values at the beginning and end of a
computational time step, Δt.
The Muskingum method traditionally has been applied to channels where
there are recorded inflow and outflow hydrographs that are used to determine
the model parameters that give a best fit between observed and simulated
outflow hydrographs. As such, the model is calibrated to the channel. It then
is applied to other inflow hydrographs to predict the outflow hydrographs.
While the Muskingum method has been used for decades, values for the coeffi-
cients, K and θ, are not readily available or easily obtained for measured hydro-
graphs. A general rule of thumb is that K can be estimated by the travel time
through the reach and θ lies between 0.2 and 0.4. The upper limit for θ is
0.5. Values less than 0.5 define dispersing hydrograph shapes. θ equal to 0.5
defines a hydrograph that does not change shape and only exhibits translation
effects. A value of 0.0 is the practical lower limit for θ and suggests inflow has
little effect. Some users have reported values for θ from 0.4 to 0.5 for natural
Muskingum–Cunge Method
To overcome the limitation with parameter estimation and extend the method
to ungaged channels, Cunge (1969) developed hydraulic-based equations to
estimate the routing parameters by matching coefficients of the advective–
diffusive transport form of the Muskingum routing model
∂Q Δx ∂Q ⎡ 1 Δx ⎤ ∂ Q
2 2
+ = ⎢Δx (1 − θ )ck − ⎥
∂t K ∂x ⎢⎣ 2 K ⎥⎦ ∂x 2
∂Q ∂Q QR ∂ Q
2
+ cK =
∂t ∂x 2TSo ∂x 2
Δx
K =
cK
1⎡ QR ⎤
θ= ⎢1 − ⎥
2 ⎢⎣ TSo cK Δx ⎥⎦
When these equations are used to estimate the routing parameters, the method
is known as the Muskingum–Cunge routing method. Numerical experience
that has shown negative values for θ can sometimes occur but do not cause
numerical problems.
An important feature and a perceived advantage of the Muskingum–Cunge model
is that it is diffusive for θ<0.5 and offers the same advantages of the diffusion wave
model. This method is included in many contemporary stormwater models.
9.8 Problems
The following diagram shows the stormwater conveyance system for an area
that underwent staged development from agricultural to suburban land use.
CM-1 was developed first, prior to any stormwater regulations. Following
standard practice at the time, stormwater runoff was collected in the sewer
system and diverted to an outfall into the creek at cross section of CS-2.
CM-2 was developed after stormwater regulations were in place. To control
the runoff peak from a 10-year design rainfall at or below the predevelopment
rate, a detention pond was constructed between the outlet from CM-2 and
the point where local runoff enters the sewer system at MH-3. During the
development period, CM-3 was built out for mixed land use—multifamily
residential, commercial, and a park. Runoff from CM-3 flows directly to the
existing creek.
Tables are included that provide data for the catchments, channel, cross sections,
conduits, manholes, and pond. The pond outflow structure is a two-stage riser
with lower stage orifice (two 6-in. diameter openings at elevation 533.25 ft MSL,
orifice coefficient = 0.60) and an upper stage rectangular, sharp crested weir
(length = 1 ft, weir coefficient = 3.10, and crest elevation at 536.75 ft MSL).
Downstream of the confluence of the storm sewer and creek, there is a road
crossing with a culvert cross-drain. Conduit CO-9 is the culvert. Note that it
connects CS-3 and CS-4. This system discharges to a free outfall.
Following build-out of CM-1, CM-2, and CM-3, the road crossing began
to flood (overtop) more frequently. Local residents noted the sewer sys-
tem between manholes MH-4 and MH-6 was surcharging with occasional
aboveground flooding during the events when the crossing overtopped. This
resulted in minor system flooding. It also was noted during these events the
water depth in the creek was above the crown of the outfall sewer. The residents
quickly concluded that the minor system flooding they were experiencing was
due to an undersized culvert at the road crossing and demanded that the city
make appropriate repairs. The city engineer analyzed the system and deter-
mined the inconvenience flooding they were experiencing was not due to an
undersized culvert, but probably to an excessive pond release rate.
Apply SewerGEMS to analyze this system. Use the default curvilinear NRCS
unit hydrograph, NRCS Type II rainfall distribution, and a 24-h rainfall depth
of 5.15 in.
Were the citizens correct about the cause of flooding in the sewer system or was
the city engineer correct?
Were both parties wrong? Was the flooding due to other conditions?
Would it be a more appropriate solution—not necessarily cheaper in terms of
material costs—to replace one or more links of sewer with a larger pipe? If so,
which one(s)?
Hint: Run animations for a profile between the detention pond and the system
outfall.
The accompanying tables provide data about the system components and tidal
data from a nearby U.S. Geological Survey gaging station. The system has a
free outfall into a tidal creek, that is, it has a tidal boundary condition. The
outfall ground elevation is 5.00 ft MSL with invert elevation of 3.80 ft MSL.
Apply SewerGEMS to analyze this system. Use the default curvilinear NRCS
unit hydrograph, NRCS Type III rainfall distribution, and a 24-h rainfall
depth of 5.30 in.
Part of the design is setting curb height to contain any surcharge overflow. Is a
curb height of 4 in. sufficient?
Does the flow in the sewer reverse during high tide? In other words does water
flow from the tidal creek into the sewer? Or does any system surcharge result
due to a “backwater effect?”
What is the duration of any about ground surcharge during the design event?
What signage should be placed along the lower portion of the parking lot?
What would be the effect of upsizing one or more of the sewer pipes? Explain
your answer in terms of the effects of in-line storage. Does this suggest in-line
storage can be a solution at other, noncoastal locations?
A nearby town proposes to add a new school onto its wastewater collection
system. The new school will discharge into MH-4 and connect with the system
trunk line at MH-3. From there, the wastewater flows to the local treatment
plant located at OF-1. There are concerns about system capacity, particularly
since this addition is immediately downstream of the location where a man-
ifold pump station system discharges into the trunk line. Specific concerns
include the potential for system surcharge with surface overflow at MH-3.
The figure and accompanying tables summarize pertinent system data. The
system has a free outfall at the head of the wastewater treatment plant with
ground elevation of 263.70 ft MSL and invert elevation of 260.00 ft MSL.
Use SewerGEMS to analyze this system. Perform simulations for a full diurnal
loading pattern at MH-1, WW-1, WW-2, and MH-4. Both pump stations are
duplex stations, so use two pumps at each station, operating in lead–lag status.
It is highly recommended that each station has two identical pumps. Data are
provided for two recommended pumps; evaluate system performance for both
pumps at each pump station.
Will MH-3 surcharge? If so, will it surcharge aboveground? If so, is the flood-
ing at MH-3 sensitive to which pump is used at the two pump stations?
What are the maximum and minimum pressures at PJ-2?
Local regulations stipulate the minimum pump runtime must be 2 min or
longer, the average number of pump cycles per hour should be in the range of
4–6, and wet well detention time should not exceed 30 min. Does each pump
station meet regulations?
What are the minimum and maximum HGL values at MH-3?
The following diagram shows the layout for a stormwater conveyance system
proposed for a 400-ft segment of a local street. Runoff from the street and adja-
cent areas (indicated as catchments CM-1, CM-2, and CM-3) flows into the
gutter and enters the underground sewer system through inlets at three catch
basins. The gutter is a triangular channel formed by the intersection of a 6-in.
barrier curb and the transversely sloped (H:V = 40:1) street surface. Pertinent
data about system components are given in the tables below.
The system has a free outfall with ground and invert elevations of 97.50 and
93.00 ft MSL, respectively.
Apply SewerGEMS to analyze the system. Use the default curvilinear NRCS
unit hydrograph, NRCS Type III rainfall distribution, and a 24-h rainfall
depth of 5.10 in. Use the HEC-22 Headloss Method to compute head losses in
the catch basins and manholes.
Will the system surcharge, that is, have HGL above the crown of the pipe, at
one or more locations if the tailwater at the free outfall is at elevation 95 ft?
If so, which pipes, manholes, and catch basins surcharge? Will the roadway
flood? If so, where?
If the system will surcharge, can the problem be resolved by upsizing one
link of pipe? If so, which link should be upsized? What is the minimum
upsized pipe diameter? What are the implications about in-line detention
storage?
Are any of the proposed catch basins undersized?
Regulations governing stormwater ponds vary from the classic statement “the
post development peak discharge cannot exceed the pre-development peak
discharge for a given design storm event” to the 10% and no net increase
rules. Under the 10% rule, the engineer must demonstrate no increase in peak
flow over predevelopment conditions at a downstream location where the land
use change area is no more than 10% of the total contributing area at that
point. At watersheds where there is substantial hydrograph attenuation due
to in-stream storage and routing effects or the presence of a pond would delay
the timing of peak outflow from a downstream pond such that it more closely
aligned with the peak runoff from upstream areas, thereby creating a greater
peak than would have occurred naturally, this criterion may actually be used
to show no pond is required on-site. Calculations to demonstrate this clearly
must include unsteady pond and channel routing.
The no net increase rule extends the requirements for stormwater ponds beyond
peak runoff control and also requires no net increase in runoff volume. Some
interpretations also hold the outflow hydrograph shape and timing must not be
altered substantially. Here again, unsteady pond and channel routing is required.
Given a 100-ac watershed with the following land use conditions: 20 ac in for-
est (good condition) over HSG C soils; 40 ac in open space (good condition)
over HSG B soils; 35 ac in single-family residential over HSG B soils; and 5 ac
in commercial land use over HSG B soils. The watershed hydraulic length is
3,750 ft and the area weighted average watershed slope is 1.6%.
At the upper end of this watershed, it is proposed to develop 10 ac of the open
space land use to single-family residential land use. Runoff from this develop-
ment will be controlled by a stormwater pond (inverted quadrilateral frustum
with 3:1 side slopes and length to width ratio of 2:1), which is to be designed
according to controlling regulations. The 10-ac site currently has an average
slope of 1.8%, which after development will be reduced to 1.5%. The hydraulic
length is 1,280 ft.
The channel connecting outflow from the 10-ac development to the outlet
of the 100-ac watershed (10 ac or 10% plus 90 ac or 90%) is 3,200 ft long
with an approximate trapezoidal cross section with the following parameters:
base width = 8 ft; side slope Hor:Vert = 2.5:1.0; longitudinal slope = 0.3%;
Manning n-value = 0.065.
The design rainfall is a 24-h event with 4.80-in. accumulation. Use the NRCS
synthetic unit hydrograph.
Apply SewerGEMS to design a stormwater pond at the outlet of the 10-ac
development site such that (a) the peak flow at the outlet of the 100-ac basin
(downstream point to check the 10% rule for the 10-ac site) does not exceed the
predevelopment rate for the design storm event and (b) the peak and volume at
the outlet of the 100 ac do not violate the “no net increase rule.”
Clark, R.M., W.M. Grayman, R.M. Males, and A.F. Hess. 1993. Modeling
Contaminant Propagation in Drinking Water Distribution Systems.
Journal of Environmental Engineering, 119(2), 349.
Cunge, J.A. 1969. On the Subject of a Flood Propagation Computation
Method (Muskingum Method). Journal of Hydraulic Research, 7(2), 205.
Debo, T.N. and A.J. Reese, 1995, Municipal Storm Water Management, CRC
Press, Boca Raton, FL.
Featherstone, R.E. and C. Nalluri, 1982, Civil Engineering Hydraulics,
Granada, New York, NY.
Federal Aviation Administration (FAA), 1970, Report A/C 150-5320-
58, Advisory Circular on Airport Drainage, U.S. Department of
Transportation, Washington, DC.
Federation, 1982, Gravity Sanitary Sewer Design and Construction, American
Society of Civil Engineers, New York, NY, ASCE Manual of Practice,
No. 60.
Frederick, R.H., V.A. Meyers, and E.P. Auciello, 1977, Five to 60-minute
Precipitation Frequency for the Eastern and Central United States, In
NOAA Technical Memorandum NEW Hydro-35. www.nws.noaa.gov/oh/
hdsc/studies/prcpfreq.html (accessed June 2004).
French, R.H., 1985, Open-Channel Hydraulics, McGraw-Hill, New York, NY.
Great Lakes-Upper Mississippi River Board of Sanitary Engineers, 1997,
Recommended Standards for Wastewater Facilities, Health Research, Inc.,
Albany, NY.
Hathaway, G.A. 1945. Design of Drainage Facilities. Transactions, American
Society of Civil Engineers, 110, 697.
Harmon, W.G. 1918. Forecasting Sewage at Toledo Under Dry Weather
Conditions. Engineering News-Record, 80.
Haestad Methods, 2000, CulvertMaster: Culvert Design & Analysis Software,
Bentley Systems.
Haestad Methods, 2000, FlowMaster: Open Channel, Pipe & Inlet Software,
Bentley Systems.
Haestad Methods, 2000, StormCAD: Storm Sewer Design & Analysis Software,
Bentley Systems.
Haestad Methods, 2001, SewerCAD: Sanitary Sewer Modeling Software,
Bentley Systems.
Haestad Methods, 2002, PondPack: Urban Hydrology & Detention Pond
Modeling Software, Bentley Systems.
Rossman, L.A. and P.F. Boulos. 1996. Numerical Methods for Modeling Water
Quality in Distribution Systems: A Comparison. Journal of Water Resources
Planning and Management 122(2), 137.
Rossman, L.A., 1993, EPANet User’s Manual (AWWA Workshop Edition), Risk
Reduction Engineering Laboratory, Office of Research and Development,
US Environmental Protection Agency, Cincinnati, OH.
Rossman, L.A., R.M. Clark, and W.M. Grayman. 1994. Modeling
Chlorine Residuals in Drinking-Water Distribution Systems. Journal of
Environmental Engineering 120(4), 803.
Rouse, H., and S. Ince. 1957. History of Hydraulics. Iowa City, Iowa. Iowa
Institute of Hydraulic Research.
Sanks, R.L., 1998, Pumping Station Design, 2nd ed., Butterworth-Heinemann,
Inc., Stoneham, MA.
Schaake, J.C. Jr., J.C. Geyer, and J.W. Knapp. 1967. Experimental Examination
of the Rational Method. Journal of the Hydraulics Division 93(HY6).
Schall, J.D., P.L. Thompson, S.M. Zerges, R.T. Kilgore, and J.L. Morris, 2012,
Hydraulic Design Of Highway Culverts, Hydraulic Design Series No. 5, 2nd
ed., Federal Highway Administration, Washington, DC.
Simon, A.L., 1976, Practical Hydraulics, John Wiley & Sons, New York, NY.
Singh, V.P., 1988, Hydrologic Systems, Volume I: Rainfall-Runoff Modeling,
Prentice Hall, Englewood Cliffs, NJ.
Singh, V.P., 1992, Elementary Hydrology, Prentice Hall, Englewood Cliffs, NJ.
Soil Conservation Service, 1969, National Engineering Handbook, Section 4,
Chapter 16, National Resources Conservation Service, Washington, DC.
www.wcc.nrcs.usda.gov/water/quality/common/neh630/4content.html
(accessed June 2004).
Soil Conservation Service, 1986, TR-55, Urban Hydrology for Small Watersheds,
Natural Resources Conservation Service, Washington, DC. http://www.
wcc.nrcs.usda.gov/water/quality/common/tr55/tr55.pdf (accessed June
2004).
Strafaci, A.M., 1998, Essential Hydraulics and Hydrology, Haestad Press,
Waterbury, CT.
Streeter, V.L. and E.B. Wylie, 1985, Fluid Mechanics, McGraw-Hill, New
York, NY.
Tchobanoglous, G., 1981, Wastewater Engineering: Collection and Pumping of
Wastewater, Metcalf & Eddy, Inc. McGraw-Hill, Inc.
Van Havern, B.P., 1986, Water Resources Measurements, American Water Works
Association, Denver, CO.
Wanielista, M.P., 1990, Hydrology and Water Quantity Control, John Wiley &
Sons, New York, NY.
Walski, T.M., 1984, Analysis of Water Distribution Systems, Van
Nostrand-Reinhold, New York, NY.
Walski, T.M., 1993, Water System Modeling Using CYBERNET®, Haestad
Methods, Waterbury, CT.
Williams, G.B. 1922. Flood Discharges and the Dimensions of Spillways in
India. Engineering (London) 134, 321.
Wylie, E.B. and V.L. Streeter, 1993, Fluid Transients in Systems, Prentice-Hall.
Zipparro, V.J. and H. Hasen, 1993, Davis’ Handbook of Applied Hydraulics,
McGraw-Hill, New York, NY.
With this code, it is possible to download a license file and install the software
according to the above article.
Online Help
Electronic documentation can be kept up-to-date more easily than paper doc-
umentation. It is accessible from within the program, and it is environmen-
tally friendly. For these reasons, Bentley OpenFlows’ programs have extensive
online Help systems to guide the user through the analysis and design process.
The Help system is used most often to:
• Obtain context-sensitive information about an area of the program
• Search for topics such as background theory and equations
• Use “How Do I” to get detailed step-by-step instructions for completing
a task
Context-Sensitive Help
The quickest way to access the Help topic related to your current task (i.e.,
context-sensitive help) is to press the F1 key, click a Help button, or right-click
the mouse and select Help from the pop-up menu. Context-sensitive help pro-
vides instant access to the help topics related to the current dialog.
How Do I?
The How Do I? section of the Help system has step-by-step instructions for
common applications of the model. You may find it beneficial to read over
these topics before using the software or at least to familiarize yourself with
them. When a related question arises, you can quickly find the answer.
Graphical Editor
Network models such as StormCAD, WaterGEMS, SewerCAD, and
PondPack have a graphical editor to aid the engineer in laying out the system.
The Bentley OpenFlows graphical editor is an “intelligent” drafting environ-
ment, meaning that it recognizes the characteristics and behavior of different
modeling elements and maintains connectivity when the user drags and drops
elements on-screen.
This intelligent graphical editor prevents problems such as incorrect pipe con-
nections and “hanging” pipes (pipes with a node defined at only one end).
Models that do not have this type of behavior can be easily “broken,” and it
may take hours (or days) to locate the source of the problem.
DXF Files
A drawing exchange file (DXF) is the standard format used for translating from
one CAD system, such as AutoCAD®, to another, such as MicroStation™.
DXF files can also be imported into StormCAD, WaterGEMS, and SewerCAD
for use as background drawings. For example, a roadway network plan can
be imported and used as a background map for laying out and designing a
proposed storm sewer system, sanitary sewer system, or water distribution
network.
If desired, drawing elements from DXF files can be automatically converted
to network components using the polyline-to-pipe import feature. DXF lines
and polylines can be converted to pipes, and DXF blocks can be converted to
node elements.
To learn more about importing and using DXF backgrounds, look in the
online help system of StormCAD, WaterGEMS, or SewerCAD. You can find
topics in the How Do I? section or perform a search with a keyword such as
“DXF” or “import.”
Scaled or Schematic?
In StormCAD, WaterGEMS, and SewerCAD, you may choose to either create
a schematic drawing or define the horizontal and vertical distance scale. The
scaled mode is typically used when a DXF or CAD background file is present.
The scale is set so that one drawing unit is equal to a unit of measure (e.g., 1
DXF unit equals 1 foot). In a scaled plan, you can also set the drawing scale
that you want to use when printing or plotting.
A schematic plan is not drawn at any defined scale, so there is no connection
between dimensions shown on the plan and real-world dimensions. There may
not even be a correlation between individual dimensions shown on a schematic
plan. For example, in a schematic plan, the length representing a pipe with a
length of 50 m may be less than, equal to, or greater than, the line representing
a pipe with a length of 20 m.
To get to the drawing scale options, choose Options from the Tools menu and
go to the Drawing tab. The available options are as follows:
• Schematic—You enter the pipe length in the Pipe Properties dialog for
each pipe. Absolute positioning in the drawing editor is not used.
• Scaled—You enter the drawing unit scale to be used to determine the
length of the drawn pipes. The pipe length will be calculated based on the
starting and ending coordinates. The length in the drawing editor corre-
sponds to a real-world length. You can also enter the horizontal and vertical
drawing scales at which output should be plotted.
Color Coding
One way to determine the performance of a system is to color code the sys-
tem according to some indicator. For hydraulic design work, pipeline veloc-
ity is often a useful indicator. For example, pipes consistently flowing below
0.15 m/s may be oversized. Pipes with velocities over 1.5 m/s are fairly heavily
stressed, and those above 2.4 m/s are often bottlenecks in the system.
FlexUnits
You can set field options—the properties of a numeric value—from almost
any field in any of Bentley OpenFlows’ programs. These properties include the
units, decimal precision, scientific notation display option, and the allowable
range of values. This feature is referred to as FlexUnits.
To set these options for a particular field type, simply right-click the field and
select Properties from the context menu. This action opens the Set Field
Properties dialog where you can select the options that you want.
If you would like to set up the units for many data types from the same loca-
tion in the program, choose FlexUnits from the Tools pull-down menu to
open the FlexUnits dialog.
Units
FlexUnit offers a wide variety of possible units for any field. In the Set Field
Options dialog, click the Units drop-down menu and select the desired unit
from the list. The value in the field will be automatically converted to the new
unit using the appropriate conversion factor. The FlexUnit feature allows you
to mix and match any units, even if you are using a combination of English
and SI units.
Once you begin changing units, you will probably notice that some units have
multiple representations (e.g., psi and lbs/in2, cfs and ft3/s). Computationally,
these units are the same; however, the additional choices allow you to control
the way the units are presented on-screen and in program output.
Display Precision
Display precision can be used to control the number of digits displayed after
a decimal point or the rounding of a numeric value. When the precision is
greater than or equal to zero, it specifies the number of digits that are displayed
after the decimal point. For example, π (3.14159265…) with a decimal preci-
sion of 4 would be presented as 3.1416.
Specifying a negative number for display precision results in a displayed value
that is rounded to the nearest power of 10. −1 rounds to 10, −2 rounds to 100,
−3 rounds to 1000, and so on. For example, the number 1,234 with a display
precision of −2 would be displayed as 1,200.
Note: The display precision affects only the way the numbers are displayed, not
their actual values. The internal values stored by the software are still carried
out to their maximum decimal precision—they are just displayed differently.
This concept is important concept to keep in mind, especially when check-
ing calculations by hand or working with values that are lower than the usual
decimal precision. For example, be aware that a 0.75-in diameter pipe with a
display precision of 0 will result in a displayed value of 1 in, even though the
calculations will be performed based on the true 0.75-in diameter.
Scientific Notation
Scientific notation displays a value as a real number multiplied by some power
of 10. It is displayed as an integer or real value followed by the letter “e” and a
positive or negative integer.
For example, 12,345 could be written in scientific notation as 1.2345 e4, and
0.12345 could be written as 1.2345 e–1.
Scientific notation follows the same display precision rules outlined previously.
To turn scientific notation on or off, just check or uncheck the box labeled
Scientific notation.
B.1 Overview
The Scenario Manager is included in many of the Bentley Systems software. This
feature can dramatically increase your productivity in the “What If?” areas of
modeling, including calibration, operations analysis, and planning. By invest-
ing a little time to understand scenario management, you can avoid unneces-
sary editing and data duplication. Tutorial 2 with Chapter 3 demonstrates how
to manage and compare scenarios using SewerCAD and SewerGEMS, which
can also be used for WaterGEMS, HAMMER, and PondPack.
In contrast to the old methods of scenario management (editing or copying
data), automated scenario management using inheritance gives you these
advantages:
• A single project file makes it possible to generate an unlimited number
of “What If?” conditions without becoming overwhelmed with numerous
modeling files and separate results.
• Because the software maintains the data for all the scenarios in a single pro-
ject, it can provide you with powerful automated tools for directly compar-
ing scenario results. Any set of results is immediately available at any time.
• The Scenario/Alternative relationship empowers you to mix and match
groups of data from existing scenarios without having to declare any data
again.
• With inheritance, you do not have to re-enter data if it remains unchanged
in a new alternative or scenario, avoiding redundant copies of the same
data. Inheritance also enables you to correct a data input error in a par-
ent scenario and automatically update the corrected attribute in all child
scenarios.
These advantages, while obvious, may not seem compelling for small projects.
It is as projects grow to hundreds or thousands of network elements that the
advantages of true scenario inheritance become clear. On a large project, being
able to maintain a collection of base and modified alternatives accurately and
• Maximize the number of scenarios you can simulate by mixing and match-
ing data from existing scenarios (data reuse).
• Minimize the amount of data that needs to be duplicated to consider con-
ditions that have a lot in common.
A single project file enables you to generate an unlimited number of “What If?”
conditions, edit only the data that needs to be changed, and quickly generate
direct comparisons of input and results for desired scenarios.
A Familiar Parallel
Although the structure of scenarios may seem a bit difficult at first, anyone
who has eaten at a restaurant should be able to relate fairly easily. A meal
(scenario) is comprised of several courses (alternatives), which might include a
salad, an entrée, and a dessert. Each course has its own attributes. For example,
the entrée may have a meat, a vegetable, and a starch. Examining the choices,
we could present a menu as in Figure B.2.
The restaurant does not have to create a new recipe for every possible meal
(combination of courses) that could be ordered. They can just assemble any
meal based on what the customer orders for each alternative course. Salad 1,
Entrée 1, and Dessert 2 might then be combined to define a complete meal.
Generalizing this concept, we see in Figure B.3 that any scenario simply refer-
ences one alternative from each category to create a “big picture” that can be
analyzed. Note that different types of alternatives may have different numbers
and types of attributes, and any category can have an unlimited number of
alternatives to choose from.
Overriding Inheritance
Overriding inheritance is the software equivalent of cosmetics. A child can
override inherited characteristics at any time by specifying a new value for that
characteristic. These overriding values do not affect the parent, and are there-
fore considered “local” to the child. Local values can also be removed at any
time, reverting the characteristic to its inherited state. The child has no choice
in the value of his inherited attributes, only in local attributes.
For example, suppose a child has inherited the attribute of blue eyes from his
parent. Now the child puts on a pair of greentinted contact lenses to hide his
natural eye color. When the contact lenses are on, we say his natural eye color
is “overridden” locally, and his eye color is green. When the child removes
the tinted lenses, his eye color instantly reverts to blue, as inherited from his
parent.
Dynamic Inheritance
Dynamic inheritance does not have a parallel in the genetic world. When
a parent’s characteristic is changed, existing children also reflect the change.
Using the eye-color example, this would be the equivalent of the parent chang-
ing eye color from blue to brown, and the children’s eyes instantly inheriting
the brown color also. Of course, if the child has already overridden a character-
istic locally, as with the green lenses, his eyes will remain green until the lenses
are removed. At this point, his eye color will revert to the inherited color, now
brown.
This dynamic inheritance has remarkable benefits for applying wide-scale
changes to a model, fixing an error, and so on. If rippling changes are not
desired, the child can override all of the parent’s values, or a copy of the parent
can be made instead of a child.
• “Entrée 2 is just like Entrée 1, except for the meat and the starch.”
• “Entrée 3 is just like Entrée 2, except for the meat.”
Note: If the vegetable of the day changes (say from green beans to peas), only
Entrée 1 needs to be updated, and the other entrées will automatically inherit
the vegetable attribute of “Peas” instead of “Green Beans”.
Note: Dessert 3 is a parent alternative and has nothing in common with the
other desserts, so it can be created as a “root” or “base” alternative. It does not
inherit its attribute data from any other alternative.
• “Meal 2 is just like Meal 1, except for the dessert.” The salad and entrée
alternatives are inherited from Meal 1.
• “Meal 3 is nothing like Meal 1 or Meal 2.” A totally new “base” or “root”
is created.
• “Meal 4 is just like Meal 3, except for the salad.” The entrée and dessert
alternatives are inherited from Meal 3.
Now we can create a child scenario from “Average Day” that inherits the phys-
ical alternative, but overrides the selected demand alternative. As a result, we
get the following scenario hierarchy:
Because no physical data (pipe diameters) have been changed, the physical
alternative hierarchy remains the same as before.
Another scenario is also created to reference these new demands, as shown here:
Note again that the physical data was not changed, so the physical alternatives
remain the same.
Correcting an Error
This analysis results in acceptable pressures, until it is discovered that the
industrial demand is not actually 500 gpm — it is 1,500 gpm! Because of the
inheritance within the demand alternatives, however, only the “Average Day”
demand for J-2 needs to be updated. The changes will ripple through to the
children. After the single change is made, the demand hierarchy is as follows:
Notice that no changes need to be made to the scenarios to reflect these cor-
rections. The three scenarios can now be calculated as a batch to update the
results.
When these results are reviewed, it is determined that the system does not have
the ability to adequately supply the system as it was originally thought. The
pressure at J-2 is too low under peak hour demand conditions.
This time, the demand alternative hierarchy remains the same since no
demands were changed. The two new scenarios (“Peak, Big P-1,” “Peak, All Big
Pipes”) can be batch run to provide results for these proposed improvements.
Also note that it would be equally effective in this case to inherit the “Avg.
Day, Big P-1” scenario from “Avg. Day” (changing the physical alternative) or
to inherit from “Peak, Big P-1” (changing the demand alternative). Likewise,
“Max. Day, Big P-1” could inherit from either “Max. Day” or “Peak, Big P-1.”
Neither the demand nor physical alternative hierarchies were changed in order
to run the last set of scenarios, so they remain as they were.
Inlet design overview, 81, 90–91 Muskingum–Cunge method, 246 rational method, 56–60
inlet and gutter network problems Muskingum routing method, time of concentration, 51–55
using StormCAD, 91 245–246 Pipe leakage, 163
Inlets in sag, 82, 88–89 Pipe networks, 170–171
Inlets on grade, 82, 84–85 Pipe outlets, 141
combination, 87–88 N Pond bottoms, 132
curb, 85–87 Nappe, 20 PondPack
Intensity–duration–frequency curves Navier–Stokes equations, 199 and tutorials, 155–156
(IDF curves), 30–32, 98 Network analysis, 171–172 for windows, 156
Interception, 39 Nonrevenue water and leaks, Precipitation, 40
Interconnected pond modeling, 130 163–164 Pressure-Reducing Valves (PRVs), 169
International System of Units (SI), 1 Nonuniform flow, 4 Pressure-Sustaining Valves
Normal depth, 4, 98 (PSVs), 169
NRCS curve number equation, 38 Pressure Breaker Valves (PBVs), 169
J NRCS (SCS) curve number Pressure flow, 3, 12, 14
Jet, 17 method, 46–47 Pressure of transients, 200–202
antecedent moisture condition, Pressure pipe friction losses, 12
48 Pressure piping systems, 161
K cover type, 48 automated optimization, 176–178
Kinematic wave model, 244 curve number tables, 48–51 control valves, 169
Kutter’s equation. See Chézy’s hydrologic condition, 48 energy gains, 166–168
equation soil group, 47 energy losses, 164–166
NRCS (SCS) peak flow network analysis, 171–172
estimation, 60–63 pipe networks, 170–171
L NRCS synthetic unit problems, 179–193
Lag, 130 hydrographs, 70–73 water demands, 161–164
Lagrangian TDM, 175 WaterGEMS and tutorials,
Laminar flow, 4–6 178–179
Land use, 42–43 O water quality analysis, 172–176
Loading, 212 Off-line pond, 129 Pressure wave celerity, 200
EFF, 215–217 Online pond, 129 Prismatic channel, 2
hydrographs, 213–214 Open-channel flow, 3, 14 Profile
load types, 213 critical flow, 14–17 classifications, 93–94
pattern loads, 214–215 specific energy, 14–17 in StormCAD, 98–99
sanitary unit loading, 215 Orifices, 17–19, 141 Proportional methods, 164
variable PF, 217–219 Oscillating pressure wave, 197 Pumps, 166
Low-Impact Development (LID), Outflow hydrograph, 130 constant horsepower pumps, 168
236–238 Outlet control hydraulics, 114–116 variable-speed pumps, 167–168
Low-level outlets, 140 Outlet structures, 132
Overflow spillways, 133, 141
R
M Rainfall, 29
Manning’s equation, 9–10, 24, 55 P abstractions, 39, 40
Mass conservation, 45 Pattern loads, 214–215 characteristics, 29
Metered demand, 161 Patterns, 162, 214 data types, 30–38
Mild slope, 93 Peaking factor (PF), 216 Horton infiltration
Minor losses, 165–166 variable, 217–219 equation, 41–43
Mixed flow profiles, 95–96 Peak runoff flow rate, 51 hyetograph, 32
Model calibration, 176–177 NRCS (SCS) peak flow NRCS (SCS) curve number
Momentum equations, 198–199, 201 estimation, 60–63 method, 46–51