Rf2D ReferenceManual
Rf2D ReferenceManual
Rf2D ReferenceManual
Reference Manual
December 2017
Hydronia LLC
RiverFlow2D
c
model and documentation produced by Hydronia, LLC, Pembroke Pines, FL. USA.
Information in this document is subject to change without notice and does not represent a commitment on part of
Hydronia, LLC. The software described in this document is furnished under a license agreement.
RiverFlow2D, OilFlow2D, RiverFlow2D, and RiverFlow2D GPU are copyrighted by Hydronia, LLC. 2011-2017.
SMS is trademark of Aquaveo, LLC.
All other products or service names mentioned herein are trademarks of their respective owners.
All rights reserved. No part of this publication may be reproduced, stored in a retrieval system or transmitted in
any form or by any means electronic, mechanical, photocopying, recording or otherwise, without the prior written
permission of Hydronia, LLC.
Last document modification date: August, 2017.
Technical Support: support@hydronia.com
Web site: www.hydronia.com
Contents
ii
5.4 Finite-Volume Numerical Solution . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
5.4.1 Numerical Optimizations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
5.4.2 Stability Region . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
5.5 Open Boundary Conditions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
5.5.1 Single Variable Boundary Condition Types (BCTYPE 1 and 6) . . . . . . . 30
5.5.2 Discharge Rating Table (BCTYPE 9) . . . . . . . . . . . . . . . . . . . . . . 31
5.5.3 “Free” Open Boundaries (BCTYPE 10, 11) . . . . . . . . . . . . . . . . . . 31
5.5.4 Uniform Flow Boundary Condition (BCTYPE 12) . . . . . . . . . . . . . . 31
5.5.5 Numerical Implementation of Open Boundaries . . . . . . . . . . . . . . . . 31
5.5.6 Closed Boundaries . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36
5.6 Dry/Wet Cell Modeling . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
5.6.1 Cell definitions Based on Dry and Wet Conditions . . . . . . . . . . . . . . 37
5.7 Volume Conservation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38
5.8 Manning’s n roughness Coefficients . . . . . . . . . . . . . . . . . . . . . . . . . . . 39
5.9 Sediment Transport Model: RiverFlow2D ST Module . . . . . . . . . . . . . . . . . 39
5.9.1 Model Equations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39
5.9.2 Sediment Transport Laws . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40
5.9.3 Boundary Conditions for the Bed Load Sediment Model . . . . . . . . . . . 42
5.9.4 Suspended Load . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 42
5.9.5 Boundary Conditions for the Suspended Sediment Model . . . . . . . . . . . 45
5.9.6 Geomorphological Collapse . . . . . . . . . . . . . . . . . . . . . . . . . . . 45
5.9.7 Sediment Transport Finite-Volume Numerical Solution . . . . . . . . . . . . 45
5.9.8 Entering Data for the Sediment Transport Model . . . . . . . . . . . . . . . 46
5.9.9 Assumptions of the Sediment Transport Model . . . . . . . . . . . . . . . . 46
5.10 Mud and Debris Flow Model: RiverFlow2D MD module . . . . . . . . . . . . . . . 47
5.10.1 Assumptions of the Mud and Debris Flow Model . . . . . . . . . . . . . . . 47
5.10.2 Rheological Formulation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 47
5.10.3 Entering Data for the Mud and Debris Flow Model . . . . . . . . . . . . . . 50
5.11 Pollutant Transport Model: RiverFlow2D PL Module . . . . . . . . . . . . . . . . 50
5.11.1 Model Equations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50
5.11.2 Pollutant Transport Finite-Volume Numerical Solution . . . . . . . . . . . . 51
5.11.3 Entering Data for the Pollutant Transport Model . . . . . . . . . . . . . . . 53
5.11.4 Assumptions of the Pollutant Transport Model . . . . . . . . . . . . . . . . 54
6 Code Parallelization . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 55
6.1 RiverFlow2D CPU . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 55
6.2 RiverFlow2D GPU . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 55
7 Hydraulic and Hydrologic Components . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57
7.1 Bridges Component . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 58
7.1.1 Bridge Geometry Data File . . . . . . . . . . . . . . . . . . . . . . . . . . . 60
7.1.2 Bridge Calculations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 62
7.2 Bridge Piers . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 66
7.2.1 Bridge Pier Calculation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 67
7.3 Culverts Component . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 68
7.3.1 Culvert Calculation using a Rating Table (CulvertType = 0) . . . . . . . . 70
7.3.2 Culvert Calculation using a Culvert Characteristics (CulvertType = 1,2) . . 70
7.3.3 Assumptions of Culvert Calculations . . . . . . . . . . . . . . . . . . . . . . 75
7.4 Gates Component . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 75
7.4.1 Gate Calculations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 78
7.5 Rainfall, Evaporation and Infiltration . . . . . . . . . . . . . . . . . . . . . . . . . . 78
7.5.1 Rainfall . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 79
7.5.2 Infiltration . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 79
7.5.3 Entering Rainfall, Evaporation and Infiltration Data . . . . . . . . . . . . . 83
7.6 Wind Component . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 84
7.7 Internal Rating Tables . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 84
7.7.1 Internal Rating Table Calculations . . . . . . . . . . . . . . . . . . . . . . . 87
7.7.2 Assumptions of Internal Rating Table Calculations . . . . . . . . . . . . . . 87
7.8 Sources and Sinks . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 87
7.9 Weirs . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 89
7.9.1 Weir Calculations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 91
7.9.2 Assumptions of Weir Calculations . . . . . . . . . . . . . . . . . . . . . . . . 92
7.10 Dam Breach Modeling . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 92
8 RiverFlow2D Data Input Program (DIP) . . . . . . . . . . . . . . . . . . . . . . . . . . 95
8.1 Control Data Panel (.DAT files) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 96
8.2 Sediment Transport Data Panel for the RiverFlow2D Plus model (.SEDS and .SEDB
Files) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 100
8.3 Mud and Debris Flow Data Panel for the RiverFlow2D Plus model (.MUD File) . 103
8.4 Oil Spill on Land for the OilFlow2D model (.OILP File) . . . . . . . . . . . . . . . 104
8.5 Graphic Output Options Tab (.PLT File) . . . . . . . . . . . . . . . . . . . . . . . 104
8.6 Profile Output Panel (.PROFILES File) . . . . . . . . . . . . . . . . . . . . . . . . 106
8.7 Cross Section Output Panel (.XSECS File) . . . . . . . . . . . . . . . . . . . . . . 107
8.8 Inflow Boundary Data Panel . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 107
8.9 Outflow Boundary Data Panel . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 109
8.10 Culverts Panel (.CULVERTS File) . . . . . . . . . . . . . . . . . . . . . . . . . . . 111
8.11 Internal Rating Tables Panel (.IRT File) . . . . . . . . . . . . . . . . . . . . . . . . 113
8.12 Weirs Panel (.WEIRS File) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 114
8.13 Sources/Sinks Panel (.SOURCES File) . . . . . . . . . . . . . . . . . . . . . . . . . 115
8.14 Open Boundary Conditions Panel (.IFL File) . . . . . . . . . . . . . . . . . . . . . 116
8.15 Bridge Piers Panel (.PIERS File) . . . . . . . . . . . . . . . . . . . . . . . . . . . . 117
8.16 Observation Points Panel (.OBS File) . . . . . . . . . . . . . . . . . . . . . . . . . 117
8.17 Tools Panel . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 119
8.17.1 Two-Dimensional Model Result Interpreter: 2dMRI . . . . . . . . . . . . . . 119
8.17.2 HEC-RAS Data Extraction Tool . . . . . . . . . . . . . . . . . . . . . . . . 119
9 Input Data File Reference . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 121
9.1 Run Control Data . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 123
9.1.1 Run Control Data File: .DAT . . . . . . . . . . . . . . . . . . . . . . . . . . 123
9.2 Mesh Data . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 130
9.2.1 Mesh Data File: .FED . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 130
9.2.2 Open Boundary Conditions Data Files: .IFL and .OBCP . . . . . . . . . . . 132
9.2.3 Mesh Boundary Data File: .TBA . . . . . . . . . . . . . . . . . . . . . . . . 137
9.3 Component data . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 138
9.3.1 Bridges . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 138
9.3.2 Bridges Data File: .BRIDGES . . . . . . . . . . . . . . . . . . . . . . . . . . 139
9.3.3 Bridge Cross Section Geometry Data File . . . . . . . . . . . . . . . . . . . 140
9.3.4 Bridge Piers Drag Forces File: .PIERS . . . . . . . . . . . . . . . . . . . . 143
9.3.5 Culverts Data File: .CULVERTS . . . . . . . . . . . . . . . . . . . . . . . . 145
9.3.6 GATES Data Files: .GATES . . . . . . . . . . . . . . . . . . . . . . . . . . 152
9.3.7 Gate Aperture Time Series File . . . . . . . . . . . . . . . . . . . . . . . . . 153
9.3.8 Internal Rating Table Data File: .IRT . . . . . . . . . . . . . . . . . . . . . 154
9.3.9 Mud and Debris Flow Data File: .MUD . . . . . . . . . . . . . . . . . . . . 155
9.3.10 Oil Spill on Land Data File: .OILP . . . . . . . . . . . . . . . . . . . . . . . 157
9.3.11 Pollutant Transport Module Data File: .SOLUTES . . . . . . . . . . . . . . 158
9.3.12 Pollutant Transport Module Initial Concentration Data File: .CINITIAL . . 160
9.3.13 Manning’s n Variable with Depth Data File: .MANNN . . . . . . . . . . . . 160
9.3.14 Rainfall And Evaporation Data File: .LRAIN . . . . . . . . . . . . . . . . . 162
9.3.15 Infiltration Data File: .LINF . . . . . . . . . . . . . . . . . . . . . . . . . . 164
9.3.16 RiverFlow2D Plus Sediment Transport Data Files: .SEDS and .SEDB . . . 166
9.3.17 Sources and Sinks Data File: .SOURCES . . . . . . . . . . . . . . . . . . . 171
9.3.18 Weirs Data File: .WEIRS . . . . . . . . . . . . . . . . . . . . . . . . . . . . 172
9.3.19 Wind Data File: .WIND . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 173
9.4 Output control data . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 175
9.4.1 Observation Points Data File: .OBS . . . . . . . . . . . . . . . . . . . . . . 175
9.4.2 Graphical Output Control Data File: .PLT . . . . . . . . . . . . . . . . . . 176
9.4.3 Data for Profile Result Output: .PROFILES . . . . . . . . . . . . . . . . . . 180
9.4.4 Cross Section Data for Result Output File: .XSECS . . . . . . . . . . . . . 180
9.5 Elevation data . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 182
9.5.1 X Y Z data with header . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 182
9.5.2 General X Y Z data . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 183
9.6 Boundary conditions data . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 183
9.6.1 One Variable Boundary Condition Files . . . . . . . . . . . . . . . . . . . . 183
9.6.2 Two Variables Boundary Condition Files . . . . . . . . . . . . . . . . . . . . 184
9.6.3 Multiple-Variable Boundary Condition Files . . . . . . . . . . . . . . . . . . 185
9.6.4 Stage-Discharge Data Files . . . . . . . . . . . . . . . . . . . . . . . . . . . 186
9.6.5 Culvert Depth-Discharge Data Files . . . . . . . . . . . . . . . . . . . . . . 187
10 Output File Reference . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 189
10.1 Output File Overview . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 189
10.2 General Output Files . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 192
10.2.1 Run Options Summary . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 192
10.2.2 Mesh Data and Mesh Metrics . . . . . . . . . . . . . . . . . . . . . . . . . . 193
10.2.3 Run Summary . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 197
10.2.4 Tabular Output Files . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 198
10.2.5 Area of Inundation Files . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 199
10.2.6 Maximum Value Tabular Files . . . . . . . . . . . . . . . . . . . . . . . . . . 200
10.2.7 Observation Point Output Files . . . . . . . . . . . . . . . . . . . . . . . . . 201
10.2.8 Hot Start File . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 202
10.3 Component Output Files . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 203
10.3.1 Culvert Output Files . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 203
10.3.2 Internal Rating Table Output Files . . . . . . . . . . . . . . . . . . . . . . . 205
10.3.3 Weir Output Files . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 206
10.4 Profile and Cross Section Output Files . . . . . . . . . . . . . . . . . . . . . . . . . 207
10.4.1 General Cross Section Output Files . . . . . . . . . . . . . . . . . . . . . . . 207
10.4.2 Cross Section Hydrograph Output Files . . . . . . . . . . . . . . . . . . . . 208
10.4.3 Profile Output Files . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 209
10.5 Output Files for GIS Post-processing . . . . . . . . . . . . . . . . . . . . . . . . . . 211
10.5.1 General Spatial Distribution Files for SMS (.HDF5) . . . . . . . . . . . . . 211
10.5.2 General Spatial Distribution Files (.EXP) . . . . . . . . . . . . . . . . . . . 212
10.5.3 Pollutant Concentration Files . . . . . . . . . . . . . . . . . . . . . . . . . . 214
10.5.4 Suspended Sediment Concentration Files . . . . . . . . . . . . . . . . . . . . 215
10.5.5 Maximum Value GIS Files . . . . . . . . . . . . . . . . . . . . . . . . . . . . 216
10.5.6 Time-to-Depth GIS File . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 216
10.5.7 Support Files to Generate ESRI Shapefiles . . . . . . . . . . . . . . . . . . . 217
10.5.8 Support Files to Generate USBR Hazard Maps in ESRI Shapefiles format
(only for Argus ONE Post Processing) . . . . . . . . . . . . . . . . . . . . . 217
10.6 Output Files for QGIS Post-processing . . . . . . . . . . . . . . . . . . . . . . . . . 218
10.6.1 General Results at Cells for QGIS . . . . . . . . . . . . . . . . . . . . . . . . 218
10.6.2 Maximum Values at Cells QGIS Files . . . . . . . . . . . . . . . . . . . . . . 219
10.6.3 Time-to-Depth at Cells QGIS File . . . . . . . . . . . . . . . . . . . . . . . 219
10.6.4 Hazard Intensity Values at Cells QGIS Files . . . . . . . . . . . . . . . . . . 220
10.7 Graphic Output Files . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 220
10.7.1 DXF Files . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 220
10.7.2 GIF Output Files . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 220
10.7.3 VTK Output Files . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 221
11 Data Input Program Tools . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 223
11.1 HEC-RAS Data Extraction Tool . . . . . . . . . . . . . . . . . . . . . . . . . . . . 223
12 References . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 225
List of Figures
2.1 SMS executable Properties dialog with the Run as administrator option set. . . . . 8
2.2 Location of the RiverFlow2D Documentation folder. . . . . . . . . . . . . . . . . . 9
vii
7.5 Bridge polyline on mesh. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 60
7.6 Front view of a bridge cross section. . . . . . . . . . . . . . . . . . . . . . . . . . . 61
7.7 Top view of a bridge showing the cross sections of interest. Only two piers are
depicted for simplicity. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 63
7.8 Simple example of A1 ,A2 , A3 and A4 used to calculate head loss in free surface bridges. 63
7.9 Simple example of A1 ,A2 , A3 and A4 used to calculate head loss in a partially
submerged bridges. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 63
7.10 Simple example of A1 ,A2 , A3 and A4 used to calculate head loss in fully submerged
bridges. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 64
7.11 Application of the scheme in triangular structured meshes. Normal bridge (left) and
oblique bridge (right). . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 64
7.12 Bridge pier proportions used to asses the influence of the structure width. . . . . . 65
7.13 Influence of the structure width on the total head change (∆H) across the bridge
as a function of the Froud number downstream. . . . . . . . . . . . . . . . . . . . . 65
7.14 RiverFlow2D Control Panel dialog with Bridge Piers component selected. . . . . . 66
7.15 Pier Data dialog. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 67
7.16 Piers inside cells. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 67
7.17 Schematic view of a rectangular pier. . . . . . . . . . . . . . . . . . . . . . . . . . . 68
7.18 RiverFlow2D Global Parameters dialog with Culverts Component selected. . . . . . 68
7.19 Schematic cut view perpendicular to a gate structure. . . . . . . . . . . . . . . . . 75
7.20 Flow modes across gates. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 76
7.21 RiverFlow2D Control Data panel with Gates Component selected. . . . . . . . . . 76
7.22 Gate polyline. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 77
7.23 Feature arc attributes. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 77
7.24 Gates dialog. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 77
7.25 Water levels for discharge under a gate in submerged conditions formulated as in
(G1). . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 78
7.26 Water levels for discharge under a gate in submerged conditions formulated as in
(G2). . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 78
7.27 RiverFlow2D Rainfall/Infiltration panel. . . . . . . . . . . . . . . . . . . . . . . . . 84
7.28 RiverFlow2D Control Data panel with IRT component selected. . . . . . . . . . . . 85
7.29 Internal Rating Table polyline. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 85
7.30 Feature arc attributes. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 86
7.31 Internal Rating Table dialog. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 86
7.32 Window to enter data pairs (discharge-WSE) for an Internal Rating Table. . . . . 86
7.33 Internal Rating Table dialog. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 87
7.34 RiverFlow2D Control Data panel with the Sources and Sinks component selected. . 88
7.35 Source/Sink dialog. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 89
7.36 Time series editor to enter discharge for a Source/Sink. . . . . . . . . . . . . . . . 89
7.37 RiverFlow2D Control Data panel with Weirs component selected. . . . . . . . . . . 90
7.38 Weir polyline. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 90
7.39 Feature arc attributes. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 91
7.40 Weirs data dialog. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 91
7.41 Coordinate system for the dam breach representation. . . . . . . . . . . . . . . . . 93
10.1 Output color coded by depth and with velocity vectors . . . . . . . . . . . . . . . . 221
7.1 Variable Descriptions for the bridge cross section geometry file. . . . . . . . . . . 62
7.2 Manning’s n roughness coefficients for various culvert materials. Adapted from
Froehlich (2003). . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 71
7.3 Entrance loss coefficients Ke . Adapted from Froehlich (2003). . . . . . . . . . . . 71
7.4 Culvert inlet control formula coefficients. Adapted from Froehlich (2003). . . . . . 72
7.5 Culvert inlet configurations. Adapted from www.xmswiki.com/xms/. . . . . . . . . 73
7.6 Horton initial infiltration for different soils. Source: Akan(1993). . . . . . . . . . . 79
7.7 Horton final infiltration for different soils. Source: Akan(1993). . . . . . . . . . . 80
7.8 Mean values and standard deviation for Green-Ampt model parameters. Source:
Rawls et al. (1983). . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 81
7.9 Antecedent Moisture Content groups (adaptation from Mishra et al. (2003) . . . 83
x
8.18 Parameters on the Outflow Boundary Data Panel. . . . . . . . . . . . . . . . . . . 110
8.19 Parameters on the Culverts Panel. . . . . . . . . . . . . . . . . . . . . . . . . . . . 111
8.20 Parameters on the Internal Rating Tables Panel. . . . . . . . . . . . . . . . . . . . 113
8.21 Parameters on the Weirs Panel. . . . . . . . . . . . . . . . . . . . . . . . . . . . . 114
8.22 Parameters on the Sources/Sinks Panel. . . . . . . . . . . . . . . . . . . . . . . . 115
8.23 Parameters on the Open Boundary Condition Panel. . . . . . . . . . . . . . . . . 116
8.24 Parameters on the Bridge Piers Panel. . . . . . . . . . . . . . . . . . . . . . . . . 117
8.25 Parameters on the Observation Points Panel. . . . . . . . . . . . . . . . . . . . . . 118
8.26 Files generated by the HEC-RAS Data Extraction Tool. . . . . . . . . . . . . . . 120
9.41 Supported image formats and their corresponding world file extensions. . . . . . . 179
9.42 Variable Descriptions for the .PROFILES File . . . . . . . . . . . . . . . . . . . . 180
9.43 Variable Descriptions for the .XSECS File . . . . . . . . . . . . . . . . . . . . . . 181
9.44 Variable Descriptions for the .EXP File. . . . . . . . . . . . . . . . . . . . . . . . 182
9.45 Variable Descriptions of Boundary Condition Files. . . . . . . . . . . . . . . . . . 184
9.47 Variable Descriptions of Two-Variable Boundary Condition Files. . . . . . . . . . 185
9.49 Variable Descriptions of Multiple-Variable Boundary Condition Files. . . . . . . . 186
9.50 Variable Descriptions of Two-Variable Boundary Condition Files. . . . . . . . . . 187
9.51 Variable Descriptions of Culvert Depth-Discharge Data Files. . . . . . . . . . . . 187
1
1.1 Summary of RiverFlow2D Features and Capabilities 2
– Triangular (TIN)
– Inverse Distance Weighting
– Sampling, and many more
• Gates
• ASCII X, Y, Z
• USGS DEM
• ESRI shapefiles
• Autodesk DXF
• Uniform flow
• Rating tables
1.1 Summary of RiverFlow2D Features and Capabilities 4
• Free outflow
• Froude number
• Sediment fluxes
• Paraview VTK
• Bed-changes (erosion-deposition)
• 10 sediment-transport formulas
Figure 1.2 – SMS main window showing a velocity field computed with RiverFlow2D.
1.1 Summary of RiverFlow2D Features and Capabilities 6
Installing RiverFlow2D SMS Edition
2
RiverFlow2D requires SMS version 11.2 or later. This section will assist you to setup the RiverFlow2D
model to enable using the it from SMS.
It is recommended to uninstall all previous RiverFlow2D SMS setups before installing this version.
1. When installing from a DVD, please insert the RiverFlow2D installation disk into your DVD
drive, double click on the RiverFlow2D SMS setup file and follow the on-screen instructions.
If you received a setup file link, please download the file and double-click on it.
2. If you are installing on a PC running Windows 7 or later, you must be logged on the PC as
an administrator before you begin the installation.
7
2.3 Activating SMS 8
4. Scroll down in the Model Executables frame till you see the Generic model.
6. Click OK to exit.
Often you may find necessary to run SMS in administration mode. To set that mode of op-
eration, assuming you have SMS v12.1, use the Windows File Explorer and go to the C:
\ProgramFiles\SMS12.164-bit folder, right-click on SMS121.exe, and then select Run as ad-
ministrator as indicated in Figure 2.1
Figure 2.1 – SMS executable Properties dialog with the Run as administrator option set.
9 2 Installing RiverFlow2D SMS Edition
Also under \RiverFlow2D_SMS, you will find example projects, videos, presentations, and other
useful resources.
11
3.1 Modules 12
• Time Steps Window: it is used to select results for specific time steps when a transient
dataset has been loaded into the project.
• Edit Window: It is used to show and/or change the coordinates of selected entities.
• Menu Bar: Allow managing data depending on the active module.
• Status Bar: it shows the mouse coordinates when the model is in plan view and the infor-
mation for selected entities.
3.1 Modules
Data is divided into modules based upon the data type they can accommodate.
3.1.6 Vectors
• Visualize vector datasets as arrows.
• Constant size or vary by magnitude.
3.1 Modules 16
3.1.7 Plots
• 2D plots to visualize results and compare to measured values.
• Profile plots view scalar data along polylines.
• Time-series plots view scalar, vector, or flux (flowrate) data at a point or across a polyline.
• Several kinds of plots can be used to compare model results with measured data.
3.1.8 Animations
• Several types of AVI animations (film loops) can be generated.
• Transient data animation shows model changes through time (contours, vectors, etc.).
• Flow trace uses vector data to generate flow paths through the geometry.
• Drogue plots use user specified starting locations and show how the particles would flow
through a vector field.
• Multiple view animations show the data while transitioning between different views.
19
20
The commands in SMS are divided based on the types of data they operate on. The dynamic
toolbar and available menu commands instantaneously change when the user switches from one
module to another. Only one module is active at any given time. The data associated with a
module (e.g. a 2D mesh) is preserved when the user switches to a different module. The different
SMS Modules used to set up a project in SMS are described below:
2. The bed shear stress is assumed to follow the depth-average velocity directions.
3. The model does not include dispersion nor turbulence terms. Turbulence dissipation and en-
ergy loses are accounted for only through the Manning’s n term in the momentum equations.
23
5.3 Unsteady Flow Model 24
∂U ∂F(U) ∂G(U)
+ + = S(U, x, y) (5.1)
∂t ∂x ∂y
T
where U = (h, qx , qy ) is the vector of conserved variables with h representing the water depth,
qx = uh and qy = vh the unit discharges, with (u, v) the depth averaged components of the velocity
vector u along the (x, y) coordinates respectively. The flux vectors are given by:
!T !T
qy2 1 qx qy qx qy qy2 1
F= qx , + gh2 , , G= qy , , + gh2 (5.2)
h 2 h h h 2
where g is the acceleration of the gravity. The terms 21 gh2 in the fluxes have been obtained after
assuming a hydrostatic pressure distribution in every water column, as usually accepted in shallow
water models. The source term vector incorporates the effect of pressure force over the bed and
the tangential forces generated by the bed stress
T
S = (0, gh(S0x − Sf x ), gh(S0y − Sf y )) (5.3)
where the bed slopes of the bottom level zb are
∂zb ∂zb
S0x = − , S0y = − (5.4)
∂x ∂y
and the bed stress contribution is modeled using the Manning friction law so that:
√ √
n2 u u2 + v 2 n2 v u2 + v 2
Sf x = , Sf y = (5.5)
h4/3 h4/3
with n the roughness coefficient.
where E = (F, G) and n = (nx , ny ) is the outward unit normal vector to the volume Ω. In order
to obtain a numerical solution of system (5.6) the domain is divided into computational cells, Ωi ,
using a fixed mesh. Assuming a piecewise representation of the conserved variables (Figure 5.1) and
an upwind and unified formulation of fluxes and source terms (Murillo & Garcı́a-Navarro 2010b)
Z NE
∂ X
UdΩ + (En − S̄)k lk = 0 (5.7)
∂t Ωi k=1
5.4 Finite-Volume Numerical Solution 25
The approximate solution can be defined using an approximate Jacobian matrix J en,k (Roe 1981)
of the non-linear normal flux En and two approximate matrices P e = (e1
e ,e2
e ,e3
e ), and Pe −1 , built
using the eigenvectors of the Jacobian, that make J
en,k diagonal
e −1 J
P en,k P
ek = Λ
ek (5.8)
k
Both the difference in vector U across the grid edge and the the source term are projected onto
the matrix eigenvectors basis
δUk = P
e k Ak (S̄)k = P
e kB (5.10)
1 2
where Ak = (α , α , α3 )Tk
contains the set of wave strengths and B = (β , β contains the 1 2
, β 3 )Tk
source strengths. Details are given in (Murillo & Garcı́a-Navarro 2010b). The complete lineariza-
tion of all terms in combination with the upwind technique allows to define the numerical flux
function (En − S̄)k as
3
X m
(En − S̄)k = Ei nk + e− θαe
λ e (5.11)
k
m=1
m
e− = 1 e e and θm = β
with λ 2 (λ − |λ|) k 1− that when inserted in (5.71) gives an explicit first
λα
e
k
order Godunov method (Godunov 1959)
NE
" 3
#
X X m lk
Un+1 = Uni − Ei nk + e−
λ θαe
e ∆t (5.12)
i
m=1
k Ai
k=1
As the quantity Ei is uniform per cell i and the following geometrical property is given at any
cell
NE
X
nk lk = 0 (5.13)
k=1
The finite-volume method can be written using a compact wave splitting formulation as follows:
NE n l
k
X
Un+1
i = Uni − δM−
i,k ∆t (5.15)
Ai
k=1
with
3
X m
δM−
i,k = e− θαe
λ e (5.16)
k
m=1
27 5 RiverFlow2D Model Equations
The use of (5.15) is efficient when dealing with boundary conditions. At the same time it ensures
conservation. In (Murillo & Garcı́a-Navarro 2010b) it was demonstrated how for a numerical scheme
written in splitting form, the total amount of contributions computed inside the domain at each
cell edge, is equal to the balance of fluxes that cross the boundary of the domain, proving exact
conservation.
min(χi , χj )
∆t ≤ CF L ∆tλ ∆tλ = (5.18)
e e
em |
max |λ
with CF L=1/2, as the construction of finite-volume schemes from direct application of one-
dimensional fluxes leads to reduced stability ranges (Toro 2001).
RiverFlow2D solution method uses variable time steps. The maximum allowed time-step is
controlled by the user-set Courant-Friederich-Lewy (CFL) number that is proportional to the
local element size, but also inversely proportional to velocity and depth. Smaller elements lead to
smaller time-steps. The maximum theoretical CFL value is 1, but in some runs it may be necessary
to reduce this number to lower values.
(see Figure 5.4). There is no restriction on the number of inlet or outlet boundaries. This section
describes the open boundary conditions.
If for some reason you need to impose open conditions on boundary segments that are adjacent,
do it in such a way that each segment is separated by a gap of one or more cells (see Figure 5.5).
Setting two or more open conditions without this separation will lead to incorrect detection of the
open boundaries.
RiverFlow2D allows having any number of inflow and outflow boundaries with various combina-
tions of imposed conditions. Proper use of these conditions is a critical component of a successful
RiverFlow2D simulation. Shallow water equation theory indicates that for two-dimensional sub-
critical flow it is required to provide at least one condition at inflow boundaries and one for outflow
boundaries. For supercritical flow all conditions must be imposed on the inflow boundaries and no
boundary condition should be imposed at outflow boundaries. The table below helps determining
which conditions to use for most applications.
The open boundary condition options are described in the table below.
5.5 Open Boundary Conditions 30
In most small slope rivers, the stage-discharge relationship is affected by hysteresis. In other words,
the stage-discharge curve is looped with higher discharges occurring on the rising limb than on the
rescission limb of the hydrograph. This is mainly caused by the depth gradient in the flow direction
that changes in sign throughout the hydrograph. In practice, this implies that there can be two
possible stages for the same discharge. Loop stage-discharge relationships are not considered in
this RiverFlow2D version.
In the cells forming the inlet discharge region the flow is characterized by the negative sign of
the following scalar product in the kΓ boundary edges
and by the state of the flow, defined commonly through the Froude number
ui · ni,kΓ
F ri = (5.20)
ci
√
with ci = ghi . When the Froude number defined as in (5.20) is greater than one, the flow is
supercritical and all the following eigenvalues are negative:
therefore the values of h, u, v, and φ must be imposed. The water solute concentration φ is inde-
pendent of the eigenvalues, and therefore has to be provided at the inlet region for all flow regimes.
w
F rs = p (5.24)
g(ST /bT )
being the cross sectional velocity w = Q/ST and defining the total wet cross section ST and total
breath as:
NB
X NB
X NB
X
ST = Sj = hj lj , bT = lj (5.25)
j=1 j=1 j=1
where N B is the number of wet boundary cells, lj is the length of each edge conforming the wet
boundary and hj is the water depth at each boundary cell.
33 5 RiverFlow2D Model Equations
The total inlet discharge at time t, QI (t), can be distributed along the inlet cross section using
a constant discharge per unit width, qI (m2 s−1 ), that can be calculated as
QI (t)
qj = qI = (5.26)
bT
In this simple case, qI is uniform along the inlet boundary and so is the resulting modulus of
the velocity, w = qI /h, with w = (u2 +v 2 )1/2 . It should be noted, that the direction of the entering
discharge is not necessarily the same as the direction normal to the inlet boundary. However, this
direction is usually chosen as the default information.
When dealing with inlet sections like that in Figure 5.8, a uniform value of qI as in (5.26) leads
to a completely unrealistic state of faster water at the section borders and slower water at the
middle of the cross section. Since the resulting velocities depend on the value of water depth h,
higher values will appear in those cells where water depth is smaller.
In order to seek a more appropriate distribution, a uniform modulus of the water velocity w
is enforced in the whole inlet boundary cross section. In this case, the unit discharge at each
boundary cell j is variable and defined depending both on the total cross section area, ST , and on
the individual cell transverse area, Sj as follows:
Sj
qj = QI (5.27)
ST lj
On the other hand, the updating of the water depth values at the inlet cells provided by
the numerical scheme leads in the general case to a set of new water depths hn+1
j (Figure 5.9)
associated, in general, to different water surface levels dj dj = hj + zj .
For our purposes a horizontal water surface level is required in that region, in order to help
in the translation between the 2D and the 1D points of view at the open boundary. The value
of that uniform cross sectional water level is fixed taking into account mass conservation, that is,
conservative redistribution of water volume. The minimum value of the water levels among all the
wet cells in the inlet boundary, dmin , is found and the water volume VS stored in the inlet section
above dmin is evaluated as
NB
X
VS = (dj − dmin )Aj |dj >dmin (5.28)
j=1
35 5 RiverFlow2D Model Equations
They are used to redistribute the volume over the inlet section, keeping constant the wet section
breadth bT . As Figure 3 shows a new uniform water level at the section, dS , is given by:
VS
ds = dmin + (5.30)
Aw
Apart from helping to decide the flow regime at the boundary, the modifications described
above make easier the treatment of supercritical inflow conditions. When modeling unsteady river
flow, high peaks in the hydrograph can be encountered. If those peaks are not correctly handled
from the numerical point of view, they can lead to local and unrealistic supercritical states in the
inlet boundary.
In that case of supercritical inlet flow, the specification of all the variables at the inlet boundary
cells is required. However, in many practical problems only the discharge hydrograph is available
as a function of time, with no data, in general, on the water level distribution or discharge direction
at the inlet boundary.
The alternative proposed is, when the inlet Froude number is bigger than 1
w
F rs = p >1 (5.31)
g(ST /bT )
to enforce a maximum Froude number, F rs,max , to the inlet flow. For that purpose, keeping
the section breadth bT , a new inlet wet cross section area, ST∗ , is computed from the F rs,max
imposed:
Q2I
ST∗ = ( 2
)1/3 (5.32)
gF rs,max /bT )
If ST∗ is greater than ST , it provides a new water surface level for the inlet section, d∗ , also
greater than ds (Figure 5.10). The associated increment in water volume is balanced by means of
a reduction in the imposed discharge QI (t) in that time step.
Occasionally, both conditions, QI (t) and d(t) are known at supercritical inlets. For those cases,
imposing both data at the inlet boundary is enough. However, due to the discrete time integration
method used, this procedure does not follow the mass conservation criterion. To guarantee that
the mass balance is preserved, one of the conditions is imposed, the other must be modified, so
that the fluxes calculated in the following step lead to mass conservation. The best solution is to
impose directly the global surface water level at the inlet boundary section,d(t), and to adapt the
discrete inlet discharge to ensure that the final volume is conserved. The imposed value of d sets
an input volume that can be transformed into discharge by means of dividing it by the time step.
This value is added to the discharge leading to a correct mass balance.
5.5 Open Boundary Conditions 36
When the boundary cell belongs to an open boundary where the inlet flow discharge is the
condition imposed and the flow is subcritical, the discharge is computed using (5.27) and imposed
in the boundary cell. Moreover, the water level is computed as a results of the contributions from
that other cell edges in (5.15) when updating the conserved values in the boundary cell at time
level n + 1 and is carefully redistributed as explained before.
In that case, the procedure to follow is well described in (Murillo & Garcı́a-Navarro 2010b).
RiverFlow2D drying and wetting algorithm is an adaptation of the the one originally proposed
by (Brufau, Garcı́a-Navarro & Vázquez-Cendón 2004) and later improved by (Murillo, Burguete,
Brufau & Garcı́a-Navarro 2007) and (Murillo & Garcı́a-Navarro 2010b) in the finite-volume context
and works as follows:
1. At the beginning of each time-step all cells are classified as wet or dry according to the
definition.
2. If a cell is dry and completely surrounded by dry cells, it is removed from the computations
and velocity components are set to zero for the ongoing time step.
5.7 Volume Conservation 38
3. All the internal cell edges are classified as active or inactive according to the definition.
4. Wet/dry cell edge contributions are computed assuming the edge is a solid boundary and
the velocities on both sides are set to zero.
5. The rest of the cell edge contributions are computed according to the numerical scheme as
described above.
6. Wet cells and dry cells surrounded by at least one wet cell are retained in the computation
and solved with the updating scheme using the contributions from the cell edges.
This method generates stable numerical solutions without spurious velocities over dry areas and
offers machine accuracy mass conservation errors allowing the use of the classical CFL condition.
where nj and nm are the directions of the flow in the inlet and in the outlet cells respectively.
The volume variation in the domain of calculation can be only due to
∆M (∆t) 6= 0 (5.35)
Therefore, the mass error of the numerical solution is measured by comparing the total amount
of water calculated at time t + ∆t
N CELLS
X
V ol(t + ∆t) = hn+1
i Si (5.36)
i=1
as follows
∂U ∂F(U) ∂G(U)
+ + = S(U) + R(U) (5.40)
∂t ∂x ∂y
where
PNp T
U= h, qx , qy , hφ1 , ..., hφNp , p=1 [zp (1 − pp )] (5.41)
are the conserved variables with h representing the water depth, qx = hu and qy = hv the
unit discharges, with (u, v) the depth averaged components of the velocity vector u along the
x and y coordinates respectively and φp , with p = 1, ..., Np representing the scalar depth-averaged
concentration of the Np different sediments transported in suspension. The term z is the bed level
5.9 Sediment Transport Model: RiverFlow2D ST Module 40
PNp
and pp is the bed porosity of each p bed sediment. The sum of p=1 [zp (1 − pp )] is motivated for
the possibility of a heterogeneous soil, where different fractions of material may coexist.
On the other hand the fluxes are given by
2
T
qx
F= qx , h + 12 gh2 qx qy /h, qx φ1 , ... qx φNp , qsx
T (5.42)
qy2
G= qy , qx qy /h, + 12 gh2 q y φ1 , ..., qy φNp , qsy
h
where g is the acceleration of the gravity and qsx and qsy are the bed load transport terms computed
by means of an empirical law. The source terms of the system are split in three kind of terms.
The term S is defined as
T
pb,x τb,x pb,y τb,y
S = 0, − , − , 0, . . . , 0, 0 (5.43)
ρw ρw ρw ρw
with pb,x , pb,y and τb,x , τb,y are the pressure force along the bottom and the bed shear stress in
the x and y direction respectively, with ρw the density of water. The former can be formulated in
terms of the bed slopes of the bottom level z
pb,x ∂z pb,y ∂z
= −gh = ghS0x , = −gh = ghS0y (5.44)
ρw ∂x ρw ∂y
and the friction losses are written in terms of the Manning’s roughness coefficient n
√ √
τb,x n2 u u2 + v 2 τb,y n2 v u2 + v 2
= ghSf x Sf x = , = ghSf y Sf y = (5.45)
ρw h4/3 ρw h4/3
and is evaluated as
where the term φ∗p is employed for defining the equilibrium concentration, which is obtained
through a solid transport discharge law. The term φp contains the information about the suspension
sediment quantity which is transported. Both concentrations take in consideration the settling
velocity, ωsp .
Mass conservation
∂(h) ∂(hu) ∂(hv)
+ + =0 (5.48)
∂t ∂x ∂y
Momentum conservation in x direction
∂(hu) ∂[hu2 + (1/2)gh2 ] ∂(huv) pbx τbx
+ + = − (5.49)
∂t ∂x ∂y ρw ρw
Momentum conservation in y direction
The bed load transport is often represented by the following dimensionless parameter,
|qs |
Φ= p (5.53)
g(s − 1)d350
where s is the solid material vs water density ratio and d50 is the median diameter.
Table 5.3 collects the formulas that are implemented in RiverFlow2D, where d90 , d50 and d30
are the grain diameter for which 90%, 50% and 30% of the weight of a non-uniform sample is finer
respectively, ρs is the sediment density, θc is the critical Shields parameter and θcs is the critical
Shield parameter as expressed by Smart (1984).
PNp
The sum of p=1 [zp (1 − pp )] in (5.51) can consider heterogeneous soil, where different fractions
of material may coexist. In order to take into account this heterogeneity, every sediment transport
discharge (qsx , qsy ) associated to each sediment size is multiplied by its corresponding soil fraction.
For every sediment particle size it is necessary to include the fraction of the material.
3
0.0 1000. 0.16 0.23
2.0 7000. 1.16 2.24
100.0 7000. 1.16 2.24
In file there are 3 times: 0, 2 and 100 hours. The first column represents the time in hours.
The second column indicates the water discharge in m3 /s or f t3 /s. The last two columns have
the sediment discharge for each given fraction in in m3 /s or f t3 /s.
Mass conservation
∂(h) ∂(hu) ∂(hv)
+ + =0 (5.54)
∂t ∂x ∂y
Momentum conservation in x direction
∂(hu) ∂[hu2 + (1/2)gh2 ] ∂(huv) pbx τbx
+ + = − (5.55)
∂t ∂x ∂y ρw ρw
settling velocity, ωsp . The settling velocity can be computed through several laws as a function of
the sediment density and the grain diameter.
Table 5.4 collects the formulas considered in RiverFlow2D. It is worth noticing that these laws
have been derived for one single particle and under steady situations. Hence, in case of considering
a transient situation with a sediment concentration in the water column, the influence of the
nearby particles must be taken into consideration, which can strongly reduce the fall velocity in
comparison with clear water. This effect, known as hindered settling (Wu 2008, Baldock, Tomkins,
Nielsen & Hughes 2004), can be determined by (Sha 1965), computing a new settling velocity ωsm
n
φ
ωsm = 1 − √ ωs (5.58)
2 d50
where ωs is the settling velocity of a single particle in clear water, φ is the volumetric concentration,
d50 is the mean diameter and n is an empirical exponent with a usual value equal to 3.
• Bagnold (Bagnold 1973) considered that the shear stress is proportional to the flow velocity
and it was established the following formula,
ρs − ρw τb |u|2
qsuspp = 0.01 (5.60)
ρs g ωs
being τb the shear stress generated at the bottom by the bed roughness which is taken into
account through the Manning’s coefficient, n
gn2 |u|2
τb = (5.61)
h1/3
• Van Rijn (Van Rijn 1984c) proposed calculating the suspended load integrating the sediment
flux within the layer where the suspension plays a key role, i.e. between the term hs and h,
Z h
qsuspp = φm udh (5.62)
hs
where, qsusp is expressed by volume per unit time, and the terms of φm and hs were defined
as follows
h
!ωs /(kU∗ )
h0
−1
φm = φs h
(5.63)
hs −1
with
1/3
ν2
0.117 θ
φs = −1 (5.64)
ds (s − 1)g θc
where θ is the non-dimensional shear stress, calculated using (5.66) and (5.67) and related
to n, the Manning roughness coefficient, θc is the Shields’s parameter, ds is the sediment
diameter, s is the rate between densities, s = ρρws and ν is the kinematic viscosity.
τb
θ= (5.66)
ρw (s − 1)gds
with
gn2 |u|2
τb = (5.67)
h1/3
√
In equation 5.63, U∗ is the friction velocity defined as U∗ = τb and κ is the Von-Karman
constant κ ≈ 0.41.
45 5 RiverFlow2D Model Equations
3
0 1000. 0.001 0.002
2 67000. 0.001 0.002
100 67000. 0.001 0.002
Z NE Z ek+1 Z Z
∂
(En)↓k dlk =
X
UdΩ + SdΩ + RdΩ (5.69)
∂t Ωi ek Ωi Ωi
k=1
with (En)↓k the value of the interface flux function through the edge k to be defined, nk = (nx , ny )
is the outward unit normal vector to the cell edge k, and N E is the number of edges in cell i. A
sketch of the fluxes is showed in Figure 5.12.
Assuming a piecewise representation of the conserved variables
Z
1
Uni = U(x, y, tn )dΩ (5.70)
Ai Ωi
5.9 Sediment Transport Model: RiverFlow2D ST Module 46
where lk is the corresponding edge length. System (5.71) is solved following the theory of Roe’s Rie-
mann solver and using the upwind discretization [(Murillo & Garcı́a-Navarro 2010a, Juez, Murillo
& P. 2014)].
• When using the bed load option, the inflow sediment transport rate is equal to the transport
capacity under equilibrium conditions.
• Each Sediment particle size is considered using a single characteristic diameter (D50 ) for all
formulas except for that of Smart (1984) that considers three sizes (D30 , D50 and D90 ).
• When assigning maximum erosion bed elevation, the model blocks further erosion when the
bed elevation reaches the given limit.
Since the formulations for the bed load discharge, the suspended load discharge and the settling
velocity are based on deterministic laws supported by experimentation, tuning parameter factors
have been considered for each one. Hence, the model makes possible the calibration of the computed
sediment transport for each particular situation.
• The flow is confined to a layer which is thin compared to the horizontal scale of interest;
• The flow is governed by equation (5.1) with the friction terms evaluated as explained below;
The mixture of water and sediments is described by using the continuum approach and
assuming the same velocity for the liquid and for the solid phase;
• The fluid is assumed to be an homogeneous single-phase mix of water and sediment and has
constant properties: e.g. density, yield stress, etc.
generated by the stresses can be lumped into the same mathematical formula despite of having a
different nature. In this way, several shear stresses can be considered to cover a wide range of mud
and debris flow conditions, including:
• Turbulent stress τt
• Dispersive stress τd
• Viscous stress τµ
The equations that describe the tangential forces generated by the stresses involve turbulent
stress τt , dispersive stress τd , Coulomb-type frictional stress τf , yield stress τy , and viscous stress
τµ . Not all stresses act simultaneously in the same location of the water column. However, since
the RiverFlow2D conceptual model is depth-averaged, all terms may actually coexist and are
mathematically lumped in the same formula. In the case of a viscous Newtonian fluid, the bed
stress τb is expressed in function of the depth averaged velocity u and depth h as
A Bingham fluid will not flow until a certain threshold value of the stress, the yield stress τy ,
is reached. Once flowing, the movement is characterized by a plastic viscosity of the mixture. The
bed stress τb can be obtained from the following cubic formula
and in cases of stress ratio τy /τb smaller than 0.5 can be reduced to
3
τb = τy + 3τµ (5.74)
2
49 5 RiverFlow2D Model Equations
If only the Coulomb-type frictional stress is considered, the bed shear stress at the bed is given
by
τb = τt = ρcf u2 (5.76)
In general the total shear stress τb can be written as a combination of stress components:
τb = f (τt , τd , τy , τµ , τf ) (5.77)
Table 5.5 shows the rheological friction laws available in RiverFlow2D MD. The Bingham-type
relations include the standard cubic formulation
5.10.3 Entering Data for the Mud and Debris Flow Model
To enter data for a mud and debris flow simulation use the Mud/Debris Flow panel in the
RiverFlow2D Data Input Program. Also make sure that the Mud/Debris Flow check box is active
in the Control Data panel.
Depending on the selected flow resistance relation, the user will need to enter different pa-
rameters. For example the Full Bingham formulation requires yield stress, Bingham viscosity and
Material density, etc..
transport equations will be expressed in a conservative form, assuming that the velocities and the
water depth may not vary smoothly in space and time.
Correspondingly, the 2D shallow water model with solute transport can be written in unique
coupled system:
∂U ∂F(U) ∂G(U)
+ + = H(U) + R(U) + D(U) (5.79)
∂t ∂x ∂y
where
T
U = (h, qx , qy , hφ)
T T
qx qy qx qy qy 2
F = qx , qxh2 + 12 gh2, h , hφu , G = qy , h , h + 21 gh2, hφv (5.80)
T
H = (0, gh(S0x − Sf x ), gh(S0y − Sf y ), 0)
and φ is the depth-averaged solute concentration. The sources terms associated to the solute
transport equation are expressed as follows:
T
→
− →
− T
R = (0, 0, 0, −Khφ) D = 0, 0, 0, ∇(Dh ∇φ) (5.81)
3 m
qk↓
X
= qi + e− γ
λ eee1 (5.82)
k
m=1
5.11 Pollutant Transport Model: RiverFlow2D PL Module 52
where qi = (hun)i and the decoupled numerical scheme for the solute transport equation is
written as:
N
E
∆t X
(hφ)n+1
i = (hφ)ni − (qφ)↓k lk (5.83)
Ai
k=1
where
qk↓ > 0
φi if
φ↓k = (5.84)
φj if qk↓ < 0
in cell i. A sketch of the fluxes is showed in Figure 5.12.
From a physical point of view, the new solute mass at a fixed cell can be seen as exchanging
water volumes with certain concentration through the neighboring walls and mixing them (finite-
volume Godunov’s type method) with the former mass existing in the previous time (Figure 5.16).
According to this philosophy, the outlet boundary cells will require a special treatment when
applying this technique in order to extract the corresponding solute mass through the boundary
walls. For this reason, it is necessary to define q ↓ = (hu · n)BC and φ↓ = φBC at the boundary
wall and to include this contribution for the updating of the boundary cell BC (see Figure 5.17).
↓ ↓
Figure 5.16 – Physical representation of solute mass exchange between cells with qi−1/2 , qi+1/2 > 0.
53 5 RiverFlow2D Model Equations
As shown, the formulation reduces to compute a class of numerical flux q ↓ using the already
computed averaged values at each edge for the hydrodynamic paPollutant Transport tabrt. Apart
from ensuring a perfect conservation and bounded free-oscillatory solutions (Murillo et al, 2012),
this simply discretization decreases substantially the number of computations that would be nec-
essary for the complete coupled system.
55
6.2 RiverFlow2D GPU 56
Figure 6.1 – Speed up using RiverFlow2D parallelized code as a function of number of processor
cores.
Figure 6.2 – Speed up using RiverFlow2D parallelized code as a function of number of processor
cores.
• Rainfall and Evaporation: accounts for spatially distributed rainfall and evaporation.
• Wind: allows incorporating the effect of spatially distributed wind stress on the water
surface.
The following hydraulic components are set over polylines (feature arcs):
• Bridges: account for general geometry bridges including pressure flow and overtopping.
• Internal Rating Tables: provide an internal relationship of water elevation and discharge..
• Weirs: represent crested structures such as weirs, levees, sound walls, etc., where there is a
unique relationship between discharge and depth.
• Culverts: one dimensional conveyance conduits where discharge can be calculated using
equations for circular or box structures, and rating tables.
• Sources and Sinks: provide a mean to enter point inflows or outflows that may vary in
time.
57
7.1 Bridges Component 58
Figure 7.1 – RiverFlow2D Control Data panel with Bridges component selected.
The bridge plan data is entered in the SMS conceptual model coverage. To create a bridge,
first make sure that the coverage where you defined the domain outline is active. Then, create an
arc that will define the bridge geometry in plan.
Now, draw a polyline clicking on each vertex until you reach the desired end when you can
double click to complete the polyline.
59 7 Hydraulic and Hydrologic Components
Then select the polyline using the Select Feature Arch icon, and right click to redistribute the
vertices on the polyline. While the polyline is selected, right click on the arc and select Attributes...
and on the Feature arc Attributes dialog, select Boundary conditions and click Options...
On the Boundary Conditions Dialog, select the Bridges component and enter the desired
BridgeID (any text), and the bridge geometry file name that contains the bridge geometry data
and is explained in the next section.
Since the bridge polyline must pass through nodes, it is essential that mesh generation engine
creates nodes along the bridge polyline. This is easily done recreating the mesh using the Feature
Objects/Map → 2D Mesh command.
7.1 Bridges Component 60
NP
Table 7.1 – Variable Descriptions for the bridge cross section geometry file.
VARIABLE TYPE RANGE UNITS DESCRIPTION
Bed elevation. Must be the lowest ele-
BEDELEV R - m or ft
vation for all polylines at a given point.
Elevation of the bridge deck. Must be
DECKELEV R - m or ft the highest elevation for all polylines at
a given point.
Elevation of the lower bridge deck.
LOWCHORD must be larger or equal
to ZLOWER and smaller or equal to
LOWCHORD R - m or ft DECKELEV for a particular point.
The space between LOWCHORD and
DECELEV is a blocked area to the
flow.
Number of points defining cross section
NP I - >1
polylines.
Distance from leftmost point defining
STATION R - m or ft cross section polyline. All polylines
points must have a common station.
Elevation of lower polyline. ZLOWER
must be larger or equal to BEDELEV
and smaller or equal to LOWCHORD
for a given point. The space between
BEDELEV and ZLOWER is a blocked
ZLOWER R - m or ft
area to the flow. The space between
ZLOWER and LOWCHORD is open
space. If the bridge has no holes,
ZLOWER must be identical to BE-
DELEV.
Note: R = Real variable. I = Integer variable. S = Text variable.
Figure 7.7 – Top view of a bridge showing the cross sections of interest. Only two piers are depicted
for simplicity.
v̄42
2 (7.5)
A4
∆He = 2g A3 − 1 + 91
where m is a typical value for the contraction coefficient, m = 0.62 (Ratia et al. 2014) and the
areas A1 to A4 refer to effective cross sectional flow area. The numbering of areas is shown in
Figures 7.7. Area 1 is a section upstream of the bridge while area 4 is a downstream section. Areas
2 and 3 are sections inside the bridge, near the entrance and exit respectively.
Figure 7.8 – Simple example of A1 ,A2 , A3 and A4 used to calculate head loss in free surface bridges.
Figure 7.9 – Simple example of A1 ,A2 , A3 and A4 used to calculate head loss in a partially submerged
bridges.
The values v̄1 and v̄4 are the cross sectional averaged velocities
Q1 Q4
v̄1 = , v̄4 = (7.6)
A1 (d1 ) A4 (d4 )
with Q1 and Q4 the total discharges in areas A1 and A4 , expressed as a function of the water
surface elevation, d = h + z. Different regimes can be described. Figure 7.8 shows a skecth of
7.1 Bridges Component 64
Figure 7.10 – Simple example of A1 ,A2 , A3 and A4 used to calculate head loss in fully submerged
bridges.
the areas considered in the free surface case, Figure 7.9 shows the equivalent areas for partially
submerged bridges and Figure 7.10 for fully submerged bridges.
with
u
en
δH = ∆H (7.8)
|e
u · n|
where ∆H is the singular loss term used to represent bridges.
Computation of ∆H in a real mesh is done as follows. The bridge is defined on cell edges (bold
line in Figure 7.11), and the cells on both sides of these edges are considered to form two cross
sections ΓL and ΓR (hatched cells in Figure 7.11). Note that it is possible to define bridges in
arbitrary orientations and in structured/unstructured meshes.
Figure 7.11 – Application of the scheme in triangular structured meshes. Normal bridge (left) and
oblique bridge (right).
In each time step, the necessary variables for the calculation of the global bridge head loss are
averaged from the cells in both upstream and downstream sections as illustrated in Figure 7.11.
The discharge is computed as
X X
QΓL = (qn)k lk QΓR = (qn)k lk (7.9)
k∈ΓL k∈ΓR
65 7 Hydraulic and Hydrologic Components
involving cells with values of h > 0. The signs of QΓL and QΓR are used to determine which section
is upstream and which downstream. If QΓL ≥ 0, the discharge across the bridge is computed as
Q = QΓL and the areas are computed using d1 = dΓL and d4 = dΓR . In case that QΓL < 0,
the discharge across the bridge is computed as Q = QΓR and the sections are reversed setting
d1 = dΓR and d4 = dΓL . Next, the different areas and the cross-sectional top width are calculated
as a function of the average water level surface. From these values the total head loss ∆H can be
evaluated.
Figure 7.12 – Bridge pier proportions used to asses the influence of the structure width.
Numerical results indicate that the changes in total head loss across the structure are very
similar for the three configurations (see Figure 7.13).
Figure 7.13 – Influence of the structure width on the total head change (∆H) across the bridge as
a function of the Froud number downstream.
7.2 Bridge Piers 66
Figure 7.14 – RiverFlow2D Control Panel dialog with Bridge Piers component selected.
Bridge pier data is entered in the SMS conceptual model coverage. To create a bridge pier,
first make sure that the Pier Component is active and then follow these steps:
• Click on the Create Feature Point icon and click on the location where you want to enter
the pier
• Now click on the Select Feature Point icon and click on the point entered in the previous
step
• Select Pier, enter the Pier Name and other data as shown:
67 7 Hydraulic and Hydrologic Components
Note: There is no limit to the number of Bridge Piers that can be used.
1
FD = CD ρU 2 AP (7.11)
2
Where CD is the pier drag coefficient, ρ is the water density, U is the water velocity, and AP
is the pier wetted area projected normal to the flow direction. Piers are assumed to be located on
elements that not necessarily conform to the pier geometry as shown on the following figure.
Piers can be circular or rectangular in plan. Rectangular piers are located on elements based
on the pier center coordinates and the angle between the axis along the largest dimension and the
X-axis as shown in the following figure.
Velocity vector magnitude and approach angle usually varies in time during unsteady flow
computations and is used to calculate the projected area. To account for the resistance force that
7.3 Culverts Component 68
the pier exerts on the flow, RiverFlow2D converts it to the distributed shear stress on the element
where the pier centroid coordinate is located. The resulting pier shear stress expressions in x and
y directions are as follows:
1 p AP
τpx = CD ρU U 2 + V 2 (7.12)
2 Ae
1 p AP
τpy = CD ρV U 2 + V 2 (7.13)
2 Ae
where Ae is the element area.
Figure 7.18 – RiverFlow2D Global Parameters dialog with Culverts Component selected.
1. Click on the Create Feature Point icon and click on the location where you want to enter
the inlet point.
2. Now click on the Select Feature Point icon and right click on the point entered in the
previous step.
5. Select Culvert, enter the Culvert ID name (e.g. Culvert1) and enter additional data depend-
ing on the culvert formulation desired.
6. Then select the Culvert Outlet Element and enter the Culvert Inlet Name:
7. Before running the model, make sure that Culverts it is selected in the Global Parameters
dialog.
There are two options to compute culvert discharge in RiverFlow2D. When the user selects
Rating Table calculation and provides a rating table on the associated file, the model determines
the discharge by interpolation as a function of the depth upstream. If the user enters Culvert
calculation using culvert characteristics, the model will calculate the discharge based on the cul-
vert geometric characteristics given in the file. Both procedures are described in more detail below.
1. If at least one of the culvert ends is wet, determine the flow direction based on the water
surface elevations at each culvert end,
2. Interpolate flow discharge from the rating table using the depth at the culvert inlet,
3. If depth at the culvert inlet is lower than minimum value in the rating table, then the
discharge is assumed to be zero.
4. If depth at entrance is higher than maximum value in the rating table, then the discharge is
assumed to be equal to that of the maximum depth.
5. The computed discharge is subtracted from the inlet element and added to the outlet element
assuming instantaneous water volume transmission.
where Hh = W SELh − Zbi is the headwater depth. Dc is the culvert diameter for circular culverts
and the base dimension for box culverts, m = 0.7 for mitered inlets and m = −0.5 for all other
inlets. For outlet control, the following formula is used to determine Cc :
−0.5
2gn2c Lc
Cc = 1 + Ke + 4/3
(7.16)
Rc
where Rc is the culvert hydraulic radius, Ke is the entrance loss coefficient that can be obtained
from Table 7.3, nc is the Manning’s n obtained from Table 7.2, Lc is the culvert length, and Y ,
K 0 , M 0 , C 0 are inlet control coefficients (see Table 7.4).
71 7 Hydraulic and Hydrologic Components
Table 7.2 – Manning’s n roughness coefficients for various culvert materials. Adapted from Froehlich
(2003).
CULVERT BAR- ENTRANCE DESCRIPTION MANNING’S N Nc
REL MATERIAL
Concrete Good joints, smooth walls 0.012
Projecting from fill, square-cut end 0.015
Poor joints, rough walls 0.017
Corrugated metal 2-2/3 inch × 1/2 inch corrugations 0.025
6 inch × 1 inch corrugations 0.024
5 inch × 1 inch corrugations 0.026
3 inch × 1 inch corrugations 0.028
6 inch × 2 inch corrugations 0.034
9 inch × 2-1/2 inch corrugations 0.035
Table 7.4 – Culvert inlet control formula coefficients. Adapted from Froehlich (2003).
BARREL BARREL INLET DESCRIP- K’ M C’ Y
MATE- SHAPE TION*
RIAL
Concrete Circular Headwall; square edge 0.3153 2.0000 1.2804 0.6700
Concrete Circular Headwall; grooved 0.2509 2.0000 0.9394 0.7400
edge
Concrete Circular Projecting; grooved 0.1448 2.0000 1.0198 0.6900
edge
Cor. metal Circular Headwall 0.2509 2.0000 1.2192 0.6900
Cor. metal Circular Mitered to slope 0.2112 1.3300 1.4895 0.7500
Cor. metal Circular Projecting 0.4593 1.5000 1.7790 0.5400
Concrete Circular Beveled ring; 45◦ bevels 0.1379 2.5000 0.9651 0.7400
Concrete Circular Beveled ring; 0.1379 2.5000 0.7817 0.8300
33.7◦ bevels
Concrete Rectangular Wingwalls; 30◦ to 0.1475 1.0000 1.2385 0.8100
75◦ flares; square edge
Concrete Rectangular Wingwalls; 90◦ and 0.2242 0.7500 1.2868 0.8000
15◦ flares; square edge
Concrete Rectangular Wingwalls; 0◦ flares; 0.2242 0.7500 1.3608 0.8200
square edge
Concrete Rectangular Wingwalls; 45◦ flare; 1.6230 0.6670 0.9941 0.8000
beveled edge
Concrete Rectangular Wingwalls; 18◦ to 1.5466 0.6670 0.8010 0.8300
33.7◦ flare; beveled
edge
Concrete Rectangular Headwall; 3/4 inch 1.6389 0.6670 1.2064 0.7900
chamfers
Concrete Rectangular Headwall; 45◦ bevels 1.5752 0.6670 1.0101 0.8200
Concrete Rectangular Headwall; 33.7◦ bevels 1.5466 0.6670 0.8107 0.8650
Concrete Rectangular Headwall; 45◦ skew; 1.6611 0.6670 1.2932 0.7300
3/4 in chamfers
Concrete Rectangular Headwall; 30◦ skew; 1.6961 0.6670 1.3672 0.7050
3/4 in chamfers
Concrete Rectangular Headwall; 15◦ skew; .7343 0.6670 1.4493 0.6800
3/4 in chamfers
Concrete Rectangular Headwall; 10-45◦ skew; 1.5848 0.6670 1.0520 0.7500
45◦ bevels
Concrete Rectangular Wingwalls; non-offset 1.5816 0.6670 1.0906 0.8030
45◦ /flares;
Concrete Rectangular Wingwalls; non-offset 1.5689 0.6670 1.1613 0.8060
18.4◦ /flares; 3/4 in
chamfers
Concrete Rectangular Wingwalls; non- 1.5752 0.6670 1.2418 0.7100
offset 18.4◦ /flares;
30◦ /skewed barrel
Concrete Rectangular Wingwalls; offset 1.5816 0.6670 0.9715 0.8350
45◦ /flares; beveled top
edge
Concrete Rectangular Wingwalls; offset 1.5752 0.6670 0.8107 0.8810
33.7◦ /flares; beveled
top edge
Concrete Rectangular Wingwalls; offset 1.5689 0.6670 0.7303 0.8870
18.4◦ /flares; top edge
bevel
Cor. metal Rectangular Headwall 0.2670 2.0000 1.2192 0.6900
Cor. metal Rectangular Projecting; thick wall 0.3023 1.7500 1.3479 0.6400
Cor. metal Rectangular Projecting; thin wall 0.4593 1.5000 1.5956 0.5700
Concrete Circular Tapered throat 1.3991 0.5550 0.6305 0.8900
Cor. metal Circular Tapered throat 1.5760 0.6400 0.9297 0.9000
Concrete Rectangular Tapered throat 1.5116 0.6670 0.5758 0.9700
Concrete Circular Headwall; square edge 0.3153 2.0000 1.2804 0.6700
Continued on next page
73 7 Hydraulic and Hydrologic Components
1. If at least one of the culvert ends is wet, Determine the flow direction based on the water
surface elevations at each culvert end.
4. Select the minimum discharge from the inlet and outlet control discharges.
5. If depth at the culvert inlet is lower than minimum value in the rating table, then the
discharge is assumed to be zero.
6. The computed discharge is subtracted from the inlet element and added to the outlet element
assuming instantaneous water volume transmission.
When using CulvertType 1 or 2, both ends of the culvert must be inside the mesh.
75 7 Hydraulic and Hydrologic Components
To run a simulation with the gates component, you need to select the option in the Control
Data panel of the RiverFlow2D Data Input Program as shown in Figure 7.21.
Figure 7.21 – RiverFlow2D Control Data panel with Gates Component selected.
The gate plan data is entered in the SMS conceptual model coverage. First create an arc that
will define the gate geometry in plan. Draw a polyline clicking on each vertex until you reach the
desired end when you can double click to complete the polyline.
77 7 Hydraulic and Hydrologic Components
Then select the polyline using the Select Feature Arch icon, and right click to redistribute the
vertices on the polyline.
While the polyline is selected, right-click on the arc and select Attributes... and on the Feature
arc Attributes dialog, select Boundary conditions and click Options...
On the Nodestring Boundary Conditions dialog, select the Gates component and enter the
desired Gate ID (any text), and gate file name that contains the aperture time series, the crest
elevation, crest height and gate Cd.
Since the gate polyline must pass through nodes, it is essential that the mesh generation en-
gine creates nodes along the polyline. This is easily done recreating the mesh using the Feature
Objects/Map → 2D Mesh command.
1/2
q = Go K2 [dr − max(zl , zr )] (7.18)
Figure 7.25 – Water levels for discharge under a gate in submerged conditions formulated as in (G1).
Figure 7.26 – Water levels for discharge under a gate in submerged conditions formulated as in (G2).
7.5.1 Rainfall
Rainfall is treated as a simple source term. It represents an additional input to the cell water depth
in the previous step of flow calculation. You can set up local rainfall events for several regions of
the watershed. This allows to simulate more realistic cases, in which data from several rain gauges
are available.
Notice that the rainfall is set before the infiltration calculation. This is an important detail because
the infiltration capacity strongly depends on the rainfall intensity, as we will see in the next section.
7.5.2 Infiltration
Infiltration represents another component of the hydrological budget and it can be defined as the
process by which surface water enters the soil. In the RiverFlow2D model, infiltration is treated as
a loss. This process is mainly governed by two forces: gravity and capillarity action. The model
includes three methods to compute the infiltration losses: Horton, Green & Ampt and SCS-Curve
Number.
Table 7.6 – Horton initial infiltration for different soils. Source: Akan(1993).
SOIL TYPE f0 (mm/hr)
Dry sandy soils with few to no vegetation 127
Dry loam soils with few to no vegetation 76.2
Dry clay soils with few to no vegetation 25.4
Dry sandy soils with dense vegetation 254
Dry loam soils with dense vegetation 152.4
Dry clay soils with dense vegetation 50.8
Moist sandy soils with few to no vegetation 43.18
Moist loam soils with few to no vegetation 25.4
Moist clay soils with few to no vegetation 7.62
Moist sandy soils with dense vegetation 83.82
Moist loam soils with dense vegetation 50.8
Moist clay soils with dense vegetation 17.78
7.5 Rainfall, Evaporation and Infiltration 80
Table 7.7 – Horton final infiltration for different soils. Source: Akan(1993).
SOIL TYPE fc (mm/hr)
Clay loam, silty clay loams 127
Sandy clay loam 1.3 - 3.8
Silt loam, loam 3.8 - 7.6
Sand, loamy sand, sandy loams 7.6 - 11.4
The equation (7.19) has to be applied after the surface ponding. In other words, we are assuming
conditions of unlimited water supply at the surface. Under this consideration, the cumulative
infiltration up to time t can be calculated by integrating the infiltration capacity:
Z t
f0 − fc
1 − e−kt
F = fp (t)dt = fc t + (7.20)
0 k
It is important to highlight the difference between the infiltration capacity fp and the infiltration
rate f. If we consider a rain event starting with a weak rainfall intensity (R ≤ fp ), then all the rain
will be infiltrated into the soil. On the other hand, if the rain exceeds the soil infiltration capacity
or if the surface becomes ponded, this magnitude will determine the infiltration rate:
fp ≤ R (7.25)
81 7 Hydraulic and Hydrologic Components
Table 7.8 – Mean values and standard deviation for Green-Ampt model parameters. Source: Rawls
et al. (1983).
SOIL TYPE θs Sf (mm/hr) Ks (mm/hr)
Sand 0.437(0.374-0.500) 4.95(0.97-25.36) the78
Loamy sand 0.437(0.363-0.506) 6.13(1.35-27.94) 2.99
Sandy loam 0.453(0.351-0.565) 11.01(2.67-45.47) 1.09
Loam 0.463(0.375-0.551) 8.89(1.33-59.38) 0.66
Silt loam 0.501(0.420-0.582) 16.68(2.92-95.39) 0.34
Sandy clay loam 0.398(0.332-0.464) 21.85(4.42-108.0) 0.15
Clay loam 0.464(0.409-0.519) 20.88(4.79-91.10) 0.10
Silty clay loam 0.471(0.418-0.524) 27.30(5.67-131.50) 0.10
Sandy clay 0.430(0.370-0.490) 23.90(4.08-140.2) 0.06
Silty clay 0.479(0.425-0.533) 29.22(6.13-139.4) 0.05
The original Green-Ampt model starts from the assumption that a ponding depth h is main-
tained over the surface. The Green-Ampt method approximates the soil infiltration capacity as
follows:
Ks (θs − θi ) Sf
fp = Ks + (7.26)
F
being Ks the effective hydraulic conductivity, Sf the suction head at the wetting front, θi the
initial uniform water content and θs the porosity. The integration of (7.26) provides the cumulative
infiltration:
dF F
fp = =⇒ Ks t = F − (θs − θi ) Sf ln 1 + (7.27)
dt (θs − θi ) Sf
Solving for the cumulative infiltration F in equation (7.27) requires an iteration procedure (e.g.
Picard iterations or Newton-Rhapson method). The effective suction head can be replaced by the
average value Ψ (Mein & Larson 1973).
Equations (7.26) and (7.27) assume that the soil is ponded from the beginning. Additional
considerations should be taken into account in order to model an unsteady storm pattern (Te Chow,
Maidment & Mays 1988). Three posibilities can occur in every timestep: 1) ponding occurs at the
beginning of the interval; 2) there is no ponding within the interval; 3) ponding occurs within the
interval. The first step consists of computing the actual infiltration capacity fp from the known
value of the cumulative infiltration F at time t. From (7.26):
Ψ∆θ
fp = K s +1 (7.28)
F
The result from eq. (7.28) is compared with the rainfall intensity i. If fp ≤ i, case 1 occurs
and the cumulative infiltration at the end of the interval is given by (7.29). Moreover, the real
infiltration f rate will be equal to the potential one fp ≤ i:
Ft+∆t + Ψ∆
Ft+∆t − F − Ψ∆θln = K∆θ (7.29)
F + Ψ∆
7.5 Rainfall, Evaporation and Infiltration 82
If fp > i, there is no ponding at the beginning of the interval. We assume that there is no
ponding during the entire interval, so the real infiltration rate is equal to the rain rate and a
tentative value for the cumulative infiltration at the end of the period can be computed as:
0
Ft+∆t = F + i∆t. (7.30)
0
From equations (7.28) and (7.30) a tentative infiltration capacity fp,t+∆t can be calculated. If
0
fp,t+∆t > i, there is no ponding during the interval, the assumption is correct and the problem
0 0
corresponds to situation number 2, so Ft+∆t = Ft+∆t . If fp,t+∆t ≤ i, there are ponding condition
within the interval (case 3). The cumulative infiltration at ponding time Fp is found by taking
fp = i and F = Fp at (7.28):
Ks Ψ∆θ
Fp = (7.31)
i − Ks
Then, the ponding time is computed as t + ∆t0 , where:
Fp − F
∆t0 = (7.32)
i
Finally, the cumulative infiltration can be found by replacing F = Fp and ∆t = ∆t − ∆t0 in
equation (7.29).
25400
S= − 254 (7.34)
CN
As seen on Eq. 7.33, runoff cannot begin until initial abstraction has been met. Initial abstrac-
tion is written in terms of its ratio α and S :
Ia = αS (7.35)
The initial abstraction ratio was traditionally taken as α = 0.2 for every watershed (USDA 1986)
but recent studies suggest that there is a wide range of values that work better than this value,
depending on the soil properties. The influence of this parameter was widely studied in Caviedes
et al. (Caviedes-Voullième, Garcı́a-Navarro & Murillo 2012).
It is important to remark that SCS-CN method was not designed to consider time. Following
(Caviedes-Voullième et al. 2012), when the method is implemented in a complex simulator, a time-
advancing methodology is used. The method is not applied to the entire catchment. Runoff is
83 7 Hydraulic and Hydrologic Components
calculated for every cell in every time step, using the cumulative rainfall since the beginning of the
storm.
Follow the guidelines provided in (USDA 1986) to determine appropriate Curve Numbers.
Another important fact that needs to be considered is the Antecedent Moisture Content (AMC)
or Antecedent Runoff Condition (ARC), which represents the preceding relative moisture of the
soil prior to the storm event (Te Chow et al. 1988). These parameters result useful in order to
attempt to account for the CN variation among storm events. Three posible assumptions can be
considered: dry conditions (AMC I), average conditions (AMC II) or wet conditions (AMC III)
(see Table 7.9).
Traditionally (Te Chow et al. 1988), the Curve Number for dry or wet conditions has been
recalculated in terms of the standard conditions according to Eqs. 7.36 and 7.37:
4.2CN (II)
CN (I) = (7.36)
10 − 0.058CN (II)
23CN (II)
CN (III) = (7.37)
10 − 0.13CN (II)
On the other hand, some newer references (USDA 2004, Mishra & Singh 2003) recommend to
use an empirical data table to compute both values.
In the infiltration without rainfall option only the Horton or the Green and Amp methods can be
used since they take into account the water depth to determine the infiltration rates. The SCS-CN
method calculates infiltration as a function of the given rainfall and does not consider the pounded
water.
Please follow the Rainfall, Evaporation and Infiltration Tutorial that explains how to setup a
hydrologic simulation from start to finish.
7.6 Wind Component 84
where
ρa ρa
Swx = Ca U|U | Swy = Ca V|U | (7.39)
ρw ρw
being U = (U, V) the wind velocity, ρa and ρw the air and water densities respectively, and Cd is
the wind stress coefficient. This source term is discretized using the same upwind technique as the
friction and bed slope source terms.
If the rating table is not fully compatible with the computed 2D flow, results can be erroneous. It
is suggested to use this condition with care to avoid inconsistencies.
To run a simulation with Internal Rating Tables, you need to select the option in the Control
Data panel of the RiverFlow2D Data Input Programas shown in Figure 7.28.
85 7 Hydraulic and Hydrologic Components
Figure 7.28 – RiverFlow2D Control Data panel with IRT component selected.
Internal Rating Table (IRT) plan data is entered in the SMS conceptual model coverage. To
create a IRT polyline, draw a polyline clicking on each vertex until you reach the desired end when
you can double click to complete the polyline.
Then select the polyline using the Select Feature Arch icon, and right click to redistribute the
vertices on the polyline. While the polyline is selected, right click on the arc and select Attributes...
and on the Feature arc Attributes dialog, select Boundary conditions and click Options...
On the Nodestring Boundary Conditions dialog, select the Interior Boundary Condition com-
ponent and enter the desired IRT ID (any text), and IRT file name that contains of water elevation
vs discharge table.
7.7 Internal Rating Tables 86
Enter the desired IRT (any text), and click Define... to enter the rating table. The following
dialog pops up and the Internal Rating Data (Discharge vs. Water Surface Elevation) can be
entered.
Figure 7.32 – Window to enter data pairs (discharge-WSE) for an Internal Rating Table.
Since the IRT must pass through nodes, it is essential that mesh generation engine creates nodes
along the IRT nodestring. This is easily done recreating the mesh using the Feature Objects/Map
→ 2D Mesh command. The resulting mesh nodes line along the IRT polyline as shown in Figure
7.33.
Note: There is no limit to the number of Internal Rating Tables that can be used.
87 7 Hydraulic and Hydrologic Components
1. For each calculation time interval, estimate an average water surface level at each side of the
IRT polyline.
2. Compute the discharge passing through the IRT polyline from the average water levels in 1
using the rating table.
3. Define an average velocity from the discharge and the cross sectional wetted area.
4. Assign a common unit discharge to every pair of cells sharing a polyline segment.
Some inappropriate IRT polyline configurations or very long polylines can over-constrain the model
and should be avoided.
Figure 7.34 – RiverFlow2D Control Data panel with the Sources and Sinks component selected.
Sources and Sinks data is entered in the SMS conceptual model coverage. To create a Source/Sink
follow these steps:
• Then, click on the Create Feature Point icon and click on the location where you want to
enter the source
• Now click on the Select Feature Point icon and click on the point entered in the previous
step
• On the Feature Objects menu, click Attributes...
• Select Source/Sink, enter the Source/Sink ID name (e.g. SOURCE1) and click Define...
89 7 Hydraulic and Hydrologic Components
• Enter the discharge time series (unsteady) or a single discharge (steady flow).
Note: There is no limit to the number of sources and sinks that can be used.
7.9 Weirs
RiverFlow2D Weirs component may be convenient when trying to simulate levee or road overtop-
ping. The tool allows defining a polyline representing the structure alignment and assigning crest
elevations that can vary along the polyline.
To run a simulation with weirs, you need to select the option in the Control Data panel of the
RiverFlow2D Data Input Program as shown in Figure 7.37.
7.9 Weirs 90
Figure 7.37 – RiverFlow2D Control Data panel with Weirs component selected.
Weir plan data is entered in the SMS conceptual model coverage. To create a weir, first make
sure that the conceptual model coverage where you defined the domain outline is active. Then,
create an arc that will define the weir geometry in plan. Now, draw a polyline clicking on each
vertex until you reach the desired end when you can double click to complete the polyline.
Then select the polyline using the Select Feature Arch icon, and right click to redistribute the
vertices on the polyline. While the polyline is selected, right click on the arc and select Attributes...
and on the Feature arc Attributes dialog, select Boundary conditions and click Options...
91 7 Hydraulic and Hydrologic Components
On the Nodestring Boundary Conditions dialog, select the Interior Boundary Condition com-
ponent and enter the desired Weir name, crest elevation and weir coefficient Cf . This version of
the software considers variable crest elevation along the weir.
Since RiverFlow2D requires that the weir passes through nodes, it is essential that mesh gen-
eration engine creates nodes along the weir polyline. To achieve this it is necessary that each weir
polyline coincides exactly with corresponding arcs on SMS Conceptual Model.
1. For each calculation time interval, the model checks for each segment defined by two pair of
opposing cells (L, R) along the weir that at least one of the opposite cells is wet and that its
water surface elevation is above the crest elevation.
2. Then the model calculates the water elevation at each weir segment as:
where hcrest is the crest elevation and dw the segment water elevation.
7.10 Dam Breach Modeling 92
3. When the water surface levels on both sides is below the weir level, M AX (dL , dR ) ≤ dw ,
the velocity component normal to the weir segment direction is set to zero.
4. Otherwise the model calculates the normal discharge for the segment according to the water
levels on both sides.
5. The discharge is imposed on both the and cells.
The weir is simulated by assuming that the discharge per unit breadth q crossing the weir is
governed by the difference between the water surface level (d = h + z) on both sides of the weir,
referred to as dl upstream and dr downstream of the weir, and by the weir crest elevation, Hw .
Several situations are accounted for. In the case that both water elevations are below the weir
crest elevation the weir behaves as a solid wall and no flow crosses it. When dl < dr , without loss
of generality, two different flow situations can occur depending on the relative values of Hw , zl , zr ,
dl and dr . When Hw + max (zl , zr ) < min (dl , dr ), the discharge is given by
3/2
q = Cf (dr − dl ) (7.41)
with Cf the discharge coefficient.
• Coordinates (x,y) of the center of the breach, assuming z = zcrest , where zcrest is the initial
dam z-coordinate.
• Value of material angle α (assumed constant).
• Table (t, b(t), Hb (t)), being t =time, b =lower breach width, Hb =breach height.
As interenting particular cases, b(t) = 0 means the generation of a triangular breach. and α = 0
means the generation of a rectangular breach.
In general, the total discharge through the breach will be calculated with a law of the type:
Qb = KB(H)3/2 (7.43)
93 7 Hydraulic and Hydrologic Components
The discharge computed in (7.43) will be distributed among the cells included in the breach
top lenght B(Hb ):
Hb (t)
B(Hb ) = b(t) + 2 (7.44)
tgα
RiverFlow2D DIP provides a data input environment with panels that include all the non-
spatial options to run the RiverFlow2D model. The left column on the main window allows you
select modules, components, output options, etc. When you click on one of the elements, the
appropriate right side panel is activated. Each panel contains the data corresponding to each of
the RiverFlow2D data files. For example, the Control Data Panel has all the data of the .DAT file.
95
8.1 Control Data Panel (.DAT files) 96
RiverFlow2D DIP lets you select different model engines. Clicking on Model Selection opens the
dialog where you can select between RiverFlow2D Plus GPU, RiverFlow2D Plus or RiverFlow2D
FE. Note that to run the GPU version you need a special license as well as the appropriate GPU
hardware.
Table 8.2 – Time Control Data Frame on the Control Data Panel.
CONTROL NAME DESCRIPTION
Simulation time (hrs.) Total simulation time in hours.
Output Interval (hrs.) Time interval for output reporting.
CFL Courant-Friederich-Lewy condition (CFL). Set this num-
ber to a value in the (0,1] interval. By default CFL is set
to 1.0 which is the recommended value for maximum per-
formance. A few rare applications may require reducing
CFL to 0.5 or so avoid model oscillations in the model
results.
Hot start Use this option to restart the model from a previously
simulation.
Table 8.4 – Computation Control Data Frame on the Control Data Panel.
CONTROL NAME DESCRIPTION
Manning’s n factor Use this factor (XNMAN) to test the sensitivity of results
to the Manning’s n. Each element Manning’s n-value will
be multiplied by XNMAN. Over several simulations, using
this factor may speed up the model calibration. Default
is XNMAN = 1.
Manning’s n variable with Select this option to set Manning’s n as a function of
depth depth. The user must enter polygons over the mesh and
each polygon should have an associated file containing the
depth vs Manning’s n table.
Minimum depth for dry bed This parameter indicates the depth below which cell veloc-
ity will be assumed 0. By default it is set to -1 which will
allow the model to dynamically set the dry cell conditions
for depths smaller than 10−6 m.
Table 8.5 – Closed Boundary Conditions Data Frame on the Control Data Panel.
CONTROL NAME DESCRIPTION
Table 8.6 – Output Options Data Frame on the Control Data Panel.
CONTROL NAME DESCRIPTION
Text file output Reports results to ASCII text output files.
Continued on next page
99 8 RiverFlow2D Data Input Program (DIP)
Table 8.9 – Parameters on the Sediment transport Mode frame and buttons of the Sediment Transport
Panel (RiverFlow2D Plus model).
CONTROL NAME DESCRIPTION
Suspended sediment When this check box is selected, the model will com-
pute sediment concentrations using the suspended sedi-
ment transport component. See comment 1.
Bed load Sediment Selecting this check box will activate the bed load sediment
transport component. See comment 1.
Buttons
Open .SED* Opens an existing .SEDS or .SEDB files.
Save .SED* Saves the sediment data to .SEDS and .SEDB files.
101 8 RiverFlow2D Data Input Program (DIP)
Table 8.10 – Parameters on the Suspended Sediment transport frame of the Sediment Transport
Panel (RiverFlow2D Plus model).
CONTROL NAME DESCRIPTION
Equilibrium Concentration When this check box is selected, the model will compute
formula sediment concentrations using one of the following sus-
pended sediment transport formulas:
1. Bagnold (1966)lb
2. Van Rijn (1984a)
Settling Velocity Formula It is a unique formula that applies for all fractions. This
drop-down list includes the following formulas:
1. Rubey (1933)
2. Zhang (1961)
3. Zanke (1977)
4. Van Rijn (1984a)
5. Raudkivi (1990)
6. Julien (1998)
7. Cheng (1997)
8. Jimenez-Madsen (2003)
9. Wu-Wong (2006)
Table 8.11 – Parameters on the Bed Load Sediment transport frame of the Sediment Transport Panel
(RiverFlow2D Plus model).
CONTROL NAME DESCRIPTION
Sediment load formula Allows selection of one of the following sediment transport
formulas:
1. Meyer-Peter & Muller (1948)
2. Ashida (1972)
3. Engelund (1976)
4. Fernandez (1976)
5. Parker fit to Einstein (1979)
6. Smart (1984)
7. Nielsen (1992)
8. Wong 1 (2003)
9. Wong 2 (2003)
10. Camenen-Larson (1966)
Table
Density Sediment density (lb/f t3 or kg/m3 ).
Diameter D30 Sediment D30 size (m). 30% of the sediment is finer than
D30. Only used for Smart Formula.
Diameter Characteristic sediment size for this fraction (m).
Diameter D90 Sediment D90 size (m). 90% of the sediment is finer than
D90. Only used for Smart Formula.
Porosity Sediment porosity.
Shields Stress Critical Shield stress.
Friction Angle Sediment friction angle (degrees).
Fraction Fraction of material in bed. All fractions must add up to
1.
Factor Transport formula factor for each fraction. This factor
multiplies the result of the transport formula selected.
Buttons
Add Bed Load Fraction Used to add a new fraction. Up to 10 fractions may be
used.
Remove Selected Fraction Deletes the selected fraction.
Table 8.12 – Parameters on the Mud/Debris Flow Panel (RiverFlow2D Plus model).
CONTROL NAME DESCRIPTION
Flow resistance relation
1. Turbulent flow
2. Full Bingham
3. Simplified Bingham
4. Turbulent and Coulomb
5. Turbulent and Yield
6. Turbulent, Coulomb and Yield
7. Quadratic
8. Granular flow
Table 8.13 – Parameters on the Oil Spill on Land Panel (OilFlow2D Plus model).
CONTROL NAME DESCRIPTION
Fluid viscosity Viscosity (Pascal s or lbs/in2 ).
Fluid density Density (kg/m3 or lb/f t3 ).
See output file section (9.4.3) for output file content description.
107 8 RiverFlow2D Data Input Program (DIP)
See output file section (4.4.1) for output file content description.
Use the Plot button to generate a graph of the time series as shown in Figure 8.13.
Figure 8.15 shows the corresponding data entry controls that appear when selecting the first
row for Culvert1 that is a circular culvert.
To access the program use first select Tools from the Output group on the left panel of the
RiverFlow2D DIP.
Explanation and guidelines on how to use the 2dMRI utility is described in the 2dMRI manual
included in the documentation folder:
...\Documents\RiverFlow2Dv4\Documentation.
To activate this panel, first select the Tools from the Output group on the left panel of the
RiverFlow2D DIP.
The tutorial: “Using HEC-RAS elevation data” in the RiverFlow2D User’s Guide, describes in
detail the procedure to use this tool.
Input Data File Reference
9
Data files for non-spatial information required to run the RiverFlow2D model are created by the
Data Input Program and the spatial data is generated based on the .2dm file. All RiverFlow2D
input data files are in ASCII free-form format, which can be opened using any text editor or
spreadsheet program. In some instances it may be convenient to directly edit the data. However,
it is recommended to edit files with extreme caution, and only after having gained a thorough
understanding of RiverFlow2D file formats. This section explains the input data files, and the
parameters included in each file.
The RiverFlow2D installation program creates a folder with several example projects that can
be consulted to review the model data files. Depending on your operating system and settings,
this folder can be found in:
.\Documents\RiverFlow2D SMS\ExampleProjects.
121
122
RiverFlow2D data files will share the same name and will use the file extensions listed in the
table below. For example a run named Run1 will have files as follows: Run1.DAT, Run1.FED, etc.
The following table summarizes the data files used by RiverFlow2D model.
Line 12: Surface detention or minimum value of flow depth for dry areas.
HMIN
2
00000
0 0.5 0.25 0.25 8
10
1
1
0.9
1
-1
0
0
0
0
4
1
1. Setting the CFL (Courant Friederich-Lewy) or Courant number is critical for adequate sta-
bility and ensure mass conservation. RiverFlow2D explicit time scheme is conditionally
stable, meaning that there is a maximum time step above which the simulations will become
unstable. This threshold can be theoretically approximated by a Courant-Frederick-Lewy
9.2 Mesh Data 130
2. For variable time step simulations, RiverFlow2D estimates the maximum DT using the theo-
retical Courant-Frederick-Lewy (CFL) condition. Sometimes, the estimated DT may be too
high, leading to instabilities, and it may be necessary to reduce CFL to with a value less than
one to adjust it. Typical CFL values range from 0.3 to 1, but may vary project to project.
3. There are three initial conditions options. If IINITIAL = 0, the initial water elevation will
be a constant horizontal surface at the elevation given as INITIAL WSE. If INITIAL WSE
is = -9999 then the program will assign a constant water elevation equal to the highest bed
elevation on the mesh. If IINITIAL = 1, the whole computational mesh will be initially dry,
except at open boundaries where discharge is prescribed and depth > 0 is assumed for the
first time step. If IINITIAL = 2, initial water surface elevations are read from the .FED data
file for each node in the mesh. These initial elevations can be set creating a dataset file (.h5)
on SMS layer, and can be arbitrary. This last option allows modeling cases where part of
the domain is flooded (e.g. a reservoir) and other parts of the mesh are dry. The initial
velocities are always assumed to be equal to 0.
4. Use the IPROFILE option to allow RiverFlow2D to generate results along a polyline. The
polyline and other required data should be given in the Profiles file (.PROFILES), which is
defined later in this document.
5. Use this option to allow RiverFlow2D to generate results along prescribed cross sections. The
cross sections and other required data should be given in Cross Section file (.XSECS) which
is defined later in this document.
6. Use the XNMAN option to test the Manning’s n value sensitivity on the results. The
prescribed Manning’s coefficient assigned to each element will be multiplied by XNMAN.
This option is useful to test model sensitivity to Manning’s n during model calibration.
7. The model will create output files with maximum values of each output variable.
8. The user can specify an initial water surface elevation setting IINITIAL = 0 and entering
INITIAL WSE.
SMS assures that the .FED file will be created error free and consistent with the boundary condi-
tions and other mesh parameters. Manually editing this file may introduce unexpected errors.
367 1 WSE97out.txt
431 1 WSE97out.txt
...
This .IFL file has 1165 nodes on the boundary. Node 365 has a BCTYPE=1 (Water Surface
Elevation) and the time series of water surface elevations vs. time is in file WSE97out.txt.
6
INFLOW1.QVT
17
2
1
...
25
26
This .OBCP file has 2 open boundaries. The first open boundary is BCTYPE=12 corresponding to
Uniform Flow outflow. The uniform flow WSE vs Discharge table is included in file UNIF1.DATP,
and there are 24 nodes on the boundary. The second open boundary is BCTYPE = 6 corresponding
to inflow hydrograph where the Discharge vs time table is given in file INFLOW1.QVT, and there are
17 nodes on the boundary.
Table 9.4 – Variable Descriptions for the .IFL and .OBCP Files.
VARIABLE TYPE RANGE UNITS DESCRIPTION
BCTYPE I - - Code to indicate type of open
boundary. See Table 9.5 and
comment 1.
Continued on next page
9.2 Mesh Data 134
2. When imposing a single variable (water surface elevation, Q, U or V), the user must provide
an ASCII file with the time series for the corresponding variable. See section Boundary
Conditions Data Files for details on the format for one-variable boundary condition files.
3. When imposing two variables (water surface elevation and discharge, etc.), it is required to
provide an ASCII file with the time series for the variables. See section Boundary Conditions
Data Files for details on the format for two-variable boundary condition files.
4. When imposing water surface elevation it is important to check that the imposed value is
higher than the bed elevation. Even though RiverFlow2D can run with that condition, it
could lead to volume conservation errors.
5. A closed boundary condition is imposed by default on all boundary nodes. In this case,
the model calculates velocities and water surface elevations for all nodes on the boundary
depending on the value of the ISPLIPBC parameter. For example ISLIPBC = 1 will impose
slip conditions setting zero-flow across the boundary. Tangential flow is free corresponding
to a slip condition.
6. When using a single valued stage-discharge condition the model first computes the discharge
on the boundary then interpolates the corresponding water surface elevation from the rating
table and imposes that value for the next time step. In case the boundary is dry, it functions
as a free condition boundary (see comment 7). Water surface elevations are imposed only
on wet nodes. This condition requires providing an ASCII file with the table values entries.
See section Boundary Conditions Data Files for details on the file format. In general it is
preferable to use stage hydrograph rather than stage-discharge condition. In most small
slope rivers, the stage-discharge relationship is affected by hysteresis. In other words, the
stage-discharge curve is looped with higher discharges occurring on the rising limb than
on the rescission limb of the hydrograph. This is mainly caused by the depth gradient
in the flow direction that changes in sign throughout the hydrograph. In practice, this
implies that there can be two possible stages for the same discharge. If the stage-discharge
relationship is not well known or if it just computed assuming steady state uniform flow,
9.2 Mesh Data 136
it may lead to considerable errors when used as downstream boundary condition. That it
is why it is often preferred to use the stage hydrograph for that purpose. However, such
hydrograph may not be available to study changes in the river and evaluating proposed
conditions. For those cases, it is useful to use a stage-discharge relationship, preferably
measured over an extensive range of discharges. When this relationship is not available, one
option would be to assume steady state flow to determine a single-value rating curve. Since
this condition may generate wave reflection that can propagate upstream, it is important to
locate the downstream boundary on a reach sufficiently far from the area of interest, therefore
minimizing artificial backwater effects. Unfortunately, there is no general way to select such
place, but numerical experimenting with the actual model will be necessary to achieve a
reasonable location.
NOTE: loop stage-discharge relationships are not implemented in this RiverFlow2D version.
7. On free outflow condition boundaries, the model calculates velocities and water surface eleva-
tions applying the full equations from the internal elements. No specific values for velocities
or depths are imposed per se on these nodes. In practice this is equivalent to assuming that
derivatives of water surface elevations and velocities are 0. In subcritical flow situations, it
is advisable to use this condition when there is at least another open boundary where WSE
or stage-discharge is imposed.
8. When using a single valued stage-discharge condition on internal sections, the model first
computes the discharge across the boundary then interpolates the corresponding water sur-
face elevation from the rating table, imposing that value for the next time step for all nodes
on the internal boundary. This condition requires providing an ASCII file with the table
values entries. See section Boundary Conditions Data Files for details on the file format.
9. When imposing a water and sediment discharge, it is required to provide an ASCII file with
the time series for water discharge and volumetric sediment discharge for each of the frac-
tions. Note that sediment discharge is always expected in volume per unit time. See section
Boundary Conditions Data Files for details on the format for multiple-variable boundary
condition files.
10. The user must provide a file with the energy slope S0 for the corresponding boundary. This
file will only contain a single value S0 . The model will use S0 , Manning’s n, and discharge to
create a rating table from which water surface elevations will be imposed as a function of the
computed outflow discharge. The rating table is calculated every 0.05 m (0.16 ft.) starting
from the lowest bed elevation in the outflow cross section up to 50 m (164 ft.) above the
highest bed elevation in the section. If S0 = −999, the model will calculate the average bed
slope perpendicular to the boundary line. Please, note than when letting the model calculate
the average bed slope, it uses the elevations on the elements adjacent to the boundary line,
which may result in adverse slopes or slopes that do not capture the general trend the reach.
11. This boundary condition is similar to the BCTYPE = 6 for inflow water discharge. However,
in this case, instead of converting the discharge into velocities that are imposed on all the
inflow nodes; the model creates sources on all the elements adjacent to the boundary line.
The condition then can be visualized as if the given discharge enters over the inflow elements.
For each time, the model evenly divides the discharge between all the inflow elements. For
example if there are Ne inflow elements and the imposed discharge is Qin, each element will
receive a discharge equal to Qin/Ne. The water volume will naturally flow away from the
inflow depending on the bed slopes, etc. Care must be taken when the inflow boundary
137 9 Input Data File Reference
elements have lower bed elevations than the surrounding elements. When imposing this
condition the user must provide an ASCII file with the discharge time series. See section
Boundary Conditions Data Files for details on the format for one-variable boundary condition
files.
NNODESBOUNDARY lines containing the list of boundary nodes in counter clockwise direction.
BOUNDARYNODE (1:NNODESBOUNDARY)
The next lines are only used if there are islands in the mesh.
Start of boundary parameter indicator for each island or internal closed contour.
IBOUNDARYID
In this example the external boundary has 132 nodes and there is one island with 34 nodes.
9.3.1 Bridges
RiverFlow2D provides four options to account for bridge piers. The most common option is to
create the pier plan geometry generating a 2D triangular-cell mesh that represents each pier as a
solid obstacle. In that case, the model will compute the flow around the pier and account for the
pier drag. This would be the preferred approach when the user needs to know the detailed flow
around the piers and the flow does not overtop the bridge deck. However, the resulting mesh may
have very small elements, leading to increasing computer times.
The second option (Bridge Piers) is a simplified formulation that does not require defining the
mesh around the piers, but will compute the pier drag force based on geometric data. This would
be the preferred approach when the flow does not overtop the bridge deck and the user does not
need to have detailed depiction of the flow around the piers but needs to account for the general
effect that the pier would have on the flow.
The third option represented in the Bridges component is a comprehensive bridge hydraulics
computation tool that does not require capturing bridge pier plan geometry in detail, therefore
allowing longer time steps, while allowing calculating the bridge hydraulics accounting for arbitrary
plan alignment, complex bridge geometry, free surface flow, pressure flow, overtopping, combined
139 9 Input Data File Reference
pressure flow and overtopping, and submergence all in 2D. This is the recommended option for
most bridges.
There is a fourth option using the Internal Rating Table component, but for most applications
it is recommended to use one of the above since they better represent the bridge hydraulics.
Bridge Id.
BRIDGE ID
Table 9.9 – Variable Descriptions for the bridge cross section geometry file.
VARIABLE TYPE RANGE UNITS DESCRIPTION
BEDELEV R - m or f t Bed elevation. Must be the low-
est elevation for all polylines at a
given point.
DECKELEV R - m or f t Elevation of the bridge deck.
Must be the highest elevation for
all polylines at a given point.
NP I - >1 Number of points defining cross
section polylines.
STATION R - m or f t Distance from leftmost point
defining cross section polyline.
All polylines points must have a
common station.
ZLOWER R - m or f t Elevation of lower polyline.
ZLOWER must be larger or
equal to BEDELEV and smaller
or equal to LOWCHORD for a
given point. The space between
BEDELEV and ZLOWER is a
blocked area to the flow. The
space between ZLOWER and
LOWCHORD is open space.
If the bridge has no holes,
ZLOWER must be identical to
BEDELEV.
LOWCHORD R - m or f t Elevation of the lower bridge
deck. LOWCHORD must be
larger or equal to ZLOWER and
smaller or equal to DECKELEV
for a particular point. The space
between LOWCHORD and DE-
CELEV is a blocked area to the
flow.
Note: R = Real variable. I = Integer variable.
143 9 Input Data File Reference
Table 9.11 – Drag Coefficients for Bridge Piers. Adapted from Froehlich (2003).
PIER PLAN SHAPE AND APPROACH VELOCITY DRAG COEFFICIENT CD
2. The drag coefficient CD is related to the drag force though the following formula:
1
FD = CD ρU 2 AP (9.2)
2
where CD is the pier drag coefficient, is the water density, U is the water velocity, and AP
is the pier wetted area projected normal to the flow direction.
To account for the drag force that the pier exerts on the flow, RiverFlow2D converts it to the
distributed shear stress on the element where the pier centroid coordinate is located. The
resulting pier shear stress expressions in x and y directions are as follows:
1 p AP
τpx = CD ρU U 2 + V 2 (9.3)
2 Ae
145 9 Input Data File Reference
1 p AP
τpy = CD ρV U 2 + V 2 (9.4)
2 Ae
INVERT Z1
INVERT Z2
Where INVERT Z1 and INVERT Z2 are the invert elevations for the inlet and outlet respectively.
7
0 0.20
0.1 1.00
1.00 36.09
2.00 60.00
3.00 84.78
4.00 110.01
100.00 110.02
5.0
1.0
Nb
Ke
nc
Kp
M
Cp
Y
m
If CulvertType=1
Hb
Base
Else if CulvertType=2
Dc
INVERT Z1
INVERT Z2
This example culvert characteristics data file indicates that the culvert one barrel (Nb =1),
Ke=0.4, nc=0.012, Kp=1, cp =1, M =1.1, Y=0.6, m=-0.5, and Dc=0.10, INVERT Z1=5.0 and
INVERT Z2 = 1.0.
Table 9.15 – Manning’s n roughness coefficients for various culvert materials. Adapted from Froehlich
(2003).
Culvert barrel material Entrance description Manning’s n
nc
Good joints, smooth walls 0.012
Concrete Projecting from fill, square-cut end 0.015
Poor joints, rough walls 0.017
2-2/3 inch × 1/2 inch corrugations 0.025
6 inch × 1 inch corrugations 0.024
5 inch × 1 inch corrugations 0.026
Corrugated metal
3 inch × 1 inch corrugations 0.028
6 inch × 2 inch corrugations 0.034
9 inch × 2 1/2 inch corrugations 0.035
Table 9.17 – Culvert inlet control formula coefficients. Adapted from Froehlich (2003).
Barrel Barrel Inlet description* K’ M c’ Y
material shape
Concrete Circular Headwall; square edge 0.3153 2.0000 1.2804 0.6700
Concrete Circular Headwall; grooved edge 0.2509 2.0000 0.9394 0.7400
Concrete Circular Projecting; grooved edge 0.1448 2.0000 1.0198 0.6900
Cor. Circular Headwall 0.2509 2.0000 1.2192 0.6900
metal
Cor. Circular Mitered to slope 0.2112 1.3300 1.4895 0.7500
metal
Cor. Circular Projecting 0.4593 1.5000 1.7790 0.5400
metal
Concrete Circular Beveled ring; 45 bevels 0.1379 2.5000 0.9651 0.7400
Concrete Circular Beveled ring; 33.7 bevels 0.1379 2.5000 0.7817 0.8300
Concrete Rectangular Wingwalls; 30 to 75 flares; 0.1475 1.0000 1.2385 0.8100
square edge
Concrete Rectangular Wingwalls; 90 and 15 flares; 0.2242 0.7500 1.2868 0.8000
square edge
Concrete Rectangular Wingwalls; 0 flares ;square 0.2242 0.7500 1.3608 0.8200
edge
Concrete Rectangular Wingwalls; 45 flare; beveled 1.6230 0.6670 0.9941 0.8000
edge
Concrete Rectangular Wingwalls; 18 to 33.7 flare; 1.5466 0.6670 0.8010 0.8300
beveled edge
Concrete Rectangular Headwall; 3/4 inch cham- 1.6389 0.6670 1.2064 0.7900
fers
Concrete Rectangular Headwall; 45 bevels 1.5752 0.6670 1.0101 0.8200
Concrete Rectangular Headwall; 33.7 bevels 1.5466 0.6670 0.8107 0.8650
Concrete Rectangular Headwall; 45 skew; 3/4 in 1.6611 0.6670 1.2932 0.7300
chamfers
Continued on next page
9.3 Component data 150
2. For CulvertTypes 1 and 2, the model will calculate culvert discharge for inlet and outlet
9.3 Component data 152
control using the FHWA procedures (Norman et al.,1985) that were later restated in dimen-
sionless form by Froehlich (2003). See User’s Guide for more details.
Gate Id
GATES ID
4031 279
4029 277
4027 156
4026 82
4024 16
Time Aperture.
TIME APERTURE
IRT NPL line groups containing the IRT polyline ID, the number of vertices defining each polyline,
the IRT boundary condition type (always equal to 19 in this version), the Rating Table file name,
followed by the list of polyline coordinate vertices as shown:
IRT ID
IRT NV IRT BCTYPE IRT FILENAME
X IRT(1) Y IRT(1)
X IRT(2) Y IRT(2)
...
X IRT(IRT NV) Y IRT(IRT NV)
This file has 2 internal rating table polylines, the ID of the first one is IRT A, which has 4
vertices, BCTYPE 19 and file name IRT A.DAT.
0.00899
3.5
2200.
In this example there are 7 entries in the Temperature Time series. For the initial time 0.0
hours the temperature is 60◦ F , and so on.
In this example there are 2 lines in the Temperature-Viscosity-Density table file. For 50◦ F the
viscosity is 1.38E-06 lb − s/in2 and the density 50 lb/f t3 .
Line 6 to NPOLLUTANT MAX+6: First order reaction rate coefficient matrix K(I,J)
K(1,1) K(1,2) ... K(1,NPOLLUTANT MAX)
K(2,1) K(2,2) ... K(2,NPOLLUTANT MAX)
...
K(NPOLLUTANT MAX,1) K(NPOLLUTANT MAX,2) ... K(NPOLLUTANT MAX,NPOLLUTANT MAX)
Line 1: Number of zones defined by polygons where Manning’s n variable with depths are defined.
NNZONES
NRZONES group of lines containing Manning’s n variable with depth data file for each zone
MANNNFILE
4
25.0 25.0
25.0 75.0
75.0 75.0
75.0 25.0
Manning2.txt
4
25.0 125.0
25.0 175.0
75.0 175.0
75.0 125.0
In this example, there are two polygons. The Manning’s n data file for the first polygon is
Manning1.txt and the polygon is defined by four vertices.
ND lines containing:
DEPTH MANNINGS N
Table 9.28 – Variable Descriptions for the Manning’s n Cariable with Depth Data File.
VARIABLE TYPE RANGE UNITS DESCRIPTION
DEPTH R ≥0 m or f t Flow depth. See comment 1.
MANNINGS N R ≥0 - Manning’s n corresponding to
DEPTH. See comment 1.
ND I - - Number values in file.
Note: I = Integer variable. R = Real variable.
9.3.13.5 Comments for the Mannign’s n variable with depth data file
1. To calculate the Manning’s n over the mesh, the model will first identify the polygon over
each cell and then will use the interpolated n value for cell depth from the table corresponding
to the polygon. In the example above, for all depth between 0.3 and 1, Manning’s n will be
obtained by linear interpolation between 0.1 and 0.03.
2. The user should provide a “DefaultManningsn.dat” file in the project folder and the program
will apply the data contained in that file to the complementary area to the polygons provided.
If the “DefaultManningsn.dat” does not exist, the model will apply a default value of 0.035
to the areas not covered by Manning’s n polygons.
Line 1: Number of zones defined by polygons where rainfall time series are defined.
NRZONES
NRZONES group of lines containing hyetograph and evaporation data file for each zone
RAINEVFILE[i]
25.0 125.0
25.0 175.0
75.0 175.0
75.0 125.0
In this example, there are two polygons. The rainfall and evaporation data file for the first
polygon is hyeto1.txt and the polygon is defined by four vertices.
Table 9.30 – Variable Descriptions for the Hyetograph and Evaporation Data File.
VARIABLE TYPE RANGE UNITS DESCRIPTION
EVAPORATION R ≥0 mm/hr Evaporation intensity. See comment
or in/hr 1.
NPRE I - - Number of times in rainfall and evap-
oration time series.
RAINFALL R ≥0 mm/hr Rainfall intensity. See comment 1.
or in/hr
TIME R >0 hours Time interval
Note: I = Integer variable. R = Real variable.
9.3 Component data 164
Line 1: Number of zones defined by polygons where infiltration parameters are defined.
NIZONES
In this example, there are two polygons. The infiltration data file for the first polygon is
inf1.inf and the polygon is defined by four vertices.
165 9 Input Data File Reference
3. If the user has a “DefaultInfiltration.dat” file in the project folder, the program will apply
the data contained in that file to the complementary area to the polygons provided.
In this example the infiltration loss method is set to 1 corresponding to the Horton model.
There are 3 parameters as follows: K = 8.3E-04, fc = 3.47E-06 and f0 = 2.22E-5.
9.3 Component data 166
0.0039370 0.0039370
0.0039000 0.0039370
0.0039370 0.0039370
0.40 0.40
0.047 0.047
35.00 35.00
0.50 0.50
1.00 1.00
NSOURCES groups of lines containing source/sink point identification text, name of the file con-
taining the discharge time series, and the coordinates of the point as follows:
SOURCEID
ISFILENAME
X S(I) Y S(I)
...
This file has 2 sources/sinks. The first one is named IntakeA located at coordinate: X = 799019.633
and Y = 309402.572 and discharge data is on file: Intake.txt.
2. To model inflows use positive discharge values and to model outflows use negative values.
NWEIRS group of lines including weir ID, number of vertices defining each weir polyline, the weir
coefficient followed by the coordinates each vertex as shown:
WEIR ID
NV CF WRCRESTELEV
X W(1) Y W(1) WRCREST(1)
X W(2) Y W(2) WRCREST(2)
...
X W(NV) Y W(NV) WRCREST(NV)
This file has 2 weirs. The first one is named WEIR A and is defined by a polyline with 4
vertices. Weir discharge coefficient is equal to 3.0.
Q = Cf LH 3/2 (9.5)
where L is the distance between nodes, H is the total head upstream of the polyline segment
and Cf is the discharge coefficient. The model checks for submergence and it occurs Cf will
be corrected according to the correction factor defined by (FHWA, 2001).
2. Weir polylines should be defined avoiding abrupt direction changes (e.g. ≥ 90 degree turns),
because such angles may create errors in the algorithm that identifies the nodes that lie over
the polyline.
Line 1: Number of zones defined by polygons where wind velocity time series are defined.
NWZONES
NWZONES group of lines containing hyetograph and evaporation data file for each zone.
WINDFILE
In this example, there are two polygons. The Cd coefficient is set to 0.009 and the wind density
to 1.225 kg/m3 . The wind velocity file for the first polygon is wind1.txt and the polygon is
defined by four vertices.
2. If the user has a “DefaultWind.dat” file in the project folder, the program will apply the
data contained in that file to the complementary area to the polygons provided.
175 9 Input Data File Reference
NOBSPOINTS groups of lines containing the observation point ID, and coordinate of each point:
ObsID
X OP(I) Y OP(I)
...
799146.926 309430.876
PointC
799721.8608 309041.615
This .OBS file has three points. The first point is named PointA and has coordinates: X=798798.380
Y=309627.950.
Line 6: Transparency.
USEBACKIMAGE
Line 7: Transparency.
TRANSP
Example:
179 9 Input Data File Reference
2.05
0.00
0.00
-2.05
795944.99
310049.73
In this example, 2.05 is the pixel size in x-direction, rotation in x and y axes is 0.00, pixel
size in y direction is 2.05 (shown in negative), x-coordinate of upper left pixel is 795944.99
and y-coordinate of upper left pixel is 310049.73.
The following table indicates the supported image formats and their corresponding world file
extensions.
Table 9.41 – Supported image formats and their corresponding world file extensions.
IMAGE FILE FOR- WORLD FILE EXTENSION
MAT
.BMP .BMPW, .BPW
.GIF .GFW, .GIFW, .WLD
.PNG .PGW, .PNGW, .WLD
.TIF, .TIFF .TFW, .WLD
3. Vector field, depth and water surface elevation text files are used for creating vector field
shapefiles. Using the IVSF = 1 option, will allow the model to generate ESRI Shapefile
support files for each report interval. This option is not used in the RiverFlow2D SMS.
9.4 Output control data 180
NPROFILES group of files including: Profile ID, number of vertices in profile I, the number of
intervals to divide each profile, and coordinates for each vertex in polyline.
PROFILEID
NVERTICES PR(I) ND PR
X PRF(I), Y PRF(I)
...
This file indicates there are 2 profiles. First profile ID is: ProfileA which is defined with a
2-vertex polyline and will be divided in 10 segments.
NCROSS SECTIONS groups of lines containing the cross section ID, the number of vertices defin-
ing the cross section (always equal to 2), the number of intervals to divide the cross section and
the list of coordinates of initial and final point in cross section:
XSECID
NPXSEC ND CS
X1 CS(I) Y1 CS(I)
X2 CS(I) Y2 CS(I)
This .XSECS file indicates there are 3 cross sections. The first one has ID = CrossSectionA and
will be divided in 40 segments.
In this example .EXP file, there are 11086 elevation data points, one parameter per point (the
elevation for each point).
the available data is in Lat/Lon format, it will need to be converted to UTM or plane Carte-
sian coordinates before importing coordinates into SMS . Presently RiverFlow2D does not
provide tools to make this conversion.
2. Elevation values should be given in the same units as the corresponding coordinates.
9.6.1.1 Example of the Boundary Condition File for One Variable Time Series
The following example shows an inflow hydrograph where NDATA is 7 and there are 7 lines with
pairs of time and discharge:
7
0. 20.
1. 30.
1.3 50.
2. 90.
4. 120.
5. 200.
7. 250.
9.6 Boundary conditions data 184
Where VARIABLE1(I) and VARIABLE2(I) depend on the boundary condition type as follows:
10
0. 20. 1420.
1. 30. 1421.5
1.3. 50. 1423.
...
7. 250. 1420.
8.1 110. 1426.
185 9 Input Data File Reference
Where VARIABLE1(I) ... VARIABLEN(I) depend on the boundary condition type as follows:
10
0. 20. 0.001 0.002
1. 30. 0.002 0.005
1.3. 50. 0.003 0.010
...
7. 250. 0.01 0.015
8.1 110. 0.005 0.009
10. 60. 0.004 0.007
9.6 Boundary conditions data 186
Where STAGE(I) is water surface elevation and Q(I) is the corresponding discharge.
3.75 103.67
4.00 110.01
7
0 0.20
0.1 1.00
1.00 36.09
2.00 60.00
3.00 84.78
4.00 110.01
100.00 110.02
189
10.1 Output File Overview 190
This file also reports the list of acute cells that have an internal angle of less than 22.5 degrees.
If there are acute cells, the model will give an error message and will not be able to execute.
For example: RESvsT_0000010.OUTI is the file name for time series results of node 10. An example
of this file is shown below.
10.2 General Output Files 202
18000 state5.out
Note that the stateN.out files are named sequentially. For instance, stateN.out corresponds to
the Nth report interval.
RiverFlow2D can be restarted from the any existing report time by reading the initial condi-
tions from the stateN.out file indicated in the 2binitialized.HOTSTART file. To restart from a
time different from the last one calculated, just edit the 2binitialized.HOTSTART file and enter
the desired time in seconds and corresponding stateN.out file name that is to be used as initial
conditions. For example, to hot start from hour 3 (10800 seconds) and assuming that the report
interval is 0.5 hours, the 2binitialized.HOTSTART file should contain the following entry:
203 10 Output File Reference
10800 state6.out
The hot start option is often useful to establishing initial conditions common to a series of
simulations for various return periods. For instance, to generate your initial state, you could run
the model with a constant discharge inflow until the model converges to a steady state. Assuming
that the final report time corresponds to the state20.out file, you can edit the 2binitialized.
HOTSTART file as shown:
0 state20.out
Then when you run the RiverFlow2D model using the hot start option, the model will start
assuming that the data in the state20.out file will define the initial conditions. You may want to
keep the 2binitialized.HOTSTART and state20.out files in a separate directory and copy them
to the project folder for each desired scenario.
Please, keep in mind that the stateN.out files are tied to the mesh you use, so if you modify
the mesh in any way, you will need to use the stateN.out corresponding to that mesh.
When running only hydrodynamics the .XSECE and .XSECI files will display the cross section
water discharge. When running sediment transport, in addition to the water discharge these files
will report the total sediment discharge in f t3 /s or m3 /s.
10.4 Profile and Cross Section Output Files 208
Mesh1_RESULTS.HDF5
Mesh1_MAXIMUM_RESULTS.HDF5
Where Mesh1 is the project name. The content of the Result at Nodes files is described on table
10.6.
TIME_ENG_0001_12_01_34.EXP
corresponds to a file in English units for time: 1 day, 12 hours, 1 minute and 34 seconds.
The format for these files is as follows. The first line indicates the number of node parameters
(19 by default). Then follows NNODES lines with results for each node in the triangular element
mesh as shown:
CONC_DDDD_HH_MM_SS.EXP
CONC_0001_12_01_34.EXP
corresponds to a file in English units for time: 1 day, 12 hours, 1 minute and 34 seconds.
The format for these files is as follows. The first line indicates the number of solutes used in
the PL run (NPOLLUTANTS MAX). Then follows NNODES lines with results for each node in
the triangular element mesh as shown:
3
1 0.000000 0.000000 0.000000
2 0.000000 0.000000 0.000000
3 0.000000 0.000000 0.000000
4 0.202378 0.000000 0.000000
5 0.326602 0.000000 0.000000
6 0.291721 0.000000 0.000000
7 0.000000 0.000000 0.000000
...
ST_DDDD_HH_MM_SS.EXP
ST_0001_12_01_34.EXP
corresponds to a file in English units for time: 1 day, 12 hours, 1 minute and 34 seconds.
The format for these files is as follows. The first line indicates the number of suspended sediment
fractions used in the ST run (NSSNFRAC). Then follows NNODES lines with results for each node
in the triangular element mesh as shown:
2
1 0.000000 0.000000
2 0.000000 0.000000
3 0.000000 0.000000
4 0.000000 0.000000
5 0.000000 0.000000
6 0.000000 0.000000
7 0.000000 0.000000
10.5 Output Files for GIS Post-processing 216
8 0.000000 0.000000
9 0.000287 0.000287
10 0.000456 0.000456
11 0.000356 0.000356
12 0.000334 0.000334
13 0.000000 0.000000
...
The format for these files is as follows. The first line indicates the number of node parameters (5
by default). There follows NNODES lines with velocity module, depth and water surface elevation
for each node as shown:
<ProjectName>_time2depths.exp
The format for these files is as follows. The first line indicates the number of node parameters
(3 by default). There follows NNODES lines with time to 0.30 m (1 ft), time to 1 m (3 ft), and
time to maximum depth for each node as shown in Table 10.13. When the node remains dry or
depth is below 0.30 m (1 ft) the reported value is -1.
Velocity field files use the .EXP extension and are named as follows:
VelField_SF_TIME_DDDD_HH_MM_SS.EXP
Water surface elevation files use the .EXP extension and are named as follows:
WSE_SF_TIME_TIME_DDDD_HH_MM_SS.EXP
For further details about the USBR Hazard classification, consult the RiverFlow2D User’s
Guide and USBR (1988).
CELL_TIME_ENG_0001_12_01_34.TEXTOUT
corresponds to a file in English units for time: 1 day, 12 hours, 1 minute and 34 seconds.
The format for these files is as follows. The first line contains the number if cells (NELEM)
and the number of cell parameters which is 16. Then it follows NELEM lines with results for each
cell in the triangular element mesh as shown:
The format for these files is as follows. The first line contains the number of cells (NELEM),
and the number of cell parameters which is 4 by default, but 8 if the run was made with the
Sediment Transport Module. There follows NELEM lines with velocity module, depth and water
surface elevation for each cell as shown:
<ProjectName>_TIME2DEPTHS_CELLS.TEXTOUT
The format for these files is as follows. The first line indicates the number of cells (NELEM)
and the number of cell parameters (5 by default). For the file in Metric Units there follows NELEM
lines with time to 0.30 m, time to 0.5 m, time to 1 m, time to maximum depth, and total inundated
time for each cell as shown in Table 10.17. When the cell remains dry or depth is below 0.30 m
the reported value is -1. For the file in English Units there follows NELEM lines with time to 1
ft, time to 2 ft, time to 3 ft, time to maximum depth, and total inundated time for each cell as
shown in Table 10.17. When the cell remains dry or depth is below 1ft the reported value is -1.
10.7 Graphic Output Files 220
<ProjectName>_CELLS_HAZARD.TEXTOUT
The format for these files is as follows. The first line indicates the number of cells (NELEM),
and the number of cell parameters (11). There follows NELEM lines with the hazard intensities
for each cell as shown:
Figure 10.1 – Output color coded by depth and with velocity vectors
10.7 Graphic Output Files 222
Data Input Program Tools
11
This section describes various utilities that are available through the RiverFlow2D Data Input
Program.
223
11.1 HEC-RAS Data Extraction Tool 224
The HEC-RAS Data Extraction Tool is accessible from the Tools panel of the Data Input Program.
225
226
Cheng, N. (1997), ‘Simplified settling velocity formula for sediment particle’, J. Hydraulic Eng.,
ASCE. 123(2), 149–152.
Cunge, J., Holly, F. & Vervey, A. (1980), Practical Aspects of Computational River Hydraulics,
Pitman: London.
De Vriend, H., Zyserman, J., Nicholson, J., Roelvink, J., Péchon, P. & Southgate, H. (1993),
‘Medium-term 2DH coastal area modelling’, Coastal Engineering 21(1-3), 193–224.
Engelund, F. & Fredsoe, J. (1976), ‘A sediment transport model for straight alluvial channels.’,
Nordic Hydrology 7, 293–306.
Engelund, F. & Hensen, E. (1967), A monograph on sediment transport to alluvial streams, Report,
Copenhagen: Teknique Vorlag.
Fernandez-Luque, R. & van Beek, R. (1976), ‘Erosion and transport of bed sediment’, Journal of
Hydraulic Research, IAHR 14(2), 127–144.
FHWA (1978), Hydraulics of bridge waterways, Report EPD-86-101, FHWA (Federal Highway
Administration).
Froehlich, D. (2003), Users Manual for FESWMS FST2DH Two-dimensional Depth-averaged Flow
and Sediment Transport Model, Report No. FHWA-RD-03-053, Federal Highway Administra-
tion, Washington, DC.
Garcı́a, R., Espinoza, R., Valera, E. & González, M. (2006), ‘An explicit two-dimensional finite
element model to simulate short- and long-term bed evolution in alluvial rivers’, Journal of
Hydraulic Research 44(6), 755–766.
Garcı́a, R., González, N. & O’Brien, J. (2009), Dam-break flood routing. Chapter 4 in: Dam-Break
Problems, Solutions and Case Studies, WIT Press, Southampton-Boston.
Garde, R. & Ranga Raju, K. (1985), Mechanics of Sediment Transportation and Alluvial Stream
Problems, John Wiley & Sons, New York.
Horton, R. (1933), ‘The role of infiltration in the hydrologic cycle’, Trans. Am. Geophys. Union
14, 446 – 460.
Jimenez, J. & O.S., M. (2003), ‘A simple formula to estimate settling velocity of natural sediments’,
J. Waterway, Port, Coastal, Ocean Eng 129(2), 70–78.
Juez, C., Murillo, J. & P., G.-N. (2014), ‘Numerical assessment of bed-load discharge formulations
for transient flow in 1d and 2d situations’, Journal of Hydroinformatics 15(4), 1234–1257.
Jury, W. & Horton, R. (2004), Soil Physics, John Wiley and Sons.
Karim, F. (1988), ‘Bed material discharge prediction for non-uniform bed sediments’, Journal of
Hydraulic Engineering, ASCE 124(6), 597–604.
227 12 References
LeVeque, R. (2002), Finite-Volume Methods for Hyperbolic problems, Cambridge University Press,
New York.
Mein, R. G. & Larson, C. L. (1973), ‘Modeling infiltration during a steady rain’, Water Resources
Research 9(2), 384–394.
Meyer-Peter, E. & Muller, R. (1948), ‘Formulas for bed-load transport’, Proc. of the Second Meet-
ing. IAHR, Stockholm, Sweden pp. 39–64.
Mishra, S. & Singh, V. (2003), Soil Conservation Service Curve Number (SCS-CN) Methodology,
Kluwer Academic Publishers.
Morvan, H., Knight, D., Wright, N., Tang, X. & Crossley, A. (2008), ‘The concept of roughness
in fluvial hydraulics and its formulation in 1-d, 2-d and 3-d numerical simulation models’,
Journal of Hydraulic Research 46(2), 191–208.
Murillo, J., Burguete, J., Brufau, P. & Garcı́a-Navarro, P. (2005), ‘Coupling between shallow water
and solute flow equations: analysis and management of source terms in 2D’, International
Journal for Numerical Methods in Fluids 49(3), 267–299.
Murillo, J., Burguete, J., Brufau, P. & Garcı́a-Navarro, P. (2007), ‘The influence of source terms
on stability, accuracy and conservation in two-dimensional shallow flow simulation using tri-
angular finite volumes’, International Journal of Numerical Methods in Fluids 54, 543–590.
Murillo, J. & Garcı́a-Navarro, P. (2010a), ‘An exner-based coupled model for two-dimensional
transient flow over erodible bed’, Journal of Computational Physics 229(23), 8704–8732.
Murillo, J. & Garcı́a-Navarro, P. (2010b), ‘Weak solutions for partial differential equations with
source terms: Application to the shallow water equations’, Journal of Computational Physics
229(11), 4327–4368.
Murillo, J., Garcı́a-Navarro, P. & Burguete, J. (2008), ‘Analysis of a second-order upwind method
for the simulation of solute transport in 2d shallow water flow’, International Journal for
Numerical Methods in Fluids 56(6), 661–686.
Murillo, J., Garcı́a-Navarro, P. & Burguete, J. (2012), ‘Wave riemann description of friction terms
in unsteady shallow flows: Application to water and mud/debris floods’, Journal of Compu-
tational Physics 231, 1963–2001.
Naef, D., Rickenmann, D., Rutschmann, P., & McArdell, B. (2006), ‘Comparison of flow resis-
tance relations for debris flows using a one-dimensional finite element simulation model’, Nat.
Hazards Earth Syst. Sci. 6, 155–156.
Nielsen, P. (1992), Coastal Bottom Boundary Layers and Sediment Transport, Advanced Series on
Ocean Engineering. World Scientific Publishing.
Norman, J. M., Houghtalen, R. J. & Johnston, W. J. (1985), Hydraulic design of highway culverts,
Hydraulic Design Series No. 5 FHWA-IP-85-15, Federal Highway Administration, Washington,
D.C.
Parker, G. (1979), ‘Hydraulic geometry of active gravel rivers’, Journal of Hydraulic Engineering,
ASCE 105(9), 1185–1201.
Parker, G., Klingeman, P. & McLean, D. (1982), ‘Bed load and size distribution in paved gravel
bed streams’, Journal of Hydraulic Engineering 108(4), 544–571.
228
Ratia, H., Murillo, J. & Garcı́a-Navarro, P. (2014), ‘Numerical modelling of bridges in 2D shallow
water flow simulations’, International Journal for Numerical Methods in Fluids 75(4), 250–
272.
Raudkivi, A. (1990), Loose boundary hydraulics, Pergamon Press, Inc., Tarrytown, N.Y.
Rawls, W. & Brakensiek, D. (1983), A procedure to predict green and ampt infiltration parameters,
in ‘Proceeding of ASAE Conferences on Advances in Infiltration’, Chicago, Illinois, pp. 102–
112.
Rawls, W., Yates, P. & Asmussen, L. (1976), Calibration of selected infiltration equation for the
georgia coastal plain, Report ARS-S-113, Agriculture Research Service.
Roe, P. (1981), ‘Approximate riemann solvers, parameter vectors, and difference schemes’, Journal
of Computational Physics 43(2), 357 – 372.
Rogers, B. & Borthwick, A.G.L. Taylor, P. (2003), ‘Mathematical balancing of flux gradient and
source terms prior to using roe’s approximate riemann solver’, J. Comput. Phys. 192(2), 422–
451.
Rosatti, G., Murillo, J. & Fraccarollo, L. (2008), ‘Generalized roe schemes for 1D two-phase,
free-surface flows over a mobile bed’, J. Comput. Phys. 227(24), 10058–10077.
Rubey, M. (1933), ‘Settling velocities of gravel, sand and silt particles’, Amer. J. Soc. 225, 325–338.
Sha, Y. (1965), Introduction to Sediment Dynamics, Industry Press, 302-310.
Smart, G. (1984), ‘Sediment transport formula for steep channels’, J. Hydr. Eng., ASCE.
111(3), 267–276.
Te Chow, V. (1959), Open Channel Hydraulics, McGraw-Hill, New York.
Te Chow, V., Maidment, D. & Mays, L. (1988), Applied Hydrology, McGraw-Hill Civil Engineering
Series, MCGRAW-HILL Higher Education.
Toro, E. (2001), Shock-Capturing Methods for Free-Surface Shallow Flows, Wiley, New York, p.
109.
USBR (1988), Downstream hazard classification guidelines, ACER Technical Memorandum 5,
United States Bureau of Reclamation (USBR).
USDA (1986), Urban hydrology for small watersheds, Environment agency report, United States
Department of Agriculture (USDA).
USDA (2004), National engineering handbook part 630 hydrology, Technical report, United States
Department of Agriculture (USDA).
Van Rijn, L. (1984a), ‘Sediment pick-up functions’, J. Hydr. Eng., ASCE. 110(10), 1494–1502.
Van Rijn, L. (1984b), ‘Sediment transport, part i: bed load transport’, J. Hydr. Eng., ASCE.
110(10), 1431–1456.
Van Rijn, L. (1984c), ‘Sediment transport, part ii: suspended load transport’, J. Hydr. Eng.,
ASCE. 110(11), 1613–1641.
Vázquez-Cendón, M. (1999), ‘Improved treatment of source terms in upwind schemes for the shal-
low water equations in channels with irregular geometry’, Journal of Computational Physics
148(2), 497–526.
229 12 References
Vreugdenhil, C. (1994), Numerical methods for shallow water flow, Water Science and Technology
Library. Kluwer Academic Publishers.
Wilcock, P. & Crowe, J. (2003), ‘Surface-based transport model for mixed-sized sediment’, Journal
of Hydraulic Engineering 129(2), 120–128.
Wong, M. (2003), ‘Does the bedload transport relation of meyer-peter and muller fits its own
data?’, Proc. 30th IAHR-Congress, Thessaloniki, Greece, 8 pp. .
Wong, M. & Parker, G. (2006), ‘Re-analysis and correction of bed load relation of meyer-peter and
mller using their own database’, Journal of Hydraulic Engineering 132(11), 1159–1168.
Yarnell, D. (1934a), ‘Bridge piers as channel obstructions’, U.S. Department of Agriculture, Tech.
Bull. 442.
Yarnell, D. (1934b), ‘Pile trestles as channel obstructions’, U.S. Department of Agriculture, Tech.
Bull. 429.
Zanke, U. (n.d.), Berechnung der sinkgeschwindigkeiten von sedimenten, Technical report, Tech-
nischen Universitt Hannover.
Zhang, R. (1961), River Dynamics (in Chinese), Industry Press.