Mike Urban: Water Distribution
Mike Urban: Water Distribution
Water Distribution
MIKE 2019
2 MIKE URBAN - © DHI
PLEASE NOTE
LIMITED LIABILITY The liability of DHI is limited as specified in your DHI Software
Licence Agreement:
3
4 MIKE URBAN - © DHI
1 EPANET Input Description . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
1.1 Network Component Editors . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
1.1.1 Junction Editor . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
Emitter . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
Air-Valve . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
Graphical Placement and Editing of Junction Nodes . . . . . . . . . . 17
1.1.2 Pipe Editor . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
Warning . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
Flow Direction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
Graphical Placement and Editing of Pipes . . . . . . . . . . . . . . . 23
Snap Tolerance . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
Defining Curved Pipes . . . . . . . . . . . . . . . . . . . . . . . . . 24
1.1.3 Pump Editor . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
Pump Energy Settings . . . . . . . . . . . . . . . . . . . . . . . . . 29
Pump Water Hammer Settings (optional) . . . . . . . . . . . . . . . . 30
Pump Variable Speed Settings (optional) . . . . . . . . . . . . . . . 31
Variable Pumps . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
Flow Direction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
Graphical Placement and Editing of Pumps . . . . . . . . . . . . . . 33
1.1.4 Valve Editor . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
Valve Water Hammer Settings (optional) . . . . . . . . . . . . . . . . 37
Variable Valves . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38
Check Valves . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38
Flow Direction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38
Graphical Placement and Editing of Valves . . . . . . . . . . . . . . 38
1.1.5 Tank Editor . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39
Graphical Placement and Editing of Tanks . . . . . . . . . . . . . . . 44
1.1.6 Air-Chamber Editor . . . . . . . . . . . . . . . . . . . . . . . . . . . 44
Graphical Placement and Editing of Air-Chambers . . . . . . . . . . . 47
1.2 Other Editors . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 48
1.2.1 Project Options Editor . . . . . . . . . . . . . . . . . . . . . . . . . 48
1.2.2 Pressure Zone Editor . . . . . . . . . . . . . . . . . . . . . . . . . . 52
1.2.3 Energy Editor . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
1.2.4 Report Editor . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
1.2.5 Curves and Relations Editor . . . . . . . . . . . . . . . . . . . . . . 56
1.2.6 Boundary Conditions Editor . . . . . . . . . . . . . . . . . . . . . . . 58
1.3 Network Demand . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 59
1.3.1 Distributed Demands . . . . . . . . . . . . . . . . . . . . . . . . . . 59
1.3.2 Demand Allocation Editor . . . . . . . . . . . . . . . . . . . . . . . . 64
Geocoding . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 66
Aggregation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 68
1.3.3 Multiple Demand Editor . . . . . . . . . . . . . . . . . . . . . . . . . 69
1.3.4 Demand Statistics . . . . . . . . . . . . . . . . . . . . . . . . . . . 71
1.4 Extended Period Simulations . . . . . . . . . . . . . . . . . . . . . . . . . . . 73
1.4.1 Simple Control Editor . . . . . . . . . . . . . . . . . . . . . . . . . . 74
Multiple Controls . . . . . . . . . . . . . . . . . . . . . . . . . . . . 76
Pipes with Check Valves . . . . . . . . . . . . . . . . . . . . . . . . 77
Available Operational Controls . . . . . . . . . . . . . . . . . . . . . 77
5
Time Based Controls . . . . . . . . . . . . . . . . . . . . . . . . . . .77
Storage Tank Water Level Based Controls . . . . . . . . . . . . . . .77
Junction Node Pressure Based Controls . . . . . . . . . . . . . . . . .77
1.4.2 Rule Based Controls . . . . . . . . . . . . . . . . . . . . . . . . . . 78
Multiple Pumps, Valves (WD-Tools) . . . . . . . . . . . . . . . . . . .81
1.4.3 Repetitive Profiles (Pattern) Editor . . . . . . . . . . . . . . . . . . . 83
Diurnal Patterns . . . . . . . . . . . . . . . . . . . . . . . . . . . . .83
Cyclic Profiles . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .86
1.4.4 Time Settings Editor . . . . . . . . . . . . . . . . . . . . . . . . . . 87
Time Units . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .88
Diurnal Curves and Patterns . . . . . . . . . . . . . . . . . . . . . . .88
1.5 Water Quality Simulations . . . . . . . . . . . . . . . . . . . . . . . . . . . . 89
1.5.1 Water Quality Analysis Selection . . . . . . . . . . . . . . . . . . . . 89
1.5.2 Water Quality Analysis Parameters . . . . . . . . . . . . . . . . . . 90
1.5.3 Initial Water Quality Editor . . . . . . . . . . . . . . . . . . . . . . . 91
1.5.4 Point Constituent Source Editor . . . . . . . . . . . . . . . . . . . . 93
Point Constituent Source Time Patterns . . . . . . . . . . . . . . . . .94
1.5.5 Reaction Rate Editor . . . . . . . . . . . . . . . . . . . . . . . . . . 95
Global Settings . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .95
1.5.6 Local Settings . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 96
Tank Reaction Rate Assignments . . . . . . . . . . . . . . . . . . . .96
Pipe Reaction Rate Assignments . . . . . . . . . . . . . . . . . . . .96
Note . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .97
1.5.7 Source Tracing . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 97
1.6 Tools . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 98
1.6.1 Send Data . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 98
1.6.2 Engineering Tables . . . . . . . . . . . . . . . . . . . . . . . . . . . 98
1.6.3 Convert Asset Valves/Pumps to Model Valves/Pumps . . . . . . . 100
1.6.4 Duplicate Pipe Attributes . . . . . . . . . . . . . . . . . . . . . . . 101
7
3.2 Pipe Roughness Calibration . . . . . . . . . . . . . . . . . . . . . . . . . . 139
3.2.1 Evolutionary Algorithms - Theoretical Background . . . . . . . . . . 139
3.2.2 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 141
3.2.3 Pipe Roughness Calibration . . . . . . . . . . . . . . . . . . . . . 141
3.3 PID Control . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 145
3.3.1 Architecture of RTC Systems . . . . . . . . . . . . . . . . . . . . 145
3.3.2 MIKE URBAN vs. Real Life . . . . . . . . . . . . . . . . . . . . . 147
3.3.3 The PID Algorithm . . . . . . . . . . . . . . . . . . . . . . . . . . 149
Control Settings . . . . . . . . . . . . . . . . . . . . . . . . . . . . 151
Set Point Settings . . . . . . . . . . . . . . . . . . . . . . . . . . . 152
PID Settings . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 152
3.4 Demand Adjusted Anlaysis . . . . . . . . . . . . . . . . . . . . . . . . . . 152
3.4.1 Low Pressure Systems . . . . . . . . . . . . . . . . . . . . . . . . 155
3.4.2 Intermittent Water Supply Systems or Modelling System Shutdown and
Maintenance 157
3.5 Variable Speed Drive Pump . . . . . . . . . . . . . . . . . . . . . . . . . . 157
3.6 Extended Rule-Based Controls . . . . . . . . . . . . . . . . . . . . . . . . 157
3.7 Pressure Zone Mapping . . . . . . . . . . . . . . . . . . . . . . . . . . . . 157
3.8 Advisor . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 159
3.8.1 Pipe Criticality . . . . . . . . . . . . . . . . . . . . . . . . . . . . 159
3.8.2 Sustainability Analysis . . . . . . . . . . . . . . . . . . . . . . . . 162
3.8.3 Cost Analysis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 165
Pump Energy Table . . . . . . . . . . . . . . . . . . . . . . . . . . 165
Pump Energy Graphs . . . . . . . . . . . . . . . . . . . . . . . . . 166
3.8.4 Alarms and Violations . . . . . . . . . . . . . . . . . . . . . . . . 167
4 MIKE URBAN Water Hammer . . . . . . . . . . . . . . . . . . . . . . . . . . 171
4.1 Water Hammer Calculation . . . . . . . . . . . . . . . . . . . . . . . . . . 171
4.2 Theoretical Background . . . . . . . . . . . . . . . . . . . . . . . . . . . . 172
4.2.1 Description of Water Hammer Model . . . . . . . . . . . . . . . . 172
Wave Speed . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 172
Rigid Conduit . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 173
Thick-Walled Elastic Conduit (D/e<=10) . . . . . . . . . . . . . . . . 173
Thin-Walled Elastic Conduit (D/e>10) . . . . . . . . . . . . . . . . . 173
Tunnels Through Solid Rock, Parmakian 1963 . . . . . . . . . . . . 174
Reinforced Concrete Pipe . . . . . . . . . . . . . . . . . . . . . . . 174
Diagrams . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 175
4.3 Numerical Scheme and Algorithm . . . . . . . . . . . . . . . . . . . . . . . 176
4.3.1 Coefficients for the Numerical Scheme . . . . . . . . . . . . . . . 177
4.3.2 Looped Network Solution Algorithm . . . . . . . . . . . . . . . . . 178
4.3.3 Hydraulic Structures . . . . . . . . . . . . . . . . . . . . . . . . . 180
4.4 Water Hammer Calculations . . . . . . . . . . . . . . . . . . . . . . . . . . 181
Water Hammer Data Preparation . . . . . . . . . . . . . . . . . . . 181
Water Hammer Result Presentation . . . . . . . . . . . . . . . . . . 181
4.4.1 Running Water Hammer Simulations . . . . . . . . . . . . . . . . 182
Initial Conditions . . . . . . . . . . . . . . . . . . . . . . . . . . . . 182
Boundary Conditions . . . . . . . . . . . . . . . . . . . . . . . . . . 182
Computational Parameters . . . . . . . . . . . . . . . . . . . . . . . 182
9
5.3 Time Patterns . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 225
5.4 Hydraulic Simulation Model . . . . . . . . . . . . . . . . . . . . . . . . . . 226
5.4.1 Water Age and Source Tracing . . . . . . . . . . . . . . . . . . . 233
5.5 Water Quality Simulation Model . . . . . . . . . . . . . . . . . . . . . . . . 234
5.5.1 Advective Transport in Pipes . . . . . . . . . . . . . . . . . . . . . 234
5.5.2 Mixing at Pipe Junctions . . . . . . . . . . . . . . . . . . . . . . . 234
5.5.3 Mixing in Storage Facilities . . . . . . . . . . . . . . . . . . . . . . 235
5.5.4 Bulk Flow Reactions . . . . . . . . . . . . . . . . . . . . . . . . . 235
5.5.5 Pipe Wall Reactions . . . . . . . . . . . . . . . . . . . . . . . . . 237
5.5.6 Lagrangian Transport Algorithm . . . . . . . . . . . . . . . . . . . 238
5.5.7 Tank Mixing Models . . . . . . . . . . . . . . . . . . . . . . . . . 240
5.6 Model Calibration . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 243
5.6.1 Accuracy Concerns . . . . . . . . . . . . . . . . . . . . . . . . . 243
Pressure Differential Method . . . . . . . . . . . . . . . . . . . . . . 244
Percent Pressure Differential Method . . . . . . . . . . . . . . . . . 244
Head Loss Differential Method . . . . . . . . . . . . . . . . . . . . . 244
5.6.2 Reasons for Calibrating a Model . . . . . . . . . . . . . . . . . . . 244
Model Credibility . . . . . . . . . . . . . . . . . . . . . . . . . . . . 245
Benchmark . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 245
Predict Potential Problems . . . . . . . . . . . . . . . . . . . . . . . 245
Understand System Operations . . . . . . . . . . . . . . . . . . . . 245
Uncover Errors . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 245
5.6.3 Calibration Model Data Requirements . . . . . . . . . . . . . . . . 246
Simulation Type Considerations . . . . . . . . . . . . . . . . . . . . 246
Data Acquisition . . . . . . . . . . . . . . . . . . . . . . . . . . . . 247
Operational Review . . . . . . . . . . . . . . . . . . . . . . . . . . 247
Recalibration Frequency . . . . . . . . . . . . . . . . . . . . . . . . 247
5.6.4 Calibration Simulations . . . . . . . . . . . . . . . . . . . . . . . . 248
Initial Simulation . . . . . . . . . . . . . . . . . . . . . . . . . . . . 248
Comparing Model Output . . . . . . . . . . . . . . . . . . . . . . . 248
5.6.5 Model Adjustments . . . . . . . . . . . . . . . . . . . . . . . . . . 250
Operational Data Adjustments . . . . . . . . . . . . . . . . . . . . . 251
Consumption Data Adjustments . . . . . . . . . . . . . . . . . . . . 251
Network Data Adjustments . . . . . . . . . . . . . . . . . . . . . . . 252
Pipe Roughness Adjustments . . . . . . . . . . . . . . . . . . . . . 252
Number of Simulation Iterations . . . . . . . . . . . . . . . . . . . . 252
The Model Manager User Guide will provide information on how to work with
the various editing tools, loading of results etc.
The network component editors are all available from the EPANET Menu, as
shown in Figure 1.1.
Figure 1.1 The EPANET Menu allow access to the network component editors
11
EPANET Input Description
The first step in defining a water distribution network is to define the junction
nodes that define the interconnection between the pipes, pumps, valves,
tanks, and reservoirs that make up the network. Junction nodes are also
placed at points of water consumption or inflow, at points where specific anal-
ysis values (e.g., pressure, concentration, etc.) are desired, and at points
where pipe attributes (e.g., diameter, roughness, etc.) change.
Junction nodes are either defined interactively on the graphical Map window
using the Junction tool, or by manual data entry using the Junction Editor dia-
log box as shown in Figure 1.2. The Junction Editor allows you to define the
junction’s ID, location, any external demand, and a description. The Junction
Editor dialog box is reached by selecting EPANET | Junctions
Figure 1.2 The Junction Editor allows you to define the junction nodes that define
the inter connectivity between the water distribution network compo-
nents
A list of the Junction Editor data entries for Figure 1.2 follows, with a short
description given for each entry.
ASSET ID (optional)
This data entry is used to specify a corresponding asset junction ID, which
uniquely identifies the junction node in the asset management system (such
as GIS, for example).
JUNCTION ID (mandatory)
This data entry is used to specify an ID which uniquely identifies the junction
node. The junction ID acts as a unique look up key that identifies the node
from all other nodes. A node can be a junction, reservoir, or tank. Therefore,
no two nodes may have the same ID. However, a node and a link (i.e., pipe,
pump, or valve) can have the same ID. The node ID value can be any string
value (up to 40 characters).
DESCRIPTION (optional)
This data entry allows you to enter a description identifying the junction node
being entered. This description can be optionally displayed on the Map win-
dow and in reports generated by the Report Generator.
STATUS (optional)
This drop down selection list data entry allows you to define whether the junc-
tion node is imported (i.e existing node was imported from the external data
source), or is inserted, modified, GIS, calibrated or similar. By default, junc-
tion node status is undefined.
Selecting «...» allows you to display the Pressure Zone selection dialog box,
where the appropriate pressure zone ID can be selected.
13
EPANET Input Description
TYPE (mandatory)
This drop down selection list data entry allows you to define whether the junc-
tion node is a simple junction, or is an emitter, or an air-valve. By default, any
junction nodes added are junctions.
DEMAND COEFFICIENT
The new field was added into the Junction editor. Demand coefficient allows
you to define the share from the whole network demand, which is taken by
the node. This filed is used only by the Demand Distribution function.
Example: the network has 3 nodes, where the demand coefficient is defined;
these values are 10, 10, and 30. For each node, the weighted coefficient is
calculated and based on it; the total network demand is distributed. The node,
where the demand coefficient is not defined will get no demand from the total
network demand: Q
QT
q i = ----------------- Ci (1.1)
i C
i = 0,n
where:
qi = node demand
Qt = total network demand
ci = demand coefficient
STATE (optional)
This drop down selection list data entry allows you to define whether the junc-
tion node is unmarked (i.e an existing node already contained in the water
distribution network), or is marked (i.e. one that is being considered for future
expansion, etc.). By default, any junction nodes added are unmarked.
ELEVATION (mandatory)
This data entry defines the elevation above a common datum for the junction
node, in units of ft. or m. This value is used to determine the pressure head at
the node during a simulation. The default elevation is zero.
Junction nodes should have their elevation specified above zero datum (i.e.,
sea level) so that pressure computations can be carried out.
Note that the baseline demand at this junction node can be (optionally) com-
puted by globally defining the demand for the entire network and then having
MIKE URBAN distribute this demand to each of the network nodes using the
Distributed Demand dialog box. See the section titled “Distributed Demands”
on page 59 for more information on computing distributed demands.
For extended period simulations, the baseline demand can be adjusted dur-
ing the simulation due to a predefined time pattern factor. See the section
titled “Repetitive Profiles (Pattern) Editor” on page 83 for more information.
15
EPANET Input Description
Emitter
The following data entries defines parameters of an emitter located at the
junction node. Emitters are needed to model flow through sprinkler systems
and irrigation networks. They can also be used to simulate leakage in a pipe
connected to the junction (if a discharge coefficient for the leading crack or
joint can be estimated).
A list of the Emitter data entries for Figure 1.3 follows, with a short description
given for each entry.
Air-Valve
The following data entries defines parameters of an air-valve located at the
junction node and are used only in case of a water hammer analysis. A list of
the Air-Valve data entries for Figure 1.3 follows, with a short description given
for each entry. See the section titled “Water Hammer Calculation” on
page 171 for more information on demand patterns.
By default this data entry is k=1.2. See the section titled “Specific Valve Data”
on page 191.
If desired, from the Map window, the user can move an existing junction node
using the ‘select’ tool. By selecting the node and holding down the left mouse
button, the node can be dragged to a new location. As the node is dragged,
its connecting links (i.e., pipes, pumps, and valves) rubber band along with it.
To delete an existing junction node graphically, select the node using the
‘Select’ tool and then press «Delete». The selected node and all its connect-
17
EPANET Input Description
ing links will then be deleted. To edit a junction node, double click on the node
using the Select tool. The Junction Editor will then be displayed allowing you
to edit the attributes of the selected node.
Pipes are used to transport water from one node to another. Pipes must
always begin and end at a node.
Pipes are either defined interactively on the Map window using the ‘Add pipe’
tool, or by manual data entry using the Pipe Editor dialog box as shown in
Figure 1.5
.
Figure 1.5 The Pipe Editor allows you to define the pipes that make up the water
distribution network
The Pipe Editor allows you to define the pipe’s ID, diameter, loss coefficients,
nodal connectivity, description, and other attributes. The Pipe Editor dialog
box is reached by selecting EPANET | Pipes.
A list of the Pipe Editor data entries for Figure 1.5 follows, with a short
description given for each entry.
ASSET ID (optional)
This data entry is used to specify a corresponding asset pipe ID, which
uniquely identifies the pipe link in the asset management system (such as
GIS, for example).
PIPE ID (mandatory)
This data entry is used to specify an ID which uniquely identifies the pipe link.
The pipe ID acts as a unique look up key that identifies the link from all other
links. A link can be a pipe, pump, or valve. Therefore, no two links may have
the same link ID. However, a link and node (i.e., junction, reservoir, or tank)
can have the same ID. The node ID value can be any string value (up to 40
characters).
When importing (or merging) multiple water distribution network models into a
single network model, MIKE URBAN will check for collisions between identi-
cal link IDs and will automatically assign a new link ID value for any link being
imported that contains the same link ID value as what already exists in the
network model. MIKE URBAN will automatically assign a new pipe ID when-
ever any pipe is split into more pipes.
The Node Type pull-down selection list allows the user to specify what type of
node is connected to the end of the pipe. Then, choosing «...» will display the
Select Node selection dialog box from which the user can select the appropri-
ate node. Or, choosing «Pick» allows the user to graphically select the node
from the Map window. To reverse the order of the nodes, choose «Advanced
- Swap Nodes».
If the computed flow is moving from the starting node to the ending node, the
computed flow value will be positive. If the computed flow is moving from the
ending node to the starting node, the computed flow value will be negative.
DESCRIPTION (optional)
This data entry allows you to enter a description identifying the pipe being
entered. This description can be optionally displayed on the Map window and
in reports generated by the Report Generator.
19
EPANET Input Description
STATUS (optional)
This drop down selection list data entry allows you to define whether the
pump is imported (i.e existing node was imported from the external data
source), or is inserted, modified, GIS, calibrated or similar. By default, pump
status is undefined.
Selecting «...» allows you to display the Pressure Zone selection dialog box,
where the appropriate pressure zone ID can be selected.
LENGTH (mandatory)
This data entry defines the pipe length, in the unit of your choice. It is possible
to define a specific pipe length—independent of the pipe network layout. Oth-
erwise, MIKE URBAN recomputes the pipe length based upon the pipe lay-
out.
DIAMETER (mandatory)
This data entry defines the internal diameter of the pipe, in inches or mm. The
second field (read-only) displays the pipe diameter as it would be used for the
hydraulic analysis. The pipe diameter is automatically adjusted when the pipe
wall is defined.
MATERIAL (optional)
This option allows the user to define the material of pipe construction. The
Pipe Material is defined by a “string” and the pipe roughness can be globally
assigned based upon the pipe material and pipe construction year, for exam-
ple.
FORMULATION (read-only)
This field is used to display the current friction loss formulation type, such as
Darcy-Weisbach, Hazen-Williams, or Manning. To change the friction formu-
lation, use the Project Options dialog box.
ROUGHNESS (mandatory)
This data entry defines the roughness of the interior surface of the pipe.
Based upon which roughness type loss coefficient has been specified, this
value is unit less for Hazen-Williams or Chezy-Manning headloss formulas,
and in millifeet or mm for the Darcy-Weisbach (or Colebrook-White) formula-
tion. The roughness type is specified by the user within the Project Options
dialog box. Choosing «...» will display the Select Pipe Roughness Coefficient
selection dialog box, allowing the user to select the appropriate roughness
value to use.
These two data entries define the demand coefficients, k1i and k2i which are
used in these two computational methods. Each method proportions the total
network demand along the pipe by splitting the computed demand between
the two nodes at each end of the pipe.
21
EPANET Input Description
The Method of Pipe Lengths computes the total water demand assigned to
the current pipe (which is then split between the starting and ending nodes)
as:
Q l i k li
q pi = -----------------
- (1.2)
k li l i
The Method of Two Coefficients computes the total water demand assigned
to the current pipe (which is then split between the starting and ending nodes)
as:
Q k li k 2i
q pi = ---------------------
- (1.3)
k li k 2i
where:
qpi=total water demand applied to the pipe, split between the two end nodes.
Q= total network water demand. Note that the total network water demand is
defined in the Distributed Demands dialog box, available from the Edit Menu.
li = pipe length
See the section titled “Distributed Demands” on page 59 for more information
on computing distributed demands.
Note that there is also Demand Coefficient 3, 4 which is available only in the
table grid.
STATE (optional)
This drop down list box data entry allows you to define whether the pipe is
unmarked (i.e. a pipe already contained in the water distribution network), or
is marked (i.e. one that is being considered for future expansion, etc.). By
default, any pipes added are unmarked.
CLOSED (optional)
This check box data entry allows the user to toggle the OPEN and CLOSED
status of the pipe. Choosing CLOSED effectively removes the pipe from the
network system.
Warning
Note that you cannot set the pipe status of a pipe containing a check valve.
Pipes with a check valve are initially open, and close only if flow within the
pipe attempts to reverse (move from the ending downstream node to the
starting upstream node).
Flow Direction
Note that the flow direction for pipes is assumed to be from the starting
(upstream) node to the ending (downstream) node. The order in which these
nodes are specified for each pipe is critical for a proper network representa-
tion. To reverse the order of the nodes, choose «Swap·Nodes».
If the computed flow is moving from the starting node to the ending node, the
computed flow value will be positive. If the computed flow is moving from the
ending node to the starting node, the computed flow value will be negative.
While selecting the starting node, if the user clicks with the ’Add pipe’ tool in
the Map window somewhere else other than a node or a pipe, MIKE URBAN
will place a starting node at the clicked location.
If the user selects an existing pipe while selecting either the starting or ending
nodes, MIKE URBAN will display a query dialog, as shown in Figure 1.6, ask-
ing the user whether to split the pipe and insert a node at the selected loca-
tion. If the user selects to insert a node into the pipe, the pipe will be broken
into two segments at the point of selection with a new node added. The
length of the two pipe segments will be adjusted accordingly so that the total
length of the two segments equals that of the original pipe.
23
EPANET Input Description
Figure 1.6 When clicking on an existing pipe with the Add Pipe tool, MIKE URBAN
will ask whether to insert a junction node into the existing pipe for con-
necting the new pipe
To delete an existing pipe graphically, select the pipe using the ‘Select’ tool
and then select «Delete Selected Features». The selected pipe will then be
deleted. To edit a pipe, double click on the pipe using the ‘Select’ tool. The
Pipe Editor will then be displayed allowing you to edit the attributes of the
selected pipe.
Snap Tolerance
As you graphically layout pipes in the Map window by clicking on existing
nodes, MIKE URBAN will snap to the nearest node if it is within a specified
snap tolerance. However, if the cursor is not placed directly on top of the
node or within the specified snap tolerance, MIKE URBAN may not snap to
the node but instead will add a new junction node.
1. Select the ‘Add Pipe’ tool from the Components toolbar and click on the
starting point in the Map window.
2. After selecting the starting point, construct the curve of the pipe by click-
ing in the Map window at the points the curvature changes to create
intermediate vertices so that an arc or a curved line is formed. Double-
click on the last position to end the curved pipe. If the last position is an
existing node, the end point of the curved pipe will snap to the existing
node.
To calculate the pipe length based upon the curved layout, press the
Advanced button on the Pipe Editor and select either ‘Recompute selected
pipe length’ or ‘Recompute all pipes length’. MIKE URBAN will then recom-
pute the pipe length for either the selected pipes or all the pipe based upon
the pipe layout.
Pumps are used to raise the hydraulic head of water. Pumps are represented
as short links of negligible length. MIKE URBAN analysis engine will automat-
ically prevent flow reversal though a pump, and will issue warning messages
when a pump operates outside of its normal operating range.
Pumps are either defined interactively on the graphical Map window using the
Add Pump tool, or by manual data entry using the Pump Editor dialog box as
shown in Figure 1.7.
The Pump Editor allows you to define the pump’s ID, pump power curve, sta-
tus, nodal connectivity, description, and other attributes. The Pump Editor
dialog box is reached by selecting EPANET | Pumps.
Figure 1.7 The Pump Editor allows you to define the pumps used in the water dis-
tribution network
25
EPANET Input Description
A list of the Pump Editor data entries for Figure 1.7 follows, with a short
description given for each entry.
ASSET ID (optional)
This data entry is used to specify a corresponding asset pump ID, which
uniquely identifies the junction node in the asset management system (such
as GIS, for example).
PUMP ID (mandatory)
This data entry is used to specify an ID which uniquely identifies the pump
link. The pump ID acts as a unique lookup key that identifies the link from all
other links. A link can be a pipe, pump, or valve. Therefore, no two links may
have the same link ID. However, a link and node (i.e., junction, reservoir, or
tank) can have the same ID. The node ID value can be any string value (up to
40 characters).
The Node Type pull-down selection list allows the user to specify what type of
node is connected to the end of the pump. Then, choosing «...» will display
the Select Node selection dialog box from which the user can select the
appropriate node. Or, choosing «Pick» allows the user to graphically select
the node from the Map window.
Pumped flow is always assumed to move from the starting (upstream) node
to the ending (downstream) node. To reverse the order of the nodes, choose
«Swap·Nodes».
DESCRIPTION (optional)
This data entry allows you to enter a description identifying the pump being
entered. This description can be optionally displayed on the Map window and
in reports generated by the Report Generator.
STATUS (optional)
This drop down selection list data entry allows you to define whether the
pump is imported (i.e existing node was imported from the external data
source), or is inserted, modified, GIS, calibrated or similar. By default, pump
status is undefined.
Selecting «...» allows you to display the Pressure Zone selection dialog box,
where the appropriate pressure zone ID can be selected.
1xQ-H Point Pump Curve pump, as shown in Figure 1.8, is used for a stand-
ard pump curve with no extended flow range, where the cutoff head is 133%
of the design head and the maximum flow is twice the design flow.
The Design Head is in units ft. or m, and are by default zero. The Design Flow
is in the user-specified units, and by default zero.
3xQ-H Point Pump Curve, as shown in Figure 1.10 can be used to describe
the flow-head relationship of the pump.
The Shutoff Head is the head value at zero flow. The Design Head is the
standard operating head. The Design Flow is corresponding flow rate. The
27
EPANET Input Description
High End Head is the head at the upper end of the normal operating flow
range. The High End Flow is the corresponding flow rate. The Maximum Flow
is the flow rate for the extended flow range.
All heads are in units of ft. or m, and flows are in the user-specified units.
Q-H Pump Curve can alternatively. The Q-H Pump Curve is created by pro-
viding either a pair of head-flow points, or four or more such points. MIKE
URBAN creates the pump curve by connecting the points with straight line
segments. The Q-H pump curve should be created in the Curve Editor found
in the Edit menu. Once the Q-H pump curve is created, it can then be
assigned to a pump in the Pump Editor.
Figure 1.9 A Q-H Table used to define the multi-point pump curve
CLOSED (optional)
This check box data entry allows the user to toggle the OPEN and CLOSED
status of the pump. Choosing CLOSED effectively removes the pump from
the network system.
29
EPANET Input Description
31
EPANET Input Description
Please note that it is possible to use this control also for the downstream
node by specifying the downstream node ID as the "Control Node".
Variable Pumps
Many times, when performing an extended period simulation, it is desirable to
model a variable pump. A variable pump can vary its speed setting and/or
change its status to open or closed during a simulation. A variable pump is
modelled by defining the pump’s initial settings using the Pump Status previ-
ously described and then modifying the pump’s operation during the
extended period simulation using the Control Editor. The Control Editor is
available from the Extended Menu. See the section titled “Simple Control Edi-
tor” on page 74 for additional information.
A special case if a variable pump is VSD (variable speed drive) pump where
the rotational speed of a pump is automatically adjusted by the program in
order to maintain the control node pressure at the required level. MIKE
URBAN WATER can model VSD pump using VSD settings in the Pump Edi-
tor or using the PID (Proportional Integral Differential Control) algorithm in the
RTC (Real-Time Control) module.
Flow Direction
Note that the flow direction for pumps is assumed to be from the starting
(upstream) node to the ending (downstream) node. The order in which these
nodes are specified is critical for a proper network representation of the
pump. To reverse the order of the nodes, choose «Swap·Nodes».
1. Split the pipe by using Split tool from the Components toolbar.
2. Split the pipe once more to create a sections for 3 pipes.
3. Select the middle section by Select Feature tool and delete it by using
the Delete Selected tool.
4. Select the Pump tool from the Components floating toolbar and click on
the beginning pipe node in the Horizontal Plan window where you want
to place the valve and double click on the ending node. This will create a
new pump.
Another method of using the ‘Add Pump’ tool to add a pump is, from the Map
window, select a starting node and then an ending node. Clicking on or near
an existing node, MIKE URBAN will then snap to the selected node and a
rubber-banding line will then be drawn from this node, representing the start-
ing node, while the user selects the ending node. After selecting the ending
node, MIKE URBAN will then place a pump between the two selected nodes.
In addition, if the user clicks with the ‘Add Pump’ tool in the Map window
somewhere else other than a node or a pipe, MIKE URBAN will place a start-
ing node at the clicked location. A rubber-banding line will then be drawn from
this node, representing the pump, while the user selects the ending node. If
the user clicks on or near an existing node, it will treat the selected node as
the ending node. If the user clicks anywhere else, it will place an ending node
33
EPANET Input Description
at the selected location. If the user selects an existing pipe, an error message
will be displayed stating that an invalid ending node was selected.
To delete an existing pump graphically, select the pump using the Select tool
and then select «Delete Selected». The selected pump will then be deleted.
To edit a pump, click on the pump using the Open Editor tool. The Pump Edi-
tor will then be displayed allowing you to edit the attributes of the selected
pump.
Valves control the flow or pressure of water from one junction node to
another. Valves are represented as short links of negligible length. Note that
valve pressure settings are pressures (e.g., psi or m) and not total head (or
hydraulic gradeline elevation). Also, valves cannot be directly connected to
reservoir or storage tank nodes.
Valves are either defined interactively on the Map window using the
Add Valve tool, or by manual data entry using the Valve Editor dialog box as
shown in Figure 1.13. The Valve Editor allows you to define the valve’s ID,
type, status, nodal connectivity, description, and other attributes. The Valve
Editor dialog box is reached by selecting EPANET|Valves.
Figure 1.13 The Valve Editor allows you to define the valves used in the water distri-
bution network
A list of the Valve Editor data entries for Figure 1.13 follows, with a short
description given for each entry.
ASSET ID (optional)
This data entry is used to specify a corresponding asset valve ID, which
uniquely identifies the valve in the asset management system (such as GIS,
for example).
VALVE ID (mandatory)
This data entry is used to specify an ID which uniquely identifies the valve
link. The valve ID acts as a unique lookup key that identifies the link from all
other links. A link can be a pipe, pump, or valve. Therefore, no two links may
have the same link ID. However, a link and node (i.e., junction, reservoir, or
tank) can have the same ID. The valve ID value can be any string value (up to
40 characters).
The Node Type pull-down selection list allows the user to specify what type of
node is connected to the end of the pump. Then, choosing «...» will display
the Select Node selection dialog box from which the user can select the
appropriate node. Or, choosing «Pick» allows the user to graphically select
the node from the Map window.
Controlled flow is always assumed from the starting (upstream) node to the
ending (downstream) node. To reverse the order of the nodes, choose
«Swap·Nodes».
DESCRIPTION (optional)
This data entry allows you to enter a description identifying the valve being
entered. This description can be optionally displayed on the Map window and
in reports generated by the Report Generator.
STATUS (optional)
This drop down selection list data entry allows you to define whether the
valve is imported (i.e existing node was imported from the external data
source), or is inserted, modified, GIS, calibrated or similar. By default, valve
status is undefined.
35
EPANET Input Description
Selecting «...» allows you to display the Pressure Zone selection dialog box,
where the appropriate pressure zone ID can be selected.
A Flow Control Valve (FCV) limits the flow through a valve to a specified
amount. The program will produce a warning message if this flow cannot be
maintained without having to add additional head at the valve.
A General Purpose Valve (GPV) provides the capability to model devices and
situations with unique headloss - flow relationships, such as reduced pres-
sure backflow prevention valves, turbines, and well drawdown behaviour. The
valve setting is the ID of a Headloss Curve.
SETTING (mandatory)
This data entry defines the pressure setting for PRVs, PSVs, and PBVs,
whose units are in psi or m. Or, this data entry defines the flow settings (in
user-defined flow units) for FCVs, or loss coefficients for TCVs.
When defining a pressure setting, the value specified is pressure (e.g., psi or
m) and not total head (or hydraulic gradeline elevation).
DIAMETER (mandatory)
This data entry is used to specify an valve diameter (inch or mm).
CLOSED (optional)
This option check box data entry allows the user to toggle the OPEN and
CLOSED status of the valve. Choosing CLOSED effectively removes the
valve from the network system.
37
EPANET Input Description
Variable Valves
Many times, when performing an extended period simulation, it is desirable to
model a variable valve. A variable valve can vary its status to open or closed
during a simulation. A variable valve is modelled by defining the valve’s initial
status using the Valve Editor and then modifying the valve’s status during the
extended period simulation using the Control Editor. The Control Editor is
available from the Extended Menu.
See the section titled Control Editor or Rule Based Control Editor on -74 for
additional information.
Check Valves
Note that check valves are defined within the Pipe Editor.
Flow Direction
Note that the flow direction for valves is assumed to be from the starting
(upstream) junction node to the ending (downstream) junction node. The
order in which junction nodes are specified for each valve is critical for a
proper network representation of the valve. To reverse the order of the junc-
tion nodes, choose «Swap·Nodes».
1. Split the pipe by using Split tool from the Components toolbar.
2. Split the pipe once more to create a sections for 3 pipes.
3. Select the middle section by Select Feature tool and delete it by using
the Delete Selected tool.
4. Select the Valve tool from the Components floating toolbar and click on
the beginning pipe node in the Horizontal Plan window where you want
to place the valve and double click on the ending node. This will create a
new Valve.
Another method of using the Add Valve tool to add a Valve is, from the Map
window, select a starting node and then an ending node. Clicking on or near
an existing node, MIKE URBAN will then snap to the selected node and a
rubber-banding line will then be drawn from this node, representing the start-
ing node, while the user selects the ending node. After selecting the ending
node, MIKE URBAN will then place a Valve between the two selected nodes.
In addition, if the user clicks with the Add Valve tool in the Map window some-
where else other than a node or a pipe, MIKE URBAN will place a starting
node at the clicked location. A rubber-banding line will then be drawn from
this node, representing the Valve, while the user selects the ending node. If
the user clicks on or near an existing node, it will treat the selected node as
the ending node. If the user clicks anywhere else, it will place an ending node
at the selected location. If the user selects an existing pipe, an error message
will be displayed stating that an invalid ending node was selected.
To delete an existing Valve graphically, select the Valve using the Select tool
and then select «Delete Selected». The selected Valve will then be deleted.
To edit a Valve, click on the Valve using the ‘Open Editor’ tool. The Valve Edi-
tor will then be displayed allowing you to edit the attributes of the selected
Valve.
Tank nodes are also placed at points in the water distribution model where a
water storage tank is located. Storage tanks can be defined as tanks with the
variable or fixed water level. The tank with the variable water level are mod-
elled as tanks where the water surface level changes with time as water flows
into and out of the tank. The tanks with the fixed water level represent places
(reservoir) within the water distribution model where a infinite source of water
(for the sake of the modelling simulation) is available. Hence, the reservoir
water level remains constant during the course of the simulation.
Tank nodes are either defined interactively on the graphical Map window
using the Add Tank tool, or by manual data entry using the Tank Editor dialog
box as shown in Figure 1.15. The Tank Editor allows you to define the reser-
voir’s ID, location, pressure zone, water surface elevation, and a description.
The Tank Editor dialog box is reached by selecting EPANET | Tanks.
39
EPANET Input Description
Figure 1.15 The Tank Editor allows you to define the storage tank nodes that supply
water to the water distribution network
A list of the Tank Editor data entries for Figure 1.15 follows, with a short
description given for each entry.
ASSET ID (optional)
This data entry is used to specify a corresponding asset tank ID, which
uniquely identifies the tank node in the asset management system (such as
GIS, for example).
TANK ID (mandatory)
This data entry is used to specify an ID which uniquely identifies the tank
node. The tank ID acts as a unique lookup key that identifies the node from all
other nodes. A node can be a junction, reservoir, tank, or air-chamber. There-
fore, no two nodes may have the same ID. However, a node and a link (i.e.,
pipe, pump, or valve) can have the same ID. The node ID value can be any
string value (up to 40 characters).
When importing (or merging) multiple water distribution network models into a
single network model, MIKE URBAN will check for collisions between identi-
cal node IDs and will automatically assign a new node ID value for any node
being imported that contains the same node ID value as what already exists
in the network model.
DESCRIPTION (optional)
This data entry allows you to enter a description identifying the tank node
being entered. This description can be optionally displayed on the Map win-
dow and in reports generated by the Report Generator.
STATUS (optional)
This drop down selection list data entry allows you to define whether the tank
node is imported (i.e existing node was imported from the external data
source), or is inserted, modified, GIS, calibrated or similar. By default, tank
node status is undefined.
Selecting «...» allows you to display the Pressure Zone selection dialog box,
where the appropriate pressure zone ID can be selected.
The Completely mixed model assumes that all water that enters a tank is
instantaneously and completely mixed with the water already in the tank. It is
the simplest form of mixing behaviour to assume, requires no extra parame-
41
EPANET Input Description
ters to describe it, and seems to apply quite well to a large number of facilities
that operate in fill-and-draw fashion.
The First-In-First-Out (FIFO) Plug Flow mixing model assumes that there is
no mixing of water at all during its residence time in a tank. Water parcels
move through the tank in a segregated fashion where the first parcel to enter
is also the first to leave. Physically speaking, this model is most appropriate
for baffled tanks that operate with simultaneous inflow and outflow. There are
no additional parameters needed to describe this mixing model.
The Last-In-First-Out (LIFO) Plug Flow mixing model assumes that there is
no mixing between parcels of water that enter a tank. However in contrast to
FIFO Plug Flow, the water parcels stack up one on top of another, where
water enters and leaves the tank on the bottom. Physically speaking this type
of model might apply to a tall, narrow standpipe with an inlet/outlet pipe at the
bottom and a low momentum inflow. It requires no additional parameters be
provided.
43
EPANET Input Description
If the user clicks with the Add Tank tool in the Map window somewhere else
other than a node or a pipe, MIKE URBAN will place a tank node at the
clicked location. If the user selects an existing pipe, a program message will
be displayed allowing to split the pipe at the cursor location.
If desired, from the Map window, the user can move an existing tank using
the Move Selected tool. By selecting the tank and holding down the left
mouse button, the tank can be dragged to a new location. As the tank is
dragged, its connecting links (i.e., pipes, pumps, and valves) rubber band
along with it.
To delete an existing tank graphically, select the tank using the Select tool
and then select «Delete Selected». The selected tank will then be deleted. To
edit a tank, click on the tank using the Open Editor tool. The Tank Editor will
then be displayed allowing you to edit the attributes of the selected tank.
Air-Chamber nodes are also placed at points in the water distribution model
where a water storage Air-Chamber is located. Storage Air-Chambers can be
defined as Air-Chambers with the variable or fixed water level. The Air-Cham-
ber with the variable water level are modelled as Air-Chambers where the
water surface level changes with time as water flows into and out of the Air-
Chamber. The Air-Chambers with the fixed water level represent places (res-
ervoir) within the water distribution model where a infinite source of water (for
the sake of the modelling simulation) is available. Hence, the reservoir water
level remains constant during the course of the simulation.
Air-Chamber nodes are either defined interactively on the graphical Map win-
dow using the ‘Add Air-Chamber’ tool, or by manual data entry using the Air-
Chamber Editor dialog box as shown in Figure 1.15. The Air-Chamber Editor
allows you to define the reservoir’s ID, location, pressure zone, water surface
elevation, and a description. The Air-Chamber Editor dialog box is reached
by selecting EPANET | Air-Chambers.
Figure 1.17 The Air-Chamber Editor allows you to define the storage tank nodes
that supply water to the water distribution network
A list of the Air-Chamber Editor data entries for Figure 1.15 follows, with a
short description given for each entry.
AIR-TANK ID (mandatory)
This data entry is used to specify an ID which uniquely identifies the air-
chamber node. The tank ID acts as a unique lookup key that identifies the
node from all other nodes. A node can be a junction, reservoir, tank, or air-
chamber. Therefore, no two nodes may have the same ID. However, a node
45
EPANET Input Description
and a link (i.e., pipe, pump, or valve) can have the same ID. The node ID
value can be any string value (up to 40 characters).
When importing (or merging) multiple water distribution network models into a
single network model, MIKE URBAN will check for collisions between identi-
cal node IDs and will automatically assign a new node ID value for any node
being imported that contains the same node ID value as what already exists
in the network model.
DESCRIPTION (optional)
This data entry allows you to enter a description identifying the air-chamber
node being entered. This description can be optionally displayed on the Map
window and in reports generated by the Report Generator.
STATE (optional)
This drop down selection list data entry allows you to define whether the air-
chamber node is imported (i.e existing node was imported from the external
data source), or is inserted, modified, GIS, calibrated or similar. By default,
tank node status is undefined.
Selecting «...» allows you to display the Pressure Zone selection dialog box,
where the appropriate pressure zone ID can be selected.
By default this data entry is k=1.2. See the section titled “Specific Valve Data”
on page 191 for more information.
If the user clicks with the Add Air-Chamber tool in the Map window some-
where else other than a node or a pipe, MIKE URBAN will place a air-cham-
ber node at the clicked location. If the user selects an existing pipe, a
47
EPANET Input Description
program message will be displayed allowing to split the pipe at the cursor
location.
If desired, from the Map window, the user can move an existing air-chamber
using the Move Selected tool. By selecting the air-chamber and holding down
the left mouse button, the air-chamber can be dragged to a new location. As
the tank is dragged, its connecting links (i.e., pipes, pumps, and valves) rub-
ber band along with it.
The project options dialog allows you to specify project configuration informa-
tion for the pipe network system to be modelled. Project configuration data
includes analysis type, friction loss formulation, simulation options, analysis
options, and flowrate units. Alternatively, to display the Project Options dialog
box, select EPANET | Project Options.
The Project Options Editor dialog is used to define the project parameters.
TITLE (Optional)
This data entry allows you to enter a describe the current project. This
description can be output in reports generated by the Report Generator.
EPANET based extended period simulation with the water quality option
such as No Water Quality Analysis (no quality analysis), Chemical con-
centrations (compute chemical concentration), Water Age (compute
water age), Source Tracing (trace flow from a specific node).
Water hammer
UNITS (Mandatory)
This selection data entry allows you to select one of the predefined project
units environments such as CFS (cubic feet per second) or LPS (liter per sec-
ond), for example.
PROPERTIES (Optional)
These data entries allow you to determine the hydraulic and water quality
behaviour of the pipe network should be analysed.
Specific Gravity
This data entry specifies the specific gravity of the fluid at the tempera-
ture condition being simulated. This data entry allows fluids other than
water to be simulated. Gravity is the weight per unit volume of the fluid
being modelled relative to water. Specific gravity is the ratio of the den-
sity of the fluid being modelled to that of water at 4 deg. C. (unitless).
Viscosity
This data entry specifies the kinematic viscosity of the fluid at the tem-
perature condition being simulated. The units of viscosity are ft2/sec (or
m2/sec for SI units). The viscosity is the kinematic viscosity of the fluid
49
EPANET Input Description
Molecular Diffusivity
This data entry specifies the molecular diffusivity of the chemical being
tracked. The diffusivity is the molecular diffusivity of the chemical being
analysed relative to that of chlorine in water. The default value is 1.0. Dif-
fusivity is only used when mass transfer limitations are considered in
pipe wall reactions. A value of 0 will cause MIKE URBAN to ignore mass
transfer limitations.
Emitter Exponent
Power to which pressure is raised when computing the flow through an
emitter device. The textbook value for nozzles and sprinklers is 0.5. This
may not apply to pipe leakage.
CONVERGENCE (Optional)
These data entries allow you to determine the hydraulic and water quality
behaviour of the pipe network should be analysed.
Accuracy
Convergence criterion used to signal that a solution has been found to
the nonlinear equations that govern network hydraulics. Trials end when
the sum of all flow changes divided by the sum of all link flows is less
than this number. Suggested value is 0.001.
Unbalanced System
Action to take if a hydraulic solution is not found within the maximum
number of trials. Choices are STOP to stop the simulation at this point or
CHECKFREQ
This sets the number of solution trials that pass during hydraulic balanc-
ing before the status of pumps, check valves, flow control valves and
pipes connected to tanks are once again updated. The default value is 2,
meaning that status checks are made every other trial. A value equal to
the maximum number of trials would mean that status checks are made
only after a system has converged. (Whenever a status change occurs
the trials must continue since the current solution may not be balanced.)
The frequency of status checks on pressure reducing and pressure sus-
taining valves (PRVs and PSVs) is determined by the DAMPLIMIT option
(see below).
MAXCHECK
MAXCHECK is the number of solution trials after which periodic status
checks on pumps, check valves flow control valves and pipes connected
to tanks are discontinued. Instead, a status check is made only after con-
vergence is achieved. The default value is 10, meaning that after 10 tri-
als, instead of checking status every CHECKFREQ trials, status is
checked only at convergence.
DAMPLIMIT
This is the accuracy value at which solution damping and status checks
on PRVs and PSVs should begin. Damping limits all flow changes to
60% of what they would otherwise be as future trials unfold. The default
is 0 which indicates that no damping should be used and that status
checks on control valves are made at every iteration. Damping might be
needed on networks that have trouble converging, in which case a limit
of 0.01 is suggested.
Theta
Numeric parameter used to centre the high order finite difference
scheme in time.
Alpha
Gravity
Gravity acceleration.
Atmospheric pressure
The pressure of the atmosphere outside the pressurized pipeline system.
51
EPANET Input Description
Temperature
Temperature of water.
Pressure zones are service areas defined by the hydraulic gradeline value of
the sources that supply them. A pressure zone has one or more sources of
supply and may have a set of closed valves that separate it from other pres-
sure zones. The Pressure Zone Editor dialog box is reached by selecting
EPANET | Pressure Zones.
The Pressure Zone Editor dialog is used to define the pressure zones for the
pipe network system.
DESCRIPTION (Optional)
This data entry allows you to enter a description identifying the pressure zone
being defined. This description can be output in reports generated by the
Report Generator.
DEMAND (optional)
This data entry is used to specify the zone demand, which can be used for
automatic demand distribution. To distribute the zone demand to junction
nodes located within the zone, use Demand Distribution dialog.
The user can compute pumping energy and cost using Edit Energy Editor.
The following parameters can be defined in Energy Editor:
53
EPANET Input Description
CONDITION (optional)
The condition data entry is used to describe the contents of the output report
produced from a simulation. The following keywords can be used:
STATUS (YES/NO/FULL)
STATUS determines whether hydraulic status messages are written to the
Report file. If YES is selected the messages will identify those network com-
ponents that change status during each time step of the simulation. If FULL is
selected, then convergence information will also be included from each trial
of each hydraulic analysis. This level of detail is only useful for de-bugging
networks that become hydraulically unstable.
SUMMARY (YES/NO)
SUMMARY determines whether a summary table of number of network com-
ponents and key analysis options is generated. The default is YES.
MESSAGES (YES/NO)
MESSAGES determines whether error and warning messages generated
during a hydraulic/water quality analysis are written to the Report file. The
default is YES.
ENERGY (YES/NO)
ENERGY determines if a table reporting average energy usage and cost for
each pump is provided. The default is NO.
Elevation
Demand
Head
Pressure
Quality
Length
Diameter
Flow
Velocity
Headloss
LinkQuality
LinkStatus
Setting (Roughness for pipes, speed for pumps, pressure/flow setting for
valves)
Reaction (reaction rate)
F-Factor (friction factor)
The default reporting variables are Demand, Head, Pressure, and Quality for
nodes and Flow, Velocity, and Headloss for links. The default precision is two
decimal places.
Example:
The following example reports on nodes N1, N2, N3, and N17 and all links
with velocity above 3.0. The standard node variables (Demand, Head, Pres-
sure, and Quality) are reported on while only Flow, Velocity, and F-Factor
(friction factor) are displayed for links.
NODES N1 N2 N3 N17
LINKS ALL
55
EPANET Input Description
FLOW YES
VELOCITY PRECISION 4
F-FACTOR PRECISION 4
The user can define data curves and their X, Y coordinate points in The
Curves and Relations Editor. The Curves and Relations Editor dialog box is
reached by selecting EPANET | Curves and Relations.
57
EPANET Input Description
Boundary conditions are used to define water level and flow changes during
the water hammer analysis.
There are in principle two types of boundary conditions, namely the piezo-
metric head, H, above a specified datum, e.g., in tanks, and the discharge, Q,
e.g., water demand. Both H and Q are given under selected names as time
series in the Curves and Relations Editor and stored in the database. These
boundary conditions may be assigned to any node in the network. The
boundary set point value for each time step is assigned from the correspond-
ing time series curve. If time step used by water hammer computation is
smaller than appropriate neighbouring values in boundary conditions time
series then linear interpolation is applied. There are nodes of the following
types:
The water hammer analysis ignores the diurnal patterns defined for the
steady state and extended period simulation. Therefore, the boundary condi-
tions using time series must be specified at nodes where the flow or water
head changes need to be defined during the water hammer simulation. The
Boundary Conditions Editor dialog box is reached by selecting EPANET |
Water Hammer | Boundary Conditions.
NODE ID (mandatory)
This data entry is used to define the ID of the boundary node the boundary
condition is being assigned to. Choosing «...» will display the Select Node
selection dialog box from which the user can select the appropriate node type
and ID. Or, choosing «Pick» allows the user to graphically select the node
from the Map window.
DESCRIPTION (optional)
This data entry allows you to enter a description identifying the boundary con-
dition being defined. This description can be output in reports generated by
the Report Generator.
Typically in large network systems, the pipe network is broken up into differ-
ent pressure zones (or distribution zones). Since pressure is related to
ground elevation, a network system covering hilly or mountainous terrain will
have more pressure zones than one covering fairly flat terrain. The section
also discusses how MIKE URBAN defines pressure zones.
59
EPANET Input Description
of Pipe Lengths and the Method of Two Coefficients. This is useful when
assigning the nodal water demand for a large network, since the software will
automatically proportion the total network demand based upon one of these
two methods. These methods are used to mimic the amount of actual
demand along a pipe, based upon the pipe length or a pre-defined demand
coefficient.
Figure 1.25 The Distributed Demands dialog box computes the demand at each
junction node automatically based upon the total demand within a pres-
sure zone or entire network system
A list of the Distributed Demands data entries for Figure 1.25 follows, with a
short description given for each entry.
Note that this total demand represents the total demand regardless into which
multiple junction demand is distributed. Multiple demands are specified in the
Junction Editor or in the Multiple Demand Editor.
The pressure zone must be specified in the provided data entry field. Select-
ing «...» displays the Pressure Zone selection dialog box, where the appropri-
ate pressure zone ID can be selected.
Selecting the Method of Two Coefficients, MIKE URBAN computes the total
water demand assigned to each pipe (which is then split between the starting
and ending nodes) as:
Q k li k 2i
q pi = ---------------------
- (1.4)
k li k 2i
Selecting the Method of Reduced Pipe Lengths, MIKE URBAN computes the
total water demand assigned to each pipe (which is then split between the
starting and ending nodes) as:
Q l i k li
q pi = -----------------
- (1.5)
k li l i
Q l i k Di
q pi = -------------------
- (1.6)
k Di l i
where:
qpi =Total water demand applied to the pipe, split between the two end nodes.
li = Pipe length
61
EPANET Input Description
These demand coefficients are defined for each pipe using the Pipe Editor, as
described on page -18. The computed demands, which are assigned once
«Compute» is selected, are stored at each individual node. These demands
are stored in the Junction Editor, as described on page -12. Selecting
«Reset» causes all of the nodal demand entries to be set to zero for the
entire network, but it will leave the additional demand entries that are defined
unaffected.
This option will only assign demand to nodes with Demand Coefficient
applied (different from 0 or NULL). In the case of an equal distribution, the
node demand coefficients have to be equal and different from zero.
Note that this total demand represents the total demand regardless into which
multiple junction demand is distributed. Multiple demands are specified in the
Junction Editor or in the Multiple Demand Editor.
Unchecking this box applies the specified water demand to the entire water
distribution network.
The pressure zone must be specified in the provided data entry field. Select-
ing «...» displays the Pressure Zone selection dialog box, where the appropri-
ate pressure zone ID can be selected.
Q
q ni = ---- (1.7)
N
where:
QT
-C
q i = ----------------- (1.8)
i
i C
i = 0,n
where:
qi = node demand
Ci = demand coefficient
63
EPANET Input Description
Junction node demands can be edited either within the junction editor for
each particular junction node or in the Multiple Demand editor, which allows
the user to display and edit all multiple demands. In addition to that, MIKE
URBAN provides another way of developing junction node demands, which is
based on geocoding the consumption data to the appropriate junctions or
pipes and aggregating their set point demand values to the junction
demands. This allows simplifies the demand development process and
allows you, for example, to import consumption data from the consumption
database systems and geocoded it based on its X, Y geographical coordi-
nates.
Figure 1.26 Demand Allocation Editor dialog box is used to store and edit the con-
sumption data defined as consumption points and link it to he appropri-
ate network junctions or pipes
A list of the Demand Allocation Editor data entries for Figure 1.26 follows,
with a short description given for each entry.
REFERENCE ID (optional)
This data entry is used to specify an ID which uniquely identifies the demand
point. The reference ID acts as a unique lookup key that identifies the
demand point from all other demand points. The reference ID value can be
any string value (up to 40 characters). It is recommended that this reference
ID corresponds to the asset ID, which uniquely identifies the demand point in
the customer information or billing database system.
X, Y COORDINATE (mandatory)
The X and Y data entries are used to define the physical (map) location of the
demand point, in units of ft. or m.
DEMAND (mandatory)
This data entry is used to specify the demand value, which will be then used
in the process of the demand aggregation. This demand set point value can
be imported from the external database systems (such as CIS Customer
Information System, for example) or it can be developed from the minimum,
average, or maximum demand values.
65
EPANET Input Description
DESCRIPTION (optional)
This data entry allows you to enter a description identifying the consumption
point being defined. This description can be output in reports generated by
the Report Generator.
ADDRESS (optional)
This data entry allows you to enter an address identifying the consumption
point being defined. This field can be output in reports generated by the
Report Generator.
OWNER (optional)
This data entry allows you to enter an owner name identifying the consump-
tion point being defined. This field can be output in reports generated by the
Report Generator.
EQUIVALENT PERSON(optional)
This data entry is used to specify the equivalent person value for each
demand point. The EP value can be used when calculating the actual junction
demand, which will be then used in the process of demand aggregation.
Geocoding
Demand geocoding is an automatic process, where each demand point is
linked to the appropriate junction node or a pipe (such as nearest node, for
example). The user can also specify the linked junction or pipe ID manually
by editing or importing from external database systems, such as GIS or CIS.
There are several different geocoding methods available. The geocoding is
done from a wizard. Press the ‘Geocode’ button on the Demand Allocation
dialog to activate the wizard.
To nearest node
To node by Pipe ID
The nearest pipe is already specified and the nearest of the pipe begging
or ending nodes is used as a demand node.
To nearest pipe
67
EPANET Input Description
Aggregation
Demand aggregation is used to develop node demands for junction nodes
based on the geocoded node or link ID/
works for demand allocations connected to pipes. Note that this does not
add to the Multiple Demands table (as is done for allocation assigned to
nodes), but sets the demand coefficients of the pipe.
TARGET (mandatory)
This data entry allows you to specify the multiple demand category, which will
be used as a target demand for the aggregated demand. If the multiple
demand with the specified category does not exist, the program will create it
and it will override the existing values in case that such demand category
already exists for each node used in the demand distribution.
Multiple demands can be edited either within the junction editor for each par-
ticular junction node or in the Multiple Demand editor, which allows the user
to display and edit all multiple demands. The Multiple Demand Editor dialog
box is reached by selecting EPANET | Multiple Demands.
69
EPANET Input Description
JUNCTION ID (mandatory)
The Junction ID identifies the selected Junctions which multiple demands are
assigned to.
DEMAND (mandatory)
The demand field shows all the demand values that are assigned to junctions
with multiple demands. The demand values must be manually entered in the
demand field.
DESCRIPTION (optional)
This data entry allows you to enter a description identifying the multiple
demand being defined. This description can be output in reports generated by
the Report Generator.
MIKE URBAN can generate statistical information for junction node demands.
Demand statistics is generated for each pressure zone as well as for the
complete network. Additionally, demand statistics dialog box allows the user
to redistribute node demands by changing the calculated statistical results.
The Demand Statistics dialog box is reached by selecting EPANET |
Demand Statistics.
A list of the Demand Statistics information for Figure 1.31 follows, with a short
description given for each entry.
Type
This data entry is used to distinguish between data (demand statistics for
the selected category), zone (demand statistics for the pressure zone),
and network (demand statistics for the whole network).
Pressure Zone
Pressure zone ID identifies the pressure zone for which the demand sta-
tistics is generated.
Category
This data entry identifies the category within the current pressure zone
for which the demand statistics is generated.
Minimum Demand
This data entry represents the minimum demand per category per pres-
sure zone. The minimum demand is calculated as minimum demand at
junction nodes at specific time level
Maximum Demand
This data entry represents the maximum demand per category per pres-
sure zone. The maximum demand is calculated as maximum demand at
junction nodes at specific time level
71
EPANET Input Description
Average Demand
This data entry represents the average demand per category per pres-
sure zone. The average demand is calculated as average demand at
junction nodes during the entire simulation period.
Sum Demand
This data entry represents the total demand per category per pressure
zone. The total demand is calculated as total demand at junction nodes
during the entire simulation period.
Category Type
This option data entry allows you to specify the category type such as
fixed or scaled, which will be used in the demand redistribute process.
Junction node demand with fixed demand category will not be scaled
during the demand redistribution.
REFRESH
Select Refresh to re-generate the demand statistics.
REDISTRIBUTE
Select Redistribute to redistribute node demands based on the new values of
Average or Total demand (zone or network). This powerful feature provides
the user with the option of specifying the new zone or a network demand from
within the Demand Statistics window and redistribute the node demand
accordingly. The process of the demand distribution is based on using the
existing node demands as coefficients - weights to calculate the new demand
values.
Figure 1.32 Demand statistics dialog box: Redistribute allows you to define the new
average or total zone (or network) demand and redistribute it to junction
nodes
Controls (sometimes called switches) allow the settings (or status) of various
network components to change at a particular time, when a specific pressure
is reached, or when a specific water tank level occurs. For example, a pump
could be turned off once a particular water tank level has been reached.
The following sections describe available editors that allow you to define an
extended period simulation.
73
EPANET Input Description
Typically during an extended period simulation, the pipes, pumps, and valves
(links) contained in a network will change their status (i.e., open or close) as
storage tanks fill and empty and pressures change throughout the network
system. Also, for a steady state simulation, network components may change
their state as the analysis model iterates to a valid solution.
Using the Control Editor, as shown in Figure 1.33, the user can specify the
operational controls for the pipe network system on a one-on-one component
basis. The Control Editor can be used for both extended period simulations
and steady state simulations.
The following situations exemplify the types of operational controls that can
be specified:
A pipe can be opened at a given time (based upon the beginning of the
network simulation). This type of operational control has no effect in a
steady state simulation.
A pump can be turned on or off depending on the water level in a speci-
fied tank.
A pump’s baseline operating speed can be adjusted using a multiplier,
thereby increasing or decreasing the pump’s output.
A valve can be opened or closed based upon the pressure in an adjacent
node.
The initial status of the network’s pipes, pumps, and valves default to an open
status, and can be changed to a closed status using the Control Editor based
upon some condition. Additionally, pumps have a default pump speed ratio of
1.0 (indicating that they are operating on their original characteristic opera-
tional curves), but which can be modified based upon a particular condition.
The Control Editor dialog box, as shown in Figure 1.33, is used to define the
operational controls for running the pipe network system. The Control Editor
dialog box is reached by selecting EPANET | Extended Period|Simple Con-
trols ....
In order to export simple controls into ASCII file, open Simple Controls Editor,
select Advanced and select "Export to ASCII file…". This will allow you to
define the name of the ASCII file and export the content of the Simple Control
Editor into the selected ASCII file. You can edit the controls within the ASCII
file and import them back to MIKE URBAN by selecting Load from ASCII file.
This is convenient in cases when you use Excel or other tools to create the
list of controls for the model.
Figure 1.33 The Control Editor allows the user to define the network operational
controls
A list of the Control Editor data entries for Figure 1.33 follows, with a short
description given for each entry.
LINK ID (mandatory)
This data entry is used to define the ID of the link to control. Choosing «...»
will display the Select Link selection dialog box from which the user can
select the appropriate link type and ID. Or, choosing «Pick» allows the user to
graphically select the link from the Map window.
DESCRIPTION (optional)
This data entry allows you to enter a description identifying the control rule
being defined. This description can be optionally included in reports gener-
ated in MIKE URBAN.
SETTING (mandatory)
This radio button selection entry is used to specify the OPEN or CLOSED sta-
tus of the link being controlled, or VALUE which then allows the user to spec-
ify a multiplier to a pump’s operating speed if a pump is being controlled. If
VALUE is selected, then an additional data entry field will be presented to
allow the user to specify the pump speed ratio.
CONDITION (mandatory)
This radio button selection entry is used to specify the control condition that
then applies the operational rule onto the link being controlled.
75
EPANET Input Description
If the user selects either IF NODE BELOW or IF NODE ABOVE, then a Con-
trol Node ID and a Control Level must be specified. Choosing «...» will dis-
play the Select Node selection dialog box from which the user can select the
appropriate node type and ID. Or, choosing «Pick» allows the user to graphi-
cally select the node from the Map window. Note that reservoirs are not
allowed to be selected as a Control Node type.
If the user selects AT TIME, then a trigger time (since the start of the simula-
tion) must be specified in the adjacent data entry field and a time unit
selected from the pull-down selection list. Note that this type of condition has
no effect for steady state simulations, although can be defined should the
user want to perform an extended period simulation at a later time.
Multiple Controls
During a simulation, a link (pipe, pump, or valve) can be operated by more
than one operational control rule. For example, a link can be opened at a
given time and then closed if a given nodal pressure is exceeded.
Not that it is necessary to re-apply the link setting (such as control valve set-
ting or pump speed) after the link is re-open by the control rules.
77
EPANET Input Description
Figure 1.34 The Control Editor allows the user to define the network operational
controls Rule Based Controls Editor Dialog Box
In order to rule based controls into ASCII file, open Rule Based Controls Edi-
tor, select Advanced and select "Export to ASCII file…". This will allow you to
define the name of the ASCII file and export the content of the Rule Based
Controls Editor into the selected ASCII file. You can edit the rules within the
ASCII file and import them back to MIKE URBAN by selecting Load from
ASCII file. This is convenient in cases when you use Excel or other tools to
create the list of rules for the model.
where
JUNCTION PIPE
RESERVOIR PUMP
TANK VALVE
DEMAND
HEAD
PRESSURE
LEVEL
FLOW
79
EPANET Input Description
= IS
<> NOT
< BELOW
> ABOVE
<=> =
where
The default values for the B (multiplier) and C (increment) are B=1, C=0 and
they do not need to be provided.
The setting value is calculated as: Setting Value = Setting Value (link A) * B +
C
Example 1:
IF SYSTEM CLOCKTIME = 8 AM
The pump 3 setting will be set equal to the settings of the pump 4 (multiplier =
1 and increment = 0) at time 8 am.
Example 2:
IF SYSTEM CLOCKTIME = 8 AM
The pump 3 setting will be higher by 10% than the settings of the pump 4
(multiplier = 1.10 and increment = 0) at time 8 am.
Example 3:
IF SYSTEM TIME >= 12 AM
The valve 10 setting (PRV setting, for example) will be lower by 10 (pressure
units) than the settings of the valve 20 (multiplier = 1 and increment = -10) at
any time (12 AM is the simulation start).
Example 4:
IF SYSTEM TIME >= 12 AM
The pump 3 status (OPEN, CLOSED) will be set equal to the status of the
pump 4 at any time (12 AM is the simulation start).
Example 5:
IF PUMP 3 SETTING LIKE PUMP 4
THEN …
81
EPANET Input Description
Control of valves
This set of rules opens and closes a Valve based on different water levels in
tanks.
Control of pumps
This set of rules opens and closes a pump based on different water levels in
tanks.
MIKE URBAN assumes that water usage rates, external water supply rates,
and constituent source concentrations (for water quality analysis) at nodes
remain constant over a fixed period of time, but that these quantities can
change from one time period to another. The default time period interval is
1 hour, but this can be set to any value. The value of any of these quantities
in a time period equals a baseline value multiplied by a time pattern factor for
that period.
Diurnal Patterns
The diurnal profiles are used to define a series of multipliers (multiplication
factors applied to a baseline value of junction node demand, constituent
source concentration, storage tank water level). The duration of such diurnal
profile curve is 1 day i.e. 24 hours. The Diurnal Profiles Editor dialog box is
reached by selecting EPANET | Extended Period|Repetitive Profiles |
Diurnal Patterns.
83
EPANET Input Description
Figure 1.35 Diurnal Profiles Editor, the Multipliers grid allows the user to define
demand and/or constituent concentration multipliers for the current pat-
tern ID
PATTERN ID (mandatory)
This data entry is used to specify the ID of the component being defined. The
pattern ID value can be any string value (up to 40 characters). There is no
limit to the number of demand patterns that can be defined.
CATEGORY (optional)
This data entry allows you to enter a description identifying the demand pat-
tern being defined. This description can be optionally included in reports.
DESCRIPTION (optional)
This data entry allows the user to enter a category that further define the
demand pattern. For example, a demand might have the description of resi-
dential, and a category of either high density, medium-density or low density
to further define what is meant by residential. This description can optionally
be included in reports generated by the Report Generator.
MULTIPLIERS (mandatory)
Selecting «Define» will display the Multipliers dialog box, as displayed in
Figure 1.35. This dialog box allows you to define the baseline demand factors
and constituent concentration levels (multipliers) for the current pattern ID.
The complete pattern is then applied to the baseline demand (or source con-
centration) specified at each node (that corresponds with this pattern ID) over
the extended period simulation.
Within each time period, the demand (or source concentration) remains con-
stant—at a level equal to the multiplier times the baseline value.
There is no limit on the number of multipliers that can be defined for a pattern.
However, since the diurnal profile duration is exactly 24 hours, the number of
multipliers is determined by the pattern time step. The typical time period
length is 1 hour, 24 multipliers are required to completely define the pattern
for the entire simulation. For more information on defining the time step inter-
val and analysis period, see the section titled Time Setting Editor which fol-
lows this section.
DISTRIBUTE
Select Distribute to let the program calculate the number of multipliers to fill
entirely the diurnal pattern duration of 24 hours with the current pattern time
step.
GRAPH
It is possible to create a GRAPH for pattern multipliers directly in the PAT-
TERN EDITOR window, not only in the MULTIPLIERS window. This function
allows the user to display the pattern graph quickly.
85
EPANET Input Description
It is possible to import and export the diurnal pattern multipliers from an ASCII
file. Select Advanced - Import or Advanced - Export to import or export the
diurnal pattern multipliers from an ASCII text file.
Cyclic Profiles
The cyclic pattern editor is used to combine 1 or more diurnal profiles into one
cyclic pattern. This allows the user to compose the cyclic patterns of two diur-
nal profiles one for e.g. working days and one for the weekend. If the total
duration of the defined cyclic pattern is less than the total duration of the sim-
ulation, then the defined pattern will be repeated. For example, a 5 day simu-
lation whose hourly demands repeat themselves on a daily cycle, only
24 multipliers would need to be defined for the cyclic pattern diurnal profile.
It is possible to display the graph for the cyclic profile, consisting of several
diurnal patterns. In order to display the graph, open Cyclic Profile Editor,
select the profile and click "Graph".
For an extended period simulation, a simulation duration and time step must
be specified. Using the Time Editor, as shown in Figure 1.38, the user can
specify the extended period simulation time parameters. Only those time
parameters that differ from the default values must be specified. The default
time units are hours.
Figure 1.38 The Time Settings Editor allows the user to define the extended period
simulation time parameters
A list of the Time Editor data entries for Figure 1.38 follows, with a short
description given for each entry.
87
EPANET Input Description
STATISTICS (optional)
Type of statistical processing used to summarize the results of an extended
period simulation. Choices are:
Statistical processing is applied to all node and link results obtained between
the Report Start Time and the Total Duration
Time Units
Note that seconds, minutes, hours, and days can be used as time units, but
that hours are generally used. It is recommended that the same time units be
used for all entries within the Time Settings Editor—that you do not mix time
units (i.e., hours and seconds).
ute the pattern multipliers. If the user selects yes, the pattern number of pat-
tern multipliers will be adjusted to fit the diurnal curve duration of 24 hrs.
Figure 1.39 When changing the pattern time step, MIKE URBAN will ask whether to
redistribute the pattern multipliers
The first step in defining a water quality simulation is to select the type of
water quality analysis to be performed. This is accomplished using the Pro-
ject Options dialog box Project Type tab, as shown in Figure 1.40. Note that
MIKE URBAN can only perform one type of water quality analysis during a
simulation.
Figure 1.40 The Project Options dialog box allows the user to define the type of
water quality analysis to be performed
89
EPANET Input Description
A list of water quality options are available from the Project Options dialog
box is shown in Figure 1.40, along with a short description given for each
option.
CHEMICAL CONCENTRATIONS
This option button entry is used to select that a chemical concentration analy-
sis is to be performed. MIKE URBAN will report the concentration of the spec-
ified chemical at each time step for each network node. This type of analysis
is typically performed to determine the amount of chlorine in the network to
maintain safe drinking water standards.
WATER AGE
This option button entry is used to select that a water age analysis is to be
performed. MIKE URBAN will report the water age at each time step for any
network node, assuming an initial age of zero at the start of the simulation.
This type of analysis is typically performed to determine dead ends (locations
of stagnant water) within a pipe network design.
SOURCE TRACING
This option button entry is used to select that a source tracing analysis is to
be performed. MIKE URBAN will report, for each time step, the percentage of
water reaching each node from a selected source (origin) node. This type of
analysis is typically performed for constituent tracking.
The Project Options dialog box Properties tab, as shown in Figure 1.41,
allows you to specify the analysis parameters used to perform a water quality
simulation.
Figure 1.41 The Project Options dialog box allows the user to define the analysis
parameters used in a water quality simulation
A list of the water quality analysis parameters from the Project Options dialog
box is shown in Figure 1.41, along with a short description given for each
parameter. These water quality modelling parameters are used when model-
ling the pipe-wall reaction mechanism.
SPECIFIC GRAVITY
This data entry specifies the specific gravity of the fluid at the temperature
condition being simulated. This data entry allows fluids other than water to be
simulated. Gravity is the weight per unit volume of the fluid being modelled
relative to water. Specific gravity is the ratio of the density of the fluid being
modelled to that of water at 4 deg. C. (unitless).
KINEMATIC VISCOSITY
This data entry specifies the kinematic viscosity of the fluid at the temperature
condition being simulated. The units of viscosity are ft2/sec (or m2/sec for SI
units). The viscosity is the kinematic viscosity of the fluid being modelled rela-
tive to that of water at 20 deg. C (1.0 centistoke). The default value is 1.0.
MOLECULAR DIFFUSIVITY
This data entry specifies the molecular diffusivity of the chemical being
tracked. The diffusivity is the molecular diffusivity of the chemical being ana-
lyzed relative to that of chlorine in water. The default value is 1.0. Diffusivity is
only used when mass transfer limitations are considered in pipe wall reac-
tions. A value of 0 will cause MIKE URBAN to ignore mass transfer limita-
tions.
The initial water quality at the start of a simulation can be assigned to individ-
ual nodes or to groups of nodes. The initial water quality can represent one of
the following:
91
EPANET Input Description
By default, all nodes are assigned with an initial water quality of zero. The Ini-
tial Water Quality Editor, as shown in Figure 1.42, is used to assign the initial
water quality levels for the pipe network system. The Initial Water Quality Edi-
tor dialog box is reached by selecting EPANET | Water Quality |
Initial Water Quality.
Figure 1.42 The Initial Water Quality Editor is used to define the initial water quality
conditions of the pipe network system
A list of the Initial Water Quality Editor data entries for Figure 1.42 follows,
with a short description given for each entry.
NODE ID
This data entry is used to specify an ID which uniquely identifies a node in
which the initial water quality is being specified for.
A new entry is automatically inserted into the list by pressing «Insert». Choos-
ing «...» will display the Select Node selection dialog box from which the user
can select the appropriate node Or, choosing «Pick» allows the user to
graphically select the node from the Map window.
QUALITY
This data entry is used to specify the initial water quality at the node (or group
of nodes) being defined.
If a water age water quality simulation is being performed, then this entry
denotes the initial age of the water in hours.
If a source tracing water quality simulation is being performed, then this entry
denotes the initial percentage of water from the source (origin) node in per-
cent.
The Point Constituent Source Editor, as shown in Figure 1.43, allows you to
specify at which nodes an external chemical constituent enters the network
system. At least one node in the network must be specified as a point source
of chemical constituent when performing a chemical concentration analysis.
Figure 1.43 The Point Constituent Source Editor is used to specify at which nodes
an external chemical constituent enters the pipe network system
A list of the Point Constituent Source Editor data entries for Figure 1.43 fol-
lows, with a short description given for each entry.
NODE TYPE
This pull-down selection list allows the user to select what type of node (i.e.,
junction, reservoir, or tank) the point constituent is being specified for.
NODE ID
This data entry is used to define the ID of the node the point constituent is
being assigned to. Choosing «...» will display the Select Node selection dia-
log box from which the user can select the appropriate node type and ID. Or,
choosing «Pick» allows the user to graphically select the node from the Map
window.
CONCENTRATION
This data entry is used to specify the baseline concentration (in mg/liter) of
the constituent entering the node as an external source.
93
EPANET Input Description
For a junction node, if there is no external inflow assigned to the node (such
as a well—denoted by a negative demand), then the resulting water quality at
the node always equals the specified concentration. This allows the user to
simulate chlorine booster stations at a node, such as is used in satellite treat-
ment in a network.
SOURCE TYPE
Water quality sources are nodes where the quality of external flow entering
the network is specified. They can represent the main treatment works, a
well-head or satellite treatment facility, or an unwanted contaminant intrusion.
Source quality can be made to vary over time by assigning it a time pattern.
MIKE URBAN can model the following types of sources:
A mass booster source adds a fixed mass flow to that entering the
node from other points in the network.
The concentration-type source is best used for nodes that represent source
water supplies or treatment works (e.g., reservoirs or nodes assigned a nega-
tive demand). The booster-type source is best used to model direct injection
of a tracer or additional disinfectant into the network or to model a contami-
nant intrusion.
CYCLIC PROFILE ID
This data entry allows you to define the ID of the constituent pattern to be
applied to the specified baseline concentration entering the node. If a pattern
ID is omitted for the specified source node, then there is no variation in the
source strength of the constituent.
Selecting «Button» allows you to display the Select Pattern selection dialog
box, where the appropriate pattern ID can be selected.
The Reaction Rate Editor, as shown in Figure 1.44, allows you to specify the
rate at which a constituent decays (or grows) by reaction as the constituent
travels through the pipe network. Reaction rates can be defined at a global
level where the same reaction rate applies to the entire network, and at a
local level allowing the user to define reaction rates based upon a particular
pipe type, etc. Note that locally defined reaction rates override globally
defined reaction rates.
The Reaction Rate Editor dialog box is reached by selecting EPANET | Water
Quality | Reaction Rates. Note that a Chemical Concentrations water quality
simulation must be specified in the Projects Options dialog box (see
Figure 1.41) in order access the Reaction Rate Editor dialog box.
Global Settings
Figure 1.44 The Reaction Rate Editor is used to specify constituent reaction rates at
both a global and local level
A list of the Reaction Rate Editor data entries for Figure 1.44 follows, with a
short description given for each entry.
95
EPANET Input Description
One method that can be used to compare the relative magnitude of the pipe
wall reaction rate with the bulk reaction rate is to divide the pipe wall reaction
rate coefficient by the hydraulic radius of the pipe (i.e., 1/2 the pipe radius).
The resulting quantity will have the same units as the bulk reaction rate coef-
ficient, days–1.
LIMITING POTENTIAL
LIMITING POTENTIAL specifies that reaction rates are proportional to the
difference between the current concentration and some limiting potential
value.
ROUGHNESS CORRELATION
ROUGHNESS CORRELATION will make all default pipe wall reaction coeffi-
cients be related to pipe roughness in the following manner
Hazen-Williams F/C
Darcy-Weisbach -F/log(e/D)
Chezy-Manning F*n
Note
Remember to use negative signs on all reaction coefficients that are to repre-
sent constituent decay (e.g., chlorine decay). Otherwise, a positive value
denotes constituent growth.
The Trace Node dialog box, as shown in Figure 1.45, allows you to track over
time what percent of water reaching any node in the network had its origin
from a specified node (i.e., junction, tank, or reservoir). Source tracing is a
useful tool for analyzing a network distribution system that draws water from
two or more different raw water supplies. It can be used to show to what
degree water from a given source blends with that from other sources, and
how the spatial pattern of this blending changes over time.
Figure 1.45 The Trace Node dialog box is used to specify a single node that acts
like a tracer in determining what percent of its water reaches any other
node in the network
A list of the Trace Node dialog box data entries for Figure 1.45 follows, with a
short description given for each entry.
97
EPANET Input Description
NODE TYPE
This pull-down selection list allows the user to select what type of node (i.e.,
junction, reservoir, or tank) the trace node is being specified for.
NODE ID
This data entry is used to define the ID of the node the trace node is being
assigned to. Choosing «...» will display the Select Node selection dialog box
from which the user can select the appropriate node type and ID. Or, choos-
ing «Pick» allows the user to graphically select the node from the Map win-
dow.
1.6 Tools
The following section describes additional capabilities in MIKE URBAN that
allow the user to work effectively with the program.
You can transfer data between the asset and the model database of MIKE
URBAN. See the section titled Transferring data between ‘Asset’ and ‘Model’
in the MIKE URBAN Model Manager book, chapter Modelling Water Distribu-
tion networks for more details on Send Data.
Included with MIKE URBAN are predefined engineering tables that include
specified pipe roughness coefficients and minor loss coefficients, as shown in
the Engineering Tables dialog box as shown in Figure 1.47. These values are
used as lookup coefficient values when defining the pipe and valves that
make up the network system. The Engineering Tables dialog box is reached
by selecting Tools | Engineering Tables. Note that additional coefficients
can be added to the provided engineering table values using this dialog box.
Figure 1.46 The Engineering Tables define lookup table values for minor loss coeffi-
cients
Figure 1.47 The Engineering Tables define lookup table values for pipe roughness
coefficients
Roughness coefficient and minor loss coefficient lookup table dialog boxes
are displayed from either the Pipe Editor and Valve Editor by clicking on the
«...» button, adjacent to the Roughness and Minor Loss data entries.
99
EPANET Input Description
The Convert Asset Valves (or Pumps) to Model Valves (or Pumps) dialog box
is tailored to convert asset valves/pumps (point feature class mwa_-
valve/mwa_pump) to the model valves/pumps (polyline feature class mw_-
valve/mw_pump). This tool is available when running MIKE URBAN in the
WD Water Distribution mode.
The Convert Asset Valves (or Pumps) to Model Valves (or Pumps) dialog box
is reached by selecting Tools | Toolbox | Convert Asset Valves (Pumps) to
Model Valves (Pumps). This tool will start converting asset valves to model
valves by finding the nearest model pipe for each asset valve (or selected
asset valve), which is located within a snap tolerance radius to the nearest
model pipe. In order to import GIS valves to the asset valves (mwa_valve fea-
ture class) use File Import and Export dialog box; it is also possible to create
asset valves by graphical editing.
A list of the Convert Asset Valves to Model Valves dialog box data fields for
Figure 34.57 follows, with a short description given for each entry.
Figure 1.49 Conversion of the asset valve to the model valve; note that the pipe is
split and two junction nodes are created nearby the asset valve loca-
tion.
It is very common that you import the model based on the CAD files, for
example, and that the pipe attributes (imported as CAD labels) were not
available for every pipe. Similar situation will also happen when you start edit-
ing pipes manually and you do not really want to assign the pipe attributes
one by one.
MIKE URBAN provides a tool, which allows you to select pipes from the Map
and duplicate their attributes. Pipe material, as illustrated in the example
below, is automatically duplicated to all pipes, which are adjacent to the
selected pipe(s) until "T" or other complex junction exists.
Figure 1.50 Step1: Select pipes which attribute(s) you want to duplicate.
Now, select attributes you would like to duplicate, such as pipe material,
diameter, and others. Click "OK" where ready.
101
EPANET Input Description
"Generate tag" tool allows you to generate unique identifier "Tag" for every
stretch of pipes within the model. This is very useful for generating unique
"GIS" tag for such "chain of pipes", which are linked together by just 1 single
node.
Figure 1.53 llustration: GIS tag is generated for each stretch of pipes
103
EPANET Input Description
105
106 MIKE URBAN - © DHI
2 Using the Model
MIKE URBAN™ is the most advanced, powerful, and comprehensive water
distribution modelling package available. MIKE URBAN WD (Water Distribu-
tion) can analyse an entire water distribution system, or selected portions,
under steady state, extended period and water hammer flow conditions, with
water quality analysis if needed. You can enter model data suited for your
choice of numerical engine EPANET (public domain EPANET engine devel-
oped by EPA, US), Water hammer engine (developed by DHI and a part of
the WD-Tools) and/or WD-Tools (DHI suit of tools such as fire flow analysis,
pipe roughness calibration, variable speed drive pump, extended controls
based on the modified EPANET engine). You can also enter ‘Asset WD’ data
suited for a water supply and water distribution system. The public domain
EPANET engine can be run from the Model Manager, while running WD-
Tools requires you to have the license for the MIKE URBAN Water Distribu-
tion system with WD-Tools.
Many times, existing water distribution systems do not have a detailed map
that can be used to graphically construct a network model. For these situa-
tions, MIKE URBAN allows the user to develop a model by simply defining
water distribution components (i.e., pipes, junction nodes, pumps, values,
tanks, and reservoirs) in interactive, easy-to-use dialog boxes. This allows
the user to define a model when an accurate map is not available for the pipe
network model.
After the pipe network has been defined, a hydraulic analysis of the network
can be performed. Detailed reports can be generated from the analysis
results and printed out. In addition, if a graphical representation of the pipe
network has been created, the computed pipe flowrates and junction node
pressures can be displayed and printed out.
107
Using the Model
2.1 Overview
MIKE URBAN Water Distribution (WD) consists of the following ‘modules’.
2.1.1 EPANET
2.1.2 WD-Tools
Specify a design fire flow rate and compute the available fire flow pres-
sure.
Specify a design fire flow pressure and compute the available fire flow
rate.
PID Control
The PID (Proportional-Integral-Differential) Control module provides
advanced Real-Time Control simulation capabilities for water supply and
water distribution systems. It permits description of various controllable
devices and makes the definition of complex operational logic for interde-
pendent regulators fully transparent and time efficient.
109
Using the Model
Water Hammer
Water hammer simulates transient (unsteady) flow in any fully pressurized
system carrying liquids. MIKE URBAN Water Hammer provides a cost effec-
tive tool for engineers seeking fast answers to questions about rapid opera-
tion of piping systems. Water hammer is based on the high-order implicit
scheme solving the continuity and momentum equation using the finite differ-
ence method. The initial conditions are modeled using MIKE URBAN
EPANET based steady state flow module.
It is generally more efficient to gather and organize the data required to define
your pipe network model before you begin input of data into MIKE URBAN
Water Distribution.
Figure 2.1 The network components that make up a typical water distribution
model
For large networks, collecting the input data to define a pipe network model
can involve considerable effort. Also, entering this input data into the com-
puter can, at times, become tedious. This process can be prolonged by errors
introduced while entering the data. To reduce the possibility of error, the fol-
lowing procedure is recommended for preparing and entering the input data
and performing the network analysis.
111
Using the Model
2. Prepare the input data for each pipe and junction node. The data for
these components should be complete. For example, pipe data includes
the pipe length, diameter, roughness coefficient, and minor loss coeffi-
cient for each pipe in the network. If the roughness coefficient is not
available, the pipe material and age can be used to interpolate a rough-
ness coefficient value using the roughness coefficient tables provided in
MIKE URBAN Water Distribution. Junction node data includes the water
demand and node elevation for each node in the network.
3. If there are control components in the network, prepare the input data for
each component. Control components include pumps, check valves, reg-
ulating valves, sustaining valves, flow control valves, and storage tanks.
4. Enter the input data which defines the network model using the MIKE
URBAN Water Distribution data input dialog boxes. The program will
report an error if invalid data is entered. This immediate checking of input
data helps prevent errors from being introduced into the pipe network
model.
5. After the input data has been entered, MIKE URBAN Water Distribution
can perform a complete check of the input data and a geometric verifica-
tion of the network connectivity. In addition, it can print a summary of the
input data to allow the user to verify the entered data. Many of the input
data dialog boxes, such as the and Pipe Editor dialog boxes, can be
placed side by side to allow the user to more easily verify model consist-
ency.
6. After the input data has been verified, the pipe network can then be ana-
lysed. Analysis results should be carefully examined to make certain that
the input data accurately defines the model and that the analysis results
appear reasonable. Of particular importance are the junction node pres-
sures. These values should lie within a reasonable range.
7. Calibrate the network model to make certain that the analysis provides
accurate modelling results. For more information on how to calibrate
your network model, see the section titled “Model Calibration” on
page 243.
Project Data
General data describing the pipe network model, such as project and output
specification settings, is categorized as project data.
Network Data
Network data consists of all components that make up the network model,
such as pipes and junction nodes. This subsection lists all components that
can be used in a steady state analysis.
113
Using the Model
Pipe Data
Pipe data includes the physical characteristics of each pipe, such as pipe
length, diameter, roughness coefficient, and minor loss due to fittings along
the pipe. In addition, check valves can be defined along a pipe to prevent flow
reversal. The pipe data provides the basic information for the hydraulic analy-
sis, and should be prepared carefully.
A text description can be defined for each pipe. Because pipes are associ-
ated with the roads and streets that they lie under, a pipe is generally labeled
by its associated street name.
Junction Data
Junction data describes the physical characteristics of each junction node,
such as external water demand and node elevation.
A text description can be defined for each junction node. Junction nodes are
often labelled with names corresponding to the buildings they are located
near.
Reservoir Data
Reservoir data specifies the reservoir's water surface elevation. In addition, a
text description can be defined for each reservoir.
Emitters Data
Emitters are devices associated with junctions that model the flow through a
nozzle or orifice. In these situations the demand (i.e. the flow rate through the
emitter) varies in proportion to the pressure at the junction raised to some
power. The constant of proportionality is termed the "discharge coefficient".
For nozzles and sprinkler heads the exponent on pressure is 0.5 and the
manufacturer usually states the value of the discharge coefficient as the flow
rate in gpm through the device at a 1 psi pressure drop. Emitters are used to
model flow through sprinkler systems and irrigation networks. They can also
be used to simulate leakage in a pipe connected to the junction (if a dis-
charge coefficient and pressure exponent for the leaking crack or joint can be
estimated) and compute a fire flow at the junction (the flow available at some
minimum residual pressure).
Pump Data
Pump data is defined by selecting one of four pump types: constant power,
standard pump curve, custom pump curve with no extended flow range, or
custom pump curve with extended flow range. Each of these pump types
Valve Data
Valve data includes the valve type, the HGL setting of the valve, and the ref-
erence junction node. In MIKE URBAN Water Distribution, there are six valve
types: pressure reducing valves (PRV), pressure sustaining valves (PSV),
pressure breaker valves (PBV), flow control valves (FCV), throttle control
valves (TCV) and general purpose valves (GPV). The GPV allows the user to
define the flow-headloss relationship through the valve. Most pipe networks
use valves to regulate flowrates and pressures. Valves greatly affect the anal-
ysis results. Therefore, valves should be defined carefully. A text description
can be defined for each valve. A check valve is not defined with the valve
data, but is defined with the pipe data.
Consumption Data
Consummation data includes junction node demands and time patterns to
additional categories of water users at a junction as well as demand alloca-
tion data, and pipe demand coefficients. For more information on how to
develop node demands, see the section titled “Distributed Demands” on
page 119.
Control Data
Control data allows pipes, pumps, and valve settings to change at specific
times or when specific pressures or tank levels are reached in the network. A
text description can be defined for each control.
115
Using the Model
Most data used to describe a water distribution model is defined using inter-
active dialog boxes, an example of which is shown in Figure 2.2. These dia-
log boxes allow the user to quickly comprehend what data input requirements
are needed to define the model.
Figure 2.2 This interactive dialog box illustrates how data is entered to define a
water distribution network model
The MIKE URBAN Water Distribution data input requirements, the methods
available to describe this data, and the dialog boxes used to enter this data
are discussed in detail in Chapter “EPANET Input Description” on page 11.
As was discussed in the previous section, MIKE URBAN allows you to enter
data interactively into the pipe network using interactive dialog boxes. This
allows you to quickly define a model if a schematic layout of the pipe network
is unavailable, or if developing a schematic layout is too expensive or time-
consuming to construct. However, a powerful capability of MIKE URBAN is to
allow the user to graphically construct the pipe network model by simply
pointing and clicking with a mouse. This graphical construction of a pipe net-
work system is shown in Figure 2.3.
117
Using the Model
Figure 2.3 MIKE URBAN allows the user to interactively develop a pipe network
model by simply pointing and clicking to graphically add the compo-
nents to construct the water distribution network
Using the tools in the Component toolbar, for example, the user can select
the junction node component tool (Add Junction icon) and then interactively
locate and place junction nodes on the Map window. A background aerial
image can be displayed in the Map window to aid in the placement of network
components. While the user is moving the cursor in the Map window, a pre-
cise X,Y location is displayed on the status line at the bottom of the program.
To delete a network component, the select component tool can be used inter-
actively to select network elements. Then, selecting «Delete Selected» will
delete the selected elements. Note that you can undelete the deleted ele-
ments.
MIKE URBAN allows you to import network data, such as ESRI Shapefiles,
DWG, DGN, DXF files and raster images (such as scanned-in quad maps or
aerial photographs), and display them as a background layer in the Map win-
dow. This enables you to import a graphical background layer to facilitate lay-
ing out the pipe network. You can also import a DXF file containing lines and
polylines representing the pipe network and MIKE URBAN will convert the
lines and polylines to equivalent pipes.
To import a graphical data file, select Add Data tool. MIKE URBAN will then
display the Add Data dialog box, as shown in Figure 2.4. Select the file format
to import and then choose «Add».
Figure 2.4 The Import dialog box provides support for importing vector files and
raster image files
For large network systems, assigning demand data can be very tedious job.
Since many times the total demand is known for a particular network pres-
sure zone or for the entire network system, MIKE URBAN Water Distribution
provides the capability to distribute this total demand among the applicable
junction nodes.
The Distributed Demands dialog box, as shown in Figure 2.5, is used to auto-
matically assign the demands at the appropriate junction nodes. The Distrib-
uted Demands dialog box is reached by selecting EPANET | Distributed
Demands.
119
Using the Model
Figure 2.5 The Distributed Demands dialog box computes the demand at each
junction node automatically based upon the total demand within a pres-
sure zone or entire network system
Selecting the Method of Two Coefficients, MIKE URBAN computes the total
water demand assigned to each pipe (which is then split between the starting
and ending nodes) as:
Q k li k 2i
q pi = ---------------------
- (2.1)
k li k 2i
Selecting the Method of Reduced Pipe Lengths, MIKE URBAN computes the
total water demand assigned to each pipe (which is then split between the
starting and ending nodes) as:
Q l i k li
q pi = -----------------
- (2.2)
k li l i
Q l i k Di
q pi = -------------------
- (2.3)
k Di l i
where:
qpi =Total water demand applied to the pipe, split between the two end nodes.
li = Pipe length
This option will only assign demand to nodes with Demand Coefficient
applied (different from 0 or NULL). In the case of an equal distribution, the
node demand coefficients have to be equal and different from zero.
121
Using the Model
Q
q ni = ---- (2.4)
N
where:
QT
-C
q i = ----------------- (2.5)
i
i C
i = 0,n
where:
qi = node demand
Ci = demand coefficient
In cases, when node demands are retrieved from the Customer Information
Systems, it is possible to calculate pipe demand coefficients in the form of
aggregated demands for streets or the counted lots
Counted Lots
The pipe demand coefficient is derived from the X,Y position of the counted
lots. In order to use X,Y position of a point-attributes for the geocoding pro-
cess, it is necessary to create a database file, such as Dbase.DBF, Microsoft
Access.MDB where X,Y and other attributes are stored, or create the X,Y
position for each point in ArcGIS, for example. Use the import features of
MIKE URBAN to import such data to the Demand Allocation editor.
Figure 2.6 The demand coefficients (represented by circles) are assigned to the
nearest pipe within the snapping radius defined in the demand alloca-
tion dialog box
123
Using the Model
Demand Allocation
Junction node demands can be edited either within the junction editor for
each particular junction node or in the Multiple Demand editor, which allows
the user to display and edit all multiple demands. In addition to that, MIKE
URBAN provides another way of developing junction node demands, which is
based on geocoding the consumption data to the appropriate junctions or
pipes and aggregating their set point demand values to the junction
demands. This allows simplifies the demand development process and
allows you, for example, to import consumption data from the consumption
database systems and geocoded it based on its X, Y geographical coordi-
nates.
See the Tutorials Guide for the tutorial on demand allocation for a complete
description on how to perform a demand allocation with MIKE URBAN Water
Distribution.
can be linked to its diurnal curve. For more details, see Junction Editor and
Multiple Demand Editor.
Global Editing
Using the Field Calculator tool can change node demands. This allows us to
select the part of the network and to increase the node demand by 20%, for
example.
To do this, open the Multiple Demand Editor, select the column Demand and
select the Field Calculator from the right-mouse click pop-up menu. Define
the following statement:
DEMAND=DEMAND * 1.2
To increase node demands associated with the specific diurnal curve (Pat-
tern), use the following statement:
First, you need to check the model for errors. MIKE URBAN contains a pro-
ject check tool that can verify the data used to define an EPANET model (see
more in the Model Manager User Guide).
Second, after the program has checked the model for errors, you can then
execute the EPANET network hydraulic analysis. The EPANET network anal-
ysis program will read the input data and will then compute the flow rates,
pressures, and water quality for the defined pipe network.
Once the project check tool has successfully checked the model for potential
problems, you are ready to perform an analysis of the water distribution net-
work model. To execute the analysis, select Simulation | Run Simulation,
as shown in Figure 2.8.
125
Using the Model
Figure 2.8 The File Menu is used to select PERFORM ANALYSIS to analyze the
water distribution network
Please note that EPANET engine does not handle special characters such as
< > ? " : | \ / * . Make sure that neither the database nor the scenario ID con-
tain such characters.
After completing the analysis, the program will report the total number of
warnings generated during the analysis and whether the analysis ended with
an error. You can then view the analysis results. If the analysis ended with an
error or if the warnings are serious, you can then make the appropriate cor-
rections to the input data.
Figure 2.9 The LF1_N and LF2_N scenarios are selected for the batch analysis
The water distribution analysis results are written to an output file with the
same filename (but with a RES file extension) as the input data file.
127
Using the Model
In addition to the standard EPANET analysis results, MIKE URBAN can dis-
play the output results in a summary table, generate custom output reports,
and display many different graphical representations of the analysis results.
The following sections describe these capabilities.
MIKE URBAN can display multiple analysis result files for a network. In order
to compare two analysis result files for a network (for example, to compare
design alternatives of the same network), MIKE URBAN can subtract the two
analysis result files from each other. MIKE URBAN will then display the differ-
ence between these two analysis results for both node and pipe values (such
as node pressure, hydraulic gradeline, water quality, demand, flow, velocity,
and hydraulic headloss). The result difference can be displayed in the same
format as the standard analysis results.
To subtract two analysis result files, select Model Results| Results Compari-
son. Note that it is only possible to subtract two analysis result files if the
number of nodes, pipes, and time-steps (if performing an extended period
simulation) are the same. The following comparisons can be performed.
Different Demands
Network demands can be different (such as comparing a work and weekend
day), or local nodal demands can increase (such as simulating fire flows or
higher demands from what was predicted from the original master plan).
When subtracting the two analysis result files, a positive resulting value
means that the first analysis result file has a larger value than the second
analysis result file.
When EPANET performs its analysis of the water distribution network, it gen-
erates two ASCII output files containing the analysis results—a summary out-
put file and a complete output file. The summary output file contains a short
description of the analysis results. The complete output file contains the con-
tents of the summary output file in addition to output results for every compo-
nent within the water distribution network.
129
Using the Model
To view the EPANET summary analysis results, open the.SUM file using
Notepad, for example or select Summary from within Run Simulation dialog.
The Result Browser window allows you to graphically select any network
component in the Map window, by simply clicking on it with the mouse, and
the program will then display that component’s input attributes and analysis
results. This allows you to quickly examine the pipe network system at the
component level (i.e., pipe, junction node, valve, pump, tank, and reservoir),
check what is defined for the model, and determine the computed analysis
results. For example, selecting a pipe from the Map window will display in the
Browser window the pipe’s ID, diameter, length, roughness coefficient, minor
loss coefficient, reaction rate coefficient, open/close status, and flowrate.
Please also refer to the Model Manager User Guide for further on presenting
simulation results.
The Results Table allows you to display the tabular list of analysis results for
junctions, tanks, pipes, pumps, and valves. It is possible to display the results
for all elements or for the selection set. The selection se can be defined
based on the rules (constraints) such as ”pressure below 15”, “velocity above
1.5”, and similar. It is also possible to display the tabular results for nodes and
links, which are graphically selected from the Map.
The Map window allows you to graphically plot the analysis results directly
onto the pipe network schematic. To add a layer to the Map window, select
Model Results | Add Results Layer and pick up the selected result item, such
as Link:Flow, from the list.
The user can change most of the display options. In the Map window, com-
plete thematic mapping of the analysis results is available, including node
elevation, HGL, pressure, demand, and any water quality constituent. This
allows you to quickly interpret the modeling results and identify any trouble
areas. And, directional flow arrows can be plotted on top of the pipes to show
the flow direction for any time-step. In addition, MIKE URBAN provides auto-
matic color-coding of pipes and nodes based upon any input or output prop-
erty, allowing the network to be color-coded based upon pipe sizes, flowrates,
velocities, headlosses, nodal pressures, nodal demands, hydraulic grades,
elevations, water age, percent source contributions, water quality concentra-
tions, and any other attribute. Numerical ranges for colors can be specified.
Furthermore, pipes can be plotted with variable width and nodes with variable
radius, allowing you to quickly identify those areas of the network experienc-
ing the most flow, headloss, water quality constituent concentration, etc.
The link arrows corresponding to the link orientation, pipe slope and/or to the
flow direction can be drawn by the selected size. The node and link labels
can be displayed for all elements or only for selected elements. It is possible
in this way to display the labels just for filtered nodes and pipes, for nodes
and pipes along a specific profile or a flow path, for the individually selected
nodes etc. This is useful especially when we want to display very quickly
results for several elements on the plan. The values of pipe diameter and
pipe length can be displayed together with the prefix.
Please also refer to the Model Manager User Guide for further on presenting
simulation results.
Profile plots allow you to graphically plot the analysis results along any pipe-
line path. To display a profile plot, a profile path must first be defined from the
pipe network horizontal plan. Once the profile path has been defined and the
profile plot displayed, the path can be saved for later re-use.
Profile plots allow you to graphically plot the analysis results along any pipe-
line path. Profile plots can have two separate vertical axes to allow plotting of
variables from two separate unit families, such as flow and pressure. Profile
plots can be plotted along any user-specified route. Profile plots can be gen-
erated as line graphs, bar graphs, or mixed—along with complete graph cus-
tomization. For example, profile plots can be plotted with an envelope to
show the minimum and maximum values reached during an extended period
simulation.
Please also refer to the Model Manager User Guide for further on presenting
simulation results.
Time series plots allow you to graphically display the analysis results for any
network element for an extended period simulation. Multiple time series plots
can be generated for the various network elements, such as pipe flow, veloc-
ity, headloss, nodal demand, pressure, hydraulic grade, water age, water
quality constituent concentration, pump characteristic operating curve, tank
water level, total and net system demand, etc.
Please also refer to the Model Manager User Guide for further on presenting
simulation results.
131
Using the Model
Time series scatter graphs can be used to display the Q-H relation for any
selected link for an extended period simulation.
The Q-H relation plots allow you to graphically display the fire flow analysis
results (calculate hydrant curve) for any fire hydrant node.
The position of the current operating points can be displayed for any pump for
an extended period simulation. In order to display this graph, load the simula-
tion results, open the Pump Editor and display the graph.
It is possible to detect pipes where the flow direction has changed during the
hydraulic simulation. The same tool can also be used for detecting links and
nodes, where the results (pressures, for example) changed more than a
specified limit.
2.6.12 Animation
See the Model Manager book for more information on viewing simulation
results.
The EPANET Errors and Warnings dialog box allows you to read the error
and warning messages from the .SUM file, generated by EPANET and store
this information in the database manner. All messages are classified based
on their type, which makes it possible to filter out errors or warning related to
pumps or to filter out errors and warning for a specific time level or an model
element, for example.
Figure 2.11 The EPANET Errors and Warnings dialog box is used to process simu-
lation messages
A list of the EPANET Errors and Warnings dialog box data fields for Figure
33.29follows, with a short description given for each entry.
MESSAGE TYPE
This list box classifies the error and warning message, such as FCV, Pump,
Pipe, Negative Pressures, System Unbalanced and similar.
LINK/NODE ID
This data entry contains the model element ID, where applicable.
LINK/NODE TYPE
This data entry contains the model element type, where applicable, such as
Junction, Tank, Pipe, Pump, or Valve.
MESSAGE
This is the error or warning message generated by EPANET engine.
SIMULATION TIME
Simulation time of the error or warning message is available in
[days:hrs:min:sec] format as well as the number of elapsed seconds from the
simulation start.
ADVANCED...
It is possible to Select or Select All messages and map highlight the corre-
sponding network elements on the Map; it is also possible to open Select by
Attribute dialog box and filter out messages based on their type, simulation
time, model element identification or a type. HTM based report can also be
generated.
133
Using the Model
It is possible to store the complete model results (for all time steps) and for
the selected part of the network into the database. This is useful when you
want to use Microsoft Access or Excel tools for analysing the model results
and for using user-defined tools for the additional results processing. In order
to store the model results to the database, load the simulation results into
MIKE URBAN, and select nodes and links you want to use for storing the sim-
ulation results to the database. You can also store all nodes and links in case
you do not use the selection set.
A fire flow is the maximum flow rate available at a specific minimum pressure,
typically 20 psi (15m). There are two basic ways to model a fire flow:
1. Specify a design fire flow rate and compute the available fire flow pres-
sure.
2. Specify a design fire flow pressure and compute the available fire flow
rate.
135
Water Distribution Tools (WD-Tools)
Figure 3.1 The Fire Flow dialog box is used to specify fire flow analysis parameters
A list of the Fire Flow Analysis dialog box data entries for Figure 3.1 follows,
with a short description given for each entry.
137
Water Distribution Tools (WD-Tools)
Maximum velocity
This data entry allows you to define the maximum velocity in pipes
Select Model Results | Results Browser to display the Results Browser win-
dow and click the junction node to view the simulated fire flow results.
The simulated fire flow results for “Calculate available flow for design pres-
sure“ and “Calculate available pressure for design flow“ are reported to the
*.FFA text file. Open this text file using any file processing program such as
Windows Notepad to see the detailed list of the fire flow results for each sim-
ulated junction node. The FFA text file format is suitable for importing into
Microsoft Excel, for example. The simulated fire flow results for “Calculate fire
hydrant Q-H curve“ are reported to the *.FQH text file.
The simulated fire flow results can be displayed using the colour coded maps.
Select Add Results Layer and select one of the fire flow result items to create
a colour coded map with the fire flow results.
Both English and Metric units are fully supported in the Pipe Roughness Cali-
bration. Pipe roughness coefficients can be calculated for the Hazen-Wil-
liams, Darcy-Weisbach (Colebrook-White), or Manning friction loss
equations.
139
Water Distribution Tools (WD-Tools)
Although different and intended for different purposes, all EAs share a com-
mon conceptual base (schematized in Figure 3.2). In principle, an initial pop-
ulation of individuals is created in a computer and allowed to evolve using the
principles of inheritance (so that offspring resemble parents), variability (the
process of offspring creation is not perfect -- some mutations occur) and
selection (more fit individuals are allowed to reproduce more often and less fit
less often so that their "genealogical" trees disappear in time). One of the
main advantages of EAs is their domain independence. EAs can evolve
almost anything, given an appropriate representation of evolving structures.
Similarly to processes observed in nature, one should distinguish between an
evolving entity's genotype and its phenotype. The genotype is basically a
code to be executed (such as a code in a DNA strand), whereas the pheno-
type represents a result of the execution of this code (such as any living
being).
3.2.2 Introduction
141
Water Distribution Tools (WD-Tools)
Pipe Groupings
Pipes can be lumped together in separate logical groups based on their
known physical characteristics such as pipe material, age, and diameter. It is
assumed that all pipes within a calibration group (where a single pipe may
constitute a group) will possess an identical roughness coefficient. Any com-
bination of pipe calibration groups can be specified and their fitness evalu-
ated to match field observations. Selected pipes can be excluded from these
groups and their roughness coefficients will remain unchanged during the
calibration process. The Group ID can be assigned for each pipe the Pipe
Editor. Pipes with undefined Group ID will be excluded from the automatic
pipe roughness calibration.
Group Definition
We can create as many pipe groups as necessary and define their minimum
and maximum roughness coefficients. Minimum and maximum pipe rough-
ness coefficient is defined for each pipe. Any group can be also displayed on
the Map.
Figure 3.6 Pipe roughness groups define the roughness coefficient limits
The measured flow and pressure values can be entered from within the cor-
responding editors.
143
Water Distribution Tools (WD-Tools)
The initial population size is the number of initial generation members (the
number of networks), which is used by the automated calibration.
After step 4 has been completed, the cycle is repeated. Three thousand
hydraulic analyses will be executed in this example.
Figure 3.9 Summary file with then calibrated pipe roughness coefficients
Coefficients to Pipes; note that this will assign the accepted pipe roughness
to each pipe based on the calibration pipe GroupID.
The control functions range from the simplest constants for the operational
variables (e.g. constant water level) to dynamic controlled variables, set in a
continuous functional relation with any of the measurable variables in the sys-
tem (e.g. pump START/STOP levels as functions of water level at strategic
location in the system).
145
Water Distribution Tools (WD-Tools)
In the latter case, the control decisions are derived by evaluating (comparing)
the current value of the controlled flow variable and the pre-defined set-point
value. The control algorithm is based on the numerical solution of the "contin-
uous control problem" equation and is usually termed as PID (Proportional-
Integral-Differential) control. The actuation signal for the regulator is gener-
ated by a PID controller, which usually appears as part of the operational
strategy programmed in a Programmable Logical Controller (PLC).
Selection of a controlled variable is, however, subject to limits set by the vari-
able's "controllability". Therefore, a controlled variable is usually selected
among the flow variables (flow, water level), preferably in the vicinity of the
regulator. As a controlled variable becomes more distant from the regulator, it
becomes more difficult to control, due to time lags, diffusion and uncontrolla-
ble interference. Control of relatively distant controlled variables is difficult
and often cannot give satisfactory results.
The global control can be extended also to include forecast data in addition to
real-time data, which is then called predictive RTC. The most comprehensive
way to obtain forecast data is to include a model in the control system. Pre-
dictive control brings additional benefits in relatively inert systems, i.e. where
the response time of an operational variable is long compared to the change
of relevant disturbance (external input or control action).
MIKE URBAN RTC simulates reactive local and global RTC systems in water
supply and water distribution networks. The software implementation is inher-
ently a conceptualisation of real life, of which the user must be fully aware.
Some conceptualisations applied in MIKE URBAN RTC are listed below.
The program does not distinguish explicitly between local and global
RTC. Per default, all elements of a modelled RTC system are assumed
available for global control.
Sensors are specified as operational devices with definition of sensor
type and position in the MIKE URBAN network. Sensors with multiple
functionality must be specified individually.
When devices (valves and pumps) are specified as controllable in the
MIKE URBAN interface, a number of additional physical parameters
about the behaviour of the structure is required to describe e.g. the
allowed change rates for the state of the structure.
The actual controllers are not specified explicitly as physical devices, but
their function (i.e. operational logic as a combination of operational con-
ditions and control functions) is associated with the respective devices.
MIKE URBAN RTC uses sampling and actuation (control loop) frequency
identical to the simulation time step.
Sensor readings are simulated as perfectly accurate and with 100%
availability.
Low-level logic of the pump START/STOP operation is built into the pro-
gram and is controlled by the START and STOP levels.
The PID control algorithm is built into the program and is controlled by the
PID constants and by factors for weighting the terms of the numerical solution
of the control equation.
147
Water Distribution Tools (WD-Tools)
Figure 3.10 through Figure 3.12 show examples of how the actual variable
(flow or water level) can fluctuate around the set-point as consequence of
various choices of the PID constants. Each figure has three different graphs
depending on whether the constant is too high, too low or adequate.
Figure 3.10 Fluctuations around the set point depending on the size of the propor-
tionality factor, K
Figure 3.11 Fluctuations around the set point depending on the size of the deriva-
tion time, TD
Figure 3.12 Fluctuations around the set point depending on the size of the integra-
tion time, Ti
At each simulation time step, the set-point is evaluated against the actual
value of the control variable (flow or water level, depending on the set-point
type).
The actual value of the control variable is determined from the following
equation:
1 i
T
de
u = K e + e dt + T d
Ti 0 dt
(3.1)
where:
149
Water Distribution Tools (WD-Tools)
1
U(s) = K 1 + + T d s E(s)
Ti s (3.2)
where:
z -1
s=
Ts z (3.3)
where Ts is the sampling period, one obtains the following discrete rep-
resentation of the PID equation:
u(n) = K 1 + T s + T d y ref (n) - y(n) -
Ti Ts
2
- K 1 + T d y ref (n - 1) - y(n - 1) +
Ts
+ K y ref (n - 2) - y(n - 2) T d +
Ts
u(n - 1) (3.4)
where:
a1..a3=factors for weighing each time step (index n, n-1 and n-2)
y=actual value of the regulated variable
Indexes n, n-1 and n-2 denote the current, the previous and the second previ-
ous time step, respectively.
The weight factors a1... a3 have been added to the PID algorithm in order to
give the user more ways of stabilising the algorithm in case of instability prob-
lems.
The following is the list of parameters that are used to specify the PID based
control.
Control Settings
Controlled element settings define the element, such as a pump or a valve
which is controlled by the PID algorithm.
151
Water Distribution Tools (WD-Tools)
Maximum decrease rate: the maximum rate at which the variable can
decrease. In units of seconds-1.
Set point type: type of the set point node; such as junction, tank.
Set point variable
Set point ID: ID of the set point node.
Set point value: the constant value.
Set point value curve: pattern ID defining how the set point value
changes in time.
PID Settings
Maximum value: the proportionality constant of the PID algorithm (KP
proportionality constant)
Max. increase rate: the integration time used in the PID algorithm (TI
integration time)
Max. decrease rate: the derivation time used in the PID algorithm (TD
derivation time)
The speed of the valve or pump control depends on the PID control settings.
The PID controls and the Control Element Settings must reflect the physical
limitations or restrictions of the real-life controlled element and the PID algo-
rithm.
The Wagner equation is used to adjust the node demand based on the avail-
able pressure.
where:
Note that nodes with negative demand i.e. inflow nodes are excluded from
this demand adjustment. In order to use Demand Adjusted Analysis, select
EPANET | WD-Tools | Demand Adjusted Analysis and define the simulation
settings.
Figure 3.14 The simulation settings for the Demand Adjusted Analysis
A list of the Demand Adjusted Analysis entries with a short description given
for each entry:
153
Water Distribution Tools (WD-Tools)
Figure 3.15 Run the Simulation with the Demand Adjusted Analysis choice
The simulation will run for any steady state or extended period model and the
simulation results will be presented in the same way as standard results. In
case that the node demands were adjusted, simulation report.INF file will pro-
vide more details. The INF-file is found in the project directory.
Where:
Figure Figure 3.17 illustrates how the original network demand (as defined in
the Multiple Demand Editor) was adjusted based on the actual pressures; this
graph was created based on the .INF file.
Figure 3.17 Comparison of the original network demand versus adjusted demands
during 24 hours
The Demand Adjusted Analysis can be used to model systems with intermit-
tent water supply or systems with low pressure conditions. We can also use
this concept in modelling system breakdowns, such as closure of parts of the
water distribution systems.
155
Water Distribution Tools (WD-Tools)
Figure 3.18 Low pressure scenario, total network demand ~ 300 l/s (given demand)
In order to see what is the network capacity, we can run the Demand
Adjusted Analysis with the settings of minimum pressure = 5 m, for example.
We can see that the network capacity is 83 l/s in order to provide minimum
pressure of 5 m, see Figure 3.19.
Figure 3.19 Pressure dependent demands, total network demand ~ 83 l/s (given
demand)
See the section titled Pump Editor on -18 for additional information.
The Pressure Zone Mapping wizard allows you to define the separators,
merge small groups of pipes into the same group, and it can automatically
generate the colour legend to display the pressure zones.
157
Water Distribution Tools (WD-Tools)
Merge zones smaller than: in case that there are many small zones (a
typical example would be small pipes located in pumping stations and
storage tanks), they will be all merged into the same zone for graphical
display instead of creating a separate zone for each of them.
This tool help you to understand the hydraulic behaviour of the network prior
to running the hydraulic simulation and it also helps you to detect possible
errors in the network connectivity.
Figure 3.21 Pressure zones are automatically mapped based on the connectivity
breakdown rules.
3.8 Advisor
Advisor is a set of tools designed to support the analysis of the water supply
and water distribution systems. It helps the user to check the results of the
hydraulic simulations, understand system deficiencies, improve design and
operation. Advisor is available under EPANET|WD-Tools|Advisor and con-
tains the following tools:
Pipe criticality
Sustainability analysis
Cost analysis
Water demand criteria (PI-1) is computed as water (in flow units) that cannot
be delivered through the pipe. The value of 1 corresponds to the total flow.
Pipe flow criteria (PI-3) is computed as the value of average hour demand
Q(hour), the value of 1 corresponds to the total value of average hour
demand.
159
Water Distribution Tools (WD-Tools)
It is recommended to use the maximum peak demand for the above simula-
tions.
In order to use Pipe Criticality, select WD-Tools from the main program menu
and then select Advisor | Pipe Criticality. The Pipe Criticality form where you
define the simulation settings will then display.
Time level: in case of Extended Period Simulation, this is the time level at
which the pipe criticality analysis will be performed
Minimum service pressure: minimum acceptable service pressure within
the network
The pipe reliability settings can be defined on the level of the whole model
area or per each specific zone. In case of Global settings:
In case of the local settings, the network demands are defined per each zone,
please note that the zone demands are defined in the Pressure Zone Editor.
Q: network demand
Qhr: maximum hour network demand
Selected: use this switch to active or deactivate zones; only selected
zones will be used for the pipe criticality analysis.
In order to run the pipe criticality analysis or to load the results of the previous
pipe criticality results, use the following controls
161
Water Distribution Tools (WD-Tools)
Figure 3.25 PIpe criticality results, Map view (pipes are colour coded by the pipe
criticality C)
In order to use Sustainability Analysis, select WD-Tools from the main pro-
gram menu and then select Advisor | Sustainability Analysis. Sustainability
Analysis dialog where you define the simulation settings will display.
Please note, that the simulation results needs to be loaded in order to start
the sustainability analysis.
Flow threshold: minimum flow criteria used for reverse flow calculation
The program will provide a detailed analysis of the simulation results and it
will create the following layers and reports:
Service pressures
– Pressure anomalies e.g. pressure too low or to high
163
Water Distribution Tools (WD-Tools)
In addition to the Map layers, the HTML report creates a detailed report for
the whole network as well as for the particular zone.
Storage tanks: reports tanks that are either drained or overflows during
the simulation, reports if the tanks are balanced within the simulation.
Balanced tank is a tank where the water level at the beginning of the sim-
ulation is the same as at the end of the simulation.
Pumps: reports pump that are operated near their maximum capacity
Unit Headloss: reports pipes with too high unit head loss
Service Pressures: reports excessive pressures
Pipe flows: reports reverse flows
Flow velocity: reports excessive flow velocity
Cost analysis allows you to review the energy consumption results on more
details, create tabular outputs, and great graphs of pump utilization, average
power consumption, and costs.
In order to use the Cost Analysis, load the simulation results and select
EPANET | WD Tools | Advisor | Cost Analysis…
165
Water Distribution Tools (WD-Tools)
Utilization: Percent utilization i.e. percent of the time that the pump was
operating (%)
Efficiency: Average efficiency of the pump (%)
kW-hr/m3 or kW-hr/gallon: Average power consumption (kW per hour)
per million gallons (or cubic meters) pumped
Average kW: Average rate of energy usage of the pump (kW)
Peak kW: Peak rate of energy usage of the pump operation (kW)
Cost/day: total cost of the pump operation per day (monetary units)
Note that the total cost and demand charge (if any) are listed at the bottom of
the table.
To define energy price, price pattern, pump efficiency use EPANET | Energy
Editor and EPANET | Pump Editor where you can define the following fields:
Alarms and Violations allows you to define critical values for various results
items anywhere within the model network such as maximum velocity, mini-
mum or maximum pressure, low or high level, high water age and let the
hydraulic model to evaluate them based on the actual simulation results.
The dialog table allows you to define or display the following parameters:
Element MUID: enter the MUID of the junction node, pipe, pump, valve,
or a tank
167
Water Distribution Tools (WD-Tools)
– Link: Flow
– Link: Absolute Flow
– Node: Water Quality
– Link: Velocity
– Link: Unit Headloss
– Link: Water Quality
– Link: Status Code
– Link: Setting
– Link: Reaction Rate
– Link: Friction Factor
– Link: Start End Head
Actual Value: this is the highest or lowest of the actual values resulting in
the alarm
Note, that the alarm will be triggered if the criteria is fulfilled i.e. if the criteria
is defined for a “tank level < 2” then the alarm will be reported once the com-
puted level is below “2”.
Similarly, if you want to alarm for a pump operation during which the pump
flow exceeds 300 l/s, enter the following criteria “pump flow > 300”.
Validate: perform the results validation for the selected results file
Clear: clear the results of the last validation
On Top: forces the dialog to remain on top of other windows
Close: close the dialog window
169
Water Distribution Tools (WD-Tools)
Q g A H
0
x a2 t
(4.1)
Q H f
gA Q Q 0
t x 2 DA (4.2)
In the governing equations the acceleration terms which are very small com-
pared to the other terms have been disregarded.
The general expression for the wave speed (only important for water hammer
computations) presented by Halliwell (1963) has been used
K
a
ρ 1 K/E ψ (4.3)
171
MIKE URBAN Water Hammer
Q g A H
0
x a2 t (4.4)
Q H f
gA Q Q 0
t x 2 DA (4.5)
In the governing equations the acceleration terms which are very small com-
pared to the other terms have been disregarded.
Wave Speed
For pure liquids Halliwell (1963) presented the general expression for the
wave speed
K
a
1 K / E (4.6)
Rigid Conduit
= 0 (4.7)
Ro 2 Ri 2 2 Ri
2
2 1 2
Ro Ri Ro 2 Ri 2
2
(4.8)
Ro 2 15
. Ri 2 Ro 2 3 Ri 2
2 2
Ro Ri Ro 2 Ri 2
2
(4.9)
Ro 2 Ri 2
2 2
Ro Ri
2
(4.10)
D
e
1 2
(4.11)
173
MIKE URBAN Water Hammer
D
e
1 0.5
(4.12)
D
e (4.13)
= 1 E G (4.14)
DE
GD Ee (4.15)
in which e is the thickness of the steel liner and E is the modulus of elasticity
of steel.
As
ee E r ec
Ls (4.16)
Diagrams
The following diagrams can be used in order to estimate the wave speed.
Table 4.1 Values of Young's Modulus of Elasticity and Poisson's Ratio for a range
of common materials
Material Young’s Modulus Poisson’s Ration (-)
(10E9 N/m2)
Aluminum 70 0.3
Glass 68 0.24
GRP 50 0.35
Polyethylene 3.1 -
175
MIKE URBAN Water Hammer
f
2 DA
Q Q
1 f
2 2 DA
Q nj 1 Q nj 11 Q nj Q nj 1 (4.17)
The coefficients for the water hammer model have been derived and have the
following form:
gA
2
a (4.18)
fric
2AD (4.19)
a t
Cr
x (4.20)
3a t
2
2 2 2 x
c 6 6 1 Cr 6 6 1 2 2
(4.21)
3a t
2
2 2 2 gA x
m 6 6 1 Cr 6 6 1 2 2
(4.22)
A1 1 x (4.23)
B1 t c (4.24)
C1 x (4.25)
D1 t c (4.26)
A2 g a t m (4.28)
177
MIKE URBAN Water Hammer
B 2 1 x fric t x Q nj1
(4.29)
C 2 g a t m (4.30)
D 2 x frict x Q nj
(4.31)
E2 1 xQ nj11 xQ nj1 gA1 t H nj1 H nj H nj11 H nj1
gA t H nj1 H nj m H nj1 2H nj H nj11 2H nj1 (4.32)
The main algorithm generates a set of grid points using a finite difference
scheme, see Cunge, Holly, Verwey (1980). The grid is introduced in time and
space, where at every point the values of H and Q are defined as the
unknown variables. Between the two successive grid points in time and
space both the continuity and the momentum equation are applied. Together
with the necessary boundary data, a sufficient number of equations are
obtained to solve H and Q at every grid point.
A2 j H nj 11 B2 j Q nj 11 C 2 j H nj 1 D2 j Q nj 1 E 2 j
(4.34)
The looped algorithm is based on the fact that a looped network contains ele-
ments known as nodes which represent the confluence of several flow paths,
some of which originate from other nodes, some from boundary points. A sys-
tem of simultaneous linear equations is developed where the piezometric
head changes at each node are the only unknowns. Solution of this system
by any matrix elimination technique yields the piezometric heads at each
node.
Suppose that there are three links, 2-1,2-3 and 2-4 and that there are b grid
points along branch 2-3 and c grid points along a link 2-4, see Figure 4.2. For
any computational grid point, equations (4.35), (4.36) may be written as
Hi L1i H1 M 1i H jj N 1i
(4.35)
Qi L2i H1 M 2i H jj N 2i
(4.36)
m
Q
k 1
n 1
1k 0
(4.37)
where n+1 indicates the (n+1)Dt time level in the solution, k is the index of the
links emanating from node 2, and m is the number of such links. These rela-
tions can be written for each from M nodes, and this leads to a system of M
179
MIKE URBAN Water Hammer
This system of linear equations may be solved by any matrix inversion tech-
niques. Once the increments of piezometric head H are known at the nodes,
it is possible to recompute Q(i) and H(i) values for all intermediate grid points
through equations (4.37) and (4.38).
The main algorithm is designed in such a way, that, after a process of lineari-
sation and discretization of the governing equations, it solves them on a pre-
scribed set of grid points using an appropriate numerical scheme. If a
hydraulic structure is present in the domain of the solution, the algorithm must
replace the governing equations by other equations defining a hydraulic
structure in order to provide the numerical solution. Various hydraulic struc-
tures can be coupled together, e.g., the closing of one valve can determine
the operating of another valve. In cases where this link exists between
hydraulic structures, communication must be maintained and controlled by
this main algorithm. This message has to be attached to the object in such a
way that it represents the reality. Object-oriented design has been applied to
create a safer interface to the numerical algorithm, since the low level opera-
tions that remain the same are hidden inside objects.
The present version will handle any number of pipes, nodes, and loops in
complex networks with various components
181
MIKE URBAN Water Hammer
In order to be able to start water hammer simulations you have to prepare the
steady state model and obtain satisfactory results. In the next step, you need
to specify Water Hammer Analysis type in the Project Options dialog, define
transient boundary conditions and computational parameters. In addition to
the specific pipe data, you need to know the wave speed in each pipe.
Initial Conditions
Initial conditions are computed with the use of the Steady State Model. The
results of the initial state are saved in the file as H, Q values at the beginning
and end of the pipes respectively and in the vicinity of hydraulic structures
such as valves, pumps, etc. There is a direct connection between the result
file from initial conditions and the water hammer execution, in spite of the fact
that the two models use different computational grids.
Boundary Conditions
There are in principle two types of boundary conditions, namely the piezo-
metric head, H, above a specified datum, e.g., in tanks, and the discharge, Q,
e.g., water demand. Both H and Q are given under selected names as time
series in the Curve Editor and stored in the database. These boundary condi-
tions may be assigned to any node in the network. Boundary for each time
step is assigned from given time series specified by the user. If time step
used by water hammer computation is smaller than appropriate neighbouring
values in boundary conditions time series then linear interpolation is applied.
There are nodes of the following types: H - boundary, Q - boundary, compati-
bility and structure (hydraulic component) description. It should be pointed
out that time patterns, used in the Steady State Model, are ignored by the
Water Hammer simulations.
For the Initial State for Water Hammer Model, the water level and/or dis-
charges are constant in time. The boundary conditions using time series must
be specified for a sufficiently long time interval.
Computational Parameters
General parameters consist of fluid density, fluid bulk modulus, absolute tem-
perature, vapour pressure and gravitational acceleration. The most important
numerical parameter is a time step. Since a numerical solution must be stable
and as accurate as possible, you have to choose a proper value of Dt. The
stability condition is given by the Courant number
a t
Cr
x (4.40)
in which a is the wave speed and Dx is the distance between two successive
grid points. In principle, an implicit, space-compact scheme is unconditionally
stable, with exact solutions generated for the Courant numbers Cr = 0.5 and
Cr = 1.0, respectively. The scheme enables us to vary the Courant number
over pipes while maintaining its high accuracy. Accurate results are produced
in the range 0< Cr < 1.1. You should try to maintain the Courant numbers
below unity, but as close as possible to Cr = 1. If you select the menu item
Geometry and Branch, you can control the values of Courant numbers. The
question how to choose the time step is dictated by the nature of the hydrau-
lic transient itself and by the shortest pipes in the system. The time step can
vary from the order of 10-3 to 10+1 seconds. The time steps must be small
enough in order to describe very fast changes of variables. It is recom-
mended to start with the shortest pipe section and to calculate the time step,
considering Cr = 1. Pipe sections with high Courant numbers are numerically
treated in MIKE URBAN Water Distribution as rigid pipelines. This simplifica-
tion enables a user to deal with a very short pipe section which would not be
important within the water hammer simulation. Once the time step has been
prescribed, you have to input the simulation time. MIKE URBAN Water Distri-
bution calculates a number of time levels, which you need to prescribe in the
Check level item. In the Project parameters' window you can also change
throughout network whether you intend to use a friction factor and/or an
absolute roughness.
183
MIKE URBAN Water Hammer
Branches
can be used to represent pipes of constant properties. In the pipe network,
branches may include hydraulic elements, for example, valves, pumps.
Nodes represent the applicable boundary conditions at the end of branches.
Nodes
are elements that represent free branch ends, branch connections or a spe-
cific storage. At nodes with one simple pipe connected, boundary conditions
are usually defined by specifying the values of piezometric head or discharge
as a constant value or as a function of time. Flow continuity and a piezometric
level compatibility is assumed at nodes connecting several branches
together.
Generally, there are these three different types of nodal boundary conditions:
Grid Points
are generated automatically by Water Hammer along the branches and they
represent the computational grid where the values of piezometric head and
discharge are solved and the input and/or output data are required. The sys-
tem requires a different computational grid for steady state and water ham-
mer computations.
Control Rules
Simple Control Rules and Rule Based Controls are not allowed during the
water hammer analysis. Valve opening and pump scheduling is handled
directly by the specific valve and pump data.
Pumps may be located inside pipeline systems (booster pumps) or they may
be connected to a suction well. Pumps are frequently used for various pipe-
line systems, and may operate during hydraulic transients with constant
pump speed. Alternatively, the pump speed can decrease and/or increase
depending on pump shut-down and/or pump start-up. The greatest difficulties
185
MIKE URBAN Water Hammer
come from hydraulic transient flows caused by turbopumps, since they may
work in four quadrants. Four quadrants pumps are currently not supported.
discharge - Q (m3/s)
total dynamic head (tdh) - H (m)
rotational speed - N (rpm)
shaft torque - T (N.m)
V 2 p V 2
tdH H d d z s p s z
2g g d 2 g g s
(4.41)
where the subscripts, d and s denote the discharge and suction flanges,
respectively. Power input P (kW) is defined as:
gQ H T 2 N
P T
60
(4.42)
where h is the pump efficiency and T (N.m) is the torque which may be calcu-
lated from this equation.
H=b + aQ2
If the pump characteristics does not satisfy parabolic relation large errors may
be produced in GGM method and in all computation modules if the pump dis-
charge is out of the Q-H curve
2 3
Q N H N P N
1 1 1 1 1 1
Q N H N P N
2 2 2 2 2 2 (4.43)
Subscripts 1 and 2 are only for corresponding points on an affinity law parab-
ola. The affinity laws for discharge and head are accurate for all types of cen-
trifugal pumps. However, large errors may be produced using the affinity law
for a power requirement. It is recommended to compute P from head, dis-
charge and efficiency and not from affinity laws.
gQH T 2 N
P T
60
(4.44)
where Iw (N.m.s) is combined moment of inertia and Dt (s) is time step used
for the calculation.
187
MIKE URBAN Water Hammer
N
Tm T t30
I
(4.45)
The relation between the pump speed and the total pump dynamic head is
described by the following equation:
tdH100% 2
tdH t 2
Nt
N100%
(4.46)
where, index (100%) represents the 100% of the pump rated speed and the
time index t represents the actual value of tdH and N during the analysis.
n
I mrj
i 1
(4.47)
The moment of inertia of a body depends on the axis about which the body is
rotated. If two axes of rotation have different distributions of mass around
them, then the body will have different moments of inertia for each of these
axes.
Torque, a twisting effort applied to an object that tends to make the object turn
about its axis of rotation. The magnitude of a torque is equal to the magnitude
of the applied force multiplied by the distance between the object's axis of
rotation and the point where the force is applied. In many ways, torque is the
rotational analogue to force. Just as a force applied to an object tends to
change the linear rate of motion of the object, a torque applied to an object
tends to change the object's rate of rotational motion.
189
MIKE URBAN Water Hammer
Pump
A Pump can operate in 3 quadrant with constant and/or variable rotational
speed. Parameters which describe a Normal Pump are as follows:
[ID]
Pump ID
[N1]
Starting node ID
[N2]
Ending node ID
[rpm], [N]
[flow direction]
[operation type]
[operational schedule]
N-time curve
[Iw]
[pump torque]
T(pump)-Q curve
[motor torque]
T(motor)-N curve
Figure 4.7 Pump Q-H curve defined in 3 quadrants for 100% rated speed
The relationship between the flow Q and the head drop DH is expressed
using a discharge coefficient Cd for:
In-line valve
1
Q Cd Av 2 gH Cd
(4.48)
where Av is the valve area and x is the valve minor loss coefficient.
Free-discharge valve
1
Q Cd Av 2 gH Cd
1
(4.49)
191
MIKE URBAN Water Hammer
Remarks:
TCV Throttle Control Valves can also be used as Isolation Valves for example
for isolation of a pipe section in case of repair, isolation of a pump, etc.
Parameters describing Control Valves are as follows:
[ID]
Valve ID
[N1]
Starting node ID
[N2]
Ending node ID
[Dd]
Valve diameter
[Valve characteristic]
or,
[Valve schedule]
Kv
Cd
3600A 2 g
(4.50)
193
MIKE URBAN Water Hammer
The relation between Cd valve coefficient and x minor loss coefficient is given
by the following equation:
1
Cd
1
(4.51)
1
Cd
(4.52)
Note, that the valve minor loss coefficient used for the steady state analysis
must correspond the initial valve opening used for the water hammer analy-
sis.
t
Cr 1 a
tx (4.53)
195
MIKE URBAN Water Hammer
Component Remark
Tank Available
Pump Available
Component Remark
Component Remark
Tanks
Surge Tanks have been widely used for hydroelectric systems in order to pro-
tect the low-pressure supply tunnel. They may also sometimes be suitable for
water supply schemes. There are various types of Surge Tanks. The sche-
matic presentation of common Surge Tanks is the same as mentioned above
for Tanks.
The governing equations describing their hydraulic behavior are the dynamic
equation and the continuity equation. Losses are disregarded at the junction,
but are taken into account for pipes. Parameters characterizing the Surge
Tank are:
Parameters:
Node ID.
Maximum water depth above datum.
Starting water depth for computation.
Tank bottom level.
Tank Type: Rectangular tank: [a] [b] right prism rectangular tank, the
base with sides a, b; Circular tank: vertical cylinder with diameter D; Var-
iable: depth versus volume curve.
Air-Chambers
Air Chambers contain compressed air which prevent very low minimum pres-
sures in the pipeline and hence column separation. They are frequently used
197
MIKE URBAN Water Hammer
behind the pumps in water supply pipelines. Mostly they are cylindrical with a
vertical and/or horizontal axis. A horizontal cylinder may be preferred for a
very long pipeline when a large volume of air is required. The analysis is sim-
ilar for both cases, but the computation of the volume of air in a horizontal cyl-
inder is more difficult. Figure 4.11. illustrates an Air Chamber with a vertical
cylindrical tank.
C H *
air air (4.54)
in which Hair and "air are the absolute pressure head and the volume of the
enclosed air, k is the exponent in the polytropic gas equation (k = 1.0 for an
isothermal expansion, k = 1.4 for adiabatic expansion). The orifice losses are
different for the inflow and outflow from the chamber.
Parameters:
Node ID. Air chamber can be located only in the "dead-end" node, i.e.
node with only one adjacent pipe.
Maximum water depth above datum.
Starting water depth for computation.
Tank bottom level.
Tank Type: Rectangular tank: [a] [b] right prism rectangular tank, the
base with sides a, b; Circular tank: vertical cylinder with diameter D; Var-
iable: depth versus volume curve.
Kappa: k = the exponent in the polytropic gas equation (default value
k=1.2).
C H *
air air (4.55)
in which Hair and "air are the absolute pressure head and the volume of the
enclosed air, k is the exponent in the polytropic gas equation (k = 1.0 for an
isothermal expansion, k = 1.4 for adiabatic expansion). m is the entrapped air
volume, pair is the air pressure inside a chamber for which the Qair (air dis-
charge) is taken from the dual-acting valve characteristics. Next chart shows
characteristics of Pont&Mousson, Ventex dual-acting valve, diameter of
150mm.
199
MIKE URBAN Water Hammer
When the pressure inside the Vented Air Chamber drops below the atmos-
pheric pressure, dual-acting valve opens and the air flow into a chamber. The
proper valve characteristics are required to set by a user. As soon as the liq-
uid starts flowing back into the dual-acting valve, valve closes.
Parameters:
Air-Valve
Air valves, similar to Vented Air Chamber contain air which prevent very low
minimum pressures in the pipeline and hence column separation. Air valves
are modelled as small Vented Air Chamber equipped by dual-acting valves
that allow air to be sucked into its chamber and to escape therefrom, while
preventing the outflow of liquid. When the pressure inside the surrounding
pipes drops below the atmospheric pressure, air-valve opens and the air flow
into a system. The proper valve characteristics are required to set by a user.
As soon as the liquid starts flowing back into the dual-acting valve, valve
closes.
Parameters:
Node ID. Air chamber can be located only in the "dead-end" node, i.e.
node with only one adjacent pipe. Air-valves are defined as junction
nodes with an air-valve type.
Valve diameter.
Dual-acting Valve characteristics: Relation between gauge pressure and
the air flow through the dual-acting valve.
201
MIKE URBAN Water Hammer
4.4.5 Tutorial
This section contains brief summary describing how to use Water Hammer
when creating a new project.
5. Add wave speed for each pipe in the model. Use the Field Calculator
from within the Pipe Editor to assign the wave speed globally. Example
of the command assigning the wave speed of 1200m/s to each pipe:
wavespeed=1200.
6. Edit specific data to pumps, valves, air-chambers, vented air-chambers,
or air valves where appropriate.
7. Run the water hammer analysis and load the analysis results.
8. Use standard way of reviewing the analysis results such as time series,
animated plans, or animated profile plots.
203
MIKE URBAN Water Hammer
The EPANET computer model used for water distribution network analysis is
composed of two parts: (1) the input data file and (2) the EPANET computer
program. The data file defines the characteristics of the pipes, the nodes
(ends of the pipe), and the control components (such as pumps and valves)
in the pipe network. The computer program solves the nonlinear energy
equations and linear mass equations for pressures at nodes and flowrates in
pipes.
205
EPANET
5.1 History
Pipe network analysis of water distribution systems has evolved from a time-
consuming process done infrequently to a quick and easy process done reg-
ularly on systems of all sizes.
Pipe network analysis initially started early in 1940. Years later, two network
analysis programs were introduced by Shamir and Howard (1968) and Epp
and Fowler (1970). Both programs used the Newton-Raphson method to line-
arize the nonlinear mass and energy equations. The major differences
between these two programs are:
Water quality models can also be used to study water retention time for reser-
voir operations, pipeline travel times, and the percentage of water supplied to
a location from multiple sources (i.e., treatment plants, wells, and reservoirs).
Additionally, water quality models can be used to develop a hydrant flushing
program to reduce water stagnation at dead ends within the pipe network.
And, site sampling locations, future rechlorination facility locations, cross-
connection locations, and reservoir operating strategies can be designed and
analyzed.
207
EPANET
Planning
MIKE URBAN WATER can be used in the planning of pipe network systems
to meet forecasted demands of the next 10 years or 20 years. For example,
the program can be used to develop long term capital-improvement plans for
the existing pipe network system. These plans can include staging, sizing,
and locating future pipe network and water chlorination facilities. The soft-
ware can also be used in the development of a main rehabilitation plan or a
system-improvement plan. And, a network analysis can provide suggestions
and recommendations to prepare for the occurrence of any unusual events.
Design
MIKE URBAN WATER can be used to design a new pipe network system or
improve on the existing pipe network system. For example, the analysis con-
ducted using MIKE URBAN WATER could help users in selecting and sizing
pipe network components, such as pipes, booster pumps, and pressure regu-
lating valves. As a part of the analysis, the performance of the pipe network
system can be analyzed to verify that the system satisfies fire-flow demand
requirements.
Operation
The operating status of a water distribution system (e.g., the pipe flowrates
and junction node pressures) can be determined by MIKE URBAN WATER.
The analysis can then be used to develop operational strategies based on the
guidelines for maximum use of available water and efficient management of
electrical energy. MIKE URBAN WATER can also be used for system trouble-
shooting, such as finding the location of a pipe break.
5.1.3 Skeletonization
In the past, water distribution models have not included all of the pipes con-
tained in the network system due to the fact that the numerical modeling
schemes used and the memory requirements required could handle only a
limited number of pipes and nodes. These limitations required skeletonization
of the pipe network system, where only a subset of all the pipes contained
within the network system was defined. However, these limitations have
eroded over the past years as enhanced programming methods and
increased computer hardware capabilities have come into being.
Skeleton network models typically include only those pipes that are consid-
ered significant to the flow and distribution of water. For example, a skeleton
model might only consider 12 inch and larger diameter mains. Smaller diame-
ter mains might also be included if they supply water to a significant area or
complete a loop in the network. Examples of a complete network model and
an equivalent skeleton network model are shown in Figure 5.1 and
Figure 5.2.
Pipe Diameter
6" and smaller
8"
12"
16"
Figure 5.1 A complete water distribution network, showing all pipes contained
within the network
Pipe Diameter
6" and smaller
8"
12"
16"
Figure 5.2 An equivalent skeleton water distribution network, showing those pipes
that are considered significant to the flow and distribution of water for a
steady state simulation
A major advantage with working with a skeleton model is that they are much
easier to define since there is less data involved, and the simulation time is
shorter since not as many pipes are involved in computing a solution. Also,
209
EPANET
the displaying of results is quicker and more readable since less pipes are
involved. In addition, the number of iterations required to converge to a solu-
tion is typically less as well.
As was stated before, a skeleton model should not be used when performing
water quality modeling. However, in other situations, the modeler must decide
when it is appropriate to skeletonize a network model to produce results with
sufficient accuracy to meet the modeling requirements. Therefore, when and
where to skeletonize a network must be decided upon a case by case basis.
Pipes
Pumps
Valves
Nodes, besides representing the connection point between pipes, can repre-
sent the following components in a network:
How the EPANET program models the hydraulic behavior of each of these
components is described in the following sections. All flow rates in this dis-
cussion will be assumed as cubic feet per second (cfs), although the program
can also accept flow rates in gallons per minute (gpm), million gallons per day
(mgd), and litres per second (L/s).
5.2.1 Pipes
Every pipe is connected to two nodes at its ends. In a pipe network system,
pipes are the channels used to convey water from one location to another.
The physical characteristics of a pipe include the length, inside diameter,
roughness coefficient, and minor loss coefficient. The pipe roughness coeffi-
cient is associated with the pipe material and age. The minor loss coefficient
is due to the fittings along the pipe.
When water is conveyed through the pipe, hydraulic energy is lost due to the
friction between the moving water and the stationary pipe surface. This fric-
211
EPANET
tion loss is a major energy loss in pipe flow and is a function of flowrate, pipe
length, diameter, and roughness coefficient.
The head lost to friction associated with flow through a pipe can be
expressed in a general fashion as:
b
h L = aq (5.1)
where
– hL = head loss, ft
– q = flow, cfs
– a = a resistance coefficient
– b = a flow exponent
EPANET can use any one of three popular forms of the headloss formula
shown in Equation 5.1: the Hazen-Williams formula, the Darcy-Weisbach for-
mula, or the Chezy-Manning formula. MIKE URBAN WATER allows the user
to choose the formulation to use.
The Hazen-Williams formula is probably the most popular head loss equation
for water distribution systems, the Darcy-Weisbach formula is more applica-
ble to laminar flow and to fluids other than water, and the Chezy-Manning for-
mula is more commonly used for open channel flow. Table 5.1 lists resistance
coefficients and flow exponents for each formula. Note that each formula
uses a different pipe roughness coefficient, which must be determined empir-
ically. Table 5.2 lists general ranges of these coefficients for different types of
new pipe materials. Be aware that a pipe's roughness coefficient can change
considerably with age.
Notes:
Pipes can contain check valves in them that restrict flow to a specific direc-
tion. They can also be made to open or close at pre-set times, when tank lev-
els fall below or above certain set-points, or when nodal pressures fall below
or above certain set-points. The normal initial condition for a pipe containing a
check valve or a pump is to be in open mode. The pipe will then switch to
closed mode only when flow is reversed.
In addition to the energy loss caused by friction between the fluid and the
pipe wall, energy losses also are caused by obstructions in the pipeline,
changes in flow direction, and changes in flow area. These losses are called
minor losses because their contribution to the reduction in energy is usually
much smaller than frictional losses. Head loss, which is the sum of friction
loss and minor losses, reduces the flowrate through the pipe.
213
EPANET
Bend points are inserted to separate two parallel pipes to allow them to be
selectable.
5.2.2 Pumps
A pump is a device that raises the hydraulic head of water. EPANET repre-
sents pumps as links of negligible length with specified upstream and down-
stream junction nodes.
Pumps are described with a pump characteristic curve. The pump curve
describes the additional head imparted to a fluid as a function of its flow rate
through the pump. EPANET is capable of modeling several types of user
defined pumps, constant energy pumps, single-point pump curves, three-
point pump curves, multiple point pump curves, and variable speed pump.
h G = 8,81Hp
------------------ (5.2)
q
where
– hG = head gain, ft
– Hp = pump horsepower
– q = flow, cfs
215
EPANET
For variable speed pumps, the pump curve shifts as the speed of the pump
changes. The relationships between flow (Q) and head (H) at speed N1 and
N2 are:
Q1 N1
-------- = ------- (5.3)
Q2 N2
2
H1 N
------- = ------1- (5.4)
H2 N2
EPANET will shut a pump down if the system demands a head higher than
the first point on a curve (i.e., the shutoff head.) A pump curve must be sup-
plied for each pump the system unless the pump is a constant energy pump.
Controlling Pumps
EPANET allows you to turn pumps on or off at pre-set times, when tank levels
fall below or rise above certain set-points, or when nodal pressures fall below
or rise above certain set-points. Variable speed pumps can also be consid-
ered by specifying using the Control Editor (see the section titled Control Edi-
tor on -74 for further information) that their speed setting be changed under
these same types of conditions. By definition, the original pump curve sup-
plied to the program has a relative speed setting of 1. If the pump speed dou-
bles, then the relative setting would be 2; if run at half speed, the relative
setting is 0.5 and so on. Figure 5.5 illustrates how changing a pump's speed
setting affects its characteristic curve.
217
EPANET
Figure 5.9 MIKE URBAN WATER can model both parallel pumps and pumps in
series
For example, as shown in Figure 5.10, if two pumps are connected in parallel
and one pump operates at a flow rate of 50 gpm with a head of 75 ft and the
other pump operates at a 60 gpm with a head of 75 ft, then the resultant sin-
gle equivalent pump curve flow rate of 110 gpm would be available with a
head of 75 ft.
200
150
Head (ft)
100 Equiv
ilent P
ump
Pu
50
mp
Pum
2
p1
0
0 25 50 75 100 125 150
Discharge (gpm)
Figure 5.10 When modeling pumps in parallel, an equivalent single pump curve can
be determined by adding together the characteristic pump curve dis-
charge values from the other individual pump curves
Similarly, as shown in Figure 5.11 if two pumps are connected in series and
one pump provides a head of 75 ft at a flow rate of 50 gpm and the other
pump provides a head of 65 feet at a flow rate of 50 gpm, then the resultant
single equivalent pump curve head of 140 ft would be available at a flow rate
of 50 gpm.
219
EPANET
200 Equiv
ilent P
ump
150
Head (ft)
Pump 1
100
Pump 2
50
0
0 25 50 75 100 125 150
Discharge (gpm)
Figure 5.11 When modeling pumps in series, an equivalent single pump curve can
be determined by adding together the characteristic pump curve head
values from the other individual pump curves
5.2.3 Valves
Aside from the valves in pipes that are either fully opened or closed (such as
check valves), EPANET can also represent valves that control either the
pressure or flow at specific points in a network. Such valves are considered
as links of negligible length with specified upstream and downstream junction
nodes. The types of valves that can be modeled are described below.
Minor head losses (also called local losses) can be associated with the added
turbulence that occurs at bends, junctions, meters, and valves. The impor-
tance of such losses will depend on the layout of the pipe network and the
degree of accuracy required. EPANET allows each pipe and valve to have a
221
EPANET
minor loss coefficient associated with it. It computes the resulting head loss
from the following formula:
2
0,0252Kq
h L = -------------------------- (5.5)
4
d
where
– hL = head loss, ft
– d = diameter in, ft
– K = loss coefficients
– q = flow rate, cfs
Table 5.3 gives values of K for several different kinds of pipe network compo-
nents.
Exit 1.0
5.2.5 Nodes
Nodes are the locations where pipes connect. Two types of nodes exist in a
pipe network system, (1) fixed nodes and (2) junction nodes. Fixed nodes are
nodes whose HGL are defined. For example, reservoirs and storage tanks
are considered fixed nodes, because their HGL are initially defined. Junction
nodes are nodes whose HGL are not yet determined and must be computed
in the pipe network analysis.
Degree of freedom, elevation, and water demand are the three important
input parameters for a node (see Figure 5.3). A node's degree of freedom is
the number of pipes that connect to that node. In EPANET, a junction node
may be connected to more than one pipe, but a fixed node (i.e., storage tank
or reservoir) must be connected to exactly one pipe. Therefore, a fixed node's
degree of freedom is always one, and a junction node's degree of freedom
may be greater than one. The elevation of a node can sometimes be obtained
from system maps or drawings. More often, it is approximated using topo-
graphic maps. Water demand at a junction node is the summation of all water
drawn from or added to the system at that node.
All nodes should have their elevation specified above sea level (i.e., greater
than zero) so that the contribution to hydraulic head due to elevation can be
computed. Any water consumption or supply rates at nodes that are not stor-
age nodes must be known for the duration of time the network is being ana-
lyzed. Storage nodes (i.e., tanks and reservoirs) are special types of nodes
where a free water surface exists and the hydraulic head is simply the eleva-
tion of water above sea level. Tanks are distinguished from reservoirs by hav-
ing their water surface level change as water flows into or out of them—
reservoirs remain at a constant water level no matter what the flow is.
EPANET models the change in water level of a storage tank with the following
equation:
q
y = ---- t (5.6)
A
where
For each storage tank, the program must know the cross-sectional area as
well as the minimum and maximum permissible water level. Reservoir-type
storage nodes are usually used to represent external water sources, such as
lakes, rivers, or groundwater wells. Storage nodes should not have an exter-
nal water consumption or supply rate associated with them.
223
EPANET
Storage Tank
Reservoir
Pump 11 12 13
11 12
Equivalent Pumped 21 22
Groundwater Well 21 23
22
121 122
31
31 32
As the storage tank empties and its water level falls to a certain amount, con-
trol rules turn on the groundwater (reservoir) pump to start refilling of the stor-
age tank. The control rules used to regulate the starting and stopping of the
groundwater pump are defined within the Control Editor (see the section titled
Control Editor on - 74 for further information).
Head
Static
ped
Pum
Discharge
Figure 5.13 The groundwater well pump curve may need to be adjusted downward
to account for groundwater well drawdown and recharge effects
225
EPANET
1.6
1.4 1.4
Usage Factor
1.2 1.2
1.1
1.0
1 0.9
0.8
0.7
0.6
;;;;;
0.5
0
0 1 2 3 4 5 6 7 8 9 10 11 12
Demand Period
Assume we have a pipe network with N junction nodes and NF fixed grade
nodes (tanks and reservoirs). Let the flow-headloss relation in a pipe between
nodes i and j be given as:
H i – H j = h ij = r Q n ij + mQ 2 ij (5.7)
below). For pumps, the headloss (negative of the head gain) can be repre-
sented by a power law of the form
n
hij 2 (h0 r (Qij / ) )
(5.8)
where h0 is the shutoff head for the pump, w is a relative speed setting, and r
and n are the pump curve coefficients. The second set of equations that must
be satisfied is flow continuity around all nodes:
Qij – Di = 0 (5.9)
j
for i = 1,... N.
where Di is the flow demand at node i and by convention, flow into a node is
positive. For a set of known heads at the fixed grade nodes, we seek a solu-
tion for all heads Hi and flows Qij that satisfy Eqs. (5.8) and (5.9).
The Gradient solution method begins with an initial estimate of flows in each
pipe that may not necessarily satisfy flow continuity. At each iteration of the
method, new nodal heads are found by solving the matrix equation.
AH = F (5.10)
Aii pij
j
(5.11)
Aij pij
(5.12)
where pij is the inverse derivative of the headloss in the link between nodes i
and j with respect to flow. For pipes,
1
pij n 1
nr Qij 2m Qij
(5.13)
227
EPANET
1
P ij = ----------------------------------
- (5.14)
Q ij n – 1
n r -------
2
Each right hand side term consists of the net flow imbalance at a node plus a
flow correction factor:
1
pij n 1
nr Qij 2m Qij
(5.15)
where the last term applies to any links connecting node i to a fixed grade
node f and the flow correction factor yij is: P
n 2
Y ij = p r Q ij + m Q ij sgn Q ij (5.16)
yij pij 2 h0 r (Qij / ) n (5.17)
for pumps, where sgn (x) is 1 if x > 0 and -1 otherwise. (Qij is always positive
for pumps.)
After new heads are computed by solving Eq. (5.10), new flows are found
from:
Q ij = Qij – y ij – p ij H i – H j (5.18)
If the sum of absolute flow changes relative to the total flow in all links is
larger than some tolerance (e.g., 0.001), then Eqs. (5.9) and (5.10) are
solved once again. The flow update formula (5.11) always results in flow con-
tinuity around each node after the first iteration.
of the analysis.
2. For the very first iteration at time 0, the flow in a pipe is chosen equal to
the flow corresponding to a velocity of 1 ft/sec. The flow through a con-
ventional pump equals the design flow specified for the pump curve. An
initial flow of 1 cfs is assumed for other types of pumps. (All computa-
tions are made with head in feet and flow in cfs).
64-
f = ------- (5.19)
Re
0.25
f 2
5.74
Ln 0.9
3.7d Re (5.20)
Cubic Interpolation From Moody Diagram for 2,000 < Re < 4,000 (Dun-
lop, 1991):
f = X1 + R X2 + R X3 + X4 (5.21)
Re
R = -------- (5.22)
64
X1 = 7FA – FB (5.23)
229
EPANET
FA = Y3 –2 (5.28)
0,00514215
FB = FA 2 – ---------------------------- (5.29)
Y2 Y3
5,74
Y2 = ---------- + ------------- (5.30)
3,7d Re 0,9
+ ----------------
Y3 = – 0,86859Ln ---------- 5,74 -
(5.31)
3,7d 0,9
4000
4. The minor loss coefficient based on velocity head (K) is converted to one
based on flow (m) with the following relation:
0,02517K
m = ----------------------- (5.32)
4
d
6. Open valves are assigned an r-value by assuming the open valve acts
as a smooth pipe (f = 0.02) whose length is twice the valve diameter.
Closed links are assumed to obey a linear headloss relation with a large
resistance factor, i.e., h = 108Q, so that p = 10-8 and y = Q. For links
where (r+m) Q < 10-7, p = 107 and y = Q/n.
7. Status checks on pumps, check valves (CVs), flow control valves, and
pipes connected to full/empty tanks are made after every other iteration,
up until the 10th iteration. After this, status checks are made only after
convergence is achieved. Status checks on pressure control valves
(PRVs and PSVs) are made after each iteration.
8. During status checks, pumps are closed if the head gain is greater than
the shutoff head (to prevent reverse flow). Similarly, check valves are
closed if the headloss through them is negative (see below). When these
conditions are not present, the link is re-opened. A similar status check is
made for links connected to empty/full tanks. Such links are closed if the
difference in head across the link would cause an empty tank to drain or
a full tank to fill. They are re-opened at the next status check if such con-
ditions no longer hold.
10. If the status check closes an open pump, pipe, or CV, its flow is set to 10-
6 cfs. If a pump is re-opened, its flow is computed by applying the current
head gain to its characteristic curve. If a pipe or CV is re-opened, its flow
is determined by solving Eq. (5.9) for Q under the current headloss h,
ignoring any minor losses.
11. Matrix coefficients for pressure breaker valves (PBVs) are set to the fol-
lowing: p = 108 and y = 108Hset, where Hset is the pressure drop setting
for the valve (in feet). Throttle control valves (TCVs) are treated as pipes
with r as described in item 6 above and m taken as the converted value
of the valve setting (see item 4 above).
12. Matrix coefficients for pressure reducing, pressure sustaining, and flow
control valves (PRVs, PSVs, and FCVs) are computed after all other
links have been analyzed. Status checks on PRVs and PSVs are made
as described in item 7 above. These valves can either be completely
open, completely closed, or active at their pressure or flow setting.
231
EPANET
– if Hi > Hj + Htol
– and Hi < Hset - Htol then new status = OPEN
– if Hi > Hj + Htol
– and Hj < Hset - Htol then new status = ACTIVE
else new status = CLOSED
where Q is the current flow through the valve, Hi is its upstream head, Hj
is its downstream head, Hset is its pressure setting converted to head,
Hml is the minor loss when the valve is open (= mQ2), and Htol and Qtol
are the same values used for check valves in item 9 above. A similar set
of tests is used for PSVs, except that when testing against Hset, the i
and j subscripts are switched as are the > and < operators.
14. Flow through an active PRV is maintained to force continuity at its down-
stream node while flow through a PSV does the same at its upstream
node. For an active PRV from node i to j:
– Pij = 0
– Fj = Fj + 108Hset
– Ajj = Ajj + 108
This forces the head at the downstream node to be at the valve setting
Hset. An equivalent assignment of coefficients is made for an active PSV
except the subscript for F and A is the upstream node i. Coefficients for
open/closed PRVs and PSVs are handled in the same way as for pipes.
15. For an active FCV from node i to j with flow setting Qset, Qset is added
to the flow leaving node i and entering node j, and is subtracted from Fi
and added to Fj. If the head at node i is less than that at node j, then the
valve cannot deliver the flow and it is treated as an open pipe.
17. For extended period simulation (EPS), the following procedure is imple-
mented:
After a solution is found for the current time period, the time step for the
next solution is the minimum of:
– The time until a new demand period begins.
– The shortest time for a tank to fill or drain.
– The shortest time until a tank level reaches a point that triggers a
change in status for some link (e.g., opens or closes a pump) as
stipulated in a simple control.
– The next time until a simple timer control on a link kicks in.
– The next time at which a rule-based control causes a status change
somewhere in the network.
In computing the times based on tank levels, the latter are assumed to
change in a linear fashion based on the current flow solution. The activa-
tion time of rule-based controls is computed as follows:
– Starting at the current time, rules are evaluated at a rule time step.
Its default value is 1/10 of the normal hydraulic time step (e.g., if
hydraulics are updated every hour, then rules are evaluated every 6
minutes).
– Over this rule time step, clock time is updated, as are the water lev-
els in storage tanks (based on the last set of pipe flows computed).
– If a rule's conditions are satisfied, then its actions are added to a list.
If an action conflicts with one for the same link already on the list
then the action from the rule with the higher priority stays on the list
and the other is removed. If the priorities are the same then the orig-
inal action stays on the list.
– After all rules are evaluated, if the list is not empty then the new
actions are taken. If this causes the status of one or more links to
change then a new hydraulic solution is computed and the process
begins anew.
– If no status changes were called for, the action list is cleared and the
next rule time step is taken unless the normal hydraulic time step
has elapsed.
Time is advanced by the computed time step, new demands are found,
tank levels are adjusted based on the current flow solution, and link con-
trol rules are checked to determine which links change status.
A new set of iterations with equations are begun at the current set of flows.
Water age is the time spent by a parcel of water in the network. It provides a
simple, non-specific measure of the overall quality of delivered drinking water.
New water entering the network from reservoirs or source nodes enters with
age of zero. As this water moves through the pipe network it splits apart and
blends together with parcels of varying age at pipe junctions and storage
facilities. EPANET provides automatic modeling of water age. Internally, it
233
EPANET
Source tracing tracks over time what percent of water reaching any node in
the network had its origin at a particular node. The source node can be any
node in the network, including storage nodes. Source tracing is a useful tool
for analyzing distribution systems drawing water from two or more different
raw water supplies. It can show to what degree water from a given source
blends with that from other sources, and how the spatial pattern of this blend-
ing changes over time. EPANET provides an automatic facility for performing
source tracing. The user need only specify which node is the source node.
Internally, EPANET treats this node as a constant source of a non-reacting
constituent that enters the network with a concentration of 100.
A dissolved substance will travel down the length of a pipe with the same
average velocity as the carrier fluid while at the same time reacting (either
growing or decaying) at some given rate. Longitudinal dispersion is usually
not an important transport mechanism under most operating conditions. This
means there is no intermixing of mass between adjacent parcels of water
traveling down a pipe. Advective transport within a pipe is represented with
the following equation:
C i C
= – ui i + r C i (5.33)
t x
At junctions receiving inflow from two or more pipes, the mixing of fluid is
taken to be complete and instantaneous. Thus the concentration of a sub-
stance in water leaving the junction is simply the flow-weighted sum of the
concentrations from the inflow pipes. For a specific node k one can write:
V C =
t s s iIsQi C i x = Li – jOsQj C s + r Cs (5.35)
r kC n (5.36)
where k = a reaction constant and n = the reaction order. When a limiting con-
centration exists on the ultimate growth or loss of a substance then the rate
expression becomes
235
EPANET
R K b (C L C )C ( n 1)
(5.37)
R K b (C C L )C ( n 1)
(5.38)
R K bC
(5.39)
R K b (C L C )
(5.40)
R K bC (C C L )
(5.41)
Kb C
R = ----------------
- (5.42)
CL – C
1,365 100UVA
K b = – 0,32UVA – -------------------------- (5.43)
DOC
CL 4.98UVA 1.91DOC
(5.44)
Note: These expressions apply only for values of Kb and CL used with
Michaelis-Menton kinetics.
R = 1.0
This special case can be used to model water age, where with each unit of
time the "concentration" (i.e., age) increases by one unit.
The relationship between the bulk rate constant seen at one temperature (T1)
to that at another temperature (T2) is often expressed using a van't Hoff -
Arrehnius equation of the form:
K b 2 K b1 T 2T 1
(5.45)
237
EPANET
ence the overall rate of this reaction. The surface area per unit volume, which
for a pipe equals 2 divided by the radius, determines the former factor. The
latter factor can be represented by a mass transfer coefficient whose value
depends on the molecular diffusivity of the reactive species and on the Reyn-
olds number of the flow (Rossman et. al, 1994). For first-order kinetics, the
rate of a pipe wall reaction can be expressed as:
2k w k f C
r = -------------------------
- (5.46)
R kw + kf
r MIN (k w , k f C )(2 / R)
(5.47)
D
k f Sh
d (5.48)
0.0668(d / L) Re Sc
Sh 3.65
1 0.04(d / L) Re Sc
2/3
(5.49)
together at junctions between fixed-length time steps (Liou and Kroon, 1987).
These water quality time steps are typically much shorter than the hydraulic
time step (e.g., minutes rather than hours) to accommodate the short times of
travel that can occur within pipes. As time progresses, the size of the most
upstream segment in a pipe increases as water enters the pipe while an
equal loss in size of the most downstream segment occurs as water leaves
the link. The size of the segments in between these remains unchanged.
The following steps occur at the end of each such time step:
1. The water quality in each segment is updated to reflect any reaction that
may have occurred over the time step.
2. The water from the leading segments of pipes with flow into each junc-
tion is blended together to compute a new water quality value at the junc-
tion. The volume contributed from each segment equals the product of its
pipe's flow rate and the time step. If this volume exceeds that of the seg-
ment then the segment is destroyed and the next one in line behind it
begins to contribute its volume.
3. Contributions from outside sources are added to the quality values at the
junctions. The quality in storage tanks is updated depending on the
method used to model mixing in the tank (see below).
4. New segments are created in pipes with flow out of each junction, reser-
voir, and tank. The segment volume equals the product of the pipe flow
and the time step. The segment's water quality equals the new quality
value computed for the node.
To cut down on the number of segments, Step 4 is only carried out if the new
node quality differs by a user-specified tolerance from that of the last seg-
ment in the outflow pipe. If the difference in quality is below the tolerance then
the size of the current last segment in the outflow pipe is simply increased by
the volume flowing into the pipe over the time step.
This process is then repeated for the next water-quality time step. At the start
of the next hydraulic time step the order of segments in any links that experi-
ence a flow reversal is switched. Initially each pipe in the network consists of
a single segment whose quality equals the initial quality assigned to the
downstream node.
239
EPANET
EPANET can use four different types of models to characterize mixing within
storage tanks: Complete Mixing, Two-Compartment Mixing, First In First Out
(FIFO) Plug Flow and Last In Fist Out (LIFO) Plug Flow. Different models can
be used with different tanks within a network.
Complete Mixing assumes that all water that enters the tank is instantane-
ously and completely mixed with the water already in the tank.
Two-Compartment Mixing divides the available storage volume into two com-
partments.
The inlet/outlet pipes of the tank are assumed to be in the first compartment.
New water that enters the tank is mixed with the first compartment. If this
compartment is full, then it sends the overflow water to the second compart-
ment where it completely mixes with the water stored there. When water
leaves the tank, it exits from the first compartment, which if full, receives an
241
EPANET
equal volume of water from the second compartment to make up the differ-
ence. The first compartment is capable of simulating short circulating
between inflow and outflow, while the second compartment represents dead
zones. The user must define a single parameter, which is the fraction of the
volume of the tank to be dedicated to the first compartment.
FIFO Plug Flow assumes that there is no mixing during the residence time of
the water in the tank.
Water parcels move through the tank in a segregated fashion, where the first
parcel to enter the tank is also the first to leave. Physically speaking, this
model best represents baffled tanks with constant inflow and outflow. There
are no additional parameters to describe in this model.
LIFO Plug Flow assumes that there is no mixing between parcels of water
that enter the tank.
However, in contrast to the FIFO Plug Flow model, the parcels stack up one
on top of another, where water enters and leaves from the bottom. In this
fashion, the first parcel to enter the tank is also the last to leave. Physically
speaking this type of model might apply to a tall narrow standpipe with an
inflow/outflow pipe at the bottom and a low momentum inflow. It requires no
additional parameters be provided
This section discusses the concepts and steps involved in calibrating a water
distribution model.
The computed pressure and measured field (actual) pressure will not exactly
match for every node contained within the network system—there will be dif-
ferences within the system. However, the maximum amount of these differ-
ences needs to be considered when performing a model calibration.
243
EPANET
Although the head loss method is better than the other methods in determin-
ing accuracy, from a practical standpoint the pressure differential method and
the percent pressure differential method are far easier to use and understand.
Therefore, since the main point is to establish a common unit of measure to
determine the degree of accuracy reached when calibrating a model, any of
these methods can be used.
During the calibration process there should be a point where the modeler has
decided that enough time has been spent on calibrating the model to the
actual network and he should move on to analyzing the pipe network system.
Model Credibility
Through the process of calibrating a model, credibility of the model is estab-
lished. The data and modeling assumptions of an uncalibrated model are
unlikely to match the actual system. On the other hand, a calibrated model is
known to simulate a network system for a range of operating conditions. Its
input data has been examined and adjusted to insure that the model can be
used as an accurate, predictive tool. Hence, use of an uncalibrated model is
not good engineering practice, since it will most likely lead to inaccurate
model results and poor engineering decisions based upon these results.
Benchmark
Once a network model has been calibrated to a known range of operating
conditions, it can be used as a benchmark. Pressure and flow rates com-
puted by this model become the benchmark from which pressure and flow
rates computed by subsequent, modified models can be compared. The dif-
ferences between the two models can then be used to analyze the changes
brought about in the modified system.
Uncover Errors
Calibration requires collecting data about the system. Many times during the
process of calibrating a model, questionable model results are investigated.
Inconsistencies between the modeled results and the actual field conditions
are examined, with additional field data being collected and analyzed. Incor-
rect pipe diameters may be determined, or even incorrectly closed valves dis-
covered.
245
EPANET
The data requirements for performing a model calibration is made up with the
fact that some data is fixed and unchanging (e.g., pipe diameter, length, etc.),
some data is variable with time (e.g., demand patterns, pump rates, dis-
charge pressures, reservoir elevations, etc.), and some data is assumed
(e.g., pipe roughness values, consumption rates, etc.). Some data, such as
consumption, is measured but sometimes also assumed.
The availability of data is also a limiting factor. For example, calibrating to the
maximum hour requires data for only one time period. Calibration should also
be performed for the minimum hour of the maximum day. Therefore, data
would be required for two time periods. Extending the calibration process to
an extended period simulation for the maximum day increases the data
requirements to 24 time periods.
What can be done to simplify the calibration process is to approach the cali-
bration process incrementally. To elaborate, first calibrate the model for the
maximum hour. This will provide a good understanding of the network sys-
tem, and can usually be accomplished with a reasonable degree of accu-
racy—even for large network systems. Then, the minimum hour can be
calibrated. Knowledge learned in calibrating the maximum hour can be used
to speed-up the calibration process for the minimum hour. Once both of these
simulations have been calibrated, calibration of the maximum day extended
period simulation can be performed if time permits and the data is readily
available.
Data Acquisition
Calibration data should be compiled and organized so that it can be used in
an efficient manner. In addition, a map of the network system is essential to
properly calibrate the system. All pipes, pumps, valves, tanks, and reservoirs
should be identified, including nodal elevations, pressure zone boundaries,
and other important information.
SCADA data and graphs are useful since they contain data that will be refer-
enced when developing the calibrated model. Even though the maximum
hour is typically of primary interest, values for the entire day are also impor-
tant since they provide an overall understanding of the operation of the net-
work system. This information, for example, allows the user to compare the
discharge pressure value at a particular location on the system during the day
with the maximum hour value.
Operational Review
Operational rules must be determined for the network system for all major
system components. Discussions with system operators should be performed
to determine the operational rules and strategy behind how specific opera-
tions are performed. For example, under what conditions do system opera-
tors turn on a pump, close a control valve, or adjust a pressure-regulating
valve? Is water pressure at particular key locations in the network act as a
flag to turn on a particular pump? Are pumping schedules ever changed to try
to minimize power costs? Are all facilities currently available for use, or are
some facilities off-line for maintenance repair?
Recalibration Frequency
Whenever changes occur in the water distribution network, such as opera-
tional changes, network configuration, or increases in water consumption, the
degree of accuracy for the calibrated model is reduced. If these changes are
severe enough, the model will need to be recalibrated. In practice this means
that the model should be recalibrated whenever major new facilities are
added to the network system, a new record for maximum hour is set, or oper-
ational procedures change significantly. Therefore, recalibration may be nec-
essary every few years.
247
EPANET
Once all of the necessary data has been gathered, the network model can be
constructed. It is important that the network model accurately represents the
physical layout of the system. Pipes and nodes must be accurately located in
the model. Pipe roughness values should be estimated—based upon the age
of the pipe.
Next, consumption values must be defined at the nodes. To save time initially,
MIKE URBAN WATER provides a method of globally applying the total net-
work demand to each of the individual network nodes using the Distributed
Demand dialog box. See the section titled Distributed Demands on - 59 for
more information on computing distributed demands. Also, operational data
must be defined for all pumps, valves, storage tanks, and reservoirs. Finally,
a network simulation can be performed with EPANET.
Initial Simulation
An initial simulation is performed to simply determine what the resulting pres-
sures and flow rates are in the pipe network. This simulation may be simpli-
fied by using single operational values rather than complete operational data.
At this stage of the network calibration, exact operational data is not required.
For example, a reservoir elevation could be entered as rounded to the near-
est whole number elevation, rather than to a tenth of a foot. Similarly, a single
value pump definition could be used rather than a pump curve for each of the
pumps defined in the network. As long as the simplified data is generally
accurate, a balanced run can be produced so that the pipe network and oper-
ational input data can be checked.
Table 5.4 Comparison of actual field flow rates and computed flow rates
249
EPANET
Once the differences between the computed and actual measured values can
be determined, adjustments to the model data can be performed to make the
computed results match more closely to that of the actual data. When per-
forming comparisons, differences in flow rates and pressures should be
examined.
The computed flow rates and pressures may be lower or higher from the
measured values. Table 5.6 allows you to quickly determine the cause for
these differences—allowing you to adjust the model input to more closely
match the measured values.
Table 5.6 Causes for differences between computed and mesured flow rates and
pressures
In adjusting the model input data, there are four areas where adjustments can
be considered. These input data adjustments are discussed in the following
sections.
251
EPANET
For a small network system with fairly reliable input data, a calibrated model
might be achieved within 10 or fewer iterations. For more complicated net-
work systems, a calibrated model may require more than 100 iterations.
253
EPANET
AWWA Manual M32, Distribution Network Analysis for Water Utilities, First
Edition, American Water Works Association, Denver, CO, 1989.
Babbitt, H. E., Plumbing, Third Edition, McGraw-Hill, Inc., New York, NY,
1960.
Bhave, P. R., Rules for Solvability of Pipe Networks, Journal of Indian Water
Works Association, January-March, 1990, pp. 7-10.
Bishop, M. M., Morgan, J. M., Cornwell, B., and Jamison, D. K., Improving the
Disinfection Detention Time of a Water Plant Clearwell, Journal of American
Water Works Association, Vol. 85, No. 3, 1993, pp. 68-75.
Boulos, P. F., Altman, T., and Liou, J. C. P., On the Solvability of Water Distri-
bution Networks with Unknown Pipe Characteristics, Journal Applied Mathe-
matical Modelling, Vol. 17, July, 1993, pp. 380-387.
Boulos, P. F., Altman, T., and Sadhal, K. S., Computer Modelling of Water
Quality in Large Multiple Source Networks, Journal Applied Mathematical
Modelling, Vol. 16, 1992, pp. 439-445.
Boulos, P. F., Altman, T., Bowcock, R. W., Dhingra, A. K., and Collevati, F., An
Explicit Algorithm for Modeling Distribution System Water Quality with Appli-
cations, 2nd BHR International Conference on Water Pipeline Systems, D. S.
Miller ed., Edinburgh, Scotland, 1994, pp. 405-423.
Boulos, P. F., Altman, T., Jarrige, P. A., and Collevati, F., An Event-Driven
Method for Modeling Contaminant Propagation in Water Networks, Journal
Applied Mathematical Modelling, Vol. 18, 1994, pp. 84-92.
Boulos, P. F., Altman, T., Jarrige, P. A., and Collevati, F., Discrete Simulation
Approach for Network Water Quality Models, Journal Water Resources Plan-
ning and Management ASCE, Vol.121 No. 1, 1995, pp. 49-60.
Boulos, P. F., and Altman, T., An Explicit Approach for Modelling Closed
Pipes in Water Networks, Journal Applied Mathematical Modelling, Vol. 17,
August, 1993, pp. 437-443.
Boulos, P. F., and Altman, T., Explicit Calculation of Water Quality Prameters
in Pipe Distribution Networks, Journal Civil Engineering Systems, Vol. 10,
1993, pp. 187-206.
Boulos, P. F., and Ormsbee, L. E., Explicit Network Calibration for Multiple
Loading Conditions, Journal Civil Engineering Systems., Vol. 8, 1991, pp.
153-160.
Boulos, P. F., and Wood, D. J., An Explicit Algorithm for Calculating Operating
Parameters for Water Networks, Journal Civil Engineering System, Vol. 8,
1991, pp. 115-122.
Clark, R. M., Grayman, W. M., Goodrich, J. A., Deininger, R. A., and Hess, A.
F., Field-Testing Distribution Water Quality Models, Journal of American
Water Works Association, Vol. 83, No. 7, 1991, pp. 67-75.
Clark, R. M., Smalley, G., Goodrich, J., Tull, R., Rossman, L. A., Vasconcelos,
J. J., and Boulos, P. F., Managing Water Quality in Distribution Systems: Sim-
ulating TTHM and Chlorine Residual Propagation, Journal Water SRT-AQUA,
Vol. 43, No. 4, 1994, pp. 182-191.
255
EPANET
Cole, G. B., Leak Detection -- Two Methods That Work, AWWA Seminar on
Operational Techniques for Distribution Systems, Los Angeles, California,
1980.
Daugherty, R.L. and Franzini, J.B. and Finnemore, E.J., "Fluid Mechanics
with Engineering Applications".
Duff, I. S., Erisman, A. M., and Reid, J. K., Direct Methods for Sparse Matri-
ces, Oxford University Press, Inc., New York, NY, 1992.
E. Cabrera and F. Martínez, "Water Supply Systems - State of the art and
future trends".
Edwards, D. K., Denny, V. E., and Mills, A. F., Transfer Processes, McGraw-
Hill, New York, NY, 1976.
Epp, R., and Fowler, A. G., Efficient Code for steady state Flows in Networks,
Journal of the Hydraulics Division, Proceedings of the American Society of
Civil Engineers, Vol. 96, No. HY1, January, 1970, pp. 43-56.
Fok, A.T.K, "Design Charts for Air Chamber on Pump Pipe Lines".
George, A., and Liu, J., Computer Solution of Large Sparse Positive Definite
Systems, Prentice-Hall, Inc., Englewood Cliffs, New Jersey, 1981.
Gofman, E. and Rodeh, M., Loop Equations with Unknown Pipe Characteris-
tics, Journal of the Hydraulics Division, ASCE, Vol. 107, No. HY9, September,
1981, pp. 1047-1060.
Grayman, W. M., and Clark, R. M., Using Computer Models to Determine the
Effect of Storage on Water Quality, Journal of American Water Works Associ-
ation, Vol. 85, No. 7, 1993, pp. 67-77.
Hamam, Y. M., and Brameller, A., Hybrid Method for the Solution of Piping
Networks, Proceedings IEE, Vol. 118, No. 11, 1971, pp. 1607-1612.
Housing and Home Finance Agency, The Uniform Plumbing Code for Hous-
ing, Section 8.3.5, February, 1948.
Chow, V. T., Maidmant, D. R., and Mays, L. W., Applied Hydrology, McGraw-
Hill, Inc., New York, NY, 1988.
Jeppson, R. W., Analysis of Flow in Pipe Networks, 5th Printing, Ann Arbor
Science, Ann Arbor, Michigan, 1981.
Journal of the Hydraulic Division, Vol. 104, No. HY9, September, 1978, pp.
1289-1303.
Journal of the Hydraulic Division, Vol. 107, No. HY10, October, 1981, pp.
1145-1161.
Lansey, K. E. and Mays, L. W., Optimization Model for Water Distribution Sys-
tem Design, Journal of the Hydraulic Engineering, Vol. 115, No. 10, October,
1989, pp. 1401-1418.
257
EPANET
Nayyar, M. L., Piping Handbook, Sixth Edition, McGraw-Hill, Inc., New York,
NY, 1992.
Nielsen, H. B., Methods for Analyzing Pipe Networks, Journal of the Hydrau-
lic Engineering, Vol. 115, February, 1989, pp. 139-157.
Ormsbee, L. E. and Wood, D. J., Hydraulic Design Algorithms for Pipe Net-
works, Journal of Hydraulic Engineering, Vol. 112, No. 12, December, 1986,
pp. 1195-1207.
Osiadacz, A. J., Simulation and Analysis of Gas Networks, E & F Spon Ltd.,
London, UK, 1991.
Pierre, D. A., Optimization Theory with Application, Wiley, New York, NY,
1969.
Press, W. H., Teukolsky, S. A., Vetterling, W. T., and Flannery, B. P., Numeri-
cal Recipes in C, Second Edition, Cambridge University, New York, NY, 1992.
325 class handout, summer semester, 1991, Alan K. Richardson is with Black
& Veatch Consulting Engineers at Kansas City, Missouri.
Rossman, L. A., and Boulos, P. F., Numerical Methods for Modeling Water
Quality in Distribution Systems: A Comparison, Journal Water Resources
Planning and Management ASCE, Vol. 122, No. 2, 1996, pp 137-146.
Rossman, L. A., Clark, R. M., and Grayman, W. M., Modeling Chlorine Resid-
uals in Drinking-Water Distribution Systems, Journal Environmental Engi-
neering ASCE, Vol. 120, No. 4, 1994, pp. 803-820.
Rossman, L. A., EPANET -An Advanced Water Quality Modeling Package for
Distribution System, Proceedings of American Water Works Association 1993
Annual Conference.
Shamir, U., and Howard, D. D., Water Distribution Systems Analysis, Journal
of the Hydraulic Division, ASCE, Vol. 94, No. HY1, January, 1968, pp. 219-
234.
Schweitzer, P. A., Handbook of Valves, Industrial Press, Inc., New York, NY,
1972.
Streeter, V. L., and Wylie, E. B., Fluid Mechanics, McGraw-Hill, Inc., Seventh
Edition, 1979.
Todini, E. and Pilati, S., "A Gradient Method for the Analysis of Pipe Net-
works".
Todini, E., and Pilati, S., A Gradient Method for the Analysis of Pipe Net-
works, Proceeding of the International Conference on Computer Applications
for Water Supply and Distribution, Leicester Polytechnic, UK, September 8-
10, 1987.
259
EPANET
Travers, K., The Mesh Method in Gas Network Analysis, Gas Journal,
November, 1967, pp. 167-177.
Verwey, A. and Yu, J.H., "A Space-Compact High-Order Scheme for Water
Hammer Simulations".
Walski, T. M., Gessler, J., and Sjostrom, J. W., Water Distribution Systems:
Simulation and Sizing, Lewis Publishers, Inc., 1990.
Warring, R. H., Handbook of Valves, Piping and Pipelines, First Edition, Trade
& Technical Press Ltd., England, 1982.
261
Index
B D
Base Elevation . . . . . . . . . . . 42, 47 Darcy-Weisbach . . . . . . . . . . . . 139
Benchmark . . . . . . . . . . . . . . . 245 Data Acquisition . . . . . . . . . . . . 247
Boundary Conditions . . . . . . . . 58, 182 Data Source . . . . . . 13, 19, 26, 35, 41
Boundary Curve . . . . . . . . . . . . 59 Defining Curved Pipes . . . . . . . . . .24
Boundary Type . . . . . . . . . . . . . 59 Defining the Model . . . . . . . . . . . 111
Branches . . . . . . . . . . . . . . . . 184 Definition of Network Layout . . . . . . 183
Browser . . . . . . . . . . . . . . . . . 130 Demand . . . . . .52, 65, 70, 71, 72, 129
Browsing Results . . . . . . . . . . . . 138 Demand Allocation . . . . . . . . . 64, 124
Bulk Flow Reactions . . . . . . . . . . 235 Demand Category . . . . . . . 15, 66, 70
Bulk Reaction Order . . . . . . . . . . 96 Demand Coefficient . . . . . . 14, 21, 70
Demand Description . . . . . . . . . . .15
C Demand Editing . . . . . . . . . . . . 124
Calculation . . . . . . . . . . . . . . . 171 Demand Pattern ID . . . . . . . . . .15, 70
Calibration Model Data Requirements . 246 Demand Processing . . . . . . . . . . 119
Calibration Simulations . . . . . . . . . 248 Demand Scenarios . . . . . . . . . . . 124
Category . . . . . . . . . . . . . . 71, 84 Demand Statistics . . . . . . . . . . . .71
Category Type . . . . . . . . . . . . . 72 Design . . . . . . . . . . . . . . . . . 208
CHeck Valve . . . . . . . . . . . . . . 22 Design Fire Flow . . . . . . . . . . . . 137
Check Valves . . . . . . . . . . . . . . 38 Design Fire Pressure . . . . . . . . . . 136
CHemical Concentrations . . . . . . . . 90 Design Flow . . . . . . . . . . . . . . .27
Closed . . . . . . . . . . . . . 23, 29, 37 Design Head . . . . . . . . . . . . . . .27
Developing Pipe Demand Coefficients . 122
Diagrams . . . . . . . . . . . . . . . 175 G
Diameter . . . . . . . . . . . 20, 37, 43, 47 General Purpose Valve . . . . 36, 115, 221
Distributed Demands . . . . . . . 59, 119 Genetic Algorithms Calibration . . . . . 143
Diurnal Curves and Patterns . . . . . . 88 Geocoding . . . . . . . . . . . . . . . 66
Diurnal Pattern . . . . . . . . . . . . . 86 Getting Started . . . . . . . . . . . . . 111
Diurnal Profiles . . . . . . . . . . . . . 83 Global Bulk Reaction Rate Coefficient . 95
Dual-Acting Valve CHaracteistics Curve 17 Global Editing . . . . . . . . . . . . . . 125
Dual-Acting Valve Characteristics . . . . 56 Global Pipe Wall Reaction Rate Coefficient
95
E GPV . . . . . . . . . . . . . . . . . . . 221
Effect of Pressure Reducing Valves . . 129 Graph . . . . . . . . . . . . . . . . . . 85
Efficiency . . . . . . . . . . . . . . . . 53 Graphical Input . . . . . . . . . . . . . 117
Efficiency Curve . . . . . . . . . . . . . 30 Graphical Placement 17, 23, 33, 38, 44, 47
Elevation . . . . . . . . . . . . . . . . 14 Gravity . . . . . . . . . . . . . . . . . 51
Emitter . . . . . . . . . . . . . . . 16, 114 Grid Points . . . . . . . . . . . . . . . 185
Emitter Exponent . . . . . . . . . . . . 50 Groundwater Wells . . . . . . . . . . . 223
Energy . . . . . . . . . . . . . . . . 53, 54 Group Definition . . . . . . . . . . . . . 142
Energy Editor . . . . . . . . . . . . . . 53
Energy Price . . . . . . . . . . . . . . . 30 H
Energy Price Pattern . . . . . . . . . . 30 Hazen-Williams . . . . . . . . . . . . . 139
Engineering Tables . . . . . . . . . . . 98 Head Loss Differential Method . . . . . 244
Entering Data . . . . . . . . . . . . . 116 Headlosses . . . . . . . . . . . . . . . 49
EPANET . . . . . . . . . . . . . . 108, 205 Hiding and Showing Graphical Data . . 119
EPANET Analysis Results . . . . . . 129 High End Flow . . . . . . . . . . . . . 28
EPANET Computer Program . . . . . 205 High End Head . . . . . . . . . . . . . 28
EPANET Input Data File . . . . . . . 205 History . . . . . . . . . . . . . . . . . 206
Equivalent Person EP . . . . . . . . . . 66 Hydraulic Simulation Model . . . . . . . 226
Equivalent Tenement ET . . . . . . . . 66 Hydraulic Structures . . . . . . . 180, 185
Estate Height . . . . . . . . . . . . . . 14 Hydraulic Time Step . . . . . . . . . . 87
Evolutionary Algorithms . . . . . . . . 139
Executing the Analysis . . . . . . . . 125 I
Extended Period Data . . . . . . . . . 115 Importing Graphical Data . . . . . . . . 118
Extended Period Hydraulics . . . . 49, 207 Inactive Volume . . . . . . . . . . . . . 43
Extended Period Simulations . . . . . . 73 Initial Conditions . . . . . . . . . . . . 182
Extended Period Water Quality . . . . . 49 Initial Level . . . . . . . . . . . . . 43, 47
Extended Rule-Based Controls . . 110, 157 Initial Simulation . . . . . . . . . . . . 248
Initial Water Quality . . . . . . . . . . . 91
F Input Data Requirements . . . . . . . . 113
FCV . . . . . . . . . . . . . . . . . . 221 Interactive Data Entry . . . . . . . . . . 116
FIFO . . . . . . . . . . . . . . . . . 41, 42 Introduction . . . . . . . . . . . . . . . 141
Fire Flow Analysis . . . . . . . . . 109, 135
Fire Flow Data . . . . . . . . . . . . . 116 J
Fixed HGL . . . . . . . . . . . . . . . . 42 Junction . . . . . . . . . . . . . . . . . 12
Flow Coefficient . . . . . . . . . . . . . 16 Junction Node Data . . . . . . . . . . . 114
Flow Control Valve . . . 36, 115, 145, 221 Junction Node Demands . . . . . . . . 185
Flow Direction . . . . . . . . . . 23, 33, 38 Junction Node Pressure Based Controls 77
Flow Paced Booster Source . . . . . . . 94
Forward Tracking . . . . . . . . . . . 100
263
Index
265
Index