Collection System
Collection System
Collection System
COLLECTION SYSTEM
Limited Liability
The liability of DHI is limited as specified in Section III of your 'DHI
Software Licence Agreement':
'IN NO EVENT SHALL DHI OR ITS REPRESENTA-TIVES (AGENTS
AND SUPPLIERS) BE LIABLE FOR ANY DAMAGES WHATSO-
EVER INCLUDING, WITHOUT LIMITATION, SPECIAL, INDIRECT,
INCIDENTAL OR CONSEQUENTIAL DAMAGES OR DAMAGES
FOR LOSS OF BUSINESS PROFITS OR SAVINGS, BUSINESS
INTERRUPTION, LOSS OF BUSINESS INFORMATION OR OTHER
PECUNIARY LOSS ARISING OUT OF THE USE OF OR THE INA-
BILITY TO USE THIS DHI SOFTWARE PRODUCT, EVEN IF DHI
HAS BEEN ADVISED OF THE POSSIBILITY OF SUCH DAMAGES.
THIS LIMITATION SHALL APPLY TO CLAIMS OF PERSONAL
INJURY TO THE EXTENT PERMITTED BY LAW. SOME COUN-
TRIES OR STATES DO NOT ALLOW THE EXCLUSION OR LIMITA-
TION OF LIABILITY FOR CONSEQUENTIAL, SPECIAL, INDIRECT,
INCIDENTAL DAMAGES AND, ACCORDINGLY, SOME PORTIONS
OF THESE LIMITATIONS MAY NOT APPLY TO YOU. BY YOUR
OPENING OF THIS SEALED PACKAGE OR INSTALLING OR
USING THE SOFTWARE, YOU HAVE ACCEPTED THAT THE
ABOVE LIMITATIONS OR THE MAXIMUM LEGALLY APPLICA-
BLE SUBSET OF THESE LIMITATIONS APPLY TO YOUR PUR-
CHASE OF THIS SOFTWARE.'
3
4 MIKE URBAN
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 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
Identification group . . . . . . . . . . . . . . . . . . . . . . . . . . 25
MOUSE model data group . . . . . . . . . . . . . . . . . . . . . . 27
Q-H relations for nodes . . . . . . . . . . . . . . . . . . . . . . . . 29
3.4 Pipes and Canals . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
Identification group . . . . . . . . . . . . . . . . . . . . . . . . . . 32
Geometrical properties . . . . . . . . . . . . . . . . . . . . . . . . 33
Hydraulic friction losses . . . . . . . . . . . . . . . . . . . . . . . . 34
Miscellaneous . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
3.5 Weirs . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36
Identification and connectivity . . . . . . . . . . . . . . . . . . . . 38
Model data . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39
3.6 Orifices . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40
Identification and connectivity . . . . . . . . . . . . . . . . . . . . 42
Model data . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 43
3.6.1 Defining a gate or a weir in an orifice . . . . . . . . . . . . . . . . 44
3.7 Stormwater Inlets . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 44
3.7.1 Curb Inlet (Lintel) . . . . . . . . . . . . . . . . . . . . . . . . . . . 46
3.7.2 On-grade Capture . . . . . . . . . . . . . . . . . . . . . . . . . . . 48
3.7.3 Capacity curves . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
3.8 Pumps . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
3.8.1 Pump types . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50
Constant flow pumps . . . . . . . . . . . . . . . . . . . . . . . . . 50
Constant speed pumps . . . . . . . . . . . . . . . . . . . . . . . . 50
Variable speed pumps . . . . . . . . . . . . . . . . . . . . . . . . 52
Identification and connectivity . . . . . . . . . . . . . . . . . . . . 57
Model data . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 58
3.9 Valves . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 59
3.10 CRS & Topography . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 62
3.11 Emptying Storage Nodes . . . . . . . . . . . . . . . . . . . . . . . . . . . . 65
4 RAINFALL-RUNOFF MODELLING WITH MOUSE . . . . . . . . . . . . . . . . . 67
4.1 Terms and Concepts . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 67
4.1.1 MIKE URBAN Catchments . . . . . . . . . . . . . . . . . . . . . . 68
5
4.1.2 Connecting Catchments to the Network . . . . . . . . . . . . . . 69
4.1.3 Hydrological Models . . . . . . . . . . . . . . . . . . . . . . . . . 69
4.1.4 Creating Hydrological Models for a Catchment . . . . . . . . . . 70
4.2 Time-Area Method (A) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 71
4.2.1 Model Data . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 71
4.2.2 Parameter Sets . . . . . . . . . . . . . . . . . . . . . . . . . . . . 74
4.2.3 Time-Area Curve Editor . . . . . . . . . . . . . . . . . . . . . . . 75
4.3 Kinematic Wave (B) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 77
4.3.1 Model Data . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 77
4.3.2 Parameter Sets . . . . . . . . . . . . . . . . . . . . . . . . . . . . 78
4.4 Linear Reservoir (C1 and C2) . . . . . . . . . . . . . . . . . . . . . . . . . 80
4.4.1 Model Data . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 81
4.4.2 Parameter Sets . . . . . . . . . . . . . . . . . . . . . . . . . . . . 84
4.5 Unit Hydrograph Method (UHM) . . . . . . . . . . . . . . . . . . . . . . . 86
4.6 Additional Flow and RDI . . . . . . . . . . . . . . . . . . . . . . . . . . . . 90
4.6.1 Model Data . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 93
4.6.2 MOUSE RDI Parameter Sets . . . . . . . . . . . . . . . . . . . . 95
4.7 Rainfall Data and Other Meteorological Variables - Boundary Conditions for
Hydrological models 98
4.8 Running the Runoff Computations . . . . . . . . . . . . . . . . . . . . . 101
4.9 MOUSE RDI - Guidelines for application . . . . . . . . . . . . . . . . . 103
4.9.1 Choice of calculation time step . . . . . . . . . . . . . . . . . . 103
4.9.2 The RDI hotstart . . . . . . . . . . . . . . . . . . . . . . . . . . . 104
4.9.3 The RDI result files . . . . . . . . . . . . . . . . . . . . . . . . . 104
4.9.4 MOUSE RDI Validation . . . . . . . . . . . . . . . . . . . . . . . 105
Surface runoff model . . . . . . . . . . . . . . . . . . . . . . . . . 106
General hydrological model - RDI . . . . . . . . . . . . . . . . . . 107
4.9.5 Overflow within the model area . . . . . . . . . . . . . . . . . . 111
4.9.6 Non-precipitation dependent flow components . . . . . . . . . 111
4.10 Using the Computed Runoff as Network Hydraulic Load . . . . . . . . 112
5 TIME SERIES . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 115
5.1 Inserting new time series . . . . . . . . . . . . . . . . . . . . . . . . . . . 115
5.1.1 Properties of time series object . . . . . . . . . . . . . . . . . . 116
5.1.2 Properties of time series item . . . . . . . . . . . . . . . . . . . 119
5.1.3 Time series plot properties . . . . . . . . . . . . . . . . . . . . . 123
5.2 Example: How to enter a rain time series . . . . . . . . . . . . . . . . . 127
5.3 Example: How to import a time series from Excel . . . . . . . . . . . . 132
6 CURVES AND RELATIONS . . . . . . . . . . . . . . . . . . . . . . . . . . . . 139
6.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 139
6.1.1 Capacity curves . . . . . . . . . . . . . . . . . . . . . . . . . . . 140
6.1.2 Pump acceleration curves . . . . . . . . . . . . . . . . . . . . . 141
6 MIKE URBAN
6.1.3 Regulation curves Qmax(H) and Qmax(dH) . . . . . . . . . . . 141
6.1.4 QH relation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 141
Manholes, Basins . . . . . . . . . . . . . . . . . . . . . . . . . . . 141
Outlets . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 141
Storage node . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 141
6.1.5 Time-Area curve . . . . . . . . . . . . . . . . . . . . . . . . . . . 141
6.1.6 Removal efficiency . . . . . . . . . . . . . . . . . . . . . . . . . 142
6.1.7 Basin geometry . . . . . . . . . . . . . . . . . . . . . . . . . . . 142
6.1.8 Valve rating curve . . . . . . . . . . . . . . . . . . . . . . . . . . 143
6.1.9 DQ and QQ relations . . . . . . . . . . . . . . . . . . . . . . . . 143
6.1.10 Capacity curve QdH & Power . . . . . . . . . . . . . . . . . . . 143
6.1.11 Undefined type . . . . . . . . . . . . . . . . . . . . . . . . . . . . 144
7 BOUNDARY CONDITIONS . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 147
7.1 Catchment Boundary Condition . . . . . . . . . . . . . . . . . . . . . . . 148
7.1.1 Visualization of catchment boundaries . . . . . . . . . . . . . . 149
7.2 Network Boundary Conditions . . . . . . . . . . . . . . . . . . . . . . . . 151
7.2.1 Visualization of network boundaries . . . . . . . . . . . . . . . . 152
7.3 Boundary Condition Editors . . . . . . . . . . . . . . . . . . . . . . . . . 154
7.3.1 Catchments Loads and Meteorological Items Editor . . . . . . 154
7.3.2 Network Loads Editor . . . . . . . . . . . . . . . . . . . . . . . . 156
7.3.3 External Water Levels Editor . . . . . . . . . . . . . . . . . . . . 157
7.3.4 Boundary Items Editor . . . . . . . . . . . . . . . . . . . . . . . 158
7.4 Examples . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 160
7.4.1 How to add a varying water level at an outlet ? . . . . . . . . . 160
7.4.2 How to add infiltration in a pipe? . . . . . . . . . . . . . . . . . . 161
7.4.3 How to add a rainfall as a boundary condition to the catchments? .
162
7.4.4 How to add a discharge to a node? . . . . . . . . . . . . . . . . 162
7.4.5 How to add runoff results as input for the network computation? . .
163
7.4.6 How do I add DWF in my network dependent on number of inhabit-
ants? 164
7.4.7 How to attach a pollutant concentration to a network load? . . 166
7.5 Repetitive Profile Editors . . . . . . . . . . . . . . . . . . . . . . . . . . . 167
7.6 Diurnal Patterns . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 168
7.7 Profiles Calendar . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 168
7.8 Cyclic Profiles . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 169
7.9 Special Days . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 170
8 MOUSE SIMULATIONS . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 173
8.1 The General Simulation Settings . . . . . . . . . . . . . . . . . . . . . . 173
8.1.1 Batch simulations . . . . . . . . . . . . . . . . . . . . . . . . . . 176
7
8.2 The Runoff Simulation Settings . . . . . . . . . . . . . . . . . . . . . . . 177
8.3 The Network Simulation Settings . . . . . . . . . . . . . . . . . . . . . . 180
8.4 The Summary Simulation Settings . . . . . . . . . . . . . . . . . . . . . 183
8.5 The 2D Overland Parameters . . . . . . . . . . . . . . . . . . . . . . . . 186
8.6 MOUSE Result Selections . . . . . . . . . . . . . . . . . . . . . . . . . . 186
9 2D OVERLAND FLOW . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 191
9.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 191
9.2 Input required - overview . . . . . . . . . . . . . . . . . . . . . . . . . . . 191
9.3 Input required - details . . . . . . . . . . . . . . . . . . . . . . . . . . . . 192
9.3.1 Defining 2D Domain, Resolution and Parameters . . . . . . . 192
9.3.2 Defining Couplings . . . . . . . . . . . . . . . . . . . . . . . . . 196
9.3.3 Flow parameters at nodes and basins . . . . . . . . . . . . . . 200
9.3.4 Outlets . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 202
9.3.5 Pumps and weirs . . . . . . . . . . . . . . . . . . . . . . . . . . 202
9.4 Running the combined 1D and 2D simulations . . . . . . . . . . . . . . 202
9.4.1 Setting the Simulation Type and Requesting 2D results . . . . 202
9.4.2 Starting Simulation . . . . . . . . . . . . . . . . . . . . . . . . . 206
9.5 Visualising simulation results . . . . . . . . . . . . . . . . . . . . . . . . 207
9.6 Understanding the files generated . . . . . . . . . . . . . . . . . . . . . 209
10 MOUSE CONTROL MODULE . . . . . . . . . . . . . . . . . . . . . . . . . . . 211
10.1 RTC in Urban Drainage and Sewer Systems . . . . . . . . . . . . . . . 211
10.2 Architecture of RTC Systems . . . . . . . . . . . . . . . . . . . . . . . . 212
10.3 MOUSE Control vs. Real Life . . . . . . . . . . . . . . . . . . . . . . . . 214
10.4 Sensors . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 215
10.5 Logical Conditions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 216
10.6 Control Actions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 217
10.7 PID parameter sets . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 219
10.7.1 Calibration of the PID constants . . . . . . . . . . . . . . . . . . 220
10.8 Controllable devices . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 222
Control Type and PID-ID . . . . . . . . . . . . . . . . . . . . . . . 224
10.8.1 Pumps . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 225
10.8.2 Weirs . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 226
10.8.3 Orifices with weirs and gates . . . . . . . . . . . . . . . . . . . 226
10.8.4 Difference between weir and orifice with weir . . . . . . . . . . 227
10.8.5 Valves . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 229
10.8.6 Control rules . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 229
10.9 MOUSE Control Computations . . . . . . . . . . . . . . . . . . . . . . . 230
10.10 User Written Control . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 231
11 LONG TERM STATISTICS . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 239
11.1 Data input . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 240
8 MIKE URBAN
11.1.1 Job list . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 240
11.1.2 Job list criteria . . . . . . . . . . . . . . . . . . . . . . . . . . . . 241
11.1.3 Initial conditions for simulated events . . . . . . . . . . . . . . . 243
11.1.4 Generating job list . . . . . . . . . . . . . . . . . . . . . . . . . . 245
11.1.5 Edit job list . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 246
11.1.6 Runtime stop criteria . . . . . . . . . . . . . . . . . . . . . . . . 247
Run-Time Stop Criteria Evaluation Matrix . . . . . . . . . . . . . 248
11.2 LTS computations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 249
11.3 Result files . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 251
11.3.1 User-Specified result files . . . . . . . . . . . . . . . . . . . . . 251
11.3.2 Statistics result file . . . . . . . . . . . . . . . . . . . . . . . . . . 252
11.4 Specification of statistical result file . . . . . . . . . . . . . . . . . . . . . 252
11.5 LTS Statistics Presentation . . . . . . . . . . . . . . . . . . . . . . . . . 258
12 AUTOMATIC PIPE DESIGN WITH MOUSE . . . . . . . . . . . . . . . . . . . 259
12.1 Design principles . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 259
12.2 Design Input . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 260
12.2.1 Example of an ADP-file . . . . . . . . . . . . . . . . . . . . . . . 262
12.2.2 Design Type . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 262
12.2.3 Design Criteria . . . . . . . . . . . . . . . . . . . . . . . . . . . . 263
12.2.4 Design Group Type . . . . . . . . . . . . . . . . . . . . . . . . . 263
12.2.5 Lower Limit . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 264
12.2.6 Commercial Diameters . . . . . . . . . . . . . . . . . . . . . . . 264
12.2.7 Creating the ADP file for the design simulation . . . . . . . . . 264
12.3 Design simulation and output . . . . . . . . . . . . . . . . . . . . . . . . 265
13 MODELLING WATER QUALITY WITH MOUSE . . . . . . . . . . . . . . . . . 271
13.1 Key features and application domain . . . . . . . . . . . . . . . . . . . . 271
13.1.1 Surface Runoff Quality (SRQ) . . . . . . . . . . . . . . . . . . . 271
13.1.2 Pipe Sediment Transport (ST) . . . . . . . . . . . . . . . . . . . 271
13.1.3 Pipe Advection-Dispersion (AD) . . . . . . . . . . . . . . . . . . 272
13.1.4 Biological Processes (BP) . . . . . . . . . . . . . . . . . . . . . 272
13.1.5 Interaction between water quality modules . . . . . . . . . . . . 273
13.2 Surface Runoff Quality (SRQ) . . . . . . . . . . . . . . . . . . . . . . . . 274
13.2.1 Surface Sediment Data Dialogs . . . . . . . . . . . . . . . . . . 275
MOUSE| Water Quality| SRQ Main Parameters . . . . . . . . . . 275
MOUSE| Water Quality| SRQ Attached Pollutants . . . . . . . . 276
MOUSE| Water Quality| SRQ Gully Pot Data . . . . . . . . . . . 277
13.3 Advection-Dispersion (AD) . . . . . . . . . . . . . . . . . . . . . . . . . . 278
13.3.1 Advection-Dispersion Data Dialogs . . . . . . . . . . . . . . . . 279
MOUSE| Water Quality| AD Components . . . . . . . . . . . . . 279
MOUSE | Water Quality | AD Dispersion . . . . . . . . . . . . . . 281
Advection-Dispersion and Open Boundary Conditions . . . . . . 282
9
13.4 Biological Processes (BP) . . . . . . . . . . . . . . . . . . . . . . . . . . 283
13.4.1 Biologoical Processes Dialog (MOUSE|Water Quality|WQ Process
Model) 284
13.5 Sediment Transport (ST) . . . . . . . . . . . . . . . . . . . . . . . . . . . 287
13.5.1 The Sediment Transport Models in MOUSE ST . . . . . . . . 288
The Explicit Sediment Transport Models . . . . . . . . . . . . . . 288
The Morphological Models . . . . . . . . . . . . . . . . . . . . . . 288
13.5.2 The Transport Formulae - Short Description . . . . . . . . . . 289
The Ackers-White formulae . . . . . . . . . . . . . . . . . . . . . 289
The Engelund-Hansen formula . . . . . . . . . . . . . . . . . . . 289
The Engelund-Fredsøe-Deigaard formulae . . . . . . . . . . . . 290
The van Rijn formulae . . . . . . . . . . . . . . . . . . . . . . . . 290
13.5.3 The Flow Resistance in Sewer Systems with Sediment Deposits .
290
13.5.4 Sediment Transport Data Dialogs . . . . . . . . . . . . . . . . . 290
MOUSE | Water Quality | ST Main Parameters . . . . . . . . . . 291
MOUSE| Water Quality | ST Sediment Fractions . . . . . . . . . 293
MOUSE | Water Quality | ST Initial Sediment Depth Local . . . 295
MOUSE| Water Quality| ST Sediment Removal Basins . . . . . 295
MOUSE | Water Quality| ST Sediment Removal Weirs . . . . . 296
13.5.5 Boundary Conditions for the Sediment Transport Model . . . . 297
13.6 Storm Water Quality . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 298
13.6.1 Assigning concentrations of pollutants to runoff and infiltrations . .
298
Cst. concentration (method 1) . . . . . . . . . . . . . . . . . . . . 301
Table concentration (method 2) . . . . . . . . . . . . . . . . . . . 301
EMC formula (method 3) . . . . . . . . . . . . . . . . . . . . . . . 301
13.6.2 Local treatment nodes . . . . . . . . . . . . . . . . . . . . . . . 302
10 MIKE URBAN
15 HYDROLOGICAL MODELLING WITH SWMM . . . . . . . . . . . . . . . . . . 345
15.1 Terms and Concepts . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 345
15.2 Catchments . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 346
15.3 Surface Routing . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 348
15.4 Infiltration . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 351
15.5 RDII . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 354
15.6 Aquifers . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 358
15.7 Groundwater . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 360
15.8 Snowpack . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 363
15.9 Climatology . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 367
15.10 Coverage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 378
15.11 LID Controls . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 379
15.12 LID Deployment . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 389
16 TIME SERIES AND BOUNDARY CONDITIONS IN SWMM . . . . . . . . . . 393
16.1 Time Series . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 395
16.2 Time Patterns . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 397
16.3 Raingauges . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 401
16.4 Dry Weather Flow . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 403
16.5 Inflow . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 406
17 PROJECT OPTIONS AND SIMULATIONS IN SWMM . . . . . . . . . . . . . . 409
17.1 The General Simulation Settings . . . . . . . . . . . . . . . . . . . . . . 409
17.2 The Runoff Simulation Settings . . . . . . . . . . . . . . . . . . . . . . . 411
17.3 The Network Simulation Settings . . . . . . . . . . . . . . . . . . . . . . 414
17.4 The Summary Simulation Settings . . . . . . . . . . . . . . . . . . . . . 418
18 WATER QUALITY MODELLING WITH SWMM . . . . . . . . . . . . . . . . . . 421
18.1 Terms and Concepts . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 421
18.2 Land Uses . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 422
18.3 Buildup . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 423
18.4 Washoff . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 425
18.5 Loading . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 428
18.6 Pollutant . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 429
18.7 Local Treatments . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 431
11
12 MIKE URBAN
MIKE URBAN CS - MOUSE
User Guide
13
14 MIKE URBAN
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.
In order to run SWMM5 simulation a Model Manager module is required,
while running MOUSE simulations require some further modules depend-
ing on type of simulation being carried out (e.g. pipeflow, rainfall runoff
simulations).
16 MIKE URBAN
2 MODELLING COLLECTION SYSTEMS WITH
MOUSE
MOUSE is a powerful and comprehensive engine for modelling complex
hydrology, advanced hydraulics in both open and closed conduits, water
quality and sediment transport for urban drainage systems, storm water
sewers and sanitary sewers.
z What are the return periods for overloading of various parts of the
existing sewer system?
z What are the main causes of that overloading - backwater or insuffi-
cient local pipe capacity?
z What are the implications of replacing critical sewers, installing new
basins, weirs, etc.?
z How is the long-term environmental impact affected by changing the
operational strategy?
z Where and why are sediments deposited in the sewer network?
z 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
18 MIKE URBAN
Introduction
3.1 Introduction
MOUSE allows for the hydrodynamic simulation of flows and water lev-
els in urban storm drainage and wastewater collection networks, thus pro-
viding an accurate information about the network functionality under a
variety of boundary 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 func-
tions. 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
20 MIKE URBAN
Definition of a MOUSE Network
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 vol-
ume. Still only a single water level computational point is located at the
node. Based on the computed water level and the description of the geom-
etry of the node the computation 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
previous time step and the flow contributions during the time step from
21
Hydraulic Network Modeling with MOUSE
each connected 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.
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 connec-
tion between two nodes.
22 MIKE URBAN
Nodes and Structures
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 dia-
log where you find fields for entering "From node:" and "To node:".
The editor organizes the related input data into the following groups:
23
Hydraulic Network Modeling with MOUSE
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.
24 MIKE URBAN
Nodes and Structures
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 approxi-
mated as a free overflow over the crest at a given level and with a "con-
ceptual" 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.
Identification group
25
Hydraulic Network Modeling with MOUSE
26 MIKE URBAN
Nodes and Structures
Table 3.2 The edit fields in the MOUSE model data group
Cover sub-group
27
Hydraulic Network Modeling with MOUSE
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
28 MIKE URBAN
Nodes and Structures
Table 3.2 The edit fields in the MOUSE model data group
29
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 click-
ing the “Advanced” button.
30 MIKE URBAN
Pipes and Canals
Depending on the selected type, a link may take the form of one of the
'standard' pipes (Circular, Rectangular, O Shaped, Egg-Shaped), or any
closed or open cross section shape (CRS) and Natural Channels. The CRS
and Natural 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 men-
tioned specified under the cross section editor. In this dialog, only the ref-
erence to the CRS ID.
31
Hydraulic Network Modeling with MOUSE
The length of a link is calculated from the shape of the line in MIKE
URBAN. The length is displayed in the 'length' field, upon executing the
'Recompute' command when clicking the ‘Advanced’ button. If a user
defined length is specified this will overwrite the calculated one during
simulation.
Identification group
32 MIKE URBAN
Pipes and Canals
Geometrical properties
Table 3.4
33
Hydraulic Network Modeling with MOUSE
Table 3.4
Table 3.5
34 MIKE URBAN
Pipes and Canals
Table 3.5
Miscellaneous
The ‘Regulation’ button provides access to inserting a regulation in the
selected link. This regulation does not require the Control module. The
regulation can be either a maximum discharge as a function of the water
level in a user specified node (Ctrl. Node A) or a maximum discharge as a
function of the water level difference between two user specified nodes
(Ctrl. Node A and Ctrl. Node B).
35
Hydraulic Network Modeling with MOUSE
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.9 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 computation will have a component from change of flow direction. If
the two nodes surrounding the device are placed exactly at the same loca-
tion then the computational engine cannot determine the direction of the
flow from the coordinates of the nodes and a default direction will be
applied. This may unintentionally introduce a change in direction and
therefore also an unexpected head loss.
36 MIKE URBAN
Weirs
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.
37
Hydraulic Network Modeling with MOUSE
Table 3.6
38 MIKE URBAN
Weirs
Model data
Table 3.7
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” AngleNo
is not speci-
is a 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 Source-
mented is
Channel channel of the weir LinkID
chosen
39
Hydraulic Network Modeling with MOUSE
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.9 The difference between real world orifice and model configuration of
orifice
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.9 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.
40 MIKE URBAN
Orifices
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 computation will have a component from change of flow direction. If
the two nodes surrounding the device are placed exactly at the same loca-
tion then the computational engine cannot determine the direction of the
flow from the coordinates of the nodes and a default direction will be
applied. This may unintentionally introduce a change in direction and
therefore also an unexpected head loss.
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.
41
Hydraulic Network Modeling with MOUSE
Table 3.8
42 MIKE URBAN
Orifices
Table 3.8
Model data
Table 3.9
43
Hydraulic Network Modeling with MOUSE
44 MIKE URBAN
Stormwater Inlets
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).
45
Hydraulic Network Modeling with MOUSE
z 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).
z ON-GRADE Type, where flow in the overland flow network can con-
tinue past the connection node. Transfer capacity of the connection is
dependent 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).
46 MIKE URBAN
Stormwater Inlets
z 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 network. As with weir flow, a crest level at least 0.01 m higher
than the connection node invert level is recommended for initial condi-
tion stability.
z 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 rela-
tionship.
z 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.
z 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
defining the transfer capacity of the connection.
47
Hydraulic Network Modeling with MOUSE
z Number of Curb Inlets, allowing multiple curb inlets of the same spec-
ified geometry (transfer capacity) applied at the same location within a
single connection.
z 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 Free-
board, including reverse flow (surcharge).
There are no limitations on the number of curb inlets specified at one loca-
tion; however, the connectivity must be ‘From’ a node in the overland
flow network ‘To’ a node in the pipe network, for correct automatic calcu-
lation of slope. Note: Link slopes must be calculated in the link dialog for
automatic calculation of slope to operate.
For calculated or user defined slopes in the Curb Inlet dialog that are out-
side the range of slopes specified in the On-grade Capture dialog, the clos-
est slope curve will be used. For intermediate calculated or user defined
slopes (lying between slope curves in the On-grade Capture dialog), linear
interpolation is applied.
48 MIKE URBAN
Pumps
The DQ relation specifies the depth based capacity curve for a SAG Type
Curb Inlet. Values must be monotonously increasing in depth and dis-
charge 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.
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 over-
land flow’.
49
Hydraulic Network Modeling with MOUSE
nection between the nodes. The pump will then be placed between the two
nodes as the flow connection.
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.9 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 computation will have a component from change of flow direction. If
the two nodes surrounding the device are placed exactly at the same loca-
tion then the computational engine cannot determine the direction of the
flow from the coordinates of the nodes and a default direction will be
applied. This may unintentionally introduce a change in direction and
therefore also an unexpected head loss.
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.
50 MIKE URBAN
Pumps
difference dH between the pump wet well and the receiving point in the
model. The pump curves for this type of pumps will in general be as
shown in the figure below.
Pump curves are defined in MIKE URBAN by tabulating the curve (see
more in section 6). 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.
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 sim-
ulation 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.
51
Hydraulic Network Modeling with MOUSE
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 infor-
mation 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 corresponding to the minimum speed and a pump curve correspond-
ing to the maximum speed.
For variable speed pumps the two pump curves are defined as 'RPMmin'
and 'RPMmax' pump curves. For constant speed pumps only the 'RPM-
max' pump curve is applied.
52 MIKE URBAN
Pumps
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.
Figure 3.21
53
Hydraulic Network Modeling with MOUSE
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 con-
sequence the water level in the wet well will rise above the defined set-
point value.
Figure 3.22
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.23
More detailed output information may be added to the standard result file
for variable speed pumps. Please refer to the documentation on the DHI-
APP.INI file found in the “DHIAPP INI and ADP Files Reference Man-
ual”.
54 MIKE URBAN
Pumps
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.
55
Hydraulic Network Modeling with MOUSE
56 MIKE URBAN
Pumps
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.10
57
Hydraulic Network Modeling with MOUSE
Model data
Table 3.11
58 MIKE URBAN
Valves
3.9 Valves
A valve is a functional relation which connects two nodes of a MOUSE
network.
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.9 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 computation will have a component from change of flow direction. If
59
Hydraulic Network Modeling with MOUSE
the two nodes surrounding the device are placed exactly at the same loca-
tion then the computational engine cannot determine the direction of the
flow from the coordinates of the nodes and a default direction will be
applied. This may unintentionally introduce a change in direction and
therefore also an unexpected head loss.
The valve is by default static, in which case the valve opening must be
specified. 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.12 Identification and connectivity edit fields of the MIKE URBAN valve
editor.
60 MIKE URBAN
Valves
Table 3.13 Geometrical and hydraulic properties, edit fields of the MIKE
URBAN valve editor
61
Hydraulic Network Modeling with MOUSE
z 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.
z Topography: The "Channel Topography" dialog is an editor for
MOUSE topography data. A topography is made up from a number of
cross sections,specified in the cross section editor.
The CRS and Topography Editors are reached by selecting MOUSE |
CRS & Topography.
62 MIKE URBAN
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
section. The markers (see Figure 3.28) are used for splitting the cross sec-
tion 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
63
Hydraulic Network Modeling with MOUSE
cross section into three parallel channels. Thus the reach of the left flood
plain is given by markers (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.
Figure 3.28 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 pro-
vides the Preissmann slot, in order to facilitate the flow computations in
pressurised conditions.
64 MIKE URBAN
Emptying Storage Nodes
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 speci-
fied at a chainage lower than the length of the pipe the last CRS will be
extrapolated to the remainder of the link.
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
controlled by the actual hydraulic situation in the system and the weirs or
pumps capacities. Weirs and pumps behave in this case as if they dis-
charge "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
controlled 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
65
Hydraulic Network Modeling with MOUSE
water level in a 'control node'. The flow according to the emptying func-
tion is possible only until the storage node is completely empty.
Table 3.14
66 MIKE URBAN
Terms and Concepts
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.
67
Rainfall-Runoff Modelling with MOUSE
68 MIKE URBAN
Terms and Concepts
relevance for the hydrological modeling will be repeated and further elab-
orated.
Details about the principles and methods for catchment connections can
be found in the chapter 13.4.
z 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
hydrographs 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 relatively 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
provide realistic results in dominantly rural catchments and for long-
term analyses involving multi-event rainfall series.
69
Rainfall-Runoff Modelling with MOUSE
70 MIKE URBAN
Time-Area Method (A)
71
Rainfall-Runoff Modelling with MOUSE
72 MIKE URBAN
Time-Area Method (A)
73
Rainfall-Runoff Modelling with MOUSE
74 MIKE URBAN
Time-Area Method (A)
75
Rainfall-Runoff Modelling with MOUSE
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).
Figure 4.6 The Curves&Relations editor with one correct user defined time-
area curve
76 MIKE URBAN
Kinematic Wave (B)
77
Rainfall-Runoff Modelling with MOUSE
78 MIKE URBAN
Kinematic Wave (B)
Figure 4.8 The Kinematic Wave Parameter Sets Editor (MOUSE|Runoff Mod-
els|Parameters Kinematic Wave)
79
Rainfall-Runoff Modelling with MOUSE
80 MIKE URBAN
Linear Reservoir (C1 and C2)
MIKE URBAN includes this model in two versions: C1 and C2. These are
variants of the same model, and are used as national standards in the Neth-
erlands and in France, respectively.
81
Rainfall-Runoff Modelling with MOUSE
82 MIKE URBAN
Linear Reservoir (C1 and C2)
83
Rainfall-Runoff Modelling with MOUSE
Figure 4.11 The Linear Reservoir Parameter Sets Editor (MOUSE|Runoff Mod-
els|Parameters Linear Reservoir) - Dutch version
84 MIKE URBAN
Linear Reservoir (C1 and C2)
Figure 4.12 The Linear Reservoir Parameter Sets Editor (MOUSE|Runoff Mod-
els|Parameters Linear Reservoir) - French version
85
Rainfall-Runoff Modelling with MOUSE
86 MIKE URBAN
Unit Hydrograph Method (UHM)
87
Rainfall-Runoff Modelling with MOUSE
88 MIKE URBAN
Unit Hydrograph Method (UHM)
89
Rainfall-Runoff Modelling with MOUSE
When studying the real flow conditions in sewer systems, flow peaks dur-
ing rain events are often found to exceed the values that can be attributed
to the contribution from participating impervious areas. This is a conse-
quence 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 precipitation, but is heavily affected by the
actual hydrological situation, i.e. the "memory" from earlier hydrological
events. For a certain rainfall event, the increase in flow will therefore dif-
fer, depending on hydrological events during the previous period. The
Rainfall Induced Infiltration is also distinguished by a slow flow response,
which takes place during several days after the rainfall event.
90 MIKE URBAN
Additional Flow and RDI
When combined with any of the MOUSE surface runoff models, MOUSE
RDI provides a platform for accurate and reliable computation of urban
runoff, liberated from the limitations inherent to "standard" urban runoff
modeling.
91
Rainfall-Runoff Modelling with MOUSE
Figure 4.14 Different catchment response under the same rainfall, due to differ-
ent soil moisture conditions at the beginning of the rainfall
92 MIKE URBAN
Additional Flow and RDI
93
Rainfall-Runoff Modelling with MOUSE
Figure 4.16 The Additional Flow and RDI Model Editor (MOUSE|Runoff Mod-
els|Additional Flow and RDI
94 MIKE URBAN
Additional Flow and RDI
Table 4.8 Overview of the Additional Flow and RDI Database Attributes (Table
‘msm_HModCRC’)
95
Rainfall-Runoff Modelling with MOUSE
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’)
96 MIKE URBAN
Additional Flow and 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’)
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-transpiration is provided as
98 MIKE URBAN
Rainfall Data and Other Meteorological Variables - Boundary Conditions for Hydrological models
MIKE URBAN MOUSE can utilize time series data stored in DHI propri-
etary format DFS0. Please refer to Chapter 5 for importing, reviewing and
editing of time series data.
All three categories of meteorological variables act as boundary condi-
tions for the hydrological model and as such must be specified through the
MOUSE model boundary system (MOUSE|Boundary Condition|Catch-
ment Loads). As with any other boundary condition, this is done in two
distinct steps:
The current boundary will apply to the catchments according to the speci-
fied "Connection Type": All, List or Individual. If "All" is selected, the
data source' (e.g. rain gauge or meteorological station) co-ordinates can
optionally be specified.
The co-ordinates are of relevance only in cases where two or more bound-
aries of the same type have been specified with connection type "All":
Instead of applying the boundary condition to the entire model, the system
will apply the specified boundary conditions using the principle of geo-
graphical proximity to the catchment's centrepoint.
Some data sources (e.g. dfs0 time series) include the information about the
data source co-ordinates. By activating the data source co-ordinates group
99
Rainfall-Runoff Modelling with MOUSE
on this dialog, the specified co-ordinates would overwrite the original co-
ordinates from the time series file.
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
must include at least one item. By the nature of the meteorological bound-
aries, they can include only the "default" item types.
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" edi-
tor 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.
Per default, the specified boundary will be applied throughout the entire
duration of the simulation or during the period covered by the time series.
Optionally, 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.
On the "General" TAB, the user must specify the essential information
about the simulation (see example in Figure 4.19), and the runoff-specific
parameters on the "Runoff parameters" TAB (see example in Figure 4.20).
101
Rainfall-Runoff Modelling with MOUSE
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.
103
Rainfall-Runoff Modelling with MOUSE
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 spec-
ified. On the contrary, the Surface Runoff simulation is carried-out only
when raining. Thus, the start time for the Surface Runoff calculation is set
as the start time for rain hydrograph. The Surface Runoff calculation con-
tinues until all the surface runoff hydrographs are regressed.
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 recon-
struction of the surface runoff component (FRC).
z *.CRF file, containing maximum five time series for each sub-catch-
ment, namely:
– discharge, calculated with the Surface Runoff Model (the FRC
component),
– discharge, calculated with the RDI model (the SRC component),
– total discharge,
– variation of water content in the surface storage for the Surface
Runoff Model,
– variation of water content in the snow storage for the Surface Run-
off Model.
The *.CRF file is used as input data for a MOUSE Hydrodynamic calcula-
tion.
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.
The available period of the measured discharge data and its resolution in
time are of major importance for the credibility of the obtained parameter
values. Ideally, for a good accuracy, a 3-5 years long time series of meas-
ured discharge data with daily values is required for the calibration of the
RDI parameters. Several months long time series with higher resolution,
i.e. minutes or hours, depending on the size of the area, are needed for the
calibration 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 hydro-
logical situations, i.e. typical wet period or dry period.
105
Rainfall-Runoff Modelling with MOUSE
and controlling 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.
Numerical 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.
To separate the Afrc component (Surface Runoff Model) and the fast part
of the SRC component (Surface Runoff Component in RDI), measured
discharge data with fairly high resolution in time (hours) is required.
107
Rainfall-Runoff Modelling with MOUSE
z The first step in the RDI calibration is usually to adjust the water bal-
ance in the system, i.e. the accuracy between the calculated and meas-
ured total volume during the observed period. This is done by
correcting the proportion of area, Asrc. An increase of Asrc propor-
tionally 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).
z 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
baseflow is done during and after wet periods with low evaporation.
Adjustment against regression is done during the start of dry periods
with high evaporation, preferably when baseflow is the only flow com-
ponent.
An adjustment of CKBF does not influence the size of the discharged
volume studied for a longer period, but displaces the volumes in time.
z 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).
Umax can be estimated from how much rainfall is required for filling-
up the surface storage, i.e. generating overland flow. The same meth-
odology can also be used for the periods with low potential evapora-
tion, but only if the rain event is preceded by a long dry period.
TG is adjusted during periods with heavy filling of the root zone stor-
age, preferably in combination with low potential evaporation and pre-
ceded 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.
109
Rainfall-Runoff Modelling with MOUSE
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-exist-
ent. Discharge from the groundwater storage exists continuously all
year around. Drawing of the surface storage is faster during summer
period since the evaporation is high, and is therefore the dominating
effect on the surface storage. During periods with low evaporation,
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 transpira-
tion, 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 calibra-
tion procedure.
111
Rainfall-Runoff Modelling with MOUSE
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.
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.
The items can be created and edited in the "Boundary Items" editor. This
can be open by "Items" button on the "Network Loads" editor (recom-
mended), or through the menu at MOUSE|Boundary Conditions|Items.
If opened from the "Network Loads" editor, the "Boundary Items" editor
displays only the items associated with the current load. Initially, for a
newly declared load, the list of items will be empty.
113
Rainfall-Runoff Modelling with MOUSE
Per default, the specified runoff result file will load the network through-
out its entire duration, corresponding to the runoff computation Start and
End times. Optionally, the user can limit the temporal validity of the spec-
ified file to any period within the file's duration, by specifying the tempo-
ral 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.
5 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 varia-
tions 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 inter-
val 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.
115
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
available.
Field Description
Title Title of the time series object
Auto Conversion When tick mark is set then the values that are in ‘Delete
of 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
without the delete value typed in being converted.
Delete Value Threshold for values that are so small that they should
be disregarded.
117
Time Series
Field Description
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
Time Type A time definition can be equidistant, meaning that the
time values difference between consecutive Timesteps
is constant through the entire time or non-equidistant,
meaning that the time values difference between con-
secutive Timesteps is not constant through the entire
time. Besides these two classifications, a time defini-
tion can be Calendar type, meaning that the time values
are specified as dates or Relative type, meaning that
time values are specified as numbers.
Field Description
Time Step Value Specifies the time step value within the TS object. If
you specify the time steps in the fields in the ‘TimeStep
Information’ group box then the time step value will be
computed automatically in the unit specified in the
‘EUM Unit’ field.
Year, Month, Day, Specification of time step. If e.g. a time step of 1 day
Hour, Minute, and 1 hour is appropriate then type ‘1’ in the ‘Day’ field
Second and ‘1’ in the ’Hour’ field
119
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
121
Time Series
Field Description
Value Type Instantaneous: means that the values are representa-
tive at one precise instant.
Reverse_Mean_Step_Accumulated: is equal to
Mean_Step_Accumulated type, but values represent the
timespan between the current Timestep and the next
Timestep. Used for forecasting purposes.
Auto Conversion If tick mark is set then the values that are entered will
of Data Values be converted on the fly if the unit specified in the ‘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 without the values typed
in being converted.
Field Description
Item Type The type of time series that is entered. E.g.if you choose
Water Level the units in the ‘Unit’ field will only dis-
play units that are relevant for that type of time series. If
‘undefined’ is chosen all the units displayed in ‘Unit’
will not be filtered but will display all possible units.
Unit Displays units relevant for the time series type chosen
in ‘Unit’
x,y The x and y co-ordinates used to connect a time series
of rainfall, evaporation and temperature to a catchment.
If x and y co-ordinates are specified in the "Catchment
loads and meterological... " dialog- then these co-ordi-
nates will overwrite the coordinates typed in with the
time series.
z, Alpha, Phi, Not used in MIKE URBAN
Theta
123
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.
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.
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 corre-
spondent cell value in Timeseries Tabular changes accordingly.
125
Time Series
Select Points
When clicking in a data value that data value is selected. The current
Timestep is set to the selected data value timestep. A red square around
the data value appears and the correspondent cell in the Timeseries Tabu-
lar (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
bookmark for the data value
3 Go to Edit|Time Series
4 Right-click on the ‘Untitled’ and choose TSObject Properties, see
Figure 5.11
127
Time Series
Figure 5.11 Right-Click to set the properties for the time series
5 On the dialog that comes up your are able to specify the various prop-
erties. The properties for this specific example are seen in Figure 5.12.
In this example we are specifying that the rain starts at the 1st of Novem-
ber 2004, that we would like 20 time steps that are a minute apart.
129
Time Series
Figure 5.13 Example of how it looks after specifying the TSObject Properties
7 After the specification of the TSItem properties, you can start to type in
the various values. A result of this is seen in Figure 5.15
131
Time Series
Figure 5.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 5.17
133
Time Series
Figure 5.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 5.19.
Figure 5.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 descrip-
tion tab page, see Figure 5.20.
135
Time Series
Figure 5.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 message 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 5.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 5.22.
You can change the properties for the display of the time series. See
Figure 5.23.
137
Time Series
Figure 5.23 In the plot the graphic settings have also been changed to make it
easier to see the difference.
6.1 Introduction
In MOUSE|Curves and Relations (Figure 6.1) a number of tabular data
used in other data dialogs are specified.
139
Curves and Relations
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
6.1.10.
6.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
specify 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
specified in the dialog "Emptying Storage Nodes". An emptying function
is a QH relation, where 'Q' is the flow from the storage node into a 'receiv-
ing 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 com-
pletely empty.
141
Curves and Relations
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 val-
ues 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)).
The H-column for the basin geometry can start at any value, e.g. 0.0 for
interpretation of H as depth in the basin. The MOUSE Engine will associ-
ate the first H-value to the bottom level of the node. This means that the
same geometry 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 open-
ing (%) and resistance (k).
143
Curves and Relations
After the simulation with a ‘QdH & Power’ the summary will contain
information on the power consumption during the simulation period.
Figure 6.4 Snapshot from the html-summary with example of output when
using the QdH & Power capacity curve
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 undefined type). The tabular data set should then hold a table
with concentrations versus normalised flow, i.e. the runoff divided by the
total catchment area. See more in section 13.6.
145
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 lev-
els, air temperature, evapo-transpiration, etc.
147
Boundary Conditions
On the top level, there are three groups of boundary condtions (see fig-
ure….below). Each of these three groups is served by a separate editor,
described in full detail further below.
Figure 7.1 Basic frame of MOUSE Boundary Conditions with three groups
z Air temperature
z Evapo-transpiration
z Rainfall
z Flow - including WQ, sediment and temperature properties
An example illustrating the catchment boundary conditions is presented in
below figure
Note that the first three boundary conditions are used by a runoff-precipi-
tation model, while the water load boundary condition can either be
appended to the computed runoff (only flow) or passed directly to the net-
work connected to the catchments.
149
Boundary Conditions
z Water level
z 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 pollutant, sediment or temperature item (constant, cyclic or a time
series).
151
Boundary Conditions
different time period. The default items (i.e. storm runoff itself) are associ-
ated with temperature and pollutants BOD and NH4.
Only one water level boundary condition is allowed at each network out-
let. 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.
153
Boundary Conditions
155
Boundary Conditions
The following describes how to define an boundary condition for the net-
work. This could e.g. be a discharge entering a manhole or pipe infiltra-
tion.
"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.
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
infiltration 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.
Applying a positive value will result in water being added to the node.
157
Boundary Conditions
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 dia-
logs 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 boundary 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 pollut-
ant 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
actually loaded into the model. Per default, the scaling factor is set to 1.
Temporal variations
Time series
Time series may come from various sources, including previous MOUSE
result file (this latter option includes MOUSE files CRF, SMR and PMR).
Temporal validity
159
Boundary Conditions
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 valid-
ity (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.
7.4 Examples
In the following paragraphs some examples on setting up boundary condi-
tions are shown step-by-step
z Choose DFS0 as the Source Type, click the list button next to the Path
field and choose the dfs0 file with the water level (item ID is only nec-
essary to fill out if the dfs0 file contains more than one item, i.e. time
series)
After completing this, the time-varying water level has been successfully
added to the outlet.
161
Boundary Conditions
z Choose Time Series/Result File in order to choose a time series for the
rain event
z This activates the Time Series part of the dialog
z 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!
7.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 exam-
ple above). Alternatively, a 'shortcut' specification can be used.
163
Boundary Conditions
z Now press the Items button. This will open the Boundary Items dialog
where the runoff result file to be used can be specified
z 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.
z After having pressed Insert an item ID will automatically be given
(hint: you can change this ID if you wish to).
z Specify that the load type is 'Storm Runoff'. This will insure a correct
categorisation of this load in the simulation summary statistics and in
the job-list generation for LTS simulations.
Now the temporal variation can be specified.
z Choose Time Series/Result File in order to choose the result file from
the catchment simulation
z 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 com-
putation.
Now the runoff results are used as input for the network computation!
165
Boundary Conditions
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 is to allow this load be directly applied to the network
model.
167
Boundary Conditions
Three checkboxes 'Week Days', 'Dates' and 'Months' are available for the
activation of three different, mutually inclusive modes of the calendar def-
inition. This means that if any day is to be included in the current calendar
definition, 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 calendar 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 unconditionally used.
169
Boundary Conditions
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.
171
Boundary Conditions
8 MOUSE SIMULATIONS
MOUSE simulations may be started from ‘Simulation | Run MOUSE’.
Simulation options for MOUSE are also defined here.
z General
z Runoff parameters
z Network parameters
z Network Summary
z 2D overland parameters
The ‘Simulation ID’ and ‘Scenario ID’ are found on all four tabsheets as a
common reference.
173
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)
175
MOUSE Simulations
Figure 8.2 Example of the dialog that comes up, when pressing the ‘Info’ but-
ton
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.
177
MOUSE Simulations
179
MOUSE Simulations
z Identification of simulation
z Simulation mode
z LTS Job list
z Runoff Input
z Dynamic simulation
z Hot Start
z Processes
z Results
181
MOUSE Simulations
z Identification of simulation
z Tables with min, max and accumulated values for
z Summary of input data
z Use selection to reduce summary tables
183
MOUSE Simulations
185
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 frequency 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
separate tables. By pressing the appropriate button a specific editor can be
accessed.
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.
187
MOUSE Simulations
chosen e.g. water level. In figure Figure 8.8 is displayed the Editor for
‘Nodes result selections’
The 'Mark List' option helps clearing or checking all items on the dis-
played 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
189
MOUSE Simulations
9 2D OVERLAND FLOW
9.1 Introduction
Surface flooding can be simulated with MIKE URBAN using the follow-
ing two approaches:
191
2D Overland Flow
1 Choosing 2D model
Select MOUSE|2D Overland... to open the dialog shown in Figure 9.1.
The underlying 2D overland flow model, MIKE 21, offers the following
different types of 2D model solvers:
z Single Grid using a rectangular cell solver. The 2D model area is cov-
ered by square calculation cells.
z 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.
z Flexible Mesh solver. The 2D model area is covered by a grid consist-
ing of a combination of triangular and quadrangular cells. Currently,
this option is not supported through the MIKE URBAN user interface.
For further information on the pros and cons of each of the 2D solvers,
please refer to the MIKE 21 documentation.
z Select the desired raster layer as the DEM to be applied in the model.
z 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 rec-
tangle’, ‘Set as extent of data set’ or ‘Set as extent of raster’ can be
used to set the coordinates corresponding to the rectangular extent of
either a manually digitized polygon, a MIKE URBAN data set or a
raster.
193
2D Overland Flow
z Enter the desired 2D model resolution. The model resolution does not
have to be the same as the resolution of the DEM. However, if chosen
to be the same and if the lower left hand corner of the model area coin-
cides with the lower left hand corner of a cell in the DEM, then the
DEM data used in the 2D model will be an exact copy of the input
DEM data. Otherwise, the DEM data used in the 2D model will be
interpolated from the DEM data provided as input.
4 Setting 2D model parameters
Select MOUSE|2D Overland... followed by clicking on the tab ‘2D model
parameters’ to open the dialog shown in Figure 9.2.
This dialog allows for the editing of the following 2D model parameters
195
2D Overland Flow
z Coupling parameters…
Selecting this will open the dialog shown in Figure 9.4.
Figure 9.4 Dialog for parameters for the 2D overland flow tools
When coupling a pipe model node to the 2D model, a coupling area cover-
ing the spatial extent of the 2D cells involved is the coupling has to be
defined. The coupling area will automatically be defined by the tool "Cou-
ple 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.
z Couple nodes
This tool applies to manholes, basins and outlets that are not already cou-
pled. 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:
1 The node will be marked as coupled and it will show-up as such in the
tab ‘2D Overland’ in the dialog MOUSE|Nodes.
2 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.
3 A line between the node and the centre of the coupling area will be de-
fined 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.
197
2D Overland Flow
Figure 9.5 below 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 loca-
tion can be manually edited using the tools in the Feature Edit toolbar.
z Un-couple nodes
If any coupled nodes are selected this tool will apply only to those - other-
wise to all coupled nodes. The tool will un-couple nodes by taking the fol-
lowing actions:
1 The node will be marked as un-coupled and show-up as such in the tab
‘2D Overland’ in the dialog MOUSE|Nodes.
2 The associated coupling area polygon will be deleted.
3 The associated coupling connection will be deleted.
After clicking this menu item, the MIKE 21 j,k coordinates at the position
of the cursor, when pointing in the map view, will appear in the status bar
alongside with the X,Y coordinates and the mapscale. See Figure 9.6.
199
2D Overland Flow
This tool will prompt the user for a j,k coordinate set and subsequently the
map view will zoom to the specified location with a zoom level such that
approximately 10x10 2D cells are visible.
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.
Inlet area
The inlet area is only used when applying the orifice equation for describ-
ing the flow exchange between the MOUSE element 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. When coupling to a basin, this value is neglected.
Qdh factor
Orifice coefficient
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
201
2D Overland Flow
certain water level difference between the MOUSE and the MIKE 21
model.
Exponent
9.3.4 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.
203
2D Overland Flow
Figure 9.9 Selection of time step and request for output in 2D model
In the tab ‘2D overland parameters’ 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 fol-
lowing has to be specified for each requested result file:
205
2D Overland Flow
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 button 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 infor-
mation 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 9.11.
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
207
2D Overland Flow
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 9.11.
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 adding raster data to the profile view. See Figure 9.12.
Figure 9.11 2D flood results shown as raster data and time series
Figure 9.12 2D flood results shown as water surface profile in longitudinal profile
209
2D Overland Flow
.dfs2
<Simulation ID> Raster file in dfs2 This file contains the DEM data for
format the 2D model area selected. This file
<ScenarioID> is only generated if the Multi Cell 2D
model has been chosen. Data in the
Fine file is extracted from the original
DEM specified in MIKE URBAN at
.dfs2 the requested model resolution
divided by the fine/coarse ratio.
z Pump
z Weir
z Orifice with gate
z Orifice with weir
z 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 defi-
nition 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 opera-
tional 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 lev-
els as functions of water level at strategic location in the system).
211
MOUSE Control module
In the latter case, the control decisions are derived by evaluating (compar-
ing) 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 “continuous control problem” equation and is usually termed as
PID (Proportional-Integral-Differential) control. The actuation signal for
the regulator is generated by a PID controller, which usually appears as
part of the operational strategy programmed in a Programmable Logical
Controller (PLC).
213
MOUSE Control module
The global control can be extended also to include forecast data in addi-
tion to real-time data, which is then called predictive RTC. The most
comprehensive way to obtain forecast data is to include a model in the
control system. 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).
z 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.
z Sensors are specified as operational devices with definition of sensor
type and position in the MOUSE network. Sensors with multiple func-
tionality must be specified individually.
z 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.
z The actual controllers are not specified explicitly as physical devices,
but their function (i.e. operational logic as a combination of operational
conditions and control actions) is associated with the respective
devices.
z MOUSE RTC uses sampling and actuation (control loop) frequency
identical to the simulation time step.
z Sensor readings are simulated as perfectly accurate and with 100%
availability.
z Low-level logic of the pump START/STOP operation is built into the
program and is controlled by the START and STOP levels.
z 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.
10.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 locality 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.
215
MOUSE Control module
Component
For the "Concentration" sensor type, a measured WQ component must be
specified.
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
evaluation 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, dis-
charge 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.
217
MOUSE Control module
Action ID
Each control action needs a unique ID, which can be used to access the
actions from other dialogs.
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.
219
MOUSE Control module
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.
-setpoint downstream:1.0
-setpoint upstream:-1.0
Weirs
-setpoint downstream:-1.0
-setpoint upstream:1.0
1.0
Pumps:1.0
1.0
Figure 10.5, Figure 10.6 and Figure 10.7 show examples of how the actual
variable (flow or water level) can fluctuate around the set-point as conse-
quence of various choices of the PID constants. Each figure has three dif-
ferent graphs depending on whether the constant is too high, too low or
adequate.
Figure 10.5 Fluctuations around the set point depending on the size of the pro-
portionality factor, K
221
MOUSE Control module
Figure 10.6 Fluctuations around the set point depending on the size of the deri-
vation time, TD.
Figure 10.7 Fluctuations around the set point depending on the size of the inte-
gration time, Ti
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).
223
MOUSE Control module
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 10.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.
10.8.1 Pumps
225
MOUSE Control module
10.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.
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 direc-
tion.
227
MOUSE Control module
Figure 10.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
completely while the ordinary weir is always open upwards. It has no ceil-
ing 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.
10.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.
229
MOUSE Control module
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
function 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 sim-
ulation”) only ONE rule can be specified. The menu allows more func-
tions 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 ful-
filled and in the case of any violation, the simulation will not start.
Figure 10.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 Initializa-
tion step is for executing code opening, preparing or reading files needed
231
MOUSE Control module
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;
MOUSETimeString : String;
begin
NodeLevelString := Floattostr(GetNodeLevel('Node_11'));
end;
procedure EndStep;
begin
CloseFile(OutFile);
end;
exports
InitStep,
TimeStep,
EndStep;
end.
233
MOUSE Control module
The example above simply writes the water level in a node and the corre-
sponding time to a file, but the possibilities are endless.
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');
//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);
235
MOUSE Control module
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.
237
MOUSE Control module
Based on the specified criteria, the system automatically selects the peri-
ods for dynamic pipe flow simulations during wet weather. This reduces
the simulation time significantly, while preserving the accuracy in the sim-
ulation of relevant 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 performance can be reliably tested and evaluated based on
the computed operational statistics. This allows the user to develop the
optimal rehabilitation / upgrade strategy, e.g. for satisfying the require-
ments of environmental regulations.
239
Long term statistics
LTS is able to answer these and similar questions and thereby appears as
an essential simulation tool for a reliable assessment of system perform-
ance.
However, use of dynamic simulation for the hydraulic and pollution trans-
port 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 simulated 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 peri-
ods from the simulated series, and to limit the dynamic simulations to a
series of disconnected events, instead of the entire simulation period
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 modi-
fied according to specific requirements
The job list criteria represent actually various thresholds for the system
loading, which are used to evaluate if the state of the system is within or
outside the "dynamic" conditions. There are separate START criteria,
which determine 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 crite-
ria 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
system, 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 potential loss of important information due to the time offset
between the loads and effects in the sewer network.
241
Long term statistics
Table 11.1
Table 11.1
If the option with hot start files is selected, the system is initialized
according to the specified hotstart information - hot start filename and hot
start date - for each event in the job list. The hot start file for an actual
event is selected according to the closest fit of the actual total loads
defined as 'runoff' and 'other', but excluding loads defined as dry weather
flow (DWF), at the start time of the event. Note that the runoff input
defined through the computation editor is also considered as 'runoff' load
and as such included into the definition of initial conditions.
If the actual inflow is outside all specified intervals, the job will start with
an empty system
243
Long term statistics
Any hot start file must cover at least 24 hours. The hot start date is user-
specified, but the hot start time is automatically set equal to the event start
time.
The hot start files for the LTS job list are selected according to the speci-
fied range of inflows. This implies that all inflows other than DWF are
constant over the entire length of the hot start file. In order to reproduce a
cyclic pattern of diurnal DWF variation, care must be taken to eliminate
the filling phase in the beginning of the hot start file. This is achieved by
running the initial simulation over two or more days, and specifying a cor-
responding date as a hot start date.
245
Long term statistics
The simulation periods for individual events represent the minimum simu-
lation 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' for-
mat), which contains all relevant information about input files for the sim-
ulation, 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 11.3
[SIMULATION_EVENT]
Dt_MaxSec = 60
Dt_MinSec = 10
EndSect // SIMULATION_EVENT
247
Long term statistics
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
evaluation matrix. The evaluation matrix consists of an arbitrary number
of lines, 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
logical operator OR. The process continues until one line has been evalu-
ated 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
evaluation repeats.
Table 11.3
The specified runtime criteria are saved in the MOUSE LTS (*.MTF) file.
249
Long term statistics
Figure 11.6
Switching the simulation mode to 'LTS' opens a new part of the dialog,
including 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 crite-
rion.
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 different phases of the project - initially using only individual
events for test simulations (or in the calibration process), and later extend-
ing 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.
251
Long term statistics
[MouseTRAP]
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.
The main parameters which control the event definition for various types
of statistics are specified in the editor 'Statistics Specification - Global'
(see Figure….).
Also, this editor allows for toggling between the annual and monthly sta-
tistics and for requesting the statistics on system emissions (flows and pol-
lutants out of the system).
The peak values statistic operates with instantaneous peak values sepa-
rated 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 statis-
tics is associated with discontinuous variables like pump or weir dis-
charges (i.e. volumes). The statistics are performed on the events
delimited with zero flows. Thus, two successive weir or pump flows are
considered as independent only if they are separated by an inactive inter-
253
Long term statistics
val (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 parame-
ters 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 accumulated 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 corresponding exactly to the number of simulated events.
Max levels:
Total discharge:
Duration discharge:
Saves ranked total accumulated event discharge duration for weirs, pumps
and outlets.
Maximum Concentration:
255
Long term statistics
Event Loads:
Saves ranked total event-accumulated mass flux for weirs, pumps and out-
lets for all simulated pollutants.
Annual Loads:
Saves total accumulated yearly/monthly mass flux for weirs, pumps and
outlets 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
system), 'List' (a number of nodes/links specified on NSE/LSE/MUS-file),
or 'Individual' (one node/link name). If 'List' or 'Individual' is selected, the
supplementary information of selection file name or node name should be
specified.
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.
257
Long term statistics
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 cal-
culate the new diameter which will then be used in next iteration.
This means that after a number of iteration - depending on how good the
initial choice of diameters was - the module will stabilise results on the
diameters needed to fulfil the design criteria.
259
Automatic pipe design with MOUSE
In the present version input data for the Pipe Design Module are specified
in the Additional Parameters File (ADP). In order to activate the design
module the simulation is started in the usual way after specification of the
ADP file name on the Start Simulation Dialog.
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 12.2.7). The format of the DESIGN section in the
ADP file is as follows:
261
Automatic pipe design with MOUSE
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 diam-
eters. 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
connected 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 corre-
sponding 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.
263
Automatic pipe design with MOUSE
When the design acts on individual pipes, each pipe will initially be auto-
matically assigned to a design node corresponding to their immediate
upstream node. Here the upstream node is defined as the node connected
to the highest end of the pipe from geometric point of view.
If more than one pipe is connected to the same upstream node they will
automatically 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
influence too much on the new pipe diameters a lower limit for designing
is specified. 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
negative value indicate a lower limit below the design level and vice versa
for a positive value.
Through the iterations the pipe design will suggest diameters without tak-
ing 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.
Figure 12.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.
Input data for the Pipe Design Module is specified in the Additional
Parameters File (ADP). When the design module is in function the simula-
tion is started the usual way (with specification of the ADP file name on
the Start Simulation Dialog), but an additional dialog for viewing and edit-
ing the design data will be available in parallel to the MOUSE Simulation
Launcher.
In the simulation launcher two new item types can be selected from the
dialog <Charts>+<Modify>: Designs and Profiles.
265
Automatic pipe design with MOUSE
z “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
during the simulation
z “Profiles” are longitudinal profiles. All profiles which are saved as
LPF files can be selected
During the simulation it will hence be possible to see the results while
simulated, see Figure 12.6.
Figure 12.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
parameters and results from the design, and can be used to adjust parame-
ters for the design between the iterations.
ID ID of the pipe
UpStrNodeID ID of the node where design criteria are evalu-
ated
Active “True“ if the pipes should actually be designed -
can be used to exclude one or more pipe itera-
tions
267
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 "Update UND file with new diameters"
which will actually change the parameters in the model, or to "Delete
design file" which will delete results of all the performed iterations and
cause the design to start from scratch.
269
Automatic pipe design with MOUSE
271
Modelling Water Quality with MOUSE
feedback due to the change in pipe area and resistance caused by sediment
deposition. The following issues can be addressed:
SRQ can be used by itself to study water quality and sediment transport
processes 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.
Problem types and processes shows which water quality modules should
be used for a list of problem types. Figure 13.1 gives an overview of the
required modules for different pollutant and sediment types.
273
Modelling Water Quality with MOUSE
The model works with two sediment fractions, a fine and a coarse, each
fraction characterised by its mean diameter. The fine fraction typically has
a particle 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.
275
Modelling Water Quality with MOUSE
Notice that if the load generated from the washoff simulation subse-
quently is supposed to be used as input in a pipe model, then it is impera-
tive that the "Save as concentration" is left un-checked.
Data for description of the build-up of pollutants for each component each
must be specified. During dry weather flow, the concentration of dissolved
pollutants will build-up according to a user-specified linear build-up func-
tion, with a maximum 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 13.4.
277
Modelling Water Quality with MOUSE
279
Modelling Water Quality with MOUSE
and water blend type. If the specification of the initial concentration for a
certain component has been omitted, a zero concentration is automatically
applied. Water blend concentration must always be given as a number
between 0 and 100 percent, and the sum of the two blend components
must add up to 100 percent.
The initial conditions are specified for individual nodes. The initial condi-
tions in the connected conduits is calculated by linear interpolation of the
concentrations specified in the upstream and downstream nodes.
dC
------- = KC (13.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 (13.2)
where:
b = a dimensionless exponent.
If the exponent is set equal to zero, then the dispersion coefficient is con-
stant and independent of the flow velocity. The unit for the dispersion fac-
tor 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.
281
Modelling Water Quality with MOUSE
The dispersion coefficient can be given either globally or locally. The glo-
bal 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
specified locally, 'overrule' the global specification and prescribe a con-
stant dispersion coefficient of 4 m2/s.
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 concen-
trations 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 transition between the calculated concentration or transport at the
boundary before the flow reversal and the specified time series of bound-
ary concentrations. The description of K-mix is given in the "MOUSE
Pollution Transport Reference" manual.
283
Modelling Water Quality with MOUSE
Group 1 includes the aspects of the classical BOD/DO problem, e.g. the
degradation 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 met-
als are modeled similarly to the nutrients, e.g. as conservative substances.
Group 4 is bacteria arising from human and animal wastes. Three groups
of bacteria can be modeled : faecal coliform, total coliform and strepto-
cocci. 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
pollutants are modeled as sediment attached. This can be any pollutant
(BOD, heavy metals etc).
285
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
coliform bacteria and streptococci. Most pathogenic microorganisms are
usually unable to multiply or survive for extensive periods in the aquatic
Eschericia coli (E. coli) is one of the dominant species in faeces from
human and warm-blooded animals. The organism itself is normally con-
sidered 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
reaction. However, the die-off constant or decay rate is highly variable due
to interaction by environmental factors on bacterial die-off. The main fac-
tors are suspected to be: light, temperature and salinity variations. In sew-
ers, the only relevant factor among these to include is the temperature.
Other less important 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.
2
u
C1 C2 ----------- u
C3
gdm
K2 = ----------------------------------------- (13.3)
dm
287
Modelling Water Quality with MOUSE
z explicit models,
z morphological models.
The main difference between the two model types is that in the explicit
models the bed levels are fixed during the simulation. In the morphologi-
cal models the bed levels / sediment deposits are updated dynamically
during the simulation.
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
z Engelund - Hansen;
z Ackers - White;
z Engelund - Fredsøe - Deigaard;
z 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.
289
Modelling Water Quality with MOUSE
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
particles 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 val-
ues were within 0.5 and 2 times the observed values, van Rijn (1984b).
are e.g. sediment grain diameter, relative density of the sediment, critical
bed shear stress, type of sediment transport model, etc.
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 computa-
tions. 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 bound-
ary conditions for sediment transport are not required. The results of an
explicit sediment transport calculation represent the potential sediment
transport.
291
Modelling Water Quality with MOUSE
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 bed shear stress from the sediment deposits can be specified in three
ways:
z the bed shear stress is calculated from the local mean grain diameter
according to:
25,4
M = ----------------------
1⁄6
(13.4)
( 2,5d )
z 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 (13.5)
else
y sediment EXP
τ C = τ bot + ( τ top – τ bot ) × ------------------------- (13.6)
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,
293
Modelling Water Quality with MOUSE
This dialog allows for the specification of the input data required for the
simulation of graded sediment transport and sediment sorting.
The 'Transport Mode' combo-box allows for the choice whether the sedi-
ment fraction is transported with ordinary non-cohesive sediment trans-
port formulations, or the fraction is transported by the Advection-
Dispersion formulations. The transport of fine organic sediment fractions
is better described by use of an Advection-Dispersion transport formula-
tion 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 coefficient must be specified. The
morphological changes resulting from the sediments transported by the
Advection-Dispersion formulations are still included in the calculation of
the bed level, but the changes in bed levels are now calculated from
source/sink terms. For further information, refer to the “MOUSE Pollution
Transport Reference” manual.
W s –n
E = 1 – 1 + ----------------------
- (13.7)
nx ( Q ⁄ A )
where:
E = efficiency
Ws = setting velocity
n = removal coefficient
295
Modelling Water Quality with MOUSE
The main limitations to the time step in the sediment transport calculation
are the boundary conditions. Both the hydrodynamic and the morphologi-
cal boundary conditions should be well resolved. Finally the calculated
changes in the bed level should be sufficiently small in order not to gener-
ate instabilities in the hydrodynamic calculations.
297
Modelling Water Quality with MOUSE
The concentrations given as input data described above come from a con-
sideration 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 buildings, traffic area, industrial area and areas with dense traffic.The
catchment processing tool, MOUSE|Catchment Tools|Catchment Process-
ing, can be used to assign concentrations based on these thematic maps.
Figure 13.17 Using the catchment processing tool for assigning concentrations
299
Modelling Water Quality with MOUSE
Table 13.3 The edit fields for the surface runoff pollutants dialog
Table 13.3 The edit fields for the surface runoff pollutants dialog
301
Modelling Water Quality with MOUSE
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 Collection System Pollutant Transport module.
303
Modelling Water Quality with MOUSE
The treatment node relations are used from the treatment relationships, see
Figure 13.20.
Edit field
Relationship ID for the relation
A1-3, B1-3, C1-3 Coefficients used in the polynomium
Low concentration Typical low concentration of incoming pollutants
to the treatment node
High concentra- Typical high concentration of incoming pollut-
tion ants to the treatment node
Normal concentra- Typical normal concentration of incoming pollut-
tion ants to the treatment node
305
Modelling Water Quality with MOUSE
Edit field
Factor low conc Correction factor for changing the computed
removal efficiency over the time period of aver-
age time is equal to the low concentration
Facor normal conc Correction factor for changing the computed
removal efficiency over the time period of aver-
age time is equal to the normal concentration
Factor hi conc Correction factor for changing the computed
removal efficiency over the time period of aver-
age time is equal to the high concentration
R = ( A1 h 2 + B1 h + C1 ) ⋅ q 2 + ( A2 h 2 + B2 h + C2 ) ⋅ q 2 + ( A3 h 2 + B3 h + C3 ) ⋅ q 2
Where
h : hydraulic efficiency (%)
q : volume load 1e −3 (m 3 / h) / m 3
Figure 13.21 The coefficients for a number of pollutants have been determined
through MIKE 21 simulations
307
Modelling Water Quality with MOUSE
To sum up the methodology, for each treatment node and for each pollut-
ant:
309
310 MIKE URBAN
Terms and Concept
z Import of external data (e.g. GIS) into MIKE URBAN CS SWMM net-
work
z Copying network data from MIKE URBAN CS Asset network into
MIKE URBAN CS SWMM network
z Copying network data from MIKE URBAN CS MOUSE network into
MIKE URBAN CS SWMM network
The last option is frequently used in a combination with one of the previ-
ous options as means for achieving a full consistency of the SWMM
model.
z Nodes
z Conduits
z Orifices
z Pumps
z Weirs
z Outlets
14.2 Nodes
Nodes are hydraulic elements such as Junctions, Outfalls, Storage Units
and Flow Dividers used in modelling a CS SWMM storm and sewer col-
lection 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 stor-
age volume. Physically they could represent something as small as a
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
hydraulic elements into the following groups:
Junction
Storage Unit
Outfall
applying Dynamic Wave routing at least one outfall in the urban drainage
system must be defined. It is possible to select between five types of out-
falls, FREE, NORMAL, FIXED, TIDAL and TIMESERIES.
z FREE outfall bases the outfall stage on the smaller of the critical and
normal depths of flow in the connecting conduit.
z NORMAL outfall uses just the normal depth of flow.
z FIXED outfall uses a constant water elevation at all times.
z TIDAL outfall has outfall stage varying in a repeating fashion over a
tidal period.
z TIMESERIES outfall allows the outfall stage varying in any specified
manner over time.
Flow Divider
Three different types of dividers can be specified for each divider node
14.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
standard open and closed geometries. Irregular natural cross-section
shapes are also supported.
The Conduits editor organizes the related input data into the following
groups:
Minor losses are only computed for the DW, Dynamic wave routing
option. The minor losses are computed as K*v2/2g.
14.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 bot-
tom or along the side of the upstream node, and have a flap gate to prevent
backflow.
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 orifice 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:
14.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
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:
z Type1 - an off-line pump with a wet well where flow increases incre-
mentally with wet well volume
z Type2 - an in-line pump where flow increases incrementally with node
depth
z Type3 - an in-line pump where flow varies continuously with head dif-
ference between the inlet and outlet nodes
z 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.
14.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:
Four varieties of weirs (Weir Type) are available, each with a slightly dif-
ferent formula for computing flow as a function of area, discharge coeffi-
cient 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.
14.7 Outlets
Outlets are flow control devices which are typically used to control out-
flows from storage units. They are used to model special flow-head rela-
tions 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:
An outlet can have a flap gate which restricts flow to only one direction
and prevents backwater flow.
14.8 Transects
Transects refer to the geometric data that describe how bottom elevation
varies 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.
Each transect are given an unique ID. When a IRREGULAR shape is cho-
sen 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
Manning'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:
z Tidal - defines how the stage at an Outfall node variates by the hour of
the day.
z Rating - relates flow through an Outlet link to the head difference
across the outlet
z Pump - relates flow through a Pump link to the depth or volume at the
upstream node or to the head delivered by the pump. Following four
types of Pump relations are available
Type 1 - an off-line pump with a wet well where flow increases incre-
mentally with the available wet well volume.
Type 3 - an in-line pump where flow varies continuously with head dif-
ference between the inlet and outlet nodes.
The Tabular Data (Curves) editor organizes the related input data into the
following groups:
14.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:
Please note there can only be either a Condition or an Action in one line.
z time-varying rainfall
z evaporation of standing surface water
z snow accumulation and melting
z rainfall interception in depression storage
z infiltration of rainfall into upper soil layers
z percolation of infiltrated water into groundwater layers
z interflow between groundwater and the conveyance system
z 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 dif-
ferent 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.
z Basin
zSub-basin
zWatershed/Sewershed
zCatchment
345
Hydrological Modelling with SWMM
15.2 Catchments
Catchments are hydrological units of land whose topography cause sur-
face runoff to drain to a single point. The user is responsible for the delin-
eation 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.
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:
z General information
z Physical Properties
z Location
Location
Physical Properties
347
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).
The CS SWMM Routing Editor organizes the related input data for the
catchment concerning the surface runoff into following groups:
349
Hydrological Modelling with SWMM
15.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
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.
This approach is adopted from the NRCS (SCS) Curve Number method
for estimating runoff. It assumes that the total infiltration capacity of a soil
351
Hydrological Modelling with SWMM
can be found from the soil's tabulated Curve Number. During a rain event
this capacity is depleted as a function of cumulative rainfall and remaining
capacity. The input parameters for this method are the curve number, the
soil's hydraulic conductivity (used to estimate a minimum separation time
for distinct rain events), and a regeneration constant that describes the res-
toration of infiltration capacity during dry periods.
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:
Infiltration Properties
353
Hydrological Modelling with SWMM
15.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 subsurface water through cracked pipes, leaky joints,
poor manhole connections, 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:
Tabular Data
355
Hydrological Modelling with SWMM
357
Hydrological Modelling with SWMM
15.6 Aquifers
Aquifers are sub-surface groundwater areas used to model the vertical
movement of water infiltrating from the catchments which lie above them.
They also permit the infiltration of groundwater into the conveyance sys-
tem, or exfiltration of surface water from the conveyance system, depend-
ing on the hydraulic gradient that exists. The same aquifer object can be
shared by several 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 eleva-
tion, and loss rate to deep groundwater. In addition, the initial water table
elevation and initial moisture content of the unsaturated zone must be sup-
plied.
The CS SWMM Aquifer Editor organizes the related input data for the
Aquifer into following groups:
Aquifer Properties
359
Hydrological Modelling with SWMM
15.7 Groundwater
A sketch of the two-zone groundwater model that is used in EPA SWMM
is shown in Figure 15.8. The upper zone is unsaturated at a varying mois-
ture 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 volume 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
computed for the next time step.
The CS SWMM Groundwater Editor organizes the related input data for
the Groundwater into following groups:
361
Hydrological Modelling with SWMM
Groundwater Properties
15.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 operations, 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 cal-
endar date.
3. Any excess snow depth on the plowable area of the pack is redistributed
according to the removal parameters established for the pack.
363
Hydrological Modelling with SWMM
5. The amount of snow in the pack that melts to liquid water is found
using:
za heat budget equation for periods with rainfall, where melt rate
increases with increasing air temperature, wind speed, and rain-
fall intensity
za degree-day equation for periods with no rainfall, where melt rate
equals the product of a melt coefficient and the difference
between the air temperature and the pack's base melt tempera-
ture.
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 avail-
able to runoff from the pack.
7. The available snow melt is then reduced by the amount of free water
holding capacity remaining in the pack. The remaining melt is treated the
same as a direct rainfall input onto the catchment.
Snowpack Properties
365
Hydrological Modelling with SWMM
15.9 Climatology
Climate-related variables used for computing runoff and snow melt are
grouped in the Climatology object:
zTemperature
zEvaporation
zWind Speed
zSnow Melt
zAreal Depletion
Temperature data are used when simulating snowfall and snow melt
processes during runoff calculations. If these processes are not being sim-
ulated then no temperature data are needed.
367
Hydrological Modelling with SWMM
Evaporation can occur for standing water on catchment surfaces, for sub-
surface water in groundwater aquifers, and from water held in storage
units. Evaporation rates can be stated as:
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.
z Temperature Tabsheet
z Evaporation Tabsheet
z Wind Speed Tabsheet
z Snow melt Tabsheet
z Areal Depletion Tabsheet
Temperature data
369
Hydrological Modelling with SWMM
Evaporation Data
371
Hydrological Modelling with SWMM
Windspeed data
373
Hydrological Modelling with SWMM
Snowmelt data
375
Hydrological Modelling with SWMM
377
Hydrological Modelling with SWMM
15.10 Coverage
The CS SWMM Coverage Editor organizes the related input data for the
coverage into following groups:
Landuse coverage
379
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.
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 peri-
ods.
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
collected runoff and allow it more time to infiltrate the native soil beneath
it.
tration into the native soil from occurring. Infiltration trenches and porous
pavement systems can also be subjected to a decrease in hydraulic con-
ductivity over time due to clogging.
Although some LID practices can also provide significant pollutant reduc-
tion benefits, at this time SWMM only models their hydrologic perform-
ance
The CS SWMM LID Controls Editor organizes the data input for the dif-
ferent types of LID practices. The input data is organized into following
group and tabs:
z Identification
z Tabs: Surface, Soil, Pavement, Storage and Drain
381
Hydrological Modelling with SWMM
383
Hydrological Modelling with SWMM
Identification
385
Hydrological Modelling with SWMM
387
Hydrological Modelling with SWMM
389
Hydrological Modelling with SWMM
391
Hydrological Modelling with SWMM
The boundary conditions can be visualized on the map. Figure 16.1 shows
how where you can toggle the visualization on/off in the TOC. Boundary
conditions are per default unchecked (one exception: catchment connec-
tions) 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.
393
Time Series and Boundary Conditions in SWMM
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.
z Temperature data
z Evaporation data
z Rainfall data
z Water stage at outfall nodes
z External inflow hydrographs at drainage system nodes
z External inflow pollutographs at drainage system nodes.
The CS SWMM Time Series editor organizes the related input data for the
different timeseries into the following groups:
395
Time Series and Boundary Conditions in SWMM
z Monthly
z Daily
z Hourly
z Weekend
The CS SWMM Time Patterns editor organizes the related input data for
the different time patterns into the following groups:
397
Time Series and Boundary Conditions in SWMM
399
Time Series and Boundary Conditions in SWMM
16.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 (spec-
ified in Time Series as Type Rainfall) or come from an external file.
Several different popular rainfall file formats currently in use are sup-
ported, 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
volume), 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.
401
Time Series and Boundary Conditions in SWMM
The CS SWMM Raingauges editor organizes the related input data for the
rainfall data into the following groups:
The CS SWMM Dry Weather Flow editor organizes the related input data
for each dry weather inflow into the following groups:
403
Time Series and Boundary Conditions in SWMM
Table 16.6 The Dry Weather Flow Identification and Connectivity Group
Table 16.8 The Pollutant Data in upper grid attached to dry weather flow
405
Time Series and Boundary Conditions in SWMM
16.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 catch-
ments are defined). A pollutant component can only be defined as inflow
to the system 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:
407
Time Series and Boundary Conditions in SWMM
The CS-SWMM Run SWMM5 editor is split into following four tab-
sheets:
z General
z Runoff
z Network
z Summary
409
Project Options and Simulations in SWMM
Simulation Period
Start Simulation
Units selection
411
Project Options and Simulations in SWMM
Infiltration Method
413
Project Options and Simulations in SWMM
Routing Simulation
415
Project Options and Simulations in SWMM
417
Project Options and Simulations in SWMM
Summary Specifications
419
Project Options and Simulations in SWMM
Subcatchment Summary
Node Summary
Link Summary
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 conduit 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:
421
Water Quality Modelling with SWMM
The following processes can be associated with each land use category:
z Pollutant Buildup
z Pollutant Washoff
z Street Cleaning
The CS SWMM Land Use editor organizes the related input data for the
land uses into the following groups:
18.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 ) (18.1)
423
Water Quality Modelling with SWMM
Exponential Function
– C 2t
B = C1 ( 1 – e ) (18.2)
Saturation Function
C1 t
B = -------------
- (18.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
pollutant buildup into the following groups:
Buildup Properties
18.4 Washoff
Pollutant Washoff from a given land use category during wet weather peri-
ods can be described in one of the following ways:
Exponential Washoff
425
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 (18.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 (18.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
pollutant washoff into the following groups:
Washoff Properties
427
Water Quality Modelling with SWMM
18.5 Loading
Input loadings of pollutants from external and dry weather inflows are
supplied through time series data associated with particular nodes of the
collection 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
18.6 Pollutant
The generation, inflow and fate of any number of user-specified pollutants
can be modelled. Concentrations in rainfall, in groundwater, in
inflow/infiltration is modelled.
There can be used following different units for the individual user defined
pollutant units (milligrams/liter, micrograms/liter, or counts/liter)
The CS SWMM Pollutant Editor organizes the related input data for the
pollutants into the following groups:
429
Water Quality Modelling with SWMM
Pollutant Properties
431
Water Quality Modelling with SWMM
For example, a first-order decay expression for BOD exiting from a stor-
age node might be:
– 0,05 ⋅ HRT
C = BOD • e (18.6)
The CS SWMM Treatment Editor organizes the related input data for the
local treatment into the following groups:
Treatment Properties
433
Water Quality Modelling with SWMM
C
Capacity Curve . . . . . . . . . . . . 139
Conduits . . . . . . . . . . . . . . . . 320
Control action . . . . . . . . . . . . . 217
Controllable devices . . . . . . . . . 224
Cross sections . . . . . . . . . . . . . 62
F
Flooding . . . . . . . . . . . . . . . . 191
Flow Dividers . . . . . . . . . . . . . 313
H
Hydrological models . . . . . . . . . . 69
I
Inlet . . . . . . . . . . . . . . . . . . . . 44
J
Job list . . . . . . . . . . . . . . . . . 241
Junctions . . . . . . . . . . . . . . . . 312
K
Kinematic Wave . . . . . . . . . . . . . 77
L
Linear Reservoir . . . . . . . . . . . . 80
Long Term Statistics . . . . . . . . . 239
N
Nodes and Structures . . . . . . . . . 23
O
Orifice . . . . . . . . . . . . . . . . . . 40
Outfalls . . . . . . . . . . . . . . . . . 312
P
PID . . . . . . . . . . . . . . . . . . . 219
Pipe Design . . . . . . . . . . . . . . 259
Pipes and Canals . . . . . . . . . . . . 30
Pump . . . . . . . . . . . . . . . . . . . 49
1
Q
QH Relation . . . . . . . . . . . . . . 139
R
RDI . . . . . . . . . . . . . . . . . . . . 90
S
Sediment transport . . . . . . . . . . 271
Sensor . . . . . . . . . . . . . . . . . 215
Storage Units . . . . . . . . . . . . . 312
T
Time series . . . . . . . . . . . . . . 115
Time-Area Curve . . . . . . . . . . . 139
Time-Area Method . . . . . . . . . . . 71
U
Unit Hydrograph Method . . . . . . . 86
W
Weirs . . . . . . . . . . . . . . . . . . . 36
2 MIKE URBAN