Collection System
Collection System
Collection System
Modelling of storm water drainage networks and sewer
collection systems
User guide
MIKE 2017
2 MIKE URBAN - © DHI
PLEASE NOTE
LIMITED LIABILITY The liability of DHI is limited as specified in Section III of your 'DHI
Software Licence Agreement':
3
4 MIKE URBAN - © DHI
MIKE URBAN CS - MOUSE User Guide . . . . . . . . . . . . . . . . . . . . . . . . . . 13
1 Modelling Collection Systems . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
2 Modelling Collection Systems with MOUSE . . . . . . . . . . . . . . . . . . . 17
3 Hydraulic Network Modeling with MOUSE . . . . . . . . . . . . . . . . . . . . 19
3.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
3.2 Definition of a MOUSE Network . . . . . . . . . . . . . . . . . . . . . . . . . 19
3.2.1 Modelling real network elements . . . . . . . . . . . . . . . . . . . . 20
3.3 Nodes and Structures . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
3.3.1 Identification group . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
3.3.2 MOUSE model data group . . . . . . . . . . . . . . . . . . . . . . . 26
3.3.3 Q-H relations for nodes . . . . . . . . . . . . . . . . . . . . . . . . . 28
3.3.4 Outlet head loss . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
3.3.5 Model Concept of Soakaway . . . . . . . . . . . . . . . . . . . . . . 29
Soakaway tab . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
3.4 Pipes and Canals . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
3.4.1 Identification group . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
3.4.2 Geometrical properties . . . . . . . . . . . . . . . . . . . . . . . . . 38
3.4.3 Hydraulic friction losses . . . . . . . . . . . . . . . . . . . . . . . . . 39
3.4.4 Miscellaneous . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40
3.5 Weirs . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41
3.5.1 Identification and connectivity . . . . . . . . . . . . . . . . . . . . . 43
3.5.2 Model data . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 44
3.6 Orifices . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 44
3.6.1 Identification and connectivity . . . . . . . . . . . . . . . . . . . . . 47
3.6.2 Model data . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 47
3.6.3 Defining a gate or a weir in an orifice . . . . . . . . . . . . . . . . . . 48
3.7 Stormwater Inlets . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
3.7.1 Curb Inlet (Lintel) . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50
3.7.2 On-grade Capture . . . . . . . . . . . . . . . . . . . . . . . . . . . 52
3.7.3 Capacity curves . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
3.8 Pumps . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54
3.8.1 Pump types . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54
Constant flow pumps . . . . . . . . . . . . . . . . . . . . . . . . . . 55
Constant speed pumps . . . . . . . . . . . . . . . . . . . . . . . . . 55
Variable speed pumps . . . . . . . . . . . . . . . . . . . . . . . . . 56
3.8.2 Identification and connectivity . . . . . . . . . . . . . . . . . . . . . 61
3.8.3 Model data . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 62
3.9 Valves . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 62
3.10 CRS & Topography . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 66
3.11 Emptying Storage Nodes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 69
4 Rainfall-Runoff Modelling with MOUSE . . . . . . . . . . . . . . . . . . . . . 71
4.1 Terms and Concepts . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 71
4.1.1 MIKE URBAN Catchments . . . . . . . . . . . . . . . . . . . . . . . 72
4.1.2 Connecting Catchments to the Network . . . . . . . . . . . . . . . . 73
4.1.3 Hydrological Models . . . . . . . . . . . . . . . . . . . . . . . . . . 73
5
4.1.4 Creating Hydrological Models for a Catchment . . . . . . . . . . . . 74
4.2 Time-Area Method (A) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 75
4.2.1 Model Data . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 75
4.2.2 Parameter Sets . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 77
4.2.3 Time-Area Curve Editor . . . . . . . . . . . . . . . . . . . . . . . . 79
4.3 Kinematic Wave (B) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 80
4.3.1 Model Data . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 80
4.3.2 Parameter Sets . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 82
4.4 Linear Reservoir (C1 and C2) . . . . . . . . . . . . . . . . . . . . . . . . . . 84
4.4.1 Model Data . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 84
4.4.2 Parameter Sets . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 87
4.5 Unit Hydrograph Method (UHM) . . . . . . . . . . . . . . . . . . . . . . . . . 89
4.6 Additional Flow and RDI . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 92
4.6.1 Model Data . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 95
4.6.2 MOUSE RDI Parameter Sets . . . . . . . . . . . . . . . . . . . . . 96
4.7 Rainfall Data and Other Meteorological Variables - Boundary Conditions for Hydro-
logical models 99
4.8 Running the Runoff Computations . . . . . . . . . . . . . . . . . . . . . . . 102
4.9 MOUSE RDI - Guidelines for application . . . . . . . . . . . . . . . . . . . 104
4.9.1 Choice of calculation time step . . . . . . . . . . . . . . . . . . . . 104
4.9.2 The RDI hotstart . . . . . . . . . . . . . . . . . . . . . . . . . . . 105
4.9.3 The RDI result files . . . . . . . . . . . . . . . . . . . . . . . . . . 105
4.9.4 MOUSE RDI Validation . . . . . . . . . . . . . . . . . . . . . . . 106
Surface runoff model . . . . . . . . . . . . . . . . . . . . . . . . . . 107
General hydrological model - RDI . . . . . . . . . . . . . . . . . . . 107
4.9.5 Overflow within the model area . . . . . . . . . . . . . . . . . . . 111
4.9.6 Non-precipitation dependent flow components . . . . . . . . . . . 112
4.10 Using the Computed Runoff as Network Hydraulic Load . . . . . . . . . . . 112
4.11 Low Impact Development (LID) . . . . . . . . . . . . . . . . . . . . . . . . 114
4.11.1 LID Controls . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 115
Bio Retention Cells . . . . . . . . . . . . . . . . . . . . . . . . . . . 116
Infiltration Trenches . . . . . . . . . . . . . . . . . . . . . . . . . . 117
Porous Pavement . . . . . . . . . . . . . . . . . . . . . . . . . . . 118
Rain Barrels . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 118
Vegetative Swales . . . . . . . . . . . . . . . . . . . . . . . . . . . 119
Rain Garden . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 119
Green Roof . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 120
4.11.2 The LID Controls Editor . . . . . . . . . . . . . . . . . . . . . . . 121
Identification . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 125
LID Control data Specification . . . . . . . . . . . . . . . . . . . . . 126
4.11.3 LID Deployment . . . . . . . . . . . . . . . . . . . . . . . . . . . 130
Identification and Connectivity . . . . . . . . . . . . . . . . . . . . . 131
LID Deployment Properties . . . . . . . . . . . . . . . . . . . . . . 132
The LID deployment result file . . . . . . . . . . . . . . . . . . . . . 132
The LID Simulation Summary . . . . . . . . . . . . . . . . . . . . . 134
5 Time Series . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 137
5.1 Inserting New Time Series . . . . . . . . . . . . . . . . . . . . . . . . . . . 137
7
8 MOUSE Simulations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 189
8.1 The General Simulation Settings . . . . . . . . . . . . . . . . . . . . . . . 189
8.1.1 Choice of numerical engine . . . . . . . . . . . . . . . . . . . . . 192
8.1.2 Hotstart files . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 192
8.1.3 Export to remote simulation . . . . . . . . . . . . . . . . . . . . . 194
8.1.4 Batch simulations . . . . . . . . . . . . . . . . . . . . . . . . . . 195
8.2 The Runoff Simulation Settings . . . . . . . . . . . . . . . . . . . . . . . . 196
8.3 The Network Simulation Settings . . . . . . . . . . . . . . . . . . . . . . . 198
8.4 The Summary Simulation Settings . . . . . . . . . . . . . . . . . . . . . . 201
8.5 The 2D Overland Parameters . . . . . . . . . . . . . . . . . . . . . . . . . 204
8.6 MOUSE Result Selections . . . . . . . . . . . . . . . . . . . . . . . . . . . 204
9 2D Overland Flow . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 207
9.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 207
9.2 Input Required - Overview . . . . . . . . . . . . . . . . . . . . . . . . . . . 207
9.3 Input Required - Details . . . . . . . . . . . . . . . . . . . . . . . . . . . . 208
9.3.1 Choosing 2D model . . . . . . . . . . . . . . . . . . . . . . . . . 208
Flood Screening Tool in MIKE URBAN . . . . . . . . . . . . . . . . 209
9.3.2 Adding the DEM to the map view . . . . . . . . . . . . . . . . . . 210
9.3.3 Defining the 2D model domain and the resolution . . . . . . . . . . 210
9.3.4 2D model settings . . . . . . . . . . . . . . . . . . . . . . . . . . 211
9.3.5 2D initial conditions . . . . . . . . . . . . . . . . . . . . . . . . . . 213
9.3.6 2D initial condition polygon layer . . . . . . . . . . . . . . . . . . . 214
9.3.7 2D bed resistance polygon layer . . . . . . . . . . . . . . . . . . . 216
9.3.8 2D boundaries . . . . . . . . . . . . . . . . . . . . . . . . . . . . 218
9.3.9 Automatic model adjustments along water level and discharge boundaries
221
Adjustment of DEM . . . . . . . . . . . . . . . . . . . . . . . . . . . 221
Adjustment of initial conditions . . . . . . . . . . . . . . . . . . . . . 222
9.3.10 Defining couplings . . . . . . . . . . . . . . . . . . . . . . . . . . 222
9.3.11 Flow parameters at manholes and basins . . . . . . . . . . . . . . 226
Calculation method . . . . . . . . . . . . . . . . . . . . . . . . . . . 228
9.3.12 Outlets . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 230
9.3.13 Pumps and weirs . . . . . . . . . . . . . . . . . . . . . . . . . . . 230
9.4 Running the Combined 1D and 2D Simulations . . . . . . . . . . . . . . . . 230
9.4.1 Setting the simulation type and requesting 2D results . . . . . . . . 230
9.4.2 Starting simulation . . . . . . . . . . . . . . . . . . . . . . . . . . 234
9.5 Visualising Simulation Results . . . . . . . . . . . . . . . . . . . . . . . . . 234
9.6 Understanding the Files Generated . . . . . . . . . . . . . . . . . . . . . . 236
10 MOUSE Control module . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 239
10.1 RTC in Urban Drainage and Sewer Systems . . . . . . . . . . . . . . . . . 239
10.2 Architecture of RTC Systems . . . . . . . . . . . . . . . . . . . . . . . . . 240
10.3 MOUSE Control vs. Real Life . . . . . . . . . . . . . . . . . . . . . . . . . 242
10.4 Sensors . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 242
10.5 Logical Conditions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 243
10.6 Control Actions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 245
9
13.2 Surface Runoff Quality (SRQ) . . . . . . . . . . . . . . . . . . . . . . . . . 298
13.2.1 Surface Sediment Data Dialogs . . . . . . . . . . . . . . . . . . . 299
MOUSE| Water Quality| SRQ Main Parameters . . . . . . . . . . . . 299
MOUSE| Water Quality| SRQ Attached Pollutants . . . . . . . . . . . 300
MOUSE| Water Quality| SRQ Gully Pot Data . . . . . . . . . . . . . 301
13.3 Advection-Dispersion (AD) . . . . . . . . . . . . . . . . . . . . . . . . . . 302
13.3.1 Advection-Dispersion Data Dialogs . . . . . . . . . . . . . . . . . 303
MOUSE| Water Quality| AD Components . . . . . . . . . . . . . . . 303
MOUSE | Water Quality | AD Dispersion . . . . . . . . . . . . . . . . 305
Advection-Dispersion and Open Boundary Conditions . . . . . . . . 306
13.4 Biological Processes (BP) . . . . . . . . . . . . . . . . . . . . . . . . . . . 307
13.4.1 Biological Processes Dialog (MOUSE|Water Quality|WQ Process Model)
307
13.5 Water quality (MIKE ECO Lab) . . . . . . . . . . . . . . . . . . . . . . . . 311
13.5.1 MIKE ECO Lab Templates and Assignments . . . . . . . . . . . . 312
13.5.2 MIKE ECO Lab State Variables . . . . . . . . . . . . . . . . . . . 314
13.5.3 MIKE ECO Lab Forcings . . . . . . . . . . . . . . . . . . . . . . . 315
13.5.4 MIKE ECO Lab Constants . . . . . . . . . . . . . . . . . . . . . . 317
13.5.5 Running MIKE ECO Lab simulation . . . . . . . . . . . . . . . . . 318
13.6 Sediment Transport (ST) . . . . . . . . . . . . . . . . . . . . . . . . . . . 319
13.6.1 The Sediment Transport Models in MOUSE ST . . . . . . . . . . . 319
The Explicit Sediment Transport Models . . . . . . . . . . . . . . . . 319
The Morphological Models . . . . . . . . . . . . . . . . . . . . . . . 320
13.6.2 The Transport Formulae - Short Description . . . . . . . . . . . . . 320
The Ackers-White formulae . . . . . . . . . . . . . . . . . . . . . . 321
The Engelund-Hansen formula . . . . . . . . . . . . . . . . . . . . . 321
The Engelund-Fredsøe-Deigaard formulae . . . . . . . . . . . . . . 321
The van Rijn formulae . . . . . . . . . . . . . . . . . . . . . . . . . 321
13.6.3 The Flow Resistance in Sewer Systems with Sediment Deposits . . 322
13.6.4 Sediment Transport Data Dialogs . . . . . . . . . . . . . . . . . . 322
MOUSE | Water Quality | ST Main Parameters . . . . . . . . . . . . 322
MOUSE| Water Quality | ST Sediment Fractions . . . . . . . . . . . 324
MOUSE | Water Quality | ST Initial Sediment Depth Local . . . . . . 325
MOUSE| Water Quality| ST Sediment Removal Basins . . . . . . . . 326
MOUSE | Water Quality| ST Sediment Removal Weirs . . . . . . . . 327
13.6.5 Boundary Conditions for the Sediment Transport Model . . . . . . . 328
13.7 Storm Water Quality . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 329
13.7.1 Assigning concentrations of pollutants to runoff and infiltrations . . . 329
Cst. concentration (method 1) . . . . . . . . . . . . . . . . . . . . . 332
Table concentration (method 2) . . . . . . . . . . . . . . . . . . . . 332
EMC formula (method 3) . . . . . . . . . . . . . . . . . . . . . . . . 332
13.7.2 Local treatment nodes . . . . . . . . . . . . . . . . . . . . . . . . 333
11
12 MIKE URBAN - © DHI
MIKE URBAN CS - MOUSE
User Guide
13
14 MIKE URBAN - © DHI
1 Modelling Collection Systems
When modelling a collection system with MIKE URBAN you can choose to
model the collection system with either the SWMM5 engine or the MOUSE
engine.
15
Modelling Collection Systems
What are the return periods for overloading of various parts of the exist-
ing sewer system?
What are the main causes of that overloading - backwater or insufficient
local pipe capacity?
What are the implications of replacing critical sewers, installing new
basins, weirs, etc.?
How is the long-term environmental impact affected by changing the
operational strategy?
Where and why are sediments deposited in the sewer network?
What are the peak concentrations of pollutants at the overflow weir or at
the treatment plant after a rainstorm?
17
Modelling Collection Systems with MOUSE
3.1 Introduction
MOUSE allows for the hydrodynamic simulation of flows and water levels in
urban storm drainage and wastewater collection networks, thus providing an
accurate information about the network functionality under a variety of bound-
ary conditions. The hydrodynamic simulations can be extended with pollution,
sediment transport and water-quality simulations. The model can also be
enhanced by the variety of real-time control functions. The simulations can be
carried out for single events or as efficient long-term simulations for longer
historical periods.
19
Hydraulic Network Modeling with MOUSE
Copying network data from MIKE URBAN CS SWMM network into MIKE
URBAN CS MOUSE network
Graphically digitizing and manual data typing within MIKE URBAN
The last option is frequently used in a combination with one of the previous
options as means for achieving a full consistency of the MOUSE model.
The nodes will only have a single computational point where the water level H
is computed. The nodes are typically circular manholes in the sewer network.
But it can also be basins or tanks with a significant volume. Still only a single
water level computational point is located at the node. Based on the com-
puted water level and the description of the geometry of the node the compu-
tation keeps track of the volume of water stored in the node.
At the nodes the water level is computed based on the water level at the pre-
vious time step and the flow contributions during the time step from each con-
nected pipe and external connected flow like a catchment runoff discharge.
When the computational grid is set up for a network of links and nodes it will
end up like shown in Figure 3.3.
21
Hydraulic Network Modeling with MOUSE
The main point to realize is the conflict between computing a discharge for
these elements and the fact that only a water level is computed at nodes.
This is why the pump, weir, orifice, valve and storm water inlet elements from
the computational and numerical point of view are links and not an element
placed in one node. All the elements are links forming a connection between
two nodes.
In MIKE URBAN we have five functional elements which from the model
building point of view are related to nodes like manholes or basins. These are
"Pumps", Weirs", "Valves", "Orifices" and “Stormwater inlets”. The concept of
elements related to nodes is reflected in the design of the dialog for editing
the parameters for these elements. Here you find a field named "Location:"
for all of the elements. The field takes the ID of a node as input. All elements
also have a field for "To:" which also takes a node ID as input.
Seen from the computational solution point of view the five elements are
actually connections from one node to another node. This is similar to how
pipes are defining the link for flow between nodes as reflected in the dialog
where you find fields for entering "From node:" and "To node:".
The editor organizes the related input data into the following groups:
Each node is geographically determined by 'x' and 'y' co-ordinates. The co-
ordinates may be specified in any local co-ordinate system.
Manholes and basins are per default considered open at the top (Cover type
equal to 'Normal'). This means, that when the water level in a node reaches
the ground level, the water spills on the ground surface. In that case, MOUSE
introduces an artificial basin on the top of the node, with a surface area 1000x
larger than the node's surface. The surcharged water is stored in the basin, to
be returned back into the sewer.
23
Hydraulic Network Modeling with MOUSE
On the other hand, a node can be specified as a 'spilling' node (Cover type
equal to 'Spilling'). In a spilling node, water escapes irreversibly from the
model, if the water level reaches and exceeds the node's ground level
(optionally set off by a 'buffer pressure level). The rate of spill is approximated
as a free overflow over the crest at a given level and with a "conceptual" crest
length. For further details, see the MOUSE Pipe Flow Reference.
In the tables given below each data variable is described shortly and if it is
required as input.
25
Hydraulic Network Modeling with MOUSE
Table 3.2 The edit fields in the MOUSE model data group
Cover sub-group
Buffer pres- Buffer pressure is only active for type Yes BufferPres-
sure = spilling. Equal to the pressure above sure
the ground level needed to cause
spills from the manhole
Spill coeff. Spill coefficient is only active for type = Yes SpillCoeff
spilling. Controls the spill capacity
Table 3.2 The edit fields in the MOUSE model data group
Use local Switch for use of local head loss data Yes LossParNo
data
27
Hydraulic Network Modeling with MOUSE
Table 3.2 The edit fields in the MOUSE model data group
Specifying a Q-H relation for an outlet controls the flow at the outlet. When
specifying a Q-H relation for a manhole or basin the Q-H relation controls
infiltration to the node. The Q-H relation specifies the relation between the
water level in the manhole (or basin) and the infiltration flow. The flow (Q)
value in the Q-H relation should be given as a positive value when water
enters the node and a negative value for specifying a loss of water from the
network model.
If you wish to make changes to the head loss parameter set that you have
made a reference to, you can change this by accessing MOUSE|Local Head
Losses. The editor is shown in Figure 3.5.
29
Hydraulic Network Modeling with MOUSE
the upper vegetated layer. Then the stormwater infiltrates vertically through
the soakaway and infiltrates out of the sides and bottom of the soakaway.
In some cases the soakaway is not connected to any drainage network and
captured runoff to the soakaway is infiltrated and in case of extreme rainfall
and exceedance of its infiltration and storage capacity storm water is sur-
charged to the surface.
The soakaway can be added graphically and/or imported from any Asset GIS
system hence existing soakaways can be illustrated in the plan view of MIKE
URBAN.
The overflow pipe, weir or orifice is modelled in MIKE URBAN. Make sure
that the up level of the pipe, weir or orifice is set correctly, and not at the bot-
tom of the basin which will be the case if no invert level is specified. If the pipe
is placed at the bottom then make sure that flow regulation is applied to the
out-going pipe. The inflow pipe is generally not modelled; instead the catch-
ment that generates the runoff to the soakaway (e.g. a road section, a roof) is
connected directly to the basin via the catchment connection tools in MIKE
URBAN.
A new point feature layer has been created to illustrate the soakaway in the
plan view and to be used for graphical digitization and graphical connection of
soakaway to existing drainage network.
31
Hydraulic Network Modeling with MOUSE
At the 'Geometry' tab in the dialog the type of node is selected in the Node
type combobox and the following attributes are set:
Ground Level
Invert Level
Basin Geometry
Cover Type (Normal, Spilling and Sealed).
Soakaway tab
At the 'Soakaway' tab in the dialog the type of infiltration, porosity and initial
water level is selected
No Infiltration
Constant Infiltration
Infiltration
The option 'No Infiltration' is used in cases when the there is no infiltration out
of the soakaway. The initial water level can be set.
33
Hydraulic Network Modeling with MOUSE
dh 1
------ = ----------------- Q in t – Q f t – Q of t (3.1)
dt lw
Q f = K l w + 2h l + w (3.2)
Vol
Q f = K fs bottom A s h = 0 + K fs side 2A c + 2 --------- (3.3)
Ac
The infiltration from the bottom can be turned off by a flag. However the infil-
tration from both side and bottom can be shut off by setting the field-saturated
hydraulic conductivity to zero.
The porosity of the fill material is used to calculate the water level in the soak-
away and the initial water level is used to set the initial water level in the soak-
away. Table values of hydraulic conductivity, Kfs, for different soil classes are
provided in Table 3.3. Within each soil type the hydraulic conductivity varies
35
Hydraulic Network Modeling with MOUSE
The respective node bottom levels are displayed in the grey areas of the
“UpLevel” and “DwLevel” fields by selecting “Recompute” when clicking the
“Advanced” button.
In case of a step-wise connection (but not allowed below node bottom level),
the elevations of both the upstream and downstream connection must be
specified in the editable “UpLevel” and “DwLevel” fields.
Depending on the selected type, a link may take the form of one of the 'stand-
ard' pipes (Circular, Rectangular, O Shaped, Egg-Shaped), or any closed or
open cross section shape (CRS) and Natural Channels. The CRS and Natu-
ral Channels are defined in the CRS and Topography Editors.
Standard pipes are defined by diameter (or cross section width and height for
non-circular pipes), the geometry of special cross sections is as mentioned
specified under the cross section editor. In this dialog, only the reference to
the CRS ID.
The length of a link is calculated from the shape of the line in MIKE URBAN.
The length is displayed in the 'Length_C' field, but is not updated until the
'Recompute' command is executing (from the ‘Advanced’ button). If a user
defined length is specified this will overwrite the calculated one during simula-
tion.
37
Hydraulic Network Modeling with MOUSE
Table 3.5
Table 3.5
Edt field Description Used or Field name in
required by datastructure
simulations
Width Width of rectangular shape Yes, if Shape Width
= Rectangu-
lar
Height Height of rectangular shape Yes, if Shape Height
= Rectangu-
lar
CRS ID ID of cross section Yes, if shape CrsID
= CRS
Topography ID of topography Yes, if Shape Topogra-
= natural phyID
channel
Max Dx Max distance between gridpoints Yes, if Shape Maxdx
= natural
channel
Length Length of link Yes Length
UpLevel Upstream invert level of link Yes UpLevel
DwLevel Downstream invert level of link Yes DwLevel
Table 3.6
Edit field Description Used or Field name in
required by datastructure
simulations
Material Material of link Yes MaterialID
Formula- Formula for calculation of the friction Yes FricTypeNo
tion loss (Manning Explicit, Manning
Implicit, Colebrook White, Hazen-Wil-
liams)
Use local Determines if roughness values from Yes FricNo
data the material are overwritten by local
values
Manning Manning roughness value Yes, if ‘Man- Manning
ning Explicit’
or ‘Manning
Implicit’ is
chosen
39
Hydraulic Network Modeling with MOUSE
Table 3.6
3.4.4 Miscellaneous
The ‘Additional’ button provides access to more advanced options for a link
such as specifying a depth-variable Manning number.
3.5 Weirs
A weir is actually a functional relation, which connects two nodes of a
MOUSE network (two-directional flow and submerged flow possible), or is
associated with only one node (free flow 'out of the system'). The latter case
is achieved if the 'To' field is left empty.
It is possible to define several weirs between the same two nodes if this is
required. This is similar to the possibility of having more than one pipe as the
link between nodes. The generation of the computational grid shown in
Figure 3.15 for the orifice is also applied for pumps, weirs and valves. The
numerical solution of the flow equations will depend on the selected device.
Please refer to the reference manual on more on this.
It is recommended not to place the two nodes in the same spot, instead place
the nodes a short distance apart. The reason is that the node head loss com-
putation will have a component from change of flow direction. If the two nodes
surrounding the device are placed exactly at the same location then the com-
putational engine cannot determine the direction of the flow from the coordi-
nates of the nodes and a default direction will be applied. This may
unintentionally introduce a change in direction and therefore also an unex-
pected head loss.
By using a small displacement of the nodes the change in flow direction will
be determined based on the coordinates and angles between the connected
pipes. Therefore consider carefully the placement of the nodes with respect
to the actual construction.
41
Hydraulic Network Modeling with MOUSE
Orientation ('degrees') plays an important role (as long as the head loss coef-
ficient is undefined), since depending on the specified orientation, kinetic
energy of the flow is included (90o) or is not included (0o) in calculations of the
weir flows.
Weirs are per default static (No Control) but can be controlled through Real
Time Control (RTC). Clicking on the “RTC” button to the right gives quick
access to the RTC specification dialog.
Table 3.7
43
Hydraulic Network Modeling with MOUSE
Table 3.8
Con-
Oper. mode No control or RTC controllable weir Yes
trolTypeNo
Yes, if dis-
Weir orientation relative to the main
charge coeff.
Orientation flow direction. “0” is Side weir, “90” is a AngleNo
is not speci-
transversal weir
fied
Yes, if weir
Discharge
Discharge coefficient formula is Coeff
coeff.
chosen
Yes, if weir
Crest width Width of rectangular weir formula is CrestWidth
chosen
Yes, if Q-H is
Q-H table Reference to tabulated Q-H funtion QHID
chosen
Yes, if frag-
Source The ID of the source or the upstream
mented is SourceLinkID
Channel channel of the weir
chosen
3.6 Orifices
An orifice is actually a functional relation, which connects two nodes of a
MOUSE network or is associated with only one node (free flow 'out of the
system'). The latter case is achieved if the 'To' field is left empty.
In the real world a flow restriction in the form of an orifice may be located in a
manhole or a similar construction which you normally would define as a node
in the model configuration. The numerical solutions for the flow equations,
however, need a model configuration with two nodes where the orifice is
defined as the connection between the nodes. The orifice will then be placed
between the two nodes as the flow connection.
Figure 3.15 The difference between real world orifice and model configuration of ori-
fice
It is possible to define several orifices between the same two nodes if this is
required. This is similar to the possibility of having more than one pipe as the
link between nodes. The generation of the computational grid shown in
Figure 3.15 for the orifice is also applied for pumps, weirs and valves. The
numerical solution of the flow equations will depend on the selected device.
Please refer to the reference manual on more on this.
It is recommended not to place the two nodes in the same spot, instead place
the nodes a short distance apart. The reason is that the node head loss com-
putation will have a component from change of flow direction. If the two nodes
surrounding the device are placed exactly at the same location then the com-
putational engine cannot determine the direction of the flow from the coordi-
nates of the nodes and a default direction will be applied. This may
unintentionally introduce a change in direction and therefore also an unex-
pected head loss.
By using a small displacement of the nodes the change in flow direction will
be determined based on the coordinates and angles between the connected
45
Hydraulic Network Modeling with MOUSE
pipes. Therefore consider carefully the placement of the nodes with respect
to the actual construction.
A discharge coefficient can be specified (default = 1.0) and a flap gate (or
non-return valve) can be specified.
Orifices are per default static (No Control) but an orifice can be controlled
through Real Time Control (RTC). Clicking on the “RTC” button to the right
gives quick access to the RTC specification dialog.
Table 3.9
Table 3.10
Edit field Description Used or Field name in
required by datastructure
simulations
Type Type of orifice according to shape Yes TypeNo
Flap Flap indicating a flap-gate built-in (i.e. Yes FlapNo
no return flow possible)
47
Hydraulic Network Modeling with MOUSE
Table 3.10
Edit field Description Used or Field name in
required by datastructure
simulations
Oper. mode Yes Con-
No control or RTC controllable orifice
trolTypeNo
Discharge Calibration coefficient. Value = 1 Yes Dis-
coeff results in the flow as determined by chargeCoeff
orifice algorithm
Invert level Absolute elevation of the orifice invert Yes InvertLevel
Height Height of rectangular orifice Yes Height
Width Width of rectangular orifice Yes Width
Diameter Diameter of circular orifice Yes Diameter
CRS ID Reference of a cross-section ID for Yes, if CRS CrsID
irregularly-shaped orifice type is cho-
sen
The orifice itself is just an opening with a static shape. In real constructions
orifices are often equipped with a controlled gate or weir which can be used in
real time control for regulating the flow through the orifice. The gate device
will move from the top of the orifice opening and downwards until the orifice is
fully closed. The weir moves from the bottom of the orifice upwards and
closes fully when the weir crest reaches the top of the orifice opening (see
Figure 3.17 for an illustration). It is possible to apply both types of movable
devices in the computations. In both cases the device is "added" to a defined
orifice. This is done from the Controllable Devices dialog (MOUSE|Con-
trol|Controllable Devices). See more on this in section 9.8.
A typical Curb Inlet/grate configuration is shown below. Flow into the pit
chamber is via both a grate and side weir (operates as an orifice for deeper
flow depths).
49
Hydraulic Network Modeling with MOUSE
SAG Type, where the connection node on the overland flow network is
located at a sag or low point where water will collect. Transfer capacity of
the connection is specified as a DQ-relation (tabular data type).
ON-GRADE Type, where flow in the overland flow network can continue
past the connection node. Transfer capacity of the connection is depend-
ent on the slope of the overland flow network, and specified as a Capture
ID (collective of QQ-relations defining the capture rate as a proportion of
approach flow).
51
Hydraulic Network Modeling with MOUSE
Invert level (m) defining the point at which spilling starts (similar to weir
crest level). The user is shown a system calculated invert level which is
the same as the invert of the connection node in the overland flow net-
work. As with weir flow, a crest level at least 0.01 m higher than the con-
nection node invert level is recommended for initial condition stability.
Freeboard (m), defining a critical water level (Invert - Freeboard) at the
connection node in the pipe network below which the defined DQ and
QQ-relations apply. For submerged and reverse flow (surcharge), the
transfer capacity of the connection reverts to a standard orifice relation-
ship.
Slope (%), representing the slope of the steepest link in the overland flow
network entering the connection node (only applies to ON-GRADE
Type). The system calculated slope is used in the calculation unless a
user defined slope is specified.
Blockage factor (%) which can be used to account for debris blockage at
the grate/inlet. This linear factor is applied to the tabular data sets defin-
ing the transfer capacity of the connection.
Number of Curb Inlets, allowing multiple curb inlets of the same specified
geometry (transfer capacity) applied at the same location within a single
connection.
Default rectangular orifice geometry, applies to those flow cases (sub-
merged and reverse flow) were the defined DQ and QQ-relations do not
apply. This generally applies when water levels at the connection node in
the pipe network exceed the critical level defined by the Freeboard,
including reverse flow (surcharge).
There are no limitations on the number of curb inlets specified at one location;
however, the connectivity must be ‘From’ a node in the overland flow network
‘To’ a node in the pipe network, for correct automatic calculation of slope.
Note: Link slopes must be calculated in the link dialog for automatic calcula-
tion of slope to operate.
The On-grade Capture dialog allows the user to group together QQ-relations
(tabular data) that comprise a single On-grade Curb Inlet geometry (similar in
function to the Topography dialog). As the transfer capacity for an On-grade
Curb Inlet is dependent on the slope in the overland flow network, a number
of QQ-relations can apply.
For calculated or user defined slopes in the Curb Inlet dialog that are outside
the range of slopes specified in the On-grade Capture dialog, the closest
slope curve will be used. For intermediate calculated or user defined slopes
(lying between slope curves in the On-grade Capture dialog), linear interpola-
tion is applied.
In the case of an On-grade Curb Inlet capacity that is not dependent on slope
of the overland flow network, the user needs to define the On-grade Capture
with a single QQ-relation. Note: In this case, the calculated or user defined
slope in the Curb Inlet dialog for ON-GRADE Type will be ignored.
Two curve types specified in the tabular data (MOUSE|Curves & Relations)
can be used with the two different types of Curb Inlets.
The DQ relation specifies the depth based capacity curve for a SAG Type
Curb Inlet. Values must be monotonously increasing in depth and discharge
and starting at (0,0). For depths in excess of the maximum value specified in
the last row of the table, the last corresponding discharge value is used.
The QQ relation specifies the relationship between approach flow in the over-
land flow network (Qapp) and the captured flow at the connection node for an
ON-GRADE Type Curb Inlet (Qcap). Values must be monotonously increasing
and starting at (0,0). For approach discharges in excess of the maximum
53
Hydraulic Network Modeling with MOUSE
value specified in the last row of the table, the last corresponding capture dis-
charge value is used.
3.8 Pumps
A pump is actually a functional relation, which connects two nodes of a
MOUSE network or is associated with only one node (free flow 'out of the
system'). The latter case is achieved if the 'To' field is left empty. In this case it
is also possible to let the water pumped out of the system enter a 2D model
(in the case of coupling MIKE URBAN to MIKE 21, i.e. MIKE FLOOD). This is
done by setting the tickmark in ‘Coupling to 2D overland flow’.
It is possible to define several pumps between the same two nodes if this is
required. This is similar to the possibility of having more than one pipe as the
link between nodes. The generation of the computational grid shown in
Figure 3.15 for the orifice is also applied for pumps, weirs and valves. The
numerical solution of the flow equations will depend on the selected device.
Please refer to the reference manual on more on this.
It is recommended not to place the two nodes in the same spot, instead place
the nodes a short distance apart. The reason is that the node head loss com-
putation will have a component from change of flow direction. If the two nodes
surrounding the device are placed exactly at the same location then the com-
putational engine cannot determine the direction of the flow from the coordi-
nates of the nodes and a default direction will be applied. This may
unintentionally introduce a change in direction and therefore also an unex-
pected head loss.
By using a small displacement of the nodes the change in flow direction will
be determined based on the coordinates and angles between the connected
pipes. Therefore consider carefully the placement of the nodes with respect
to the actual construction. For pumps the distance between the nodes will
often be larger.
When specifying a constant flow pump, set the Speed to ‘Constant’ and
choose ‘Constant’ in the Cap. curve type field and type in the constant flow in
the Constant flow field.
Notice: 'Constant flow' pump is different from 'Constant speed' pump. The
later may have varying discharge.
Pump curves are defined in MIKE URBAN by tabulating the curve (see more
in section 5). This can be done by selecting representative points on the
pump curve as shown in the next figure. Information about pump curves are
typically provide by the manufacturer of the pumps.
55
Hydraulic Network Modeling with MOUSE
When a pump curve is tabulated like this and used in the simulations exe-
cuted with the MOUSE Engine, then the pump is allowed to operate within
the range of the dH values in the table. At any time during the simulation the
MOUSE Engine will compute the actual dH and determine the corresponding
Q value from the table, resulting in the actual duty point position for the pump.
As a standard feature the MOUSE Engine will stop the simulation with an
error message if the hydrodynamic conditions result in an actual dH value
outside the range of the pump curve table.
It is also found that variable speed pumps cannot be regulated over the full
range between 0% and 100%. Instead the regulation may be in the range of
70% to 100% as indicated in the figure above. Search for specific information
available for the pumps applied at the pumping stations being modeled.
For simplicity in the modeling the actual method of varying the pump capacity
is not considered. It is not directly related to the actual number of rotations
per minute (RPM) or the actual electrical power input. In the modeling we
only consider the pump capacity as varying between a pump curve corre-
sponding to the minimum speed and a pump curve corresponding to the max-
imum speed.
For variable speed pumps the two pump curves are defined as 'RPMmin' and
'RPMmax' pump curves. For constant speed pumps only the 'RPMmax' pump
curve is applied.
When the MOUSE simulation is executed it is the standard condition that the
simulation will stop and give an error message if the conditions exceed the
range for the dH value in the pump curve table. This applies both for the
'RPMmax' and the 'RPMmin' pump curve.
When a variable speed pump is in operation during the simulation with the
MOUSE Engine the duty point will move between the two pump curves
depending on the actual conditions. The actual flow is determined by using a
PID function which will attempt to control the water level at the wet well to
stay at the specified set point value.
The discharge Q determined by the PID function and the dH found by the
actual hydraulic conditions defines the duty point. This may vary in the area
between the two pump curves.
57
Hydraulic Network Modeling with MOUSE
Figure 3.27
Eventually the conditions may result in the duty point ending up on one of the
pump curves. If the PID regulation sets the discharge Q to a value higher
than the limitation by the 'RPMmax' pump curve at the given dH condition,
then the discharge will be defined by the pump curve. As a consequence the
water level in the wet well will rise above the defined set-point value.
Figure 3.28
If the inflow to the pumping station is low, then the operation of the pump may
result in an actual duty point located at the 'RPMmin' pump curve. In this case
the water level at the wet well will drop below the defined set-point value.
Eventually the water level will reach the stop level defined for the pump and
the operation is switched off
Figure 3.29
More detailed output information may be added to the standard result file for
variable speed pumps. Please refer to the documentation on the DHIAPP.INI
file found in the “DHIAPP INI and ADP Files Reference Manual”.
There can be various reasons for situations where pumps are running dry.
When this happens in the model simulation it indicates a mismatch in the
configuration. E.g. too small wet well volume, errors in the applied pump
capacity curves or deceleration times.
Pumps will be switched off when the water level drops to the defined stop
level. But the pump flow continues during the specified deceleration period by
a linear decrease to zero flow. During this time interval the water level in wet
well may eventually drop to the bottom of the wet well during a single time
step. In this situation the MOUSE Engine will force an immediate 'dry stop' of
the pump discharge.
59
Hydraulic Network Modeling with MOUSE
A pump is characterised by the 'Start Level' and 'Stop Level', an offset, accel-
eration and deceleration time and a capacity curve. The capacity curve is
specified in the Curves and Relations dialog. The capacity curve can be
specified as a 'Capacity Curve QH' relation (for screw pumps) or as 'Capacity
Curve QdH' relation (for differential head pumps), where 'H' is the absolute
water level in the pump's wet well (at 'Location'), and 'dH' is the water level
difference between the 'To' and the 'Location' nodes. A pump type with a
'Capacity Curve QH' relation is named a screw pump, while a pump type with
a 'dH-Q' relation is named a differential head pump.
Pumps are per default static (No Control) but can be controlled through Real
Time Control (RTC). Clicking on the “RTC” button to the right gives quick
access to the RTC specification dialog.
Table 3.11
61
Hydraulic Network Modeling with MOUSE
Table 3.12
3.9 Valves
A valve is a functional relation which connects two nodes of a MOUSE net-
work.
It is possible to define several valves between the same two nodes if this is
required. This is similar to the possibility of having more than one pipe as the
link between nodes. The generation of the computational grid shown in
Figure 3.15 for the orifice is also applied for pumps, weirs and valves. The
numerical solution of the flow equations will depend on the selected device.
Please refer to the reference manual on more on this.
It is recommended not to place the two nodes in the same spot, instead place
the nodes a short distance apart. The reason is that the node head loss com-
putation will have a component from change of flow direction. If the two nodes
surrounding the device are placed exactly at the same location then the com-
putational engine cannot determine the direction of the flow from the coordi-
nates of the nodes and a default direction will be applied. This may
unintentionally introduce a change in direction and therefore also an unex-
pected head loss.
By using a small displacement of the nodes the change in flow direction will
be determined based on the coordinates and angles between the connected
63
Hydraulic Network Modeling with MOUSE
pipes. Therefore consider carefully the placement of the nodes with respect
to the actual construction.
The valve is by default static, in which case the valve opening must be speci-
fied. It is also possible to define a valve to be the RTC controlled and then the
valve opening will be controlled by any specified control rules.
Table 3.13 Identification and connectivity edit fields of the MIKE URBAN valve edi-
tor.
Table 3.14 Geometrical and hydraulic properties, edit fields of the MIKE URBAN
valve editor
Table 3.15
65
Hydraulic Network Modeling with MOUSE
Table 3.15
CRS editor: Open and closed Cross sections. The X, Z types are appro-
priate for irregular cross sections, while H, W are best for symmetric
cross sections.
Topography: The "Channel Topography" dialog is an editor for MOUSE
topography data. A topography is made up from a number of cross sec-
tions,specified in the cross section editor.
The CRS and Topography Editors are reached by selecting MOUSE | CRS &
Topography.
Cross sections are classified in seven types: three of them are closed cross
sections, and four of them are open cross sections. Each of the types has
three sub-types, defined by the way how the CRS geometry is described.
Thus, the following CRS types are supported:
A series of lateral resistance factors can hence be specified for the cross sec-
tion. The markers (see Figure 3.34) are used for splitting the cross section
into parallel channels defined by the markers. The marker (1) and the marker
(3) define where the cross section is ‘cut off’, i.e. outside of these markers the
cross section will not be taken into account. If marker (1) and (3) are not
specified the entire cross section will be used. Marker (2) is not used for
MOUSE simulations. Marker (8) and (9) indicate the left and right channel
banks respectively. The markers are used for splitting the cross section into
three parallel channels. Thus the reach of the left flood plain is given by mark-
ers (1) and (8), the reach of the main channel is defined by markers (8) and
(9) and finally the right flood plain is defined through markers (9) and (3). The
calculation of the hydraulic parameters is carried out for each of the three
channels.
67
Hydraulic Network Modeling with MOUSE
Figure 3.34 The marker options for the cross section type X-Z-R-M.
The X, Z types are appropriate for irregular cross sections, while H, W are
best for symmetric cross sections.
For X,Z and H,W closed cross section types, MOUSE automatically provides
the Preissmann slot, in order to facilitate the flow computations in pressurised
conditions.
The 'Chainage' is the location for which the CRS is valid and 'Bottom Level' is
the elevation of the CRS at the specific chainage. MOUSE will interpolate
between each of the inserted CRS according to chainage. A CRS must
always be specified at chainage = 0.0. If the last CRS is specified at a chain-
age lower than the length of the pipe the last CRS will be extrapolated to the
remainder of the link.
Storage nodes are defined only by their name (ID-string), a receiving node, a
control node and Emptying function (QH-relation).
Storage nodes are not connected to the MOUSE network by links, which is
the case with all other types of nodes.
Water arrives into a storage node over a weir or a pump. This process is con-
trolled by the actual hydraulic situation in the system and the weirs or pumps
capacities. Weirs and pumps behave in this case as if they discharge "out of
the system". The volume of water that can be stored in a storage node is
unlimited.
Return of water from storage basins back into the MOUSE network is con-
trolled by the emptying functions specified in the dialog "Emptying Storage
Nodes". An emptying function is actually a Q-H function, where 'Q' is the flow
from the storage node into a 'receiving node', and 'H' is the water level in a
'control node'. The flow according to the emptying function is possible only
until the storage node is completely empty.
Table 3.16
69
Hydraulic Network Modeling with MOUSE
The user must be aware that runoff computation and its subsequent use as a
network load are, in principle two distinct steps in the modeling process.
71
Rainfall-Runoff Modelling with MOUSE
Catchments are essential for any hydrological model. In MIKE URBAN, the
geographical extent of a catchment is determined by the catchment polygon
perimeter. MIKE URBAN provides upon user request the total catchment
area.
In order to transfer the runoff generated on the catchment surfaces into the
collection network, the MIKE URBAN MOUSE model must include the infor-
mation about the connection of the catchment outlet to the collection network.
One or multiple catchments can be connected to one node.
Details about the principles and methods for catchment connections can be
found in the chapter 13.4.
Surface runoff model: These are the most common type in urban runoff
analysis. The common characteristic of all the models in this class is that
only surface runoff is computed. This implies discontinuous runoff hydro-
graphs where flow starts as a result of rainfall and ceased back to zero
again after the end of rainfall. As such, these models are suitable for rel-
atively densely urbanized catchments with dominant amount of runoff
generated on impervious surfaces, and for single-event analyses (e.g.
design rainfall of certain recurrence interval). These models fail to pro-
vide realistic results in dominantly rural catchments and for long-term
analyses involving multi-event rainfall series.
73
Rainfall-Runoff Modelling with MOUSE
MIKE URBAN MOUSE includes a series of surface runoff models and one
continuous hydrological model. The surface runoff models available are:
Any of these models can be used on its own, or in a combination of one sur-
face runoff and MOUSE RDI.
Each time a catchment is inserted a record is also inserted under the hydrol-
ogy tabs. Then the user can add hydrology data for each of the runoff meth-
ods available. See e.g. Figure 4.2.
75
Rainfall-Runoff Modelling with MOUSE
The Time-Area model uses several parameters. For practical reasons, these
parameters have been grouped in parameter sets, which, in turn, can be
associated with certain catchment. By these means the entire model setup
can be established with a very small amount of information, while still allow-
ing for full spatial variation of model parameters for individual catchments
through the application of local values.
New parameter sets can be inserted and values of individual parameters can
be edited in the "Parameters Time-Area" editor. This can be activated at
MOUSE|Runoff Models|Parameters Time-Area.
77
Rainfall-Runoff Modelling with MOUSE
MIKE URBAN comes with a default parameter set (-DEFAULT-). User can
insert any number of parameter sets and edit them according to the needs.
Table 4.2 Overview of the Time-Area Parameter Set Database Attributes (Table
‘msm_HParA’)
Table 4.2 Overview of the Time-Area Parameter Set Database Attributes (Table
‘msm_HParA’)
You can define any number of custom time-area curves. These can be
inserted and edited in the editor for the curves and relations
(MOUSE|Curves&Relations). Each Time-area table must start with a pair of
values (0,0) and must end with a pair of values representing the whole (per
default, MIKE URBAN maintains T-A curves in percent (%), i.e. the last pair of
values in the table must be (100,100).
79
Rainfall-Runoff Modelling with MOUSE
Figure 4.6 The Curves&Relations editor with one correct user defined time-area
curve
Table 4.3 Overview of the Kinematic Wave Database Attributes (Table ‘msm_H-
ModB)
81
Rainfall-Runoff Modelling with MOUSE
Table 4.3 Overview of the Kinematic Wave Database Attributes (Table ‘msm_H-
ModB)
The Kinematic Wave model uses relatively large number of parameters. For
practical reasons, these parameters have been grouped in parameter sets,
which, in turn, can be associated with certain catchment. By these means the
entire model setup can be established with a very small amount of informa-
tion, while still allowing for full spatial variation of model parameters for indi-
vidual catchments through the application of local values.
New parameter sets can be inserted and values of individual parameters can
be edited in the "Parameters Kinematic Wave" editor. This can be activated at
MOUSE|Runoff Models|Parameters Kinematic Wave.
MIKE URBAN comes with a default parameter set (-DEFAULT-). User can
insert any number of parameter sets and edit them according to the needs.
Figure 4.8 The Kinematic Wave Parameter Sets Editor (MOUSE|Runoff Mod-
els|Parameters Kinematic Wave)
Table 4.4 Overview of the Kinematic Wave Parameter Set Database Attributes
(Table ‘msm_HParB’
83
Rainfall-Runoff Modelling with MOUSE
Table 4.4 Overview of the Kinematic Wave Parameter Set Database Attributes
(Table ‘msm_HParB’
MIKE URBAN includes this model in two versions: C1 and C2. These are var-
iants of the same model, and are used as national standards in the Nether-
lands and in France, respectively.
Figure 4.9 The Linear Reservoir Model Editor (MOUSE|Runoff Models|Linear Res-
ervoir (C)), set up for the work with model C1 (Dutch)
85
Rainfall-Runoff Modelling with MOUSE
Figure 4.10 The Linear Reservoir Model Editor (MOUSE|Runoff Models|Linear Res-
ervoir (C)), set up for the work with model C2 (French)
Table 4.5 Overview of the Linear Reservoir Database Attributes (Table ‘msm_H-
ModC’)
Table 4.5 Overview of the Linear Reservoir Database Attributes (Table ‘msm_H-
ModC’)
The Linear Reservoir model uses a number of parameters. For practical rea-
sons, these parameters have been grouped in parameter sets, which, in turn,
can be associated with certain catchment. By these means the entire model
setup can be established with a very small amount of information, while still
87
Rainfall-Runoff Modelling with MOUSE
allowing for full spatial variation of model parameters for individual catch-
ments through the application of local values.
New parameter sets can be inserted and values of individual parameters can
be edited in the "Parameters Linear Reservoir" editor. This can be activated
at MOUSE|Runoff Models|Parameters Linear Reservoir.
MIKE URBAN comes with a default parameter set (-DEFAULT-). User can
insert any number of parameter sets and edit them according to the needs.
Figure 4.11 The Linear Reservoir Parameter Sets Editor (MOUSE|Runoff Mod-
els|Parameters Linear Reservoir) - Dutch version
Figure 4.12 The Linear Reservoir Parameter Sets Editor (MOUSE|Runoff Mod-
els|Parameters Linear Reservoir) - French version
Table 4.6 Overview of the Linear Reservoir Parameter Set Database Attributes
(Table ‘msm_HParC’)
89
Rainfall-Runoff Modelling with MOUSE
91
Rainfall-Runoff Modelling with MOUSE
Instead of performing hydrological load analysis of the sewer system only for
short periods of high intensity rainstorms, a continuous, long-term analysis is
applied to look at periods of both wet and dry weather, as well as inflows and
infiltration to the sewer network. This provides a more accurate picture of
actual loads on treatment plants and combined sewer overflows.
When studying the real flow conditions in sewer systems, flow peaks during
rain events are often found to exceed the values that can be attributed to the
contribution from participating impervious areas. This is a consequence of the
phenomenon, usually named Rainfall Induced Infiltration. This differs from the
Rainfall Induced Inflow by the fact that it does not depend only on the actual
From a hydrological point of view, parts of the infiltration behave in the same
way as the inflow. Therefore, classification of total hydrological loads to infil-
tration and inflow is not suitable for modeling approach. Rather, to describe
appropriately the constitutive components of flow hydrographs, distinguished
by their hydrological behaviour, the following concept is used instead:
Distinctive for the FRC component is that it is not influenced by the previous
hydrological situation, i.e. high or low soil moisture content. It occurs as a
direct consequence of a rainfall. The FRC component consists of the inflow to
the sewer system and the fast flow component of the infiltration, not depend-
ent on previous hydrological conditions.
When combined with any of the MOUSE surface runoff models, MOUSE RDI
provides a platform for accurate and reliable computation of urban runoff, lib-
erated from the limitations inherent to "standard" urban runoff modeling.
93
Rainfall-Runoff Modelling with MOUSE
Figure 4.14 Different catchment response under the same rainfall, due to different
soil moisture conditions at the beginning of the rainfall
The model data for additional catchment flow and RDI are stored in the data-
base table msm_HmodCRC.
Figure 4.16 The Additional Flow and RDI Model Editor (MOUSE|Runoff Mod-
els|Additional Flow and RDI
95
Rainfall-Runoff Modelling with MOUSE
Table 4.8 Overview of the Additional Flow and RDI Database Attributes (Table
‘msm_HModCRC’)
The MOUSE RDI model uses a large number of parameters. For practical
reasons, these parameters have been grouped in parameter sets, which, in
turn, can be associated with certain catchment. By these means the entire
model setup can be established with a very small amount of information.
New parameter sets can be inserted and values of individual parameters can
be edited in the "Parameters RDI" editor. This can be activated at
MOUSE|Runoff Models|Parameters RDI.
MIKE URBAN comes with a default RDI parameter set (-DEFAULT-). User
can insert any number of parameter sets and edit them according to the
needs.
Figure 4.17 The MOUSE RDI Parameter Sets Editor (MOUSE|Runoff Mod-
els|Parameters RDI)
Table 4.9 Overview of the RDI Parameter Set Database Attributes (Table
‘msm_HparRDII’)
97
Rainfall-Runoff Modelling with MOUSE
Table 4.9 Overview of the RDI Parameter Set Database Attributes (Table
‘msm_HparRDII’)
Table 4.9 Overview of the RDI Parameter Set Database Attributes (Table
‘msm_HparRDII’)
MOUSE RDI also uses air temperature (used for snowmelt process)
and potential evapo-transpiration data. Air temperatures are given
as time series of instantaneous values. Potential evapo-transpira-
tion is provided as time series of accumulated evapo-transpiration
(e.g. on monthly basis). Equally as rainfall, both air temperature
and potential evapo-transpiration can be specified as constant for
the entire simulation period.
MIKE URBAN MOUSE can utilize time series data stored in DHI pro-
prietary format DFS0. Please refer to Chapter 4 for importing,
reviewing and editing of time series data.
All three categories of meteorological variables act as boundary conditions
for the hydrological model and as such must be specified through the
99
Rainfall-Runoff Modelling with MOUSE
Some data sources (e.g. dfs0 time series) include the information
about the data source co-ordinates. By activating the data source
Figure 4.18 The "Catchment Loads…" editor with a rainfall boundary condition
applied for all catchments in the model. The co-ordinates specified
under the "Data Source Location" represent the rain gauge location
At this point, the declared catchment boundary is just a "placeholder" for the
actual rainfall data. These have to be specified in the second step.
The items can be created and edited in the "Boundary Items" editor. This can
be open by "Items" button on the "Catchment Loads…" editor (recom-
mended), or through the menu at MOUSE|Boundary Conditions|Items.
If opened from the "Boundary Loads…" editor, the "Boundary Items" editor
displays only the items associated with the current boundary. Initially, for a
newly declared boundary, the list of items will be empty. A new item is created
by the "Insert" button.
101
Rainfall-Runoff Modelling with MOUSE
Per default, the specified boundary will be applied throughout the entire dura-
tion of the simulation or during the period covered by the time series. Option-
ally, the user can limit the temporal validity of the specified file to any period
within the DFS0 file's duration, by specifying the temporal validity. For the rest
of the simulation period, the boundary conditions must be "extended" by
some other time series or constant value. This is achieved by inserting and
specifying another boundary item.
The computation is set-up and initiated from the MOUSE Computation dialog
(Simulation|Run MOUSE) - see Chapter 7 for full detail).
On the "General" TAB, the user must specify the essential information about
the simulation (see example in Figure 4.19), and the runoff-specific parame-
ters on the "Runoff parameters" TAB (see example in Figure 4.20).
Figure 4.19 The "MOUSE Computation dialog (Simulation|Run MOUSE), with gen-
eral computation parameters
103
Rainfall-Runoff Modelling with MOUSE
Figure 4.20 The "MOUSE Computation dialog (Simulation|Run MOUSE), with run-
off-specific (T-A Curve) computation parameters
When calculating with MOUSE RDI, time steps are given separately for the
Surface Runoff Model and for the rain dependent infiltration part.
The RDI calculation can often be performed with a relatively long time step
(several hours), while calculation with the Surface Runoff Model is typically
performed with a time step in order of value of several minutes.
The time step for Surface Runoff computations is primarily concerned about
the sufficient resolution of the runoff process in time.
Generally, the RDI simulation time step should be chosen in accordance with
the resolution of precipitation data, e.g. a time step of 24 hours could be suit-
able if only daily precipitation data is available. However, in case when pre-
cipitation data with high resolution of e.g. few minutes are available, the RDI
time step should be chosen in accordance with the response of the discharge
when raining. E.g., an RDI time step of 2 - 4 hours should be chosen, if the
time constant CKOF is given a value of 8 hours.
To minimize the calculation time as well as the size of the result files the RDI
calculations are performed according to the following principle:
The RDI simulation is carried out continuously for the whole period specified.
On the contrary, the Surface Runoff simulation is carried-out only when rain-
ing. Thus, the start time for the Surface Runoff calculation is set as the start
time for rain hydrograph. The Surface Runoff calculation continues until all
the surface runoff hydrographs are regressed.
There is a HOTSTART facility for MOUSE RDI, i.e. the initial conditions for
the various storages can be automatically taken from a former result file, at a
simulation start time.
The structure and contents of the result file used as a HOTSTART file
requires that the time series in the boundary connection start at least for the
maximum specified concentration time Tc earlier than the start time for the
HOTSTART is specified. This is required for the correct reconstruction of the
surface runoff component (FRC).
Two result files are generated by a MOUSE RDI calculation. These are:
*.CRF file, containing maximum five time series for each sub-catchment,
namely:
– discharge, calculated with the Surface Runoff Model (the FRC com-
ponent),
– discharge, calculated with the RDI model (the SRC component),
– total discharge,
– variation of water content in the surface storage for the Surface Run-
off Model,
– variation of water content in the snow storage for the Surface Runoff
Model.
The *.CRF file is used as input data for a MOUSE Hydrodynamic calculation.
105
Rainfall-Runoff Modelling with MOUSE
In the *.CRF file the time series are saved with two various intervals, the
shorter one for the periods when the Surface Runoff Model is used, and a
larger one in the remaining periods. In the *.NOF file the time series are
saved with the larger time interval, which is equal to the time step used for the
RDI calculation.
Some of the parameters in MOUSE RDI (here meaning both for the rain
dependent inflow and the infiltration part) are related to actual physical data.
However, the final choice of parameter values must be based upon a compar-
ison with historical measured discharges, since a number of the parameters
have an empirical character.
The available period of the measured discharge data and its resolution in
time are of major importance for the credibility of the obtained parameter val-
ues. Ideally, for a good accuracy, a 3-5 years long time series of measured
discharge data with daily values is required for the calibration of the RDI
parameters. Several months long time series with higher resolution, i.e. min-
utes or hours, depending on the size of the area, are needed for the calibra-
tion of the surface runoff model. Measured time series with shorter duration
are also useful, although not securing optimal parameter values. In such case
it is important that the time series represents different hydrological situations,
i.e. typical wet period or dry period.
MOUSE RDI calculates the total generated discharge from a catchment, i.e.
overflow within the sub-catchment will also be included in the calculated dis-
charge. Therefore, when comparing with measured peak flows and con-
trolling the water balance (total volume) this has to be taken into
consideration.
The choice of criteria is important since it may affect the final choice of
parameter values and by that the behaviour of the calibrated model. Numeri-
cal criteria are, however, limited and therefore a visual comparison between
the hydrographs is indispensable.
MOUSE supports visual comparison of the calculated time series with any
time series of the same type contained in the time series database. E.g.,
when validating the model, the calculated discharge can be plotted on the
same graph with the measured discharge and compared.
In the example related to the illustrations, overflow occurs within the model
area. MOUSE RDI can not describe this kind of processes, which compli-
cates the choice of validation criteria.
To separate the Afrc component (Surface Runoff Model) and the fast part of
the SRC component (Surface Runoff Component in RDI), measured dis-
charge data with fairly high resolution in time (hours) is required.
107
Rainfall-Runoff Modelling with MOUSE
The effects of changing each particular parameter are discussed below. Also,
the most suitable hydrological periods for calibrating certain parameters are
identified, which implies that a certain parameter affects the model behaviour
more during periods with specific hydrological conditions. Usually, effects will
also be obtained during other periods, why these should also be studied
when adjusting a parameter.
The first step in the RDI calibration is usually to adjust the water balance
in the system, i.e. the accuracy between the calculated and measured
total volume during the observed period. This is done by correcting the
proportion of area, Asrc. An increase of Asrc proportionally increases
every flow component at each time step.
The total volume generally also contains the runoff from impervious
areas (Surface Runoff Model) - see Chapter 4.
The measured flow peaks generally also contain the runoff from imper-
vious areas (Surface Runoff Model).
CKBF is adjusted against the response of the baseflow, i.e. the build-up
and regression of the baseflow. Adjustment against the built-up of base-
flow is done during and after wet periods with low evaporation. Adjust-
ment against regression is done during the start of dry periods with high
evaporation, preferably when baseflow is the only flow component.
An adjustment of CKBF does not influence the size of the discharged
volume studied for a longer period, but displaces the volumes in time.
CKOF is adjusted against the response, i.e. the shape of the peak flows.
This is done during periods with heavy rainfall, preferably after a wet
period.
The measured flow peaks generally also contain the runoff from imper-
vious areas (Surface Runoff Model).
The threshold values indicate at which relative water content in the root
zone, L/Lmax, overland flow, interflow and baseflow respectively will be
generated. Therefore, the threshold values can be estimated from the
time of filling the root zone storage when each flow component starts dis-
charging.
109
Rainfall-Runoff Modelling with MOUSE
The threshold values have no effect during periods when the root zone
storage is full, L = Lmax.
TG is adjusted during periods with heavy filling of the root zone stor-
age, preferably in combination with low potential evaporation and
preceded by a dry period. TG is therefore an important parameter for
adjusting the increase of the groundwater level in the beginning of wet
periods.
TOF is adjusted after a dry period at events with heavy filling of the root
zone storage. For example adjustment can be done for events where
even larger rainfall volumes does not generate overland flow.
TIF is adjusted after a dry period when filling of the root zone storage,
preferably in combination with low potential evaporation. However, TIF
is one of the less important parameters.
Since the variation of water contents in the surface and root zone stor-
age controls many of the other processes, they should be studied con-
tinuously throughout the calibration procedure. Figure 9 shows an
example of the variation of water content in the surface storage, root
zone storage and groundwater storage. The example comes from the
catchment of Rya treatment plant, Göteborg, Sweden. It appears that
the root zone storage is emptied only during the summer period,
because the evaporation during the rest of the year is almost non-
existent. Discharge from the groundwater storage exists continuously
all year around. Drawing of the surface storage is faster during sum-
mer period since the evaporation is high, and is therefore the dominat-
ing effect on the surface storage. During periods with low evaporation,
drawing of the surface storage is controlled by the given time constant
for interflow, CKIF.
The example also shows that filling of the root zone and groundwater
storage only occurs when the surface storage is completely filled, i.e.
when precipitation has filled-up the surface storage. A larger surface
storage, i.e. a larger Umax, will therefore imply that this happens more
rarely and at a smaller extent, allowing a larger part of the precipitation
to evaporate.
A smaller root zone storage, i.e. a smaller Lmax, would have led to an
increased relative variation in the storage. Furthermore, the actual
evaporation will decrease in case of smaller root zone storage,
because less water is available for the vegetation to draw water for
transpiration, mainly during summer period.
Monthly and yearly values for the different processes, e.g. precipitation
volume, real evaporation and total discharge, are written to an ASCII
file, NAMSTAT.TXT after every RDI calculation. It is recommended that
the content of this file is studied now and then during the calibration
procedure.
In those cases when overflow occurs in the studied model area, e.g.
when simulating the discharge to the treatment plant, this has to be
considered when calibrating the peak flows during rainfall. MOUSE
RDI calculates the total generated discharge in the catchment area
and is therefore not able to describe hydraulic processes like e.g. over-
flow ("loss of water"). Calibration of parameters affecting the volume in
the peak flows should therefore be performed for rain events, when
overflow is unlikely to occur. Model parameters affecting the response
111
Rainfall-Runoff Modelling with MOUSE
of the discharge, for rain events when overflow occur, can be calibrated
against the peak flows base or width.
The foul flow is preferably estimated through daily values from pro-
duced water volumes weighted with yearly charged water volumes.
This will however only give a rough estimation, why departure from this
methodology may be necessary, e.g. for areas where a large amount
of freshwater is used for irrigation.
Concerning the network connectivity of the storm runoff, the option "Defined
by Catchment Connections" must be selected. This selection means that the
model will utilize the catchment connection information from the runoff result
file (*.CRF) used. I.e. no explicit declaration of the runoff connection point is
needed.
Figure 4.21 The "Network Loads" editor with runoff defined as a load. "Defined by
catchment connections" implies that the MOUSE runoff result file
(*.CRF) will be used as a source of connection information
At this point, the declared network load is just a "placeholder" for the actual
data. These have to be specified in the second step.
113
Rainfall-Runoff Modelling with MOUSE
The items can be created and edited in the "Boundary Items" editor. This can
be open by "Items" button on the "Network Loads" editor (recommended), or
through the menu at MOUSE|Boundary Conditions|Items.
If opened from the "Network Loads" editor, the "Boundary Items" editor dis-
plays only the items associated with the current load. Initially, for a newly
declared load, the list of items will be empty.
Per default, the specified runoff result file will load the network throughout its
entire duration, corresponding to the runoff computation Start and End times.
Optionally, the user can limit the temporal validity of the specified file to any
period within the file's duration, by specifying the temporal validity. For the
rest of the simulation period, the network can be loaded by some other runoff
result file. This is achieved by inserting and specifying another flow item.
Many practices have been used to adhere to these principles such as bio-
retention facilities, rain gardens, vegetated rooftops, rainwater tanks, and per-
meable pavements. By implementing LID principles and practices, storm
water can be managed in a way that reduces the impact of impervious built-
up areas and supports the natural movement of water within an urban eco-
system or catchment.
The ability to assess the benefit of the LID practices installed in urban catch-
ments is required by Water Utilities and other stakeholders responsible for the
urban drainage. Essential part of this requirement is ability to model the vari-
ous LID practices at both hydrological screening level and at a detailed
hydraulic level in order to evaluate the effect of installing LID practices as
parts of a storm water drainage system.
115
Rainfall-Runoff Modelling with MOUSE
Bio-retention cells, infiltration trenches, and porous pavement systems can all
be supplied with optional underdrain systems in their gravel storage beds to
convey storage runoff off the site rather than infiltrate it all. They can also
have an impermeable floor of liner that prevents any infiltration into the native
soil from occurring. Infiltration trenches and porous pavement systems can
also be subject to a decrease in hydraulic conductivity over time due to clog-
ging.
While some LID practices can also provide important pollutant reduction, the
current LIDS implementation in MIKE URBAN can only be used to model the
LIDS’ quantitative hydrologic performance.
The different structural layers of this feature are presented in Figure 4.22.
There are one or several relevant hydrological processes associated with
each layer.
In the surface occur surface storage, surface infiltration, the collecting area
run-on, surface evaporation and overflow.
In the soil layer occur percolation (vertical water movement to the storage
layer), evapotranspiration (loss due to the plant root action) and storage in
the soil's voids.
In the storage layer, storm water provided by the percolation through the soil
layer is detained in the storage layer's voids. Infiltration (leakage) to the
native soil through the storage bottom, controlled by the characteristics of the
surrounding soil, restores the storage capacity.
Infiltration Trenches
Infiltration Trenches are storage cells filled with gravel that capture runoff
from upstream impervious areas. They provide storage capacity and the pos-
sibility for captured runoff to infiltrate the soil underneath.
The processes simulated for the infiltration trench are similar as for the bio-
retention cell, except for the missing soil layer: storm water from the surface
enters the storage directly, i.e. without detention in the soil layer.
117
Rainfall-Runoff Modelling with MOUSE
Porous Pavement
Porous Pavement systems are excavated areas filled with gravel and paved
over with a porous concrete or asphalt mix. Porous pavements are imple-
mented in places where surface has to provide a firm support for vehicle traf-
fic, such as streets and parking lots.
Normally, all rainfall will immediately pass through the pavement into the
gravel storage layer below it where it can infiltrate at natural rates into the
site´s native soil.
Rain Barrels
Rain Barrels are containers that collect roof runoff during storm events and
can either release or re-use the rain water during dry periods.
A rain barrel includes storage (with total porosity, i.e. empty volume) and
underdrain.
The processes simulated include water detention in the storage, overflow and
drainage (i.e. storage recovery). The action of the underdrain is timed: it can
be specified the time offset between the end of rain and the drain activation.
This feature allows for the simulation of rainwater re-use for garden watering,
which normally takes place several hours or days after the rain has stopped.
Vegetative Swales
Vegetative Swales are waterways or depressed areas with sloping sides cov-
ered with grass and other vegetation. They slow down the conveyance of col-
lected runoff and allow it more time to infiltrate to inherent soil beneath it. The
only relevant layers regarding vegetative swales are surface and surrounding
soil.
Rain Garden
Rain Garden, as a simplified form of bio-retention cell is a garden, which
takes advantage of rainfall and storm water runoff in its design, simultane-
ously contributing to the reduction of runoff volume and pollutants released
from the site. Typically, it is designed as a small garden with selected types of
119
Rainfall-Runoff Modelling with MOUSE
The different layers comprising rain gardens are presented in Figure 4.27.
These include surface and soil layer
In the surface occur surface storage, surface infiltration, the collecting area
run-on, surface evaporation and overflow.
In the soil layer occur storage in the soil's voids, evapotranspiration (loss due
to the plant root action) and infiltration (leakage) to the native surrounding soil
through the bottom. The infiltration is controlled by the characteristics of the
surrounding soil.
Green Roof
Green Roofs are structural components that reduce the negative effects of
urbanization on water quality and rainfall runoff by absorbing or detaining,
and by filtering runoff. They are built of a planted lightweight soil media and a
drainage layer.
Green roofs are equipped with a drainage layer named drainage mat, placed
beneath the soil media. The purpose of the drain mat is to conduct the sur-
plus water percolated through the soil layer from the roof to the drainage sys-
tem.
The MIKE URBAN LID Controls editor organizes the data input for the differ-
ent layers and functional elements of Low Impact Development structures.
The input data is organized into following group and tabs:
Tabs: Surface, Soil, Pavement, Storage, Drain and Drainage Mat. These
tabs represent different layers and functional elements defining the struc-
ture of the LID control. The properties of each layer and functional ele-
ment, which constitute the actual LID control type, are entered in these
tabs. Only the actually relevant data fields for the current LID control type
are activated and available for the data input.
Figure 4.29 through Figure 4.34 illustrate the six tabs containing the data for
various functional elements of the LID structure. Table 4.10 and Table 4.11
provide a full reference for the attributes of the msm_LIDControl database
table.
121
Rainfall-Runoff Modelling with MOUSE
Figure 4.29 The MIKE URBAN LID Control Editor - Surface Tab
Figure 4.30 The MIKE URBAN LID Control Editor - Soil Tab
Figure 4.31 The MIKE URBAN LID Control Editor - Pavement Tab
Figure 4.32 The MIKE URBAN LID Control Editor - Storage Tab
123
Rainfall-Runoff Modelling with MOUSE
Figure 4.33 The MIKE URBAN LID Control Editor - Drain Tab
Figure 4.34 The MIKE URBAN LID Control Editor - Drainage Mat Tab
Identification
125
Rainfall-Runoff Modelling with MOUSE
127
Rainfall-Runoff Modelling with MOUSE
In Table 4.12 are presented examples of Manning' values (M) for different
types of surfaces. Table 4.13 shows hydraulic conductivity properties for vari-
ous porous media.
129
Rainfall-Runoff Modelling with MOUSE
The LID Controls are assigned to a catchment by means of the LID deploy-
ment editor. In this editor, it is possible to specify the size (i.e. area) of the
individual control, the size of the catchment´s impervious area that is fun-
neled into it and the number of the current LID control units.
Furthermore, level of initial saturation and routing width for overland flow can
be specified.
131
Rainfall-Runoff Modelling with MOUSE
If the user does not activate this option, the only visible outputs from the run-
off simulation including LIDs are the changed (reduced) runoff hydrographs,
caused by infiltration loss and storage in the LID control, and the LID sum-
mary table in the simulation summary file.
The contents of the dfs0 file depends on the actual LID control type. In
Figure 4.36, an example of the dfs0 file for porous pavement is presented.
Thirteen columns contain time series for the processes occurring inside the
porous pavement.
NOTE: The flow inside the LID control and the drain flow are reported as
intensities based on the LID control area. In cases where collecting area
is bigger than the LID control area, the reported intensities will not be
comparable with rainfall and evapotranspiration intensities, which are
given as the model boundary conditions. In order to make the compari-
son possible, the reported flow intensities must be scaled down by the
ratio between the LID control area and the collecting area.
Figure 4.36 Example of a result report dfs0 file per deployment (a porous pavement
control)
133
Rainfall-Runoff Modelling with MOUSE
3. Rain (mm/h): input rain for the catchment containing the LID control unit.
This rain loads the LID area directly.
4. Evaporation (mm/h): this is a given parameter specified by the user by
means of boundary conditions valid for the catchment as well as for the
LID unit.
5. Infiltration (mm/h): infiltration from the LID unit to the surrounding
native soil.
6. Surface flow (mm/h): surface water flow. Represents overland flow for
vegetative swale; overflow for other LID controls
7. Drain flow (mm/h): water flow through the drain.
8. Surface depth (m): accumulated water on the surface expressed as the
water height on the surface. Maximum value is the specified surface stor-
age height.
9. Soil moisture (): presented as a fraction, relative to the total volume of
the soil layer. Corresponds to the water that is held in the spaces
between soil particles. Maximum soil moisture (i.e. full saturation) is
equivalent to the specified soil porosity.
10. Pavement moisture(): presented as a fraction, relative to the total vol-
ume of the porous pavement layer. Corresponds to the water that is held
in the pavement's voids. Maximum pavement moisture (i.e. full satura-
tion) is equivalent to the specified pavement porosity.
11. Storage depth (m): accumulated water in the storage layer expressed
as the water height in the storage. Maximum value is the storage height.
12. Surface to soil (mm/h): infiltration from the surface layer to the soil layer
expressed in mm/h.
13. Surface to pavement (mm/h): infiltration from the surface layer to the
pavement layer expressed in mm/h.
14. Soil to storage (mm/h): ): flow from the soil layer to the storage layer
expressed in mm/h.
15. Pavement to storage (mm/h): flow from the pavement layer to the stor-
age layer mm/h.
16. MassChecksum (m3): regards to the error check in the mass water bal-
ance per time step computed, a low value in the order of 1*E-015 indi-
cates an accurate internal flow estimation.
135
Rainfall-Runoff Modelling with MOUSE
4 Time Series
Time series can represent both input and result time series. Examples of
input time series are e.g. rain event, a varying water level at an outlet etc.
These input time series can be entered in MIKE URBAN through the time
series editor. Time series in the time series editor are saved in dfs0 format. In
the time series editor it is possible so save several time series with the same
time axis in different time series items, these various items are then saved in
one time series object. When using the time series as tidal variations of e.g. a
specific load in the simulation, reference to the individual time series items
within the time series objects are made.
If the data are to be typed in anyway, the TS data input starts by specifying
the basic TS properties: number of values, start date and time, time interval
between successive values (assuming a uniform time resolution of the TS)
etc. When data are imported it may also be desirable to change some of the
properties, please refer to the example later in this section.
71
Time Series
After specifying the properties for the TS object, you can now continue to
specify the properties for the TS items. One item is per default already availa-
ble.
The properties of the time series object can be set in the dialog seen in
Figure 4.3. A description of the input fields is found in Table 4.1.
Field Description
Title Title of the time series object
Auto Conversion of When tick mark is set then the values that are in ‘Delete
Delete Values Value’ will be converted on the fly if the unit specified in
the ‘EUM Unit’ field is changed. Only conversion
between units of the same type can be done. If the tick
mark is not set then it is possible to change the unit with-
out the delete value typed in being converted.
Delete Value Threshold for values that are so small that they should
be disregarded.
Nr Days in Year If a time step of months is used, the number of days is
used to compute the number of days in a month. Other-
wise not used
TS Data Type Not relevant
73
Time Series
The properties of the time series item can be set in the dialog seen in
Figure 4.4. A description of the input fields is found in Table 4.2.
75
Time Series
Field Description
Item Name Name of the item
Data Type Type_Float: 4 byte floating point number,
Type_Double: 8 byte floating point number ,
Type_Char: 1 byte signed integer char, Type_Int: 4
byte signed integer int, Type_Unsigned: 4 byte
unsigned integer, Type_Short: 2 byte signed
integer short Type_Unsigned_Short: 2 byte
unsigned integer
Field Description
Value Type Instantaneous: means that the values are repre-
sentative at one precise instant.
Reverse_Mean_Step_Accumulated: is equal to
Mean_Step_Accumulated type, but values repre-
sent the timespan between the current Timestep
and the next Timestep. Used for forecasting pur-
poses.
77
Time Series
Properties
Allows to change properties on the x-axis, y-axis and text annotations, e.g.
display and placement of legends.
Graphic Settings
Allows to change the display of the graphs, e.g. line thickness, line colour,
line type (solid, dotted etc.).
Font
Allows to change the font used for the title, legends and the axis.
Zoom in, Zoom Out, Zoom to Selected Item, Next Zoom, Previous
Zoom
Allows to access various zoom functionality.
Draw Time Axis Label, Draw Left Axis Label, Draw Right Axis Label,
Draw Left Legend Label, Draw Right Legend Label
Controls the settings of the axis and legends related to the axis.
79
Time Series
Mode
Delete Points
When this mode is selected, when the mouse pointer is near a data value, the
pointer becomes a rubber and clicking on the left button of the mouse,
deletes the data value. The data value is not really deleted from the data set
but its value is set to delete value.
Move Points
When the mouse pointer is near a data value, the mouse pointer becomes a
moving cross and you can move the data value by moving the mouse while
keeping the left button of the mouse down (drag). Before the data value drag
begins, the data value is selected. When moving data values, if the selected
item is also being shown in a Timeseries Tabular, the correspondent cell
value in Timeseries Tabular changes accordingly.
Select Points
When clicking in a data value that data value is selected. The current Time-
step is set to the selected data value timestep. A red square around the data
value appears and the correspondent cell in the Timeseries Tabular (if the
selected item is also being shown in a Timeseries Tabular) is selected
Set Bookmarks
When the mouse pointer is near a data value not bookmarked, the pointer
becomes a pushpin. Clicking on the left button of the mouse bookmarks the
data value.
Clear Bookmarks
When the mouse pointer is near a bookmarked data value, the pointer
becomes a rubber. Clicking on the left button of the mouse removes the book-
mark for the data value
81
Time Series
profile plot, e.g. if you place the line at a specific time, the horizontal plot (if it
displays results layers) will display the results for the same time.
1. Go to Edit|Time Series
2. Right-click on the ‘Untitled’ and choose TSObject Properties, see
Figure 4.11
Figure 4.11 Right-Click to set the properties for the time series
3. On the dialog that comes up your are able to specify the various proper-
ties. The properties for this specific example are seen in Figure 4.12.
In this example we are specifying that the rain starts at the 1st of November
2004, that we would like 20 time steps that are a minute apart.
Figure 4.13 Example of how it looks after specifying the TSObject Properties
83
Time Series
Now right-click on the TS item (per default named ‘Undefined’ before renam-
ing) and specify the properties for the TS item. This is done as shown in
Figure 4.14.
5. After the specification of the TSItem properties, you can start to type in
the various values. A result of this is seen in Figure 4.15
Figure 4.15 Example of how the rainfall is displayed after typing in the values
1. Go to Edit|Time Series
2. Choose File|Open
3. On the dialog that comes up choose between available bridges for time
series, see e.g. Figure 4.17
85
Time Series
Figure 4.18 Dialog for specifying the import from Excel and ASCII files
6. Now you need to specify how data should be interpreted, e.g. is data a
rainfall in mm - or an intensity. First, choose the file and then specify how
time should be interpreted. In this case the time and date column is
found in one column. Hence the time description tab page will look as
seen in Figure 4.19.
Figure 4.19 The time description for the excel file in this example
7. In this example both time series are rainfall depth measured in mm. The
values are mean step accumulated values (see previous section). As the
first row in the Excel sheet is used for descriptive headings, data will firs
start in row two. All this is then specified on the data description tab
page, see Figure 4.20.
87
Time Series
Figure 4.20 The data description for the excel file in this example
8. Then go to the Preview tab page. Here you can see how your time series
is interpreted before you import the data into MIKE URBAN. If an error
has been made on the previous pages, e.g. you have specified that data
will start in row one instead of row two you will now receive an error mes-
sage and data cannot be previewed. You can now choose to either
import the data straight ahead or to save a template file first that you can
reuse when importing data of the same format.
9. The result of the import then looks like seen in Figure 4.21.
10. The save the data in the native storage format used in MIKE URBAN for
time series, i.e. the dfs0 format. This is done from the file menu of the
time series dialog. See Figure 4.22.
You can change the properties for the display of the time series. See
Figure 4.23.
Figure 4.23 In the plot the graphic settings have also been changed to make it eas-
ier to see the difference.
89
Time Series
5.1 Introduction
In MOUSE|Curves and Relations (Figure 5.1) a number of tabular data used
in other data dialogs are specified.
91
Curves and Relations
It is possible to define two different capacity curves - both used for specifying
pump operation. The capacity curve can be specified either as a 'Capacity
Curve QH' relation (for screw pumps) or as 'Capacity Curve QdH' relation (for
differential head pumps), where 'H' is the absolute water level in the pump's
wet well (at 'Location'), and 'dH' is the water level difference between the 'To'
and the 'Location' nodes. A pump type with a 'Capacity Curve QH' relation is
named a screw pump, while a pump type with a 'Capacity Curve Q-dH' rela-
tion is named a differential head pump.
Please also note that you can specify a pump capacity curve with the energy
consumption as well - for more on this type please refer to section 5.1.10.
Pumps can be RTC controlled. For those that are PID controlled, the acceler-
ation of a pump can be specified as dependent on the actual flow. This pump
acceleration curve is then specified as a number of ‘Q, dQ/dt’ data sets.
The regulation curves Qmax(H) and Qmax(dH) are used in the regulation of
the maximum discharge in links. The regulation can be either a maximum dis-
charge as a function of the water level in user specified node or a maximum
discharge as a function of the water level difference between two user speci-
fied nodes. Please refer to section 3.4.43.4 for regulation in links.
5.1.4 QH relation
QH relations can be used for various types of manholes, outlets, basins and
storage nodes.
Manholes, Basins
When using a QH relation in either a manhole or a basin it is possible to spec-
ify an infiltration (Q) based on the water level in the manhole/basin.
Outlets
Using a QH relation in an outlet means that you specify the discharge out of
the outlet based on the water level in the outlet.
Storage node
Return of water from storage basins back into the MOUSE network is speci-
fied in the dialog "Emptying Storage Nodes". An emptying function is a QH
relation, where 'Q' is the flow from the storage node into a 'receiving node',
and 'H' is the water level in a 'control node'. The flow according to the empty-
ing function is possible only until the storage node is completely empty.
The Time-Area curve is used in the Time-Area runoff model (see section 4.2).
A Time-Area curve represents the contributing part of the catchment surface
as function of time.
You can define any number of custom time-area curves. Each Time-area
table must start with a pair of values (0,0) and must end with a pair of values
representing the whole (per default, MIKE URBAN maintains T-A curves in
percent (%), i.e. the last pair of values in the table must be (100,100)).
93
Curves and Relations
There are three methods available for removal of sediments in weirs (see
also section 12.6). In one of these methods you specify the relation between
discharge towards the weir and the effiency, i.e. the efficiency curve. The
removal efficency is hence a function of Q and the efficiency (unit less 1/1).
The H-column for the basin geometry can start at any value, e.g. 0.0 for inter-
pretation of H as depth in the basin. The MOUSE Engine will associate the
first H-value to the bottom level of the node. This means that the same geom-
etry can be reused several places in the model. If we let Hmax indicate the
maximum level before flooding the highest value of the top of the specified
basin or the ground level is applied. Therefore if the top of the basin geometry
is below the ground level the specified basin geometry will be extended with
additional points to allow for flooding.
A rating curve is specified to define the relation between the valve opening
(%) and resistance (k).
Two curve types can be specified for the two different types of Curb Inlets
(see 3.7):
The QQ relation specifies the relationship between approach flow in the over-
land flow network (Qapp) and the captured flow at the connection node for an
ON-GRADE Type Curb Inlet (Qcap). Values must be monotonously increasing
and starting at (0,0). For approach discharges in excess of the maximum
value specified in the last row of the table, the last corresponding capture dis-
charge value is used.
If you have the knowledge of the power consumption of the pump at various
pump levels it is possible to add this to the pump capacity curve by using the
‘QdH & Power’ capacity curve (see Figure 5.3).
95
Curves and Relations
After the simulation with a ‘QdH & Power’ the summary will contain informa-
tion on the power consumption during the simulation period.
Figure 5.4 Snapshot from the html-summary with example of output when using
the QdH & Power capacity curve
The undefined type is used for surface runoff pollutants - SRC, FRC and TOT
dataset.
The concentrations of a given pollutant can be given for each specific catch-
ment. For one catchment it is possible to give a concentration on the slow
runoff component (SRC) and the fast runoff component (FRC) or override
these settings stating a concentration for the total runoff (TOT). Instead of
97
Curves and Relations
Boundary conditions other than water loads are "simple" in a sense that they
are fully defined by the boundary condition variable itself, i.e. no additional
item is possible to be associated with them. These are water levels, air tem-
perature, evapo-transpiration, etc.
99
Boundary Conditions
On the top level, there are three groups of boundary conditions (see fig-
ure….below). Each of these three groups is served by a separate editor,
described in full detail further below.
Figure 6.1 Basic frame of MOUSE Boundary Conditions with three groups
Air temperature
Evapo-transpiration
Rainfall
Flow - including WQ, sediment and temperature properties
An example illustrating the catchment boundary conditions is presented in
below figure
"Rainfall" is a boundary condition of the type "rainfall", and always has rainfall
(constant or a time series) associated as a default item. This type of boundary
is used for precipitation-runoff hydrological calculations. Association of WQ
properties (e.g. temperature, acidity, etc.) with rainfall is not supported. When
using a time series the dfs0-file needs to be one of the following Item Types:
‘Rainfall Intensity’, ‘Rainfall’ or ‘Rainfall Depth’. Delete values will be consid-
ered equal to 0.0.
Note that the first three boundary conditions are used by a runoff-precipitation
model, while the water load boundary condition can either be appended to
the computed runoff (only flow) or passed directly to the network connected
to the catchments.
101
Boundary Conditions
The following items can be toggled on/off under ‘MOUSE Catchment Bound-
aries’ in the TOC:
Water level
Flow (discharge) - incl. WQ and sediment properties
The discharge boundary conditions are termed as "Network Loads". An
example illustrating the network loads is presented (see illustration below)
Network Loads represent all kinds of hydraulic loads, such as DWF point-
inflows, individual lateral inflows, infiltration, storm runoff, etc. This type of
boundary condition always has discharge (constant, cyclic or a time series)
associated as a default item. The discharge can be associated with any pol-
lutant, sediment or temperature item (constant, cyclic or a time series).
"Water Load.Storm Runoff" represent calculated storm runoff loads. Note that
in the example above, two default items are specified, each covering different
time period. The default items (i.e. storm runoff itself) are associated with
temperature and pollutants BOD and NH4.
103
Boundary Conditions
"catchment load" (in this case "Water Load.Resident DWF") and a signal to
the model to pass the catchment load to the network model. This type of load
definition does not need any item. Please note that this is only valid for the
MOUSE engine. The MIKE 1D engine does not require this step.
Only one water level boundary condition is allowed at each network outlet. If
nothing is specified for an outlet, the model assumes a constant water level
equal to the outlet invert level. I.e. a free outlet is assumed.
This type of boundary condition always has water level (constant or a time
series) associated as a default item.
In the example above, two external water level boundary conditions have
been defined.
The network loads are visualized by toogling the display of ‘MOUSE Network
Boundaries’ on in the TOC.
The following items can be toggled on/off under ‘MOUSE Network Bounda-
ries’ in the TOC:
To sort the different loads (catchment loads, network loads and external water
levels) in each editor, use the right click option "Sort" on top of the "Boundary
ID" column. Double clicking on the column will display all boundaries condi-
tions loaded in the database indistinctly of the editor type where the user is
positioned.
105
Boundary Conditions
As the first action after inserting a new boundary conditions, a proper name
(ID) must be specified. We recommend that you use a descriptive ID.
To sort the different loads use the right click option "Sort" on top of the
"Boundary ID" column. Double clicking on the column will display all bounda-
ries conditions loaded in the database.
The following describes how to define an boundary condition for the network.
This could e.g. be a discharge entering a manhole or pipe infiltration.
As the first action after inserting a new boundary condition, a proper name
(ID) must be specified. We recommend that you use a descriptive ID.
The "Connection Type" determines the mode of definition of the network loca-
tion(s) to which this boundary applies. "Individual" applies for a single node or
pipe, and requires a specification of the Node ID or Link ID.
"List" applies for a set of node or links, contained in the selection files (*.NSE,
LSE or *.MUS). "Defined by Catchment Connections" applies to the storm
runoff loads from MOUSE *.CRF file.
107
Boundary Conditions
the node or loss of water from the node to the surrounding ground can be
used.
The loss can be defined for any type of node, but it is primarily intended for
use with basins or nodes with geometry representing small ponds. The infil-
tration can be defined as a constant rate (velocity) which is multiplied by the
actual wet area in the node. Note that the rate value should be given as a
negative value for specifying a loss of water from the network model. Apply-
ing a positive value will result in water being added to the node.
Please note that when you use the boundary type “Defined by catchment”
and afterwards specify a time series ID in the boundary items dialog then this
item is always used regardless of type (water, component) or extension (dfs0,
CRF). If you have not specified a time series ID then the catchment name will
be used. This is only applicable for the boundary type “Defined by catchment”
for other types the time series ID is always used.
To sort the different loads use the right click option "Sort" on top of the
"Boundary ID" column. Double clicking on the column will display all bounda-
ries conditions loaded in the database.
To sort the different loads use the right click option "Sort" on top of the
"Boundary ID" column. Double clicking on the column will display all bounda-
ries conditions loaded in the database.
Boundary item editor allows for a full definition of a boundary item in terms of
type, temporal variation, quantity/quality and temporal validity
109
Boundary Conditions
The 'Boundary Items' editor has been designed to operate in conjunction with
the three boundary condition editors. When open simultaneously with one of
these, this editor displays only the boundary items associated with the cur-
rently active boundary condition. Otherwise, displays all boundary items.
'Boundary Items' editor allows for a full definition of a boundary item in terms
of type, temporal variation, quantity/quality and temporal validity.
A new boundary item is created by pressing "Insert". The new item is auto-
matically associated with the currently active boundary condition and is
appended to the list of associated items. Alternatively, already existing item
can be re-used and associates it with the currently active boundary. "Re-use
Item" actually makes a copy of a selected item and adds it on the list of items
for the current boundary condition. By these means user can assign the same
item (e.g. concentration of some pollutant in wastewater) to many boundary
conditions quickly.
A new item can also be inserted while none of the boundary condition dialogs
is open. In that case, user must "manually" connect the item to a boundary
condition. There is no automatic control over the match between the bound-
ary condition type and the actual item definition, so due care is needed.
Each new item is set as "default". For some types of boundary conditions, this
can be changed to "Pollutant" or Sediment". In the later case, a pollutant
component or a sediment fraction must be selected from the respective lists.
The pollutant items can be defined in terms of concentrations or loads.
The scaling factor determines the fraction of the specified load which is actu-
ally loaded into the model. Per default, the scaling factor is set to 1.
Temporal variations
Temporal variation of boundary items may take one of the following forms:
Time series may come from various sources, including previous MOUSE
result file (this latter option includes MOUSE files CRF, SMR and PMR).
Temporal validity
Each boundary item may be specified with limited temporal validity.
Per default, each item is considered valid during the entire simulation period.
Indeed, in case of TS items, the validity is also limited by the length of TS.
If multiple items are specified they may all be with default temporal validity
(i.e. over the entire simulation period), or with overlapping validity periods
(e.g. when mixing de-fault and limited validity) or with periods with missing
coverage.
6.4 Examples
In the following paragraphs some examples on setting up boundary condi-
tions are shown step-by-step
This boundary condition falls into the category 'External Water Levels'. The
definition involves the specification of a boundary condition itself and the
water level (default) item.
111
Boundary Conditions
This boundary condition falls into the category 'Network Loads'. The definition
involves the specification of a boundary condition itself and the infiltration
item.
Press the Items button to open the Boundary items dialog to add the infiltra-
tion value. Note that the Boundary ID is automatically filled in on the Bound-
ary items dialog.
This boundary condition falls into the category 'Catchment Loads And Mete-
orological Items'. The definition involves the specification of a boundary con-
dition itself and the rainfall (default) item.
Choose Time Series/Result File in order to choose a time series for the
rain event
This activates the Time Series part of the dialog
113
Boundary Conditions
Choose DFS0 as the Source Type, and select the dfs0 file with the rain
event (item ID is only necessary to fill out if the dfs0 file contains more
than one item, i.e. time series)
Now the rain event has successfully been connected to the catchments!
This boundary condition falls into the category 'Network Loads'. The definition
involves the specification of a boundary condition itself and the discharge
(default) item.
6.4.5 How to add runoff results as input for the network computation?
The runoff results are treated as network loads in MIKE URBAN, and can be
specified following the same principle as any other load (see the example
above). Alternatively, a 'shortcut' specification can be used.
Now press the Items button. This will open the Boundary Items dialog
where the runoff result file to be used can be specified
Press Insert button on the Boundary Items dialog. The Boundary ID field
is automatically filled in by the current boundary ID. This ensures a bond
between the Boundary condition and this item.
After having pressed Insert an item ID will automatically be given (hint:
you can change this ID if you wish to).
Specify that the load type is 'Storm Runoff'. This will insure a correct cat-
egorisation of this load in the simulation summary statistics and in the
job-list generation for LTS simulations.
Now the temporal variation can be specified.
Choose Time Series/Result File in order to choose the result file from the
catchment simulation
This activates the Time Series part of the dialog
Choose MOUSE as the Source Type, click the list button next to the Path field
and select the *.CRF file containing the results from the runoff computation.
Now the runoff results are used as input for the network computation!
115
Boundary Conditions
Please note that the steps are different depending on your choice of numeri-
cal engine (MOUSE or MIKE1D).
The PE-based and catchment area-based DWF fall into the category of
catchment loads. This is because both the input parameters and the results
of the computation are catchment-related. However, in most applications, this
load is supposed to be passed into the network through the catchment con-
nection points. The following example describes the entire definition process.
After having pressed Insert an item ID will automatically be given (hint: you
can change this ID if you wish to)
The scaling factor and item type are per default set properly, so no action is
required. Now the temporal variation can be specified:
This completes a full definition of the catchment DWF load based on PEs.
The remaining step (which is only needed in case of using the MOUSE
engine) is to allow this load be directly applied to the network model.
117
Boundary Conditions
Press the 'Items' button. This will open the Boundary Items dialog where
items associated with the actual load will be displayed. There will be at
least one discharge item and, possibly, some pollutant and sediment
items.
Press Insert button on the Boundary Items dialog. The Boundary ID field
is automatically filled in by the current boundary ID. This ensures a bond
between the Boundary condition and the new item.
After having pressed Insert an item ID will automatically be given (hint:
you should change this ID in order to be able to identify it easier)
Specify the new item type as 'Pollutant'. This will activate the part of the
editor where a pollutant must be selected from the list. Note that pollutant
can be specified as concentration (i.e. mass per unit volume of water) or
as load (i.e. mass per unit of time). This specification is not effective for
the cases where actual pollutant data include the information of actual
unit (e.g. dfs0 time series).
Now the temporal variation can be specified. Note that specification of
cyclic variation in a combination with "concentration" affects the mass
balance, i.e. the simulated pollutant mass changes with changes in cyclic
pattern. This must be compensated by appropriate adjustment of the pol-
lutant concentration.
Now a pollutant (concentration or load) has been successfully connected to
your network load!
An unlimited number of diurnal patterns can be entered. E.g. one for week-
days and one for weekends, but also one for weekdays belonging to one
group of catchments and one for weekdays belonging to a different group of
catchments. Clicking on ‘Graph’ gives a graphical display of the profile.
Three checkboxes 'Week Days', 'Dates' and 'Months' are available for the
activation of three different, mutually inclusive modes of the calendar defini-
tion. This means that if any day is to be included in the current calendar defi-
nition, it must be included in all three definition modes. If any of the definition
modes is not activated (i.e. the checkbox not ticked), it does not affect the cal-
endar definition. If only a calendar's ID is given, but no tick marks are set, the
diurnal pattern that the profiles calendar is assigned to would be uncondition-
ally used.
119
Boundary Conditions
For each pattern the 'Interpolation' method is chosen to either 'Linear Interpo-
lation' or 'No Interpolation'. The 'Linear Interpolation' performs a linear inter-
polation between the values given in the diurnal profile while the 'No
Interpolation' will apply a step function. For use with DWF, the 'Linear Interpo-
lation’ is recommended
E.g. the 1st of January each year (the New year holiday) even though it
comes every year on a different day of the week could be coupled with the
diurnal profile of e.g. Sunday, since the diurnal pattern of a holiday matches
better with Sunday than by any other week day.
The special days are divided in two categories: Special days that are to be
considered every year and special days that are only considered once
(Unique date). An example of this could be the Easter holiday which will not
be on the same date every year.
121
Boundary Conditions
7 MOUSE Simulations
MOUSE simulations may be started from ‘Simulation | Run MOUSE’. Simula-
tion options for MOUSE are also defined here.
General
Runoff parameters
Network parameters
Network Summary
2D overland parameters
The ‘Simulation ID’ and ‘Scenario ID’ are found on all four tabsheets as a
common reference.
123
MOUSE Simulations
‘Info’ button
The ‘Info’ button opens a dialog containing:
1. Start and end time of the simulation similar to the values found under
max time
2. Start and end time of a hotstart file (if specified) (Table: msm_Project,
Attribute: PRFHotstartFileName)
Figure 7.2 Example of the dialog that comes up, when pressing the ‘Info’ button
Each boundary contains a number of items which can cover different parts of
the simulation.
If a validity is specified for an item this specifies the start and end time. If a
validity is not specified only items specified as timeseries have a start and
end time. If either a constant or cyclic value is given without validity the item is
not included in the evaluation.
125
MOUSE Simulations
For each simulation you can choose to use either the classic MOUSE engine
or the MIKE 1D engine. Please refer to seperate documentation on the differ-
ence between these two engines and the possibilities provided with them.
It is important to note that when running the numerical engine of you choice
you need to keep the input for the simulation within the same numerical
engine. I.e. when using results from a runoff simulation as input to a pipe flow
simulation then both the runoff and the pipeflow simulation needs to be done
with the same numerical engine. If runoff result file is computed with e.g. the
MOUSE engine and you try to use it as input for MIKE 1D pipeflow simulation
you will hence get an error message. You can, however, create a number of
simulation records and then choose to have some done with MOUSE and
some with MIKE 1D.
The MIKE 1D result files all have the same extension (*res1d), while the
MOUSE result files have different extensions based on simulation type. Sev-
eral places you will hence see both result types appear as a filter for file types
you can choose between:
Computation dialog - RDII hotstart file. The filter says: "Hotstart Files
(*.nof, *.res1d)
Computation dialog - runoff input. The filter says: "Runoff Files (*.crf,
*.res1d)
Computation dialog - hotstart file. The filter says: "Hotstart files (*.prf,
res1d)"
Boundary Items - result file/time series. Both MOUSE Result Files and
res1d files can be used
LTS initial conditions - hotstart file. The filter says: "Hotstart Files (*.prf,
res1d)"
Relevant information from a hot start file for a network simulation includes:
For the RDII model, relevant information in the hot start file includes all rele-
vant hydrological state variables, such as surface storage, groundwater
depth, etc.
A hot start file for a network simulation must be generated with exactly the
same model setup as the model used in a hot start simulation run. A result file
used as a hot start file must be a complete ("standard") result file. I.e. user-
specified result file with reduced contents cannot be used as a hot start file.
For the RDII runoff model, the only parameter is the Hot start Filename. The
system automatically detects the values of the relevant variables at the simu-
lation start time and initiates the actual simulation run accordingly.
For the pipe flow model, the following parameters are needed:
Hot start Filename: This is the name of an existing result file from which
initial conditions for the simulation will be loaded.
Hot start Date and Time: The date and time at which the initial conditions
are taken from the hot start file. The hot start time must be within the hot
start file time coverage.
127
MOUSE Simulations
If the Hot Start Filename specified is identical with the current simulation
result file name (concatenated by simulation MUID, ScenarioID and the
extension .PRF),
The results of the current simulation will be appended to the hot start file,
starting at the simulation start time. In this case, hot start date and time are
ignored, as the simulation hot starts at the simulation start time. Any informa-
tion in the hot start file beyond the simulation start date and time (the actual
hot start date and time) will be lost. This part of the file will be replaced by the
new simulation results.
To create an engine specific input file, open the ‘General’ tab of the ‘Compu-
tation’ dialog and click ‘Advanced’ . This opens a menu where there is an
option called ‘Export to Remote Simulation’:
Selecting ‘Export to Remote Simulation’ will export the setup to a *.m1dx file
or a *.mex file, depending on the choice of simulation engine. The full name
of the exported file will be composed of the Simulation ID, the Scenario ID
and the extension.
129
MOUSE Simulations
131
MOUSE Simulations
Identification of simulation
Simulation mode
LTS Job list
Runoff Input
Dynamic simulation
Hot Start
Processes
Results
133
MOUSE Simulations
The additional parameter file (the ADP-file) is a separate file with additional
settings for the simulation. Please refer to the seperate documentation on this
file for further information.
Identification of simulation
Tables with min, max and accumulated values for
Summary of input data
Use selection to reduce summary tables
135
MOUSE Simulations
137
MOUSE Simulations
By ticking the 'Save volume data' checkbox, data regarding mass balance in
the system will be saved in the result file. If ticked, 'Save time step data' will
save information about the used time steps in the simulation.
It is possible to save data in a time period that does not necessarily cover the
whole period of the simulation. Make sure that the interval specified by 'start'
and 'end' dates is at least partly covered by the simulation period. Saving fre-
quency specified for the 'Reduced' result file is not dependent on the time
interval chosen for the complete result file, if such a file has been specified.
Results can be saved from four categories; Nodes, Pipes, Weirs, Orifices and
Pumps. Consequently, the result save specification is divided into four sepa-
It should be noted that if a set of nodes is selected, and the selection (from
the horizontal plot or from the selection file) is loaded into the pumps or weirs
list, only nodes where pumps or weirs do actually exist would turn up on the
list.
In cases where several pumps or weirs are associated with one node, the
results will be saved for all these pumps or weirs.
Presence of a certain item on the list is not sufficient for it to be included into
the user specified result file. This is first ensured after a result type is chosen
e.g. water level. In figure Figure 7.9 is displayed the Editor for ‘Nodes result
selections’
The 'Mark List' option helps clearing or checking all items on the displayed
list. If a queried list is displayed, 'Mark List' will operate only on the currently
displayed list. The remaining part of the list will retain earlier settings.
When a result selection file has been specified then it is activated in the ‘Run
MOUSE | Network parameters’ editor by selecting the Selection ID in the
‘Results’group in the field labeled ‘Contents:’ - see Figure 7.10
139
MOUSE Simulations
8 2D Overland Flow
8.1 Introduction
Surface flooding can be simulated with MIKE URBAN using the following two
approaches:
141
2D Overland Flow
format or DEMs provided in the .dfs2 file format. The latter is the DHI file
format for raster data sets.
3. A number of couplings between the 1D model and the 2D model. Such
couplings are connections between the two models allowing for
exchange of water between the two model domains. Couplings can, for
example, represent manholes or curb inlets.
Select MOUSE | 2D Overland Flow | 2D Model Settings ... to open the dialog
shown in Figure 8.1.
The underlying 2D overland flow model, MIKE 21, offers the following differ-
ent types of 2D model solvers:
1. Single Grid using a rectangular cell solver. The 2D model area is covered
by square calculation cells.
2. Single Grid using a rectangular multi-cell solver. Again, the 2D model
area is covered by a grid of square calculation cells. This is referred to as
the coarse grid. Each coarse grid cell is further subdivided into N by N
cells called the fine grid. N is an integer number - the ratio between the
linear cell dimension of the coarse and the fine grid.
3. MIKE 21 Flood Screening Tool (FST) is designed for inland flooding
applications. By means of FST the user can run 2D overland flood simu-
lations faster than the regular MIKE 21 overland flow simulations. The
Flood Screening Tool (FST) provides a simple hydrodynamic solver for
computations regarding inland flooding, is based on MIKE 21 Flow
Model but uses simpler equations with some parameters from MIKE 21
Flow Model being excluded. In FST solver the convective terms are omit-
ted from the model, the viscous terms are omitted in the governing equa-
tions, and the influence of wind, atmospheric pressure variation, Coriolis
and wave radiation stresses cannot be included.
For further information on the pros and cons of each of the 2D solvers, please
refer to the MIKE 21 documentation. For more information regarding the
equations employed in MIKE 21 FST please consult MIKE 21 Flow Model &
MIKE 21 Flood Screening Tool Hydrodynamic Module - Scientific Documen-
tation.
Enabling MIKE 21 Flood Screening Tool automatically disable the Eddy Vis-
cosity settings which are consequently greyed out in the "2D model settings",
see Figure 8.2.
143
2D Overland Flow
The DEM to be applied is loaded to the map view using the Insert Layer but-
ton in the Layer Management toolbar or by choosing Insert Layer in Back-
ground Layers within the TOC. Both these options will accept DEM raster
data stored in ESRI Grid format or in DHI dfs2 format. Depending on your
GIS installation other raster formats could be available.
1. Select the desired raster layer as the DEM to be applied in the model.
2. Enter the X- and Y-coordinates for the lower left hand and upper right
hand corner of the rectangular area to be modelled with the 2D model.
Instead of manually entering the coordinates the buttons ‘Digitize rectan-
gle’, ‘Set as extent of data set’ or ‘Set as extent of raster’ can be used to
This dialog allows for the editing of the following 2D model settings:
145
2D Overland Flow
Eddy viscosity. The eddy viscosity influences how eddies are simulated
as a result of spatial velocity or flux gradients. Higher values of the eddy
viscosity may also dampen numerical instabilities, and increasing the
eddy viscosity is thus sometimes used as a method for stabilizing mod-
els. The eddy viscosity can be specified in three ways:
– It can be included in terms of an automatically calculated, flux based
value equal to 0.2.x.y/t, where x and y are the model cell size
and t is the simulation time step.
– It can be excluded.
– It can be included as a user defined value, using either a flux based
or velocity based equation.
Miscellaneous. Read Section 8.3.9.
The initial water level (water surface elevation) applied in the 2D model can
be specified as a constant value applied throughout the 2D model area or as
spatially varying for which two options are available:
147
2D Overland Flow
The specified initial water level (spatially constant or varying) will only apply
to 2D cells with a ground elevation less than the initial water level. 2D cells
with higher ground elevation will initially be dry.
Spatially varying initial conditions for the 2D model are supported in two
ways:
The form and grid view dialog for the table holding the polygons opens when
choosing MOUSE | 2D Overland Flow | 2D Initial Conditions … as shown in
Figure 8.5.
Figure 8.5 Form and grid view for 2D Initial Conditions dialog
Table 8.1 Editor fields in form and grid view for 2D Initial Conditions dialog
149
2D Overland Flow
Table 8.1 Editor fields in form and grid view for 2D Initial Conditions dialog
Spatially varying bed resistance for the 2D model are supported in two ways:
The form and grid view dialog for the table holding the polygons opens when
choosing MOUSE | 2D Overland Flow | 2D Bed Resistance … as shown in
Figure 8.6.
Figure 8.6 Form and grid view for 2D Bed Resistance dialog
Bed resistance polygons are applied to represent land use or other ground
surface classification with a strong relation to bed resistance. The polygon
table is typically populated by importing land use polygons for instance from a
number of Shape files each holding multiple polygons with the same land
use. The description field in the table can be used to store the land use and
thereby, with the selection and field calculator tools, ease bulk changes of the
bed resistance for all polygons of the same land use.
151
2D Overland Flow
Table 8.2 Editor fields in form and grid view for 2D Bed Resistance dialog
8.3.8 2D boundaries
By default the 2D model has a closed boundary and water will not enter or
leave the 2D model area across the perimeter of the 2D model areas. How-
ever, it is possible to define the following types of 2D boundaries:
Rain: This applies in cases where rainfall should enter each 2D compu-
tational cell at each time step. Both constant rain and rain with spatial
and temporal variation are supported. Rain directly on 2D cells is used
as an alternative to using one of the rainfall-runoff models in MIKE
URBAN.
Water level: This applies to a stretch of the perimeter of the 2D model
area and is typically applied where the 2D model perimeter crosses a
larger water body such as a lake or the sea where the water level is
known (constant or time varying). However, a water level boundary can
also be applied as a free outlet allowing water to leave the 2D model
area rather than backing up against the closed boundary.
Discharge: This applies to a stretch of the perimeter of the 2D model
area and is typically applied where the 2D model perimeter crosses a
river or similar where the discharge is known (constant or time varying).
Boundaries for the 2D model is defined with the form and grid view which
opens when selecting MOUSE | 2D Overland Flow | 2D Boundaries ... , see
Figure 8.7.
Figure 8.7 Form and grid view for 2D boundaries dialog. The arrow points to the
graphical line showing the location of the boundary
153
2D Overland Flow
points”. If the user manually types in the coordinates then snapping will
occur when choosing Advanced | Recompute, when the boundary is
shown on the map or when a simulation is run. The location of 2D bound-
aries can be visualized in the map view by means of a red line between
the snapped boundary end points. See Figure 8.7.
3. Value: The value of the boundary condition can be defined as constant
or time varying with a .dfs0 file. Rain boundaries can additionally have a
spatial variation in which case the data has to be provided in a .dfs2 file.
The discharge boundary is to be specified with positive numbers for flow in
the same direction as the X and Y axis. This means that flow into the model
area on the northern and eastern boundaries has to be specified as a nega-
tive number.
It should further be noted that 2D boundary files (dfs0 and dfs2) should be
specified as floats and with an equidistant calendar axis see Chapter 4 ‘Time
Series’ for configuration. For a rainfall boundary the item type should be set
to “Rainfall rate” and must only include one item in the file. For a spatially var-
ying rain specified with a dfs2 file the file must have the same number of cells
and cell size as the DEM.
Table 8.3 Editor fields in form and grid view for 2D Boundaries dialog
Table 8.3 Editor fields in form and grid view for 2D Boundaries dialog
Adjustment of DEM
Along the 2D model perimeter the ground elevation in the first row/column of
2D cells is automatically raised to a high elevation (elevation threshold for
active cells) to ensure zero flow across the perimeter. However, this does not
apply along stretches where a 2D water level or discharge boundary is
defined. Instead, the ground elevation is lowered automatically to ensure suit-
able numerics of the simulation at the boundary, because a boundary must
be wet all the time. By default the lowering is 0.5 meter and it applies to the
first two rows/columns parallel to the 2D boundary. All 2D cells are lowered to
the same value which equals the lowering value subtracted from the lowest
cell in the first non-adjusted row/column (by default the 3rd) along the 2D
boundary. The lowering and the width of the indentation can be adjusted by
clicking the “Additional” button in MOUSE | 2D Overland Flow | 2D Model Set-
tings ... Figure 8.8 shows an example of how the DEM is adjusted along a
boundary.
155
2D Overland Flow
Defining the couplings between the pipe flow model and the 2D overland flow
model involves the following steps:
The following set of tools found in the sub-menu MOUSE|2D Overland Tools
(see Figure 8.9) are available for defining the couplings:
Coupling parameters
Selecting this will open the dialog shown in Figure 8.10.
When coupling a pipe model node to the 2D model, a coupling area
covering the spatial extent of the 2D cells involved in the coupling has
to be defined. The coupling area will automatically be defined by the
tool "Couple nodes". The selection in this dialog defines how many 2D
model cells the coupling area will span. There are two choices. 1)
Either a squared coupling area with the width/height as the specified
number of cells. 2) A coupling area defined by all 2D cells within the
specified search radius from the node in question. If no 2D cells will
have a cell centre within the specified coupling area, then the 2D cell
containing the node in question will be selected as the coupling area.
157
2D Overland Flow
Figure 8.10 Dialog for parameters for the 2D overland flow tools
Couple nodes
This tool applies to manholes, basins and outlets that are not already
coupled. If any nodes are selected the tool will apply only to those
selected nodes - otherwise to all. Once a node is coupled the following
will happen:
– The node will be marked as coupled and it will show-up as such in
the tab ‘2D Overland’ in the dialog MOUSE|Nodes.
– A coupling area polygon will be defined and an entry will be added to
the table ‘Coupled 2D Areas’. If selected for drawing the coupling
area will show up in the map view.
– A line between the node and the centre of the coupling area will be
defined and an entry will be added to the table ‘Coupled 2D Connec-
tions’. If selected for drawing the connection line will show up in the
map view.
Figure 8.11 shows an example of three nodes which have been cou-
pled to respectively a single 2D cell, to 3x3 2D cells and to the 2D cells
within a circular search radius.
Once a coupling area has been defined automatically, its shape and
location can be manually edited using the tools in the Feature Edit tool-
bar.
Un-coupled nodes
If any coupled nodes are selected this tool will apply only to those - oth-
erwise to all coupled nodes. The tool will un-couple nodes by taking the
following actions:
– The node will be marked as un-coupled and show-up as such in the
tab ‘2D Overland’ in the dialog MOUSE|Nodes.
– The associated coupling area polygon will be deleted.
– The associated coupling connection will be deleted.
159
2D Overland Flow
Once a manhole or basin has been coupled, the node attribute ‘Coupling to
2D Overland’ will be selected. See Figure 8.13. On this tab the user can set
various parameters for how the exchange flow between the pipe flow model
and the overland flow model is calculated. For more information, please refer
to MIKE FLOOD User Manual.
Max flow
This parameters gives an upper limit to the discharge that can flow
through the link. If the discharge calculated from the flow equation
exceeds this value, then the discharge is reduced to this value.
Qdh factor
A parameter indicating the water level at which the calculated discharge
should be surpressed. May be used for stabilizing setups where the
instability is caused by a small difference in the water level in the two
coupled models. The suppression is calculated as the water level differ-
ence divided by the Qdh factor to the power of 2.
161
2D Overland Flow
Calculation method
The exchange flow between the 1D and 2D models may be calculated using
4 different methods:
Inlet area
The inlet area is only used when applying the orifice equation for describ-
ing the flow exchange between the 1D model and the 2D model. The
greater the cross sectional area the greater the conveyance capacity of
the coupling. This parameter corresponds physically to the area of the
manhole cover. For a manhole the inlet area used in the calculation will
be the smallest value of the specified inlet area and the inlet area calcu-
lated from the manhole diameter. For a basin the specified inlet area will
always be used.
Orifice coefficient
The orifice coefficient is a non-dimensional factor that may be used to
scale the orifice flow.
Crest width
The crest width is used when the flow exchange is described through a
weir equation. The weir crest should typically be the circumference of the
manhole cover.
Scaling factor
This factor is only applicable if the simple exponential expression is used
to describe the flow exchange. The scaling factor has a linear effect on
the discharge. Thus, an increase in this value will generate a larger flow
for a certain water level difference between the 1D and 2D model.
Exponent
Only applicable if the simple exponential expression is used to describe
the flow exchange. An increase in the exponent factor has a strong
impact on the discharge. An increase in this value will generate a larger
flow for a specific water level difference between the 1D and 2D model.
Freeboard
The freeboard (see Figure 8.14) defines a critical water level (Invert
level- Freeboard) at the connection node in the pipe network below
which the defined DQ relation apply. For submerged and reverse flow
(surcharge), the transfer capacity of the connection reverts to a standard
orifice relationship.
The freeboard describes the distance (depth) from the Grade Invert Level to
the Critical Water Level (i.e. Freeboard = Grade Invert Level - Critical Water
Level). The Critical Water Level is the level above which the basin is consid-
ered to be 'drowned' such that the DQ relationship for transferring water from
the surface to the node no longer applies. For submerged and reverse flow
conditions the exchange of water across the coupling reverts to a standard
orifice equation.
163
2D Overland Flow
8.3.12 Outlets
An outlet coupled to the 2D overland flow model does not require other
parameters than the definition of the coupling area itself. The average water
levels in the 2D cells to which the coupling is made will then be applied as a
water level boundary. The same algorithm will apply as when a water level
boundary is specified in an outlet in the pipe flow model. This will result in an
outflow through the outlet (positive or negative) which will correspondingly
apply as a point inflow (positive or negative) distributed to the 2D cells in the
coupling area. For more information, please refer to MIKE FLOOD User Man-
ual.
To run the combined pipe flow and overland flow simulation three issues
needs to be addressed in the simulation dialog as shown in Figure 8.15.
Firstly, the simulation type needs to be "Network + 2D overland". Secondly,
the time step in the 2D model and the request for output from the 2D model
have to be defined in the dialog "2D Overland Results" as shown in
Figure 8.16. This dialog can be opened from the "2D Results" button located
in the "Network parameters" tab.
165
2D Overland Flow
Figure 8.16 Selection of time step and request for output in 2D model
In the tab ‘2D Overland Results’ the 2D time step is specified. The request for
output is specified as a single 2D output file with maximum flood depth and in
the table below, a number of output files with the flood depth at multiple time
steps through the simulation. For the latter the following has to be specified
for each requested result file:
Type The choice between the following output files both containing
multiple time steps:
Flood maps: A 2D output file in .dfs2 raster format.
Time series: A time series output file in .dfs0 format.
Model area This is only enabled for flood map output files. The user can
choose be-tween:
Model area: The output file will then cover the full 2D model
area.
Selected area: The output file will then cover only a sub-
region of the full 2D model area. The sub-region is a rectan-
gle defined by the following two set of X,Y coordinates:
X0, Y0 Type = Flood Maps: Only enabled if the model area is cho-
X1, Y1 sen as selected area. The lower left and upper right corners
of the sub-region from which 2D output will be generated is
defined by the two coordinate sets (X0, Y0) and (X1, Y1)
Type = Time series: Results are saved from the location (X0,
Y0). (X1, Y1) is not used and thus disabled.
FileName Name of the result file to be saved. The file name cannot be
changed by the user.
Depth, Level, Option for saving water depth, level (water surface eleva-
Velocity & Flux tion), velocity and flux (P, Q) respectively
167
2D Overland Flow
The standard MIKE 21 results (H,P,Q) will be generated when Depth and Flux
are chosen. If you would like more information on the result type flux (P and
Q), you can check the MIKE 21 documentation.
After pushing the button Start Simulation in the dialog shown in Figure 8.15
MIKE URBAN data will be exported as for a regular MOUSE simulation. Fol-
lowing that export, the MIKE FLOOD simulation progress window shown in
Figure 8.17 will appear.
The .dfs2 raster files may contain a single time step or multiple time steps.
The former is the case for files containing the maximum flood extend. Such
should be added as a raster layer to the map view using the Insert Layer but-
ton in the Layer Management toolbar. Dfs2 files with multiple time steps are
visualized as explained in the following.
In the main menu bar, choose Model Result > Load Results > Overland
Results in order to select the 2D overland result file. This will add information
about the file content to the MIKE URBAN session, but as .dfs2 files are
potentially very large the data from the individual time steps are not loaded to
memory. Thus, nothing will be added to the TOC or the map view when a
.dfs2 files is selected. There are three ways to subsequently visualize the
content of a .dfs2 file:
1. As a raster layer in the map view. To do so, select Model Results > Add
Results Layer. Subsequently, select the desired.dfs2 file and data type
which then will be added to the TOC and map view. The map view will
show data from the time step selected with the Animation toolbar. See
Figure 8.18.
2. As a time series plot extracted at a particular location. To do so, in the
Simulation toolbar click on the Time Series button and then select the
desired .dfs2 file and data type. Subsequently, move the cursor to the
location in the map view from where a time series should be extracted
from the .dfs2 file. See Figure 8.18.
3. As a water surface elevation line in a longitudinal profile view. To do so,
firstly a longitudinal profile view should be opened. In here right-click and
choose Properties. In the Properties dialog the Raster tab allows for add-
ing raster data to the profile view. See Figure 8.19.
Figure 8.18 2D flood results shown as raster data and time series
169
2D Overland Flow
Figure 8.19 2D flood results shown as water surface profile in longitudinal profile
<Simulation ID> MIKE FLOOD This file is an input file for the MIKE
<ScenarioID> input file FLOOD simulation and it can be edited
.couple further with the MIKE FLOOD editor in
MIKE ZERO. In short, the file contains
reference to the MIKE 21 input file and
the MIKE URBAN .mex file as well as
how these two models are coupled.
<Simulation ID> MIKE 21 input file This file is an input file for the MIKE 21
<ScenarioID> component of the MIKE FLOOD simula-
.m21 tion and it can be edited further with the
MIKE 21 editor in MIKE ZERO. In short,
the file contains the parameters for MIKE
21 as specified in the MIKE URBAN dia-
log for 2D Overland Flow.
<Simulation ID> Raster file in dfs2 This file contains the DEM data for the
<ScenarioID> format 2D model area selected. Data in the file
Coarse is extracted from the original DEM speci-
fied in MIKE URBAN at the requested
.dfs2
model resolution.
<Simulation ID> Raster file in dfs2 This file contains the DEM data for the
<ScenarioID> format 2D model area selected. This file is only
Fine generated if the Multi Cell 2D model has
been chosen. Data in the file is extracted
.dfs2
from the original DEM specified in MIKE
URBAN at the requested model resolu-
tion divided by the fine/coarse ratio.
<Simulation ID> Raster file in dfs2 This file contains spatially varying 2D ini-
<Scenario ID> format tial conditions.
InitialCondi-
tions.dfs2
<Simulation ID> Raster file in dfs2 This file contains spatially varying 2D
<Scenario ID> format bed resistance
Bed_Resist.dfs2
171
2D Overland Flow
Pump
Weir
Orifice with gate
Orifice with weir
Valve
The devices may be specified as settings or PID-controlled, with control
function selection based on a global system analysis. Each regulator or pump
operates under the control logic encapsulated into a set of simple logical
rules and control functions. The system allows a schematised definition of
any form and size of decision tree, featuring logical operators AND, OR, NOT
and NOR, in any combination.
The control functions range from the simplest constants for the operational
variables (e.g. constant weir crest setting or constant flow set-point) to
dynamic controlled variables, set in a continuous functional relation with any
of the measurable variables in the system (e.g. CSO discharge set-point as a
function of flow concentration or a pump START/STOP levels as functions of
water level at strategic location in the system).
173
MOUSE Control module
where typical rainfall patterns over the catchment area exhibit high
degrees of spatial variability resulting in some parts of the system
becoming overloaded whilst others are under utilised;
where the urban wastewater system includes treatment processes
whose performance is amenable to active, short term control;
where the assimilative capacity of the receiving waters is variable over
time.
Usually, RTC is implemented as an integral part of a rehabilitation/upgrade
scheme also involving significant civil upgrading works to increase the trans-
port, treatment or assimilative capacity of the urban wastewater system. In
such circumstances, the role of RTC is to optimise the operation of both the
new and the existing facilities, thereby maximising the benefit in performance
terms. Where the overall objective is to achieve compliance with specified
performance targets, RTC serves to minimise the scale and extent of the nec-
essary works.
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 comprehen-
sive way to obtain forecast data is to include a model in the control system.
Predictive 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).
175
MOUSE Control module
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 MOUSE network. Sensors with multiple function-
ality must be specified individually.
When devices (weirs, gates and pumps) are specified as controllable in
the MOUSE 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 actions) is associated with the respective devices.
MOUSE RTC uses sampling and actuation (control loop) frequency iden-
tical 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.
9.4 Sensors
A sensor is a physical device positioned somewhere in the system, which
provides information of the actual value of a monitored variable. A sensor can
only monitor one variable – if more variables are measured at the same local-
ity a corresponding number of sensors has to be described.
Sensor ID
Each sensor needs a unique ID, which can be used to access the sensor
information from other dialogs.
Type
The type of the sensor defines which variable the sensor measures.
Component
For the "Concentration" sensor type, a measured WQ component must be
specified.
177
MOUSE Control module
Condition ID
Each Logical Condition needs a unique ID, which can be used to access the
logical condition information from other dialogs.
For each of the tests included in a logical condition the following should be
specified.
Type and ID
The type and identification of the variable which should be evaluated.
Depending on the type, it might be necessary with two IDs –e.g. for the eval-
uation of difference between two sensor values.
For the type "Diff of 2 Sensor Values" the sensors specified under ID1 and
ID2 must be of the same type and sensors measuring either level, discharge
or concentration.
For the type "Sum of 2 Sensor Values" the sensors specified under ID1 and
ID2 must both be sensors measuring discharge.
For the type "Change In Sensor Value" the sensor specified under ID1 must
be a sensor measuring level. The threshold will correspond to a speed of
change (dH/dt).
For the type "Accumulated Value" the sensor specified under ID1 must be a
sensor measuring mass flux.
Action ID
Each control action needs a unique ID, which can be used to access the
actions from other dialogs.
179
MOUSE Control module
Setpoint Sensor
For actions of PID-type, it must be specified where in the system the set-point
is located. For this purpose a Sensor ID must be specified.
PID ID
Each set of PID parameter sets is identified with a unique ID which can be
used to access the information from other dialogs.
Tuning of the PID constants (Ti, Td and Kd) is not a straightforward task.
Understanding of the theoretical background and the numerical solution of
the control equation would be beneficial in this process. The following may be
used as an elementary guideline.
-setpoint downstream:1.0
-setpoint upstream:-1.0
Weirs
-setpoint downstream:-1.0
-setpoint upstream:1.0
1.0
181
MOUSE Control module
Pumps:1.0
1.0
NOTE: The sign on the K-factor is very important - if it is wrong it will cause
the control function not to work at all since the device will typically move to
one of the extreme positions and stay there till the end of the simulation.
Figure 9.5, Figure 9.6 and Figure 9.7 show examples of how the actual varia-
ble (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 9.5 Fluctuations around the set point depending on the size of the propor-
tionality factor, K
Figure 9.6 Fluctuations around the set point depending on the size of the deriva-
tion time, TD.
Figure 9.7 Fluctuations around the set point depending on the size of the integra-
tion time, Ti
183
MOUSE Control module
For 'Orifice with gate' and 'Orifice with weir' you need to define an orifice first.
Then the 'Controllable devices' dialog is used for adding either a movable
gate or a movable weir to the orifice.
For 'Orifice with weir' you only need to make the selection in this dialog. You
should not add the weir in the 'Weir' dialog (please refer later section on the
difference between an ‘ordinary’ weir and a weir in a orifice).
When the 'Device type' has been selected the specific parameters for the
configuration of the device can be entered. For this purpose the dialog has
three separate sections. The section to the left covers pumps while the sec-
tion to the right covers valves.The center section covers three different device
types: 'Weir', 'Orifice with gate' and 'Orifice with weir'. See Figure 9.10.
If e.g. the device type is 'Weir' then you have access to the 'Weir ID' field for
selecting the weir to be controlled.
For the device types 'Orifice with gate' and 'Orifice with weir' you will have
access to the 'Orifice ID' field where you can select the orifice for which you
want to add either a gate or a weir as defined by the device type selection.
The input for each for these types is described in later sections.
185
MOUSE Control module
9.8.1 Pumps
9.8.2 Weirs
Max Level
The maximum elevation of the movable weir crest.
Min Level
The minimum elevation of the movable weir crest.
The fixed weir crest level is not used for a RTC weir.
Max Speed Up
The maximum velocity for movement of the weir in upward direction.
If the device type is either ‘Orifice with gate’ or ‘Orifice with weir’ the mid-sec-
tion of the dialog becomes active.
187
MOUSE Control module
Max Level
The maximum elevation of the movable gate/weir.
Min Level
The minimum elevation of the movable gate/weir.
Max Speed Up
The maximum velocity for movement of the gate/weir in upward direction.
It is possible to define weirs in two different ways. Either through the 'Weir'
dialog alone which we will refer to as an 'ordinary weir'. The second way is to
use the combination of an orifice and the controllable device type 'Orifice with
weir'.
Figure 9.14 Difference between a ordinary weir and an orifice with a weir.
The main difference is that the 'Orifice with weir' can close the orifice com-
pletely while the ordinary weir is always open upwards. It has no ceiling and
in principle the flow can always pass over the weir if the water level is higher
than the weir crest level.
For the 'Orifice with weir' once the weir is fully closing the orifice then no flow
will pass the weir even if the water level is above the crest level.
189
MOUSE Control module
9.8.5 Valves
Max Opening
The maximum opening of the valve in percentage.
Min Opening
The minimum opening of the valve in percentage.
Max Speed
The maximum velocity for movement of the valve in percentages/sec.
Any number of rules can be specified to control the device. The statements
are evaluated sequentially starting from the top. This means that appropriate
sequence of rules is essential for the achievement of the desired control
logic. Evaluation of a logical condition belonging to a rule as "TRUE", leads to
the selection of the specified control action. If a logical conditions is "FALSE",
the evaluation proceeds to the next rule on the list.
The blocking time ensures that a certain rule is applied for a minimum period.
This means that the system is locked to certain rule in a period equal to the
blocking period after it was activated.
All control actions used to control a specific device must be of the same func-
tion type – corresponding to the specified Device Type and Control Type.
For PID-control, all control actions must refer to the same set-point sensor –
change of the set-point sensor during simulation is not allowed
For time dependent control (Input sensor of type “Time since start of simula-
tion”) only ONE rule can be specified. The menu allows more functions to be
inserted, but only the first one will be used during the simulation.
When starting a simulation, the system checks if these conditions are fulfilled
and in the case of any violation, the simulation will not start.
191
MOUSE Control module
Figure 9.16 The control simulation is activated on the network parameters tab page
The user can execute code at 3 different steps in the code. The Initialization
step is for executing code opening, preparing or reading files needed in the
simulation. The Time step is the step in which the control strategies are run.
The End step is the step in which files opened in the Initialization step can be
closed. The three steps are defined under the mouse604.exe Simulation |
Options dialog, where a filename and a procedure name are specified for
each step under UWC. The compiled DLL with the UWC must be located in
the same folder as the mouse_hd.dll file.
There are currently over 140 different procedures and functions exported
from MOUSE ranging from setting PID constants runtime to retrieving the
flow in a specific grid point in a link. Using the methods require the user to
include the MOUSEDLL.pas file into the DLL Delphi project. An example of
the use of this is listed below.
library Test;
uses
SysUtils,
Classes,
{$R *.res}
var
OutFile : TextFile;
procedure InitStep;
begin
AssignFile(OutFile, 'C:\MOUSEOutput.txt');
Rewrite(OutFile);
end;
procedure TimeStep;
var
NodeLevelString : String;
193
MOUSE Control module
MOUSETimeString : String;
begin
NodeLevelString := Floattostr(GetNodeLevel('Node_11'));
end;
procedure EndStep;
begin
CloseFile(OutFile);
end;
exports
InitStep,
TimeStep,
EndStep;
end.
The example above simply writes the water level in a node and the corre-
sponding time to a file, but the possibilities are endless.
The next example is the RTCExample provided in the installation where the
actual control is handled from UWC. The RTC control in the UND file is not
changed for this, which means the UWC RTC overrides the RTC in the UND.
This provides a way of overriding standard RTC in e.g. emergency situations
in the network. Alternatively the UWC RTC can be used alone. In this case
the Control Functions and Logical Conditions are not needed for the Con-
trolled Devices, since the actual control is handled externally. In addition to
this, the control type number for the 3 PID controlled devices is changed to 5.
This change is done using a text editor.
library RTCExample;
uses
SysUtils,
Classes,
MOUSEDLL;
{$R *.res}
procedure TimeStep;
begin
LevelA := GetSensorValue('Level_A');
Level02 := GetSensorValue('Level_02');
Level12 := GetSensorValue('Level_12');
Level22 := GetSensorValue('Level_22');
195
MOUSE Control module
//Pump_2_3
SetPoint1 := -0.6*LevelA+67;
SetPumpPIDSetpoint('Pump_2_3', SetPoint1);
Measured1 := Level02;
SetPumpPIDMeasurement('Pump_2_3', Measured1);
//Weir_12_3
SetPoint2 := -0.6*LevelA+67;
SetWeirPIDSetpoint('Weir_12_3', SetPoint2);
Measured2 := Level12;
SetWeirPIDMeasurement('Weir_12_3', Measured2);
//Orifice_22_3
SetPoint3 := -0.6*LevelA+67;
SetGatePIDSetpoint('Orifice_22_3', SetPoint3);
Measured3 := Level22;
SetGatePIDMeasurement('Orifice_22_3', Measured3);
//Pump_2a_3
StartLevel := LevelA+0.2;
StopLevel := LevelA;
SetPumpStartLevel('Pump_2a_3', StartLevel);
SetPumpStopLevel('Pump_2a_3', StopLevel);
//Weir_12a_3
WeirPosition := LevelA;
SetWeirPosition('Weir_12a_3', WeirPosition);
//Orifice_22a_3
GatePosition := LevelA;
SetGatePosition('Orifice_22a_3', GatePosition);
end;
exports
TimeStep;
begin
end.
197
MOUSE Control module
Based on the specified criteria, the system automatically selects the periods
for dynamic pipe flow simulations during wet weather. This reduces the simu-
lation time significantly, while preserving the accuracy in the simulation of rel-
evant dynamic effects in the system - CSOs, surcharges, etc. For example,
by running successive LTS simulations with the existing drainage system
configuration and the planned upgrade, impacts of the planned investments
(e.g. new sewers, retention tanks, control schemes) on the system perfor-
mance can be reliably tested and evaluated based on the computed opera-
tional statistics. This allows the user to develop the optimal rehabilitation /
upgrade strategy, e.g. for satisfying the requirements of environmental regu-
lations.
199
Long term statistics
However, use of dynamic simulation for the hydraulic and pollution transport
analysis of drainage networks is only beneficial in periods of high-flow
dynamics. Otherwise, simpler and much faster simulation tools can be used
with a similar level of accuracy. Since high-flow dynamics in urban drainage
networks are normally associated with rainfall events, the dynamically simu-
lated periods can be reduced to the wet periods, i.e. to only a few percent of
the certain historical period.
The fundamental principle behind the LTS concept is therefore the ability to
reduce the simulation time without losing significant information. LTS is able,
based on the user-specified criteria, to eliminate the irrelevant periods from
the simulated series, and to limit the dynamic simulations to a series of dis-
connected events, instead of the entire simulation period
For the purpose of isolating the intervals relevant for dynamic simulation, LTS
generates a job list on the basis of system information, the knowledge of time
series for dry weather flow (i.e. wastewater), rainfall and other loads and the
set of specified 'Job-List Criteria'. The simulation jobs from the list, limited by
the simulation start time end the simulation stop time, may optionally be sub-
ject to further evaluation by a set of 'Runtime Criteria' during the simulation.
Runtime Criteria are founded on the evaluation of the operational variables
within the network itself which can potentially extend the simulations beyond
the time intervals defined in the job list. Since the 'runtime criteria' are evalu-
ated during the simulation, it is not possible to determine the exact duration of
the dynamic simulation in advance.
Upon its creation, the job list can be edited, so that the actual start time and
date and the estimated stop time for the individual jobs may be modified
according to specific requirements
Job list criteria are used to select the intervals which (at least) should be cal-
culated with the hydrodynamic model during a computation with LTS.
The job list criteria represent actually various thresholds for the system load-
ing, which are used to evaluate if the state of the system is within or outside
the "dynamic" conditions. There are separate START criteria, which deter-
mine the start dates and times for the dynamic simulations, and a single
STOP criterion, which is used for the estimation of the ending dates and
times for dynamic simulations.
The START criteria are evaluated against the loading time series - dry
weather flow, rainfall events, etc., for the system as a whole, for a part of the
system, or for an individual sub-catchment. A number of START criteria can
be defined, and if at least one of them is fulfilled, a certain "event" will be
included into the job list.
The STOP criterion includes a threshold for the total load for the entire sys-
tem, and an additional time condition which extends the simulation for the
specified time after the loading has fallen below the threshold. The objective
of the STOP criterion is to delimit the dynamic simulation, but without a poten-
tial loss of important information due to the time offset between the loads and
effects in the sewer network.
201
Long term statistics
Table 10.1
Appropriate initial conditions for the individual jobs are essential for achieving
realistic statistics in hydrodynamic LTS simulations.
LTS provides two different methods for initializing the system. The system
can either be initialized in an empty state before each job (default behaviour),
or a set of hotstart files can be provided. If the latter option is used, then the
specific hotstart file used for a job is selected based on the inflow to the sys-
tem at the start of the job. A detailed description of the method is presented
below.
If the first option is selected (empty initial system), then the pipe system is
emptied of water before each job in the LTS job list is started. Water levels in
outlets will still be extended backward into the system similar to a 'normal'
simulation without hotstart files (see MOUSE Pipe Flow reference manual).
This option is valid only for storm drainage systems where there is no water in
the system during dry periods. It is not recommended to use this option in any
other situation.
If a set of hotstart files is provided, then each hotstart file must be provided
with an inflow interval. This interval is the "validity interval", meaning that the
specified hotstart file is valid (the file will be used) for all jobs in the job list
where the inflow to the network at the beginning of the job is in that interval.
Note that the inflow to the system when selecting a hotstart file is calculated
from inflows to the system that are marked "runoff" or "other". Inflows marked
as "dry weather flow" (DWF) are excluded from the total inflow when select-
ing a hotstart file. Inflow from the rainfall runoff simulation is considered "run-
off" and therefore included in the total inflow calculation.
If the total inflow to the system at the beginning of job is outside all the validity
intervals provided for the hotstart files, then the system will be initialized with
empty pipes before the job is started.
In order to take daily variations in inflow into account, the time-of-day at which
the values from the hotstart file should be used is set equal to the time-of-day
when the job starts. The date however, is provided by the user. This means
that care must be taken to eliminate the filling phase in the beginning of the
hotstart files. This is achieved by running simulation that produce the hotstart
files over two or more days before the date provided as hotstart date. It also
means that the hotstart file must cover a period of at least 24 hours. If all 24
hours of the hotstart date are not covered by the hotstart file, then the simula-
tion might fail.
This method of using the starting time of the job to account for daily variations
has the drawback that weekly and seasonal variations cannot be taken into
account.
203
Long term statistics
Total inflow Defines a lower threshold for inter- Yes, if Hot- InitFrom
interval - vals for total inflow (excluding start
From DWF) which is used at the start of
each job in the job list to evaluate
what hotstart parameters to use.
Total inflow Defines an upper threshold for Yes, if Hot- InitTo
interval - To intervals for total inflow (excluding start
DWF) which is used at the start of
each job in the job list to evaluate
what the hotstart parameters to
use
On the basis of user-specified Job List criteria, the simulation input files, the
simulation period (specified in the Start Computation dialog) and optionally,
the initiation hot start files, a job list is automatically created.
The simulation periods for individual events represent the minimum simula-
tion time, which may be extended in run-time according to the fulfilment of
run-time criteria.
The job list is written into an ASCII file (*.MJL, written in the 'PFS' format),
which contains all relevant information about input files for the simulation,
definition of the time periods to be simulated (for the individual jobs), and
optionally information about the hot start files and hotstart date/time. An
example of a job list file is presented in Figure 10.3
205
Long term statistics
Optionally, the job list file can be edited before starting the LTS calculation.
There is no dedicated dialog available for editing the job list.. Instead, Win-
dows Notepad is used. Alternatively, any ASCII editor can be used. When
editing the job list before starting the computation it is also possible to change
the minimum and maximum default time step for each job. Per default, each
job has assigned the same minimum and maximum time step, according to
the global definition. This ca be changed by editing the values for Dt_MaxSec
and Dt_MinSec (see below).
[SIMULATION_EVENT]
Dt_MaxSec = 60
Dt_MinSec = 10
EndSect // SIMULATION_EVENT
At every time step during computation (but only after the end time in the job
list), each criterion is subject to logical evaluation as an element of an evalua-
tion matrix. The evaluation matrix consists of an arbitrary number of lines,
207
Long term statistics
each containing up to five individual criteria. Each line of the evaluation matrix
represents a combination of criteria which have to be fulfilled simultaneously
(i.e. criteria are linked with the logical operator AND) if the combination is to
be evaluated as TRUE and the simulation stopped.
If one line has been evaluated as FALSE, the evaluation algorithm jumps to
the next line. This means that the successive lines are connected by the logi-
cal operator OR. The process continues until one line has been evaluated as
TRUE (the simulation stops).
If none of the lines has been evaluated as TRUE, the simulation continues for
one time step. Then, a new evaluation cycle is initiated and the entire evalua-
tion repeats.
If no evaluation matrix is specified, periods specified in the job list will be sim-
ulated exactly, i.e. no runtime criteria will be evaluated.
Table 10.3
Table 10.3
Edit field Description Used or Field name
required by in data struc-
simulation ture
Threshold The residual value of the variable Yes StopValue
defined by the ‘Type’ that must be
achieved in order to evaluate the
criterion as true
Duration A period in which the parameter Yes StopTime
value must be below the threshold
in order to evaluate the criterion as
true.
The specified runtime criteria are saved in the MOUSE LTS (*.MTF) file.
209
Long term statistics
Figure 10.6
Switching the simulation mode to 'LTS' opens a new part of the dialog, includ-
ing fields for selecting, generating and editing a job list file.
Prior to the computation, a job list must be generated. Several different job
lists can be generated and the appropriate one selected from the list.
While generating a job list, MOUSE scans the boundary inputs, and evalu-
ates the job list start criteria in order to identify the events to be simulated.
Each event is preliminarily delimited by the evaluation of the stop criterion.
Simultaneously, total loads (excluding DWF) at the simulation start for each
event are determined, and the appropriate initialization file is selected. Please
note that MOUSE LTS engine only handles a unique AD component. If sev-
eral AD components are specified then MIKE 1D engine has to be used.
When using a job list in a LTS simulation, the specified simulation period may
include the entire job list or only a small fraction - ultimately only a single
event. This feature opens the possibility to use the same job list through dif-
ferent phases of the project - initially using only individual events for test sim-
ulations (or in the calibration process), and later extending the simulation to
the entire job list. The advantage of using the job list instead of an "ordinary"
continuous mode simulation is access to useful statistics available only with
MOUSE LTS.
The user-specified result file contains only selected time series of particular
interest. This makes it especially practical for keeping the long-term results
sufficiently resolved, and the standard result file may often be completely
abandoned. Detailed information on the specification of the contents of user-
specified result files is provided in the Chapter entitled ‘MOUSE Simulations’.
The specification of the contents of the user-specified result file for hydraulic
variables (i.e. *.PRF file) applies also for the results of pollution transport
model (AD), i.e. defines the user-specified *.TRF file. This practically means
that at the specified locations, all computed water quality variables will be
saved. Selective saving of AD variables can be further controlled through the
DHIAPP.INI file, located in MIKE URBAN/Bin directory (section MouseTRAP,
see below).
[MouseTRAP]
211
Long term statistics
PRINT_WEIR_AD=1
PRINT_PUMP_AD=1
PRINT_NODE_AD=1
PRINT_GRID_AD=1
PRINT_GRID_MASS_AD=1
PRINT_MASS_BALANCE_AD=1
If the parameter value is changed to 0 (zero), then the item will not be written
in the result file.
When running any simulation over a longer period, focus of interest is usually
turned to the operational statistics rather than to the raw time series. There-
fore, MOUSE LTS provides an option to save an additional result file (*.ERF)
containing various statistics (event-based and annual) over the individual
computed time series and the system performance as the whole.
The main parameters which control the event definition for various types of
statistics are specified in the editor 'Statistics Specification - Global' (see Fig-
ure….).
Also, this editor allows for toggling between the annual and monthly statistics
and for requesting the statistics on system emissions (flows and pollutants
out of the system).
The peak values statistic operates with instantaneous peak values separated
by at least dT. If,for instance, in peak value statistics a continuous series of
water levels includes 2 peaks within an interval shorter than that specified for
Max(H) events, only the larger of the two values will be counted in the peaks
statistics. Accumulated discharge and duration statistics is associated with
discontinuous variables like pump or weir discharges (i.e. volumes). The sta-
tistics are performed on the events delimited with zero flows. Thus, two suc-
cessive weir or pump flows are considered as independent only if they are
separated by an inactive interval (i.e. with zero flow) of at least the specified
length dT. On the contrary, if the two flow occurrences are close to each
other, the statistical parameters will be computed from the sum of the two.
For the continuous flow variables (e.g. pipe flow), which in principle never
drop to zero within the simulation period, it means practically that the accu-
mulated discharge and duration statistics are calculated for the entire length
of the simulated event.
Both for peak statistics and for accumulated flow statistics, two events are
always considered as independent if the computation has been stopped in
between by the specified stop criteria. Thus, specification of a very long dT -
longer than any individual event to be simulated, would result in the number
of peak values and the number of accumulated discharges and duration cor-
responding exactly to the number of simulated events.
213
Long term statistics
Max levels:
Total discharge:
Saves ranked total event-accumulated discharges for weirs, pumps and out-
lets.
Saves total accumulated yearly/monthly discharge for weirs, pumps and out-
lets in chronological order.
Duration discharge:
Saves ranked total accumulated event discharge duration for weirs, pumps
and outlets.
Maximum Concentration:
Event Loads:
Saves ranked total event-accumulated mass flux for weirs, pumps and outlets
for all simulated pollutants.
Annual Loads:
Saves total accumulated yearly/monthly mass flux for weirs, pumps and out-
lets for all simulated pollutants, in chronological order.
Max levels:
Max flow:
Max velocity:
Location
For all result types it must be specified for which location(s) the desired
results are to be saved. Location type can be either 'General' (the whole sys-
tem), 'List' (a number of nodes/links specified on NSE/LSE/MUS-file), or 'Indi-
vidual' (one node/link name). If 'List' or 'Individual' is selected, the
supplementary information of selection file name or node name should be
specified.
215
Long term statistics
Save Type
For links, it must be specified where in the link data should be saved. There
are four saving possibilities:
The maximum number of events must be specified within the rand 1-5000.
In both cases the design module will be an excellent help to speed up the
process by designing the exact pipe sizes which is needed to fulfil the design
criteria and to finally verify that the design system is correctly designed.
In the present version the pipes are always designed from the principle of
increasing the diameter until the maximum level upstream in the system is
below the design criteria. Thus the module cannot be used to design pipes
from downstream criteria.
After end of the simulation the maximum Manning number is used to calcu-
late the new diameter which will then be used in next iteration.
This means that after a number of iteration - depending on how good the ini-
tial choice of diameters was - the module will stabilise results on the diame-
ters needed to fulfil the design criteria.
Each pipe (group) is designed independently of the others, which means that
the process can be long if there are many pipes, since changing of one diam-
eter can often influence the design of both upstream and downstream pipes.
217
Automatic pipe design with MOUSE
The ADP file can either be written directly in a text editor, e.g. Notepad, or a
wizard like tool can be utilized for writing the ADP-file for the design module
(see further in 11.2.7). The format of the DESIGN section in the ADP file is as
follows:
219
Automatic pipe design with MOUSE
[DESIGN]
DESIGNRESULTFILE = ‘MyOwnResults.txt’
[Commercial_Pipes]
Diameter = 50
Diameter = 100
Diameter = 200
EndSect // Commercial_Pipes
EndSect // DESIGN
There are two design types available for designing pipes – fixed invert levels
of pipes and fixed top levels of pipes.
When fixed invert levels are selected the network geometry is updated only
with pipe diameters hence all invert levels are kept as in the original setup.
When fixed top levels are selected the top of each pipe is fixed, which
means that all invert levels need to be updated according to the new diame-
ters. All adjacent pipes are unchanged, hence when a nodes invert level is
changed any attached pipe which was previously connected at the invert
level of the manhole is now connected above invert of the manhole - at the
original invert level of the manhole. A node invert level is never changed to a
higher level than for the original setup.
Note: In the present version fixed top levels can only be used for pipes con-
nected to manholes in both ends. If either upstream of downstream node is
specified as a basin, the pipe will be designed according to the principle of
fixed invert levels, regardless of design type. This is necessary in order to
accommodate for the complexity of changed basin geometry, weir levels,
pump start-/stop levels etc. which will often be present in basin structures.
The design tool accommodates two types of design criteria – critical levels
and filling degree. When the Critical Levels are applied, the pipes will be
designed according to the critical level specified in the MOUSE nodes dialog.
The Filling Degree is in this case interpreted as the maximum water level in
upstream node compared to the upstream pipe levels. E.g. if the filling degree
is set to 0.8 the design criteria will be set as a water level corresponding to
the pipe being filled 80% at the design node.
If more than one pipe is connected to the design node the lowest value for all
connected pipes is applied in order to ensure that the criteria is fulfilled for all
pipes.
If the pipes are designed as a group their sizes will decrease or increase with
the same factor in order to obtain the design criteria for one specified
upstream node. This means that the pipe sizes will change with same per-
centage for all pipes.
When the design acts on individual pipes, each pipe will initially be automat-
ically assigned to a design node corresponding to their immediate upstream
node. Here the upstream node is defined as the node connected to the high-
est end of the pipe from geometric point of view.
If more than one pipe is connected to the same upstream node they will auto-
matically be grouped together and thereby maintain their relative sizes.
During the designing simulation there will be long periods where the water
levels are below the design level. In order to avoid that these periods influ-
ence too much on the new pipe diameters a lower limit for designing is speci-
fied. Hence when the water level is below this limit the design module will not
try to resize the current pipe size until the level rises above the limit again.
The limit is defined relative to water level of the design criteria. Thus a nega-
tive value indicate a lower limit below the design level and vice versa for a
positive value.
221
Automatic pipe design with MOUSE
Through the iterations the pipe design will suggest diameters without taking
into account that only some diameters are actually manufactured. Specifying
the available pipe diameters can be used to let the module select the most
appropriate standard pipe after end of simulation.
The standard pipe diameters should be given in increasing order and in milli-
metres (or inches for US-units).
The ADP file can either be written directly in a text editor, e.g. Notepad, or a
wizardlike tool can be utilized for writting the ADP-file for the design module.
This tool can be activated under Tools|User Tools as seen in Figure 11.2.
Figure 11.2 The tool for creating the ADP-file is now located under User Tools
When the tool is launched you are guided step-by-step setting up the ADP-
file.
In the simulation launcher two new item types can be selected from the dialog
<Charts>+<Modify>: Designs and Profiles.
Figure 11.4 For changing the charts displayed during simulation go to Charts|Modify
“Designs” are time series plots for the nodes where design criteria has
been specified. The default selection will include actual level, critical
level and calculated factor of the manning number which is changed dur-
ing the simulation
“Profiles” are longitudinal profiles. All profiles which are saved as LPF
files can be selected
223
Automatic pipe design with MOUSE
During the simulation it will hence be possible to see the results while simu-
lated, see Figure 11.6.
Figure 11.6 The results are displayed during the simulation. The critical levels used
for the design are clearly marked on the profiles.
The additional dialog for pipe design is dominated by a grid with input param-
eters and results from the design, and can be used to adjust parameters for
the design between the iterations.
ID ID of the pipe
UpStrNodeID ID of the node where design criteria are evaluated
Init Factor Initial factor of Manning number for next iteration - val-
ues smaller than 1.0 will allow the pipe to decrease in
size if necessary, however too small values can cause
numerical problems. Values bigger than 1.0 are usu-
ally not relevant since it will cause the pipe to increase
in size even if it's not necessary. Default is 0.98.
225
Automatic pipe design with MOUSE
If changes are made in the grid these must be saved in order to be used for
the next iteration. Saving can be done by pressing "Apply changes for next
iteration", but if data has been changed and a new iteration started, a mes-
sage will occur asking if the changes should be used or lost.
Finally there is the possibility to "Create MEX-file With New Diameters" which
will change the parameters for the specified elements and create a MEX file.
This MEX file can be used to import and update the elements in MIKE
URBAN once the optimal set of values has been identified. In order to import
and update the elements in the network based on this MEX file you need to
go to File|Import/Export and switch to the ‘Use Advanced’ mode. Under the
list of pre-defined import jobs you will find the job ‘CS-Update From Pipe
Design’. Choosing this option and the MEX file in question (see example in
Figure 11.8) will update the elements in your MIKE URBAN model after
pressing ‘OK’.
The option "Delete Design File" which will delete results of all the performed
iterations and cause the design to start from scratch.
Figure 11.8 Choosing to import and update the elements based on the result of the
pipe design
227
Automatic pipe design with MOUSE
The primary role of the Surface Runoff Quality (SRQ) is to provide a physi-
cally-based description of the relevant processes associated with sediments
and pollutants due to surface runoff, and then provide surface runoff sedi-
ment and pollutant data for the other pipe sewer network sediment transport
and water quality modules. The following processes can be accounted for:
Sediment deposits can greatly reduce the hydraulic capacity of sewer pipes
by restricting their flow area and increasing the bed friction resistance. Pipe
Sediment Transport (ST) can be used to account for these problems, by sim-
ulating pipe sewer network sediment transport-including deposition and ero-
sion from non-uniform (graded) sediments. Contributions from rainstorm
wash-off and dry-weather wastewater flow can be included. The ST process
runs in conjunction with the dynamic flow routing, thereby simulating dynamic
deposition of sediment and providing feedback due to the change in pipe
area and resistance caused by sediment deposition. The following issues can
be addressed:
229
Modelling Water Quality with MOUSE
SRQ can be used by itself to study water quality and sediment transport pro-
cesses on catchments. The result can be used as input to a river model
(MIKE11) or to other water quality modules which deal with processes in the
pipe model.
Finally, if the sediment transport is the only modelling purpose the ST module
can be used separately.
Problem types and processes shows which water quality modules should be
used for a list of problem types. Figure 12.1 gives an overview of the required
modules for different pollutant and sediment types.
231
Modelling Water Quality with MOUSE
The model works with two sediment fractions, a fine and a coarse, each frac-
tion characterised by its mean diameter. The fine fraction typically has a parti-
cle size less than 0.1 mm and is limited in supply. On the contrary, the coarse
fraction is typically unlimited in supply. Thus, for the coarse fraction, transport
capacity is always satisfied and, therefore, the build-up description is only
applicable to the fine sediment fraction to determine the mass available for
transport.
The data for the surface sediment transport calculation is entered and edited
through dedicated editors accessible under the 'MOUSE|Water Quality'-
menu. The data for the surface sediment model consists of data for the Build-
up/Wash-off model, data for the Sediment and Pollutant model and data for
the Gully Pot model.
The parameters for the build-up/wash-off model can be specified in the dialog
displayed in the SRQ Main Parameters dialog Figure 12.2.The 'Build Up
Method' defines build up method for the fine sediment on the surface. A linear
build-up function, or an exponential function, can be chosen.
233
Modelling Water Quality with MOUSE
Notice that if the load generated from the washoff simulation subsequently is
supposed to be used as input in a pipe model, then it is imperative that the
"Save as concentration" is left un-checked.
Data for description of the build-up of pollutants for each component must be
specified. During dry weather flow, the concentration of dissolved pollutants
will build-up according to a user-specified linear build-up function, with a max-
imum threshold value. During storms, all sediment is routed straight through
the gully pots and a simple mixing model is used to mix the incoming water
with the gully pot liquor.
The data for the gully pots are specified and edited under 'MOUSE|Water
Quality| SRQ Gully Pot Data - see Figure 12.4.
235
Modelling Water Quality with MOUSE
Please note that when using MOUSE LTS engine only a single AD compo-
nent should be specified. If several AD components have to be included in
the simulation then use MIKE 1D engine.
237
Modelling Water Quality with MOUSE
The initial conditions are specified for individual nodes. The initial conditions
in the connected conduits is calculated by linear interpolation of the concen-
trations specified in the upstream and downstream nodes.
dC
-------- = KC (12.1)
dt
where:
C = the concentration
The decay constant is defined as an uniform decay over the entire model.
The AD model can be run with the components specifications only. In this
case all model specific parameters (decay constant, dispersion coefficient,
initial concentration) as well as boundary conditions are set to zero.
Please not that the sum of the two blend components must always add up to
100.
b
D = au (12.2)
where:
b = a dimensionless exponent.
If the exponent is set equal to zero, then the dispersion coefficient is constant
and independent of the flow velocity. The unit for the dispersion factor will
then be m2/s. If the exponent is 1, i.e. the dispersion coefficient is a linear
function of the flow velocity, then the unit of the dispersion factor will be
meter, and the dispersion factor will in this case be equal to what is generally
termed the dispersivity. It is possible to specify values of the minimum and the
maximum dispersion coefficients, in order to limit the range of the dispersion
coefficient calculated during the simulation.
The dispersion coefficient can be given either globally or locally. The global
description will be used at all locations except for those pipes where local
conditions have been specified in the menu. In the example above the global
statement indicates that a constant dispersion coefficient of 2.00 m2/s is
applicable, and independent of the flow velocity (the exponent is zero).
For two conduits (see the example below), the dispersion coefficients speci-
fied locally, 'overrule' the global specification and prescribe a constant disper-
sion coefficient of 4 m2/s.
239
Modelling Water Quality with MOUSE
If the flow is directed out of the model, the concentrations or transport at the
boundary are computed from the concentrations or transport within the model
area. In this case the concentration or transport outside the model area is
assumed not to influence the conditions in the model area.
If a flow reversal takes place so that flow is directed from the receiving waters
into the model area, the specified time series of boundary concentrations or
transport are used and assumed to be unaffected by the previous outflow
from the model area. A parameter K-mix is used to ensure a smooth transi-
tion between the calculated concentration or transport at the boundary before
the flow reversal and the specified time series of boundary concentrations.
The description of K-mix is given in the "MOUSE Pollution Transport Refer-
ence" manual.
Group 1 includes the aspects of the classical BOD/DO problem, e.g. the deg-
radation of organic material and subsequent consumption of oxygen. The
oxygen demand can be specified as BOD or as COD. The processes
included in the simulation of these parameters are the degradation of organic
material and the resulting oxygen consumption and the reaeration of oxygen.
Group 2 includes nutrients relevant for sewer systems, e.g. the reduced form
of nitrogen - ammonia/ammonium and phosphate. The nutrients are modeled
as conservative substances.
Group 3 is the heavy metals. The model does not a priori determines which
metals to be included. This is defined by the user. The Heavy metals are
modeled similarly to the nutrients, e.g. as conservative substances.
241
Modelling Water Quality with MOUSE
Group 4 is bacteria arising from human and animal wastes. Three groups of
bacteria can be modeled: faecal coliform, total coliform and streptococci. The
coliform bacteria are assumed to be subject to a decay (death) that depends
on the temperature.
Group 5 and 6 are facilities for specifying the parameters needed when pol-
lutants are modeled as sediment attached. This can be any pollutant (BOD,
heavy metals etc).
The determinants included in this group is the BOD/COD and dissolved oxy-
gen. The BOD is modeled as two fractions: the dissolved and readily degra-
dable BOD and the suspended BOD that is degraded slower and of which a
part is the heterotrophic bacteria responsible for the degradation of a part of
the dissolved BOD. The dissolved BOD can also in some cases besides the
degradation by the bacteria, be degraded in a biofilm at the inside of the
sewer pipe and at the water surface.
Figure 12.9 The processes included in the modelling of BOD/COD - dissolved oxy-
gen problems in sewer pipes modeled with MOUSE
243
Modelling Water Quality with MOUSE
The constants in this expression can be changed by the user (totally three
constants). Besides these constants a temperature coefficient can be speci-
fied.
The Group 4 of determinants is the bacteria that can be faecal and total coli-
form bacteria and streptococci. Most pathogenic microorganisms are usually
unable to multiply or survive for extensive periods in the aquatic environment.
Sedimentation, starvation, pH changes, temperature plus competition with
and predation by other naturally occurring microorganisms are factors
involved in the decay of pathogenic bacteria in the marine environment.
Escherichia coli (E. coli) is one of the dominant species in faeces from human
and warm-blooded animals. The organism itself is normally considered non-
pathogenic, but is very often used as indicator organism for faecal pollution
and hence a potential of real pathogenic organisms (e.g. other bacteria and
vira). Infections obtained from bathing waters are most often caused by vira,
protozoas, Shigella and Campylobacter.
Enteric bacteria die-off can be modeled very well by a first order decay reac-
tion. However, the die-off constant or decay rate is highly variable due to
interaction by environmental factors on bacterial die-off. The main factors are
suspected to be: light, temperature and salinity variations. In sewers, the only
relevant factor among these to include is the temperature. Other less impor-
tant factors are, e.g. predation,sedimentation and the amount of available
nutrients.
For each bacteria group, a first order decay rate at 20ºC has to be specified
together with a temperature coefficient.
MIKE ECO Lab is a highly flexible and open framework for formulation of
water quality models. MIKE ECO Lab utilises a concept of templates where
water quality models are defined transparently: Complete contents of the pre-
defined templates supplied by the MU installation is accessible for review and
can be modified without any need to involve the software vendor. Also, com-
pletely new templates can be developed by the user.
The MIKE URBAN installation includes two pre-defined water quality tem-
plates that can be applied “as is” - or they can be adjusted by the user to con-
form to the specific project requirements.
A MIKE ECO Lab template is as ASCII file which can be accessed and edited
by the MIKE ECO Lab editor. Note that MIKE ECO Lab is part of MIKE Zero
and it is necessary install MIKE Zero and acquire a license for MIKE ECO
Lab in order to use MIKE ECO Lab with MIKE URBAN.
Please consult the MIKE Zero ECO Lab documentation on how to create or
modify model templates, i.e. how to work with the MIKE ECO Lab editor.
Please also note that MIKE ECO Lab couplings are only possible with the
MIKE 1D engine. The older MOUSE engine does not support MIKE ECO
Lab couplings.
245
Modelling Water Quality with MOUSE
A simulation job with the coupled HD/AD and MIKE ECO Lab model is set up
and launched from within MIKE URBAN in a usual way.
MIKE ECO Lab in MIKE URBAN is supported by a set of four MIKE ECO Lab
editors in MIKE URBAN:
These four editors and the database tables behind them are described below.
This editor is used for loading (import) MIKE ECO Lab templates into MIKE
URBAN and assigning them to the model elements – nodes and links. The
data for this editor is stored in the msm_EcoLabTemplate MU database table.
A template can be assigned to the entire model (i.e. to all points where
HD/AD calculation is performed or to an arbitrary sub-set. A sub-set is
defined by a MIKE URBAN Selection file (*.mus) or as an individual element
(a node or a link).
Note that there can only be one MIKE ECO Lab template on a node or link. If
a template is used for All, then a template on a single or multiple nodes/links
will override the template on all elements.
Table 12.2 Field description for MIKE ECO Lab Templates editor
Pressing <Import> button after the template has been already loaded (im-
ported) into MIKE URBAN will re-load the template and any associations of
the MIKE ECO Lab state variables with MIKE URBAN AD components, as
well as any modifications of the MIKE ECO Lab constants will be lost.
247
Modelling Water Quality with MOUSE
This editor is used for coupling MIKE URBAN AD components to MIKE ECO
Lab state variables and for possible reconciliation of differences in units used
in MIKE URBAN and in MIKE ECO Lab, respectively. The data for this editor
is stored in the msm_EcoLab.Component MU database table.
The import operation described in the MIKE ECO Lab Templates and Assign-
ments section will automatically create the required entries in this table. The
user must couple each of these MIKE ECO Lab state variables (components)
to the relevant AD component.
The "Conversion Factor" entry (default value 1, i.e. no conversion) can be set
different from 1.0 if the units for the AD components and MIKE ECO Lab state
variables are different. Please note that the AD components require a specific
unit while MIKE ECO Lab does not require a specific unit for a specific com-
ponent type.
Table 12.3 Field description for MIKE ECO Lab State Variables editor
Please note that the only active fields in this editor are those for coupling the
MU AD components with MIKE ECO Lab state variables, and for specification
of the conversion factor. The remaining fields are read-only.
Forcings are external variables associated with computational points (i.e. are
spatially fixed), constant or variable in time, which affect some of the pro-
cesses on the state variables.
Forcings can be “Built-In” or “User Defined”. “Built-In” means that forcings are
provided by the hydraulic model through coupling with MIKE ECO Lab. “User
De-fined” forcings must be provided by the user, either as constant values or
as time series (dfs0 file name and item in the file).
This editor is used for setting the User Defined forcings associated with the
loaded MIKE ECO Lab templates. Built-In forcings are also imported into
249
Modelling Water Quality with MOUSE
The data table in MIKE URBAN database containing MIKE ECO Lab forcings
is msm_EcoLabForcing. The database records are created automatically dur-
ing the template import. The table contains as many forcings as there are in
the MIKE ECO Lab template.
Table 12.4 Field description for MIKE ECO Lab Forcings editor
Edit field label Database field Description
Forcing ID msm_EcoLabForcing.Forc- MIKE ECO Lab forcing identi-
ingID fier
MIKE ECO Lab msm_EcoLabForcing.Tem- MIKE URBAN name of the
Template pID template the forcing belongs to
ECO Lab Unit msm_EcoLabForcing.Eco- Unit for the MIKE ECO Lab
LabUnit forcing, as defined in the tem-
plate
Value msm_EcoLabForc- Constant value for the forcing
ing.DefaultValue
File msm_EcoLabForcing.File- Path/filename for dfs0 time
Name series file containing forcing
TS data
Item msm_EcoLabForcing.Item- dfs0 item reference
Label
Please note that the only active fields are those for the specification of the
forcing’s constant value or the reference to dfs0 TS file. The remaining fields
are read-only.
This editor is used for setting the values of constants used in the MIKE ECO
Lab template, typically in relation to the model calibration.
The data table in MIKE URBAN database containing MIKE ECO Lab con-
stants is msm_EcoLabCoeff. The database records are created automatically
during the template import. The table contains as many constants as there
are in the MIKE ECO Lab template.
251
Modelling Water Quality with MOUSE
Table 12.5 Field description for MIKE ECO Lab Constants editor
Please note that the only active field in the one containing the constant value.
The remaining fields are read-only.
To run a MIKE ECO Lab simulation, the network model must be prepared in
consistency with the MIKE ECO Lab template which is to be simulated. This
means that it should contain at least those AD components which are to be
coupled with MIKE ECO Lab state variables.
Further, the wanted MIKE ECO Lab template should be imported to MIKE
URBAN. Possible unit inconsistencies must be considered and resolved by
providing proper conversion factors. Values of the constant should be
reviewed and modified, as appropriate. User-defined forcings should be
specified, either as constants or as time series.
The network simulation with MIKE ECO Lab is activated by checking “Trans-
port(AD)” and “ECO Lab” checkboxes, as shown in Figure 12.15. Default
integration method (EULER) can be applied in most cases.
The simulation will result in three files containing all the hydrodynamic infor-
mation (*.res1d), advection and dispersion results for the state variables
(*AD.res1D) and the auxiliary variables, processes and derived outputs
included in the MIKE ECO Lab results file (ADEL.res1D).
explicit models,
morphological models.
The main difference between the two model types is that in the explicit mod-
els the bed levels are fixed during the simulation. In the morphological mod-
els the bed levels / sediment deposits are updated dynamically during the
simulation.
253
Modelling Water Quality with MOUSE
The morphological model for uniform sediment uses one grain diameter at
each point in the sewer system to calculate the sediment transport. In the
morphological model for non-uniform sediment it is possible to calculate the
sediment transport from a range of grain diameters. For the transport of fine
organic particles some of the non-uniform sediment fractions can be chosen
to be transported by use of the Advection-Dispersion transport formulation
instead of the non-cohesive sediment transport formula, as the transport of
particles transported in suspension is better described by used of the advec-
tion-dispersion equation. The fractions transported by the advection-disper-
sion module are still included in the morphological calculation, but the
description of erosion/deposition is changed.
Note, the morphological models require sediment or bed level boundary con-
ditions at all inflow boundaries. In order to chose the morphological sediment
transport model type, check the checkbox "Apply Morphological Modeling" on
the dialog 'MOUSE | Water Quality | ST Main Parameters'.
Engelund - Hansen;
Ackers - White;
Engelund - Fredsøe - Deigaard;
van Rijn.
Short descriptions of the sediment transport formulae are given below. No
general guidelines can be given for the preference of one formulation over
another, as the applicability of each depends on a number of factors.
All these formulae demonstrate that the sediment transport is a highly non-
linear function of the flow velocity: depending on the formulation, the sedi-
ment transport is proportional to the velocity raised to the power from 3 to 5.
Hence, the user should be very careful in the description of the hydrodynam-
ics in the model setup.
The suspended load is calculated as the depth integration of the local con-
centration and flow velocity. The method uses the reference concentration
computed from the bed load transport. The formula has been verified for par-
ticles in the range 200 - 2000 mm. The verification based on 600 data sets,
showed that 77% of the predicted bed load rates were within 0.5 and 2 times
the observed values, van Rijn (1984a). The verification for the suspended
load, using 800 data sets showed that 76% of the predicted values were
within 0.5 and 2 times the observed values, van Rijn (1984b).
255
Modelling Water Quality with MOUSE
The hydraulic resistance in the sewer originates from the pipe wall and from
the sediment deposits on the bottom of the sewer. The resistance from the
sediment deposits consists of two contributions, one part originates from the
grain friction and the other part originates from the expansion loss behind the
bed forms. The dimensions of the bed forms are determined by the sediment
transport and the flow, the resistance from the bed forms can be described
through sediment transport formulae.
The average shear stress in a pipe with sediment deposits is calculated from
the Einstein side-wall elimination procedure. The calculation is based on the
pipe roughness and the bed shear stress calculated from the sediment
deposits.
The data dialogs under the Menu option "MOUSE|Water Quality" contain
input parameters to the sediment transport model - MOUSE ST. The data are
e.g. sediment grain diameter, relative density of the sediment, critical bed
shear stress, type of sediment transport model, etc.
Selection of the sediment transport formulae and type of model takes place in
the Data Dialog "MOUSE| Water Quality| ST Main Parameters".
In a case when it has been chosen to calculate the bed shear stress during
the simulation, the updated values will be used in the hydraulic computations.
This implies that the Manning number specified for the bed may vary during
the simulation. A detailed description of how the shear stress is calculated is
given in the “MOUSE Pollution Transport Reference” manual. If it has been
chosen to update the bed level then a morphological model is run during the
simulation. Hence, boundary conditions in terms of time series of sediment
transport, bed levels and/or available fraction percentages, must be specified
at every point where inflow occurs in the sewer network. If it has been chosen
not to update the bed level, then an explicit sediment transport model is run
during the simulation and boundary conditions for sediment transport are not
required. The results of an explicit sediment transport calculation represent
the potential sediment transport.
The initial depth of the sediment deposits, the bed shear stress caused by the
sediment deposits and variation of the critical bed shear stress as a function
of the initial sediment depth.
The depth of the initial sediment deposits should be specified relative to the
bottom level of the pipe.
The bed shear stress from the sediment deposits can be specified in three
ways:
257
Modelling Water Quality with MOUSE
the bed shear stress is calculated from the local mean grain diameter
according to:
25,4
M = ----------------------- (12.3)
16
2,5d
the bed shear stress is the Manning's M taken from the dialog.
The variation of the critical bed shear stress can be given as a function of the
initial sediment depth as long as only erosion occurs. If deposition occurs, the
critical bed shear stress is set equal to the Shields critical bed shear stress.
This value of the critical bed shear stress will be used as long as the depth of
the sediment deposit is larger than the maximum depth to which the initial
sediment deposit has been eroded.
If the actual sediment depth > maximum eroded sediment depth then
c = c Shields (12.4)
else
y sediment EXP
C = bot + top – bot -------------------------- (12.5)
y initialdepth
where:
bot = the critical bed shear stress at the bottom of the pipe,
bot = the critical bed shear stress at top of the initial sediment deposit,
This dialog allows for the specification of the input data required for the simu-
lation of graded sediment transport and sediment sorting.
The 'Transport Mode' combo-box allows for the choice whether the sediment
fraction is transported with ordinary non-cohesive sediment transport formu-
lations, or the fraction is transported by the Advection-Dispersion formula-
tions. The transport of fine organic sediment fractions is better described by
use of an Advection-Dispersion transport formulation instead of the non-
cohesive sediment transport formulae. If it has been selected to transport the
fraction with the Advection-Dispersion model then the critical shear stress for
deposition cd, the critical shear stress for erosion ce and the erodibility coef-
ficient must be specified. The morphological changes resulting from the sedi-
ments transported by the Advection-Dispersion formulations are still included
in the calculation of the bed level, but the changes in bed levels are now cal-
culated from source/sink terms. For further information, refer to the “MOUSE
Pollution Transport Reference” manual.
259
Modelling Water Quality with MOUSE
Ws –n
E = 1 – 1 + -----------------------
- (12.6)
nx Q A
where:
E = efficiency
Ws = setting velocity
n = removal coefficient
261
Modelling Water Quality with MOUSE
Boundary conditions for the ST module can be specified via sediment trans-
port results from a surface calculation (.SMR-file) and/or via the boundary
connection system. It is possible to connect sediment loads or concentrations
to nodes in this system either as time series, constant values or cyclic pat-
terns. Boundary conditions must be specified at all inflow points in the sedi-
ment transport model. Hence, if a water level at an outlet results in inflow to
the model, a boundary condition must be specified at this location for the
morphological model. An appropriate boundary condition at such a location
will often be to specify the sediment transport equal to zero.
The main limitations to the time step in the sediment transport calculation are
the boundary conditions. Both the hydrodynamic and the morphological
boundary conditions should be well resolved. Finally the calculated changes
in the bed level should be sufficiently small in order not to generate instabili-
ties in the hydrodynamic calculations.
The functionality is aimed at modelling the pollutant loads from runoff and
infiltration water and is implemented in all the runoff models of the MOUSE
engine.
The concentrations of a given pollutant can be given for each specific catch-
ment. For one catchment it is possible to give a concentration on the slow
runoff component (SRC) and the fast runoff component (FRC) or override
these settings stating a concentration for the total runoff (TOT). Instead of
fixed concentrations it is possible to state a name of a tabular data set. The
tabular data set should then hold a table with concentrations versus normal-
ised flow, i.e. the runoff divided by the total catchment area.
263
Modelling Water Quality with MOUSE
The concentrations given as input data described above come from a consid-
eration of pollutant loads and land use types for the various catchments.
These can be difficult to assess. The requirements for a GIS analysis, using
MIKE URBAN, are thematic maps of different land-use types such as build-
ings, traffic area, industrial area and areas with dense traffic.The catchment
processing tool, MOUSE|Catchment Tools|Catchment Processing, can be
used to assign concentrations based on these thematic maps.
Figure 12.22 Using the catchment processing tool for assigning concentrations
Table 12.7 The edit fields for the surface runoff pollutants dialog
Edit field Description
Catchment ID The ID of the catchment. It is recommeended
to add pollutants only to catchments with a
surface runoff areas
Pollutant compo- The ID name of the pollutant component in the
nent AD model
Method Method = Cst concentration (method 1), Table
concentration (method 2), EMC formula
(method 3)
For all methods values in the fields for total
concentrations (TOT) overrides settings speci-
fied for SRC and FRC.
SRC concentration Constant concentration value for the slow run-
off component. Only used in RDI models.
FRC concentration Constant concentration value for the fast run-
off component. This is also the component
used in the model A, B, C and UHM.
TOT concentration Constant concentration value for the total run-
off component. Can be used in all models.
SRC dataset ID Only applied with method 2. ID identifier for the
tabular dataset applied to the slow runoff com-
ponent.
FRC dataset ID Only applied with method 2. ID identifier for the
tabular data set applied to the fast runoff com-
ponent
265
Modelling Water Quality with MOUSE
Table 12.7 The edit fields for the surface runoff pollutants dialog
Edit field Description
TOT dataset ID Only applied with method 2. ID identifier for the
tabular data set applied to the total runoff com-
ponent
K Only applied with method 3. K factor used in
the formula.
Event threshold Threshold that defines when a runoff event
starts and stops. Events with flows less than
the treshold are seen as dry weather when cal-
culating the antcendent dry weather period
(ADWP)
ADWP initial The initial ADWP used before the first runoff
event is detected in the data.
ADWP minimum If the duration between two events is smaller
than the ADWP minimum then the two events
are perceived as only one event.
Exponent Only applied with method 3. Exponent used in
the formula.
exp onent
EMC = SMC K ADWP (12.7)
input data and the antecedent dry weather period (ADWP) is determined from
the runoff from the given catchment. The method only applies for the FRC
runoff component.
Based on the layout of the basin or pond, the hydraulic load and the type of
pollutant the removal efficiency is computed and the corresponding amount
of pollutant is removed during the simulation with the MIKE URBAN Collec-
tion System Pollutant Transport module.
267
Modelling Water Quality with MOUSE
The treatment node relations are used from the treatment relationships, see
Figure 12.25.
269
Modelling Water Quality with MOUSE
2 2 2 2
R = A1 h + B1 H + C1 q + A2 h + B2 H + C2 q 2 + A3 h + B3 H + C3 q 2
where (12.8)
h:hydraulic efficiency (%)
3 3
q:volume load 0.001 m h m
271
Modelling Water Quality with MOUSE
Figure 12.26 The coefficients for a number of pollutants have been determined
through MIKE 21 simulations
To sum up the methodology, for each treatment node and for each pollutant:
273
Modelling Water Quality with MOUSE
275
276 MIKE URBAN - © DHI
Terms and Concept
A SWMM network within MIKE URBAN can be defined in one of the following
ways. Please note that MIKE URBAN supports the EPA SWMM5 format (if
you have a SWMM44 model you can locate conversion tools on the US EPA
website).
Import of external data (e.g. GIS) into MIKE URBAN CS SWMM network
Copying network data from MIKE URBAN CS Asset network into MIKE
URBAN CS SWMM network
Copying network data from MIKE URBAN CS MOUSE network into
MIKE URBAN CS SWMM network
Graphical data digitalisation and manual data typing within MIKE URBAN
The last option is frequently used in a combination with one of the previous
options as means for achieving a full consistency of the SWMM model.
277
Hydraulic Network Modelling with SWMM
Nodes
Conduits
Orifices
Pumps
Weirs
Outlets
13.2 Nodes
Nodes are hydraulic elements such as Junctions, Outfalls, Storage Units
and Flow Dividers used in modelling a CS SWMM storm and sewer collec-
tion system.
Outfalls are defined as terminal nodes of the conveyance system and used
to define final downstream boundaries under Dynamic Wave flow routing. For
other types of flow routing they behave as a junction. Only one link can be
incident on an outfall node.
Storage Units are defined as conveyance system nodes that provide storage
volume. Physically they could represent something as small as a catchbasin
or as large as a lake. They are allowed to have any shape as described by a
function or table of surface area versus height.
Flow Dividers are defined as conveyance system nodes used to split the
total flow between exactly to two outflow conduits in a prescribed manner.
Flow dividers are only active under Uniform or Kinematic Wave routing and
are treated as simple junctions under Dynamic Wave routing.
The CS SWMM Nodes editor organizes the related input data for the hydrau-
lic elements into the following groups:
279
Hydraulic Network Modelling with SWMM
Junction
Storage Unit
281
Hydraulic Network Modelling with SWMM
Outfall
FREE outfall bases the outfall stage on the smaller of the critical and nor-
mal depths of flow in the connecting conduit.
NORMAL outfall uses just the normal depth of flow.
FIXED outfall uses a constant water elevation at all times.
TIDAL outfall has outfall stage varying in a repeating fashion over a tidal
period.
TIMESERIES outfall allows the outfall stage varying in any specified
manner over time.
283
Hydraulic Network Modelling with SWMM
Flow Divider
Three different types of dividers can be specified for each divider node
OVERFLOW divider, all flow above the non-diverted conduits full flow
depth is diverted
CUTOFF divider, all flow above the cutoff level is diverted.
TABULAR divider, a curve is supplied that specifies the amount of
diverted flow to total flow.
WEIR divider the diverted flow is linearly proportional to the total flow in
excess of some minimum flow
285
Hydraulic Network Modelling with SWMM
13.3 Conduits
Conduits are pipes or channels that move water from one node (Junctions,
Storage Units, Flow Dividers and Outfalls) to another node in the drainage
system. Their cross-sectional shapes can be selected from a variety of stand-
ard open and closed geometries. Irregular natural cross-section shapes are
also supported.
The Conduits editor organizes the related input data into the following groups:
287
Hydraulic Network Modelling with SWMM
Minor losses are only computed for the DW, Dynamic wave routing option.
The minor losses are computed as K*v2/2g.
13.4 Orifices
Orifices are openings in the wall of a manhole junction or storage unit used to
model outlet and diversion structures in drainage systems. Orifices belong to
the group of Link components connecting two nodes. An orifice can have
either a circular or rectangular shape, be located either at the bottom or along
the side of the upstream node, and have a flap gate to prevent backflow.
289
Hydraulic Network Modelling with SWMM
Orifices attached to storage units are active under all types of flow routing.
Those attached to other types of nodes are active only under Dynamic Wave
flow routing, otherwise the flow through them is zero. The flow through an ori-
fice is computed based on the area of its opening, its discharge coefficient,
and the head difference across the orifice. The area of an orifice's opening
can be controlled dynamically through user-defined control rules.
The Orifices editor organizes the related input data into the following groups:
291
Hydraulic Network Modelling with SWMM
13.5 Pumps
Pumps are drainage system links used to lift water to higher elevations. The
relation between a pump's flow rate and conditions at its end nodes is
described by a pump curve. MIKE URBAN CS SWMM represents a pump as
a link connecting two nodes, where the pump itself is placed at the upstream
node.
The Pumps editor organizes the related input data into the following groups:
Pump Properties
293
Hydraulic Network Modelling with SWMM
For each individual pump a pump curve must be defined. Pump curves are
defined in Tabular Data (Curves). Following four different types of pump
curves are supported:
Type1 - an off-line pump with a wet well where flow increases incremen-
tally with wet well volume
Type2 - an in-line pump where flow increases incrementally with node
depth
Type3 - an in-line pump where flow varies continuously with head differ-
ence between the inlet and outlet nodes
Type4 - a variable speed in-line pump where flow varies continuously
with node depth
The on/off status of pumps can be controlled dynamically through user-
defined Control Rules.
13.6 Weirs
Weirs, like orifices, are used to model outlet and diversion structures in a
drainage system. Weirs are typically located in a chamber of a manhole,
along the side of a channel, or along the top of a storage unit. MIKE URBAN
CS SWMM represents a weir as a link connecting two nodes, where the weir
itself is placed at the upstream node. A flap gate can be included to prevent
backflow.
The Weir editor organizes the related input data into the following groups:
295
Hydraulic Network Modelling with SWMM
Four varieties of weirs (Weir Type) are available, each with a slightly different
formula for computing flow as a function of area, discharge coefficient and
head difference across the weir:
The height of the weir crest above the inlet node invert can be controlled
dynamically through user-defined Control Rules. This feature can be used to
model RTC controlled Weirs and inflatable dams.
13.7 Outlets
Outlets are flow control devices which are typically used to control outflows
from storage units. They are used to model special flow-head relations that
cannot be characterized by either orifices or weirs. MIKE URBAN CS SWMM
represents an Outlet as a link connecting two nodes, where the outlet itself is
placed at the upstream node.
Outlets attached to storage units are active under all types of flow routing.
Those attached to other types of nodes are active only under Dynamic Wave
flow routing, otherwise the flow through them is zero.
The Outlet editor organizes the related input data into the following groups:
297
Hydraulic Network Modelling with SWMM
An outlet can have a flap gate which restricts flow to only one direction and
prevents backwater flow.
13.8 Transects
Transects refer to the geometric data that describe how bottom elevation var-
ies with horizontal distance over the cross section of a natural channel or
irregular-shaped conduit. The figure below displays an example of a transect
for a natural channel.
299
Hydraulic Network Modelling with SWMM
Each transect are given an unique ID. When a IRREGULAR shape is chosen
as the shape type for a conduit then a reference to one of the user defined
Transects must be made for representing the crossectional shape.
The SWMM engine converts these data into tables of area, top width, and
hydraulic radius versus channel depth. In addition, as shown in the figure
above, each Transect can have a left and right overbank section whose Man-
ning's roughness can be different from that of the main channel. This feature
can provide more realistic estimates of channel conveyance under high flow
conditions.
The Transects editor organizes the related input data into the following
groups:
The above Transect can be viewed graphically by clicking on the Graph but-
ton
301
Hydraulic Network Modelling with SWMM
303
Hydraulic Network Modelling with SWMM
The Tabular Data (Curves) editor organizes the related input data into the fol-
lowing groups:
305
Hydraulic Network Modelling with SWMM
13.10 Controls
Controls is used to determine pumps and regulators will be controlled based
on simulation time or conditions at specific nodes and links.
The Controls editor organizes the related input data into the following groups:
307
Hydraulic Network Modelling with SWMM
Please note there can only be either a Condition or an Action in one line.
time-varying rainfall
evaporation of standing surface water
snow accumulation and melting
rainfall interception in depression storage
infiltration of rainfall into upper soil layers
percolation of infiltrated water into groundwater layers
interflow between groundwater and the conveyance system
nonlinear reservoir routing of overland flow
Spatial variability in all of these processes is achieved by dividing a study
area into a smaller collection of homogeneous areas. In the literature different
terms are used for describing the delineation of large study areas into smaller
sub-areas each containing its own fraction of pervious and impervious areas.
Basin
Sub-basin
Watershed/Sewershed
Catchment
Throughout this documentation the term Catchments will be used as term for
the smallest hydrological unit of land.
14.2 Catchments
Catchments are hydrological units of land whose topography cause surface
runoff to drain to a single point. The user is responsible for the delineation of
a watershed/sewershed into a compatible number of catchments, and for
identifying the outlet point of each catchment. Outlet points can be either
nodes of the conveyance system or other catchments.
309
Hydrological Modelling with SWMM
The CS SWMM Catchments Editor organizes the related input data for the
catchment into the following groups. The catchment data in the Catchment
Editor is shared by CS models:
General information
Physical Properties
Location
Location
Physical Properties
311
Hydrological Modelling with SWMM
General Information
Surface runoff per unit area, Q, occurs only when the depth of water in the
"reservoir" exceeds the maximum depression storage, dp, in which case the
outflow is given by Manning's equation (Kinematic Wave).
Depth of water over the catchment is continuously updated with time by solv-
ing numerically a water balance equation over the catchment.
The CS SWMM Routing Editor organizes the related input data for the catch-
ment concerning the surface runoff into following groups:
313
Hydrological Modelling with SWMM
14.4 Infiltration
Infiltration is the process of rainfall penetrating the ground surface into the soil
over the pervious areas of catchments. CS SWMM supports three choices for
modelling infiltration:
Horton's Equation
This method is based on empirical observations showing that infiltration
decreases exponentially from an initial maximum rate to some minimum rate
over the course of a long rainfall event. Input parameters required by this
method include the maximum and minimum infiltration rates, a decay coeffi-
cient that describes how fast the rate decreases over time, and a regenera-
tion constant that describes the restoration of infiltration rate during dry
periods.
Green-Ampt Method
This method for modelling infiltration assumes that a sharp wetting front
exists in the soil column, separating soil with some initial moisture content
below from saturated soil above. The input parameters required are the initial
moisture deficit of the soil, the soil's hydraulic conductivity, and the suction
head at the wetting front.
Only one of the above infiltration methods can be used in one simulation - the
option for which method is used is set in Run SWMM5, tabsheet RUNOFF.
The CS SWMM Infiltration Editor organizes the related input data for these
three infiltration approaches into following groups:
315
Hydrological Modelling with SWMM
Infiltration Properties
14.5 RDII
Rainfall Dependent Inflow/Infiltration (RDII) - These are stormwater flows that
enter sanitary or combined sewers due to "inflow" from direct connections of
downspouts, sump pumps, foundation drains, etc. as well as "infiltration" of
317
Hydrological Modelling with SWMM
subsurface water through cracked pipes, leaky joints, poor manhole connec-
tions, etc.
RDII can be computed for a given rainfall record based on set of triangular
unit hydrographs (UH) that determine a short-term, intermediate-term, and
long-term inflow response for each interval of rainfall. Each unit hydrograph is
defined by three parameters:
The CS SWMM RDII Hydrographs Editor organizes the related input data for
the RDII Unit Hydrograph approach into following groups:
319
Hydrological Modelling with SWMM
Tabular Data
321
Hydrological Modelling with SWMM
14.6 Aquifers
Aquifers are sub-surface groundwater areas used to model the vertical move-
ment of water infiltrating from the catchments which lie above them. They
also permit the infiltration of groundwater into the conveyance system, or
exfiltration of surface water from the conveyance system, depending on the
hydraulic gradient that exists. The same aquifer object can be shared by sev-
eral catchments.
Aquifers are represented using two zones - an unsaturated zone and a satu-
rated zone. Their behaviour is characterized using such parameters as soil
porosity, hydraulic conductivity, evapotranspiration depth, bottom elevation,
and loss rate to deep groundwater. In addition, the initial water table elevation
and initial moisture content of the unsaturated zone must be supplied.
The CS SWMM Aquifer Editor organizes the related input data for the Aquifer
into following groups:
Aquifer Properties
323
Hydrological Modelling with SWMM
14.7 Groundwater
A sketch of the two-zone groundwater model that is used in EPA SWMM is
shown in Figure 14.8. The upper zone is unsaturated at a varying moisture
content of q. The lower zone is saturated and therefore its moisture content is
fixed at the soil porosity h. The fluxes shown in the figure, expressed as vol-
ume per unit area per unit time, consist of the following:
After computing the water fluxes that exist at a given point in time, a mass
balance is written for the change in water volume stored in each zone so that
a new water table depth and unsaturated zone moisture content can be com-
puted for the next time step.
The CS SWMM Groundwater Editor organizes the related input data for the
Groundwater into following groups:
325
Hydrological Modelling with SWMM
Groundwater Properties
14.8 Snowpack
The snow melt routine is part of the runoff modelling process. It updates the
state of the snow packs associated with each catchment by accounting for
snow accumulation, snow redistribution by areal depletion and removal oper-
ations, and snow melt via heat budget accounting. Any snow melt coming off
the pack is treated as an additional rainfall input into the runoff process.
1. Air temperature and melt coefficients are updated according to the calen-
dar date.
3. Any excess snow depth on the plowable area of the pack is redistributed
according to the removal parameters established for the pack.
4. Areal coverages of snow on the impervious and pervious areas of the pack
are reduced according to the Areal Depletion Curves defined for the study
area.
5. The amount of snow in the pack that melts to liquid water is found using:
a heat budget equation for periods with rainfall, where melt rate
increases with increasing air temperature, wind speed, and rainfall inten-
sity
a degree-day equation for periods with no rainfall, where melt rate equals
the product of a melt coefficient and the difference between the air tem-
perature and the pack's base melt temperature.
6. If no melting occurs, the pack temperature is adjusted up or down based
on the product of the difference between current and past air temperatures
and an adjusted melt coefficient. If melting occurs, the temperature of the
pack is increased by the equivalent heat content of the melted snow, up to the
base melt temperature. Any remaining melt liquid beyond this is available to
runoff from the pack.
7. The available snow melt is then reduced by the amount of free water hold-
ing capacity remaining in the pack. The remaining melt is treated the same as
a direct rainfall input onto the catchment.
327
Hydrological Modelling with SWMM
Snowpack Properties
329
Hydrological Modelling with SWMM
14.9 Climatology
Climate-related variables used for computing runoff and snow melt are
grouped in the Climatology object:
Temperature
Evaporation
Wind Speed
Snow Melt
Areal Depletion
Temperature data are used when simulating snowfall and snow melt pro-
cesses during runoff calculations. If these processes are not being simulated
then no temperature data are needed.
Evaporation can occur for standing water on catchment surfaces, for subsur-
face water in groundwater aquifers, and from water held in storage units.
Evaporation rates can be stated as:
If a NCDC file is used then a set of monthly pan coefficients should also be
supplied to convert the pan evaporation data to free-water-surface values.
Wind speed is an optional climatological variable that is only used for snow
melt calculations. There can be use either a set of monthly average speeds
or wind speed data contained in the same NCDC climatological file used for
daily min/max temperatures.
Snow melt parameters are climatological variables that apply across the
entire model area when simulating snowfall and snow melt. They include:
Two such curves can be supplied, one for impervious areas and another for
pervious areas.
The CS SWMM Climatology Editor organizes these data into following tab-
sheets and groups:
Temperature Tabsheet
Evaporation Tabsheet
Wind Speed Tabsheet
331
Hydrological Modelling with SWMM
Temperature data
Evaporation Data
333
Hydrological Modelling with SWMM
Windspeed data
335
Hydrological Modelling with SWMM
Snowmelt data
337
Hydrological Modelling with SWMM
339
Hydrological Modelling with SWMM
14.10 Coverage
The CS SWMM Coverage Editor organizes the related input data for the cov-
erage into following groups:
Landuse coverage
341
Hydrological Modelling with SWMM
Bio-retention Cells
Infiltration Trenches
Infiltration Trenches are narrow ditches filled with gravel that intercept runoff
from upslope impervious areas. They provide storage volume and additional
time for captured runoff to infiltrate the native soil below.
Continuous Porous Pavement systems are excavated areas filled with gravel
and paved over with a porous concrete or asphalt mix. Normally all rainfall will
immediately pass through the pavement into the gravel storage layer below it
where it can infiltrate at natural rates into the site's native soil. Block Paver
systems consist of impervious paver blocks placed on a sand or pea gravel
bed with a gravel storage layer below. Rainfall is captured in the open spaces
between the blocks and conveyed to the storage zone and native soil below.
Rain Barrels
Rain Barrels (or Cisterns) are containers that collect roof runoff during storm
events and can either release or re-use the rainwater during dry periods.
Vegetative Swales
Vegetative Swales are channels or depressed areas with sloping sides cov-
ered with grass and other vegetation. They slow down the conveyance of col-
lected runoff and allow it more time to infiltrate the native soil beneath it.
Bio-retention cells, infiltration trenches, and porous pavement systems can all
contain optional underdrain systems in their gravel storage beds to convey
captured runoff off of the site rather than letting it all infiltrate. They can also
have an impermeable floor or liner that prevents any infiltration into the native
soil from occurring. Infiltration trenches and porous pavement systems can
also be subjected to a decrease in hydraulic conductivity over time due to
clogging.
Although some LID practices can also provide significant pollutant reduction
benefits, at this time SWMM only models their hydrologic performance
The CS SWMM LID Controls Editor organizes the data input for the different
types of LID practices. The input data is organized into following group and
tabs:
Identification
Tabs: Surface, Soil, Pavement, Storage and Drain
343
Hydrological Modelling with SWMM
Identification
345
Hydrological Modelling with SWMM
347
Hydrological Modelling with SWMM
349
Hydrological Modelling with SWMM
351
Hydrological Modelling with SWMM
The boundary conditions can be visualized on the map. Figure 15.1 shows
how where you can toggle the visualization on/off in the TOC. Boundary con-
ditions are per default unchecked (one exception: catchment connections)
and not displayed in the horizontal view but they can be activated at any time
during the project.
The boundaries that can be seen on the map are explained in below.
353
Time Series and Boundary Conditions in SWMM
SWMM5 Outfall shows outfalls where a water level has been assigned.
The boundary condition will only be displayed if the Outfall type is set to
one of the following type: Fixed (3), Tidal (4) or Timeseries (5).
SWMM5 Catchment Connection shows the link between the center of
the catchment and the node it is connected to
SWMM5 Raingauge to Catchment displays the link between the rain-
gauge and the center of the catchment when a raingauge is specified in
the Routing dialog
SWMM5 Catchment to Catchment displays the link between the
centers of two connected catchments
SWMM5 Rainfall shows the catchments where a raingauge has been
assigned (this one is similar to the SWMM5 Raingauge to Catchment but
will show all catchments where a raingauge has been assigned to)
An example on how this may look is shown in Figure 15.2.
The SWMM engine supports two time series formats. Time series values can
either be in date / time / value format or in time / value format, where each
entry is separated by one or more spaces or tab characters.
For the date / time / value format which is the only format supported in MIKE
URBAN interface, dates are entered as month/day/year (e.g., 7/21/2004) and
times as 24-hour military time (e.g., 8:30 pm is 20:30).
For the time / value format, times are entered as hours since simulation start.
If the time / value format has been used outside MIKE URBAN interface, then
it is converted during import to the date / time /value format given the start of
simulation date as initial start date.
Temperature data
Evaporation data
Rainfall data
Water stage at outfall nodes
External inflow hydrographs at drainage system nodes
External inflow pollutographs at drainage system nodes.
The CS SWMM Time Series editor organizes the related input data for the dif-
ferent timeseries into the following groups:
355
Time Series and Boundary Conditions in SWMM
Monthly
Daily
Hourly
Weekend
The CS SWMM Time Patterns editor organizes the related input data for the
different time patterns into the following groups:
357
Time Series and Boundary Conditions in SWMM
359
Time Series and Boundary Conditions in SWMM
15.3 Raingauges
Rain Gauges supply precipitation data for one or more catchments in a study
area. The rainfall data can be either a user-defined time series (specified in
Time Series as Type Rainfall) or come from an external file.
Several different popular rainfall file formats currently in use are supported,
as well as a standard user-defined format. These are:
When a rain gauge is designated as receiving its rainfall data from a file, the
user must supply the name of the file and the name of the recording station
referenced in the file.
For the standard user-prepared format, the rainfall type (e.g., intensity or vol-
ume), recording time interval, and depth units must also be supplied as rain
gauge properties. For the other file types these properties are defined by their
respective file format and are automatically recognized by SWMM.
The CS SWMM Raingauges editor organizes the related input data for the
rainfall data into the following groups:
361
Time Series and Boundary Conditions in SWMM
stream channels. They are represented by an average inflow rate that can be
periodically adjusted on a monthly, daily, and hourly basis by applying Time
Patterns multipliers to this average value.
The CS SWMM Dry Weather Flow editor organizes the related input data for
each dry weather inflow into the following groups:
Table 15.6 The Dry Weather Flow Identification and Connectivity Group
363
Time Series and Boundary Conditions in SWMM
Table 15.8 The Pollutant Data in upper grid attached to dry weather flow
Edit field Description Used or required Field name in
by simulations datastructure
PollutantID ID of Pollutant Yes, if pollutant PollutantID
attached
BValue Base value of pollutant Yes, if pollutant BValue
attached
Monthly ID of Pattern of Monthly type Yes, if pollutant PatternMonthID
attached
Daily ID of Pattern of Daily type Yes, if pollutant PatternWeekID
attached
Week hourly ID of Pattern of Hourly type Yes, if pollutant PatternWeek-
attached HourlyID
Weekend ID of Pattern of Weekend type Yes, if pollutant PatternWeek-
hourly attached endHourlyID
15.5 Inflow
Direct Inflows are user-defined time series of inflows added directly into a
node. They can be used to perform flow and water quality routing in the
absence of any runoff computations (as in a study area where no catchments
are defined). A pollutant component can only be defined as inflow to the sys-
tem if it is attached to a flow component.
The CS SWMM Inflow editor organizes the related input data for each inflow
into the following groups:
365
Time Series and Boundary Conditions in SWMM
The CS-SWMM Run SWMM5 editor is split into following four tabsheets:
General
Runoff
Network
Summary
367
Project Options and Simulations in SWMM
Simulation Period
Start Simulation
Units selection
369
Project Options and Simulations in SWMM
Infiltration Method
371
Project Options and Simulations in SWMM
Routing Simulation
373
Project Options and Simulations in SWMM
375
Project Options and Simulations in SWMM
Summary Specifications
Subcatchment Summary
Node Summary
Link Summary
377
Project Options and Simulations in SWMM
Water quality routing within conduit links assumes that the conduit behaves
as a continuously stirred tank reactor (CSTR). Although a plug flow reactor
assumption might be more realistic, the differences will be small if the travel
time through the conduit is on the same order as the routing time step. The
concentration of a constituent exiting the conduit at the end of a time step is
found by integrating the conservation of mass equation, using average values
for quantities that might change over the time step such as flow rate and con-
duit volume.
Water quality modelling within storage unit nodes follows the same approach
used for conduits. For other types of nodes that have no volume, the quality
of water exiting the node is simply the mixture concentration of all water
entering the node.
The production and fate of pollutant loads associated with this runoff is
also modelled. The following processes can be modelled for any number of
user-defined water quality constituents:
379
Water Quality Modelling with SWMM
include roof tops, lawns, paved roads, undisturbed soils, etc. Land uses are
used solely to allow spatial variation in pollutant buildup and washoff rates.
The following processes can be associated with each land use category:
Pollutant Buildup
Pollutant Washoff
Street Cleaning
The CS SWMM Land Use editor organizes the related input data for the land
uses into the following groups:
17.3 Buildup
Pollutant Buildup that accumulates over a category of Land Use is described
by either a mass per unit of catchment area or per unit of curb length. The
amount of buildup as a function of days of dry weather can be computed
using one of the following functions:
Power Function
C3
B = Min C 1 C 2 t (17.1)
where C1 = maximum buildup possible (mass per unit of area or curb length),
C2 = buildup rate constant, and C3 = time exponent.
Exponential Function
381
Water Quality Modelling with SWMM
– C 2t
B = C1 1 – e (17.2)
where C1 = maximum buildup possible (mass per unit of area or curb length)
and C2 = buildup rate constant (1/days).
Saturation Function
Buildup begins at a linear rate which proceeds to decline constantly over time
until a saturation value is reached,
C1 t
B = -------------
- (17.3)
C2 + t
where C1 = maximum buildup possible (mass per unit area or curb length)
and C2 = half-saturation constant (days to reach half of the maximum
buildup).
The CS SWMM Buildup Editor organizes the related input data for the pollut-
ant buildup into the following groups:
Buildup Properties
17.4 Washoff
Pollutant Washoff from a given land use category during wet weather periods
can be described in one of the following ways:
Exponential Washoff
383
Water Quality Modelling with SWMM
The washoff load (W) in units of mass per hour is proportional to the product
of runoff raised to some power and to the amount of buildup remaining, i.e.,
C2
W = C1 q B (17.4)
The rate of washoff W in mass per second is proportional to the runoff rate
raised to some power, i.e.,
C2
W = C1 Q (17.5)
This is a special case of Rating Curve Washoff where the exponent is 1.0 and
the coefficient C1 represents the concentration of any and all runoff in mass
per liter (the conversion between user-specified flow units used for runoff and
liters is handled internally by SWMM).
Washoff loads for a given pollutant and land use category can be reduced by
a fixed percentage by specifying a BMP Removal Efficiency which reflects the
effectiveness of any BMP controls associated with the land use.
The CS SWMM Buildup Editor organizes the related input data for the pollut-
ant washoff into the following groups:
Washoff Properties
385
Water Quality Modelling with SWMM
17.5 Loading
Input loadings of pollutants from external and dry weather inflows are sup-
plied through time series data associated with particular nodes of the collec-
tion system.The definition of the boundary connection of external timeseries
to a node is done in the CS SWMM Loading Editor.
The CS SWMM Loading Editor organizes the related data for setting up the
boundary connection of the pollutant inflow timeseries to the model into the
following groups:
Loading Properties
17.6 Pollutant
The generation, inflow and fate of any number of user-specified pollutants
can be modelled. Concentrations in rainfall, in groundwater, in inflow/infiltra-
tion is modelled.
There can be used following different units for the individual user defined pol-
lutant units (milligrams/liter, micrograms/liter, or counts/liter)
Co-pollutants can also be modelled. For example, pollutant X can have a co-
pollutant Y, meaning that the runoff concentration of X will have some fixed
fraction of the runoff concentration of Y added to it.
The pollutant buildup and washoff on catchment areas are determined by the
Land Uses assigned to those areas.
The CS SWMM Pollutant Editor organizes the related input data for the pol-
lutants into the following groups:
387
Water Quality Modelling with SWMM
Pollutant Properties
For example, a first-order decay expression for BOD exiting from a storage
node might be:
– 0,05 HRT
C = BOD e (17.6)
The CS SWMM Treatment Editor organizes the related input data for the local
treatment into the following groups:
389
Water Quality Modelling with SWMM
Treatment Properties
457
Index
B R
Basin geometry . . . . . . . . . . . . . 157 RDI . . . . . . . . . . . . . . . . . . . .92
Biological Processes . . . . . . . . . . 296
S
C Sediment transport . . . . . . . . . . . 295
Capacity Curve . . . . . . . . . . . . . 157 Sensor . . . . . . . . . . . . . . . . . 242
Conduits . . . . . . . . . . . . . . . . 352 Soakaway. . . . . . . . . . . . . . . . .29
Control action . . . . . . . . . . . . . . 245 Storage Units . . . . . . . . . . . . . . 344
Controllable devices . . . . . . . . . . 251
Cross sections . . . . . . . . . . . . . 66 T
Time series . . . . . . . . . . . . . . . 137
F Time-Area Curve . . . . . . . . . . . . 157
Flooding . . . . . . . . . . . . . . . . 207 Time-Area Method . . . . . . . . . . . .75
Flow Dividers . . . . . . . . . . . . . . 344
FST . . . . . . . . . . . . . . . . . . . 208 U
Unit Hydrograph Method . . . . . . . . .89
H
Hydrological models . . . . . . . . . . 73 W
Weirs . . . . . . . . . . . . . . . . . . .41
I
Inlet . . . . . . . . . . . . . . . . . . . 49
J
Job list . . . . . . . . . . . . . . . . . 266
Junctions . . . . . . . . . . . . . . . . 344
K
Kinematic Wave . . . . . . . . . . . . 80
L
Linear Reservoir . . . . . . . . . . . . 84
Long Term Statistics . . . . . . . . . . 265
N
Nodes and Structures . . . . . . . . . 22
O
Orifice . . . . . . . . . . . . . . . . . . 44
Outfalls . . . . . . . . . . . . . . . . . 344
P
PID . . . . . . . . . . . . . . . . . . . 246
Pipe Design . . . . . . . . . . . . . . . 283
Pipes and Canals . . . . . . . . . . . . 36
Pump . . . . . . . . . . . . . . . . . . 54
Q
QH Relation . . . . . . . . . . . . . . . 157