ADVANCED WATER

DISTRIBUTION MODELING
AND MANAGEMENT

Authors
Thomas M. Walski
Donald V. Chase
Dragan A. Savic
Walter Grayman
Stephen Beckwith
Edmundo Koelle

Contributing Authors
Scott Cattran, Rick Hammond, Kevin Laptos, Steven G. Lowry,
Robert F. Mankowski, Stan Plante, John Przybyla, Barbara Schmitz

Peer Review Board
Lee Cesario (Denver Water), Robert M. Clark (U.S. EPA),
Jack Dangermond (ESRI), Allen L. Davis (CH2M Hill),
Paul DeBarry (Borton-Lawson), Frank DeFazio (Franklin G. DeFazio Corp.),
Kevin Finnan (Bristol Babcock), Wayne Hartell (Bentley Systems),
Brian Hoefer (ESRI), Bassam Kassab (Santa Clara Valley Water District),
James W. Male (University of Portland), William M. Richards
(WMR Engineering), Zheng Wu (Bentley Systems ),
and E. Benjamin Wylie (University of Michigan)



Click here to visit the Bentley Institute
Press Web page for more information

C H A P T E R
1
Introduction to Water Distribution
Modeling
Water distribution modeling is the latest technology in a process of advancement that
began two millennia ago when the Minoans constructed the first piped water convey-
ance system. Today, water distribution modeling is a critical part of designing and
operating water distribution systems that are capable of serving communities reliably,
efficiently, and safely, both now and in the future. The availability of increasingly
sophisticated and accessible models allows these goals to be realized more fully than
ever before.
This book is structured to take the engineer through the entire modeling process, from
gathering system data and understanding how a computer model works, through con-
structing and calibrating the model, to implementing the model in system design and
operations. The text is designed to be a first course for the novice modeler or
engineering student, as well as a reference for those more experienced with
distribution system simulations.
This chapter introduces the reader to water distribution modeling by giving an over-
view of the basic distribution system components, defining the nature and purposes of
distribution system simulations, and outlining the basic steps in the modeling process.
The last section of the chapter presents a chronology of advancements in water
distribution.
1.1 ANATOMY OF A WATER DISTRIBUTION SYSTEM
Although the size and complexity of water distribution systems vary dramatically,
they all have the same basic purposeto deliver water from the source (or treatment
facility) to the customer.
2 Introduction to Water Distribution Modeling Chapter 1
Sources of Potable Water
Untreated water (also called raw water) may come from groundwater sources or sur-
face waters such as lakes, reservoirs, and rivers. The raw water is usually transported
to a water treatment plant, where it is processed to produce treated water (also known
as potable or finished water). The degree to which the raw water is processed to
achieve potability depends on the characteristics of the raw water, relevant drinking
water standards, treatment processes used, and the characteristics of the distribution
system.
Before leaving the plant and entering the water distribution system, treated surface
water usually enters a unit called a clearwell. The clearwell serves three main pur-
poses in water treatment. First, it provides contact time for disinfectants such as
chlorine that are added near the end of the treatment process. Adequate contact time is
required to achieve acceptable levels of disinfection.
Second, the clearwell provides storage that acts as a buffer between the treatment
plant and the distribution system. Distribution systems naturally fluctuate between
periods of high and low water usage, thus the clearwell stores excess treated water
during periods of low demand and delivers it during periods of peak demand. Not
only does this storage make it possible for the treatment plant to operate at a more sta-
ble rate, but it also means that the plant does not need to be designed to handle peak
demands. Rather, it can be built to handle more moderate treatment rates, which
means lower construction and operational costs.
Third, the clearwell can serve as a source for backwash water for cleaning plant filters
that, when needed, is used at a high rate for a short period of time.
In the case of groundwater, many sources offer up consistently high quality water that
could be consumed without disinfection. However, the practice of maintaining a dis-
infectant residual is almost always adhered to for protection against accidental
contamination and microbial regrowth in the distribution system. Disinfection at
groundwater sources differs from sources influenced by surface water in that it is usu-
ally applied at the well itself.
Customers of Potable Water
Customers of a water supply system are easily identified they are the reason that
the system exists in the first place. Homeowners, factories, hospitals, restaurants, golf
courses, and thousands of other types of customers depend on water systems to pro-
vide everything from safe drinking water to irrigation. As demonstrated throughout
the book, customers and the nature in which they use water are the driving mechanism
behind how a water distribution system behaves. Water use can vary over time both in
the long-term (seasonally) and the short-term (daily), and over space. Good knowledge
of how water use is distributed across the system is critical to accurate modeling.
Transport Facilities
Moving water from the source to the customer requires a network of pipes, pumps,
valves, and other appurtenances. Storing water to accommodate fluctuations in
demand due to varying rates of usage or fire protection needs requires storage facili-
Section 1.1 Anatomy of a Water Distribution System 3
ties such as tanks and reservoirs. Piping, storage, and the supporting infrastructure are
together referred to as the water distribution system (WDS).
Transmission and Distribution Mains. This system of piping is often cate-
gorized into transmission/trunk mains and distribution mains. Transmission mains
consist of components that are designed to convey large amounts of water over great
distances, typically between major facilities within the system. For example, a trans-
mission main may be used to transport water from a treatment facility to storage tanks
throughout several cities and towns. Individual customers are usually not served from
transmission mains.
Distribution mains are an intermediate step toward delivering water to the end cus-
tomers. Distribution mains are smaller in diameter than transmission mains, and
typically follow the general topology and alignment of the city streets. Elbows, tees,
wyes, crosses, and numerous other fittings are used to connect and redirect sections of
pipe. Fire hydrants, isolation valves, control valves, blow-offs, and other maintenance
and operational appurtenances are frequently connected directly to the distribution
mains. Services, also called service lines, transmit the water from the distribution
mains to the end customers.
Homes, businesses, and industries have their own internal plumbing systems to trans-
port water to sinks, washing machines, hose bibbs, and so forth. Typically, the internal
plumbing of a customer is not included in a WDS model; however, in some cases,
such as sprinkler systems, internal plumbing may be modeled.
System Configurations. Transmission and distribution systems can be either
looped or branched, as shown in Figure 1.1. As the name suggests, in looped systems
there may be several different paths that the water can follow to get from the source to
a particular customer. In a branched system, also called a tree or dendritic system, the
water has only one possible path from the source to a customer.
Figure 1.1
Looped and branched
networks
Looped Branched
4 Introduction to Water Distribution Modeling Chapter 1
Looped systems are generally more desirable than branched systems because, coupled
with sufficient valving, they can provide an additional level of reliability. For exam-
ple, consider a main break occurring near the reservoir in each system depicted in
Figure 1.2. In the looped system, that break can be isolated and repaired with little
impact on customers outside of that immediate area. In the branched system, however,
all the customers downstream from the break will have their water service interrupted
until the repairs are finished. Another advantage of a looped configuration is that,
because there is more than one path for water to reach the user, the velocities will be
lower, and system capacity greater.
Figure 1.2
Looped and branched
networks after
network failure
Looped Branched
Customers
Without
Service
Customers
Without
Service
Pipe Break
Most water supply systems are a complex combination of loops and branches, with a
trade-off between loops for reliability (redundancy) and branches for infrastructure
cost savings. In systems such as rural distribution networks, the low density of cus-
tomers may make interconnecting the branches of the system prohibitive from both
monetary and logistical standpoints.
1.2 WHAT IS A WATER DISTRIBUTION SYSTEM
SIMULATION?
The term simulation generally refers to the process of imitating the behavior of one
system through the functions of another. In this book, the term simulation refers to the
process of using a mathematical representation of the real system, called a model.
Network simulations, which replicate the dynamics of an existing or proposed sys-
tem, are commonly performed when it is not practical for the real system to be
directly subjected to experimentation, or for the purpose of evaluating a system before
it is actually built. In addition, for situations in which water quality is an issue,
directly testing a system may be costly and a potentially hazardous risk to public
health.
Section 1.2 What Is a Water Distribution System Simulation? 5
Simulations can be used to predict system responses to events under a wide range of
conditions without disrupting the actual system. Using simulations, problems can be
anticipated in proposed or existing systems, and solutions can be evaluated before
time, money, and materials are invested in a real-world project.
For example, a water utility might want to verify that a new subdivision can be pro-
vided with enough water to fight a fire without compromising the level of service to
existing customers. The system could be built and tested directly, but if any problems
were to be discovered, the cost of correction would be enormous. Regardless of
project size, model-based simulation can provide valuable information to assist an
engineer in making well-informed decisions.
Simulations can either be steady-state or extended-period. Steady-state simulations
represent a snapshot in time and are used to determine the operating behavior of a sys-
tem under static conditions. This type of analysis can be useful in determining the
short-term effect of fire flows or average demand conditions on the system. Extended-
period simulations (EPS) are used to evaluate system performance over time. This
type of analysis allows the user to model tanks filling 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.
Modern simulation software packages use a graphical user interface (GUI) that
makes it easier to create models and visualize the results of simulations. Older-
generation software relied exclusively on tabular input and output. A typical modern
software interface with an annotated model drawing is shown in Figure 1.3.
Figure 1.3
Software interface and
annotated model
drawing
6 Introduction to Water Distribution Modeling Chapter 1
1.3 APPLICATIONS OF WATER DISTRIBUTION MODELS
Most water distribution models (WDMs) can be used to analyze a variety of other
pressure piping systems, such as industrial cooling systems, oil pipelines, or any net-
work carrying an incompressible, single-phase, Newtonian fluid in full pipes. Munic-
ipal water utilities, however, are by far the most common application of these models.
Models are especially important for WDSs due to their complex topology, frequent
growth and change, and sheer size. It is not uncommon for a system to supply hun-
dreds of thousands of people (large networks supply millions); thus, the potential
impact of a utility decision can be tremendous.
Water distribution network simulations are used for a variety of purposes, such as
Long-range master planning, including both new development and rehabili-
tation
Fire protection studies
Water quality investigations
Energy management
System design
Daily operational uses including operator training, emergency response, and
troubleshooting
Long-Range Master Planning
Planners carefully research all aspects of a water distribution system and try to deter-
mine which major capital improvement projects are necessary to ensure the quality of
service for the future. This process, called master planning (also referred to as capital
improvement planning or comprehensive planning), may be used to project system
growth and water usage for the next 5, 10, or 20 years. System growth may occur
because of population growth, annexation, acquisition, or wholesale agreements
between water supply utilities. The capability of the hydraulic network to adequately
serve its customers must be evaluated whenever system growth is anticipated.
Not only can a model be used to identify potential problem areas (such as future low
pressure areas or areas with water quality problems), but it can also be used to size
and locate new transmission mains, pumping stations, and storage facilities to ensure
that the predicted problems never occur. Maintaining a system at an acceptable level of
service is preferable to having to rehabilitate a system that has become problematic.
Rehabilitation
As with all engineered systems, the wear and tear on a water distribution system may
lead to the eventual need to rehabilitate portions of the system such as pipes, pumps,
valves, and reservoirs. Pipes, especially older, unlined, metal pipes, may experience
an internal buildup of deposits due to mineral deposits and chemical reactions within
the water. This buildup can result in loss of carrying capacity, reduced pressures, and
Section 1.3 Applications of Water Distribution Models 7
poorer water quality. To counter these effects of aging, a utility may choose to clean
and reline a pipe. Alternatively, the pipe may be replaced with a new (possibly larger)
pipe, or another pipe may be installed in parallel. Hydraulic simulations can be used
to assess the impacts of such rehabilitation efforts, and to determine the most econom-
ical improvements.
Fire Protection Studies
Water distribution systems are often required to provide water for fire fighting pur-
poses. Designing the system to meet the fire protection requirements is essential and
normally has a large impact on the design of the entire network. The engineer deter-
mines the fire protection requirements and then uses a model to test whether the
system can meet those requirements. If the system cannot provide certain flows and
maintain adequate pressures, the model may also be used for sizing hydraulic
elements (pipes, pumps, etc.) to correct the problem.
Water Quality Investigations
Some models provide water quality modeling in addition to hydraulic simulation
capabilities. Using a water quality model, the user can model water age, source trac-
ing, and constituent concentration analyses throughout a network. For example, chlo-
rine residual maintenance can be studied and planned more effectively, disinfection
by-product formation (DBP) in a network can be analyzed, or the impact of storage
tanks on water quality can be evaluated. Water quality models are also used to study
the modification of hydraulic operations to improve water quality.
Energy Management
Next to infrastructure maintenance and repair costs, energy usage for pumping is the
largest operating expense of many water utilities (Figure 1.4). Hydraulic simulations
can be used to study the operating characteristics and energy usage of pumps, along
with the behavior of the system. By developing and testing different pumping strate-
gies, the effects on energy consumption can be evaluated, and the utility can make an
educated effort to save on energy costs.
Daily Operations
Individuals who operate water distribution systems are generally responsible for mak-
ing sure that system-wide pressures, flows, and tank water levels remain within
acceptable limits. The operator must monitor these indicators and take action when a
value falls outside the acceptable range. By turning on a pump or adjusting a valve,
for example, the operator can adjust the system so that it functions at an appropriate
level of service. A hydraulic simulation can be used in daily operations to determine
the impact of various possible actions, providing the operator with better information
for decision-making.
8 Introduction to Water Distribution Modeling Chapter 1
Figure 1.4
Pumping is one of the
largest operating
expenses of many
utilities
Operator Training. Most water distribution system operators do their jobs very
well. As testimony to this fact, the majority of systems experience very few water out-
ages, and those that do occur are rarely caused by operator error. Many operators,
however, gain experience and confidence in their ability to operate the system only
over a long period of time, and sometimes the most critical experience is gained under
conditions of extreme duress. Hydraulic simulations offer an excellent opportunity to
train system operators in how their system will behave under different loading condi-
tions, with various control strategies, and in emergency situations.
Emergency Response. Emergencies are a very real part of operating a water
distribution system, and operators need to be prepared to handle everything from main
breaks to power failures. Planning ahead for these emergencies by using a model may
prevent service from being compromised, or may at least minimize the extent to
which customers are affected. Modeling is an excellent tool for emergency response
planning and contingency.
System Troubleshooting. When hydraulic or water quality characteristics in
an existing system are not up to standard, a model simulation can be used to identify
probable causes. A series of simulations for a neighborhood that suffers from chronic
low pressure, for example, may point toward the likelihood of a closed valve in the
area. A field crew can then be dispatched to this area to check nearby valves.
1.4 THE MODELING PROCESS
Assembling, calibrating, and using a water distribution system model can seem like a
foreboding task to someone confronted with a new program and stacks of data and
maps of the actual system. As with any large task, the way to complete it is to break it
down into its components and work through each step. Some tasks can be done in par-
allel while others must be done in series. The tasks that make up the modeling process
are illustrated in Figure 1.5. Note that modeling is an iterative process.
Section 1.4 The Modeling Process 9
Figure 1.5
Flowchart of the
modeling process
10 Introduction to Water Distribution Modeling Chapter 1
The first step in undertaking any modeling project is to develop a consensus within
the water utility regarding the need for the model and the purposes for which the
model will be used in both the short- and long-term. It is important to have utility per-
sonnel, from upper management and engineering to operations and maintenance,
commit to the model in terms of human resources, time, and funding. Modeling
should not be viewed as an isolated endeavor by a single modeler, but rather a utility-
wide effort with the modeler as the key worker. After the vision of the model has been
accepted by the utility, decisions on such issues as extent of model skeletonization
and accuracy of calibration will naturally follow.
Figure 1.5 shows that most of the work in modeling must be done before the model
can be used to solve real problems. Therefore, it is important to budget sufficient time
to use the model once it has been developed and calibrated. Too many modeling
projects fall short of their goals for usage because the model-building process takes up
all of the allotted time and resources. There is not enough time left to use the model to
understand the full range of alternative solutions to the problems.
Modeling involves a series of abstractions. First, the real pipes and pumps in the sys-
tem are represented in maps and drawings of those facilities. Then, the maps are
converted to a model that represents the facilities as links and nodes. Another layer of
abstraction is introduced as the behaviors of the links and nodes are described
mathematically. The model equations are then solved, and the solutions are typically
displayed on maps of the system or as tabular output. A models value stems from the
usefulness of these abstractions in facilitating efficient design of system improve-
ments or better operation of an existing system.
1.5 A BRIEF HISTORY OF WATER DISTRIBUTION
TECHNOLOGY
The practice of transporting water for human consumption has been around for sev-
eral millennia. From the first pipes in Crete some 3,500 years ago, to todays complex
hydraulic models, the history of water distribution technology is quite a story. The fol-
lowing highlights some of the key historical events that have shaped the field since its
beginnings.
1500 B.C. First water distribution pipes used in Crete. The Minoan civilization
flourishes on the island of Crete. The City of Knossos develops an aqueduct system
that uses tubular conduits to convey water. Other ancient civilizations have had sur-
face water canals, but these are probably the first pipes.
250 B.C. Archimedes principle developed. Archimedes, best known for his dis-
covery of S and for devising exponents, develops one of the earliest laws of fluids
when he notices that any object in water displaces its own volume. Using this princi-
ple, he proves that a crown belonging to King Hiero of Syracuse is not made of gold.
A legend will develop that he discovered this principle while bathing and became so
excited that he ran naked through the streets shouting Eureka (Ive found it).
100 A.D. Roman aqueducts. The Romans bring water from great distances to
their cities through aqueducts (Figure 1.6). While many of the aqueducts are above-
Section 1.5 A Brief History of Water Distribution Technology 11
ground, there are also enclosed conduits to supply public fountains and baths. Sextus
Julius Frontinus, water commissioner of Rome, writes two books on the Roman water
supply.
Figure 1.6
Roman aqueduct
1455 First cast iron pipe. Casting of iron for pipe becomes practical, and the first
installation of cast iron pipe, manufactured in Siegerland, Germany, occurs at Dillen-
burg Castle.
1652 Piped water in Boston. The first water pipes in the U.S. are laid in Boston to
bring water from springs to what is now the Quincy Market area.
1664 Palace of Versailles. King Louis XIV of France orders the construction of a
15-mile cast iron water main from Marly-on-Seine to the Palace of Versailles. This is
the longest pipeline of its kind at this time, and portions of it remain in service into the
21st century. A section of the line, after being taken out of service, was shipped in the
1960s from France to the United States (Figure 1.7) where it is still on display.
Figure 1.7
King Louis XIV of
France and a section
of the Palace of
Versailles pipeline
Courtesy of the Ductile Iron Pipe Research Association
12 Introduction to Water Distribution Modeling Chapter 1
1732 Pitot invents a velocity-measuring device. Henri Pitot is tasked with mea-
suring the velocity of water in the Seine River. He finds that by placing an L-shaped
tube into the flow, water rises in the tube proportionally to the velocity squared, and
the Pitot tube is born.
1738 Bernoulli publishes Hydrodynamica. The Swiss Bernoulli family extends
the early mathematics and physics discoveries of Newton and Leibniz to fluid sys-
tems. Daniel Bernoulli publishes Hydrodynamica while in St. Petersburg and
Strasbourg, but there is a rivalry with his father Johann regarding who actually
developed some of the principles presented in the book. These principles will become
the key to energy principles used in hydraulic models and the basis for numerous
devices such as the Venturi meter and, most notably, the airplane wing. In 1752,
however, it will actually be their colleague, Leonard Euler, who develops the forms of
the energy equations that will live on in years to come.
1754 First U.S. water systems built. The earliest water distribution systems in the
United States are constructed in Pennsylvania. The Moravian community in Bethle-
hem, Pennsylvania claims to have the first water system, and it is followed quickly by
systems in Schaefferstown and Philadelphia, Pennsylvania. Horses drive the pumps in
the Philadelphia system, and the pipes are made of bored logs. They will later be
replaced with wood stave pipes made with iron hoops to withstand higher pressures.
The first steam driven pumps will be used in Bethlehem ten years later.
1770 Chezy develops head loss relationship. While previous investigators real-
ized that energy was lost in moving water, it is Antoine Chezy who realizes that V
2
/
RS is reasonably constant for certain situations. This relationship will serve as the
basis for head loss equations to be used for centuries.
1785 Bell and spigot joint developed. The Chelsea Water Company in London
begins using the first bell and spigot joints. The joint is first packed with yarn or hemp
and is then sealed with lead. Sir Thomas Simpson is credited with inventing this joint,
which replaced the crude flanged joints used previously.
1839 Hagen-Poiseuille equation developed. Gotthilf Hagen and Jean Louis Poi-
seuille independently develop the head loss equations for laminar flow in small tubes.
Their work is experimental, and it is not until 1856 that Franz Neuman and Eduard
Hagenbach will theoretically derive the Hagen-Poiseuille equation.
1843 St. Venant develops equations of motion. Several researchers, including
Louis Navier, George Stokes, Augustin de Cauchy, and Simeon Poisson, work toward
the development of the fundamental differential equations describing the motion of
fluids. They become known as the Navier-Stokes equations. Jean-Claude Barre de
Saint Venant develops the most general form of these equations, but the term St.
Venant equations will be used to refer to the vertically and laterally averaged (that is,
one-dimensional flow) form of equations.
1845 Darcy-Weisbach head loss equation developed. Julius Weisbach publishes
a three-volume set on engineering mechanics that includes the results of his experi-
ments. The Darcy-Weisbach equation comes from this work, which is essentially an
extension of Chezys work, as Chezys C is related to Darcy-Weisbachs f by C
2
=8g/f.
Section 1.5 A Brief History of Water Distribution Technology 13
Darcys name is also associated with Darcys law for flow through porous media,
widely used in groundwater analysis.
1878 First automatic sprinklers used. The first Parmelee sprinklers are installed.
These are the first automatic sprinklers for fire protection.
1879 Lambs Hydrodynamics published. Sir Horace Lamb publishes his Treatise
on the Mathematical Theory of the Motion of Fluids. Subsequent editions will be pub-
lished under the title Hydrodynamics, with the last edition published in 1932.
1881 AWWA formed. The 22 original members create the American Water Works
Association. The first president is Jacob Foster from Illinois.
1883 Laminar/turbulent flow distinction explained. While earlier engineers
such as Hagen observed the differences between laminar and turbulent flow, Osborne
Reynolds is the first to conduct the experiments that clearly define the two flow
regimes. He identifies the dimensionless number, later referred to as the Reynolds
number, for quantifying the conditions under which each type of flow exists. He pub-
lishes An Experimental Investigation of the Circumstances which Determine
whether the Motion of Water shall be Direct or Sinuous and the Law of Resistance in
Parallel Channels.
1896 Cole invents Pitot tube for pressure pipe. Although numerous attempts
were made to extend Henri Pitots velocity measuring device to pressure pipes,
Edward Cole develops the first practical apparatus using a Pitot tube with two tips
connected to a manometer. The Cole Pitometer will be widely used for years to come,
and Coles company, Pitometer Associates, will perform flow measurement studies
(among many other services) into the 21st century.
1906 Hazen-Williams equation developed. A. Hazen and G.S. Williams develop
an empirical formula for head loss in water pipes. Although not as general or precise
in rough, turbulent flow as the Darcy-Weisbach equation, the Hazen-Williams equa-
tion proves easy to use and will be widely applied in North America.
1900 1930 Boundary Layer Theory developed. The interactions between flu-
ids and solids are studied extensively by a series of German scientists lead by Ludwig
Prandtl and his students Theodor von Karman, Johan Nikuradse, Heinrich Blasius,
and Thomas Stanton. As a result of their research, they are able to theoretically
explain and experimentally verify the nature of drag between pipe walls and a fluid.
In particular, the experiments of Nikuradse, who glues uniform sand grains inside
pipes and measures head loss, lead to a better understanding of the calculation of the f
coefficient in the Darcy-Weisbach equation. Stanton develops the first graphical rep-
resentation of the relationship between f, pipe roughness, and the Reynolds number,
which later leads to the Moody diagram. This work is summarized in H. Schichtings
book, Boundary Layer Theory.
1914 First U.S. drinking water standards established. The U.S. Public Health
Service publishes the first drinking water standards, which will continually evolve.
The U.S. Environmental Protection Agency (U.S. EPA) will eventually assume the
role of setting the water quality standards in the United States.
14 Introduction to Water Distribution Modeling Chapter 1
1920s Cement-mortar lining of water mains. Cement mortar lining of water
mains is used to minimize corrosion and tuberculation. Procedures for cleaning and
lining existing pipes in place will be developed by the 1930s.
1921 First Hydraulic Institute Standards published. The first edition of Trade
Standards in the Pump Industry is published as a 19-page pamphlet. These standards
become the primary reference for pump nomenclature, testing, and rating.
1936 Hardy Cross method developed. Hardy Cross, a structural engineering pro-
fessor at the University of Illinois, publishes the Hardy Cross method for solving head
loss equations in complex networks. This method is widely used for manual calcula-
tions and will serve as the basis for early digital computer programs for pipe network
analysis.
1938 Colebrook-White equation developed. Cyril Colebrook and Cedric White
of Imperial College in London build upon the work of Prandtl and his students to
develop the Colebrook-White equation for determining the Darcy-Weisbach f in com-
mercial pipes.
1940 Hunter curves published. During the 1920s and 30s, Roy Hunter of the
National Bureau of Standards conducts research on water use in a variety of build-
ings. His fixture unit method will become the basis for estimating building water
use, even though plumbing fixtures will change over the years. His probabilistic anal-
ysis captured the mathematics of the concept that the more fixtures in a building, the
less likely they are to be used simultaneously.
1944 Moody diagram published. Lewis Moody of Princeton University publishes
the Moody diagram, which is essentially a graphical representation of the Colebrook-
White equation in the turbulent flow range and the Hagen-Poisseuille equation in the
laminar range. This diagram is especially useful because, at the time, no explicit solu-
tion exists for the Colebrook-White equation. Stanton had developed a similar chart
30 years earlier.
1950 McIlroy network analyzer developed. The McIlroy network analyzer, an
electrical analog computer, is developed to simulate the behavior of water distribution
systems using electricity instead of water. The analyzer uses special elements called
fluistors to reproduce head loss in pipes, because in the Hazen-Williams equation,
head loss varies with flow raised to the 1.85 power, while normal resistors comply
with Ohms law, in which voltage drop varies linearly with current.
1950s Earliest digital computers developed. The Electronic Numerical Integra-
tor and Computer (ENIAC) is assembled at the University of Pennsylvania. It con-
tains approximately 18,000 vacuum tubes and fills a 30 x 50 ft (9 x 15 m) room.
Digital computers such as the ENIAC and Univac show that computers can carry out
numerical calculations quickly, opening the door for programs to solve complex
hydraulic problems.
1956 Push-on joint developed. The push-on pipe joint using a rubber gasket is
developed. This type of assembly helps speed the construction of piping.
1960s and 70s Earliest pipe network digital models created. With the coming
of age of digital computers and the establishment of the FORTRAN programming
Section 1.5 A Brief History of Water Distribution Technology 15
language, researchers at universities begin to develop pipe network models and make
them available to practicing engineers. Don Wood at the University of Kentucky, Al
Fowler at the University of British Columbia, Roland Jeppson of Utah State Univer-
sity, Chuck Howard and Uri Shamir at MIT, and Simsek Sarikelle at the University of
Akron all write pipe network models.
Figure 1.8
A computer punch
card
1963 First U.S. PVC pipe standards. The National Bureau of Standards accepts
CS256-63 Commercial Standard for PVC Plastic Pipes (SDR-PR and Class T),
which is the first U.S. standard for polyvinyl chloride water pipe.
1963 URISA is founded. The Urban and Regional Information Systems Associa-
tion is founded by Dr. Edgar Horwood. URISA becomes the premier organization for
the use and integration of spatial information technology to improve the quality of life
in urban and regional environments.
1960s and 70s Water system contamination. Chemicals that can result in health
problems when ingested or inhaled are dumped on the ground or stored in leaky
ponds because of lack of awareness of their environmental impacts. Over the years,
these chemicals will make their way into water distribution systems and lead to
alleged contamination of water systems in places like Woburn, Massachusetts; Phoe-
nix/Scottsdale, Arizona; and Dover Township, New Jersey. Water quality models of
distribution systems will be used to attempt to recreate the dosages of chemicals
received by customers. These situations lead to popular movies like A Civil Action
and Erin Brockovich.
1970s Early attempts to optimize water distribution design. Dennis Lai and
John Schaake at MIT develop the first approach to optimize water system design.
Numerous papers will follow by researchers such as Arun Deb, Ian Goulter, Uri
Shamir, Downey Brill, Larry Mays, and Kevin Lansey.
1970s Models become more powerful. Although the earliest pipe network models
could only solve steady-state equations for simple systems, the 70s bring modeling
features such as pressure regulating valves and extended-period simulations.
16 Introduction to Water Distribution Modeling Chapter 1
1975 Data files replace input cards. Modelers are able to remotely create data
files on time-share terminals instead of using punched cards.
1975 AWWA C-900 approved. The AWWA approves its first standard for PVC
water distribution piping. C900 pipe is made to match old cast iron pipe outer diame-
ters.
1976 Swamee-Jain equation published. Dozens of approximations to the
Colebrook-White equations have been published in an attempt to arrive at an explicit
equation that would give the same results without the need for an iterative solution.
Indian engineers P. K. Swamee and Akalnank Jain publish the most popular form of
these approximations. The use of an explicit equation results in faster numerical solu-
tions of pipe network problems.
1976 Jeppson publishes Analysis of Flow in Pipe Networks. Roland Jeppson
authors the book Analysis of Flow in Pipe Networks, which presents a summary of the
numerical techniques used to solve network problems.
1980 Personal computers introduced. Early personal computers make it possible
to move hydraulic analysis to desktop systems. Initially, these desktop models are
slow, but their power will grow exponentially over the next two decades.
Figure 1.9
Time-share
terminal
Early 1980s Water Quality Modeling First Developed. The concept of modeling
water quality in distribution systems is first developed, and steady state formulations
are proposed by Don Wood at the University of Kentucky and USEPA researchers in
Cincinnati, Ohio.
Section 1.5 A Brief History of Water Distribution Technology 17
1985 Battle of the Network Models. A series of sessions is held at the ASCE
Water Resources Planning and Management Division Conference in Buffalo, New
York, where researchers are given a realistic system called Anytown and are asked
to optimize the design of that network. Comparison of results shows the strengths and
weaknesses of the various models.
1986 Introduction of Dynamic Water Quality Models. At the AWWA Distribu-
tion System Symposium, three groups independently introduce dynamic water quality
models of distribution systems.
1988 Gradient Algorithm. Ezio Todini and S. Pilati publish A Gradient Algo-
rithm for the Analysis of Pipe Networks, and R. Salgado, Todini, and P. O'Connell
publish Comparison of the Gradient Method with some Traditional Methods of the
Analysis of Water Supply Distribution Networks. The gradient algorithm serves as
the basis for the WaterCAD model.
1989 AWWA holds specialty conference. AWWA holds the Computers and Auto-
mation in the Water Industry conference. This conference will later grow into the
popular IMTech event (Information Management and Technology).
1990s Privatization of water utilities. The privatization of water utilities
increases significantly as other utilities experience a greater push toward deregulation.
1991 Water Quality Modeling in Distribution Systems Conference. The
USEPA and the AWWA Research Foundation bring together researchers from around
18 Introduction to Water Distribution Modeling Chapter 1
the world for a two-day meeting in Cincinnati. This meeting is a milestone in the
establishment of water quality modeling as a recognized tool for investigators.
1991 GPS technology becomes affordable. The cost of global positioning sys-
tems (GPS) drops to the point where a GPS can be an economical tool for determining
coordinates of points in hydraulic models.
1993 Introduction of water quality modeling tool. Water quality modeling
comes of age with the development of EPANET by Lewis Rossman of the USEPA.
Intended as a research tool, EPANET provides the basis for several commercial-grade
models.
1990 through present. Several commercial software developers release water distri-
bution modeling packages. Each release brings new enhancements for data manage-
ment and new abilities to interoperate with other existing computer systems.
2001 Automated calibration. Automated calibration of distribution models
moves from being a research tool to a standard modeling feature with the use of
Genetic Algorithms.
2001 Security awareness. Water system security increases in importance and util-
ities realize the value of water quality modeling as a tool for protecting a water system.
2002 Integration with GIS. Water modeling and GIS software become highly
integrated with the release of WaterGEMS, software that combines the functionality
of both tools.
What Next?
Predicting the future is difficult, especially with rapidly changing fields such as the
software industry. However, there are definite trends as data sharing continues to gain
popularity, modeling spreads into operations, and automated design tools add to the
modelers arsenal.
The next logical question is, When will network models eliminate the need for engi-
neers? The answer is, never. Though a word processor can reduce the number of
spelling and grammar mistakes, it cannot write a best-selling novel. Even as technol-
ogy advances, an essential need still exists for living, breathing, thinking human
beings. A network model is just another tool (albeit a very powerful, multi-purpose
tool) for an experienced engineer or technician. It is still the responsibility of the user
to understand the real system, understand the model, and make decisions based on
sound engineering judgement.
REFERENCES
Mays, L. W. (2000). Introduction. Water Distribution System Handbook, Mays, L. W., ed., McGraw Hill,
New York, New York.




ADVANCED WATER
DISTRIBUTION MODELING
AND MANAGEMENT

Authors
Thomas M. Walski
Donald V. Chase
Dragan A. Savic
Walter Grayman
Stephen Beckwith
Edmundo Koelle

Contributing Authors
Scott Cattran, Rick Hammond, Kevin Laptos, Steven G. Lowry,
Robert F. Mankowski, Stan Plante, John Przybyla, Barbara Schmitz

Peer Review Board
Lee Cesario (Denver Water), Robert M. Clark (U.S. EPA),
Jack Dangermond (ESRI), Allen L. Davis (CH2M Hill),
Paul DeBarry (Borton-Lawson), Frank DeFazio (Franklin G. DeFazio Corp.),
Kevin Finnan (Bristol Babcock), Wayne Hartell (Bentley Systems),
Brian Hoefer (ESRI), Bassam Kassab (Santa Clara Valley Water District),
James W. Male (University of Portland), William M. Richards
(WMR Engineering), Zheng Wu (Bentley Systems ),
and E. Benjamin Wylie (University of Michigan)



Click here to visit the Bentley Institute
Press Web page for more information

C H A P T E R
2
Modeling Theory
Model-based simulation is a method for mathematically approximating the behavior
of real water distribution systems. To effectively utilize the capabilities of distribution
system simulation software and interpret the results produced, the engineer or mod-
eler must understand the mathematical principles involved. This chapter reviews the
principles of hydraulics and water quality analysis that are frequently employed in
water distribution network modeling software.
2.1 FLUID PROPERTIES
Fluids can be categorized as gases or liquids. The most notable differences between
the two states are that liquids are far denser than gases, and gases are highly com-
pressible compared to liquids (liquids are relatively incompressible). The most impor-
tant fluid properties taken into consideration in a water distribution simulation are
specific weight, fluid viscosity, and (to a lesser degree) compressibility.
Density and Specific Weight
The density of a fluid is the mass of the fluid per unit volume. The density of water is
1.94 slugs/ft
3
(1000 kg/m
3
) at standard pressure of 1 atm (1.013 bar) and standard tem-
perature of 32.0
o
F (0.0
o
C). A change in temperature or pressure will affect the density,
although the effects of minor changes are generally insignificant for water modeling
purposes.
The property that describes the weight of a fluid per unit volume is called specific
weight and is related to density by gravitational acceleration:
g =
20 Modeling Theory Chapter 2
(2.1)
where = fluid specific weight (M/L
2
/T
2
)
= fluid density (M/L
3
)
g = gravitational acceleration constant (L/T
2
)
The specific weight of water, , at standard pressure and temperature is 62.4 lb/ft
3
(9,806 N/m
3
).
Viscosity
Fluid viscosity is the property that describes the ability of a fluid to resist deformation
due to shear stress. For many fluids, most notably water, viscosity is a proportionality
factor relating the velocity gradient to the shear stress, as described by Newtons law
of viscosity:
t
V d
y d
------ = (2.2)
where t = shear stress (M/L/T
2
)
= absolute (dynamic) viscosity (M/L/T)
dV
dy
------- = time rate of strain (1/T)
The physical meaning of this equation can be illustrated by considering the two paral-
lel plates shown in Figure 2.1. The space between the plates is filled with a fluid, and
the area of the plates is large enough that edge effects can be neglected. The plates are
separated by a distance y, and the top plate is moving at a constant velocity V relative
to the bottom plate. Liquids exhibit an attribute known as the no-slip condition, mean-
ing that they adhere to surfaces they contact. Therefore, if the magnitude of V and y
are not too large, then the velocity distribution between the two plates is linear.
From Newtons second law of motion, for an object to move at a constant velocity, the
net external force acting on the object must equal zero. Thus, the fluid must be exert-
ing a force equal and opposite to the force F on the top plate. This force within the
fluid is a result of the shear stress between the fluid and the plate. The velocity at
which these forces balance is a function of the velocity gradient normal to the plate
and the fluid viscosity, as described by Newtons law of viscosity.
Thick fluids, such as syrup and molasses, have high viscosities. Thin fluids, such as
water and gasoline, have low viscosities. For most fluids, the viscosity remains con-
stant regardless of the magnitude of the shear stress that is applied to it.
Returning to Figure 2.1, as the velocity of the top plate increases, the shear stresses in
the fluid increase at the same rate. Fluids that exhibit this property conform to New-
tons law of viscosity and are called Newtonian fluids. Water and air are examples of
Newtonian fluids. Some types of fluids, such as inks and sludge, undergo changes in
viscosity when the shear stress changes. Fluids exhibiting this type of behavior are
called pseudo-plastic fluids.
Section 2.1 Fluid Properties 21
Figure 2.1
Physical interpretation
of Newtons law of
viscosity
F
V
y
dV
dy
Relationships between the shear stress and the velocity gradient for typical Newtonian
and non-Newtonian fluids are shown in Figure 2.2. Since most distribution system
models are intended to simulate water, many of the equations used consider Newto-
nian fluids only.
Figure 2.2
Stress versus strain for
plastics and fluids
I
d
e
a
l
P
l
a
s
t
i
c
P
s
e
u
d
o
-
p
la
s
t
ic
F
lu
id
N
e
w
to
n
ia
n
F
lu
id
E
l
a
s
t
i
c
S
o
l
i
d
Ideal Fluid
dV/dy
t
=
S
h
e
a
r
S
t
r
e
s
s
Viscosity is a function of temperature, but this relationship is different for liquids and
gases. In general, viscosity decreases as temperature increases for liquids, and viscos-
ity increases as temperature increases for gases. The temperature variation within
22 Modeling Theory Chapter 2
water distribution systems, however, is usually quite small, and thus changes in water
viscosity are considered negligible for this application. Generally, water distribution
system modeling software treats viscosity as a constant [assuming a temperature of
68
o
F (20
o
C)].
The viscosity derived in Equation 2.2 is referred to as the absolute viscosity (or
dynamic viscosity). For hydraulic formulas related to fluid motion, the relationship
between fluid viscosity and fluid density is often expressed as a single variable. This
relationship, called the kinematic viscosity, is expressed as follows:
v

--- = (2.3)
where v = kinematic viscosity (L
2
/T)
Just as there are shear stresses between the plate and the fluid in Figure 2.1, there are
shear stresses between the wall of a pipe and the fluid moving through the pipe. The
higher the fluid viscosity, the greater the shear stresses that will develop within the
fluid, and, consequently, the greater the friction losses along the pipe. Distribution
system modeling software packages use fluid viscosity as a factor in estimating the
friction losses along a pipes length. Packages that can handle any fluid require the
viscosity and density to be input by the modeler, while models that are developed only
for water usually account for the appropriate value automatically.
Fluid Compressibility
Compressibility is a physical property of fluids that relates the volume occupied by a
fixed mass of fluid to its pressure. In general, gases are much more compressible than
liquids. An air compressor is a simple device that utilizes the compressibility of air to
store energy. The compressor is essentially a pump that intermittently forces air mole-
cules into the fixed volume tank attached to it. Each time the compressor turns on, the
mass of air, and therefore the pressure within the tank, increases. Thus a relationship
exists between fluid mass, volume, and pressure.
This relationship can be simplified by considering a fixed mass of a fluid. Compress-
ibility is then described by defining the fluids bulk modulus of elasticity:
E
v
V
dP
dV
------- = (2.4)
where E
v
= bulk modulus of elasticity (M/L/T
2
)
P = pressure (M/L/T
2
)
V = volume of fluid (L
3
)
All fluids are compressible to some extent. The effects of compression in a water dis-
tribution system are very small, and thus the equations used in hydraulic simulations
are based on the assumption that the liquids involved are incompressible. With a bulk
modulus of elasticity of 410,000 psi (2.83 10
6
kPa) at 68
o
F (20
o
C), water can safely
be treated as incompressible. For instance, a pressure change of over 2,000 psi
(1.379 10
4
kPa) results in only a 0.5 percent change in volume.
Section 2.2 Fluid Statics and Dynamics 23
Although the assumption of incompressibility is justifiable under most conditions,
certain hydraulic phenomena are capable of generating pressures high enough that the
compressibility of water becomes important. During field operations, a phenomenon
known as water hammer can develop due to extremely rapid changes in flow (when,
for instance, a valve suddenly closes, or a power failure occurs and pumps stop oper-
ating). The momentum of the moving fluid can generate pressures large enough that
fluid compression and pipe wall expansion can occur, which in turn causes destructive
transient pressure fluctuations to propagate throughout the network. Specialized net-
work simulation software is necessary to analyze these transient pressure effects. For
complete coverage of transient flow, see Chapter 13.
Vapor Pressure
Consider a closed container that is partly filled with water. The pressure in the con-
tainer is measured when the water is first added, and again after some time has
elapsed. These readings show that the pressure in the container increases during this
period. The increase in pressure is due to the evaporation of the water, and the result-
ing increase in vapor pressure above the liquid.
Assuming that temperature remains constant, the pressure will eventually reach a con-
stant value that corresponds to the equilibrium or saturation vapor pressure of water
at that temperature. At this point, the rates of evaporation and condensation are equal.
The saturation vapor pressure increases with increasing temperature. This relationship
demonstrates, for example, why the air in humid climates typically feels moister in
summer than in winter, and why the boiling temperature of water is lower at higher
elevations.
If a sample of water at a pressure of 1 atm and room temperature is heated to 212
o
F
(100
o
C), the water will begin to boil since the vapor pressure of water at that tempera-
ture is equal to 1 atm. In a similar vein, if water is held at a temperature of 68
o
F
(20
o
C), and the pressure is decreased to 0.023 atm, the water will also boil.
This concept can be applied to water distribution in cases in which the ambient pres-
sure drops very low. Pump cavitation occurs when the fluid being pumped flashes
into a vapor pocket and then quickly collapses. For this to happen, the pressure in the
pipeline must be equal to or less than the vapor pressure of the fluid. When cavitation
occurs, it sounds as if gravel is being pumped, and severe damage to pipe walls and
pump components can result. For complete coverage of cavitation, see Chapter 13.
2.2 FLUID STATICS AND DYNAMICS
Static Pressure
Pressure can be thought of as a force applied normal, or perpendicular, to a body that
is in contact with a fluid. In the English system of units, pressure is expressed in
pounds per square foot (lb/ft
2
), but the water industry generally uses lb/in
2
, typically
abbreviated as psi. In the SI system, pressure has units of N/m
2
, also called a Pascal.
24 Modeling Theory Chapter 2
However, because of the magnitude of pressures occurring in distribution systems,
pressure is typically reported in kilo-Pascals (kPa), or 1,000 Pascals.
Pressure varies with depth, as illustrated in Figure 2.3. For fluids at rest, the variation
of pressure over depth is linear and is called the hydrostatic pressure distribution.
P h = (2.5)
where P = pressure (M/L/T
2
)
h = depth of fluid above datum (L)
= fluid specific weight (M/L
2
/T
2
)
Figure 2.3
Static pressure in a
standing water
column
D
e
p
t
h
Pressure = (Depth) g
P
This equation can be rewritten to find the height of a column of water that can be sup-
ported by a given pressure:
h
P

--- = (2.6)
The quantity P/ is called the pressure head, which is the energy resulting from water
pressure. Recognizing that the specific weight of water in English units is 62.4 lb/ft
3
, a
convenient conversion factor can be established for water as 1 psi = 2.31 ft (1 kPa =
0.102 m) of pressure head.
Section 2.2 Fluid Statics and Dynamics 25
Example Pressure Calculation. Consider the storage tank in Figure 2.4 in which the
water surface elevation is 120 ft above a pressure gage. The pressure at the base of the tank is due to
the weight of the column of water directly above it, and can be calculated as follows:
P h
62.4
lb
ft
3
------ 120ft ( )
144
in
2
ft
2
-------
----------------------------------- = =
P 52 = psi
Figure 2.4
Storage tank
120 ft
P = 52 psi
base
Absolute Pressure and Gage Pressure. Pressure at a given point is due to
the weight of the fluid above that point. The weight of the earths atmosphere pro-
duces a pressure, referred to as atmospheric pressure. Although the actual atmo-
spheric pressure depends on elevation and weather, standard atmospheric pressure at
sea level is 1 atm (14.7 psi or 101 kPa).
Two types of pressure are commonly used in hydraulics: absolute pressure and gage
pressure. Absolute pressure is the pressure measured with absolute zero (a perfect
vacuum) as its datum, and gage pressure is the pressure measured with atmospheric
pressure as its datum. The two are related to one another as shown in Equation 2.7 and
as illustrated in Figure 2.5. Note that when a pressure gage located at the earths sur-
face is open to the atmosphere, it registers zero on its dial. If the gage pressure is neg-
ative (that is, the pressure is below atmospheric), then the negative pressure is called a
vacuum.
P
abs
P
gage
P
atm
+ =
26 Modeling Theory Chapter 2
(2.7)
where P
abs
= absolute pressure (M/L/T
2
)
P
gage
= gage pressure (M/L/T
2
)
P
atm
= atmospheric pressure (M/L/T
2
)
In most hydraulic applications, including water distribution systems analysis, gage
pressure is used. Using absolute pressure has little value, since doing so would simply
result in all the gage pressures being incremented by atmospheric pressure. Addition-
ally, gage pressure is often more intuitive because people do not typically consider
atmospheric effects when thinking about pressure.
Figure 2.5
Gage versus absolute
pressure
H = 20 ft
P = 0 psi
gage
P = 14.6 psi
abs
P
gage
= 8.7 psi
P
abs
= 23.3 psi
Water
Velocity and Flow Regime
The velocity profile of a fluid as it flows through a pipe is not constant across the
diameter. Rather, the velocity of a fluid particle depends on where the fluid particle is
located with respect to the pipe wall. In most cases, hydraulic models deal with the
average velocity in a cross-section of pipeline, which can be found by using the fol-
lowing formula:
V
Q
A
---- = (2.8)
where V = average fluid velocity (L/T)
Q = pipeline flow rate (L
3
/T)
A = cross-sectional area of pipeline (L
2
)
The cross-sectional area of a circular pipe can be directly computed from the diameter
D, so the velocity equation can be rewritten as:
V
4Q
tD
2
---------- = (2.9)
where D = diameter (L)
Section 2.2 Fluid Statics and Dynamics 27
For water distribution systems in which diameter is measured in inches and flow is
measured in gallons per minute, the equation simplifies to:
V 0.41
Q
D
2
------ = (2.10)
where V = average fluid velocity (ft/s)
Q = pipeline flow rate (gpm)
D = diameter (in.)
Reynolds Number. In the late 1800s, an English scientist named Osborne Rey-
nolds conducted experiments on fluid passing through a glass tube. His experimental
setup looked much like the one in Figure 2.6 (Streeter, Wylie, and Bedford, 1998).
The experimental apparatus was designed to establish the flow rate through a long
glass tube (meant to simulate a pipeline) and to allow dye (from a smaller tank) to
flow into the liquid. He noticed that at very low flow rates, the dye stream remained
intact with a distinct interface between the dye stream and the fluid surrounding it.
Reynolds referred to this condition as laminar flow. At slightly higher flow rates, the
dye stream began to waver a bit, and there was some blurring between the dye stream
and the surrounding fluid. He called this condition transitional flow. At even higher
flows, the dye stream was completely broken up, and the dye mixed completely with
the surrounding fluid. Reynolds referred to this regime as turbulent flow.
When Reynolds conducted the same experiment using different fluids, he noticed that
the condition under which the dye stream remained intact not only varied with the
flow rate through the tube, but also with the fluid density, fluid viscosity, and the
diameter of the tube.
Figure 2.6
Experimental
apparatus used to
determine Reynolds
number
28 Modeling Theory Chapter 2
Based on experimental evidence gathered by Reynolds and dimensional analysis, a
dimensionless number can be computed and used to characterize flow regime. Con-
ceptually, the Reynolds number can be thought of as the ratio between inertial and vis-
cous forces in a fluid. The Reynolds number for full flowing circular pipes can be
found using the following equation:
Re
VD

------------
VD
v
-------- = = (2.11)
where Re = Reynolds number
D = pipeline diameter (L)
= fluid density (M/L
3
)
= absolute viscosity (M/L/T)
v = kinematic viscosity (L
2
/T)
The ranges of the Reynolds number that define the three flow regimes are shown in
Table 2.1. The flow of water through municipal water systems is almost always turbu-
lent, except in the periphery where water demand is low and intermittent, and may
result in laminar and stagnant flow conditions.
Table 2.1 Reynolds number for various flow regimes
Flow Regime Reynolds Number
Laminar < 2000
Transitional 20004000
Turbulent > 4000
Velocity Profiles. Due to the shear stresses along the walls of a pipe, the velocity
in a pipeline is not uniform over the pipe diameter. Rather, the fluid velocity is zero at
the pipe wall. Fluid velocity increases with distance from the pipe wall, with the max-
imum occurring along the centerline of the pipe. Figure 2.7 illustrates the variation of
fluid velocity within a pipe, also called the velocity profile.
The shape of the velocity profile will vary depending on whether the flow regime is
laminar or turbulent. In laminar flow, the fluid particles travel in parallel layers or
lamina, producing very strong shear stresses between adjacent layers, and causing the
dye streak in Reynolds experiment to remain intact. Mathematically, the velocity
profile in laminar flow is shaped like a parabola as shown in Figure 2.7. In laminar
flow, the head loss through a pipe segment is primarily a function of the fluid viscos-
ity, not the internal pipe roughness.
Turbulent flow is characterized by eddies that produce random variations in the veloc-
ity profiles. Although the velocity profile of turbulent flow is more erratic than that of
laminar flow, the mean velocity profile actually exhibits less variation across the pipe.
The velocity profiles for both turbulent and laminar flows are shown in Figure 2.7.
Section 2.3 Energy Concepts 29
Figure 2.7
Velocity profiles for
different flow regimes
v
Laminar Profile
v
UniformVelocity Profile
v
Turbulent Profile
2.3 ENERGY CONCEPTS
Fluids possess energy in three forms. The amount of energy depends on the fluids
movement (kinetic energy), elevation (potential energy), and pressure (pressure
energy). In a hydraulic system, a fluid can have all three types of energy associated
with it simultaneously. The total energy associated with a fluid per unit weight of the
fluid is called head. The kinetic energy is called velocity head (V
2
/2g), the potential
energy is called elevation head (Z), and the internal pressure energy is called pressure
head (P/ J ). While typical units for energy are foot-pounds (Joules), the units of total
head are feet (meters).
H Z
P
J
---
V
2
2g
------ + + = (2.12)
where H = total head (L)
Z = elevation above datum (L)
P = pressure (M/L/T
2
)
J = fluid specific weight (M/L
2
/T
2
)
V = velocity (L/T)
g = gravitational acceleration constant (L/T
2
)
Each point in the system has a unique head associated with it. A line plotted of total
head versus distance through a system is called the energy grade line (EGL). The sum
of the elevation head and pressure head yields the hydraulic grade line (HGL), which
corresponds to the height that water will rise vertically in a tube attached to the pipe
and open to the atmosphere. Figure 2.8 shows the EGL and HGL for a simple pipe-
line.
In most water distribution applications, the elevation and pressure head terms are
much greater than the velocity head term. For this reason, velocity head is often
ignored, and modelers work in terms of hydraulic grades rather than energy grades.
Therefore, given a datum elevation and a hydraulic grade line, the pressure can be
determined as
P J HGL Z = (2.13)
30 Modeling Theory Chapter 2
Figure 2.8
Energy and hydraulic
grade lines
Datum
H
G
L
E
G
L
Elevation Head, Z
Pressure Head,
P
V
2
2g
Velocity Head,
Flow
Head Loss, h
L
Energy Losses
Energy losses, also called head losses, are generally the result of two mechanisms:
Friction along the pipe walls
Turbulence due to changes in streamlines through fittings and appurtenances
Head losses along the pipe wall are called friction losses or head losses due to fric-
tion, while losses due to turbulence within the bulk fluid are called minor losses.
2.4 FRICTION LOSSES
When a liquid flows through a pipeline, shear stresses develop between the liquid and
the pipe wall. This shear stress is a result of friction, and its magnitude is dependent
on the properties of the fluid that is passing through the pipe, the speed at which it is
moving, the internal roughness of the pipe, and the length and diameter of the pipe.
Consider, for example, the pipe segment shown in Figure 2.9. A force balance on the
fluid element contained within a pipe section can be used to form a general expression
describing the head loss due to friction. Note the forces in action:
Pressure difference between sections 1 and 2
The weight of the fluid volume contained between sections 1 and 2
The shear at the pipe walls between sections 1 and 2
Assuming the flow in the pipeline has a constant velocity (that is, acceleration is equal
to zero), the system can be balanced based on the pressure difference, gravitational
forces, and shear forces.
P
1
A
1
P
2
A
2
AL o ( ) sin t
o
NL 0 =
Section 2.4 Friction Losses 31
(2.14)
where P
1
= pressure at section 1(M/L/T
2
)
A
1
= cross-sectional area of section 1(L
2
)
P
2
= pressure at section 2 (M/L/T
2
)
A
2
= cross-sectional area of section 2 (L
2
)
A = average area between section 1 and section 2 (L
2
)
L = distance between section 1 and section 2 (L)
= fluid specific weight (M/L
2
/T
2
)
o = angle of the pipe to horizontal
t
o
= shear stress along pipe wall (M/L/T
2
)
N = perimeter of pipeline cross-section (L)
Figure 2.9
Free body diagram of
water flowing in an
inclined pipe
The last term on the left side of Equation 2.14 represents the friction losses along the
pipe wall between the two sections. By recognizing that sin( o ) = (Z
2
-Z
1
)/L, the equa-
tion for head loss due to friction can be rewritten to obtain the following equation.
(Note that the velocity head is not considered in this case because the pipe diameters,
and therefore the velocity heads, are the same.)
h
L
t
o
NL
A
-------
P
1

------ Z
1
+
\ .
| |
P
2

------ Z
2
+
\ .
| |
= = (2.15)
32 Modeling Theory Chapter 2
where h
L
= head loss due to friction (L)
Z
1
= elevation of centroid of section 1 (L)
Z
2
= elevation of centroid of section 2 (L)
Recall that the shear stresses in a fluid can be found analytically for laminar flow
using Newtons law of viscosity. Shear stress is a function of the viscosity and veloc-
ity gradient of the fluid, the fluid specific weight (or density), and the diameter of the
pipeline. The roughness of the pipe wall is also a factor (that is, the rougher the pipe
wall, the larger the shear stress). Combining all these factors, it can be seen that
t
o
F V D c , , , , ( ) = (2.16)
where = fluid density (M/L
3
)
= absolute viscosity (M/L/T)
V = average fluid velocity (L/T)
D = diameter (L)
c = index of internal pipe roughness (L)
Darcy-Weisbach Formula
Using dimensional analysis, the Darcy-Weisbach formula was developed. The for-
mula is an equation for head loss expressed in terms of the variables listed in Equation
2.16, as follows (note that head loss is expressed with units of length):
h
L
f
LV
2
D2g
-----------
8fLQ
2
gD
5
t
2
---------------- = = (2.17)
where f = Darcy-Weisbach friction factor
g = gravitational acceleration constant (L/T
2
)
Q = pipeline flow rate (L
3
/T)
The Darcy-Weisbach friction factor, f, is a function of the same variables as wall shear
stress (Equation 2.16). Again using dimensional analysis, a functional relationship for
the friction factor can be developed:
f F
VD

------------
c
D
---- ,
\ .
| |
F Re
c
D
---- ,
\ .
| |
= = (2.18)
where Re = Reynolds number
The Darcy-Weisbach friction factor is dependent on the velocity, density, and viscos-
ity of the fluid; the size of the pipe in which the fluid is flowing; and the internal
roughness of the pipe. The fluid velocity, density, viscosity, and pipe size are
expressed in terms of the Reynolds number. The internal roughness is expressed in
terms of a variable called the relative roughness, which is the internal pipe roughness
( c ) divided by the pipe diameter (D).
Section 2.4 Friction Losses 33
In the early 1930s, the German researcher Nikuradse performed an experiment that
would become fundamental in head loss determination (Nikuradse, 1932). He glued
uniformly sized sand grains to the insides of three pipes of different sizes. His experi-
ments showed that the curve of f versus Re is smooth for the same values of c /D.
Partly because of Nikuradses sand grain experiments, the quantity c is called the
equivalent sand grain roughness of the pipe. Table 2.2 provides values of c for
various materials.
Other researchers conducted experiments on artificially roughened pipes to generate
data describing pipe friction factors for a wide range of relative roughness values.
Table 2.2 Equivalent sand grain roughness for various pipe materials
Equivalent Sand Roughness, c
Material (ft) (mm)
Copper, brass 1x10
-4
3x10
-3
3.05x10
-2
0.9
Wrought iron, steel 1.5x10
-4
8x10
-3
4.6x10
-2
2.4
Asphalted cast iron 4x10
-4
7x10
-3
0.1 2.1
Galvanized iron 3.3x10
-4
1.5x10
-2
0.102 4.6
Cast iron 8x10
-4
1.8x10
-2
0.2 5.5
Concrete 10
-3
10
-2
0.3 3.0
Uncoated cast iron 7.4x10
-4
0.226
Coated cast iron 3.3x10
-4
0.102
Coated spun iron 1.8x10
-4
5.6x10
-2
Cement 1.3x10
-3
4x10
-3
0.4 1.2
Wrought iron 1.7x10
-4
5x10
-2
Uncoated steel 9.2x10
-5
2.8x10
-2
Coated steel 1.8x10
-4
5.8x10
-2
Wood stave 6x10
-4
3x10
-3
0.2 0.9
PVC 5x10
-6
1.5x10
-3
Compiled from Lamont (1981), Moody (1944), and Mays (1999)
Colebrook-White Equation and the Moody Diagram. Numerous for-
mulas exist that relate the friction factor to the Reynolds number and relative rough-
ness. One of the earliest and most popular of these formulas is the Colebrook-
White equation:
1
f
----- 0.86
c
3.7D
------------
2.51
Re f
------------ +
\ .
| |
ln = (2.19)
The difficulty with using the Colebrook-White equation is that it is an implicit func-
tion of the friction factor (f is found on both sides of the equation). Typically, the
equation is solved by iterating through assumed values of f until both sides are equal.
34 Modeling Theory Chapter 2
The Moody diagram, shown in Figure 2.10, was developed from the Colebrook-White
equation as a graphical solution for the Darcy-Weisbach friction factor.
It is interesting to note that for laminar flow (low Re) the friction factor is a linear
function of the Reynolds number, while in the fully turbulent range (high c /D and
high Re) the friction factor is only a function of the relative roughness. This difference
occurs because the effect of roughness is negligible for laminar flow, while for very
turbulent flow the viscous forces become negligible.
Swamee-Jain Formula. Much easier to solve than the iterative Colebrook-
White formula, the formula developed by Swamee and Jain (1976) also approximates
the Darcy-Weisbach friction factor. This equation is an explicit function of the Rey-
nolds number and the relative roughness, and is accurate to within about one percent
of the Colebrook-White equation over a range of
4 10
3
Re 1 10
8
s s and
1 10
6
c D 1 10
2
s s
f
1.325
c
3.7D
------------
5.74
Re
0.9
------------ +
\ .
| |
ln
2
-------------------------------------------------- = (2.20)
Because of its relative simplicity and reasonable accuracy, most water distribution
system modeling software packages use the Swamee-Jain formula to compute the
friction factor.
Hazen-Williams
Another frequently used head loss expression, particularly in North America, is the
Hazen-Williams formula (Williams and Hazen, 1920; ASCE, 1992):
h
L
C
f
L
C
1.852
D
4.87
---------------------------Q
1.852
= (2.21)
where h
L
= head loss due to friction (ft, m)
L = distance between sections 1 and 2 (ft, m)
C = Hazen-Williams C-factor
D = diameter (ft, m)
Q = pipeline flow rate (cfs, m
3
/s)
C
f
= unit conversion factor (4.73 English, 10.7 SI)
The Hazen-Williams formula uses many of the same variables as Darcy-Weisbach,
but instead of using a friction factor, the Hazen-Williams formula uses a pipe carrying
capacity factor, C. Higher C-factors represent smoother pipes (with higher carrying
capacities) and lower C-factors describe rougher pipes. Table 2.3 shows typical C-
factors for various pipe materials, based on Lamont (1981).
Section 2.4 Friction Losses 35
Figure 2.10
Moody diagram
V
D
V
a
l
u
e
s
o
f
(
)
f
o
r
w
a
t
e
r
a
t
6
0
F
(
d
i
a
m
e
t
e
r
i
n
i
n
c
h
e
s
,
v
e
l
o
c
i
t
y
i
n
f
t
.
/
s
e
c
.
)
0
.
1
0
.
2
0
.
4
0
.
6
1
2
4
6
1
0
2
0
4
0
6
0
1
0
0
2
0
0
4
0
0
6
0
0
1
0
0
0
2
0
0
0
4
0
0
0
1
0
,
0
0
0
0
.
1
0
0
.
0
9
0
.
0
8
0
.
0
7
0
.
0
6
0
.
0
5
0
.
0
4
0
.
0
3
0
.
0
2
5
0
.
0
2
0
0
.
0
1
5
0
.
0
1
0
0
.
0
0
9
0
.
0
0
8
f F r i c t i o n f a c t o r , =
h
L
L
D
V
2
2 g
1
0
3
2
(
1
0
)
3
4
6
8
1
0
4
2
(
1
0
)
4
4
6
8
1
0
5
2
(
1
0
)
5
4
6
8
1
0
6
2
(
1
0
)
6
4
6
8
1
0
7
2
(
1
0
)
7
4
6
8
1
0
8
R
e
R
e
y
n
o
l
d
s
n
u
m
b
e
r
,
=
V
D
v
/ D R e l a t i v e r o u g h n e s s ,
0
.
0
5
0
.
0
3
0
.
0
2
0
.
0
1
0
.
0
0
6
0
.
0
0
4
0
.
0
0
2
0
.
0
0
1
0
.
0
0
0
6
0
.
0
0
0
4
0
.
0
0
0
2
0
.
0
0
0
0
5
f
=
L
a
m
i
n
a
r
f
l
o
w
,
R
e
6
4
D
r
a
w
n
T
u
b
i
n
g
S
t
e
e
l
o
r
w
r
o
u
g
h
t
i
r
o
n
A
s
p
h
a
l
t
e
d
c
a
s
t
i
r
o
n
G
a
l
v
a
n
i
z
e
d
i
r
o
n
C
a
s
t
i
r
o
n
W
o
o
d
s
t
a
v
e
C
o
n
c
r
e
t
e
R
i
v
e
t
e
d
s
t
e
e
l
0
.
0
0
0
0
0
5
0
.
0
0
0
1
5
0
.
0
0
0
4
0
.
0
0
0
5
0
.
0
0
0
8
5
0
.
0
0
0
6
-
0
.
0
0
3
0
.
0
0
1
-
0
.
0
1
0
.
0
0
3
-
0
.
0
3
f
t
.
m
m
0
.
0
0
1
5
0
.
0
4
5
0
.
1
2
0
.
1
5
0
.
2
5
0
.
1
8
-
0
.
9
0
.
3
-
3
0
.
9
-
9
S
m
o
o
t
h
p
i
p
e
s
F
r
o
m

L
.

F
.

M
o
o
d
y
,

F
r
i
c
t
i
o
n

F
a
c
t
o
r
s

f
o
r

P
i
p
e

F
l
o
w
,

T
r
a
n
s
.

A
.
S
.
M
.
E
.
,

V
o
l
.

6
6
,

1
9
4
4
,

u
s
e
d

w
i
t
h

p
e
r
m
i
s
s
i
o
n
.

36 Modeling Theory Chapter 2
Lamont found that it was not possible to develop a single correlation between pipe
age and C-factor and that, instead, the decrease in C-factor also depended heavily on
the corrosiveness of the water being carried. He developed four separate trends in
carrying capacity loss depending on the attack of the water on the pipe. Trend 1,
slight attack, corresponded to water that was only mildly corrosive. Trend 4, severe
attack, corresponded to water that would rapidly attack cast iron pipe. As can be seen
from Table 2.3, the extent of attack can significantly affect the C-factor. Testing pipes
to determine the loss of carrying capacity is discussed further on page 196.
Table 2.3 C-factors for various pipe materials
C-factor Values for Discrete Pipe Diameters
Type of Pipe
1.0 in.
(2.5 cm)
3.0 in.
(7.6 cm)
6.0 in.
(15.2 cm)
12 in.
(30 cm)
24 in.
(61 cm)
48 in.
(122 cm)
Uncoated cast iron - smooth and
new
121 125 130 132 134
Coated cast iron - smooth and
new
129 133 138 140 141
30 years old
Trend 1 - slight attack 100 106 112 117 120
Trend 2 - moderate attack 83 90 97 102 107
Trend 3 - appreciable
attack
59 70 78 83 89
Trend 4 - severe attack 41 50 58 66 73
60 years old
Trend 1 - slight attack 90 97 102 107 112
Trend 2 - moderate attack 69 79 85 92 96
Trend 3 - appreciable
attack
49 58 66 72 78
Trend 4 - severe attack 30 39 48 56 62
100 years old
Trend 1 - slight attack 81 89 95 100 104
Trend 2 - moderate attack 61 70 78 83 89
Trend 3 - appreciable
attack
40 49 57 64 71
Trend 4 - severe attack 21 30 39 46 54
Miscellaneous
Newly scraped mains 109 116 121 125 127
Newly brushed mains 97 104 108 112 115
Coated spun iron - smooth and
new
137 142 145 148 148
Old - take as coated cast iron
of same age
Galvanized iron - smooth and
new
120 129 133
Wrought iron - smooth and new 129 137 142
Coated steel - smooth and new 129 137 142 145 148 148
Uncoated steel - smooth and
new
134 142 145 147 150 150
Section 2.4 Friction Losses 37
From a purely theoretical standpoint, the C-factor of a pipe should vary with the flow
velocity under turbulent conditions. Equation 2.22 can be used to adjust the C-factor
for different velocities, but the effects of this correction are usually minimal. A two-
fold increase in the flow velocity correlates to an apparent five percent decrease in the
roughness factor. This difference is usually within the error range for the roughness
estimate in the first place, so most engineers assume the C-factor remains constant
regardless of flow (Walski, 1984). However, if C-factor tests are done at very high
velocities (i.e., >10 ft/s), then a significant error can result when the resulting C-
factors are used to predict head loss at low velocities.
C C
o
V
o
V
------
\ .
| |
0.081
= (2.22)
where C = velocity adjusted C-factor
C
o
= reference C-factor
V
o
= reference value of velocity at which C
0
was determined (L/T)
Manning Equation
Another head loss expression more typically associated with open channel flow is the
Manning equation:
Coated asbestos cement - clean 147 149 150 152
Uncoated asbestos cement -
clean
142 145 147 150
Spun cement-lined and spun
bitumen- lined - clean
147 149 150 152 153
Smooth pipe (including lead,
brass, copper, polyethylene,
and PVC) - clean
140 147 149 150 152 153
PVC wavy - clean 134 142 145 147 150 150
Concrete - Scobey
Class 1 - Cs = 0.27; clean 69 79 84 90 95
Class 2 - Cs = 0.31; clean 95 102 106 110 113
Class 3 - Cs = 0.345; clean 109 116 121 125 127
Class 4 - Cs = 0.37; clean 121 125 130 132 134
Best - Cs = 0.40; clean 129 133 138 140 141
Tate relined pipes - clean 109 116 121 125 127
Prestressed concrete pipes -
clean
147 150 150
Lamont (1981)
Table 2.3 (cont.) C-factors for various pipe materials
C-factor Values for Discrete Pipe Diameters
Type of Pipe
1.0 in.
(2.5 cm)
3.0 in.
(7.6 cm)
6.0 in.
(15.2 cm)
12 in.
(30 cm)
24 in.
(61 cm)
48 in.
(122 cm)
h
L
C
f
L nQ ( )
2
D
5.33
------------------------ =
38 Modeling Theory Chapter 2
(2.23)
where n = Manning roughness coefficient
C
f
= unit conversion factor (4.66 English, 10.29 SI)
As with the previous head loss expressions, the head loss computed using Manning
equation is dependent on the pipe length and diameter, the discharge or flow through
the pipe, and a roughness coefficient. In this case, a higher value of n represents a
higher internal pipe roughness. Table 2.4 provides typical Mannings roughness coef-
ficients for commonly used pipe materials.
Table 2.4 Mannings roughness values
Material
Manning
Coefficient
Material
Manning
Coefficient
Asbestos cement .011 Corrugated metal .022
Brass .011 Galvanized iron .016
Brick .015 Lead .011
Cast iron, new .012 Plastic .009
Concrete Steel
Steel forms .011 Coal-tar enamel .010
Wooden forms .015 New unlined .011
Centrifugally spun .013 Riveted .019
Copper .011 Wood stave .012
Comparison of Friction Loss Methods
Most hydraulic models have features that allow the user to select from the Darcy-
Weisbach, Hazen-Williams, or Manning head loss formulas, depending on the nature
of the problem and the users preferences.
The Darcy-Weisbach formula is a more physically-based equation, derived from the
basic governing equations of Newtons Second Law. With appropriate fluid viscosi-
ties and densities, Darcy-Weisbach can be used to find the head loss in a pipe for any
Newtonian fluid in any flow regime.
The Hazen-Williams and Manning formulas, on the other hand, are empirically-based
expressions (meaning that they were developed from experimental data), and gener-
ally only apply to water under turbulent flow conditions.
The Hazen-Williams formula is the predominant equation used in the United States,
and Darcy-Weisbach is predominant in Europe. The Manning formula is not typically
used for water distribution modeling; however, it is sometimes used in Australia.
Table 2.5 presents these three equations in several common unit configurations. These
equations solve for the friction slope (S
f
), which is the head loss per unit length of
pipe.
Table 2.5 Friction loss equations in typical units
Equation Q (m
3
/s); D (m) Q (cfs); D (ft) Q (gpm); D (in.)
Darcy-Weisbach
S
f
0.083f Q
2
D
5
---------------------- = S
f
0.025f Q
2
D
5
---------------------- = S
f
0.031f Q
2
D
5
---------------------- =
Hazen-Williams
S
f
10.7
D
4.87
------------
Q
C
----
\ .
| |
1.852
= S
f
4.73
D
4.87
------------
Q
C
----
\ .
| |
1.852
= S
f
10.5
D
4.87
------------
Q
C
----
\ .
| |
1.852
=
Manning
S
f
10.3 nQ ( )
2
D
5.33
------------------------- = S
f
4.66 nQ ( )
2
D
5.33
------------------------- = S
f
13.2 nQ ( )
2
D
5.33
------------------------- =
Compiled from ASCE (1975) and ASCE/WEF (1982)
Section 2.5 Minor Losses 39
2.5 MINOR LOSSES
Head losses also occur at valves, tees, bends (see Figure 2.11), reducers, and other
appurtenances within the piping system. These losses, called minor losses, are due to
turbulence within the bulk flow as it moves through fittings and bends. Figure 2.12
illustrates the turbulent eddies that develop within the bulk flow as it travels through a
valve and a 90-degree bend.
Figure 2.11
48-in. elbow fitting
Head loss due to minor losses can be computed by multiplying a minor loss coeffi-
cient by the velocity head, as shown in Equation 2.24.
40 Modeling Theory Chapter 2
Figure 2.12
Valve and bend cross-
sections generating
minor losses
h
m
K
L
V
2
2g
------ K
L
Q
2
2gA
2
------------ = = (2.24)
where h
m
= head loss due to minor losses (L)
K
L
= minor loss coefficient
V = velocity (L/T)
g = gravitational acceleration constant (L/T
2
)
A = cross-sectional area (L
2
)
Q = flow rate (L
3
/T)
Minor loss coefficients are found experimentally, and data are available for many dif-
ferent types of fittings and appurtenances. Table 2.6 provides a list of minor loss coef-
ficients associated with several of the most commonly used fittings. More thorough
treatments of minor loss coefficients can be found in Crane (1972), Miller (1978), and
Idelchik (1999).
For water distribution systems, minor losses are generally much smaller than the head
losses due to friction (hence the term "minor" loss). For this reason, many modelers
frequently choose to neglect minor losses. In some cases, however, such as at pump
stations or valve manifolds where there may be more fittings and higher velocities,
minor losses can play a significant role in the piping system under consideration.
Like pipe roughness coefficients, minor head loss coefficients will vary somewhat
with velocity. For most practical network problems, however, the minor loss coeffi-
cient is treated as constant.
Valve Coefficient
Most valve manufacturers can provide a chart of percent opening versus valve coeffi-
cent (C
v
), which can be related to the minor loss (K
L
) by using Equation 2.25.
Table 2.6 Minor loss coefficients
Fitting K
L
Fitting K
L
Pipe entrance 90
o
smooth bend
Bellmouth 0.03-0.05 Bend radius/D = 4 0.16-0.18
Rounded 0.12-0.25 Bend radius/D = 2 0.19-0.25
Sharp-edged 0.50 Bend radius/D = 1 0.35-0.40
Projecting 0.78 Mitered bend
Contraction sudden
u = 15
o 0.05
D
2
/D
1
=0.80 0.18
u = 30
o 0.10
D
2
/D
1
=0.50 0.37
u = 45
o 0.20
D
2
/D
1
=0.20 0.49
u = 60
o 0.35
Contraction conical
u = 90
o 0.80
D
2
/D
1
=0.80 0.05 Tee
D
2
/D
1
=0.50 0.07 Line flow 0.30-0.40
D
2
/D
1
=0.20 0.08 Branch flow 0.75-1.80
Expansion sudden Tapping T Branch
D
2
/D
1
=0.80 0.16 d = tapping hole diameter
D = main line diameter
1.97/(d/D)
4
D
2
/D
1
=0.50 0.57 Cross
D
2
/D
1
=0.20 0.92 Line flow 0.50
Expansion conical Branch flow 0.75
D
2
/D
1
=0.80 0.03 45
o
Wye
D
2
/D
1
=0.50 0.08 Line flow 0.30
D
2
/D
1
=0.20 0.13 Branch flow 0.50
Gate valve open 0.39 Check valve conventional 4.0
3/4 open 1.10 Check valve clearway 1.5
1/2 open 4.8 Check valve ball 4.5
1/4 open 27 Cock straight through 0.5
Globe valve open 10 Foot valve hinged 2.2
Angle valve open 4.3 Foot valve poppet 12.5
Butterfly valve open 1.2
Walski (1984)
K
L
C
f
D
4
C
v
2
=
Section 2.5 Minor Losses 41
(2.25)
where D = diameter (in., m)
C
v
= valve coefficient [gpm/(psi)
0.5
, (m
3
/s)/(kPa)
0.5
]
C
f
= unit conversion factor (880 English, 1.22 SI)
42 Modeling Theory Chapter 2
Equivalent Pipe Length
Rather than including minor loss coefficients directly, a modeler may choose to adjust
the modeled pipe length to account for minor losses by adding an equivalent length of
pipe for each minor loss. Given the minor loss coefficient for a valve or fitting, the
equivalent length of pipe to give the same head loss can be calculated as:
L
e
K
L
D
f
----------- = (2.26)
where L
e
= equivalent length of pipe (L)
D = diameter of equivalent pipe (L)
f = Darcy-Weisbach friction factor
The practice of assigning equivalent pipe lengths was typically used when hand calcu-
lations were more common because it could save time in the overall analysis of a
pipeline. With modern computer modeling techniques, this is no longer a widespread
practice. Because it is now so easy to use minor loss coefficients directly within a
hydraulic model, the process of determining equivalent lengths is actually less effi-
cient. In addition, use of equivalent pipe lengths can unfavorably affect the travel time
predictions that are important in many water quality calculations.
2.6 RESISTANCE COEFFICIENTS
Many related expressions for head loss have been developed. They can be mathemati-
cally generalized with the introduction of a variable referred to as a resistance coeffi-
cient. This format allows the equation to remain essentially the same regardless of
which friction method is used, making it ideal for hydraulic modeling.
h
L
K
P
Q
z
= (2.27)
where h
L
= head loss due to friction (L)
K
P
= pipe resistance coefficient (T
z
/L
3z - 1
)
Q = pipeline flow rate (L
3
/T)
z = exponent on flow term
Equations for computing K
P
with the various head loss methods are given next.
Darcy-Weisbach
K
P
f
L
2gA
z
D
----------------- = (2.28)
where f = Darcy-Weisbach friction factor
L = length of pipe (L)
D = pipe diameter (L)
A = cross-sectional area of pipeline (L
2
)
z = 2
Section 2.6 Resistance Coefficients 43
Hazen-Williams
K
P
C
f
L
C
z
D
4.87
------------------- = (2.29)
where K
P
= pipe resistance coefficient (s
z
/ft
3z-1
, s
z
/m
3z-1
)
L = length of pipe (ft, m)
C = C-factor with velocity adjustment
z = 1.852
D = pipe diameter (ft, m)
C
f
= unit conversion factor (4.73 English, 10.7 SI)
Manning
K
P
C
f
Ln
z
D
5.33
-------------- = (2.30)
where n = Mannings roughness coefficient
z = 2
C
f
= unit conversion factor [4.64 English, 10.3 SI (ASCE/WEF, 1982)]
Minor Losses
A resistance coefficient can also be defined for minor losses, as shown in the equation
below. Like the pipe resistance coefficient, the resistance coefficient for minor losses
is a function of the physical characteristics of the fitting or appurtenance and
the discharge.
h
m
K
M
Q
2
= (2.31)
where h
m
= head loss due to minor losses (L)
K
M
= minor loss resistance coefficient (T
2
/L
5
)
Q = pipeline flow rate (L
3
/T)
Solving for the minor loss resistance coefficient by substituting Equation 2.24
results in:
K
M
K
L _
2gA
2
-------------- = (2.32)
where K
L _
= sum of individual minor loss coefficients
44 Modeling Theory Chapter 2
2.7 ENERGY GAINS PUMPS
On many occasions, energy needs to be added to a hydraulic system to overcome ele-
vation differences, friction losses, and minor losses. A pump is a device to which
mechanical energy is applied and transferred to the water as total head. The head
added is called pump head and is a function of the flow rate through the pump. The
following discussion is oriented toward centrifugal pumps because they are the most
frequently used pumps in water distribution systems. Additional information about
pumps can be found in Bosserman (2000), Hydraulic Institute Standards (2000),
Karassik (1976), and Sanks (1998).
Pump Head-Discharge Relationship
The relationship between pump head and pump discharge is given in the form of a
head versus discharge curve (also called a head characteristic curve) similar to the
one shown in Figure 2.13. This curve defines the relationship between the head that
the pump adds and the amount of flow that the pump passes. The pump head versus
discharge relationship is nonlinear, and as one would expect, the more water the pump
passes, the less head it can add. The head that is plotted in the head characteristic
curve is the head difference across the pump, called the total dynamic head (TDH).
This curve must be described as a mathematical function to be used in a hydraulic
simulation. Some models fit a polynomial curve to selected data points, but a more
common approach is to describe the curve by using a power function in the following
form:
h
P
h
o
c Q
P
m
= (2.33)
where h
P
= pump head (L)
h
o
= cutoff (shutoff) head (pump head at zero flow) (L)
Q
P
= pump discharge (L
3
/T)
c, m = coefficients describing pump curve shape
More information on pump performance testing is available in Chapter 5 (see page
197).
Affinity Laws for Variable-Speed Pumps. A centrifugal pumps charac-
teristic curve is fixed for a given motor speed and impeller diameter, but it can be
determined for any speed and any diameter by applying relationships called the affin-
ity laws. For variable-speed pumps, these affinity laws are presented as follows:
Q
P1
Q
P2
n
1
n
2
= (2.34)
h
P1
h
P2
n
1
n
2
( )
2
= (2.35)
where Q
P1
= pump flow at speed 1 (L
3
/T)
n
1
= pump speed 1 (1/T)
h
P1
= pump head at speed 1 (L)
Section 2.7 Energy Gains Pumps 45
Figure 2.13
Pump head
characteristic curve
0.0
50.0
100.0
150.0
200.0
0.0 250 500 750 1000 1250 1500 1750
Shutoff Head
Flow, gpm
H
e
a
d
,
f
t
Design Point
Maximum Flow
Thus, pump discharge rate is directly proportional to pump speed, and pump dis-
charge head is proportional to the square of the speed. Using this relationship, once
the pump curve at any one speed is known, then the curve at another speed can be pre-
dicted. Figure 2.14 illustrates the affinity laws for variable-speed pumps where the
line through the pump head characteristic curves represents the locus of best effi-
ciency points.
Inserting Equations 2.34 and 2.35 into Equation 2.33 and solving for h gives a general
equation for adjusting pump head curves for speed:
h
P2
n
2
h
o
cn
2 m
Q
P2
m
= (2.36)
where n = n
2
/n
1
System Head Curves
The purpose of a pump is to overcome elevation differences and head losses due to
pipe friction and fittings. The amount of head the pump must add to overcome eleva-
tion differences is dependent on system characteristics and topology (and independent
of the pump discharge rate), and is referred to as static head or static lift. Friction and
minor losses, however, are highly dependent on the rate of discharge through the
pump. When these losses are added to the static head for a series of discharge rates,
the resulting plot is called a system head curve (see Figure 2.15).
46 Modeling Theory Chapter 2
Figure 2.14
Relative speed factors
for variable-speed
pumps
Figure 2.15
A family of system
head curves
Section 2.7 Energy Gains Pumps 47
The pump characteristic curve is a function of the pump and independent of the sys-
tem, while the system head curve is dependent on the system and is independent of
the pump. Unlike the pump curve, which is fixed for a given pump at a given speed,
the system head curve is continually sliding up and down as tank water levels change
and demands change. Rather than there being a unique system head curve, a family of
system head curves forms a band on the graph.
For the case of a single pipeline between two points, the system head curve can be
described in equation form as follows:
H h
l
K
P
Q
z
K
M
Q
2

+ = (2.37)
where H = total head (L)
h
l
= static lift (L)
K
P
= pipe resistance coefficient (T
z
/L
3z - 1
)
Q = pipe discharge (L
3
/T)
z = coefficient
K
M
= minor loss resistance coefficient (T
2
/L
5
)
Thus, the head losses and minor losses associated with each segment of pipe are
summed along the total length of the pipeline. When the system is more complex, the
interdependencies of the hydraulic network make it impossible to write a single equa-
tion to describe a point on the system curve. In these cases, hydraulic analysis using a
hydraulic model may be needed. It is helpful to visualize the hydraulic grade line as
increasing abruptly at a pump and sloping downward as the water flows through pipes
and valves (see Figure 2.16).
Figure 2.16
Schematic of
hydraulic grade line
for a pumped system
48 Modeling Theory Chapter 2
Pump Operating Point
When the pump head discharge curve and the system head curve are plotted on the
same axes (as shown in Figure 2.17), only one point lies on both the pump character-
istic curve and the system head curve. This intersection defines the pump operating
point, which represents the discharge that will pass through the pump and the head
that the pump will add. This head is equal to the head needed to overcome the static
head and other losses in the system.
Figure 2.17
System operating
point
Pump Head Curve
Pump Operating Point
Head Losses
Flow, gpm
Static Lift
H
e
a
d
,
f
t
S
ystem
H
ead
C
u
rve
Other Uses of Pump Curves
In addition to the pump head-discharge curve, other curves representing pump behav-
ior describe power, water horsepower, and efficiency (see Figure 2.18), and are dis-
cussed further in Chapter 3 (see page 95) and Chapter 5 (see page 197). Since utilities
want to minimize the amount of energy necessary for system operation, the engineer
should select pumps that run as efficiently as possible. Pump operating costs are dis-
cussed further in Chapter 10 (see page 436).
Another issue when designing a pump is the net positive suction head (NPSH)
required (see page 324). NPSH is the head that is present at the suction side of the
pump. Each pump requires that the available NPSH exceed the required NPSH to
ensure that local pressures within the pump do not drop below the vapor pressure of
the fluid, causing cavitation. As discussed on page 23, cavitation is essentially a boil-
ing of the liquid within the pump, and it can cause tremendous damage. The NPSH
required is unique for each pump model, and is a function of flow rate. The use of a
calibrated hydraulic model in determining available net positive suction head is dis-
cussed further on page 324.
Section 2.8 Network Hydraulics 49
Figure 2.18
Pump efficiency curve
2.8 NETWORK HYDRAULICS
In networks of interconnected hydraulic elements, every element is influenced by
each of its neighbors; the entire system is interrelated in such a way that the condition
of one element must be consistent with the condition of all other elements. Two con-
cepts define these interconnections:
Conservation of mass
Conservation of energy
Conservation of Mass
The principle of conservation of mass (shown in Figure 2.19) dictates that the fluid
mass entering any pipe will be equal to the mass leaving the pipe (since fluid is typi-
cally neither created nor destroyed in hydraulic systems). In network modeling, all
outflows are lumped at the nodes or junctions.
Q
i
U
pipes
_
0 = (2.38)
where Q
i
= inflow to node in i-th pipe (L
3
/T)
U = water used at node (L
3
/T)
50 Modeling Theory Chapter 2
Figure 2.19
Conservation of mass
principle
Q
1
Q
2
Q
3
U
Note that for pipe outflows from the node, the value of Q is negative.
When extended-period simulations are considered, water can be stored and
withdrawn from tanks, thus a term is needed to describe the accumulation of water at
certain nodes:
Q
i
U
dS
dt
------
pipes
_
0 = (2.39)
where
dS
dt
------ = change in storage (L
3
/T)
The conservation of mass equation is applied to all junction nodes and tanks in a net-
work, and one equation is written for each of them.
Conservation of Energy
The principle of conservation of energy dictates that the difference in energy between
two points must be the same regardless of the path that is taken (Bernoulli, 1738). For
convenience within a hydraulic analysis, the equation is written in terms of head as
follows:
1
P
1

------
V
1
2
2g
------ h
P _
+ + + Z
2
P
2

------
V
2
2
2g
------ h
L
h
m _
+
_
+ + + = (2.40)
where Z = elevation (L)
P = pressure (M/L/T
2
)
= fluid specific weight (M/L
2
/T
2
)
V = velocity (L/T)
Section 2.8 Network Hydraulics 51
g = gravitational acceleration constant (L/T
2
)
h
P
= head added at pumps (L)
h
L
= head loss in pipes (L)
h
m
= head loss due to minor losses (L)
Thus the difference in energy at any two points connected in a network is equal to the
energy gains from pumps and energy losses in pipes and fittings that occur in the path
between them. This equation can be written for any open path between any two
points. Of particular interest are paths between reservoirs or tanks (where the differ-
ence in head is known), or paths around loops because the changes in energy must
sum to zero, as illustrated in Figure 2.20.
Figure 2.20
The sum of head
losses around a pipe
loop is equal to zero
2 Loss
1

L
o
s
s
3

L
o
s
s
A B
C
A to B to C to A = 0
+ 2 + 1 - 3 = 0
Solving Network Problems
Real water distribution systems do not consist of a single pipe and cannot be
described by a single set of continuity and energy equations. Instead, one continuity
equation must be developed for each node in the system, and one energy equation
must be developed for each pipe (or loop), depending on the method used. For real
systems, these equations can number in the thousands.
The first systematic approach for solving these equations was developed by Hardy
Cross (1936). The invention of digital computers, however, allowed more powerful
numerical techniques to be developed. These techniques set up and solve the system
of equations describing the hydraulics of the network in matrix form. Because the
energy equations are nonlinear in terms of flow and head, they cannot be solved
directly. Instead, these techniques estimate a solution and then iteratively improve it
until the difference between solutions falls within a specified tolerance. At this point,
the hydraulic equations are considered solved.
Some of the methods used in network analysis are described in Bhave (1991); Lansey
and Mays (2000); Larock, Jeppson, and Watters (1999); and Todini and Pilati (1987).
52 Modeling Theory Chapter 2
2.9 WATER QUALITY MODELING
Water quality modeling is a direct extension of hydraulic network modeling and can
be used to perform many useful analyses. Developers of hydraulic network simulation
models recognized the potential for water quality analysis and began adding water
quality calculation features to their models in the mid 1980s. Transport, mixing, and
decay are the fundamental physical and chemical processes typically represented in
water quality models. Water quality simulations also use the network hydraulic solu-
tion as part of their computations. Flow rates in pipes and the flow paths that define
how water travels through the network are used to determine mixing, residence times,
and other hydraulic characteristics affecting disinfectant transport and decay. The
results of an extended period hydraulic simulation can be used as a starting point in
performing a water quality analysis.
The equations describing transport through pipes, mixing at nodes, chemical forma-
tion and decay reactions, and storage and mixing in tanks are adapted from Grayman,
Rossman, and Geldreich (2000). Additional information on water quality models can
be found in Clark and Grayman (1998).
Transport in Pipes
Most water quality models make use of one-dimensional advective-reactive transport
to predict the changes in constituent concentrations due to transport through a pipe,
and to account for formation and decay reactions. Equation 2.41 shows concentration
within a pipe i as a function of distance along its length (x) and time (t).
cC
i
ct
--------
Q
i
A
i
-----
cC
i
cx
-------- u C
i
( ) i , + 1...P = = (2.41)
where C
i
= concentration in pipe i (M/L
3
)
Q
i
= flow rate in pipe i (L
3
/T)
A
i
= cross-sectional area of pipe i (L
2
)
u C
i
( ) = reaction term (M/L
3
/T)
Equation 2.41 must be combined with two boundary condition equations (concentra-
tion at x = 0 and t = 0) to obtain a solution. Solution methods are described later in this
section.
The equation for advective transport is a function of the flow rate in the pipe divided
by the cross-sectional area, which is equal to the mean velocity of the fluid. Thus, the
bulk fluid is transported down the length of the pipe with a velocity that is directly
proportional to the average flow rate. The equation is based on the assumption that
longitudinal dispersion in pipes is negligible and that the bulk fluid is completely
mixed (a valid assumption under turbulent conditions). Furthermore, the equation can
also account for the formation or decay of a substance during transport with the sub-
stitution of a suitable equation into the reaction term. Such an equation will be devel-
oped later. First, however, the nodal mixing equation is presented.
Section 2.9 Water Quality Modeling 53
Mixing at Nodes
Water quality simulation uses a nodal mixing equation to combine concentrations
from individual pipes described by the advective transport equation, and to define the
boundary conditions for each pipe as mentioned previously. The equation is written
by performing a mass balance on concentrations entering a junction node.
C
OUT
j
Q
i
C
i n
i
,
U
j
+
i IN
j
e
_
Q
i
i OUT
j
e
_
------------------------------------------- =
\ .
|
|
|
|
|
| |
(2.42)
where C
OUT
j
= concentration leaving the junction node j (M/L
3
)
OUT
j
= set of pipes leaving node j
IN
j
= set of pipes entering node j
Q
i
= flow rate entering the junction node from pipe i (L
3
/T)
C
i n
i
,
= concentration entering junction node from pipe i (M/L
3
)
U
j
= concentration source at junction node j (M/T)
The nodal mixing equation describes the concentration leaving a network node (either
by advective transport into an adjoining pipe or by removal from the network as a
demand) as a function of the concentrations that enter it. The equation describes the
flow-weighted average of the incoming concentrations. If a source is located at a junc-
tion, constituent mass can also be added and combined in the mixing equation with
the incoming concentrations. Figure 2.21 illustrates how the nodal mixing equation is
used at a pipe junction. Concentrations enter the node with pipe flows. The incoming
concentrations are mixed according to Equation 2.42, and the resulting concentration
is transported through the outgoing pipes modeled as demand leaving the system. The
nodal mixing equation assumes that incoming flows are completely and instanta-
neously mixed. The basis for the assumption is that turbulence occurs at the junction
node, which is usually sufficient for good mixing.
Mixing in Tanks
Pipes are sometimes connected to reservoirs and tanks as opposed to junction nodes.
Again, a mass balance of concentrations entering or leaving the tank or reservoir can
be performed.
dC
k
dt
---------
Q
i
V
k
----- C
i np ,
t ( ) C
k
( ) u C
k
( ) + = (2.43)
where C
k
= concentration within tank or reservoir k (M/L
3
)
Q
i
= flow entering the tank or reservoir from pipe i (L
3
/T)
V
k
= volume in tank or reservoir k (L
3
)
u (C
k
) = reaction term (M/L
3
/T)
54 Modeling Theory Chapter 2
Figure 2.21
Nodal mixing
Demand
Equation 2.43 applies when a tank is filling. During a hydraulic time step in which the
tank is filling, the water entering from upstream pipes mixes with water that is already
in storage. If the concentrations are different, blending occurs. The tank mixing equa-
tion accounts for blending and any reactions that occur within the tank volume during
the hydraulic step. During a hydraulic step in which draining occurs, terms can be
dropped and the equation simplified.
dC
k
dt
--------- u C
k
( ) = (2.44)
Specifically, the dilution term can be dropped because it does not occur. Thus, the
concentration within the volume is subject only to chemical reactions. Furthermore,
the concentration draining from the tank becomes a boundary condition for the advec-
tive transport equation written for the pipe connected to it.
Equations 2.43 and 2.44 assume that concentrations within the tank or reservoir are
completely and instantaneously mixed. This assumption is frequently applied in water
quality models. There are, however, other useful mixing models for simulating flow
processes in tanks and reservoirs (Grayman et al., 1996). For example, contact basins
or clearwells designed to provide sufficient contact time for disinfectants are fre-
quently represented as simple plug-flow reactors using a first in first out (FIFO)
model. In a FIFO model, the first volume of water to enter the tank as inflow is the
first to leave as outflow.
If severe short-circuiting is occurring within the tank, a last in first out (LIFO)
model may be applied, in which the first volume entering the tank during filling is the
Section 2.9 Water Quality Modeling 55
last to leave while draining. More complex tank mixing behavior can be captured
using more generalized compartment models. Compartment models have the ability
to represent mixing processes and time delays within tanks more accurately.
Many water distribution models offer a simple two-compartment model, as shown in
Figure 2.22 (Rossman, 2000). In this type of model, water enters or exits the tank
through a completely mixed inlet-outlet compartment, and if the first compartment is
completely full, the overflow is exchanged with a completely mixed second main
compartment. The inlet-outlet compartment can represent short-circuiting with the
last flow in becoming the first out (LIFO). The main compartment can represent a
stagnant or dead zone that will contain older water than the first compartment. The
only parameter for this model is the fraction of the total tank volume in the first
compartment. Selection of an appropriate value for this fraction is generally done by
comparing model results to field measurements of a tracer or chlorine residual.
Figure 2.22
Two-compartment
mixing model
Inlet -
Outlet
Main Zone
Figure 2.23 illustrates a more complex, three-compartment model for a tank with a
single pipe for filling and draining. This example illustrates a tank that is stratified.
New (good quality) water entering the tank occupies the first compartment and is then
transferred to a mixing compartment containing older water, and finally, to a third
dead-zone compartment that contains much older, poorer quality water. The model
simulates the exchange of water between different compartments, and in doing so,
mimics complex tank mixing dynamics. CompTank, a model that can be used to sim-
ulate the three-compartment model, as well as the other models described, is available
as part of an AWWA Research Foundation report (Grayman et al., 2000).
All the models mentioned in this section can be used to simulate a non-reactive (con-
servative) constituent, as well as decay or formation reactions for substances that react
over time. The models can also be used to represent tanks that either operate in fill
and draw mode or operate with simultaneous inflow and outflow.
Chemical Reaction Terms
Equations 2.42, 2.43, and 2.44 compose the linked system of first-order differential
equations solved by typical water quality simulation algorithms. This set of equations
and the algorithms for solving them can be used to model different chemical reactions
56 Modeling Theory Chapter 2
known to impact water quality in distribution systems. Chemical reaction terms are
present in Equations 2.43 and 2.44. Concentrations within pipes, storage tanks, and
reservoirs are a function of these reaction terms. After water leaves the treatment plant
and enters the distribution system, it is subject to many complex physical and chemi-
cal processes, some of which are poorly understood, and most of which are not mod-
eled. Three chemical processes that are frequently modeled, however, are bulk fluid
reactions, reactions that occur on a surface (typically the pipe wall), and formation
reactions involving a limiting reactant. First, an expression for bulk fluid reactions is
presented, and then a reaction expression that incorporates both bulk and pipe wall
reactions is developed.
Figure 2.23
Three-compartment
tank mixing model
Fill Drain
Good Quality
Good Quality
Mixing Zone
Mixed
Poor Quality
Poor Quality
Air
Bulk Reactions. Bulk fluid reactions occur within the fluid volume and are a
function of constituent concentrations, reaction rate and order, and concentrations of
the formation products. A generalized expression for n
th
order bulk fluid reactions is
developed in Equation 2.45 (Rossman, 2000).
u C ( ) kC
n
= (2.45)
where u C ( ) = reaction term (M/L
3
/T)
k = reaction rate coefficient [(L
3
/M)
n-1
/T]
C = concentration (M/L
3
)
n = reaction rate order constant
Equation 2.45 is the generalized bulk reaction term most frequently used in water
quality simulation models. The rate expression accounts for only a single reactant
Section 2.9 Water Quality Modeling 57
concentration, tacitly assuming that any other reactants (if they participate in the reac-
tion) are available in excess of the concentration necessary to sustain the reaction. The
sign of the reaction rate coefficient, k, signifies that a formation reaction (positive) or
a decay reaction (negative) is occurring. The units of the reaction rate coefficient
depend on the order of the reaction. The order of the reaction depends on the compo-
sition of the reactants and products that are involved in the reaction. The reaction rate
order is frequently determined experimentally.
Zero-, first-, and second-order decay reactions are commonly used to model chemical
processes that occur in distribution systems. Figure 2.24 is a conceptual illustration
showing the change in concentration versus time for these three most common reac-
tion rate orders. Using the generalized expression in Equation 2.45, these reactions
can be modeled by allowing n to equal 0, 1, or 2 and then performing a regression
analysis to experimentally determine the rate coefficient.
The most commonly used reaction model is the first order decay model. This has been
applied to chlorine decay, radon decay, and other decay processes. A first order decay
is equivalent to an exponential decay, represented by Equation 2.46.
C
t
C
o
e
kt
= (2.46)
where C
t
= concentration at time t (M/L
3
)
C
o
= initial concentration (at time zero)
k = reaction rate (1/T)
58 Modeling Theory Chapter 2
For first order reactions, the units of k are (1/T) with values generally expressed in 1/
days or 1/hours. Another way of expressing the speed of the reaction is the concept of
half-life that is frequently used when describing the decay rate for radioactive materi-
als. The half-life is the time it takes for the concentration of a substance to decrease to
50 percent of its original concentration. For example, the half-life of radon is approx-
imately 3.8 days, and the half-life of chlorine can vary from hours to many days. The
relationship between the decay rate, k, and half-life is easily calculated by solving
Equation 2.45 for the time t when C
t
/C
o
is equal to a value of 0.5. This results in
Equation 2.47.
T
0.693
k
------------- = (2.47)
For example, if the decay rate k is 1.0, the half-life is 0.693 days.
Figure 2.24
Conceptual
illustration of
concentration versus
time for zero, first-,
and second-order
decay reactions
2nd Order
Conservative
Time
C
o
n
c
e
n
t
r
a
t
i
o
n
1st Order
0 Order
Bulk and Wall Reactions. Disinfectants are the most frequently modeled con-
stituents in water distribution systems. Upon leaving the plant and entering the distri-
bution system, disinfectants are subject to a poorly characterized set of potential
chemical reactions. Figure 2.25 illustrates the flow of water through a pipe and the
types of chemical reactions with disinfectants that can occur along its length. Chlo-
rine (the most common disinfectant) is shown reacting in the bulk fluid with natural
organic matter (NOM), and at the pipe wall, where oxidation reactions with biofilms
and the pipe material (a cause of corrosion) can occur.
Many disinfectant decay models have been developed to account for these reactions.
The first-order decay model has been shown to be sufficiently accurate for most distri-
bution system modeling applications and is well established. Rossman, Clark, and
Grayman (1994) proposed a mathematical framework for combining the complex
reactions occurring within distribution system pipes. This framework accounts for the
physical transport of the disinfectant from the bulk fluid to the pipe wall (mass trans-
fer effects) and the chemical reactions occurring there.
Section 2.9 Water Quality Modeling 59
Figure 2.25
Disinfectant reactions
occurring within a
typical distribution
system pipe
Q
i
Biofilm
Corrosion
Cl
Cl
Cl
Cl
Cl
Cl
N.O.M.
N.O.M.
N.O.M.
Cl
Cl
Cl
u C ( ) KC = (2.48)
where K = overall reaction rate constant (1/T)
Equation 2.48 is a simple first-order reaction (n = 1). The reaction rate coefficient K,
however, is now a function of the bulk reaction coefficient and the wall reaction coef-
ficient, as indicated in the following equation.
K k
b
k
w
k
f
R
H
k
w
k
f
+ ( )
----------------------------- + = (2.49)
where k
b
= bulk reaction coefficient (1/T)
k
w
= wall reaction coefficient (L/T)
k
f
= mass transfer coefficient, bulk fluid to pipe wall (L/T)
R
H
= hydraulic radius of pipeline (L)
The rate that disinfectant decays at the pipe wall depends on how quickly disinfectant
is transported to the pipe wall and the speed of the reaction once it is there. The mass
transfer coefficient is used to determine the rate at which disinfectant is transported
using the dimensionless Sherwood number, along with the molecular diffusivity coef-
ficient (of the constituent in water) and the pipeline diameter.
k
f
S
H
d
D
--------- = (2.50)
where S
H
= Sherwood number
d = molecular diffusivity of constituent in bulk fluid (L
2
/T)
D = pipeline diameter (L)
For stagnant flow conditions (Re < 1), the Sherwood number, S
H
, is equal to 2.0. For
turbulent flow (Re > 2,300), the Sherwood number is computed using Equation 2.51.
S
H
0.023Re
0.83 v
d
---
\ .
| |
0.333
= (2.51)
where Re = Reynolds number
v = kinematic viscosity of fluid (L
2
/T)
60 Modeling Theory Chapter 2
For laminar flow conditions (1 < Re < 2,300), the average Sherwood number along
the length of the pipe can be used. To have laminar flow in a 6-in. (150-mm) pipe, the
flow would need to be less than 5 gpm (0.3 l/s) with a velocity of 0.056 ft/s (0.017 m/
s). At such flows, head loss would be negligible.
S
H
3.65
0.0668
D
L
----
\ .
| |
Re ( )
v
d
---
\ .
| |
1 0.04
D
L
----
\ .
| |
Re
v
d
---
\ .
| |
2 3
+
------------------------------------------------------------ + = (2.52)
where L = pipe length (L)
Using the first-order reaction framework developed immediately above, both bulk
fluid and pipe wall disinfectant decay reactions can be accounted for. Bulk decay
coefficients can be determined experimentally. Wall decay coefficients, however, are
more difficult to measure and are frequently estimated using disinfectant concentra-
tion field measurements and water quality simulation results.
Formation Reactions. One shortcoming of the first-order reaction model is that
it accounts for the concentration of only one reactant. This model is sufficient if only
one reactant is being considered. For example, when chlorine residual concentrations
are modeled, chlorine is assumed to be the limiting reactant and the other reactants
material at the pipe walls and natural organic matter (NOM)are assumed to be
present in excess. The behavior of some disinfection by-product (DBP) formation
reactions, however, differs from this assumption. NOM, not chlorine, is frequently the
limiting reactant. DBP formation is just one example of a generalized class of reac-
tions that can be modeled using a limiting reactant. The reaction term for this class of
formation and decay reactions as proposed by Rossman (2000) is shown in Equation
2.53.
u C ( ) k C
li m
C ( )C
n 1
= (2.53)
where C
lim
= limiting concentration of the reaction (M/L
3
)
A first-order growth rate to a limiting value has been used to represent the formation
of trihalomethanes, a common form of DBP, in distribution systems (Vasconcelos et
al., 1996). Mathematically this is represented by Equation 2.54 and is shown graphi-
cally in Figure 2.26.
THM t ( ) C
o
FP C
o
| | 1 e
kt
| | + = (2.54)
where THM(t) = THM concentration at time t
C
o
= initial THM concentration
FP = formation potential (concentration)
k = reaction rate (a positive value)
Section 2.9 Water Quality Modeling 61
Figure 2.26
First-order growth
rate to a limiting value
Other Types of Water Quality Simulations
Although the water quality features of individual software packages vary, the most
common types of water quality simulations, in addition to the constituent analysis
already described, are source trace and water age analyses. The solution methods used
in both of these simulations are actually specific applications of the method used in
constituent analysis.
Source Trace Analysis. For the sake of reliability, or to simply provide suffi-
cient quantities of water to customers, a utility often uses more than one water supply
source. Suppose, for instance, that two treatment plants serve the same distribution
system. One plant draws water from a surface source, and the other pulls from an
underground aquifer. The raw water qualities from these sources are likely to differ
significantly, resulting in quality differences in the finished water as well.
Using a source trace analysis, the areas within the distribution system influenced by a
particular source can be determined, and, more important, areas where mixing of
water from different sources has occurred can be identified. The significance of
source mixing depends on the quality characteristics of the waters. Sometimes, mix-
ing can reduce the aesthetic qualities of the water (for example, creating cloudiness as
solids precipitate, or causing taste and odor problems to develop), and can contribute
to disinfectant residual maintenance problems. Source trace analyses are also useful
in tracking water quality problems related to storage tanks by tracing water from stor-
age as it is transported through the network.
A source trace analysis is a useful tool for better management of these situations. Spe-
cifically, it can be used to determine the percentage of water originating from a partic-
ular source for each junction node, tank, and reservoir in the distribution system
model. The procedure the software uses for this calculation is a special case of
constituent analysis in which the trace originates from the source as a conservative
constituent with an output concentration of 100 units. The constituent transport and
mixing equations introduced in the beginning of this section are then used to simulate
the transport pathways through the network and the influence of transport delays and
dilution on the trace constituent concentration. The values computed by the
simulation are then read directly as the percentage of water arriving from the source
location.
62 Modeling Theory Chapter 2
Water Age Analysis. The chemical processes that can affect distribution sys-
tem water quality are a function of water chemistry and the physical characteristics of
the distribution system itself (for example, pipe material and age). More generally,
however, these processes occur over time, making residence time in the distribution
system a critical factor influencing water quality. The cumulative residence time of
water in the system, or water age, has come to be regarded as a reliable surrogate for
water quality. Water age is of particular concern when quantifying the effect of stor-
age tank turnover on water quality. It is also beneficial for evaluating the loss of disin-
fectant residual and the formation of disinfection by-products in distribution systems.
The chief advantage of a water age analysis when compared to a constituent analysis
is that once the hydraulic model has been calibrated, no additional water quality cali-
bration procedures are required. The water age analysis, however, will not be as pre-
cise as a constituent analysis in determining water quality; nevertheless, it is an easy
way to leverage the information embedded in the calibrated hydraulic model. Con-
sider a project in which a utility is analyzing mixing in a tank and its effect on water
quality in an area of a network experiencing water quality problems. If a hydraulic
model has been developed and adequately calibrated, it can immediately be used to
Section 2.9 Water Quality Modeling 63
evaluate water age. The water age analysis may indicate that excessively long resi-
dence times within the tank are contributing to water quality degradation. Using this
information, a more precise analysis can be planned (such as an evaluation of tank
hydraulic dynamics and mixing characteristics, or a constituent analysis to determine
the impact on disinfectant residuals), and preliminary changes in design or operation
can be evaluated.
The water age analysis reports the cumulative residence time for each parcel of water
moving through the network. Again, the algorithm the software uses to perform the
analysis is a specialized case of constituent analysis. Water entering a network from a
source is considered to have an age of zero. The constituent analysis is performed
assuming a zero-order reaction with a k value equal to +1 [(mg/l)/s]. Thus, constituent
concentration growth is directly proportional to time, and the cumulative residence
time along the transport pathways in the network is numerically summed.
Using the descriptions of water quality transport and reaction dynamics provided
here, and the different types of water quality-related simulations available in modern
software packages, water quality in the distribution system can be accurately pre-
dicted. Water quality modeling can be used to help improve the performance of distri-
bution system modifications meant to reduce hydraulic residence times, and as a tool
for improving the management of disinfectant residuals and other water quality-
related operations. Continuing advancements in technology combined with more
stringent regulations on quality at the customers tap are motivating an increasing
number of utilities to begin using the powerful water quality modeling capabilities
already available to them.
Solution Methods
The earliest water quality models of distribution systems were steady-state models
(Wood, 1980 and Males, Clark, Wehrman, and Gates, 1985). These models used
simultaneous equation or marching out solution methods to determine the steady-
state water quality concentrations throughout the distribution systems. However, it
quickly became apparent that steady-state water quality models were of limited use in
representing actual systems due to the temporal variability in distribution system
operation, the impacts of tanks on water quality, and temporal changes in source con-
centrations. This led to the development of several dynamic water quality models dur-
ing the mid to late 1980s (Clark, Grayman, Males, and Coyle, 1986; Hart, Meader,
and Chiang, 1986; Liou and Kroon, 1987; and Grayman, Clark, and Males, 1988).
Two methods are available to solve the dynamic water quality equations used in water
quality models. One method is based on an Eulerian approach that divides each sepa-
rate pipe into a series of equal length sub-links. The other method is a Lagrangian
approach that tracks parcels of water of homogeneous water quality concentrations as
they move through the pipe system. These solution methods are shown graphically in
Figures 2.27 and 2.28 and explained in detail in the paragraphs that follow. In both
solution methods, a hydraulic model must first be applied in extended-period simula-
tion (EPS) mode to determine the flow, flow direction, and velocity in each pipe at all
times during the simulation.
64 Modeling Theory Chapter 2
The Eulerian Approach. With an Eulerian approach (illustrated in Figure
2.27), an observer located at a fixed location watches water as it flows by. Grayman,
Clark, and Males (1988) developed an Eulerian solution method for water quality
modeling in distribution systems, and Rossman, Boulos, and Altman (1993) formal-
ized this method and named it the Discrete Volume Method (DVM).
In DVM, for each time period, a pipe is divided into a series of sub-links with the sub-
link length selected so that the time of travel through each sub-link is equal to a user-
selected water quality time step that remains constant throughout the simulation. As a
result, water moves from one sub-link to the next adjacent downstream sub-link in
one water quality time step.
In order to meet this constraint, the sub-link length varies from pipe to pipe and within
a pipe as flow changes. If the constituent being simulated is reactive, then the water
quality concentration is adjusted according to the appropriate reaction method during
each water quality time step.
At a junction at the downstream end of one or more pipes, the water quality concen-
tration in the junction is calculated by taking a flow-weighted average of incoming
inflows, as described previously by Equation 2.42. Water moves instantaneously
through pumps and valves without a change in water quality. At the end of a hydraulic
time step, if there is a change in flow or direction, then the sub-link gridding is
changed, and the water quality concentrations at the end of the previous time step are
used to define the initial water quality in each of the new sub-links. There are special
numerical assumptions made to accommodate problem situations such as very short
pipes (with travel time less than the water quality time step) and very long pipes (with
a very large number of sub-links).
Figure 2.27
Eulerian solution
method
The Lagrangian Approach. In the Lagrangian approach (illustrated in Figure
2.28), rather than observing the flow from the sidelines, the observer moves with
the flow. Additionally, rather than having a fixed grid, parcels of water with homoge-
neous concentrations are tracked through the pipe. New parcels are added when water
quality changes occur due to changes in source quality or when parcels are combined
at junctions. In order to reduce the number of parcels, algorithms have been devel-
oped for combining adjacent parcels when the difference in concentrations are less
than a user-defined tolerance. Liou and Kroon (1987) and Hart, Meader, and Chiang
(1986) developed Lagrangian solution methods for water distribution system water
quality models.
Lagrangian solutions can be either time-driven or event-driven. In a time-driven
method, conditions are updated at a fixed time step. In an event-driven model, condi-
tions are updated when the source water quality changes or when the front of a parcel
References 65
reaches a junction. In comparing the methods, Rossman and Boulos (1996) found that
the Lagrangian time-driven method is the most efficient and versatile of the solution
methods for water distribution system quality models.
Figure 2.28
Lagrangian solution
method
REFERENCES
ASCE. (1975). Pressure Pipeline Design for Water and Wastewater. ASCE, New York, New York.
ASCE/WEF. (1982). Gravity Sanitary Sewer Design and Construction. ASCE, Reston, Virginia.
ASCE Committee on Pipeline Planning. (1992). Pressure Pipeline Design for Water and Wastewater.
ASCE, Reston, Virginia.
Bernoulli, D. (1738). Hydrodynamica. Argentorati.
Bhave, P. R. (1991). Analysis of Flow in Water Distribution Networks. Technomics, Lancaster, Pennsylva-
nia.
Bosserman, B. E., (2000). Pump System Hydraulic Design. Water Distribution System Handbook, Mays,
L. W., ed., McGraw-Hill, New York, New York.
Clark, R. M., and Grayman, W. M. (1998). Modeling Water Quality in Distribution Systems. AWWA, Den-
ver, Colorado.
Clark, R. M., Grayman, W. M., Males, R. M., and Coyle, J. A. (1986). Predicting Water Quality in Distri-
bution Systems. Proceedings of the AWWA Distribution System Symposium, American Water Works
Association, Denver, Colorado.
Crane Company (1972). Flow of Fluids through Valves and Fittings. Crane Co., New York, New York.
Cross, H. (1936). Analysis of Flow in Networks of Conduits or Conductors. University of Illinois Experi-
ment Station Bulletin No. 286, Department of Civil Engineering, University of Illinois, Champaign
Urbana, Illinois.
Grayman, W. M., Clark, R. M., and Males, R. M. (1988). Modeling Distribution System Water Quality:
Dynamic Approach. Journal of Water Resources Planning and Management, ASCE, 114(3).
Grayman, W. M., Deininger, R. A., Green, A., Boulos, P. F., Bowcock, R. W., and Godwin, C. C. (1996).
Water Quality and Mixing Models for Tanks and Reservoirs. Journal of the American Water Works
Association, 88(7).
Grayman, W. M., Rossman, L. A., and Geldreich, E. E. (2000). Water Quality. Water Distribution Sys-
tems Handbook, Mays, L. W., ed., McGraw-Hill, New York, New York.
Grayman, W. M., Rossman, L. A., Arnold, C., Deininger, R. A., Smith, C., Smith, J. F., and Schnipke, R.
(2000). Water Quality Modeling of Distribution System Storage Facilities. AWWA.
66 Modeling Theory Chapter 2
Hart, F. F., Meader, J. L., and Chiang, S. N. (1986). CLNETA Simulation Model for Tracing Chlorine
Residuals in a Potable Water Distribution Network. Proceedings of the AWWA Distribution System
Symposium, American Water Works Association, Denver, Colorado.
Hydraulic Institute (2000). Pump Standards. Parsippany, New Jersey.
Idelchik, I. E. (1999). Handbook of Hydraulic Resistance. 3rd edition, Begell House, New York, New York.
Lamont, P. A. (1981). Common Pipe Flow Formulas Compared with the Theory of Roughness. Journal
of the American Water Works Association, 73(5), 274.
Lansey, K., and Mays, L. W. (2000). Hydraulics of Water Distribution Systems. Water Distribution Sys-
tems Handbook, Mays, L. W., ed., McGraw-Hill, New York, New York.
Larock, B. E., Jeppson, R. W., and Watters, G. Z. (1999). Handbook of Pipeline Systems. CRC Press, Boca
Raton, Florida.
Liou, C. P. and Kroon, J. R. (1987). Modeling the Propagation of Waterborne Substances in Distribution
Networks. Journal of the American Water Works Association, 79(11), 54.
Karassik, I. J., ed. (1976). Pump Handbook. McGraw-Hill, New York, New York.
Males, R. M., Clark, R. M., Wehrman, P. J., and Gates, W. E. (1985). An Algorithm for Mixing Problems
in Water Systems. Journal of Hydraulic Engineering, ASCE, 111(2).
Mays, L. W., ed. (1999). Hydraulic Design Handbook. McGraw-Hill, New York, New York.
Miller, D. S. (1978). Internal Flow Systems. BHRA Fluid Engineering, Bedford, United Kingdom.
Moody, L. F. (1944). Friction Factors for Pipe Flow. Transactions of the American Society of Mechanical
Engineers, Vol. 66.
Nikuradse (1932). Gestezmassigkeiten der Turbulenten Stromung in Glatten Rohren. VDI-Forschungsh,
No. 356 (in German).
Rossman, L.A. (2000). EPANET Users Manual. Risk Reduction Engineering Laboratory, U.S. Environ-
mental Protection Agency, Cincinnati, Ohio.
Rossman, L. A., and Boulos, P. F. (1996). Numerical Methods for Modeling Water Quality in Distribution
Systems: A Comparison. Journal of Water Resources Planning and Management, ASCE, 122(2), 137.
Rossman, L. A., Boulos, P. F., and Altman, T. (1993). Discrete Volume-Element Method for Network
Water-Quality Models. Journal of Water Resources Planning and Management, ASCE, 119(5), 505.
Rossman, L. A., Clark, R. M., and Grayman, W. M. (1994). Modeling Chlorine Residuals in Drinking
Water Distribution Systems. Journal of Environmental Engineering, ASCE, 1210(4), 803.
Sanks, R. L., ed. (1998). Pumping Station Design. 2nd edition, Butterworth, London, UK.
Streeter, V. L., Wylie, B. E., and Bedford, K. W. (1998). Fluid Mechanics. 9th edition, WCB/McGraw-Hill,
Boston, Massachusetts.
Swamee, P. K., and Jain, A. K. (1976). Explicit Equations for Pipe Flow Problems. Journal of Hydraulic
Engineering, ASCE, 102(5), 657.
Todini, E., and Pilati, S. (1987). A Gradient Method for the Analysis of Pipe Networks. Proceedings of
the International Conference on Computer Applications for Water Supply and Distribution, Leicester
Polytechnic, UK.
Vasconcelos, J. J., Boulos, P. F., Grayman, W. M., Kiene, L., Wable, O., Biswas, P., Bhari, A., Rossman, L.,
Clark, R., and Goodrich, J. (1996). Characterization and Modeling of Chlorine Decay in Distribution
Systems. AWWA, Denver, Colorado.
Walski, T. M. (1984). Analysis of Water Distribution Systems. Van Nostrand Reinhold, New York, New
York.
Williams, G. S., and Hazen, A. (1920). Hydraulic Tables. John Wiley & Sons, New York, New York.
Wood, D. J. (1980). Slurry Flow in Pipe Networks. Journal of Hydraulics, ASCE, 106(1), 57.
Discussion Topics and Problems 67
DISCUSSION TOPICS AND PROBLEMS
Read the chapter and complete the problems. Submit your work to Haestad Methods
and earn up to 11.0 CEUs. See Continuing Education Units on page xxix for more
information, or visit www.haestad.com/awdm-ceus/.
2.1 Find the viscosity of the fluid contained between the two square plates shown in the figure. The top
plate is moving at a velocity of 3 ft/s.
t = 0.5 in.
F = 50 lb
6 ft
2.2 Find the force P required to pull the 150 mm circular shaft shown in the figure through the sleeve at
a velocity of 1.5 m/s. The fluid between the shaft and the sleeve is water at a temperature of 15
o
C.
P
150 mm
2 mm
75 mm
2.3 Find the pressure at the base of a container of water having a depth of 15 m.
2.4 How high is the water level from the base of an elevated storage tank if the pressure at the base of the
tank is 45 psi?
2.5 Water having a temperature of 65
o
F is flowing through a 6-in. ductile iron main at a rate of 300 gpm.
Is the flow laminar, turbulent, or transitional?
2.6 What type of flow do you think normally exists in water distribution systems: laminar, turbulent, or
transitional? Justify your selection with sound reasoning.
68 Modeling Theory Chapter 2
2.7 What is the total head at point A in the system shown in the figure if the flow through the pipeline is
1,000 gpm? What is the head loss in feet between point A and point B?
Q 8 in.
Datum
550 ft 550 ft
P = 62 psi
A
P = 48 psi
B
A B
2.8 For the piping system shown in Problem 2.7, what would the elevation at point B have to be in order
for the reading on the two pressure gages to be the same?
2.9 Assuming that there are no head losses through the Venturi meter shown in the figure, what is the
pressure reading in the throat section of the Venturi? Assume that the discharge through the meter is
158 l/s.
P = 497kPa
P = ? kPa
400 mm
150 mm
2.10 What is the head loss through a 10-in. diameter concrete water main 2,500 ft in length if water at
60
o
F is flowing through the line at a rate of 1,250 gpm? Solve using the Darcy-Weisbach formula.
2.11 For Problem 2.10, what is the flow through the line if the head loss is 32 ft? Solve using the Darcy-
Weisbach formula.
2.12 Find the length of a pipeline that has the following characteristics: Q=41 l/s, D=150 mm, Hazen-
Williams C=110, H
L
=7.6 m.
2.13 For the pipeline shown in Problem 2.7, what is the Hazen-Williams C-factor if the distance between
the two pressure gages is 725 ft and the flow is 1,000 gpm?
Discussion Topics and Problems 69
2.14 English Units: Compute the pipe resistance coefficient, K
P
, for the following pipelines.
Length
(ft)
Diameter
(in.)
Hazen-Williams
C-factor
Pipe Resistance
Coefficient (K
p
)
1,200 12 120
500 4 90
75 3 75
3,500 10 110
1,750 8 105
SI Units: Compute the pipe resistance coefficient, K
P
, for the following pipelines.
Length
(m)
Diameter
(mm)
Hazen-Williams
C-factor
Pipe Resistance
Coefficient (K
p
)
366 305 120
152 102 90
23 76 75
1067 254 110
533 203 105
2.15 English Units: Compute the minor loss term, K
M
, for the fittings shown in the table below.
Type of Fitting/Flow Condition
Minor Loss
Coefficient
Pipe Size
(in.)
Minor Loss Term
(K
M
)
Gate Valve - 50% Open 4.8 8
Tee - Line Flow 0.4 12
90
o
Mitered Bend 0.8 10
Fire Hydrant 4.5 6
SI Units: Compute the minor loss term, K
M
, for the fittings shown in the table below.
Type of Fitting/Flow Condition
Minor Loss
Coefficient
Pipe Size
(mm)
Minor Loss Term
(K
M
)
Gate Valve - 50% Open 4.8 200
Tee - Line Flow 0.4 300
90
o
Mitered Bend 0.8 250
Fire Hydrant 4.5 150
2.16 English Units: Determine the pressures at the following locations in a water distribution system,
assuming that the HGL and ground elevations at the locations are known.
Node Label
HGL
(ft)
Elevation
(ft)
Pressure
(psi)
J-1 550.6 423.5
J-6 485.3 300.5
J-23 532.6 500.0
J-5 521.5 423.3
J-12 515.0 284.0
70 Modeling Theory Chapter 2
SI Units: Determine the pressures at the following locations in a water distribution system, assuming
that the HGL and ground elevations at the locations are known.
Node Label
HGL
(m)
Elevation
(m)
Pressure
(kPa)
J-1 167.8 129.1
J-6 147.9 91.6
J-23 162.3 152.4
J-5 159.0 129.0
J-12 157.0 86.6
2.17 Using the concept of conservation of mass, is continuity maintained at the junction node shown in
the following figure?
J-10 P -8
P -12
P -10
P -9
Q = 38 gpm
P-10
Q = 120 gpm
P-12
D = 45 gpm
J-10
Q = 72 gpm
P-9
Q = 55 gpm
P-8
2.18 Find the magnitude and direction of the flow through pipe P-9 so that continuity is maintained at
node J-10 in the following figure.
J-10 P -8
P -12
P -10
P -9
Q = 3.0 L/s
P-10
Q = ? L/s
P-9
D = 0.9 L/s
J-10
Q = 1.9 L/s
P-12
Q = 1.4 L/s
P-8
Discussion Topics and Problems 71
2.19 Does conservation of energy around a loop apply to the loop shown in the following figure? Why or
why not? The total head loss (sum of friction losses and minor losses) in each pipe and the direction
of flow are shown in the figure.
P-23
P-32
P-27
P-25
H = 7.73 ft
L,P-23
H = 3.76 ft
L,P-32
H = 2.10 ft
L,P-27
H = 1.87 ft
L,P-25
2.20 Does Conservation of Energy apply to the system shown in the figure? Data describing the physical
characteristics of each pipe are presented in the table below. Assume that there are no minor losses in
this loop.
P-23
P-32
P-27
P-25
Q = 22.7 L/s
P-23
Q = 4.7 L/s
P-32
Q = 27.1 L/s
P-27
Q = 7.6 L/s
P-25
Pipe Label
Length
(m)
Diameter
(mm)
Hazen-Williams
C-factor
P-23 381.0 305 120
P-25 228.6 203 115
P-27 342.9 254 120
P-32 253.0 152 105
72 Modeling Theory Chapter 2
2.21 Find the discharge through the system shown in the figure. Compute friction loss using the Hazen-
Williams equation.
R-1
R-2
Pipe 2:
L=5,000 ft
D=4 in.
C=85
125 ft
225 ft
2.22 Find the pump head needed to deliver water from reservoir R-1 to reservoir R-2 in the following fig-
ure at a rate of 70.8 l/s. Compute friction losses using the Hazen-Williams equation.
R-1
R-2
Pipe 2:
L=1524 m
D=102 mm
C=85
Pipe 1:
L=152 m
D=305 mm
C=120
38.1 m
68.6 m
2.23 Compute the age of water at the end of a 12-in. pipe that is 1,500 ft in length and has a flow of 900
gpm. The age of the water when it enters the pipeline is 7.2 hours.
2.24 Suppose that a 102-mm pipe is used to serve a small cluster of homes at the end of a long street. If
the length of the pipe is 975 m, what is the age of water leaving it if the water had an age of 6.3 hours
when entering the line? Assume that the water use is 1.6 l/s.
Discussion Topics and Problems 73
2.25 Given the data in the tables below, what is the average age of the water leaving junction node J-4
shown in the figure? What is the flow rate through pipe P-4? What is the average age of the water
arriving at node J-5 through pipe P-4? Fill in your answers in the tables provided.
J -1
J -2
J-3
J-4
J-5
P -1
P -3
P -2
P -4
35 gpm
Pipe Label
Flow
(gpm)
Length
(ft)
Diameter
(in.)
P-1 75 1,650 10
P-2 18 755 8
P-3 23 820 6
P-4 2,340 10
Node Label
Average Age
(hours)
J-1 5.2
J-2 24.3
J-3 12.5
J-4
J-5
2.26 What will be the concentration of chlorine in water samples taken from a swimming pool after 7
days if the initial chlorine concentration in the pool was 1.5 mg/l? Bottle tests performed on the pool
water indicate that the first-order reaction rate is -0.134 day
-1
.
2.27 Do you think that the actual reaction rate coefficient for water in the swimming pool described above
(i.e., the water being considered remains in the pool, and is not being stored under laboratory condi-
tions) would be equal to -0.134 day
-1
? Suggest some factors that might cause the actual reaction rate
to differ. Would these factors most likely cause the actual reaction rate to be greater than or less than
-0.134 day
-1
?
74 Modeling Theory Chapter 2
2.28 For the system presented in Problem 2.25, what is the concentration of a constituent leaving node J-
4 (assume it is a conservative constituent)? The constituent concentration is 0.85 mg/L in pipe P-1,
0.50 mg/l in pipe P-2, and 1.2 mg/l in pipe P-3.
2.29 What is the fluoride concentration at the end of a 152-mm diameter pipeline 762 m in length if the
fluoride concentration at the start of the line is 1.3 mg/l? Fluoride is a conservative species; that is, it
does not decay over time. Ignoring dispersion, if there is initially no fluoride in the pipe and it is
introduced at the upstream end at a 2.0 mg/l concentration, when will this concentration be reached
at the end of the line if the flow through the pipe is 15.8 l/s? Assume that there are no other junction
nodes along the length of this pipe.
2.30 A community has found high radon levels in one of their wells. Because radon decays relatively
quickly, they are exploring the use of a baffled clearwell at the well in order to provide some time
delay before the water enters the distribution system. Their goal is to reduce the radon levels by 80
percent in the clearwell. The half-life for radon is 3.8 days. What minimum residence time is
required in the clearwell to meet this goal? Is the use of the clearwell as a means of meeting their
goal reasonable?
2.31 Trihalomethane levels leaving a treatment plant are 20 ug/l. Based on bottle tests, the formation was
found to follow a first-order reaction with a growth rate of 2 l/day and an ultimate formation poten-
tial of 100 ug/l. What would the THM concentration be after 1 day? How long would it take for the
THM concentration to reach 0.99 ug/l?
2.32 A city is building a 6-in. diameter, 5000-ft pipeline to serve a small community with an average con-
sumption of 20 gpm. Based on tests, they estimate that the bulk decay rate for chlorine is 0.5 1/day
and the wall decay rate will be 1 feet per day. If the chlorine residual at the start of the pipe is
1 mg/l, what will the chlorine residual be in the water delivered to the community? If you only con-
sidered the effects of bulk decay, what would the chlorine residual be? If you only considered wall
decay, what would the chlorine residual be? (Hint: Using a distribution system model, build a net-
work model composed of a reservoir, a single pipe, and a junction representing the community.)




ADVANCED WATER
DISTRIBUTION MODELING
AND MANAGEMENT

Authors
Thomas M. Walski
Donald V. Chase
Dragan A. Savic
Walter Grayman
Stephen Beckwith
Edmundo Koelle

Contributing Authors
Scott Cattran, Rick Hammond, Kevin Laptos, Steven G. Lowry,
Robert F. Mankowski, Stan Plante, John Przybyla, Barbara Schmitz

Peer Review Board
Lee Cesario (Denver Water), Robert M. Clark (U.S. EPA),
Jack Dangermond (ESRI), Allen L. Davis (CH2M Hill),
Paul DeBarry (Borton-Lawson), Frank DeFazio (Franklin G. DeFazio Corp.),
Kevin Finnan (Bristol Babcock), Wayne Hartell (Bentley Systems),
Brian Hoefer (ESRI), Bassam Kassab (Santa Clara Valley Water District),
James W. Male (University of Portland), William M. Richards
(WMR Engineering), Zheng Wu (Bentley Systems ),
and E. Benjamin Wylie (University of Michigan)



Click here to visit the Bentley Institute
Press Web page for more information

C H A P T E R
3
Assembling a Model
As Chapter 1 discusses, a water distribution model is a mathematical description of a
real-world system. Before building a model, it is necessary to gather information
describing the network. In this chapter, we introduce and discuss sources of data used
in constructing models.
The latter part of the chapter covers model skeletonization. Skeletonization is the pro-
cess of simplifying the real system for model representation, and it involves making
decisions about the level of detail to be included.
3.1 MAPS AND RECORDS
Many potential sources are available for obtaining the data required to generate a
water distribution model, and the availability of these sources varies dramatically
from utility to utility. The following sections discuss some of the most commonly
used resources, including system maps, as-built drawings, and electronic data files.
System Maps
System maps are typically the most useful documents for gaining an overall under-
standing of a water distribution system because they illustrate a wide variety of
valuable system characteristics. System maps may include such information as
Pipe alignment, connectivity, material, diameter, and so on
The locations of other system components, such as tanks and valves
Pressure zone boundaries
Elevations
Miscellaneous notes or references for tank characteristics
Background information, such as the locations of roadways, streams, plan-
ning zones, and so on
Other utilities
76 Assembling a Model Chapter 3
Topographic Maps
A topographic map uses sets of lines called contours to indicate elevations of the
ground surface. Contour lines represent a contiguous set of points that are at the same
elevation and can be thought of as the outline of a horizontal slice of the ground sur-
face. Figure 3.1 illustrates the cross-sectional and topographic views of a sphere, and
Figure 3.2 shows a portion of an actual topographic map. Topographic maps are often
referred to by the contour interval that they present, such as a 20-foot topographic
map or a 1-meter contour map.
By superimposing a topographic map on a map of the network model, it is possible to
interpolate the ground elevations at junction nodes and other locations throughout the
system. Of course, the smaller the contour interval, the more precisely the elevations
can be estimated. If available topographic maps cannot provide the level of precision
needed, other sources of elevation data need to be considered.
Topographic maps are also available in the form of Digital Elevation Models (DEMs),
which can be used to electronically interpolate elevations. The results of the DEM are
only as accurate as the underlying topographic data on which they are based; thus, it is
possible to calculate elevations to a large display precision but with no additional
accuracy.
Figure 3.1
Topographic
representation of a
hemisphere
Looking Down
Generates
Plan View
Profile
Plan
Looking from
the Side
Generates
Profile View
0
100
200
2
0
0
1
0
0
0
4
0
0
3
0
0
300
400
As-Built Drawings
Site restrictions and on-the-fly changes often result in differences between original
design plans and the actual constructed system. As a result, most utilities perform
post-construction surveys and generate a set of as-built or record drawings for the
purpose of documenting the system exactly as it was built. In some cases, an inspec-
tor's notes may even be used as a supplemental form of documentation. As-built
drawings can be especially helpful in areas where a fine level of precision is required
for pipe lengths, fitting types and locations, elevations, and so forth.
Section 3.1 Maps and Records 77
Figure 3.2
Typical topographic
map
As-built drawings can also provide reliable descriptions of other system components
such as storage tanks and pumping stations. There may be a complete set of drawings
for a single tank, or the tank plans could be included as part of a larger construction
project.
Electronic Maps and Records
Many water distribution utilities have some form of electronic representation of their
systems in formats that may vary from a nongraphical database, to a graphics-only
Computer-Aided Drafting (CAD) drawing, to a Geographic Information System (GIS)
that combines graphics and data.
Nongraphical Data. It is common to find at least some electronic data in non-
graphical formats, such as a tracking and inventory database, or even a legacy text-
based model. These sources of data can be quite helpful in expediting the process of
model construction. Even so, care needs to be taken to ensure that the network topol-
ogy is correct, because a simple typographic error in a nongraphical network can be
difficult to detect.
Computer-Aided Drafting. The rise of computer technology has led to many
improvements in all aspects of managing a water distribution utility, and mapping is
no exception. CAD systems make it much easier to plug in survey data, combine data
from different sources, and otherwise maintain and update maps faster and more reli-
ably than ever before.
78 Assembling a Model Chapter 3
Even for systems having only paper maps, many utilities digitize those maps to con-
vert them to an electronic drawing format. Traditionally, digitizing has been a process
of tracing over paper maps with special computer peripherals, called a digitizing tab-
let and puck (see Figure 3.3). A paper map is attached to the tablet, and the drafts-
person uses crosshairs on the puck to point at locations on the paper. Through
magnetic or optical techniques, the tablet creates an equivalent point at the
appropriate location in the CAD drawing. As long as the tablet is calibrated correctly,
it will automatically account for rotation, skew, and scale.
Figure 3.3
A typical digitizing
tablet
Another form of digitizing is called heads-up digitizing (see Figure 3.4). This method
involves scanning a paper map into a raster electronic format (such as a bitmap),
bringing it into the background of a CAD system, and electronically tracing over it on
a different layer. The term heads-up is used because the draftsperson remains focused
on the computer screen rather than on a digitizing tablet.
Geographic Information Systems. A Geographic information system
(GIS) is a computer-based tool for mapping and analyzing objects and events that
happen on earth. GIS technology integrates common database operations such as
query and statistical analysis with the unique visualization and geographic analysis
benefits offered by maps (ESRI, 2001). Because a GIS stores data on thematic layers
linked together geographically, disparate data sources can be combined to determine
relationships between data and to synthesize new information.
Section 3.2 Model Representation 79
Figure 3.4
Network model
overlaid on an aerial
photograph
GIS can be used for tasks such as proximity analysis (identifying customers within a
certain distance of a particular node), overlay analysis (determining all junctions that
are completely within a particular zoning area), network analysis (identifying all
households impacted by a water-main break), and visualization (displaying and com-
municating master plans graphically). With a hydraulic model that links closely to a
GIS, the benefits can extend well beyond just the process of building the model and
can include skeletonization, demand generalization, and numerous other operations.
3.2 MODEL REPRESENTATION
The concept of a network is fundamental to a water distribution model. The network
contains all of the various components of the system, and defines how those elements
are interconnected. Networks are comprised of nodes, which represent features at spe-
cific locations within the system, and links, which define relationships between nodes.
80 Assembling a Model Chapter 3
Network Elements
Water distribution models have many types of nodal elements, including junction
nodes where pipes connect, storage tank and reservoir nodes, pump nodes, and con-
trol valve nodes. Models use link elements to describe the pipes connecting these
nodes. Also, elements such as valves and pumps are sometimes classified as links
rather than nodes. Table 3.1 lists each model element, the type of element used to rep-
resent it in the model, and the primary modeling purpose.
Table 3.1 Common network modeling elements
Element Type Primary Modeling Purpose
Reservoir Node Provides water to the system
Tank Node Stores excess water within the system and releases that water
at times of high usage
Junction Node Removes (demand) or adds (inflow) water from/to
the system
Pipe Link Conveys water from one node to another
Pump Node
or link
Raises the hydraulic grade to overcome elevation differences
and friction losses
Control
Valve
Node
or link
Controls flow or pressure in the system based on specified
criteria
Naming Conventions (Element Labels). Because models may contain tens
of thousands of elements, naming conventions are an important consideration in mak-
ing the relationship between real-world components and model elements as obvious
as possible (see Figure 3.5). Some models allow only numeric numbering of ele-
ments, but most modern models support at least some level of alphanumeric labeling
(for example, J-1, Tank 5, or West Side Pump A).
Figure 3.5
Schematic junction
with naming
convention
J5- 115- ElmStreet
Description
115 = Sequential Number
5 = Zone 5
J = Junction Node
Naming conventions should mirror the way the modeler thinks about the particular
network by using a mixture of prefixes, suffixes, numbers, and descriptive text. In
general, labels should be as short as possible to avoid cluttering a drawing or report,
but they should include enough information to identify the element. For example, a
naming convention might include a prefix for the element type, another prefix to indi-
cate the pressure zone or map sheet, a sequential number, and a descriptive suffix.
Section 3.2 Model Representation 81
Of course, modelers can choose to use some creativity, but it is important to realize
that a name that seems obvious today may be baffling to future users. Intelligent use
of element labeling can make it much easier for users to query tabular displays of
model data with filtering and sorting commands. In some cases, such as automated
calibration, it may be very helpful to group pipes with like characteristics to make cal-
ibration easier. If pipe labels have been set up such that like pipes have similar labels,
this grouping becomes easy.
Rather than starting pipe labeling at a random node, it is best to start from the water
source and number outward along each pipeline. In addition, just as pipe elements
were not laid randomly, a pipe labeling scheme should be developed to reflect that.
For example, consider the pipes in Figure 3.6 (Network A), which shows that the
pipes were laid in four separate projects in four different years. By labeling the pipes
as shown in Figure 3.6 (Network B), the user will be able to more rationally group,
filter, and sort pipes. For example, pipes laid during the 1974 construction project
were labeled P-21, P-22, and so on so that those pipes could be grouped together. This
can have major time-saving benefits in working with a large system.
Figure 3.6
Logical element
labeling schemes
82 Assembling a Model Chapter 3
Boundary Nodes. A boundary node is a network element used to represent
locations with known hydraulic grade elevations. A boundary condition imposes a
requirement within the network that simulated flows entering or exiting the system
agree with that hydraulic grade. Reservoirs (also called fixed grade nodes) and tanks
are common examples of boundary nodes.
Every model must have at least one boundary node so that there is a reference point
for the hydraulic grade. In addition, every node must maintain at least one path back
to a boundary node so that its hydraulic grade can be calculated. When a node
becomes disconnected from a boundary (as when pipes and valves are closed), it can
result in an error condition that needs to be addressed by the modeler.
Network Topology
The most fundamental data requirement is to have an accurate representation of the
network topology, which details what the elements are and how they are intercon-
nected. If a model does not faithfully duplicate real-world layout (for example, the
model pipe connects two nodes that are not really connected), then the model will
never accurately depict real-world performance, regardless of the quality of the
remaining data.
System maps are generally good sources of topological information, typically includ-
ing data on pipe diameters, lengths, materials, and connections with other pipes.
There are situations in which the modeler must use caution, however, because maps
may be imperfect or unclear.
False Intersections. Just because mains appear to cross on a map does not nec-
essarily mean that a hydraulic connection exists at that location. As illustrated in Fig-
ure 3.7, it is possible for one main to pass over the other (called a crossover).
Modeling this location as an intersecting junction node would be incorrect, and could
result in serious model inaccuracies. Note that some GISs automatically assign nodes
where pipes cross, which may not be hydraulically correct.
When pipes are connected in the field via a bypass (as illustrated in Figure 3.7), the
junction node should only be included in the model if the bypass line is open. Since
the choice to include or omit a junction in the model based on the open or closed sta-
tus of a bypass in the field is somewhat difficult to control, it is recommended that the
bypass itself be included in the model. As a result, the modeler can more easily open
or close the bypass in accordance with the real system.
Figure 3.7
Pipe crossover and
crossover with bypass
Cross Crossover
Crossover w/
Bypass Line
Section 3.2 Model Representation 83
Converting CAD Drawings into Models. Although paper maps can some-
times falsely make it appear as though there is a pipe intersection, CAD maps can
have the opposite problem. CAD drawings are often not created with a hydraulic
model in mind; thus, lines representing pipes may visually appear to be connected on
a large-scale plot, but upon closer inspection of the CAD drawing, the lines are not
actually touching. Consider Figure 3.8, which demonstrates three distinct conditions
that may result in a misinterpretation of the topology:
T-intersections: Are there supposed to be three intersecting pipes or two
non-intersecting pipes? The drawing indicates that there is no intersection,
but this could easily be a drafting error.
Crossing pipes: Are there supposed to be four intersecting pipes or two non-
intersecting pipes?
Nearly connecting line endpoints: Are the two pipes truly non-intersecting?
Automated conversion from CAD drawing elements to model elements can save time,
but (as with any automated process) the modeler needs to be aware of the potential
pitfalls involved and should review the end result. Some models assist in the review
process by highlighting areas with potential connectivity errors. The possibility of
difficult-to-detect errors still remains, however, persuading some modelers to trace
over CAD drawings when creating model elements.
84 Assembling a Model Chapter 3
Figure 3.8
Common CAD
conversion errors
3.3 RESERVOIRS
The term reservoir has a specific meaning with regard to water distribution system
modeling that may differ slightly from the use of the word in normal water distribu-
tion construction and operation. A reservoir represents a boundary node in a model
that can supply or accept water with such a large capacity that the hydraulic grade of
the reservoir is unaffected and remains constant. It is an infinite source, which means
that it can theoretically handle any inflow or outflow rate, for any length of time,
without running dry or overflowing. In reality, there is no such thing as a true infinite
source. For modeling purposes, however, there are situations where inflows and out-
flows have little or no effect on the hydraulic grade at a node.
Reservoirs are used to model any source of water where the hydraulic grade is con-
trolled by factors other than the water usage rate. Lakes, groundwater wells, and
clearwells at water treatment plants are often represented as reservoirs in water distri-
bution models. For modeling purposes, a municipal system that purchases water from
a bulk water vendor may model the connection to the vendors supply as a reservoir
(most current simulation software includes this functionality).
For a reservoir, the two pieces of information required are the hydraulic grade line
(water surface elevation) and the water quality. By model definition, storage is not a
concern for reservoirs, so no volumetric storage data is needed.
3.4 TANKS
A storage tank (see Figure 3.9) is also a boundary node, but unlike a reservoir, the
hydraulic grade line of a tank fluctuates according to the inflow and outflow of water.
Tanks have a finite storage volume, and it is possible to completely fill or completely
exhaust that storage (although most real systems are designed and operated to avoid
such occurrences). Storage tanks are present in most real-world distribution systems,
and the relationship between an actual tank and its model counterpart is typically
straightforward.
Section 3.4 Tanks 85
Figure 3.9
Storage tanks
For steady-state runs, the tank is viewed as a known hydraulic grade elevation, and the
model calculates how fast water is flowing into or out of the tank given that HGL.
Given the same HGL setting, the tank is hydraulically identical to a reservoir for a
steady-state run. In extended-period simulation (EPS) models, the water level in the
tank is allowed to vary over time. To track how a tanks HGL changes, the relationship
between water surface elevation and storage volume must be defined. Figure 3.10 illus-
trates this relationship for various tank shapes. For cylindrical tanks, developing this
relationship is a simple matter of identifying the diameter of the tank, but for non-
cylindrical tanks it can be more challenging to express the tanks characteristics.
Some models do not support noncylindrical tanks, forcing the modeler to approximate
the tank by determining an equivalent diameter based on the tanks height and capac-
ity. This approximation, of course, has the potential to introduce significant errors in
hydraulic grade. Fortunately, most models do support non-cylindrical tanks, although
the exact set of data required varies from model to model.
Regardless of the shape of the tank, several elevations are important for modeling pur-
poses. The maximum elevation represents the highest fill level of the tank, and is
usually determined by the setting of the altitude valve if the tank is equipped with one.
The overflow elevation, the elevation at which the tank begins to overflow, is slightly
higher. Similarly, the minimum elevation is the lowest the water level in the tank
should ever be. A base or reference elevation is a datum from which tank levels are
measured.
The HGL in a tank can be referred to as an absolute elevation or a relative level,
depending on the datum used. For example, a modeler working near the Mile High
city of Denver, Colorado, could specify a tanks base elevation as the datum, and then
work with HGLs that are relative to that datum. Alternatively, the modeler could work
with absolute elevations that are in the thousands of feet. The choice of whether to use
absolute elevations or relative tank levels is a matter of personal preference. Figure
86 Assembling a Model Chapter 3
3.11 illustrates these important tank elevation conventions for modeling tanks. Notice
that when using relative tank levels, it is possible to have different values for the same
level, depending on the datum selected.
Figure 3.10
Volume versus level
curves for various
tank shapes
V
o
l
u
m
e
R
a
t
i
o
Depth Ratio
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0
1.0
0.9
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0
Spherical
45 Deg. Cone
Cylindrical/
Rectangular
Water storage tanks can be classified by construction material (welded steel, bolted
steel, reinforced concrete, prestressed concrete), shape (cylindrical, spherical, torroi-
dal, rectangular), style (elevated, standpipe, ground, buried), and ownership (utility,
private) (Walski, 2000). However, for pipe network modeling, the most important
classification is whether or not the tank floats on the system. A tank is said to float
on the system if the hydraulic grade elevation inside the tank is the same as the HGL
in the water distribution system immediately outside of the tank. With tanks, there are
really three situations that a modeler can encounter:
1. Tank that floats on the system with a free surface
2. Pressure (hydropneumatic) tank that floats on the system
3. Pumped storage in which water must be pumped from a tank
Section 3.4 Tanks 87
Figure 3.11
Important tank
elevations
Figure 3.12 shows that elevated tanks, standpipes, and hydropneumatic tanks float on
the system because their HGL is the same as that of the system. Ground tanks and
buried tanks may or may not float on the system, depending on their elevation. If the
HGL in one of these tanks is below the HGL in the system, water must be pumped
from the tank, resulting in pumped storage.
A tank with a free surface floating on the system is the simplest and most common
type of tank. The pumped storage tank needs a pump to deliver water from the tank to
the distribution system and a control valve (usually modeled as a pressure sustaining
valve) to gradually fill the tank without seriously affecting pressure in the surrounding
system.
Figure 3.12
Relationship between
floating, pressurized,
and pumped tanks
HGL
Elevated Standpipe
Ground
Buried
Ground
Buried
Hydro-
Pneumatic
//=//=
/
/
=
/
/
=
/
/
=
/
/
=
Pumped Storage Floating on System
88 Assembling a Model Chapter 3
Hydropneumatic Tanks. In most tanks, the water surface elevation in the tank
equals the HGL in the tank. In the case of a pressure tank, however, the HGL is higher
than the tanks water surface. Pressure tanks, also called hydropneumatic tanks, are
partly full of compressed air. Because the water in the tank is pressurized, the HGL is
higher than the water surface elevation, as reflected in Equation 3.1.
HGL C
f
P Z + = (3.1)
where HGL = HGL of water in tank (ft, m)
P = pressure recorded at tank (psi, kPa)
Z = elevation of pressure gage (ft, m)
C
f
= unit conversion factor (2.31 English, 0.102 SI)
In steady-state models, a hydropneumatic tank can be represented by a tank or reser-
voir having this HGL. In EPS models, the tank must be represented by an equivalent
free-surface tank floating on the system. Because of the air in the tank, a hydropneu-
matic tank has an effective volume that is less than 30 to 50 percent of the total
volume of the tank. Modeling the tank involves first determining the minimum and
maximum pressures occurring in the tank and converting them to HGL values using
Equation 3.1. The cross-sectional area (or diameter) of this equivalent tank can be
determined by using Equation 3.2.
A
eq
V
eff
HGL
max
HGL
min

------------------------------------------------ = (3.2)
where A
eq
= area of equivalent tank (ft
2
, m
2
)
V
eff
= effective volume of tank (ft
3
, m
3
)
HGL
max
= maximum HGL in tank (ft, m)
HGL
min
= minimum HGL in tank (ft, m)
The relationship between the actual hydropneumatic tank and the model tank is
shown in Figure 3.13.
Using this technique, the EPS model of the tank will track HGL at the tank and vol-
ume of water in the tank, but not the actual water level.
3.5 JUNCTIONS
As the term implies, one of the primary uses of a junction node is to provide a loca-
tion for two or more pipes to meet. Junctions, however, do not need to be elemental
intersections, as a junction node may exist at the end of a single pipe (typically
referred to as a dead-end). The other chief role of a junction node is to provide a loca-
tion to withdraw water demanded from the system or inject inflows (sometimes
referred to as negative demands) into the system.
Junction nodes typically do not directly relate to real-world distribution components,
since pipes are usually joined with fittings, and flows are extracted from the system at
any number of customer connections along a pipe. From a modeling standpoint, the
Section 3.5 Junctions 89
importance of these distinctions varies, as discussed in the section on skeletonization
on page 112. Most water users have such a small individual impact that their with-
drawals can be assigned to nearby nodes without adversely affecting a model.
Figure 3.13
Relationship between
a hydropneumatic
tank and a model tank
V
eff
HGL
max
HGL
min
Pressure Tank
Equivalent
Model
Tank
V
eff
Pump off
Pump on
Pump
Junction Elevation
Generally, the only physical characteristic defined at a junction node is its elevation.
This attribute may seem simple to define, but there are some considerations that need
to be taken into account before assigning elevations to junction nodes. Because pres-
sure is determined by the difference between calculated hydraulic grade and eleva-
tion, the most important consideration is, at what elevation is the pressure most
important?
Selecting an Elevation. Figure 3.14 represents a typical junction node, illus-
trating that at least four possible choices for elevation exist that can be used in the
model. The elevation could be taken as point A, the centerline of the pipe. Alterna-
tively, the ground elevation above the pipe (point B), or the elevation of the hydrant
(point C), may be selected. As a final option, the ground elevation at the highest ser-
vice point, point D, could be used. Each of these possibilities has associated benefits,
so the determination of which elevation to use needs to be made on a case-by-case
basis. Regardless of which elevation is selected, it is good practice to be consistent
within a given model to avoid confusion.
90 Assembling a Model Chapter 3
Figure 3.14
Elevation choices for
a junction node
//=
//=
//=
//=
//=
//=
/
/
=
/
/
=
/
/
=
/
/
=
//=//= //=//= //=//=
Pipe - A (622)
Ground - B (630)
High Service - D (650)
Hydrant
Elevation - C (635)
Service
Line
The elevation of the centerline of the pipe may be useful for determining pressure for
leakage studies, or it may be appropriate when modeling above-ground piping sys-
tems (such as systems used in chemical processing). Ground elevations may be the
easiest data to obtain and will also overlay more easily onto mapping systems that use
ground elevations. They are frequently used for models of municipal water distribu-
tion systems. Both methods, however, have the potential to overlook poor service
pressures because the model could incorrectly indicate acceptable pressures for a cus-
tomer who is notably higher than the ground or pipe centerline. In such cases, it may
be more appropriate to select the elevation based on the highest service elevation
required.
In the process of model calibration (see Chapter 7 for more about calibration), accu-
rate node elevations are crucial. If the elevation chosen for the modeled junction is not
the same as the elevation associated with recorded field measurements, then direct
pressure comparisons are meaningless. Methods for obtaining good node elevation
data are described in Walski (1999).
3.6 PIPES
A pipe conveys flow as it moves from one junction node to another in a network. In
the real world, individual pipes are usually manufactured in lengths of around 18 or
20 feet (6 meters), which are then assembled in series as a pipeline. Real-world pipe-
lines may also have various fittings, such as elbows, to handle abrupt changes in
direction, or isolation valves to close off flow through a particular section of pipe.
Figure 3.15 shows ductile iron pipe sections.
For modeling purposes, individual segments of pipe and associated fittings can all be
combined into a single pipe element. A model pipe should have the same characteris-
tics (size, material, etc.) throughout its length.
Section 3.6 Pipes 91
Figure 3.15
Ductile iron pipe
sections
Length
The length assigned to a pipe should represent the full distance that water flows from
one node to the next, not necessarily the straight-line distance between the end nodes
of the pipe.
Scaled versus Schematic. Most simulation software enables the user to indi-
cate either a scaled length or a user-defined length for pipes. Scaled lengths are auto-
matically determined by the software, or scaled from the alignment along an
electronic background map. User-defined lengths, applied when scaled electronic
maps are not available, require the user to manually enter pipe lengths based on some
other measurement method, such as use of a map wheel (see Figure 3.16). A model
using user-defined lengths is a schematic model. The overall connectivity of a sche-
matic model should be identical to that of a scaled model, but the quality of the plani-
metric representation is more similar to a caricature than a photograph.
Even in some scaled models, there may be areas where there are simply too many
nodes in close proximity to work with them easily at the model scale (such as at a
pump station). In these cases, the modeler may want to selectively depict that portion
of the system schematically, as shown in Figure 3.17.
Diameter
As with junction elevations, determining a pipes diameter is not as straightforward as
it might seem. A pipes nominal diameter refers to its common name, such as a 16-in.
(400-mm) pipe. The pipes internal diameter, the distance from one inner wall of the
pipe to the opposite wall, may differ from the nominal diameter because of manufac-
turing standards. Most new pipes have internal diameters that are actually larger than
the nominal diameters, although the exact measurements depend on the class (pres-
sure rating) of pipe.
92 Assembling a Model Chapter 3
Figure 3.16
Use of a map
measuring wheel for
measuring pipe
lengths
Figure 3.17
Scaled system with a
schematic of a pump
station
Scaled System
Pump Station Schematic
(not to scale)
For example, Figure 3.18 depicts a new ductile iron pipe with a 16-in. nominal diam-
eter (ND) and a 250-psi pressure rating that has an outside diameter (OD) of 17.40 in.
and a wall thickness (Th) of 0.30 in., resulting in an internal diameter (ID) of 16.80 in.
(AWWA, 1996).
To add to the confusion, the ID may change over time as corrosion, tuberculation, and
scaling occur within the pipe (see Figure 3.19). Corrosion and tuberculation are
related in iron pipes. As corrosion reactions occur on the inner surface of the pipe, the
reaction by-products expand to form an uneven pattern of lumps (or tubercules) in a
process called tuberculation. Scaling is a chemical deposition process that forms a
material build-up along the pipe walls due to chemical conditions in the water. For
Section 3.6 Pipes 93
example, lime scaling is caused by the precipitation of calcium carbonate. Scaling can
actually be used to control corrosion, but when it occurs in an uncontrolled manner it
can significantly reduce the ID of the pipe.
Figure 3.18
Cross-section of a
16-in. pipe
Figure 3.19
Pipe corrosion and
tuberculation
Courtesy of Donald V. Chase, Department of Civil Engineering, University of Dayton
Of course, no one is going to refer to a pipe as a 16.80-in. (426.72-mm) pipe, and
because of the process just described, it is difficult to measure a pipes actual internal
diameter. As a result, a pipes nominal diameter is commonly used in modeling, in
combination with a roughness value that accounts for the diameter discrepancy. How-
ever, using nominal rather than actual diameters can cause significant differences
94 Assembling a Model Chapter 3
when water quality modeling is performed. Because flow velocity is related to flow
rate by the internal diameter of a pipe, the transport characteristics of a pipe are
affected. Chapter 7 discusses these calibration issues further (see page 255). Typical
roughness values can be found in Section 2.4.
Minor Losses
Including separate modeling elements to represent every fitting and appurtenance
present in a real-world system would be an unnecessarily tedious task. Instead, the
minor losses caused by those fittings are typically associated with pipes (that is, minor
losses are assigned as a pipe property).
In many hydraulic simulations, minor losses are ignored because they do not contrib-
ute substantially to the overall head loss throughout the system. In some cases,
however, flow velocities within a pipe and the configuration of fittings can cause
minor losses to be considerable (for example, at a pump station). The term minor is
relative, so the impact of these losses varies for different situations.
Red Water
Distribution systems with unlined iron or steel
pipes can be subject to water quality problems
related to corrosion, referred to as red water. Red
water is treated water containing a colloidal sus-
pension of very small, oxidized iron particles that
originated from the surface of the pipe wall. Over
a long period of time, this form of corrosion weak-
ens the pipe wall and leads to the formation of
tubercles. The most obvious and immediate
impact, however, is that the oxidized iron particles
give the water a murky, reddish-brown color. This
reduction in the aesthetic quality of the water
prompts numerous customer complaints.
Several alternative methods are available to con-
trol the pipe corrosion that causes red water. The
most traditional approach is to produce water that
is slightly supersaturated with calcium carbonate.
When the water enters the distribution system, the
dissolved calcium carbonate slowly precipitates
on the pipe walls, forming a thin, protective scale
(Caldwell and Lawrence, 1953; Merrill and Sanks,
1978). The Langelier Index (an index of the corro-
sive potential of water) can be used as an indica-
tion of the potential of the water to precipitate
calcium carbonate, allowing better management
of the precipitation rate (Langelier, 1936).
A positive saturation index indicates that the pipe
should be protected, provided that sufficient alka-
linity is present.
More recently, corrosion inhibitors such as zinc
orthophosphate and hexametaphosphate have
become popular in red water prevention (Ben-
jamin, Reiber, Ferguson, Vanderwerff, and Miller,
1990; Mullen and Ritter, 1974; Volk, Dundore,
Schiermann, and LeChevallier, 2000). Several
theories exist concerning the predominant mecha-
nism by which these inhibitors prevent corrosion.
The effectiveness of corrosion control measures
can be dependent on the hydraulic flow regime
occurring in the pipe. Several researchers have
reported that corrosion inhibitors and carbonate
films do not work well in pipes with low velocities
(Maddison and Gagnon, 1999; McNeil and
Edwards, 2000). Water distribution models pro-
vide a way to identify pipes with chronic low veloc-
ities, and therefore more potential for red water
problems. The effect of field operations meant to
control red water (for example, flushing and blow-
offs) can also be investigated using hydraulic
model simulations.
Section 3.7 Pumps 95
Composite Minor Losses. At any instant in time, velocity in the model is con-
stant throughout the length of a particular pipe. Since individual minor losses are
related to a coefficient multiplied by a velocity term, the overall head loss from sev-
eral minor losses is mathematically equivalent to having a single composite minor
loss coefficient. This composite coefficient is equal to the simple sum of the individual
coefficients.
3.7 PUMPS
A pump is an element that adds energy to the system in the form of an increased
hydraulic grade. Since water flows downhill (that is, from higher energy to lower
energy), pumps are used to boost the head at desired locations to overcome piping
head losses and physical elevation differences. Unless a system is entirely operated by
gravity, pumps are an integral part of the distribution system.
In water distribution systems, the most frequently used type of pump is the centrifugal
pump. A centrifugal pump has a motor that spins a piece within the pump called an
impeller. The mechanical energy of the rotating impeller is imparted to the water,
resulting in an increase in head. Figure 3.20 illustrates a cross-section of a centrifugal
pump and the flow path water takes through it. Water from the intake pipe enters the
pump through the eye of the spinning impeller (1) where it is then thrown outward
between vanes and into the discharge piping (2).
Figure 3.20
Cross-section of a
centrifugal pump
Frank M. White, Fluid Mechanics, 1994, McGraw-Hill, Inc. Reproduced by permission of the McGraw-Hill Companies.
Pump Characteristic Curves
With centrifugal pumps, pump performance is a function of flow rate. The perfor-
mance is described by the following four parameters, which are plotted versus
discharge.
Head: Total dynamic head added by pump in units of length (see page 44)
Efficiency: Overall pump efficiency (wire-to-water efficiency) in units of
percent (see pages 199 and 442)
96 Assembling a Model Chapter 3
Brake horsepower: Power needed to turn pump (in power units)
Net positive suction head (NPSH) required: Head above vacuum (in units
of length) required to prevent cavitation (see page 48)
Only the head curve is an energy equation necessary for solving pipe network prob-
lems. The other curves are used once the network has been solved to identify power
consumption (energy), motor requirements (brake horsepower), and suction piping
(NPSH).
Fixed-Speed and Variable-Speed Pumps. A pump characteristic curve is
related to the speed at which the pump motor is operating. With fixed-speed pumps,
the motor remains at a constant speed regardless of other factors. Variable-speed
pumps, on the other hand, have a motor or other device that can change the pump
speed in response to the system conditions.
A variable-speed pump is not really a special type of pump, but rather a pump con-
nected to a variable-speed drive or controller. The most common type of variable-
speed drive controls the flow of electricity to the pump motor, and therefore controls
the rate at which the pump rotates. The difference in pump speed, in turn, produces
different head and discharge characteristics. Variable-speed pumps are useful in appli-
cations requiring operational flexibility, such as when flow rates change rapidly, but
the desired pressure remains constant. An example of such a situation would be a net-
work with little or no storage available.
Power and Efficiency. The term power may have one of several meanings
when dealing with a pump. These possible meanings are listed below:
Input power: The amount of power that is delivered to the motor, usually in
electric form
Brake power: The amount of power that is delivered to the pump from the
motor
Water power: The amount of power that is delivered to the water from the
pump
Of course, there are losses as energy is converted from one form to another (electricity
to motor, motor to pump, pump to water), and every transfer has an efficiency associ-
ated with it. The efficiencies associated with these transfers may be expressed either
as percentages (100 percent is perfectly efficient) or as decimal values (1.00 is per-
fectly efficient), and are typically defined as follows:
Motor efficiency: The ratio of brake power to input power
Pump efficiency: The ratio of water power to brake power
Wire-to-water (overall) efficiency: The ratio of water power to input power
Pump efficiency tends to vary significantly with flow, while motor efficiency remains
relatively constant over the range of loads imposed by most pumps. Note that there
Section 3.7 Pumps 97
may also be an additional efficiency associated with a variable-speed drive. Some
engineers refer to the combination of the motor and any speed controls as the driver.
Figure 3.21 shows input power and wire-to-water efficiency curves overlaid on a typi-
cal pump head curve. Notice that the input power increases as discharge increases,
and head decreases as discharge increases. For each impeller size, there is a flow rate
corresponding to maximum efficiency. At higher or lower flows, the efficiency
decreases. This maximum point on the efficiency curve is called the best efficiency
point (BEP).
Figure 3.21
Pump curves with
efficiency, NPSH, and
horsepower overlays
Courtesy of Peerless Pumps
Obtaining Pump Data. Ideally, a water utility will have pump operating
curves on file for every pump in the system. These are usually furnished to the utility
with the shop drawings of the pump stations or as part of the manufacturers submit-
tals when replacing pumps. If the pump curve cannot be located, a copy of the curve
can usually be obtained from the manufacturer (provided the model and serial num-
bers for the pump are available).
To perform energy cost calculations, pump efficiency curves should also be obtained.
Note that the various power and efficiency definitions can be confusing, and it is
important to distinguish which terms are being referred to in any particular document.
Every pump differs slightly from its catalog model, and normal wear and tear will
cause a pump's performance to change over time. Thus, pumps should be checked to
verify that the characteristic curves on record are in agreement with field perfor-
mance. If an operating point does not agree with a characteristic curve, a new curve
can be developed to reflect the actual behavior. More information is available on this
subject in Chapter 5 (see page 199).
98 Assembling a Model Chapter 3
Even though a pump curve on record may not perfectly match the actual pump char-
acteristics, many utilities accept that the catalogued values for the pump curve are
sufficiently accurate for the purposes of the model, and forgo any performance testing
or field verification. This decision is dependent on the specific situation.
Model Representation
In order to model a pumps behavior, some mathematical expression describing its
pump head curve must be defined. Different models support different definitions, but
most are centered on the same basic concept, furnishing the model with sufficient
sample points to define the characteristic head curve.
Selecting Representative Points. As discussed previously, the relationship
between pump head and discharge is nonlinear. For most pumps, three points along
the curve are usually enough to represent the normal operating range of the pump.
These three points include
The zero-discharge point, also known as the cutoff or shutoff point
The normal operating point, which should typically be close to the best effi-
ciency point of the pump
The point at the maximum expected discharge value
Positive Displacement Pumps
Virtually all water distribution system pumps are
centrifugal pumps. However, pipe network models
are used in other applicationssuch as chemical
feeds, low-pressure sanitary sewer collection sys-
tems, and sludge pumpingin which positive dis-
placement pumps (for example, diaphragm,
piston, plunger, lobe, and progressive cavity
pumps) are used. Unlike centrifugal pumps, these
pumps produce a constant flow, regardless of the
head supplied, up to a very high pressure.
The standard approximations to pump curves
used in most models do not adequately address
positive displacement pumps because the head
characteristic curve for such pumps consists of a
virtually straight, vertical line. Depending on the
model, forcing a pump curve to fit this shape usu-
ally results in warning messages.
An easy way to approximate a positive displace-
ment pump in a model is to not include a pump at
all but rather to use two nodesa suction node
and a discharge nodethat are not connected.
The suction side node would have a demand set
equal to the pump flow, while the discharge node
would have an inflow set equal to this flow. The
model will then give the suction and discharge
HGLs and pressures at the nodes. (Custom
extended curve options can also be used.)
Because the suction and discharge systems are
separated, it is important for the modeler to
include a tank or reservoir on both the suction and
discharge sides of the pump. Otherwise, the
model will not be able to satisfy the law of conser-
vation of mass. For example, if the demands on
the discharge side do not equal the inflow to the
discharge side, the model may not give a valid
solution. Because most models assume demands
as independent of pressure, inflows must equal
system demands, plus or minus any storage
effects. If no storage is present, the model cannot
solve unless inflows and demands are equal.
Section 3.7 Pumps 99
It is also possible to provide some models with additional points along the pump
curve, but not all models treat these additional data points in the same way. Some
models perform linear interpolation between points, some fit a polynomial curve
between points, and others determine an overall polynomial or exponential curve that
fits the entire data set.
Constant Power Pumps. Many models also support the concept of a constant
power pump. With this type of pump, the water power produced by the pump remains
constant, regardless of how little or how much flow the pump passes.
Water power is a product of discharge and head, which means that a curve depicting
constant water power is asymptotic to both the discharge and head axes, as shown in
Figure 3.22.
Figure 3.22
Characteristic pump
curve for a constant
power pump
Discharge
H
e
a
d Actual
pump
curve
Equivalent
pump curve
Some modelers use a constant power pump definition to define a curve simply
because it is easier than providing several points from the characteristic curve, or
because the characteristic curve is not available. The results generated using this defi-
nition, however, can be unreliable and sometimes counter-intuitive. As shown in Fig-
ure 3.22, the constant power approximation will be accurate for a specific range of
flows, but not at very high or low flows. For very preliminary studies when all the
modeler knows is the approximate size of the pump, this approximation can be used
to get into pipe sizing quickly. However, it should not be used for pump selection.
The modeler must remember that the power entered for the constant power pump is
not the rated power of the motor but the water power added. For example, a 50 hp
motor that is 90 percent efficient, running at 80 percent of its rated power, and con-
nected to a pump that is operating at 70 percent efficiency will result in a water power
of roughly 25 hp (that is, 50 0.9 u 0.8 u 0.7 u ). The value 25 hp, not 50 hp, should be
entered into the model.
Node versus Link Representation. A pump can be represented as a node or
a link element, depending on the software package. In software that symbolizes
pumps as links, the pump connects upstream and downstream nodes in a system the
same way a pipe would. A link symbolization more closely reflects the internal math-
ematical representation of the pump, but it can introduce inaccuracies. For example,
Figure 3.23 illustrates how the pump intake and discharge piping may be ignored and
the head losses occurring in them neglected.
100 Assembling a Model Chapter 3
Figure 3.23
Comparison of an
actual pump and a
pump modeled as a
link element
Other models represent pumps as nodes, typically with special connectivity rules (for
example, only allowing a single downstream pipe). This nodal representation is less
error-prone, more realistic, and easier for the modeler to implement. Nodal represen-
tation may also be more intuitive, since a real-world pump is usually thought of as
being in a single location with two distinct hydraulic grades (one on the intake side
and one on the discharge side). Figure 3.24 illustrates a nodal representation of a
pump.
Figure 3.24
Comparison of an
actual pump and a
pump modeled as a
node element
Real World Model (Node)
HGL HGL
3.8 VALVES
A valve is an element that can be opened and closed to different extents (called throt-
tling) to vary its resistance to flow, thereby controlling the movement of water
through a pipeline (see Figure 3.25). Valves can be classified into the following five
general categories:
Isolation valves
Directional valves
Altitude valves
Air release and vacuum breaking valves
Control valves
Section 3.8 Valves 101
Figure 3.25
Different valve types
Gate Valve Butterfly Valve Check Valve
Courtesy of Crane Co. All Rights Reserved.
Some valves are intended to automatically restrict the flow of water based on pres-
sures or flows, and others are operated manually and used to completely turn off por-
tions of the system. The behaviors of different valve types vary significantly
depending on the software used. This section provides an introduction to some of the
most common valve types and applications.
Isolation Valves
Perhaps the most common type of valve in water distribution systems is the isolation
valve, which can be manually closed to block the flow of water. As the term isola-
tion implies, the primary purpose of these valves is to provide a field crew with a
means of turning off a portion of the system to, for example, replace a broken pipe or
a leaky joint. Well-designed water distribution systems have isolation valves through-
out the network, so that maintenance and emergencies affect as few customers as pos-
sible. In some systems, isolation valves may be intentionally kept in a closed position
to control pressure zone boundaries, for example.
There are several types of isolation valves that may be used, including gate valves (the
most popular type), butterfly valves, globe valves, and plug valves.
In most hydraulic models, the inclusion of each and every isolation valve would be an
unnecessary level of detail. Instead, the intended behavior of the isolation valve
(minor loss, the ability to open and close, and so on) can be defined as part of a pipe.
A common question in constructing a model is whether to explicitly include minor
losses due to open gate valves, or to account for the effect of such losses in the Hazen-
Williams C-factor. If the C-factor for the pipe with no minor losses is known, an
equivalent C-factor that accounts for the minor losses is given by:
C
e
C
L
L D
6K
L
f
----------


+
------------------------------





0.54
= (3.3)
102 Assembling a Model Chapter 3
where C
e
= equivalent Hazen-Williams C-factor accounting for minor losses
C = Hazen-Williams C-factor
L = length of pipe segment (ft, m)
D = diameter (ft, m)
f = Darcy-Weisbach friction factor
6 K
L
= sum of minor loss coefficients in pipe
For example, consider a 400-ft (122-m) segment of 6-in. (152-mm) pipe with a C-
factor of 120 and an f of 0.02. From Equation 3.3, the equivalent C-factor for the pipe
including a single open gate valve (K
L
= 0.39) is 118.4. For two open gate valves, the
equivalent C-factor is 116.9. Given that C-factors are seldom known to within plus or
minus 5, these differences are generally negligible. Note that if a model is calibrated
without explicitly accounting for many minor losses, then the C-factor resulting from
the calibration is the equivalent C-factor, and no further adjustment is needed.
Directional Valves
Directional valves, also called check valves, are used to ensure that water can flow in
one direction through the pipeline, but cannot flow in the opposite direction (back-
flow). Any water flowing backwards through the valve causes it to close, and it
remains closed until the flow once again begins to go through the valve in the forward
direction.
Simple check valves commonly use a hinged disk or flap to prevent flow from travel-
ing in the undesired direction. For example, the discharge piping from a pump may
include a check valve to prevent flow from passing through the pump backwards
(which could damage the pump). Most models automatically assume that every pump
has a built-in check valve, so there is no need to explicitly include one (see Figure
3.26). If a pump does not have a check valve on its discharge side, water can flow
backwards through the pump when the power is off. This situation can be modeled
with a pipe parallel to the pump that only opens when the pump is off. The pipe must
have an equivalent length and minor loss coefficient that will generate the same head
loss as the pump running backwards.
Figure 3.26
A check valve
operating at a pump
Pump Off
Demand
Check Valve
HGL
HGL
Demand
Check Valve
Pump On
Section 3.8 Valves 103
Mechanically, some check valves require a certain differential in head before they will
seat fully and seal off any backflow. They may allow small amounts of reverse flow,
which may or may not have noteworthy consequences. When potable water systems
are hydraulically connected to nonpotable water uses, a reversal of flow could be
disastrous. These situations, called cross-connections, are a serious danger for water
distributors, and the possibility of such occurrences warrants the use of higher quality
check valves. Figure 3.27 illustrates a seemingly harmless situation that is a potential
cross-connection. A device called a backflow preventer is designed to be highly sensi-
tive to flow reversal, and frequently incorporates one or more check valves in series to
prevent backflow.
Figure 3.27
A potential cross-
connection
As far as most modeling software is concerned, there is no difference in sensitivity
between different types of check valves (all are assumed to close completely even for
the smallest of attempted reverse flows). As long as the check valve can be repre-
sented using a minor loss coefficient, the majority of software packages allow them to
be modeled as an attribute associated with a pipe, instead of requiring that a separate
valve element be created.
Altitude Valves
Many water utilities employ devices called altitude valves at the point where a pipe-
line enters a tank (see Figure 3.28). When the tank level rises to a specified upper
limit, the valve closes to prevent any further flow from entering, thus eliminating
overflow. When the flow trend reverses, the valve reopens and allows the tank to drain
to supply the usage demands of the system.
Most software packages, in one form or another, automatically incorporate the behav-
ior of altitude valves at both the minimum and maximum tank levels and do not
require explicit inclusion of them. If, however, an altitude valve does not exist at a
tank, tank overflow is possible, and steps must be taken to include this behavior in the
model.
104 Assembling a Model Chapter 3
Figure 3.28
Altitude valve
controlling the
maximum fill level of
a tank
Filling/Draining
Max. Fill Level
Altitude Valve Open
(Upper Limit)
No Flow
//=//= //=//=
Vault
Altitude Valve Closed
Vault
//=//= //=//= //=//= //=//= //=//= //=//=
Air Release Valves and Vacuum Breaking Valves
Most systems include special air release valves to release trapped air during system
operation, and air/vacuum valves that discharge air upon system start-up and admit
air into the system in response to negative gage pressures (see Figure 3.29). These
types of valves are often found at system high points, where trapped air settles, and at
changes in grade, where pressures are most likely to drop below ambient or atmo-
spheric conditions. Combination air valves that perform the functions of both valve
types are often used as well.
Air release and air/vacuum valves are typically not included in standard water distri-
bution system modeling. The importance of such elements is significant, however, for
advanced studies such as transient analyses.
Figure 3.29
Air release and air/
vacuum valves
Air Release Valve Vacuum Breaking Valve
Courtesy of Val-Matic Valve and Manufacturing Corporation, Elmhurst, Illinois.
Control Valves
For any control valve, also called regulating valve, the setting is of primary impor-
tance. For a flow control valve, this setting refers to the flow setting, and for a throttle
Section 3.8 Valves 105
control valve, it refers to a minor loss coefficient. For pressure-based controls, how-
ever, the setting may be either the hydraulic grade or the pressure that the valve tries
to maintain. Models are driven by hydraulic grade, so if a pressure setting is used, it is
critically important to have not only the correct pressure setting, but also the correct
valve elevation.
Given the setting for the valve, the model calculates the flow through the valve and
the inlet and outlet HGL (and pressures). A control valve is complicated in that,
unlike a pump, which is either on or off, it can be in any one of the several states
described in the following list. Note that the terminology may vary slightly between
models.
Active: Automatically controlling flow
- Open: Opened fully
- Closed (1): Closed fully
- Throttling: Throttling flow and pressure
Closed (2): Manually shut, as when an isolating valve located at the control
valve is closed
Inactive: ignored
Because of the many possible control valve states, valves are often points where
model convergence problems exist.
Pressure Reducing Valves (PRVs). Pressure reducing valves (PRVs) throttle
automatically to prevent the downstream hydraulic grade from exceeding a set value,
and are used in situations where high downstream pressures could cause damage. For
example, Figure 3.30 illustrates a connection between pressure zones. Without a PRV,
the hydraulic grade in the upper zone could cause pressures in the lower zone to be
high enough to burst pipes or cause relief valves to open.
Figure 3.30
Schematic network
illustrating the use of
a pressure reducing
valve
PRV
Target Maximum Grade
Tank
on
Hill
Higher
Service
Area
Lower
Service
Area
Without PRV
With PRV
HGL
106 Assembling a Model Chapter 3
Unlike the isolation valves discussed earlier, PRVs are not associated with a pipe but
are explicitly represented within a hydraulic model. A PRV is characterized in a
model by the downstream hydraulic grade that it attempts to maintain, its controlling
status, and its minor loss coefficient. Because the valve intentionally introduces losses
to meet the required grade, a PRV's minor loss coefficient is really only a concern
when the valve is wide open (not throttling).
Like pumps, PRVs connect two pressure zones and have two associated hydraulic
grades, so some models represent them as links and some represent them as nodes.
The pitfalls of link characterization of PRVs are the same as those described previ-
ously for pumps (see page 99).
Pressure Sustaining Valves (PSVs). A pressure sustaining valve (PSV)
throttles the flow automatically to prevent the upstream hydraulic grade from drop-
ping below a set value. This type of valve can be used in situations in which unregu-
lated flow would result in inadequate pressures for the upstream portion of the system
(see Figure 3.31). They are frequently used to model pressure relief valves (see page
313).
Like PRVs, a PSV is typically represented explicitly within a hydraulic model and is
characterized by the upstream pressure it tries to maintain, its status, and its minor
loss coefficient.
Figure 3.31
Schematic network
illustrating the use of
a pressure sustaining
valve
PSV
Target Minimum Grade
Without PSV
With PSV
Tank
on
Hill
Higher
Service
Area
Lower
Service
Area
HGL
HGL
Flow Control Valves (FCVs). Flow control valves (FCVs) automatically throt-
tle to limit the rate of flow passing through the valve to a user-specified value. This
type of valve can be employed anywhere that flow-based regulation is appropriate,
such as when a water distributor has an agreement with a customer regarding maxi-
mum usage rates. FCVs do not guarantee that the flow will not be less than the setting
value, only that the flow will not exceed the setting value. If the flow does not equal
the setting, modeling packages will typically indicate so with a warning.
Similar to PRVs and PSVs, most models directly support FCVs, which are character-
ized by their maximum flow setting, status, and minor loss coefficient.
Section 3.9 Controls (Switches) 107
Throttle Control Valves (TCVs). Unlike an FCV where the flow is specified
directly, a throttle control valve (TCV) throttles to adjust its minor loss coefficient
based on the value of some other attribute of the system (such as the pressure at a crit-
ical node or a tank water level). Often the throttling effect of a particular valve posi-
tion is known, but the minor loss coefficients as a function of position are unknown.
This relationship can frequently be provided by the manufacturer.
Valve Books
Many water utilities maintain valve books, which are sets of records that provide
details pertaining to the location, type, and status of isolation valves and other fittings
throughout a system. From a modeling perspective, valve books can provide valuable
insight into the pipe connectivity at hydraulically complex intersections, especially in
areas where system maps may not show all of the details.
3.9 CONTROLS (SWITCHES)
Operational controls, such as pressure switches, are used to automatically change the
status or setting of an element based on the time of day, or in response to conditions
within the network. For example, a switch may be set to turn on a pump when pres-
sures within the system drop below a desired value. Or a pump may be programmed
to turn on and refill a tank in the early hours of the morning.
Without operational controls, conditions would have to be monitored and controlled
manually. This type of operation would be expensive, mistake-prone, and sometimes
impractical. Automated controls enable operators to take a more supervisory role,
focusing on issues larger than the everyday process of turning on a pump at a given
time or changing a control valve setting to accommodate changes in demand. Conse-
quently, the system can be run more affordably, predictably, and practically.
Models can represent controls in different ways. Some consider controls to be sepa-
rate modeling elements, and others consider them to be an attribute of the pipe, pump,
or valve being controlled.
Pipe Controls
For a pipe, the only status that can really change is whether the pipe (or, more accu-
rately, an isolation valve associated with the pipe) is open or closed. Most pipes will
always be open, but some pipes may be opened or closed to model a valve that auto-
matically or manually changes based on the state of the system. If a valve in the pipe
is being throttled, it should be handled either through the use of a minor loss directly
applied to the pipe or by inserting a throttle control valve in the pipe and adjusting it.
Pump Controls
The simplest type of pump control turns a pump on or off. For variable-speed pumps,
controls can also be used to adjust the pumps relative speed factor to raise or lower
108 Assembling a Model Chapter 3
the pressures and flow rates that it delivers. For more information about pump relative
speed factors, see Chapter 2 (page 44).
The most common way to control a pump is by tank water level. Pumps are classified
as either lead pumps, which are the first to turn on, or lag pumps, the second to
turn on. Lead pumps are set to activate when tanks drain to a specified minimum level
and to shut off when tanks refill to a specified maximum level, usually just below the
tank overflow point. Lag pumps turn on only when the tank continues to drain below
the minimum level, even with the lead pump still running. They turn off when the tank
fills to a point below the shut off level for the lead pump. Controls get much more
complicated when there are other considerations such as time of day control rules or
parallel pumps that are not identical.
Regulating Valve Controls
Similar to a pump, a control valve can change both its status (open, closed, or active)
and its setting. For example, an operator may want a flow control valve to restrict flow
more when upstream pressures are poor, or a pressure reducing valve to open com-
pletely to accommodate high flow demands during a fire event.
Indicators of Control Settings
If a pressure switch setting is unknown, tank level charts and pumping logs may pro-
vide a clue. As shown in Figure 3.32, pressure switch settings can be determined by
looking at tank level charts and correlating them to the times when pumps are placed
into or taken out of service. Operations staff can also be helpful in the process of
determining pressure switch settings.
Figure 3.32
Correlation between
tank levels and pump
operation
Time
T
a
n
k
L
e
v
e
l
P
u
m
p
F
l
o
w
R
a
t
e
Time
Section 3.10 Types of Simulations 109
3.10 TYPES OF SIMULATIONS
After the basic elements and the network topology are defined, further refinement of
the model can be done depending on its intended purpose. There are various types of
simulations that a model may perform, depending on what the modeler is trying to
observe or predict. The two most basic types are
Steady-state simulation: Computes the state of the system (flows, pres-
sures, pump operating attributes, valve position, and so on) assuming that
hydraulic demands and boundary conditions do not change with respect to
time.
Extended-period simulation (EPS): Determines the quasi-dynamic behav-
ior of a system over a period of time, computing the state of the system as a
series of steady-state simulations in which hydraulic demands and boundary
conditions do change with respect to time.
Steady-State Simulation
As the term implies, steady-state refers to a state of a system that is unchanging in
time, essentially the long-term behavior of a system that has achieved equilibrium.
Tank and reservoir levels, hydraulic demands, and pump and valve operation remain
constant and define the boundary conditions of the simulation. A steady-state simula-
tion provides information regarding the equilibrium flows, pressures, and other vari-
ables defining the state of the network for a unique set of hydraulic demands and
boundary conditions.
Real water distribution systems are seldom in a true steady state. Therefore, the
notion of a steady state is a mathematical construct. Demands and tank water levels
are continuously changing, and pumps are routinely cycling on and off. A steady-state
hydraulic model is more like a blurred photograph of a moving object than a sharp
photo of a still one. However, by enabling designers to predict the response to a
unique set of hydraulic conditions (for example, peak hour demands or a fire at a par-
ticular node), the mathematical construct of a steady state can be a very useful tool.
Steady-state simulations are the building blocks for other types of simulations. Once
the steady-state concept is mastered, it is easier to understand more advanced topics
such as extended-period simulation, water quality analysis, and fire protection studies
(these topics are discussed in later chapters).
Steady-state models are generally used to analyze specific worst-case conditions such
as peak demand times, fire protection usage, and system component failures in which
the effects of time are not particularly significant.
Extended-Period Simulation
The results provided by a steady-state analysis can be extremely useful for a wide
range of applications in hydraulic modeling. There are many cases, however, for
which assumptions of a steady-state simulation are not valid, or a simulation is
required that allows the system to change over time. For example, to understand the
110 Assembling a Model Chapter 3
effects of changing water usage over time, fill and drain cycles of tanks, or the
response of pumps and valves to system changes, an extended-period simulation
(EPS) is needed.
It is important to note that there are many inputs required for an extended-period sim-
ulation. Due to the volume of data and the number of possible actions that a modeler
can take during calibration, analysis, and design, it is highly recommended that a
model be examined under steady-state situations prior to working with extended-
period simulations. Once satisfactory steady-state performance is achieved, it is much
easier to proceed into EPSs.
EPS Calculation Process. Similar to the way a film projector flashes a series of
still images in sequence to create a moving picture, the hydraulic time steps of an
extended-period simulation are actually steady-state simulations that are strung
together in sequence. After each steady-state step, the system boundary conditions are
reevaluated and updated to reflect changes in junction demands, tank levels, pump
operations, and so on. Then, another hydraulic time step is taken, and the process con-
tinues until the end of the simulation.
Simulation Duration. An extended-period simulation can be run for any length
of time, depending on the purpose of the analysis. The most common simulation dura-
tion is typically a multiple of 24 hours, because the most recognizable pattern for
demands and operations is a daily one. When modeling emergencies or disruptions
that occur over the short-term, however, it may be desirable to model only a few hours
into the future to predict immediate changes in tank level and system pressures. For
water quality applications, it may be more appropriate to model a duration of several
days in order for quality levels to stabilize.
Even with established daily patterns, a modeler may want to look at a simulation
duration of a week or more. For example, consider a storage tank with inadequate
capacity operating within a system. The water level in the tank may be only slightly
less at the end of each day than it was at the end of the previous day, which may go
unnoticed when reviewing model results. If a duration of one or two weeks is used,
the trend of the tank level dropping more and more each day will be more evident.
Even in systems that have adequate storage capacity, a simulation duration of 48
hours or longer can be helpful in better determining the tank draining and filling char-
acteristics.
Hydraulic Time Step. An important decision when running an extended-period
simulation is the selection of the hydraulic time step. The time step is the length of
time for one steady-state portion of an EPS, and it should be selected such that
changes in system hydraulics from one increment to the next are gradual. A time step
that is too large may cause abrupt hydraulic changes to occur, making it difficult for
the model to give good results.
For any given system, predicting how small the time increment should be is difficult,
although experience is certainly beneficial in this area. Typically, modelers begin by
assuming one-hour time steps, unless there are considerations that point to the need
for a different time step.
Section 3.10 Types of Simulations 111
When junction demands and tank inflow/outflow rates are highly variable, decreasing
the time step can improve the accuracy of the simulation. The sensitivity of a model to
time increment changes can be explored by comparing the results of the same analysis
using different increments. This sensitivity can also be evaluated during the calibra-
tion process. Ultimately, finding the correct balance between calculation time and
accuracy is up to the modeler.
Intermediate Changes. Of course, changes within a system dont always
occur at even time increments. When it is determined that an elements status changes
between time steps (such as a tank completely filling or draining, or a control condi-
tion being triggered), many models will automatically report a status change and
results at that intermediate point in time. The model then steps ahead in time to the
next even increment until another intermediate time step is required. If calculations
are frequently required at intermediate times, the modeler should consider decreasing
the time increment.
Why Use a Scenario Manager?
When water distribution models were first created,
data were input into the computer program by
using punch cards, which were submitted and
processed as a batch run. In this type of run, a
separate set of input data was required to gener-
ate each set of results. Because a typical model-
ing project requires analysis of many alternative
situations, large amounts of time were spent cre-
ating and debugging multiple sets of input cards.
When data files replaced punch cards, the batch
approach to data entry was carried over. The
modeler could now edit and copy input files more
easily, but there was still the problem of trying to
manage a large number of model runs. Working
with many data files or a single data file with doz-
ens of edits was confusing, inefficient, and error-
prone.
The solution to this problem is to keep alternative
data sets within a single model data file. For
example, data for current average day demands,
maximum day demands with a fire flow at node
37, and peak hour demands in 2020 can be cre-
ated, managed, and stored in a central database.
Once this structure is in place, the user can then
create many runs, or scenarios, by piecing
together alternative data sets.
For example, a scenario may consist of the peak
hour demands in 2020 paired with infrastructure
data that includes a proposed tank on Washington
Hill and a new 16-in. (400 mm) pipe along North
Street. This idea of building model runs from alter-
native data sets created by the user is more intui-
tive than the batch run concept, and is consistent
with the object-oriented paradigm found in mod-
ern programs. Further, descriptive naming of sce-
narios and alternative data sets provides internal
documentation of the users actions.
Because alternative plans in water modeling tend
to grow out of previous alternatives, a good sce-
nario manager will use the concept of inheritance
to create new child alternatives from existing par-
ent alternatives. Combining this idea of inherit-
ance with construction of scenarios from
alternative data sets gives the model user a self-
documenting way to quickly create new and better
solutions based on the results of previous model
runs.
A user accustomed to performing batch runs may
find some of the terminology and concepts
employed in scenario management a bit of a chal-
lenge at first. But, with a little practice, it becomes
difficult to imagine building or maintaining a model
without this versatile feature.
112 Assembling a Model Chapter 3
Other Types of Simulations
Using the fundamental concepts of steady-state and extended-period simulations,
more advanced simulations can be built. Water quality simulations are used to ascer-
tain chemical or biological constituent levels within a system or to determine the age
or source of water (see page 61). Automated fire flow analyses establish the suitability
of a system for fire protection needs. Cost analyses are used for looking at the mone-
tary impact of operations and improvements. Transient analyses are used to investi-
gate the short-term fluctuations in flow and pressure due to sudden changes in the
status of pumps or valves (see page 573).
With every advance in computer technology and each improvement in software meth-
ods, hydraulic models become a more integral part of designing and operating safe
and reliable water distribution systems.
3.11 SKELETONIZATION
Skeletonization is the process of selecting for inclusion in the model only the parts of
the hydraulic network that have a significant impact on the behavior of the system.
Attempting to include each individual service connection, gate valve, and every other
component of a large system in a model could be a huge undertaking without a signif-
icant impact on the model results. Capturing every feature of a system would also
result in tremendous amounts of data; enough to make managing, using, and trouble-
shooting the model an overwhelming and error-prone task. Skeletonization is a more
practical approach to modeling that allows the modeler to produce reliable, accurate
results without investing unnecessary time and money.
Eggener and Polkowski (1976) did the first study of skeletonization when they sys-
tematically removed pipes from a model of Menomonie, Wisconsin, to test the sensi-
tivity of model results. They found that under normal demands, they could remove a
large number of pipes and still not affect pressure significantly. Shamir and Hamberg
(1988a, 1988b) investigated rigorous rules for reducing the size of models.
Skeletonization should not be confused with the omission of data. The portions of the
system that are not modeled during the skeletonization process are not discarded; rather,
their effects are accounted for within parts of the system that are included in the model.
Skeletonization Example
Consider the following proposed subdivision, which is tied into an existing water sys-
tem model. Figures 3.33, 3.34, 3.35, and 3.36 show how demands can be aggregated
from individual customers to nodes with larger and larger nodal service areas.
Although a modeler would almost never include the individual connections as shown
in Figure 3.33, this example, which can be extrapolated to much larger networks,
shows the steps that are followed to achieve various levels of skeletonization.
As depicted in the network segment in Figure 3.33, it is possible to not skeletonize at
all. In this case, there is a junction at each service tap, with a pipe and junction at each
house. There are also junctions at the main intersections, resulting in a total of nearly
50 junctions (not including those required for fire hydrants).
Section 3.11 Skeletonization 113
Figure 3.33
An all-link network
0.3 0.2 0.4 0.3 0.2 0.2 0.3
0.4
0.4
0.3
0.3
0.6
0.3 0.3 0.2 0.2 0.4
0.3
0.3
0.3 0.3 0.2 0.2 0.3
0.3 0.2 0.3 0.3 0.4
Demand in gpm
The same subdivision could be modeled again, but slightly more skeletonized. Instead
of explicitly including each household, only the tie-ins and main intersections are
included. This level of detail results in a junction count of less than 20 (Figure 3.34).
Note that in this level of skeletonization, hydraulic results for the customer service
lines would not be available since they were not included in the model. If results for
service lines are not important, then the skeletal model shown in Figure 3.34 repre-
sents an adequate level of detail.
Figure 3.34
Minimal
skeletonization
The system can be skeletonized even more, modeling only the ends of the main piping
and the major intersections (Figure 3.35). Attributing the demands to the junctions
becomes a little trickier since a junction is not being modeled at each tap location.
The demands for this model are attributed to the junction nearest to the service (fol-
lowing the pipeline). The dashed boundary areas indicate the contributing area for
each model junction. For example, the junction in the upper right will be assigned the
demand for eight houses, while the lower right junction has demands for ten houses,
and so on.
114 Assembling a Model Chapter 3
Figure 3.35
Moderate
skeletonization
3.5 gpm
2.5 gpm
2.7 gpm
An even greater level of skeletonization can be achieved using just a single junction
node where the subdivision feeds from the existing system. The piping within the
entire subdivision has been removed, with all demands being attributed to the remain-
ing junction (see Figure 3.36). In this case, the model will indicate the impact of the
demands associated with the subdivision on the overall hydraulic network. However,
the modeler will not be able to determine how pressures and flows vary within the
subdivision.
Figure 3.36
Maximum
skeletonization
8.7 gpm
An even broader level of skeletonization is possible in which even the junction node
where the subdivision piping ties into the main line is excluded. The subdivision
demands would simply be added to a nearby junction, where other effects may be
combined with those from several other subdivisions that also have not been included
in detail. As this example demonstrates, the extent of skeletonization depends on the
intended use of the model and, to a large degree, is subject to the modelers discretion.
Section 3.11 Skeletonization 115
Skeletonization Guidelines
There are no absolute criteria for determining whether a pipe should be included in
the model, but it is safe to say that all models are most likely skeletonized to some
degree. Water distribution networks vary drastically from one system to another, and
modeling judgment plays a large role in the creation of a solution. For a small-
diameter system, such as household plumbing or a fire sprinkler system, small
differences in estimated flow rate may have perceptible effects on the system head
losses. For a large city system, however, the effects of water demanded by an entire
subdivision may be insignificant for the large-transmission main system.
Opposing Philosophies. There are definitely opposing philosophies regarding
skeletonization that stem from different modeling perspectives. Some modelers assert
that a model should never be bigger than a few hundred elements, because no one can
possibly digest all of the data that pours out of a larger model. Others contend that a
model should include all the pipes, so that data-entry can be done by less skilled per-
sonnel, who will not need to exercise judgment about whether or not an element
should be included. Followers of this approach then use database queries, automated
consolidation algorithms, and demand allocation procedures (see page 136) to gener-
ate skeletonized models for individual applications.
Somewhere in the Middle. Most network models, however, fall somewhere
between the two extremes. The level of skeletonization used depends on the intended
use of the model. At one extreme, energy operation studies require minimal detail,
while determining available fire flow at individual hydrants requires the most. For
master planning or regional water studies, a broader level of skeletonization will typi-
116 Assembling a Model Chapter 3
cally suffice. For detailed design work or water quality studies, however, much more
of the system needs to be included to accurately model the real-world system.
The responsibility really comes back to the modeler, who must have a good under-
standing of the models intended use and must select a level of detail appropriate for
that purpose. Most modelers choose to develop their own skeletonization guidelines.
Elements of High Importance
Any elements that are important to the system or can potentially influence system
behavior should be included in the model. For most models this criterion includes
Large water consumers
Points of known conditions, such as sampling points
Critical points with unknown conditions
Large-diameter pipes
Pipes that complete important loops
Pumps, control valves, tanks, and other controlling elements
Elements of Unknown Importance
If the modeler is unsure what the effects of including or excluding specific elements
may be, there is a very simple method that can be used to find out exactly what the
effects are on the system. Run the model and see what happens.
A base skeleton can be created using experience and judgment, with pipes of ques-
tionable importance included. The model should be run over a range of study condi-
tions and the results noted. One or more questionable pipes can then be closed
(preventing them from conveying water) and the model run again. If the modeler
determines that the results from the two analyses are essentially the same, then the
pipes apparently did not have a significant effect on the system and can be removed
from the skeleton.
If a pipes level of significance cannot be determined or is questionable, it is usually
better to leave the pipe in the model. With older, nongraphical interfaces, it was often
desirable to limit the number of pipes as much as possible to prevent becoming lost in
the data. With the advanced computers and easy-to-use software tools of today, how-
ever, there are fewer reasons to exclude pipes from the model.
Automated Skeletonization
An increasing number of water utilities are linking their models to GIS systems and
even creating models from scratch by importing data from their GIS. However, there
are generally far more GIS elements than the user would want pipes in the model pre-
senting an obstacle for a smooth data conversion process. For example, Figure 3.37
shows how a single pipe link from a model can correspond to a large number of GIS
Section 3.11 Skeletonization 117
elements. The number of pipes in the GIS is even greater when each hydrant lateral
and service line is included in the GIS. Of course, the modeler can manually eliminate
pipes from the model, but this task can be extremely tedious and error prone, espe-
cially if it must be repeated for several time periods or planning scenarios. Thus, it is
highly desirable and clearly more efficient to automate the process of model skeleton-
ization.
Figure 3.37
GIS pipes versus
model pipes
Simply removing pipes and nodes from a model based on a rule, such as pipe size, is a
straightforward process. The process becomes complicated, however, when it is nec-
essary to also keep track of the demands (and associated demand patterns) and emitter
coefficients that were assigned to the nodes being removed, and it becomes even more
complicated when one tries to account for the hydraulic capacity of the pipes being
removed.
Skeletonization is not a single process but several different low-level element removal
processes that must be applied in series to ensure that the demands are logically
brought back to their source of supply. The skeletonization process also involves
developing rules for pumps, tanks, and valves, and deciding which pipes and nodes
should be identified as nonremovable.
As with manual skeletonization, the degree to which a system is skeletonized depends
on the type of raw data and the ultimate purpose of the model. If the raw data are a
complete GIS of the system including service lines and hydrant laterals and the model
is going to be used to set up pump controls or study energy costs, it may be possible to
remove the overwhelming majority of the pipes. On the other hand, if the model was
built manually from distribution maps and the model is to be used to determine avail-
able fire flow at every hydrant, then there may be little room for skeletonization.
118 Assembling a Model Chapter 3
The individual processes involved with skeletonization are discussed in the subsec-
tions that follow.
Simple Pipe Removal. The simplest type of pipe removal is when pipes are
simply removed from the system based on size or other criteria without any consider-
ation of their effects on demand loading or hydraulic capacity. This can be useful
when importing data from a GIS if the dataset contains service lines and hydrant later-
als. This type of pipe removal is usually practiced before demands are assigned to
model nodes (as a preprocessing step), although that is not always the case. Some
models that claim to perform automated skeletonization only perform this type of
skeletonization process.
Removing Branch Pipes. The next simplest type of skeletonization consists of
removing dead-end branches that do not contain tanks at the end. This process is
referred to as branch trimming, or branch collapsing, and the user needs to determine
whether some finite number of branches should be trimmed or if the network should
be trimmed back to a pipe that is part of a loop. Figure 3.38 shows how a branch is
trimmed back to a node that is part of a loop. When dead-end branch pipes are
removed, the removal has no effect on the carrying capacity of the remainder of the
system.
Figure 3.38
Branch pipe removal
(branch trimming) J-12
Q=4
J-11
Q=5
J-13
Q=8
J-10
Q=10
Before Trimming
After Trimming
J-10
Q=27
Removing Pipes in Series (with no other pipes connected to the
common node). In most cases, removing pipes in series (sometimes called pipe
merging) has a negligible effect on model performance. For example, in Figure 3.39,
pipes P-121 and P-122 can be combined to form a new P-121. In this example, the
Section 3.11 Skeletonization 119
demand (Q) at J-12 is split evenly between the two nodes at the ends of the resulting
pipe. Depending on the situation, however, other rules regarding demands can be
applied. For example, either of the nodes could receive all the demand, or the demand
could be split according to user-specified rules.
Figure 3.39
Series pipe removal
similar attributes
Before Series Pipe Removal
After Series Pipe Removal
J-11
Q=5
J-12
Q=8
J-13
Q=5
P-121
L=350
D=8
C=120
P-122
L=250
D=8
C=120
J-11
Q=9
P-121
L=600
D=8
C=120
J-13
Q=9
If the node between two pipes in series has a large demand, removing it may
adversely impact the model results. To prevent such situations, the modeler may con-
sider setting a limit on flows such that nodes that exceed the limit cannot be elimi-
nated.
A key issue in combining two pipes into one lies in determining the attributes of the
resulting pipe. In Figure 3.39, the length of the resulting pipe is equal to the sum of
the lengths of the two pipes being combined and because the two pipes have the same
diameter and C-factor, the resulting pipe also has the same diameter and C-factor.
The problem becomes more complicated when the two pipes have different attributes,
as shown in Figure 3.40. In this case, the length is still the sum of the length of the two
pipes, but now there are an infinite number of combinations of diameter and C-factor
that would produce the same head loss through the pipe. As an option, the modeler can
choose to use the diameter and C-factor of one of the pipes as the attributes for the
resulting pipe. Or, the modeler can pick either the C-factor or the diameter for the result-
ing pipe and then calculate the other property. For example, if the modeler specifies
the diameter, then the C-factor of the resulting pipe can be given by Equation 3.4.
C
r
L
r
D
r
4.87
------------



0.54
L
i
D
i
4.87
C
i
1.85
-------------------------
i

0.54
=
120 Assembling a Model Chapter 3
(3.4)
where L = length (ft, m)
D = diameter (in., m)
C = Hazen-Williams C-factor
r = subscript referring to resulting pipe
i = subscript referring to the i-th pipe being combined
The mathematics are considerably more complicated when using the Darcy-Weisbach
equation.
In Figure 3.40, the length of the resulting pipe is 600 ft, so that if an 8-in. diameter
pipe was used, the Hazen-Williams C-factor of that pipe would be 55, and if a 6-in.
diameter was used, the C-factor would be 118. Either of these values will give the cor-
rect head loss. Minor loss coefficients and check valves can then be assigned to the
resulting pipe if needed.
Figure 3.40
Series pipe removal
different attributes
Removing Parallel Pipes. Another way to skeletonize a system is to remove
parallel pipes. (Two pipes are considered to be in parallel if they have the same begin-
ning and ending nodes.) When removing parallel pipes, one of the pipes is considered
to be the dominant pipe and the length and either the diameter or C-factor from that
pipe is used for the new equivalent pipe. Depending on whether the diameter or C-
factor is used from the dominant pipe, the other parameter is calculated using equiva-
Section 3.11 Skeletonization 121
lent pipe formulas. For example, if the diameter of the dominant pipe is used, then the
C-factor is given by the following equation:
C
r
L
r
0.54
D
r
2.63
------------
C
i
D
i
2.63
L
i
0.54
------------------
i

= (3.5)
In Figure 3.41, the length and diameter of P-40 are kept, but to account for the
removal of P-41, the capacity of P-40 is increased by increasing the C-factor.
Other factors to consider when removing parallel pipes are check valves and minor
losses. If both pipes have check valves, then the resulting pipe should also have a
check valve. Accurately assigning minor loss coefficients when determining equiva-
lent pipes, however, can be more difficult. In most cases, assigning some average
value does not cause a significant error.
Figure 3.41
Removing parallel
pipes
Removing Pipes to Break Loops. The types of pipe removal described in
the preceding sections can reduce the complexity of a model somewhat, but to dra-
matically reduce system size for typical water distribution systems, it is necessary to
actually break loops. Although two parallel pipes are considered a loop, they can be
handled with the basic action described in the previous section, and the hydraulic
capacity can be accounted for using Equation 3.5. This section applies to loops with
more than two attachments to the remainder of the system.
122 Assembling a Model Chapter 3
Consider the three-pipe loop in Figure 3.42 made up of pipes P-31, P-32, and P-33.
Removing any pipe in the loop can possibly result in a branch system that can be fur-
ther skeletonized using the methods described previously. However, in contrast to the
unique solutions that result from the pipe removal operations described in the preced-
ing sections, the results of breaking this loop by removing a pipe are different depend-
ing on which pipe is removed because the removal can have an impact on the carrying
capacity of the rest of the system.
Figure 3.42
Removing pipes in
loop
Therefore, there needs to be a rule for determining which pipe should be removed
first. Usually, it is best to remove the pipe with the least carrying capacity, which may
be defined as the smallest or the pipe with the minimum value of the quantity
CD
2.63
L
0.54
----------------- (3.6)
In Figure 3.42, pipe P-33 has the lowest carrying capacity so its removal should have
the least adverse impact to the carrying capacity of the system.
It is important to note that removing the pipe with the least carrying capacity does not
always do the least harm to the overall accuracy of the model. In some cases, a pipe
with very little capacity may be very important in some scenario and may need to be
kept in spite of its low carrying capacity.
Summary of Basic Pipe Removals. The results of the possible pipe
removal actions can be summarized as shown in Table 3.2. The first three actions are
fairly simple in that the system will end up with the correct flows and head loss. With
the fourth removal action, however, some carrying capacity is lost and removing one
pipe from a loop will give a different carrying capacity than removal of a different
pipe.
Table 3.2 Summary of pipe removal actions
Action Effect on Node Loss of System Capacity
Remove branch pipe Removes node No
Remove pipe in series Removes node No
Remove pipe in parallel No nodes removed No
Remove pipe from loop No nodes removed Yes
Section 3.11 Skeletonization 123
Removing Nonpipe Elements. Removing link-type elements other than fully
open pipes can be problematic, thus special rules must be developed for handling the
skeletonization of other network elements including pumps, tanks, closed pipes, and
valves.
A closed pipe or a pump that is not running has essentially already been skeletonized
out of the system and any effort to skeletonize it is trivial. If the element may be open,
however, then it should be treated as being open during the removal process.
Some other rules regarding the skeletonization of other network elements include the
following:
When loads are being aggregated from removed nodes, they cannot be
passed through pumps, control valves, check valves, or closed valves.
Pumps, control valves, and check valves in a branch can be trimmed and rep-
resented as an outflow from the remaining upstream system.
Pumps, control valves, and check valves can be removed from series, paral-
lel, or looped systems only if their effect can be accounted for, which is usu-
ally difficult.
If there is a check valve on a pipe in series, the resulting pipe must also have
the check valve.
Tanks are usually too important to be removed during skeletonization and no
pipes connected to tanks should be removed.
Complex Skeletonization. Skeletonizing a real system involves applying the
basic removal actions in a sequence. In general, it is best to perform the skeletonizing
actions in the order given in Table 3.2. First, remove all dead-ends or branch pipes,
then remove series pipes, then combine parallel pipes, and finally, remove loops.
After each action, it is necessary to review the network because the previous action
may have created a dead-end or a parallel pipe that did not exist previously.
Figure 3.43, which shows a network being reduced, illustrates these actions. The net-
work looks like a dead-end branch and if one were doing the skeletonization manu-
ally, a modeler would simply add together the demands and place them on node J-10.
However, it is difficult for a computer to recognize that this is a branch, and it must
first eliminate series pipes and loops to identify the branch.
124 Assembling a Model Chapter 3
Figure 3.43
Steps for applying
automatic
skeletonization to
complex pipe systems
P-53
P-54
P-52
P-51
J-10
Original System
P-53 P-52
P-51
J-10
Remove Pipe P-54 in Series
P-52
P-51
J-10
Remove Parallel Pipe P-53
J-10
Remove Branch Pipes P-51 and P-52
Stopping Criteria. Using the basic steps described in the preceding sections,
automated skeletonization can reduce any network to a handful of tanks and pumps.
In most cases, however, a user would not want this much reduction. The key to stop-
ping the skeletonization lies in defining criteria for links and nodes not to be skeleton-
ized.
Usually, the user will specify that all pipes with a certain diameter or larger will not be
skeletonized. This preserves the larger pipes in the system. The user can also specify
that certain pipes, especially those that close loops, are not to be removed (or that the
basic action of removing pipes from loops will not be carried out at all). The user can
also specify that if a pipe removal removes a node with greater than a specified
demand, then that removal action will not be carried out.
After these limits are set, the skeletonization process continues until it results in a sys-
tem that is skeletonized to the level specified by the user.
Skeletonization Conclusions
No hard and fast rules exist regarding skeletonization. It all depends on perspective
and the intended use of the model. For a utility that operates large transmission mains
and sells water to community networks, a model may be skeletonized to include only
Section 3.12 Model Maintenance 125
the source and large-diameter pipes. For a community that receives water from that
utility, the opposite may be true. Although most planning and analysis activities can
be performed successfully with a moderately skeletonized model, local fire flow eval-
uations and water quality analyses call for little to no skeletonization.
3.12 MODEL MAINTENANCE
Once a water distribution model is constructed and calibrated, it can be modified to
simulate and predict system behavior under a range of conditions. The model repre-
sents a significant investment on the part of the utility, and that investment should be
maximized by carefully maintaining the model for use well into the future.
Good record-keeping that documents model runs and history is necessary to ensure
that the model is used correctly by others or at a later date, and that time is not wasted
in deciphering and reconstructing what was done previously. There should be notes in
the model files or paper records indicating the state of the system in the various model
versions. These explanations will help subsequent users determine the best model run
to use as a starting point in future analyses.
Although the initial calibrated model reflects conditions in the current system, the
model is frequently used to test future conditions and alternative piping systems. The
scenario manager features in modeling software (see page 111) enable the user to
maintain the original model while keeping track of numerous proposed changes to the
system, some of which are never constructed. Eventually, a model file may contain
many proposed facilities and demands that fall into the following categories:
Installed
Under design or construction
To be installed later
Never to be installed
The user needs to periodically update the model file so that installed piping is accu-
rately distinguished from proposed facilities, and that facilities that will most likely
never be installed are removed from the model. The modeler also needs to be in regu-
lar contact with operations personnel to determine when new piping is placed into ser-
vice. Note that there may be a substantial lag between the time that a pipe or other
facility is placed into service, and the time that facility shows up in the system map or
GIS.
Once a master plan or comprehensive planning study is completed, model use typi-
cally becomes sporadic, though the model will still be used to respond to developer
inquiries, address operations problems, and verify project designs. Each of these spe-
cial studies involves creating and running additional scenarios. A single model even-
tually becomes cluttered with extraneous data on alternatives not selected.
A good practice in addressing these special studies is to start from the existing model
and create a new data file that will be used to study alternative plans. Once the project
design is complete, the facilities and demands associated with the selected plan
126 Assembling a Model Chapter 3
should be placed into the main model file as future facilities and demands. The ver-
sion of the model used for operational studies should not be updated until the facilities
are actually placed into service.
REFERENCES
American Water Works Association (1996). Ductile Iron Pipe and Fitting. AWWA Manual M-41, Denver,
Colorado.
Benjamin, M. M., Reiber, S. H., Ferguson, J. F., Vanderwerff, E. A., and Miller, M. W. (1990). Chemistry of
corrosion inhibitors in potable water. AWWARF, Denver, Colorado.
Caldwell, D. H., and Lawrence, W. B. (1953). Water Softening and Conditioning Problems. Industrial
Engineering Chemistry, 45(3), 535.
Eggener, C. L., and Polkowski, L. (1976). Network Modeling and the Impact of Modeling Assumptions.
Journal of the American Water Works Association, 68(4), 189.
ESRI. (2001). What is a GIS? http://www.esri.com/library/gis/abtgis/what_gis.html.
Langelier, W. F. (1936). The Analytical Control of Anti-Corrosion in Water Treatment. Journal of the
American Water Works Association, 28(10), 1500.
Maddison, L. A., and Gagnon, G. A. (1999). Evaluating Corrosion Control Strategies for a Pilot-Scale
Distribution System. Proceedings of the Water Quality Technology Conference, American Water
Works Association, Denver, Colorado.
McNeil, L. S., and Edwards, M. (2000). Phosphate Inhibitors and Red Water in Stagnant Iron Pipes.
Journal of Environmental Engineering, ASCE, 126(12), 1096.
Merrill, D. T. and Sanks, R. L. (1978). Corrosion Control by Deposition of CaCO
3
Films. AWWA, Denver,
Colorado.
Mullen, E. D., and Ritter, J. A. (1974). Potable-Water Corrosion Control. Journal of the American Water
Works Association, 66(8), 473.
Shamir, U. and Hamberg, D. (1988a). Schematic Models for Distribution Systems Design I: Combination
Concept. Journal of Water Resources Planning and Management, ASCE, 114(2), 129.
Shamir, U. and Hamberg, D. (1988b). Schematic Models for Distribution Systems Design II: Continuum
Approach. Journal of Water Resources Planning and Management, ASCE, 114(2), 141.
Volk, C., Dundore, E., Schiermann, J., LeChevallier, M. (2000). Practical Evaluation of Iron Corrosion
Control in a Drinking Water Distribution System. Water Research, 34(6), 1967.
Walski, T. M. (1999). Importance and Accuracy of Node Elevation Data. Essential Hydraulics and
Hydrology, Haestad Press, Waterbury, Connecticut.
Walski, T. M. (2000). Hydraulic Design of Water Distribution Storage Tanks. Water Distribution System
Handbook, Mays L. W., ed., McGraw Hill, New York, New York.
Discussion Topics and Problems 127
DISCUSSION TOPICS AND PROBLEMS
Read the chapter and complete the problems. Submit your work to Haestad Methods
and earn up to 11.0 CEUs. See Continuing Education Units on page xxix for more
information, or visit www.haestad.com/awdm-ceus/.
3.1 Manually find the flow rate through the system shown in the figure and compute the pressure at node
J-1. Also, find the suction and discharge pressures of the pump if it is at an elevation of 115 ft. Use
the Hazen-Williams equation to compute friction losses. Assume h
P
is in ft and Q is in cfs.
R- A
R-B
J-1
Elev = 150 ft
Pipe 3:
L=1,000 ft
D=12 in.
C=120
Pipe 2:
L=2,200 ft
D=12 in.
C=120
Pipe 1:
L=220 ft
D=16 in.
C=120
h = 225
P
- 10Q
1.50
125 ft
300 ft
128 Assembling a Model Chapter 3
3.2 Manually find the flow in each pipeline and the pressure at node J-1 for the system shown in the fig-
ure. Assume that h
P
is in m and Q is in m
3
/s and note the demand at junction J-1 of 21.2 l/s. Use the
Hazen-Williams equation to compute friction losses.
Hint: Express the flow in Pipe 3 in terms of the flow in Pipe 1 or Pipe 2.
J-1
Elev = 45.7 m
Pipe 3:
L=304.8 m
D=305 mm
C=120
Pipe 2:
L=670.6 m
D=305 mm
C=120
Pipe 1:
L=67.1 m
D=406 mm
C=120
h = 68.58
P
- 639.66Q
1.50
38.1 m
91.4 m
Q = 21.2 l/s
R - B
R - A
3.3 English Units: Manually find the discharge through each pipeline and the pressure at each junction
node of the rural water system shown in the figure. Physical data for this system are given in the
tables that follow. Fill in the tables at the end of the problem.
P-1
J-1
J-3
J-2
J-7
J-6
J-5
J-4
J-11
J-12
J-8
J-9
J-10
P-2
P-6 P-4
P-3
P-5
P-8
P-7
P-9
P-12
P-11
P-10
(Not To Scale)
R-1
Pipe Label
Length
(ft)
Diameter
(in.)
Hazen-Williams
C-factor
P-1 500 10 120
P-2 1,200 6 120
P-3 4,200 10 120
P-4 600 6 110
P-5 250 4 110
P-6 500 4 100
P-7 5,200 8 120
P-8 4,500 4 100
P-9 5,500 3 90
P-10 3,000 6 75
P-11 570 6 120
P-12 550 4 80
Node Label
Elevation
(ft)
Demand
(gpm)
R-1 1050 N/A
J-1 860 40
J-2 865 15
J-3 870 30
J-4 875 25
J-5 880 5
J-6 885 12
J-7 880 75
J-8 850 25
J-9 860 0
J-10 860 18
J-11 850 15
J-12 845 10
Pipe Label
Flow
(gpm)
Head loss
(ft)
P-1
P-2
P-3
P-4
P-5
P-6
P-7
P-8
P-9
P-10
P-11
P-12
Discussion Topics and Problems 129
Node Label
HGL
(ft)
Pressure
(psi)
J-1
J-2
J-3
J-4
J-5
J-6
J-7
J-8
J-9
J-10
J-11
J-12
130 Assembling a Model Chapter 3
SI Units: Manually find the discharge through each pipeline and the pressure at each junction node
of the rural water system shown in the figure. Physical data for this system are given in the tables
that follow. Fill in the tables at the end of the problem.
Pipe Label
Length
(m)
Diameter
(mm)
Hazen-Williams
C-factor
P-1 152.4 254 120
P-2 365.8 152 120
P-3 1,280.2 254 120
P-4 182.9 152 110
P-5 76.2 102 110
P-6 152.5 102 100
P-7 1,585.0 203 120
P-8 1,371.6 102 100
P-9 1,676.4 76 90
P-10 914.4 152 75
P-11 173.7 152 120
P-12 167.6 102 80
Node Label
Elevation
(m)
Demand
(l/s)
R-1 320.0 N/A
J-1 262.1 2.5
J-2 263.7 0.9
J-3 265.2 1.9
J-4 266.7 1.6
J-5 268.2 0.3
J-6 269.7 0.8
J-7 268.2 4.7
J-8 259.1 1.6
J-9 262.1 0
J-10 262.1 1.1
J-11 259.1 0.9
J-12 257.6 0.6
Pipe Label
Flow
(l/s)
Head loss
(m)
P-1
P-2
P-3
P-4
P-5
P-6
P-7
P-8
P-9
P-10
P-11
P-12
Node Label
HGL
(m)
Pressure
(kPa)
J-1
J-2
J-3
J-4
J-5
J-6
J-7
J-8
J-9
J-10
J-11
J-12
Discussion Topics and Problems 131
132 Assembling a Model Chapter 3
3.4 Determine the effect of placing demands at points along a pipe rather than at the end node (point D)
for the 300-m long pipe segment A-D shown in the figure. The pipe has a diameter of 150 mm and a
roughness height of 0.0001 m, and the kinematic viscosity of water at the temperature of interest is
1x10
-6
m
2
/s. The total head at Point A is 200 m, and the ground elevation along the pipe is 120 m.
The flow past point A is 9 l/s. Points A, B, C, and D are equidistant from each other.
Upstream
Point A
Intermediate
Point B
Intermediate
Point C
End
Point D
a) Assume that there is no water use along the pipe (that is, flow is 9 l/s in all segments). Determine
the head loss in each segment and the pressure head (in meters) at points B, C, and D.
b) Assume that a small amount of water is used at points B and C (typical of a pipe in a residential
neighborhood), such that the flow in the second and third segments decreases to 8 and 7 l/s,
respectively. Determine the pressures at points B, C, and D.
c) Assume that the water is withdrawn evenly along the pipe, such that the flows in the second and
third segments are 6 and 3 l/s, respectively. Find the pressures at points B, C, and D.
d) At these flows, do the pressures in the pipe vary significantly when the water use is lumped at the
endpoint versus being accounted for along the length of the pipe? Would you expect a similar
outcome at much higher flows?
Pressure in meters of water
Point Part (a) Part (b) Part (c)
B
C
D




ADVANCED WATER
DISTRIBUTION MODELING
AND MANAGEMENT

Authors
Thomas M. Walski
Donald V. Chase
Dragan A. Savic
Walter Grayman
Stephen Beckwith
Edmundo Koelle

Contributing Authors
Scott Cattran, Rick Hammond, Kevin Laptos, Steven G. Lowry,
Robert F. Mankowski, Stan Plante, John Przybyla, Barbara Schmitz

Peer Review Board
Lee Cesario (Denver Water), Robert M. Clark (U.S. EPA),
Jack Dangermond (ESRI), Allen L. Davis (CH2M Hill),
Paul DeBarry (Borton-Lawson), Frank DeFazio (Franklin G. DeFazio Corp.),
Kevin Finnan (Bristol Babcock), Wayne Hartell (Bentley Systems),
Brian Hoefer (ESRI), Bassam Kassab (Santa Clara Valley Water District),
James W. Male (University of Portland), William M. Richards
(WMR Engineering), Zheng Wu (Bentley Systems ),
and E. Benjamin Wylie (University of Michigan)



Click here to visit the Bentley Institute
Press Web page for more information

C H A P T E R
4
Water Consumption
The consumption or use of water, also known as water demand, is the driving force
behind the hydraulic dynamics occurring in water distribution systems. Anywhere
that water can leave the system represents a point of consumption, including a cus-
tomers faucet, a leaky main, or an open fire hydrant.
Three questions related to water consumption must be answered when building a
hydraulic model: (1) How much water is being used? (2) Where are the points of con-
sumption located? and (3) How does the usage change as a function of time? This
chapter addresses these questions for each of the three basic demand types described
below.
Customer demand is the water required to meet the non-emergency needs of
users in the system. This demand type typically represents the metered por-
tion of the total water consumption.
Unaccounted-for water (UFW) is the portion of total consumption that is
lost due to system leakage, theft, unmetered services, or other causes.
Fire flow demand is a computed system capacity requirement for ensuring
adequate protection is provided during fire emergencies.
Determining demands is not a straightforward process like collecting data on the
physical characteristics of a system. Some data, such as billing and production
records, can be collected directly from the utility but are usually not in a form that can
be directly entered into the model. For example, metering data are not grouped by
node. Once this information has been collected, establishing consumption rates is a
process requiring study of past and present usage trends and, in some cases, the pro-
jection of future ones.
After consumption rates are determined, the water use is spatially distributed as
demands, or loads, assigned to model nodes. This process is referred to as loading the
model. Loading is usually a multistep process that may vary depending on the prob-
lem being considered. The following steps outline a typical example of the process
the modeler might follow.
134 Water Consumption Chapter 4
1. Allocate average-day demands to nodes.
2. Develop peaking factors for steady-state runs (page 153) or diurnal curves for
EPS runs (page 155).
3. Estimate fire and other special demands.
4. Project demands under future conditions for planning and design.
This chapter presents some of the methods to follow when undertaking the process of
loading a water distribution system model.
4.1 BASELINE DEMANDS
Most modelers start by determining baseline demands to which a variety of peaking
factors and demand multipliers can be applied, or to which new land developments
and customers can be added. Baseline demands typically include both customer
demands and unaccounted-for water. Usually, the average day demand in the current
year is the baseline from which other demand distributions are built.
Data Sources
Pre-Existing Compiled Data. The first step in finding demand information
for a specific utility should always be to research the utilitys existing data. Previous
studies, and possibly even existing models, may have a wealth of background infor-
mation that can save many hours of investigation.
However, many utilities do not have existing studies or models, or may have only lim-
ited resources to collect this type of information. Likewise, models that do exist may
be outdated and may not reflect recent expansion and growth.
System Operational Records. Various types of operational records are avail-
able that can offer insight into the demand characteristics of a given system. Treat-
ment facility logs may provide data regarding long-term usage trends such as seasonal
pattern changes or general growth indications. Pumping logs and tank level charts
(such as the one shown in Figure 4.1) contain data on daily system usage, as well as
the changing pattern of demand and storage levels over time.
Water distribution systems may measure and record water usage in a variety of forms,
including
Flow information, such as the rate of production of a treatment or well
facility
Volumetric information, such as the quantity of water consumed by a
customer
Hydraulic grade information, such as the water level within a tank
The data described above are frequently collected in differing formats and require
conversion before they can be used. For example, tank physical characteristics can be
used to convert tank level data to volumes. If data describing the temporal changes in
tank levels are incorporated, volumes can be directly related to flow rates.
Section 4.1 Baseline Demands 135
Figure 4.1
Tank level chart
Courtesy of the City of Waterbury, CT Bureau of Water
Customer Meters and Billing Records. If meters are employed throughout
a system, they can be the best source of data for determining customer demands. Cus-
tomers are typically billed based on a volumetric measure of usage, with meter read-
ings taken on a monthly or quarterly cycle. Using these periodically recorded usage
volumes, customers average usage rates can be computed. Billing records, therefore,
provide enough information to determine a customers baseline demand, but not
enough to determine fluctuations in demand on a finer time scale such as that required
for extended-period simulations.
Ideally, the process of loading demand data into a model from another source would
be relatively automatic. Cesario and Lee (1980) describe an early approach to auto-
mate model loading. Coote and Johnson (1995) developed a system in Valparaiso,
Indiana in which each customer account was tied to a node in their hydraulic model.
With the increasing popularity of geographic information systems (GIS) among water
utilities, more modelers are turning to GIS to store and manipulate demand data to be
imported into the model. Stern (1995) described how Cybernet data were loaded from
a GIS in Los Angeles, and Basford and Sevier (1995) and Buyens, Bizier, and
Combee (1996) describe similar applications in Newport News, Virginia, and Lake-
136 Water Consumption Chapter 4
land, Florida, respectively. As GIS usage becomes more widespread, more utilities
will construct automated links between customer data, GIS, and hydraulic models.
Spatial Allocation of Demands
Although water utilities make a large number of flow measurements, such as those at
customer meters for billing and at treatment plants and wells for production monitor-
ing, data are usually not compiled on the node-by-node basis needed for modeling.
The modeler is thus faced with the task of spatially aggregating data in a useful way
and assigning the appropriate usage to model nodes.
The most common method of allocating baseline demands is a simple unit loading
method. This method involves counting the number of customers [or acres (hectares)
of a given land use, number of fixture units, or number of equivalent dwelling units]
that contribute to the demand at a certain node, and then multiplying that number by
the unit demand [for instance, number of gallons (liters) per capita per day] for the
applicable load classification. For example, if a junction node represents a population
of 200, and the average usage is 100 gal/day/person (380 l/day/person), the total base-
line demand for the node would be 20,000 gal/day (75,710 l/day).
In applying unit demands, the user must be careful to understand what is accounted
for by that measure. Equations 4.1, 4.2, and 4.3 show three different unit demands that
can be determined for a utility (Male and Walski, 1990).
Section 4.1 Baseline Demands 137
system-wide use = (production) / (domestic customers) (4.1)
non-industrial use = (production industrial use) / (domestic customers) (4.2)
domestic use = (domestic metered consumption) / (domestic customers) (4.3)
All three unit demands can be determined on a per capita or per account basis.
Although all three can be referred to as unit demands, each yields a different result,
and it is important that the modeler understand which unit demand is being used. The
first unit demand includes all nonemergency uses and is the largest numerical value;
the second excludes industrial uses; the third excludes industrial use and
unaccounted-for water and is the smallest. If the third unit demand is used, then
unaccounted-for water and industrial use must be handled separately from the unit
demands. This approach may be advantageous where industrial use is concentrated in
one portion of the network.
Another approach to determining the baseline demand for individual customers
involves the use of billing records. However, rarely does a system have enough
recorded information to directly define all aspects of customer usage. Even in cases
where both production records and full billing records are available, disagreements
between the two may exist that need to be resolved.
Demands in the United Kingdom
Not all water systems are universally metered as
is customary in North America. For example, in
the United Kingdom, only roughly 10 percent of
the domestic customers are metered.
Instead of metering individual customers, distribu-
tion systems in the UK are divided into smaller
zones, called District Metered Areas (DMAs),
which are isolated by valving and are fed through
a smaller number of inlet and outlet meters (WRc,
1985). The number of properties in a DMA is
known fairly precisely, and usually varies from 500
to 5,000 properties but can go as high as 10,000.
The flows are recorded using data-logging tech-
nology or telemetered to a central location.
Per capita consumption at the unmetered resi-
dences is estimated to be on the order of 150
liters per capita per day, although there is consid-
erable variation (Ofwat, 1998). Some of the
variation is attributed to different socioeconomic
classes as accounted for by ACORN (A Classifi-
cation of Residential Neighborhoods), which
classifies properties in England and Wales into
categories such as modern family housing with
higher income to poorest council estates.
Demand patterns in the UK are similar to most
other developed nations, and the patterns can be
established by DMA or groups of DMAs. Data log-
ging is used to determine individual demand pat-
terns only for the largest users.
Because most residences are not metered, unac-
counted-for water in the UK is large, but most of
this water is delivered to legitimate users and can
be estimated fairly reasonably. The amount of
actual leakage depends on pressure, burst fre-
quency, leakage control policy, and age of pipes.
Despite the differences in metering practices
between the UK and North America, loading of
the model still involves many of the same steps,
and the system metering data collected in the UK
can make calibration easier than in locations with-
out pervasive distribution metering.
138 Water Consumption Chapter 4
Two basic approaches exist for filling in the data gaps between water production and
computed customer usage: top-down and bottom-up. Both of these methods are based
on general mass-balance concepts and are shown schematically in Figure 4.2.
Figure 4.2
Approaches to model
loading
Nodes Nodes Nodes Nodes Nodes
Large Users Nodes / Service Area
Production
Meter
Routes
Billing Records Unaccounted-For
Nodes
Start: Production Data
Goal: Nodal Demands
Start: Meter Records
Top-down demand determination involves starting from the water sources (at the
top) and working down to the nodal demands. With knowledge about the produc-
tion of water and any large individual water customers, the remainder of the demand
is disaggregated among the rest of the customers. Bottom-up demand determination is
exactly the opposite, starting with individual customer billing records and summing
their influences using meter routes as an intermediate level of aggregation to deter-
mine the nodal demands.
Most methods for loading models are some variation or combination of the top-down
and bottom-up approaches, and tend to be system-specific depending on the availabil-
ity of data, the resources for data-entry, and the need for accuracy in demands. For
some systems, the decision to use top-down or bottom-up methods can be made on a
zone-by-zone basis.
Cesario (1995) uses the terms estimated consumption method and actual consumption
method to describe these two approaches. However, both methods involve a certain
level of estimation. An intermediate level of detail can be achieved by applying the
top-down approach with usage data on a meter-route-by-meter-route basis (AWWA,
1989).
Most design decisions, especially for smaller pipes, are controlled by fire flows, so
modest errors in loading have little impact on pipe sizing. The case in which loading
becomes critical is in the tracking of water quality constituents through a system,
because fire flows are not typically considered in such cases.
Section 4.1 Baseline Demands 139
Example Top-Down Demand Determination. Consider a system that serves a
community of 1,000 people and a single factory, which is metered. Over the course of a year, the total
production of potable water is 30,000,000 gallons (114,000 m
3
). The factory meter registered a usage
of 10,000,000 gallons (38,000 m
3
). Determining the average per capita residential usage in this case is
straightforward:
Total volume of residential usage = (Total usage) (Non-residential usage)
= 30,000,000 gallons 10,000,000 gallons
= 20,000,000 gallons
Residential volume usage per capita = (Total volume of residential usage) / (no. of residents)
= 20,000,000 gallons / 1,000 capita
= 20,000 gallons/capita
Residential usage rate per capita (given that prior volume calculations were for a period of one year)
= (Residential volume usage per capita) / time
= (20,000 gallons/capita/year) (1 year / 365 days)
= 55 gallons/capita/day = 210 liters/capita/day
Models usually require demands in gallons per minute or liters per second, which gives
0.038 gpm/capita = 0.0024 l/s/capita
Next, the approximate number of people (or houses) per node (e.g., 25 houses with 2.5 residents per
house = 62.5 residents/node) is determined to give average nodal demand of
2.37 gpm/node =0.15 l/s/node
These average residential nodal demands can be adjusted for different parts of town based on popula-
tion density, amount of lawn irrigation, and other factors.
Example Bottom-Up Demand Determination. Each customer account is assigned
an x-y coordinate in a GIS. Then, each account can be assigned to a node in the model based on poly-
gons around each node in the GIS. (If a GIS is not available, customer accounts can be directly
assigned to a node in the customer service information system used for billing purposes.) Then, each
account in the customer information database records can be assigned to a model node. By querying
the customer information database, the average demand at each node for any billing period can be
determined.
The billing data must now be corrected for unaccounted-for water. Consider a user who decides to
allocate unaccounted-for water uniformly to each node. The daily production is 82,000 gpd, and
metered sales are 65,000 gpd. For each node, the demand must be corrected for unaccounted-for
water. One approach is to assign unaccounted-for water in proportion to the demand at a node using:
Corrected demand = (Node consumption) [(Production) / (Metered Sales)]
For a node with a consumption of 4.2 gpm, the corrected demand is:
Corrected demand = (4.2 gpm) (82,000/65,000) = 5.3 gpm = 0.33 l/s
As can be seen in the preceding examples, bottom-up demand allocation requires a
great deal of initial effort to set up links between accounts and nodes, but after this
work is done, the loads can be recalculated easily. Of course, the corrections due to
unaccounted-for water and the fact that instantaneous demands are most likely not
equal to average demands suggest that both approaches are subject to error.
140 Water Consumption Chapter 4
Example Demand Allocation. In a detailed demand allocation, a key step is determin-
ing the customers assigned to each node. Figure 4.3 demonstrates the allocation of customer demands
to modeled junction nodes. The dashed lines represent the boundaries between junction associations.
For example, the junction labeled J-1 should have demands that represent nine homes and two com-
mercial establishments. Likewise, J-4 represents the school, six homes, and one commercial building.
Figure 4.3
Allocating demands
to network junctions
Commercial
Establishments
Homes
School
J-1 J-2
J-3 J-4 J-5
Node Service
Area Boundary (typ.)
Following demand allocation, the modeler must ensure that demands have been assigned to junction
nodes in such a way that (1) the sums of the nodal demands system-wide and in each pressure zone
are in agreement with production records, and (2) the spatial allocation of demands closely approxi-
mates actual demands.
When working with high-quality GIS data, the modeler can much more precisely
assign demands to nodes. Nodal demands can be loaded using several GIS-related
methodologies, ranging from a simple inverse-pipe-diameter allocation model to a
comprehensive polygon overlay. The inverse-pipe-diameter approach assumes that
demands are associated with small-diameter pipes, whereas large-diameter pipes are
mainly used for transmission and thus should have less weight associated with
them. More detailed methods make use of extensive statistical data analysis and GIS
processing by combining layers of data that account for variables such as population
changes over time, land use, seasonal changes, planning, and future development
rates. Davis and Brawn (2000) describe an approach they employed to allocate
demands using a GIS.
Using GIS for Demand Allocation
As discussed previously, an integral part of creating a water distribution model is the
accurate allocation of demands to the node elements within the model. The spatial
Section 4.1 Baseline Demands 141
analysis capabilities of GIS make it a logical tool for the automation of the demand
allocation process.
The following sections provide descriptions of some of the automated allocation strat-
egies that have been used successfully.
Meter Assignment. This allocation strategy uses the spatial analysis capabili-
ties of GIS to assign geocoded (possessing coordinate data based on physical location,
such as an x-y coordinate) customer meters to the nearest demand node. Therefore,
this type of model loading is a point-to-point demand allocation technique, meaning
that known point demands (customer meters) are assigned to network demand points
(demand nodes). Meter assignment is the simplest technique in terms of required data,
because there is no need for service polygons to be applied (see Figure 4.4).
However, meter assignment can prove less accurate than the more complex allocation
strategies because nearest is determined by straight-line proximity between the
demand node and the consumption meter. Piping routes are not considered, so the
nearest demand node may not be the location from which the meter actually receives
its flow. In addition, the actual location of the service meter may not be known. Ide-
ally, these meter points should be placed at the location of the tap, but the centroid of
the building or land parcel may be all that is known about a customer account.
Figure 4.4
Meter assignment
142 Water Consumption Chapter 4
Meter Aggregation. Meter aggregation is the technique of assigning all meters
within a service polygon to a specified demand node (see Figure 4.5). Service poly-
gons define the service area for each of the demand junctions.
Meter aggregation is a polygon-to-point allocation technique because the service
areas are contained in a GIS polygon layer and the demand junctions are contained in
a point layer. The demands associated with each of the service-area polygons are
assigned to the respective demand node points.
Because of the need for service polygons, the initial setup for this approach is more
involved than for the simpler meter assignment strategy, with the tradeoff being
greater control over the assignment of meters to demand nodes. Automated construc-
tion of the service polygons may not produce the desired results, so it may be neces-
sary to manually adjust the polygon boundaries, especially at the edges of the
drawing.
Figure 4.5
Meter aggregation
Flow Distribution. This strategy involves distributing a lump-sum demand
among a number of service polygons (service areas) and, by extension, their associ-
ated demand nodes. The lump-sum area is a polygon for which the total (lump-sum)
demand of all of the service areas (and their demand nodes) is known (metered), but
Section 4.1 Baseline Demands 143
the distribution of the total demand among the individual nodes is not. Lump-sum
areas can be based on pressure zones, meter routes, or other criteria.
The known demand within the lump-sum area is divided among the service polygons
within the area using one of two techniques: equal distribution or proportional distri-
bution. The equal distribution option simply divides the known demand evenly
between the demand nodes. For example, in Figure 4.6, the total demand in meter
route A may be 55 gpm (3.48 l/s), and the total demand in meter route B may be 72
gpm (4.55 l/s). Since there are 11 nodes in meter route A and 8 nodes in meter route
B, the demand at each node will be 5 gpm (0.32 l/s) and 9 gpm (0.57 l/s), respectively.
Figure 4.6
Equal flow
distribution
The proportional distribution option divides the lump-sum demand among the service
polygons based on one of two attributes of the service polygons: the area or the popu-
lation. That is, the greater the percentage of the lump-sum area or population that a
service polygon contains, the greater the percentage of total demand that will be
assigned to that service polygon.
Each service polygon has an associated demand node, and the demand that is calcu-
lated for each service polygon is assigned to this demand node. For example, if a ser-
vice polygon makes up 50 percent of the lump-sum polygons area, then 50 percent of
the demands associated with the lump-sum polygon will be assigned to the demand
node associated with that service polygon.
144 Water Consumption Chapter 4
Flow distribution strategies require the definition of lump-sum area or population
polygons, service polygons, and their related demand nodes.
Sometimes, a combination of demand allocation methods is recommended. One case
where this technique is particularly helpful is in accounting for unaccounted-for
water. A meter assignment or meter aggregation method can be used to distribute the
normal demands, and a flow distribution technique can be used in addition to assign
the unaccounted-for water.
Point Demand Assignment. A point demand assignment technique is used to
directly assign a demand to a demand node. This strategy is primarily a manual oper-
ation, and is used to assign large (generally industrial or commercial) water users to
the demand node that serves the consumer in question. This technique is unnecessary
if all demands are accounted for by using one of the other allocation strategies.
Projection of Future Demands. Automated techniques have also been devel-
oped to assist in the assignment of future demands to nodes. These are similar to flow
distribution allocation except that the type of base layer that is used to intersect with
the service layer may contain information other than average-day demands.
Demand projection relies on a polygon layer that contains data regarding expected
future conditions. Some data types that can be used for this include future land use
and projected population, in combination with a demand density (for example, gal-
lons per capita per day), with the polygons based on traffic analysis zones, census
tracts, planning districts, or another classification. Many of these data types do not
include demand information, so demand density is required to translate the informa-
tion contained in the future-condition polygons into projected demand values.
Methods of using water-use data based on population or land use involve overlaying
those polygons on node service-area polygons and are described in more detail in
Chapter 12.
Categorizing Demands
Sometimes water users at a single node fall into several categories, and the modeler
would like to keep track of these categories within the model. Composite demands
enable the modeler to do this type of tracking. The modeler can selectively search for
all demands of a certain type (for example, residential or industrial) and make adjust-
ments. The modeler can also make changes to the characteristics of an entire category,
and all of the customers of that type will automatically be modified.
Composite Demands. Whether a unit-loading-based or a billing-record-based
method is used to generate the baseline demand, the user may need to convert it into a
composite demand at a particular node. This conversion is necessary since a junction
node does not always supply a single customer type. When more than one demand
type is served by a particular junction, the demand is said to be a composite. Deter-
mining the total rate of consumption for a junction node with a composite demand is a
simple matter of summing the individual components. Composite demands are also a
way of keeping track of unaccounted-for water independent of the other demands at a
node.
Section 4.1 Baseline Demands 145
When temporal patterns are applied to composite demands, the total demand for a
junction at any given time is equal to the sum of each baseline demand times its
respective pattern multiplier. It is also possible with most software packages to assign
a different pattern to the different components of the composite demand.
Q
i t ,
B
i j ,
P
i j t , ,
j
_
= (4.4)
where Q
i, t
= total demand at junction i at time t (cfs, m
3
/s)
B
i, j
= baseline demand for demand type j at junction i (cfs, m
3
/s)
P
i, j, t
= pattern multiplier for demand type j at junction i at time t
Nomenclature. Depending on the scale of the model, the demand type may con-
sist of such broad categories as residential, commercial, and industrial, or be
broken down into a finer level of detail with categories such as school, restaurant,
multifamily dwelling, and so on.
An issue that arises when discussing demands is that each utility classifies customers
differently. For example, an apartment may be a residential account at one utility, a
commercial account at another, and a multifamily residential account at yet
another. Schools may be classified as institutional, commercial, public, or sim-
ply schools. A modeler working for a utility can easily adapt to the naming conven-
tions, but a consultant who works with many utilities may have a difficult time
keeping track of the nomenclature when moving from one system to another.
Mass Balance Technique
Regardless of whether a modeler is studying the entire system, one particular pressure
zone, or an individual customer, mass balance techniques are useful for determining
changes in demand occurring on a finer time scale than a monthly billing cycle. For a
water distribution system, a mass balance simply indicates that what goes into the sys-
tem must be equal to what comes out of the system or zone (accounting for changes in
storage). In equation form, this can be stated as follows:
Q
demand
Q
inf low
Q
outf low
AV
st orage
At + = (4.5)
where Q
inflow
= average rate of production (cfs, m
3
/s)
Q
demand
= average rate of demand (cfs, m
3
/s)
Q
outflow
= average outflow rate (cfs, m
3
/s)
AV
st orage
= change in storage within the system (ft
3
, m
3
)
At = time between volume measurements (s)
Note that the rates of production and demand in the above equation are representative
of the average flow rates over the time period. The change in storage, however, is
146 Water Consumption Chapter 4
found by taking the difference between storage volumes at the beginning and end of
the time period for each tank, as follows:
AV
st orage
V
i t At + ,
V
i t ,
( )
i
_
= (4.6)
where V
i t At + ,
= storage volume of tank i at time t+A t (ft
3
, m
3
)
V
i, t
= storage volume of tank i at time t (ft
3
, m
3
)
When calculating volume changes in storage, a sign convention must apply. If the vol-
ume in storage decreased during the time interval, then that volume is added to the
inflows, and if it increased over the time period, then it is subtracted from inflows.
For upright cylindrical tanks (or any tank with vertical sides), the change in storage
can be determined directly from the change in tank level, as follows:
AV
st orage
H
i t At + ,
H
i t ,
( )A
i t ,
i
_
= (4.7)
where H
i t At + ,
= water level at beginning of times step t+A t at tank i (ft, m)
H
i,t
= water level at beginning of times step t at tank i (ft, m)
A
i,t
= surface area of tank i during time step t (ft
2
, m
2
)
Example Mass Balance. Consider a pressure zone with a single cylindrical tank having
a diameter of 40 feet. At the beginning of a daily monitoring interval, the water level is at 28.3 feet,
and at the beginning of the next day it is 29.1 feet. During that time, the total flow into that zone is
determined to be 455 gallons per minute, and there is no outflow to other zones. What is the total
average daily demand within the zone?
Knowing the tanks diameter, its area is found to be
A
tD
2
4
----------
t 40 ( )
2
4
----------------- 1256 = = = ft
2
The change in storage is then found as
AV = A(H
i+1
- H
i
) = 1256 ft
2
(29.1 ft 28.3 ft) 7.48 gal/ft
3
= 7,516 gal
Storage in the tank increased over the monitoring period, thus the sign convention dictates that the
flows to storage be subtracted from the total inflow. With the change in storage and the average
inflow, the average zone demand occurring over the hourly monitoring period is
Q 455
7516
60 24
------------------ 449.8 = = gpm
This answer makes sense because the average zone demand must be smaller than the average inflow
for the tank to fill during the monitoring period.
Section 4.1 Baseline Demands 147
Using Unit Demands
In the case of new water customers, flows can usually be estimated based on similar
customers in the community. Numerous investigators have compiled typical water
consumption for different types of facilities. To use this data, the modeler needs to
determine the number of units (for example, number of rooms in a hotel or number of
seats in a restaurant) and multiply by the typical unit flow to determine the average
daily flow from that establishment.
Table 4.1 provides typical unit loads for a number of different types of users. Ranges
are given because there is considerable variation between establishments within a
given category.
Table 4.1 Typical rates of water use for various establishments
Range of Flow
User (l/person or unit/day) (gal/person or unit/day)
Airport, per passenger 1020 35
Assembly hall, per seat 610 23
Bowling alley, per alley 60100 1626
Camp
Pioneer type 80120 2132
Childrens, central toilet and bath 160200 4253
Day, no meals 4070 1118
Luxury, private bath 300400 79106
Labor 140200 3753
Trailer with private toilet and bath,
per unit (2 1/2 persons)
500600 132159
Country clubs
Resident type 300600 79159
Transient type serving meals 60100 1626
Dwelling unit, residential
Apartment house on individual well 300400 79106
Apartment house on public water supply,
unmetered
300500 79132
Boardinghouse 150220 4058
Hotel 200400 53106
Lodging house and tourist home 120200 3253
Motel 400600 106159
Private dwelling on individual well or
metered supply
200600 53159
Private dwelling on public water
supply, unmetered
400800 106211
Factory, sanitary wastes, per shift 40100 1126
Table extracted from Ysuni, 2000 based on Metcalf and Eddy, 1979
148 Water Consumption Chapter 4
Other investigators have linked water use in nonresidential facilities to the Standard
Industrial Classification (SIC) codes for each industry as shown in Table 4.2. To use
this table, the modeler determines the number of employees in an industry and multi-
plies the number by the use rate given in the table. As is the case with Table 4.1, there
will be considerable variation about the typical values given Table 4.2.
Fairground (based on daily attendance) 26 12
Institution
Average type 400600 106159
Hospital 7001200 185317
Office 4060 1116
Picnic park, with flush toilets 2040 511
Restaurant (including toilet)
Average 2540 711
Kitchen wastes only 1020 35
Short order 1020 35
Short order, paper service 48 12
Bar and cocktail lounge 812 23
Average type, per seat 120180 3248
Average type, 24 h, per seat 160220 4258
Tavern, per seat 60100 1626
Service area, per counter seat (toll road) 10001600 264423
Service area, per table seat (toll road) 600800 159211
School
Day, with cafeteria or lunchroom 4060 1116
Day, with cafeteria and showers 6080 1621
Boarding 200400 53106
Self-service laundry, per machine 10003000 264793
Store
First 7.5 m (25 ft) of frontage 16002000 423528
Each additional 7.5 m of frontage 14001600 370423
Swimming pool and beach, toilet and shower 4060 1116
Theater
Indoor, per seat, two showings per day 1020 35
Outdoor, including food stand, per car
(3 1/3 persons)
1020 35
Table 4.1 (cont.) Typical rates of water use for various establishments
Range of Flow
User (l/person or unit/day) (gal/person or unit/day)
Table extracted from Ysuni, 2000 based on Metcalf and Eddy, 1979
Section 4.1 Baseline Demands 149
Table 4.2 Average rates of nonresidential water use from establishment-level data
Category SIC Code
Use Rate
(gal/employee/day)
Sample
Size
Construction 31 246
General building contractors 15 118 66
Heavy construction 16 20 30
Special trade contractors 17 25 150
Manufacturing 164 2790
Food and kindred products 20 469 252
Textile mill products 22 784 20
Apparel and other textile products 23 26 91
Lumber and wood products 24 49 62
Furniture and fixtures 25 36 83
Paper and allied products 26 2614 93
Printing and publishing 27 37 174
Chemicals and allied products 28 267 211
Petroleum and coal products 29 1045 23
Rubber and miscellaneous plastics products 30 119 116
Leather and leather products 31 148 10
Stone, clay, and glass products 32 202 83
Primary metal industries 33 178 80
Fabricated metal products 34 194 395
Industrial machinery and equipment 35 68 304
Electronic and other electrical equipment 36 95 409
Transportation equipment 37 84 182
Instruments and related products 38 66 147
Miscellaneous manufacturing industries 39 36 55
Transportation and public utilities 50 226
Railroad transportation 40 68 3
Local and interurban passenger transit 41 26 32
Trucking and warehousing 42 85 100
U.S. Postal Service 43 5 1
Water transportation 44 353 10
Transportation by air 45 171 17
Transportation services 47 40 13
Communications 48 55 31
Electric, gas, and sanitary services 49 51 19
Wholesale trade 53 751
Wholesale tradedurable goods 50 46 518
Wholesale tradenondurable goods 51 87 233
Table from Dziegielweski, Opitz, and Maidment, 1996
150 Water Consumption Chapter 4
Retail trade 93 1044
Building materials and garden supplies 52 35 56
General merchandise stores 53 45 50
Food stores 54 100 90
Automotive dealers and service stations 55 49 498
Apparel and accessory stores 56 68 48
Furniture and home furnishings stores 57 42 100
Eating and drinking places 58 156 341
Miscellaneous retail 59 132 161
Finance, insurance, and real estate 192 238
Depository institutions 60 62 77
Nondepository institutions 61 361 36
Security and commodity brokers 62 1240 2
Insurance carriers 63 136 9
Insurance agents, brokers, and service 64 89 24
Real estate 65 609 84
Holding and other investment offices 67 290 5
Services 137 1878
Hotels and other lodging places 70 230 197
Personal services 72 462 300
Business services 73 73 243
Auto repair, services, and parking 75 217 108
Miscellaneous repair services 76 69 42
Motion pictures 78 110 40
Amusement and recreation services 79 429 105
Health services 80 91 353
Legal services 81 821 15
Educational services 82 110 300
Social service 83 106 55
Museums, botanical, zoological gardens 84 208 9
Membership organizations 86 212 45
Engineering and management services 87 58 5
Services, NEC 89 73 60
Public administration 106 25
Executive, legislative, and general 91 155 2
Justice, public order, and safety 92 18 4
Administration of human resources 94 87 6
Table 4.2 (cont.) Average rates of nonresidential water use from establishment-level data
Category SIC Code
Use Rate
(gal/employee/day)
Sample
Size
Table from Dziegielweski, Opitz, and Maidment, 1996
Section 4.1 Baseline Demands 151
Unaccounted-For Water
Ideally, if individual meter readings are taken for every customer, they should exactly
equal the amount of water that is measured leaving the treatment facility. In practice,
however, this is not the case. Although inflow does indeed equal outflow, not all of the
outflows are metered. These lost flows are referred to as unaccounted-for water
(UFW).
There are many possible reasons why the sum of all metered customer usage may be
less than the total amount of water produced by the utility. The most common reasons
for discrepancies are leakage, errors in measurement, and unmetered usage. Ideally,
customer demands and unaccounted-for water should be estimated separately. In this
way, a utility can analyze the benefits of reducing unaccounted-for water.
Unaccounted-for water must be loaded into the model just like any other demand.
However, the fact that it is unaccounted-for means that the user does not know where
to place it. Usually, the user simply calculates total unaccounted-for water and divides
that quantity equally among all nodes. If the modeler knows that one portion of a sys-
tem has a greater likelihood of leakage because of age, then more unaccounted-for
water can be placed within that section.
Leakage. Leakage is frequently the largest component of UFW and includes dis-
tribution losses from supply pipes, distribution and trunk mains, services up to the
meter, and tanks. The amount of leakage varies from system to system, but there is a
general correlation between the age of a system and the amount of UFW. Newer sys-
tems may have as little as 5 percent leakage, while older systems may have 40 percent
leakage or higher. Leakage tends to increase over time unless a leak detection and
repair program is in place. Use of acoustic detection equipment to listen for leaks is
shown in Figure 4.7.
Other factors affecting leakage include system pressure (the higher the pressure, the
more leakage), burst frequencies of mains and service pipes, and leakage detection
and control policies. These factors make leakage very difficult to estimate, even with-
out the complexity of approximating other UFW causes. If better information is not
available, UFW is usually assigned uniformly around the system.
Environmental quality and housing 95 101 6
Administration of economic programs 96 274 5
National security and international affairs 97 445 2
Table 4.2 (cont.) Average rates of nonresidential water use from establishment-level data
Category SIC Code
Use Rate
(gal/employee/day)
Sample
Size
Table from Dziegielweski, Opitz, and Maidment, 1996
152 Water Consumption Chapter 4
Figure 4.7
Use of leak detection
equipment
Meter Under Registration. Flow measurement errors also contribute to UFW.
Flow measurements are not always exact, and thus metered customer usage may con-
tain inaccuracies. Some flow meters under-register usage at low flow rates, especially
as they get older.
Unmetered Usage. Systems may have illegal connections or other types of
unmetered usage. Not all unmetered usage is indicative of water theft. Fire hydrants,
blow-offs, and other maintenance appurtenances are typically not metered.
4.2 DEMAND MULTIPLIERS
By definition, baseline demands during a steady-state simulation do not change over
time. However, in reality, water demand varies continuously over time according to
several time scales:
Daily. Water use varies with activities over the course of a day.
Weekly. Weekend patterns are different from weekdays.
Seasonal. Depending on the extent of outdoor water use or seasonal
changes, such as tourism, consumption can vary significantly from one sea-
son to another.
Section 4.2 Demand Multipliers 153
Long-term. Demands can grow due to increases in population and the
industrial base, changes in unaccounted-for water, annexation of areas previ-
ously without service, and regionalization of neighboring water systems.
The modeler needs to be cognizant of the impacts of temporal changes on all of these
scales. These time-varying demands are handled in the model by either
Steady-state runs for a particular condition, or
Extended-period model runs
For extended-period simulations, the model requires both baseline demand data and
information on how demands vary over time. Modeling of these temporal variations is
described in the next section.
In steady-state runs, the user can build on the baseline demand by using multipliers
and/or assigning different demands to specific nodes. Fortunately, the entire demand
allocation need not be redone.
The following are some examples of demand events frequently considered:
Average-day demand: The average rate of demand for an average day (past,
present, or future)
Maximum-day demand: The average rate of use on the maximum usage
day (past, present, or future)
Peak-hour demand: The average rate of usage during the maximum hour of
usage (past, present, or future)
Maximum day of record: The highest average rate of demand for the his-
torical record
Peaking Factors
For some consumption conditions (especially predicted consumption conditions),
demands can be determined by applying a multiplication factor or a peaking factor.
For example, a modeler might determine that future maximum day demands will be
double the average-day demands for a particular system. The peaking factor is calcu-
lated as the ratio of discharges for the various conditions. For example, the peaking
factor applied to average-day demands to obtain maximum day demands can be found
by using Equation 4.8.
PF Q
max
Q
avg
= (4.8)
where PF = peaking factor between maximum day and average-day demands
Q
max
= maximum day demands (cfs, m
3
/s)
Q
avg
= average-day demands (cfs, m
3
/s)
Determining system-wide peaking factors is fairly easy because most utilities keep
good records on production and tank levels. However, peaking factors for different
types of demands applied at individual nodes are more difficult to determine, because
154 Water Consumption Chapter 4
individual nodes do not necessarily follow the same demand pattern as the system as a
whole.
Peaking factors from average day to maximum day tend to range from 1.2 to 3.0, and
factors from average day to peak hour are typically between 3.0 and 6.0. Of course,
these values are system-specific, so they must be determined based on the demand
characteristics of the system at hand.
Fire flows represent a special type of peaking condition, and they are described on
page 165. Fire flows are usually added to maximum day flows when evaluating the
capacity of the system for fire fighting.
Demands in Systems with High Unaccounted-For Water. Using glo-
bal demand multipliers for projections in systems with high unaccounted-for water is
based on the assumption that the relative amount of unaccounted-for water will
remain constant in the future. Unaccounted-for water can also be treated as one of the
parts of a composite demand, as discussed on page 144. If unaccounted-for water is
reduced, then the utility will see higher peaking factors because unaccounted-for
water tends to flatten out the diurnal demand curve. Walski (1999) describes a method
for correcting demand multipliers for systems where leakage is expected to change
over time.
M
A
---- -
M
A
-----
\ .
| |
c
Q
c
L +
Q
c
L +
----------------------------- = (4.9)
where M/A = corrected multiplier
(M/A)
c
= multiplier for consumptive users only
Q
c
= water use through customer meters in future (cfs, m
3
/s)
L = leakage in future (cfs, m
3
/s)
Example Peaking Factors. If the multiplier for metered customers (M/A)
c
is 2.1, and
the metered demand (Q
c
) is projected to be 2.4 MGD in a future condition, then the overall multiplier
can be determined based on estimated future leakage as shown in the following table.
Leakage
(MGD)
M/A
0.0 2.1
0.5 1.9
1.0 1.8
Because leakage contributes the same to average and peak demands, the peak demand multipliers
increase as leakage decreases. The numerical value of (M/A)
c
can be calculated using current year
data and Equation 4.10.
M
A
-----
\ .
| |
c
M
A
---- - Q
c
L + ( ) L
Q
c
------------------------------------ = (4.10)
Section 4.3 Time-Varying Demands 155
The L and Q values are based on current year actual values. For example, in this problem, say that the
current year overall multiplier is 1.8, the metered demand is 1.5 MGD, and the leakage is 0.6 MGD.
The multiplier for metered consumption is then
M
A
-----
\ .
| |
c
1.8 1.5 0.6 + ( ) 0.6
1.5
----------------------------------------------- 2.1 = =
Commercial Building Demands. A means of estimating design demands for
proposed commercial buildings is called the Fixture Unit Method. If the nature of the
customer/building is known, and the number and types of water fixtures (toilets, dish-
washers, drinking fountains, and so on) can be calculated, then the peak design flow
can be determined. The fixture unit method accounts for the fact that it is very
unlikely that all of the fixtures in a building will be operated simultaneously. Chapter
9 contains more information on using this method (see page 399).
4.3 TIME-VARYING DEMANDS
Water usage in municipal water distribution systems is inherently unsteady due to
continuously varying demands. In order for an extended period simulation to accu-
rately reflect the dynamics of the real system, these demand fluctuations must be
incorporated into the model.
The temporal variations in water usage for municipal water systems typically follow a
24-hour cycle called a diurnal demand pattern. However, system flows experience
changes not only on a daily basis, but also weekly and annually. As one might expect,
weekend usage patterns often differ from weekday patterns. Seasonal differences in
water usage have been related to climatic variables such as temperature and precipita-
tion, and also to the changing habits of customers, such as outdoor recreational and
agricultural activities occurring in the summer months.
Diurnal Curves
Each city has its own unique level of usage that is a function of recent climatic condi-
tions and the time of day. (Economic growth also influences demands, but its effect
occurs over periods longer than the typical modeling time horizon, and it is accounted
for using future demand projections.) Figure 4.8 illustrates a typical diurnal curve for
a residential area. There is relatively low usage at night when most people sleep,
increased usage during the early morning hours as people wake up and prepare for the
day, decreased usage during the middle of the day, and finally, increased usage again
in the early evening as people return home.
For other water utilities and other types of demands, the usage pattern may be very
different. For example, in some areas, residential irrigation occurs overnight to mini-
mize evaporation, which may cause peak usage to occur during the predawn hours.
For small towns that are highly influenced by a single industry, the diurnal pattern
may be much more pronounced because the majority of the population follows a sim-
ilar schedule. For example, if a large water-using industry runs 24 hours per day, the
overall demand pattern for the system may appear relatively flat because the steady
industrial usage is much larger than peaks in the residential patterns.
156 Water Consumption Chapter 4
Figure 4.8
A typical diurnal
curve
Developing System-Wide Diurnal Curves
A system-wide diurnal curve can be constructed using the same mass balance tech-
niques discussed earlier in this chapter. The only elaboration is that the mass balance
is performed as a series of calculations, one for each hydraulic step of an EPS simula-
tion.
Time Increments. The amount of time between measurements has a direct corre-
lation to the resolution and precision of the constructed diurnal curve. If measure-
ments are only available once per day, then only a daily average can be calculated.
Likewise, if measurements are available in hourly increments, then hourly averages
can be used to define the pattern over the entire day.
If the modeler tries to use a time step that is too small, small errors in tank water level
can lead to large errors in water-use calculations. This type of error is explained fur-
ther in Walski, Lowry, and Rhee (2000). Modeling of hydraulic time steps smaller
than one hour is usually only justified in situations in which tank water levels change
rapidly. Even if facility operations (such as pump cycling) occur frequently, it may
still be acceptable for the demand pattern time interval to be longer than the hydraulic
time step.
The modeler should be aware that incremental measurements can still overlook a peak
event. For example, consider something as simple as determining the peak-hour
demand (the highest average demand over any continuous one-hour period). If mea-
surements are taken every hour on the hour, then the determination of the computed
peak hour will only be accurate if the actual peak begins and ends right on an even
Section 4.3 Time-Varying Demands 157
hour increment (such as a peak occurring from 7:00 to 8:00 a.m.). The modeler will
underestimate peak hour usage if the true peak occurs, for example, from 7:15 to 8:15
a.m. The diurnal demand curve in Figure 4.9 illustrates this point. As can also be seen
in Figure 4.9, as time increments become smaller, peak flows become higher (for
instance, the 15-minute peak is higher than the one-hour peak).
Figure 4.9
Missed peak on a
diurnal curve due to
model time step
U
s
a
g
e
Time
Diurnal Curve Peak Missed
5am 6am 7am 8am 9am 10am 11am 11pm 12am
Hourly Measurements
Developing Customer Diurnal Curves
Frequently, developing a diurnal curve for a specific customer requires more informa-
tion than can be extracted from typical billing records. In these situations, more inten-
sive data collection methods are needed to portray the time-variant nature of the
demands.
Data Logging for Customer Usage. Manually reading a customers water
meter at frequent intervals would obviously be a tedious and expensive undertaking.
The process of data logging refers to the automated gathering of raw data in the field.
These data are later compiled and analyzed for a variety of purposes, among them the
creation of diurnal demand curves. Various applications of data logging are described
in papers by Brainard (1994); Rhoades (1995); and DeOreo, Heaney, and Mayer
(1996).
There have been many recent advances in data-logging technology, making it a reli-
able and fairly inexpensive way to record customer water usage. Figure 4.10 illus-
158 Water Consumption Chapter 4
trates a typical meter/data logger setup. Utilities can now easily place a data logger on
a customers meter and determine that customers consumption pattern. While it
would be nice to have such a detailed level of data on all customers, the cost of obtain-
ing the information is only justifiable for large customers and a sampling of smaller
ones.
Figure 4.10
A typical meter/data
logger setup
Meter-Master Model 100EL Flow Recorder manufactured by F.S. Brainard & Company
Representative Customers. Although it is possible to study a few customers
in detail and extend the conclusions of that study to the rest of the system, this type of
data extrapolation has some inherent dangers. The probability of selecting the per-
fect average customer is small, and any deviation from the norm or error in measure-
ment will be compounded when it is applied to an entire community. As with all
statistical data collection methods, the smaller the sample size, the less confidence
there can be in the results.
There are also applications in which use of a representative customer is inappropriate
under any circumstances. With large industries, for example, there may be no rela-
tionship at all between the volumes and patterns of usage even though they share a
similar zoning classification. Therefore, demands for large consumers (industries,
hospitals, hotels, and so on) and their diurnal variations should be individually deter-
mined.
Even if data logging cannot be applied to all customers, studying the demands of large
consumers and applying the top-down demand determination concept to the smaller
consumers can still yield reasonable demand calculations. The large customer data is
subtracted from the overall system or zone usage, and the difference in demand is
attributable to the smaller customers.
It is impossible to know with absolute certainty when water will be used or how much
water is used in a short period of time, even though usage per billing period is known
exactly. Bowen, Harp, Baxter, and Shull (1993) collected data from single and multi-
Section 4.3 Time-Varying Demands 159
family residential customers in several U.S. cities. These demand patterns can be used
as a starting point for assigning demand patterns to residential nodes.
Buchberger and Wu (1995) and Buchberger and Wells (1996) developed a stochastic
model for residential water demands and verified it by collecting extensive data on
individual residential customers. The model is particularly useful for evaluating the
hydraulics of dead-ends and looped systems in the periphery of distribution networks.
The researchers found that the demand at an individual house cannot simply be multi-
plied by the number of houses to determine the demand in a larger area. The methods
that they developed provide a way of combining the individual stochastic demands
from individual customers who are brushing their teeth or running their washing
machines, dishwashers, and so on into the aggregate for use in a larger area over a
longer time interval.
In general, hotels and apartments have demand patterns similar to those of residential
customers, office buildings have demand patterns corresponding to 8 a.m. to 5 p.m.
operations, and retail area demand patterns reflect 9 a.m. to 9 p.m. operations. Every
large industry that uses more than a few percentage points of total system production
should have an individual demand pattern developed for it.
Defining Usage Patterns within a Model
Usage could be defined directly by describing a series of actual flow versus time
points for each junction in the system. One shortcoming of this type of definition is
that it does not offer much data reuse for nodes with similar usage patterns. Conse-
quently, most hydraulic models express demands by using a constant baseline demand
multiplied by a dimensionless demand pattern factor at each time increment.
A demand multiplier is defined as
Mult
i
Q
i
Q
base
= (4.11)
where Multi
i
= demand multiplier at the i
th
time step
Q
i
= demand in i
th
time step (gpm, m
3
/s)
Q
base
= base demand (gpm, m
3
/s)
The series of demand pattern multipliers models the diurnal variation in demand and
can be reused at nodes with similar usage characteristics. The baseline demand is
often chosen to be the average daily demand (although peak day demand or some
other value can be used). Assuming a baseline demand of 200 gpm, Table 4.3 illus-
trates how nodal demands are computed using a base demand and pattern multipliers.
Table 4.3 Calculation of nodal demands using pattern multipliers
Time
Pattern
Multiplier
Demand
0:00 0.7
200 gpm 0.7 = 140 gpm
160 Water Consumption Chapter 4
As one can imagine, usage patterns are as diverse as the customers themselves. Figure
4.11 illustrates just how different diurnal demand curves for various classifications
can be. A broad zoning classification, such as commercial, may contain differences
significant enough to warrant the further definition of subcategories for the different
types of businesses being served. For instance, a hotel may have a demand pattern that
resembles that of a residential customer. A dinner restaurant may have its peak usage
during the late afternoon and evening. A clothing store may use very little water,
regardless of the time of day. Water usage in an office setting may coincide with cof-
fee breaks and lunch hours.
Figure 4.11
Diurnal curve for
different user
categories
Single Family Businesses
Restaurant Factory
D
e
m
a
n
d
M
u
l
t
i
p
l
i
e
r
D
e
m
a
n
d
M
u
l
t
i
p
l
i
e
r
D
e
m
a
n
d
M
u
l
t
i
p
l
i
e
r
D
e
m
a
n
d
M
u
l
t
i
p
l
i
e
r
Time Time
Time Time
There will sometimes be customers within a demand classification whose individual
demand patterns differ significantly from the typical demand pattern assigned to the
classification as a whole. For most types of customers, the impact such differences
have on the model is insignificant. For other customers, such as industrial users, errors
in the usage pattern may have a large impact on the model. In general, the larger the
individual usage of a customer, the more important it is to ensure the accuracy of the
consumption data.
Stepwise and Continuous Patterns. In a stepwise demand pattern, demand
multipliers are assumed to remain constant over the duration of the pattern time step.
1:00 1.1
200 gpm 1.1 = 220 gpm
2:00 1.8
200 gpm 1.8 = 360 gpm
Table 4.3 Calculation of nodal demands using pattern multipliers
Time
Pattern
Multiplier
Demand
Section 4.3 Time-Varying Demands 161
A continuous pattern, on the other hand, refers to a pattern that is defined indepen-
dently of the pattern time step. Interpolation methods are used to compute multiplier
values at intermediate time steps. If the pattern time step is reset to a smaller or larger
value, the pattern multipliers are automatically recalculated. The pattern multiplier
value is updated by linearly interpolating between values occurring along the continu-
ous curve at the new time step interval. The result is a more precise curve fit that is
independent of the time step specified, as shown in Figure 4.12.
Figure 4.12
Stepwise and
continuous pattern
variation
P
a
t
t
e
r
n
M
u
l
t
i
p
l
i
e
r
Time of Day
Continuous
Stepwise
1.5
1.0
0.5
0.0
0 8 12 18 24
Average
For example, the pattern from Table 4.3 can be extended to show how a typical model
might determine multipliers after the time step had been changed from 1 hour to 15
minutes over the time period 0:00 to 1:00. As Table 4.4 shows, a pattern multiplier for
an intermediate time increment in a continuous pattern can differ significantly from
its stepwise pattern counterpart.
Table 4.4 Interpolated stepwise and continuous pattern multipliers
Time Pattern Multiplier Stepwise Multiplier Continuous Multiplier
0:00 0.7 0.7 0.7
0:15 0.7 0.7 0.8
0:30 0.7 0.7 0.9
0:45 0.7 0.7 1.0
1:00 1.1 1.1 1.1
Pattern Start Time and Repetition. When defining and working with pat-
terns, it is important to understand how the pattern start time is referenced. Does pat-
tern hour 2 refer to 2:00 a.m., or does it refer to the second hour from the beginning of
162 Water Consumption Chapter 4
a simulation? If a model simulation begins at midnight, then there is no difference
between military time and time step number. If the model is intended to start at some
other time (such as 6:00 a.m., when many systems have refilled all their tanks), then
the patterns may need to be adjusted, advancing or retarding them in time accordingly.
Most modelers accept that demand patterns repeat every 24 hours with only negligible
differences, and are willing to use the same pattern each day in such a way that hours
25 and 49 use the same demands as the first hour. For a factory with three shifts, a pat-
tern may repeat every eight hours. Other patterns may not repeat at all. Each software
package handles pattern repetition in its own way; thus, some research and experi-
mentation may be required to produce the desired behavior for a particular applica-
tion.
4.4 PROJECTING FUTURE DEMANDS
Water distribution models are created not only to solve the problems of today, but also
to prevent problems in the future. With almost any endeavor, the future holds a lot of
uncertainty, and demand projection is no exception. Long-range planning may
include the analysis of a system for 5-, 10-, and 20-year time frames. When perform-
ing long-term planning analyses, estimating future demands is an important factor
influencing the quality of information provided by the model.
Section 4.4 Projecting Future Demands 163
The uncertainty of this process puts the modeler in the difficult position of trying to
predict the future. The complexity of such analyses, however, can be reduced to some
extent with software that supports the creation and comparison of a series of possible
alternative futures. Testing alternative future projections provides a way for the mod-
eler to understand the sensitivity of decisions regarding demand projections. Scenario
management tools in models help make this process easier. Even the most compre-
hensive scenario management, however, is just another tool that needs to be applied
intelligently to obtain reasonable results.
Historical Trends
Since the growth of cities and industries is hard to predict, it follows that it is also dif-
ficult to predict future water demands. Demand projections are only as accurate as the
assumptions made and the methods used to extrapolate development. Some cities
have relatively stagnant demands, but others experience volatile growth that chal-
lenges engineers designing water systems.
How will the economy affect local industries? Will growth rates continue at their cur-
rent rate, or will they level off? Will regulations requiring low-flow fixtures actually
result in a drop in water usage? What will be the combined result of increased popula-
tion and greater interest in water conservation? These questions are all difficult to
answer, and no method exists that can answer them with absolute certainty.
In general, the decision about which alternative future projection should be used is not
so much a modeling decision as a utility-wide planning decision. The modeler alone
should not try to predict the future, but rather facilitate the utility decision-makers
process of coming to a consensus on likely future demands.
Figure 4.13 illustrates some possible alternative futures given a historical demand pat-
tern. In spite of its shortcomings, the most commonly used method for predicting
demands is to examine historical demand trends and to extrapolate them into the
future under the assumption that they will continue.
Spatial Allocation of Future Demands
Planning departments and other groups may provide population projections for future
years and associate these population estimates with census tracts, traffic analysis
zones, planning districts, or other areas. The data must then be manipulated to deter-
mine the spatial allocation of nodal demands for the water model.
This manipulation requires a good deal of judgment on the modelers part, reflecting
the uncertainty of the process. Predictions concerning the future, by their nature, con-
tain varying degrees of uncertainty. If the significant factors affecting community
growth have been identified, the modeler can usually save time by making good judg-
ments about how the current baseline demand allocation can be modified and reused
for planning purposes. It is also important for the modeler to consider the future fire
protection requirements. Because fire protection demands are often much larger than
baseline demands, they are usually a major factor in future pipe-sizing decisions.
164 Water Consumption Chapter 4
Figure 4.13
Several methods for
projecting future
demands
Time, year
P
e
a
k
D
a
y
D
e
m
a
n
d
,
M
G
D
Annual Demand Data
2040 2030 2020 2010 2000 1990 1980 1970 1960
1.5
2.0
2.5
3.0
3.5
4.0
4.5
5.0
Constant Percent Growth
Growth to Buildout
Economic Downturn
Linear Growth
Disaggregated Projections
Rather than basing projections on extrapolation of flow rate data, it is somewhat more
rational to examine the causes of demand changes and then project that data into the
future. This technique is called disaggregated projection. Instead of predicting
demands, the user predicts such things as industrial production, number of hotel
rooms, and cost of water, and then uses a forecasting model to predict demand.
The simplest type of disaggregated demand projection involves projecting population
and per capita demand separately. In this way, the modeler can, for example, separate
the effects of population growth from the effects of a decrease in per capita consump-
tion due to low-volume fixtures and other water conservation measures.
These types of approaches attempt to account for many variables that influence future
demands, including population projections, water pricing, land use, industrial growth,
and the effects of water conservation (Vickers, 1991; and Macy, 1991). The IWR-
Main model (Opitz et al., 1998; Dziegielewski and Boland, 1989) is a sophisticated
model that uses highly disaggregated projections to forecast demands.
The most difficult factor to predict when performing a projection is drastic change in
the economy of an area (for example, a military base closure or the construction of a
factory). Using disaggregated projections, population projections can be modified
Section 4.5 Fire Protection Demands 165
more rationally than can flow projections when developing demand forecasts that
reflect these types of events.
Population Estimates. Planning commissions often have population studies
and estimates that predict the future growth of a city or town. Though population esti-
mates usually contain uncertainties, they can be used as a common starting point for
any model requiring future estimates, such as water distribution models, sewer plans,
and traffic models.
Starting with current per capita usage rates or projections of per capita usage trends,
future demands can be estimated by taking the product of the future population and
the future per capita usage. In areas that are already densely populated, the growth
may be only slightly positive, or even negative.
The United States Geological Survey (USGS) publishes per capita water consumption
rates for each state, but these values include nonmunicipal uses such as power genera-
tion and agriculture. A per capita consumption rate developed in this manner cannot
be widely applied because there are large differences in water consumption among
customers in different areas within a particular state.
Land Use. Sometimes, water demands can be estimated based on land use desig-
nations such as single-family residential, high-density residential, commercial, light
industrial, heavy industrial, and so on. Information regarding a representative water
usage rate based on land use can then aid in planning for other areas that are in the
same category.
As with population estimates, using land use designation requires some level of pre-
diction regarding future growth in every area from residential land use to industrial
and commercial operations. For example, the loss or gain of a single large industry
can have a tremendous effect on the overall consumption in the system.
4.5 FIRE PROTECTION DEMANDS
When a fire is in progress, fire protection demands can represent a huge fraction of
the total demand for the system. The effects of fire demands are difficult to derive pre-
cisely since fires occur with random frequency in different areas, with each area hav-
ing unique fire protection requirements. Generally, the amount of water needed to
adequately fight a fire depends on the size of the burning structure, its construction
materials, the combustibility of its contents, and the proximity of adjacent buildings.
For some systems, fire protection is a lower priority than water quality or construction
costs. To reduce costs in situations in which customers are very spread out, such as in
rural areas, the network may not be designed to provide fire protection. Instead, the
fire departments rely on water tanker trucks or other sources for water to combat fires
(for example, ponds constructed specifically for that purpose).
One of the primary benefits of providing water for fire protection is a reduction in the
insurance rates of residents and businesses in the community. In the United States,
community fire protection infrastructure (the fire-fighting capabilities of the fire
department and the capacity of the water distribution network) is audited and rated by
the Insurance Services Office (ISO) using the Fire Protection Rating System (ISO,
166 Water Consumption Chapter 4
1998). In Canada, the Insurers Advisory Organization (IAO) evaluates water supply
systems using the Grading Schedule for Municipal Fire Protection (IAO, 1974). The
ISO evaluation process is summarized in AWWA M-31 (1998).
In Europe, no fire prevention standards exist that apply to all European countries;
therefore, each country must develop or adopt its own fire flow requirements. For
example, the flow rates that the UK Fire Services ideally require to fight fires are
based on the national guidance document on the provision of water for fire fighting
(Water UK and LGA, 1998). Similarly, the German standards (DVGW, 1978), the
French standards (Circulaire, 1951, 1957, and 1967), the Russian standards (SNIP,
1985), and others, are based on fire risk categories that assign the level of risk accord-
ing to the type of premises to be protected, fire-spread risk, installed fire proofing, or
any combination of these factors.
Because systems will be evaluated using ISO methods, engineers in the United States
usually base design of fire protection systems on the ISO rating system, which
includes determining fire flow demands according to the ISO approach. Although the
actual water needed to fight a fire depends on the structure and the fire itself, the ISO
method yields a Needed Fire Flow (NFF) that can be used for design and evaluation
of the system. Different calculation methods are used for different building types,
such as residential, commercial, or industrial.
For one- and two-family residences, the needed fire flow is determined based on the
distance between structures, as shown in Table 4.5.
Table 4.5 Needed fire flow for residences two stories and less
Distance Between Buildings
(ft)
Fire Flow
(gpm)
More than 100 500
31-100 750
11-30 1,000
Less than 11 1,500
For commercial and industrial structures, the needed fire flow is based on building
area, construction class (that is, frame or masonry construction), occupancy (such as
a department store or chemical manufacturing plant), exposure (distance to and type
of nearest building), and communication (types and locations of doors and walls). The
formula can be summarized as:
NFF 18FA
0.5
O X P + ( ) = (4.12)
where NFF = needed fire flow (gpm)
F = class of construction coefficient
A = effective area (ft
2
)
O = occupancy factor
X = exposure factor
P = communication factor
References 167
The procedure for determining NFF is documented in the Fire Protection Rating Sys-
tem (1998) and AWWA M-31 (1998). The minimum needed fire flow is not less than
500 gpm (32 l/s), and the maximum is no more than 12,000 gpm (757 l/s). Most fre-
quently, the procedure produces values less than 3,500 gpm (221 l/s). Values are
rounded to the nearest 250 gpm (16 l/s) for NFFs less than 2,500 gpm (158 l/s), and to
the nearest 500 gpm (32 l/s) for values greater than 2,500 gpm (158 l/s). Values are
also adjusted if a building is equipped with sprinklers.
In addition to a flow rate requirement, a requirement exists for the duration over
which the flow can be supplied. According to ISO (1998), fires requiring 3,500 gpm
(221 l/s) or less are referred to as receiving Public Fire Suppression, and those
requiring greater than 3,500 gpm (221 l/s) are classified as receiving Individual
Property Fire Suppression. For fires requiring 2,500 gpm (158 l/s) or less, a two-hour
duration is sufficient; for fires needing 3,000 to 3,500 gpm (190 to 221 l/s), a three-
hour duration is used; and for fires needing more than 3,500 gpm (221 l/s), a four-
hour duration is used along with slightly different rules for evaluation.
Methods for estimating sprinkler demands are based on the area covered and a flow
density in gpm/ft
2
as described in NFPA 13 (1999) for commercial and industrial
structures, and in NFPA 13D (1999) for single- and two-family residential dwellings.
For residences, the sprinklers shall provide at least 18 gpm (1.14 l/s) when one sprin-
kler operates and no less than 13 gpm (0.82 l/s) per sprinkler when more than one
operates. For commercial and industrial buildings, the flow density can vary from
0.05 to 0.35 gpm/ft
2
(2 to 14 l/min/m
2
), depending on the hazard class associated with
the building and the floor area. Sprinkler design is covered in detail on page 401.
NFPA 13 provides a chart for determining flow density based on whether occupancy
is light, ordinary hazard, or extra hazard. A hose stream requirement exists as well for
water used to supplement the sprinkler flows. These values range from 100 to 1000
gpm (6.3 to 63 l/s), depending on the hazard classification.
REFERENCES
American Water Works Association (1989). Distribution Network Analysis for Water Utilities. AWWA
Manual M-32, Denver, Colorado.
American Water Works Association (1998). Distribution System Requirements for Fire Protection.
AWWA Manual M-31, Denver, Colorado.
Basford, C., and Sevier, C. (1995). Automating the Maintenance of Hydraulic Network Model Demand
Database Utilizing GIS and Customer Billing Records. Proceedings of the AWWA Computer Confer-
ence, American Water Works Association, Norfolk, Virginia.
Bowen, P. T., Harp, J., Baxter, J., and Shull, R. (1993). Residential Water Use Patterns. AWWARF, Denver,
Colorado.
Brainard, B. (1994). Using Electronic Rate of Flow Recorders. Proceedings of the AWWA Distribution
System Symposium, American Water Works Association, Omaha, Nebraska.
Buchberger, S. G., and Wu, L. (1995). A Model for Instantaneous Residential Water Demands. Journal of
Hydraulic Engineering, ASCE, 121(3), 232.
Buchberger, S. G., and Wells, G. J. (1996). Intensity, Duration, and Frequency of Residential Water
Demands. Journal of Water Resources Planning and Management, ASCE, 122(1), 11.
168 Water Consumption Chapter 4
Buyens, D. J., Bizier, P. A., and Combee, C. W. (1996). Using a Geographical Information System to
Determine Water Distribution Model Demands. Proceedings of the AWWA Annual Conference, Ameri-
can Water Works Association, Toronto, Canada.
Cesario, A. L., and Lee T. K. (1980). A Computer Method for Loading Model Networks. Journal of the
American Water Works Association, 72(4), 208.
Cesario, A. L. (1995). Modeling, Analysis, and Design of Water Distribution Systems. American Water
Works Association, Denver, Colorado.
Circulaire des Ministreres de lIntriur et de lAgriculture du Fvrier (1957). Protection contre lincendie
dans les communes rurales. Paris, France.
Circulaire du Ministrere de lAgriculture du Auout (1967). Rserve deau potable. Protection contre
lincendie dans les communes rurales. Paris, France.
Circulaire Interministrielle du Dcembre (1951). Alimentation des communes en eau potable - Lutte con-
tre lincendie. Paris, France.
Coote, P. A., and Johnson, T. J. (1995). Hydraulic Model for the Mid-Size Utility. Proceedings of the
AWWA Computer Conference, American Water Works Association, Norfolk, Virginia.
Davis, A. L., and Brawn, R. C. (2000). General Purpose Demand Allocator (DALLOC). Proceedings of
the Environmental and Water Resources Institute Conference, American Society of Civil Engineers,
Minneapolis, Minnesota.
DeOreo, W. B., Heaney, J. P., and Mayer, P. W. (1996). Flow Trace Analysis to Assess Water Use. Jour-
nal of the American Water Works Association, 88(1), 79.
DVGW. (1978). DVGW W405 Bereitstellung von Lschwasser durch die ffentliche Trinkwasserversor-
gung. Deutscher Verein des Gas und Wasserfaches, Franfurt, Germany.
Dziegielewski, B., and Boland J. J. (1989). Forecasting Urban Water Use: the IWR-MAIN Model. Water
Resource Bulletin, 25(1), 101 119.
Dziegielewski, B., Opitz, E. M., and Maidment, D. (1996). Water Demand Analysis. Water Resources
Handbook, Mays, L. W., ed., McGraw-Hill, New York, New York.
Insurance Advisory Organization (IAO) (1974). Grading Schedule for Municipal Fire Protection. Toronto,
Canada.
Insurance Services Office (ISO) (1998). Fire Suppression Rating Schedule. New York, New York.
Macy, P. P. (1991). Integrating Construction and Water Master Planning. Journal of the American Water
Works Association, 83(10), 44 47.
Male, J. W., and Walski, T. M. (1990). Water Distribution: A Troubleshooting Manual. Lewis Publishers,
Chelsea, Florida.
Metcalf & Eddy, Inc. (1979). Water Resources and Environmental Engineering. 2nd Edition, McGraw-Hill,
New York, New York.
National Fire Protection Association (NFPA) (1999). Sprinkler Systems in One- and Two-Family Dwell-
ings and Manufactured Homes. NFPA 13D, Quincy, Massachusetts.
National Fire Protection Association (NFPA) (1999). Standard for Installation of Sprinkler Systems.
NFPA 13, Quincy, Massachusetts.
Office of Water Services (Ofwat) (1998). 1997-98 Report on Leakage and Water Efficiency. http://
www.open.gov.uk/ofwat/leak97.pdf, United Kingdom.
Opitz, E. M., et al. (1998). Forecasting Urban Water Use: Models and Application. Urban Water Demand
Management and Planning, Baumann D., Boland, J. and Hanemann, W. H., eds., McGraw Hill. New
York, New York, 350.
Rhoades, S. D. (1995). Hourly Monitoring of Single-Family Residential Areas. Journal of the American
Water Works Association, 87(8), 43.
References 169
SNIP (1985). Water Supply Standards (in Russian). 2.04.02-84, Moscow, Russia.
Stern, C. T. (1995). The Los Angeles Department of Water and Power Hydraulic Modeling Project. Pro-
ceedings of the AWWA Computer Conference, American Water Works Association, Norfolk, Virginia.
Vickers, A. L. (1991). The Emerging Demand Side Era in Water Conservation. Journal of the American
Water Works Association, 83(10), 38.
Walski, T. M. (1999). Peaking Factors for Systems with Leakage. Essential Hydraulics and Hydrology,
Haestad Press, Waterbury, Connecticut.
Walski, T. M., Lowry, S. G., and Rhee, H. (2000). Pitfalls in Calibrating an EPS Model. Proceedings of
the Environmental and Water Resource Institute Conference, American Society of Civil Engineers, Min-
neapolis, Minnesota.
Water Research Centre (WRc) (1985). District Metering, Part I - System Design and Installation. Report
ER180E, United Kingdom.
Water UK and Local Government Association (1998). National Guidance Document on the Provision of
Water for Fire Fighting. London, United Kingdom.
Ysuni, M. A. (2000). System Design: An Overview. Water Distribution Systems Handbook, Mays, L. W.,
ed., McGraw-Hill, New York, New York.
170 Water Consumption Chapter 4
DISCUSSION TOPICS AND PROBLEMS
Read the chapter and complete the problems. Submit your work to Haestad Methods
and earn up to 11.0 CEUs. See Continuing Education Units on page xxix for more
information, or visit www.haestad.com/awdm-ceus/.
4.1 Develop a steady-state model of the water distribution system shown in the figure. Data describing
the system and average daily demands are provided in the tables that follow.
P -1
P -7 P -6 P -5
P -4 P -3
P -2
P -9
P -8
J-1
J-5
J-4
J-3
J-2
R-1
Pipe Label
Length
(ft)
Diameter
(in.)
Hazen-
Williams
C-factor
Minor Loss
Coefficient
P-1 500 12 120 10
P-2 2,600 10 120 0
P-3 860 8 120 0
P-4 840 8 120 5
P-5 710 6 120 0
P-6 1,110 4 120 0
P-7 1,110 4 120 0
P-8 710 6 120 0
P-9 1,700 6 120 0
Node Label
Elevation
(ft)
Demand
(gpm)
R-1 750 N/A
J-1 550 250
J-2 520 75
J-3 580 125
J-4 590 50
J-5 595 0
Discussion Topics and Problems 171
a) Fill in the tables below with the pipe and junction node results.
Pipe Label
Flow
(gpm)
Hydraulic Gradient
(ft/1000 ft)
P-1
P-2
P-3
P-4
P-5
P-6
P-7
P-8
P-9
Node Label
Hydraulic Grade
(ft)
Pressure
(psi)
J-1
J-2
J-3
J-4
J-5
172 Water Consumption Chapter 4
b) Complete the tables below assuming that all demands are increased to 225 percent of average-
day demands.
Pipe Label
Flow
(gpm)
Hydraulic Gradient
(ft/1000 ft)
P-1
P-2
P-3
P-4
P-5
P-6
P-7
P-8
P-9
Node Label
Hydraulic Grade
(ft)
Pressure
(psi)
J-1
J-2
J-3
J-4
J-5
c) Complete the tables below assuming that, in addition to average-day demands, there is a fire flow
demand of 1,850 gpm added at node J-3.
Pipe Label
Flow
(gpm)
Hydraulic Gradient
(ft/1000 ft)
P-1
P-2
P-3
P-4
P-5
P-6
P-7
P-8
P-9
Node Label
Hydraulic Grade
(ft)
Pressure
(psi)
J-1
J-2
J-3
J-4
J-5
Discussion Topics and Problems 173
4.2 English Units: Perform a 24-hour extended-period simulation with a one-hour time step for the sys-
tem shown in the figure. Data necessary to conduct the simulation are provided in the tables that fol-
low. Alternatively, the pipe and junction node data has already been entered into Prob4-02.wcd. Use
a stepwise format for the diurnal demand pattern. Answer the questions presented at the end of this
problem.
J-1
J-9
J-8
J-6
J-5
J-4
J-3 J-2
J-7
J-10
Suction
P-5
P-1
P-12 P-13
P-4
P-2
P-9
P-8
P-7 P-6
Discharge
P-10 P-11
P-14
P-15
P-16
P-3
Crystal Lake
Miamisburg
Tank
West Carrolton
Tank
(Not To Scale)
Pipe Label
Length
(ft)
Diameter
(in.)
Hazen-Williams
C-factor
Suction 25 24 120
Discharge 220 21 120
P-1 1,250 6 110
P-2 835 6 110
P-3 550 8 130
P-4 1,010 6 110
P-5 425 8 130
P-6 990 8 125
P-7 2,100 8 105
P-8 560 6 110
P-9 745 8 100
P-10 1,100 10 115
P-11 1,330 8 110
P-12 890 10 115
P-13 825 10 115
P-14 450 6 120
P-15 690 6 120
P-16 500 6 120
Node Label
Elevation
(ft)
Demand
(gpm)
Crystal Lake 320 N/A
J-1 390 120
J-2 420 75
J-3 425 35
J-4 430 50
J-5 450 0
J-6 445 155
J-7 420 65
J-8 415 0
J-9 420 55
J-10 420 20
Pump Curve Data
Head
(ft)
Flow
(gpm)
Shutoff 245 0
Design 230 1,100
Max Operating 210 1,600
Elevated Tank Information
Miamisburg
Tank
West Carrolton
Tank
Base Elevation (ft) 0 0
Minimum Elevation (ft) 535 525
Initial Elevation (ft) 550 545
Maximum Elevation (ft) 570 565
Tank Diameter (ft) 49.3 35.7
Diurnal Demand Pattern
Time of Day Multiplication Factor
Midnight 1.00
6:00 am 0.75
Noon 1.00
6:00 pm 1.20
Midnight 1.00
174 Water Consumption Chapter 4
a) Produce a plot of the HGL in the Miamisburg and West Carrolton tanks as a function of time.
b) Produce a plot of the pressures at node J-3 versus time.
Discussion Topics and Problems 175
SI Units: Perform a 24-hour extended-period simulation with a one-hour time step for the system
shown in the figure. Data necessary to conduct the simulation are provided in the tables below. Alter-
natively, the pipe and junction node data has already been entered into Prob4-02m.wcd. Use a step-
wise format for the diurnal demand pattern. Answer the questions presented at the end of this
problem.
Pipe Label
Length
(m)
Diameter
(mm)
Hazen-Williams
C-factor
Suction 7.6 610 120
Discharge 67.1 533 120
P-1 381.0 152 110
P-2 254.5 152 110
P-3 167.6 203 130
P-4 307.8 152 110
P-5 129.5 203 130
P-6 301.8 203 125
P-7 640.1 203 105
P-8 170.7 152 110
P-9 227.1 203 100
P-10 335.3 254 115
P-11 405.4 203 110
P-12 271.3 254 115
P-13 251.5 254 115
P-14 137.2 152 120
P-15 210.3 152 120
P-16 152.4 152 120
Node Label
Elevation
(m)
Demand
(l/s)
Crystal Lake 97.5 N/A
J-1 118.9 7.6
J-2 128.0 4.7
J-3 129.5 2.2
J-4 131.1 3.2
J-5 137.2 0
J-6 135.6 9.8
J-7 128.0 4.1
J-8 126.5 0
J-9 128.0 3.5
J-10 128.0 1.3
Pump Curve Data
Head
(m)
Flow
(l/s)
Shutoff 74.6 0
Design 70.1 69
Max Operating 64.0 101
Elevated Tank Information
Miamisburg
Tank
West Carrolton
Tank
Base Elevation (m) 0 0
Minimum Elevation (m) 163.1 160.0
Initial Elevation (m) 167.6 166.1
Maximum Elevation (m) 173.7 172.2
Tank Diameter (m) 15.0 10.9
Diurnal Demand Pattern
Time of
Day
Multiplication
Factor
Midnight 1.00
6:00 a.m. 0.75
Noon 1.00
6:00 p.m. 1.20
Midnight 1.00
176 Water Consumption Chapter 4
a) Produce a plot of the HGL in the Miamisburg and West Carrolton tanks as a function of time.
b) Produce a plot of the pressures at node J-3 versus time.
4.3 Develop a steady-state model for the system shown in the figure and answer the questions that fol-
low. Data necessary to conduct the simulation are provided in the following tables. Alternatively, the
pipe and junction node data has already been entered into Prob4-03.wcd. Note that there are no
minor losses in this system. The PRV setting is 74 psi.
Central Tank
Newtown Res.
High Field Res.
Lower Pressure Zone
Pressure Zone Boundary
Higher Pressure Zone
J-10
P-16
J-9
P-17
J-11 P-18
J-3
P-19
P-15 J-8
P-5 P-4
PRV-1 P-14
J-2 P-3
P-2
P-1
PMP-1
J-1
P-9
J-7
P-13
J-6 P-8
P-12
P-20
P-11
J-5
PMP-2
P-10
P-7
J-4
P-6
Pipe Label
Length
(ft)
Diameter
(in.)
Hazen-Williams
C-factor
P-1 120 24 120
P-2 435 16 120
P-3 2,300 12 120
P-4 600 10 110
P-5 550 10 110
P-6 1,250 12 110
P-7 850 12 110
P-8 4,250 12 120
P-9 2,100 12 120
P-10 50 24 105
P-11 250 16 105
P-12 1,650 10 115
P-13 835 8 110
P-14 800 8 100
P-15 1,300 6 95
P-16 1,230 6 95
P-17 750 6 95
P-18 1,225 8 95
P-19 725 6 100
P-20 155 4 75
Node Label
Elevation
(ft)
Demand
(gpm)
High Field Reservoir 1,230 N/A
Newtown Reservoir 1,050 N/A
Central Tank 1,525 N/A
J-1 1,230 0
J-2 1,275 0
J-3 1,235 120
J-4 1,250 35
J-5 1,300 55
J-6 1,250 325
J-7 1,260 0
J-8 1,220 100
J-9 1,210 25
J-10 1,210 30
J-11 1,220 45
PRV-1 1,180 N/A
PMP-1 1,045 N/A
PMP-2 1,225 N/A
Discussion Topics and Problems 177
Pump Curve Data
PMP-1 PMP-2
Head
(ft)
Flow
(gpm)
Head
(ft)
Flow
(gpm)
Shutoff 550 0 320 0
Design 525 750 305 1,250
Max Operating 480 1,650 275 2,600
178 Water Consumption Chapter 4
a) Fill in the tables for pipe and junction node results.
Pipe Label
Flow
(gpm)
Hydraulic Gradient
(ft/1000 ft)
P-1
P-2
P-3
P-4
P-5
P-6
P-7
P-8
P-9
P-10
P-11
P-12
P-13
P-14
P-15
P-16
P-17
P-18
P-19
P-20
Node Label
Hydraulic Grade
(ft)
Pressure
(psi)
J-1
J-2
J-3
J-4
J-5
J-6
J-7
J-8
J-9
J-10
J-11
Discussion Topics and Problems 179
Analyze the following demand conditions for this system by using the average-day demands as your
base demands.
b) Increase all demands to 150 percent of average-day demands. What are the pressures at nodes J-2
and J-10?
c) Add a fire flow demand of 1,200 gpm to node J-4. What is the discharge from the Newtown
pump station? What is the pressure at node J-4?
d) Replace the demand of 120 gpm at node J-3 with a demand of 225 gpm. How does the pressure at
node J-3 change between the two demand cases?
e) Replace the existing demands at nodes J-3, J-9, J-10, and J-11 with 200 gpm, 50 gpm, 90 gpm,
and 75 gpm, respectively. Is Central Tank filling or draining? How does the tank condition com-
pare with the original simulation before demands were changed?
4.4 Perform an extended-period simulation on the system from part (a) of Problem 4.3. However, first
add a PRV to pipe P-6 and close pipe P-14. Note that pipe P-6 must split into two pipes when the
PRV is inserted. Specify the elevation of the PRV as 1,180 ft and the setting as 74 psi.
The simulation duration is 24 hours and starts at midnight. The hydraulic time step is 1 hour. The
capacity and geometry of the elevated storage tank and the diurnal demand pattern are provided
below. Assume that the diurnal demand pattern applies to each junction node and that the demand
pattern follows a continuous format. Assume that the High Field pump station does not operate.
Central Tank Information
Base Elevation (ft) 1,260
Minimum Elevation (ft) 1,505
Initial Elevation (ft) 1,525
Maximum Elevation (ft) 1,545
Tank Diameter (ft) 46.1
Diurnal Demand Pattern
Time of Day Multiplication Factor
Midnight 0.60
3:00 a.m. 0.75
6:00 a.m. 1.20
9:00 a.m. 1.10
Noon 1.15
3:00 p.m. 1.20
6:00 p.m. 1.33
9:00 p.m. 0.80
Midnight 0.60
a) Produce a plot of HGL versus time for Central Tank.
b) Produce a plot of the discharge from the Newtown pump station versus time.
c) Produce a plot of the pressure at node J-3 versus time.
d) Does Central Tank fill completely? If so, at what time does the tank completely fill? What hap-
pens to a tank when it becomes completely full or completely empty?
180 Water Consumption Chapter 4
e) Why does the discharge from the Newtown pump station increase between midnight and 6:00
a.m.? Why does the discharge from the pump station decrease, particularly after 3:00 p.m.?
f) Does the pressure at node J-3 vary significantly over time?
4.5 Given a pressure zone with one pump station pumping into it and a smaller one pumping out of it,
and a single 40-ft diameter cylindrical tank, develop a diurnal demand pattern. The pumping rates
and tank water levels are given in the table below. The pumping rates are the average rates during the
hour, and the tank levels are the values at the beginning of the hour.
What is the average use in this pressure zone? What is the average flow to the higher pressure zone?
Hour
Pump In
(gpm)
Pump Out
(gpm)
Tank Level
(ft)
0 650 0 35.2
1 645 210 38.5
2 645 255 40.4
3 652 255 42.1
4 310 255 43.5
5 0 255 42.8
6 0 255 39.6
7 0 0 36.0
8 0 0 33.4
9 225 0 30.3
10 650 0 28.9
11 650 0 30.5
12 650 0 32.1
13 650 0 33.8
14 650 45 35.8
15 645 265 37.5
16 645 260 37.5
17 645 260 37.2
18 645 260 36.5
19 645 260 36.4
20 645 255 36.7
21 645 150 37.2
22 115 0 38.7
23 0 0 38.3
24 0 0 37.1




ADVANCED WATER
DISTRIBUTION MODELING
AND MANAGEMENT

Authors
Thomas M. Walski
Donald V. Chase
Dragan A. Savic
Walter Grayman
Stephen Beckwith
Edmundo Koelle

Contributing Authors
Scott Cattran, Rick Hammond, Kevin Laptos, Steven G. Lowry,
Robert F. Mankowski, Stan Plante, John Przybyla, Barbara Schmitz

Peer Review Board
Lee Cesario (Denver Water), Robert M. Clark (U.S. EPA),
Jack Dangermond (ESRI), Allen L. Davis (CH2M Hill),
Paul DeBarry (Borton-Lawson), Frank DeFazio (Franklin G. DeFazio Corp.),
Kevin Finnan (Bristol Babcock), Wayne Hartell (Bentley Systems),
Brian Hoefer (ESRI), Bassam Kassab (Santa Clara Valley Water District),
James W. Male (University of Portland), William M. Richards
(WMR Engineering), Zheng Wu (Bentley Systems ),
and E. Benjamin Wylie (University of Michigan)



Click here to visit the Bentley Institute
Press Web page for more information

C H A P T E R
5
Testing Water Distribution Systems
Verifying that a water distribution model replicates field conditions requires an inti-
mate knowledge of how the system performs over a wide range of operating condi-
tions. For example, can the model reproduce the flow patterns and pressures that
occur during periods of peak summertime usage, or can the model accurately simulate
chlorine decay? Collecting water distribution system data in the field provides valu-
able insight into system performance and is an essential part of calibration.
Data collection, the first step in the model calibration process, is discussed in depth
within this chapter. The chapter begins with a brief discussion of system testing,
including descriptions of some simple tests for measuring flow and pressure, as well
as some of the pitfalls that may be encountered. The details of performing fire hydrant
flow tests, head loss tests, pump performance tests, and water quality tests are dis-
cussed as well. The chapter concludes with a discussion of the importance of data
quality, particularly when automated calibration methods are used.
5.1 TESTING FUNDAMENTALS
Pressure Measurement
Pressures are measured throughout the water distribution system to monitor the level
of service and to collect data for use in model calibration. Pressure readings are com-
monly taken at fire hydrants (see Figure 5.1) but can also be read at hose bibs (also
called spigots); home faucets; pump stations (both suction and discharge sides);
tanks; reservoirs; and blow-off, air release, and other types of valves.
If the measurements are taken at a location other than a direct connection to a water
main (for example, at a house hose bib), the head loss between the supply main and
the site where pressure is measured must be considered. Of course, the best solution is
to have no flow (and hence no head loss) between the main and the gage. To check if
flow into the building is occurring, listen at the hose bib for the sound of rushing water.
182 Testing Water Distribution Systems Chapter 5
Figure 5.1
Pressure gages on a
fire hydrant
When measuring pressure, slight fluctuations may be seen on the gage due to chang-
ing flows in the system. Devices such as pressure snubbers and liquid-filled pressure
gages can be used to dampen the pressure fluctuations, unless the fluctuations them-
selves are a source of interest.
Pressure gages are most accurate when measuring pressures within 50 to 75 percent
of the maximum value on the scale. Using several pressure gages of varying pressure
ranges is advisable when working with a water distribution system. A pressure gage
with a range of 0 to 100 psi (690 kPa) is commonly used; however, a pressure gage
that can read up to 200 psi (1,380 kPa) may be necessary for measurements taken at a
pump discharge or at a low elevation. If pressure measurements are taken on the suc-
tion side of a pump, then a pressure gage capable of reading negative pressures, called
a pressure-vacuum gage, may be required. Remember that it is the elevation of the
gage, not the elevation of the node, that is used in calculating the elevation of the HGL
(see page 252).
Flow Measurement
Flow is measured at key locations throughout a system to provide insight into flow
patterns and system performance, develop consumption data, and determine flow
rates for calibration.
Many of the tests described in this chapter require measuring flow in pipes. A variety
of flow meters are available for this purpose, including Venturi meters, magnetic flow-
meters, and ultrasonic meters. Pressure and flow metering and recording equipment
should be calibrated regularly and undergo routine performance checks to ensure that
it is in good working order. Furthermore, even if a flow meter is accurate and cali-
brated, the monitoring station may use an analog gage or dial readout that has a coarse
level of precision, which limits the overall precision.
Section 5.1 Testing Fundamentals 183
The extent of flow measurement employed varies from system to system. Usually,
flow is measured continuously at only a few key locations in the distribution system
such as treatment plants and pump stations. Data from these sites should be used to
the greatest extent possible in system calibration. Flow from higher to lower pressure
zones can also be measured at pressure zone boundaries using combination pressure
reducing valve/flow meters (Walski, Gangemi, Kaufman, and Malos, 2001). More
rarely, systems employ in-line flow meters at key points throughout the network and
transmit the flow rates back to a control center using Supervisory Control and Data
Acquisition (SCADA) systems and telemetry (See Chapter 6). This type of comprehen-
sive flow monitoring is not typically done in the United States; however, more utility
managers and operators are starting to see the value of in-line flow information.
Temporary flow metering may be a cost-effective option to check pump discharges or
to see if in-line flow measurements are required throughout the system. Field mea-
surement using a Pitot rod is shown in Figure 5.2. The rod is inserted into the pipe to
measure total head and pressure head, which can then be converted into velocity
(Walski, 1984a). The Pitot rod should not be confused with the Pitot gage, which
measures velocity head only. Clamp-on or insertion electromagnetic or ultrasonic
meters may also be used.
Placement of the flow-measuring device is important. To be sure that disturbances
caused by any bends or obstructions do not influence the readings, the device should
be placed far enough downstream of the disturbance (usually at a distance of approxi-
mately 10 times the pipe diameter) that the effects will have completely dissipated.
In certain cases it may be desirable to isolate one end of the pipe such that all of the
flow through the pipe is diverted through a hydrant for measurement. The hydrant
flow can then be measured with a hydrant Pitot gage as described in Section 5.2.
Net flow in and out of a tank during a time period can be measured by monitoring
water level in the tank and then calculating the flow based on cross-sectional area in
the tank.
184 Testing Water Distribution Systems Chapter 5
Figure 5.2
Tip of Pitot rod
inserted into clear
pipe
Potential Pitfalls in System Measurements
Flow measurement tests can be beneficial, but there are potential drawbacks to keep
in mind. Testing may result in disruption of service to some customers. For example,
fire flow tests typically cause lower than normal pressures and higher than normal
velocities, particularly in residential areas. Higher velocities can entrain sediments in
pipes or shear against tuberculation on pipe walls, causing customers to experience
discolored water.
Customers may, either by accident or necessity, be disconnected from the system
when valves are operated to facilitate flow tests. As described in the following sec-
tions, head loss tests require the operation of system valves to isolate sections of water
main. Valve operation needs to be carefully planned when conducting such tests to
avoid inadvertently disconnecting customers from the system. To avoid surprises, cus-
tomers should be notified prior to the tests.
5.2 FIRE HYDRANT FLOW TESTS
Obtaining data for a wide range of operating conditions, including peak (high)
demand periods, would be difficult without fire hydrant flow tests. These tests can be
used to simulate high flow conditions (see page 218) and allow the system behavior to
be analyzed under extreme conditions. Fire hydrant flow tests are primarily used to
measure the fire flow capacity of the system. They also provide data on pressures
within the system under static conditions (no hydrants flowing) and stressed condi-
tions (high flows occurring at the hydrants) and can be used in conjunction with the
hydraulic model to calibrate parameters such as pipe roughness (Walski, 1988). Pro-
cedures for conducting fire hydrant flow tests are described in AWWA (1989) and
ISO (1963).
Section 5.2 Fire Hydrant Flow Tests 185
Two or more hydrants are required to perform a fire hydrant flow test, as illustrated in
Figure 5.3. One hydrant is identified as the residual hydrant(s), where all pressure
measurements are taken, and the other is identified as the flowed hydrant(s), where all
flow measurements are taken. When the flowed hydrant(s) is closed, referred to as
static conditions, the pressure at the residual hydrant is called the static pressure.
When one or more of the flowed hydrants are open, referred to as flowed conditions,
the pressure at the residual hydrant is called the residual pressure.
Figure 5.3
Hydrant flow test
Residual
Hydrant
Flowed Hydrants
Pressure
Gage
Q
2
Q
1
Q
3
//=//= //=//= //=//=
Conducting a fire hydrant flow test is a simple procedure, and a number of these tests
can be conducted throughout the system in a days time. Although not essential, many
utilities have a policy requiring that the residual hydrant be opened and allowed to
flow prior to connecting the pressure gage. This precaution helps remove any particles
that have accumulated in the hydrant lateral and barrel since it was last exercised.
After that, a pressure gage is connected to the residual hydrant and a static pressure
reading taken.
Next, the first of the flowed hydrants is opened and flowed. Once the readings stabi-
lize, a reading is taken at the flowed hydrant using a hand-held or clamp-on Pitot gage
(shown in Figure 5.4) or a Pitot diffuser (shown in Figure 5.5). Meanwhile, another
pressure reading is taken at the residual hydrant. Once the residual pressure is taken
and the discharge rate of the flowed hydrant is recorded, the same procedure can be
repeated for additional hydrants if needed.
The number of hydrants that should be flowed during a test is determined by the pres-
sure drop observed at the residual hydrant. Usually, a drop of at least 10 psi (70 kPa)
is needed to give good results. In a 6- to 8-in. pipe (150 to 200 mm), flowing a single
hydrant is sufficient. For larger pipes, more hydrants may need to be flowed.
Pitot Gages and Diffusers
Because a Pitot gage (shown in Figure 5.4) converts virtually all of the velocity head
associated with the flow stream to pressure head, the Pitot gage pressure reading can
be converted to a hydrant discharge rate using the orifice relationship in Equation 5.1.
Q C
f
C
d
D
2
P = (5.1)
186 Testing Water Distribution Systems Chapter 5
Figure 5.4
Hand-held Pitot gage
where Q = hydrant discharge (gpm, l/s)
C
d
= discharge coefficient
D = outlet diameter (in., cm)
P = pressure reading from Pitot gage (psi, kPa)
C
f
= unit conversion factor (29.8 English, 0.111 SI)
For a typical 2.5-in. (64 mm) outlet with a discharge coefficient of 0.9, Equation 5.1
can be reduced to:
Q 167 P =
The discharge coefficient in Equation 5.1 accounts for the decrease in the diameter of
flow that occurs between the hydrant opening and the end of the Pitot gage, as well as
the head losses through the opening. The coefficient depends on the geometry of the
inside of the hydrant opening and can be determined by feeling the inside of the
hydrant nozzle (see Figure 5.6).
The Pitot diffuser is similar to a Pitot gage except that it incorporates a nozzle that
redirects the flow from the hydrant, reducing its momentum and thus the potential for
erosion. Because the velocity head sensor is measuring inside the diffuser at a point
where the pressure is not equal to zero, a slightly modified formula is required to
compute flow. This formula varies with the manufacturer of the diffuser (Walski and
Lutes, 1990; and Morin and Rajaratnam, 2000). For example, for the Pitot diffuser
shown in Figure 5.5, the coefficient of 167 given previously reduces to 140.
Section 5.2 Fire Hydrant Flow Tests 187
Figure 5.5
Pitot diffuser
Figure 5.6
Discharge coefficients
at hydrant openings
Rounded
C = 0.9
d
Square and Sharp
C = 0.8
d
Projecting
C = 0.7
d
To briefly review, the procedure for conducting a fire hydrant flow test is as follows:
1. Place a pressure gage on the residual hydrant and record static pressure.
2. Take the 2 -in. (64 mm) cap off of the flowed hydrant.
3. Feel the inside of the hydrant opening to determine its geometry.
4. Slowly start the flow.
5. Once readings stabilize, take a Pitot gage reading at the flowed hydrant(s).
6. Simultaneously measure the residual pressure(s) at the residual hydrant(s).
7. Slowly close the hydrants.
188 Testing Water Distribution Systems Chapter 5
8. Assign the discharge coefficient according to the geometry of the hydrant open-
ing.
9. Determine the hydrant discharge rate by using Equation 5.1 or the equation pro-
vided by the Pitot diffuser manufacturer.
Once all of the data have been collected, a table similar to Table 5.1 can be con-
structed to present the results of the fire hydrant flow test.
Table 5.1 Results of fire hydrant flow test
Number of
Hydrants Flowing
Residual
Pressure
(psi)
Hydrant #1
Discharge
(gpm)
Hydrant #2
Discharge
(gpm)
Hydrant #3
Discharge
(gpm)
Total
Discharge
(gpm)
0 78 N/A N/A N/A 0
1 72 1,360 N/A N/A 1360
2 64 1,150 975 N/A 2125
3 49 850 745 600 2195
When sufficient resources are available, additional residual pressure measurements
can be taken during the fire hydrant flow test at various locations throughout the sys-
tem. Taking these additional pressure readings will provide more information on how
the hydraulic grade changes across the system. Depending on the nature of the water
distribution system, the pressure drop may be localized to the vicinity of the flowed
hydrants.
If the hydrant flow test is conducted to provide data for model calibration, it is
extremely important to note the boundary conditions at the time of the test. Recall that
boundary conditions reflect the water levels in tanks and reservoirs, as well as the
operational status of any high-service pumps, booster pumps, or control valves (for
example, pressure reducing valves) for both static and flowed conditions. As will be
discussed in Chapter 7 (see page 261), these boundary conditions must also be
defined in the hydraulic model.
In addition, system demands in place at the time of the test need to be replicated in the
model. It is important to note the time of day and the weather conditions when the test
was performed to assist in establishing the demands and boundary conditions.
Potential Problems with Fire Flow Tests
Fire hydrant flow tests are a useful tool. They do, however, present some areas of con-
cern. Because the discharges from fire hydrants can be quite large, the following sug-
gestions can reduce potential problems associated with these flows.
1. Minimize the period of time over which hydrants are flowed to limit flooding
potential. (In some locations it may be necessary to dechlorinate water before it
can be discharged into receiving waters.)
2. Direct the flow through the 2 -in. (64 mm) nozzle opening instead of the 4 -in.
(115 mm) opening. This will help to reduce street flooding while still producing
flow velocities sufficient for calibration.
Section 5.2 Fire Hydrant Flow Tests 189
3. Use hydrant diffusors to reduce the high velocity of the hydrant stream. This will
help to avoid erosion problems and damage to vegetation.
4. Conduct fire hydrant flow tests during warm weather to avoid ice problems.
5. Notify customers who may be impacted by the test beforehand. In some systems,
hydrant flow tests can stir up sediments and rust, causing temporary water quality
problems.
6. Make sure to open and close the hydrants gradually, as sudden changes in flow
can induce dangerous pressure surges in the system.
7. Make sure that the residual and flowed hydrants are hydraulically close to one
another. It is possible to have two hydrants that are near each other at the street
but are fed by different mains that may not be hydraulically connected for several
blocks. Ideally, the flowed and residual hydrants would be located side-by-side
on the same pipeline, but because this will almost never be the case, accuracy can
instead be improved by minimizing the flow between the hydrants. (The flow can
often be reduced by bracketing the residual hydrant between two flowed
hydrants.)
Evaluating Distribution Capacity
with Hydrant Tests
The results of hydrant flow tests described in this
chapter are used primarily to evaluate the distribu-
tion systems capacity to provide water for fighting
fires. The standard formula for converting the test
flow to the distribution capacity at some desired
residual pressureusually 20 psi (135 kPa)was
developed by the Insurance Services Office
(1963), and is given in AWWA M-17 (1989) as:
where Q
r
= fire flow at residual pressure P
r
(gpm, l/s)
Q
t
= hydrant discharge during test (gpm, l/s)
P
s
= static pressure (psi, kPa)
P
r
= desired residual pressure (psi, kPa)
P
t
= residual pressure during test (psi, kPa)
The value of Q
r
is referred to as the distribution
main capacity in that location, and is used in eval-
uation of water systems for insurance purposes.
Assumptions made when using the above equa-
tion are as follows:
1. Head loss is negligible during static conditions.
2. Demands correspond to maximum day
demands.
3. All pumps and regulating valves that would
open during an actual fire are open and operat-
ing during the test.
4. There is sufficient water quantity to supply the
fire throughout the duration of the fire event.
5. Tank level is at normal day low level.
6. The residual and flowed hydrants are close to
one another (Walski, 1984b).
Water system models can explicitly account for
these factors and are a more accurate and flexible
way of assessing available fire flow at a given
residual pressure. However, this equation is still
widely used.
The previous equation can also be rearranged to
provide a rough estimate of residual pressure for
some future flow, given hydrant flow test results,
according to
In this case, Q
r
is the estimated flow, and P
r
is the
pressure that will exist at that flow rate, given that
all other conditions remain the same.
Q
r
Q
t
P
s
P
r