SMS:SRH-2D: Difference between revisions
AlanZundel (talk | contribs) |
|||
(44 intermediate revisions by 4 users not shown) | |||
Line 1: | Line 1: | ||
{{SMS Infobox Model | | {{SMS Infobox Model | | ||
|name= SRH-2D | |name= SRH-2D | ||
|screenshot= File:SMS SRH-2D 1.PNG | |||
|model_type= Two-dimensional (2D) hydraulic, sediment, temperature, and vegetation model for river systems | |model_type= Two-dimensional (2D) hydraulic, sediment, temperature, and vegetation model for river systems | ||
|developer= | |developer= | ||
Yong Lai <br /> | Yong Lai <br /> | ||
Bureau of Reclamation | Bureau of Reclamation | ||
|web_site= [ | |web_site= [https://www.usbr.gov/tsc/techreferences/computer%20software/models/srh2d/index.html Bureau of Reclamation] | ||
|tutorials= | |tutorials= | ||
General Section | '''General Section''' | ||
* Data Visualization | * Data Visualization | ||
* Mesh Editing | * Mesh Editing | ||
* Observation | * Observation | ||
Models Section | '''Models Section''' | ||
* SRH-2D | * SRH-2D | ||
* Boundary Conditions | |||
* Post Processing | |||
* Simulations | |||
* FHWA Tutorials | |||
}} | }} | ||
Line 21: | Line 26: | ||
== Features == | == Features == | ||
SRH-2D is a 2D model | SRH-2D is a 2D model that is particularly useful for problems where 2D effects are important. Examples include flows with in-stream structures, through bends, with perched rivers, with multiple channel systems, and with complex floodplains. A 2D model may also be needed for local flow velocities, eddy patterns and flow recirculation, lateral variations, flow spills over banks and levees, and flow diversion and bifurcation. | ||
<!-- check statements as accuracy has been called into question by customers | |||
* SRH-2D solves the 2D depth-averaged form of the diffusive wave or the dynamic wave equations. The dynamic wave equations are the standard St. Venant depth-averaged shallow water equations. | * SRH-2D solves the 2D depth-averaged form of the diffusive wave or the dynamic wave equations. The dynamic wave equations are the standard St. Venant depth-averaged shallow water equations. | ||
* Both the diffusive wave and dynamic wave solvers use the implicit scheme to achieve solution robustness and efficiency. | * Both the diffusive wave and dynamic wave solvers use the implicit scheme to achieve solution robustness and efficiency.--> | ||
* Both steady or unsteady flows may be simulated. | * Both steady or unsteady flows may be simulated. | ||
* All flow regimes | * All flow regimes (i.e. subcritical, transcritical, and supercritical flows) may be simulated simultaneously without the need of a special treatment. | ||
* Solution domain may include a combination of main channels, side channels, floodplains, and overland flow areas. | * Solution domain may include a combination of main channels, side channels, floodplains, and overland flow areas. | ||
* Solved variables include water surface elevation, water depth | * Solved variables include water surface elevation, water depth and depth averaged velocity, magnitude, and vector. Output information includes above variables, plus Froude number, and bed shear stress. | ||
* SRH-2D can also be used to simulate sediment transport conditions, including erosion and deposition of cohesive and non-cohesive sediments. | * SRH-2D can also be used to simulate sediment transport conditions, including erosion and deposition of cohesive and non-cohesive sediments. Simulations can include bed load, suspended load, and total load using a variety of sediment transport and erosion equations. | ||
** Limitation: SRH-2D does not simulate bank erosion and channel migration processes other than what would be computed based only on the sediment transport relationships. | ** Limitation: SRH-2D does not simulate bank erosion and channel migration processes other than what would be computed based only on the sediment transport relationships. | ||
** Limitation: SRH-2D sediment transport does not work with | ** Limitation: SRH-2D sediment transport does not work with 1D hydraulic structures (e.g. culverts and 1D weir over topping), nor does the current version work with bridge pressure flow internal boundary conditions. | ||
* The additional solved variables for a sediment transport simulation include sediment transport by size fraction, bed elevation, bed elevation change (erosion and deposition), sediment concentration, changes to bed gradation, and D50 particle size. | * The additional solved variables for a sediment transport simulation include sediment transport by size fraction, bed elevation, bed elevation change (erosion and deposition), sediment concentration, changes to bed gradation, and D50 particle size. | ||
Line 36: | Line 42: | ||
==Graphical Interface== | ==Graphical Interface== | ||
SRH-2D uses a custom interface to specify boundary conditions, model parameters, model control and material parameters. The interface includes the following: | SRH-2D uses a custom interface to specify boundary conditions, model parameters, model control, and material parameters. The interface includes the following: | ||
* [[SMS:SRH-2D Simulation|SRH-2D Simulation]] | * [[SMS:SRH-2D Simulation|SRH-2D Simulation]] | ||
Line 42: | Line 48: | ||
* [[SMS:SRH-2D Model Control|SRH-2D Model Control]] | * [[SMS:SRH-2D Model Control|SRH-2D Model Control]] | ||
* [[SMS:SRH-2D Material Properties|SRH-2D Material Properties]] – SRH-2D is able to work with a number of material zones. Materials may be created from a Materials coverage, an SRH-2D coverage, or directly from the mesh. | * [[SMS:SRH-2D Material Properties|SRH-2D Material Properties]] – SRH-2D is able to work with a number of material zones. Materials may be created from a Materials coverage, an SRH-2D coverage, or directly from the mesh. | ||
* SRH-2D Sediment Material Properties – Similar to ''Material Properties'', a number of sediment material zones are required to describe the bed material thickness, particle sizes, and bulk densities for multiple sediment layers. | * [[SMS:SRH-2D Sediment Material Properties|SRH-2D Sediment Material Properties]] – Similar to ''Material Properties'', a number of sediment material zones are required to describe the bed material thickness, particle sizes, and bulk densities for multiple sediment layers. | ||
* [[SMS:SRH-2D Coverages|SRH-2D Coverages]] | * [[SMS:SRH-2D Coverages|SRH-2D Coverages]] | ||
* [[SMS:SRH-2D Files|SRH-2D input and output files]] – When wanting to execute an SRH-2D model, export the model native files using the '''Export SRH-2D Files''' or '''Save, Export, and Launch SRH-2D''' commands in the ''SRH-2D'' menu. The native files include: | * [[SMS:SRH-2D Files|SRH-2D input and output files]] – When wanting to execute an SRH-2D model, export the model native files using the '''Export SRH-2D Files''' or '''Save, Export, and Launch SRH-2D''' commands in the ''SRH-2D'' menu. The native files include: | ||
** [[SMS:SRH-2D Files#SRHHYDRO File|SRHHYDRO]] – Contains key information about the simulation while acting as a directory to other files for SRH-2D to use. | |||
* [[SMS:SRH-2D Files#SRHHYDRO File|SRHHYDRO]] – Contains key information about the simulation while acting as a directory to other files for SRH-2D to use. | ** [[SMS:SRH-2D Files#SRHGEOM File|SRHGEOM]] – Tells SRH-2D where each element is located and the characteristics of that element. | ||
* [[SMS:SRH-2D Files#SRHGEOM File|SRHGEOM]] – Tells SRH-2D where each element is located and the characteristics of that element. | ** [[SMS:SRH-2D Files#SRHMAT File|SRHMAT]] – Gives each element a material type. | ||
* [[SMS:SRH-2D Files#SRHMAT File|SRHMAT]] – Gives each element a material type. | ** SRHSEDMAT – Gives each element a sediment material type. | ||
* SRHSEDMAT – Gives each element a sediment material type. | ** [[SMS:SRH-2D Files#SRHMONITORPTS File|SRHMPOINTS]] – Tells SRH-2D that there are monitor points to watch and where those points are located. | ||
* [[SMS:SRH-2D Files#SRHMONITORPTS File|SRHMPOINTS]] – Tells SRH-2D that there are monitor points to watch and where those points are located. | |||
In the past, this model has been utilized through the [[SMS:Generic_Model_Graphical_Interface|Generic Model | In the past, this model has been utilized through the [[SMS:Generic_Model_Graphical_Interface|Generic Model graphical interface]]. The [http://www.usbr.gov/tsc/sediment/model/srh2d/Downloads/SRH2D-v2-Distribution-Package-Jan2009.zip SRH-2D version 2.0 Distribution] included SRH-2D template files for both SMS 8.0 and SMS 10.0. These are no longer needed with the custom interface. | ||
==Steps to Create an SRH-2D Model== | ==Steps to Create an SRH-2D Model== | ||
[[File:SRH-2D model dev.png|thumb|500 px|SRH-2D Model Development]] | |||
To create an SRH-2D model, the following general steps should be followed: | To create an SRH-2D model, the following general steps should be followed: | ||
# Gather data pertinent to the project and location. This should include bathymetry data, roughness data (Manning's n value), coordinate system corresponding to the data, and flow data. | # Gather data pertinent to the project and location. This should include bathymetry data, roughness data (Manning's n value), coordinate system corresponding to the data, and flow data. | ||
Line 61: | Line 66: | ||
# Add bathymetry data. This may come as survey data, Lidar data, or Raster DEM data to name a few. | # Add bathymetry data. This may come as survey data, Lidar data, or Raster DEM data to name a few. | ||
# Check the triangulation or raw data display. It is important to make sure that SMS is reading the data the same way that it was measured. Turning on contours will allow the user to view what SMS sees and make adjustments as needed. Contours may be turned on using the '''Display Options''' command. Optionally, the user may use the tools available in SMS for refinement of the data. | # Check the triangulation or raw data display. It is important to make sure that SMS is reading the data the same way that it was measured. Turning on contours will allow the user to view what SMS sees and make adjustments as needed. Contours may be turned on using the '''Display Options''' command. Optionally, the user may use the tools available in SMS for refinement of the data. | ||
# [[SMS:SRH-2D Coverages|Create coverages]]. A simple SRH-2D project would likely include a | # [[SMS:SRH-2D Coverages|Create coverages]]. A simple SRH-2D project would likely include a Mesh Generator coverage which contains mesh type and bathymetry data, an SRH-2D Boundary Conditions coverage which holds boundary arcs and flow data, an SRH-2D Materials coverage which maps material types defined for a region to each cell/element defined in the mesh, and an SRH-2D Monitor coverage which specifies locations where results are gathered. Coverages may be created by right-clicking on the map data folder from the data tree and selecting '''New Coverage''' or by duplicating an existing coverage. The coverage is assigned a type upon creation which can be changed at any time. Changing a coverage type may necessitate other modifications to a simulation. For an SRH-2D model, select the SRH-2D coverages as they relate to the data that will correspond to that type. | ||
#* For sediment transport, | #* For sediment transport, an SRH-2D Sediment Materials coverage will be necessary. | ||
#* For sediment transport, set ''BC Type'' parameters for the | #* For sediment transport, set ''BC Type'' parameters for the Boundary Conditions coverage. | ||
# Outline the workspace with arcs. Here the user is defining regions of the model location that will have unique features. For example, locations of more water interaction will need more detail which equates to more nodes; locations with different roughness values will need to be separated for material type assignments. Create polygons for areas of similar characteristics. Keep in mind that SMS has a variety of tools available to adjust the arcs to meet the modeling needs of the project. | # Outline the workspace with arcs. Here the user is defining regions of the model location that will have unique features. For example, locations of more water interaction will need more detail which equates to more nodes; locations with different roughness values will need to be separated for material type assignments. Create polygons for areas of similar characteristics. Keep in mind that SMS has a variety of tools available to adjust the arcs to meet the modeling needs of the project. | ||
# Build polygons in the | # Build polygons in the Mesh Generator, Materials, and Sediment Transport coverages. | ||
# Assign attributes to the polygons. Direct SMS to what materials and what mesh type is to be built over that polygon. Specify how SMS should assign elevation data by selecting the bathymetry source. | # Assign attributes to the polygons. Direct SMS to what materials and what mesh type is to be built over that polygon. Specify how SMS should assign elevation data by selecting the bathymetry source. | ||
#* Additional polygon attributes must be set for sediment transport in the | #* Additional polygon attributes must be set for sediment transport in the Sediment Materials coverage. | ||
# In the | # In the Boundary Conditions coverage, assign attributes to the arcs by giving the arcs [[SMS:SRH-2D Coverages#Boundary Conditions Coverage|boundary conditions]]. For interior arcs, the only option is a monitor line. For exterior arcs the user may choose from a variety of inlet conditions, whether subcritical or supercritical, as well as outlet or water surface elevation options. For no-flow boundaries, the options of a wall or symmetry are available. | ||
# Prepare to build the mesh. Review the information given to SMS to ensure that the field data matches what is represented virtually for the region. After a review of the inputs for the model, the mesh is ready to be built. The mesh is built by converting the Map coverage to a 2D Mesh. | # Prepare to build the mesh. Review the information given to SMS to ensure that the field data matches what is represented virtually for the region. After a review of the inputs for the model, the mesh is ready to be built. The mesh is built by converting the Map coverage to a 2D Mesh. | ||
# Build the mesh from the | # Build the mesh from the Mesh Generator coverage. If bathymetry data was assigned to the polygons on the Mesh Generator coverage, the mesh will contain the correct data for the project. Otherwise, the bathymetry data must be interpolated to the mesh. | ||
# Create an SRH-2D simulation. | # Create an SRH-2D simulation. | ||
# Add components to the simulation. Components usually include all coverages and the 2D mesh. | # Add components to the simulation. Components usually include all coverages and the 2D mesh. | ||
Line 78: | Line 83: | ||
==Sediment Transport== | ==Sediment Transport== | ||
Sediment transport modeling in SRH-2D follows the same steps outlined above. | Sediment transport modeling in SRH-2D follows the same steps outlined above. However, sediment transport modeling requires the following additional data: | ||
* Bed surface and subsurface samples to determine particle gradation which could vary through the model domain. | * Bed surface and subsurface samples to determine particle gradation which could vary through the model domain. | ||
* Sediment transport data of suspended and bed loads over a range of flows. | * Sediment transport data of suspended and bed loads over a range of flows. | ||
* Samples of floodplain materials. | * Samples of floodplain materials. | ||
* Results of erodibility testing | * Results of erodibility testing on floodplain materials. | ||
To perform a sediment transport simulation, follow these general steps: | To perform a sediment transport simulation, follow these general steps: | ||
#Set up and | #Set up and run an SRH-2D hydraulic model. The following should be kept in mind: | ||
## Due to the computational expense of sediment transport modeling, if the ultimate goal of a modeling project is to conduct a sediment transport study the number of elements should be limited to less than 40,000. In fact, we recommend that the number of elements be kept to less than 30,000. When converting an existing hydraulic model into a sediment transport model the same limits apply meaning the mesh may require coarsening. | ## Due to the computational expense of sediment transport modeling, if the ultimate goal of a modeling project is to conduct a sediment transport study, the number of elements should be limited to less than 40,000. In fact, we recommend that the number of elements be kept to less than 30,000. When converting an existing hydraulic model into a sediment transport model the same limits apply, meaning the mesh may require coarsening. This can be accomplished by changing the parameters of the Mesh Generation coverage and generating a coarser mesh. Hydraulic results of the coarse mesh should be reviewed to ensure they still retain comparable and adequate results. | ||
##A sediment transport simulation must have a representative condition at the beginning of the simulation to avoid unrealistic sediment transport, erosion, and deposition as the model transitions from an initial condition to a convergence condition. | ##A sediment transport simulation must have a representative condition at the beginning of the simulation to avoid unrealistic sediment transport, erosion, and deposition as the model transitions from an initial condition to a convergence condition. | ||
#Specify sediment transport properties in a [[SMS:SRH-2D_Coverages#BC_Types_Dialog| | #Specify sediment transport properties in a [[SMS:SRH-2D_Coverages#BC_Types_Dialog|Boundary Condition coverage]]. Typically this would be created by copying the Boundary Condition coverage and changing the ''Run Type'' of the new coverage to "Mobile". This is accessed by right-clicking on the coverage and selecting the ''BC Types'' dialog. This allows for the specification of the sediment parameters. | ||
#Specify sediment data on the boundary conditions which include inflow. This includes Inlet_Q and Inlet_SC boundary types. | #Specify sediment data on the boundary conditions which include inflow. This includes Inlet_Q and Inlet_SC boundary types. | ||
#Additional monitor lines (in the | #Additional monitor lines (in the Boundary Condition coverage) and monitor points (in the Monitor coverage) may be desired in the sediment transport simulation to monitor the sediment transport rates. In a sediment transport run two output files are create for each point. One focuses on detailed hydraulic results and the other contains detailed sediment transport results. In the "quickwin" version of SRH-2D, bed elevation plots are shown during the sediment transport simulation for the first two monitoring points (ID 1 and ID 2). | ||
#A sediment transport simulation also requires the creation of a [[SMS:SRH-2D_Coverages#Sediment_Materials_Coverage| | #A sediment transport simulation also requires the creation of a [[SMS:SRH-2D_Coverages#Sediment_Materials_Coverage|Sediment Materials coverage]]. Sediment transport material properties include sediment gradations, bulk densities, and sediment layer thicknesses for each sediment material type. SRH-2D guidance indicates that a minimum of two layers is recommended even if the bed material is vertically uniform. The top layer interacts with the sediment in the water column and is part of the active layer used for computing sediment transport. A zero-thickness layer can also be specified to set the area as non-erodible.” | ||
#At the simulation level, the model control parameters for a sediment transport simulation are identical to those of a hydraulic simulation including run times, timestep, restart file location, etc. Note that sediment transport simulations take much longer than hydraulic simulations. A 12-hour sediment transport simulation can take several hours of computer time, so simulation times are often limited to hours or days. A one-year sediment transport simulation can take multiple weeks of computer time. | #At the simulation level, the model control parameters for a sediment transport simulation are identical to those of a hydraulic simulation including run times, timestep, restart file location, etc. Note that sediment transport simulations take much longer than hydraulic simulations. A 12-hour sediment transport simulation can take several hours of computer time, so simulation times are often limited to hours or days. A one-year sediment transport simulation can take multiple weeks of computer time. | ||
Line 99: | Line 104: | ||
Some common errors in running SRH-2D include: | Some common errors in running SRH-2D include: | ||
* The [[SMS:Mesh Quality|quality]] of the 2D mesh and the placement of the boundary conditions arcs will significantly impact how SRH-2D processes the project. Errors frequently occur when a mesh is poorly generated (uneven patches, narrow triangles, extreme size differences between elements, etc.) or boundary arcs are misplaced. | |||
* When importing SRH-2D projects created in earlier versions of SMS, the project may be slow to open in a current version of SMS. This is common with SRH-2D projects created in SMS 11.2. After the project has been opened in a newer version of SMS, it needs to be saved to prevent the project from having a slow open. | |||
* Running SRH-2D requires both the pre-SRH-2D executable and the SRH-2D executable. If SMS cannot find these executables, the model run will encounter an error. The [[SMS:Preferences#File Locations|''File Locations'']] tab in the [[SMS:Preferences|''Preferences'']] dialog can be used to specify the path to these executables if SMS cannot find them. | * Running SRH-2D requires both the pre-SRH-2D executable and the SRH-2D executable. If SMS cannot find these executables, the model run will encounter an error. The [[SMS:Preferences#File Locations|''File Locations'']] tab in the [[SMS:Preferences|''Preferences'']] dialog can be used to specify the path to these executables if SMS cannot find them. | ||
* While SRH-2D can handle a large variety of geometries, using a 2D mesh with a large number of elements, poor transition between elements, and other mesh quality issues can cause SRH-2D to become unstable. It is advised that the mesh quality be reviewed before running SRH-2D. When using sediment transport options in a simulation, it is advised that the 2D mesh contain fewer than 50,000 elements. | |||
*Starting in SMS version 12.3.4, SRH-2D solutions are cell-centered. A post-processor was added to interpolate the solution to the mesh nodes after running SRH-2D. | |||
*Starting in SMS version 13.0, monitor lines are no longer set on the Boundary Conditions coverage. Because of this change, when importing SRH-2D projects create in earlier version several items will happen: | |||
**SMS will separate the monitoring lines from other arcs in the boundary conditions coverage. The monitoring lines will be placed in their own Monitor coverage. | |||
**The former Boundary Condition coverage will be made into two coverages, one for the monitor lines and the other for the existing boundary conditions. The names of these coverages will have either "Boundary Conditions" or "Monitor" attached to designate which coverage contains the boundary condition arcs and which contains the monitor lines. | |||
**If a Monitor Points coverage existed in the project, this coverage will be changed to a Monitor coverage and the coverage name will have "Monitor" attached. | |||
**If the imported project had both monitor points and monitor lines, their will now be two Monitor coverages in the project. | |||
**By default, the Monitor coverage with the monitor lines will be included in any simulations. | |||
**The monitor points and monitor lines from the previous project will need to be merged into one Monitor coverage and linked to the simulation in order for both points and lines to be included in the simulation run. | |||
**The default boundary condition for arcs on an SRH-2D boundary condition coverage is now the "Wall" type. The boundary condition arcs should be reviewed to ensure no parameters were lost in the conversion. | |||
* Note that iTIME is in seconds. If the restart times are in hours, then iTIME values will half the number of seconds for those hours. The TSO file will identify the times corresponding to the restart files based on their ID number. If not needed it can be safely deleted. | |||
* There is sometimes a disparity between the depth+elevation and the water surface elevation because of how PostSRH-2D processes the data. SRH-2D calculates the values at the centroid, and then PostSRH-2D goes back through and averages to the nodes, which can cause the depth to appear a little bumpy. This generally does not cause an issue for most projects. | |||
See [[SMS:SRH-2D Errors|SRH-2D Error Messages]] for more information. | See [[SMS:SRH-2D Errors|SRH-2D Error Messages]] for more information. | ||
Line 141: | Line 159: | ||
==SRH-2D Files== | ==SRH-2D Files== | ||
See the article [[SMS:SRH-2D Files|SRH-2D Files]] | See the article [[SMS:SRH-2D Files|SRH-2D Files]] | ||
== External Links – SRH-2D Version 3.2 == | |||
SRH-2D version 3 uses the same user manual as SRH-2D version 2 with some additional papers on sediment transport. | |||
* [http://smsdocs.aquaveo.com/SRH2D-v3-Readme.pdf SRH-2D version 3.2.0 Read Me] | |||
* [https://www.usbr.gov/tsc/techreferences/computer%20software/models/srh2d/downloads/Version_3.2.0_2017_0306.zip SRH-2D version 3.2 Distribution] – Includes software, manual, template file for integration into SMS interface, and tutorials | |||
* [http://smsdocs.aquaveo.com/Final_Report_FHWA_2016.pdf FWHA 2016 report] | |||
* [http://smsdocs.aquaveo.com/Giude_on_Unsteady_Flow_Modeling_with_SRH-2D.pdf Guide on Unsteady Flow Modeling with SRH-2D] | |||
== External Links – SRH-2D Version 2.0 == | == External Links – SRH-2D Version 2.0 == | ||
Line 148: | Line 174: | ||
=== Papers / Presentations === | === Papers / Presentations === | ||
* [ | * [https://www.usbr.gov/tsc/techreferences/computer%20software/models/srh2d/downloads/Theory_Paper-SRH2D-v2-2009.pdf SRH-2D Theory Paper] | ||
* [ | * [https://www.usbr.gov/tsc/techreferences/computer%20software/models/srh2d/downloads/SRH2D-v2-Training-Dec2007-Taiwan.pdf SRH-2D Training Presentation] | ||
* [http://smsdocs.aquaveo.com/Session_9A-1_Bountry.pdf 2006 FISC Paper on Savage Rapids Dam Removal Project "Comparison of Numerical Hydraulic Models Applied To The Removal of Savage Rapids Dam Near Grants Pass, Oregon"] | * [http://smsdocs.aquaveo.com/Session_9A-1_Bountry.pdf 2006 FISC Paper on Savage Rapids Dam Removal Project "Comparison of Numerical Hydraulic Models Applied To The Removal of Savage Rapids Dam Near Grants Pass, Oregon"] | ||
* [ | * [https://shira.lawr.ucdavis.edu/publications List of journal articles using SMS by Prof. Greg Pasternack, UC Davis] | ||
* [https://www.mdpi.com/2073-4441/12/4/1032/htm#cite Lai, Y. G. (2020). A Two-Dimensional Depth-Averaged Sediment Transport Mobile-Bed Model with Polygonal Meshes. Water, 12(4), 1032.] | |||
=== Project Reports === | === Project Reports === | ||
* [http://www.usbr.gov/tsc/sediment/model/srh2d/Downloads/Final%20Report%20Savage%20Rapids%20Jan%2023%202006.pdf Bountry J.A. and Lai, Y.G. (2006)."Numerical modeling of flow hydraulics in support of the Savage Rapids Dam removal."] | <!-- broken link * [http://www.usbr.gov/tsc/sediment/model/srh2d/Downloads/Final%20Report%20Savage%20Rapids%20Jan%2023%202006.pdf Bountry J.A. and Lai, Y.G. (2006)."Numerical modeling of flow hydraulics in support of the Savage Rapids Dam removal."]--> | ||
* [http://smsdocs.aquaveo.com/SandyRiverDamRemovalSedimentModel.pdf Lai, Y.G., Holburn, E.R., and Bauer, T.R. (2006)."Analysis of sediment transport following removal of the Sandy River Delta Dam."] | * [http://smsdocs.aquaveo.com/SandyRiverDamRemovalSedimentModel.pdf Lai, Y.G., Holburn, E.R., and Bauer, T.R. (2006)."Analysis of sediment transport following removal of the Sandy River Delta Dam."] | ||
* [http://www.usbr.gov/tsc/sediment/model/srh2d/Downloads/Final%20Report-Elwha%20Nov%2014%202006.pdf Lai, Y.G. and Bountry, J.A. (2006). "Numerical hydraulic modeling and assessment in support of Elwha Surface Diversion Project."] | <!-- broken links * [http://www.usbr.gov/tsc/sediment/model/srh2d/Downloads/Final%20Report-Elwha%20Nov%2014%202006.pdf Lai, Y.G. and Bountry, J.A. (2006). "Numerical hydraulic modeling and assessment in support of Elwha Surface Diversion Project."] | ||
* [http://www.usbr.gov/tsc/sediment/model/srh2d/Downloads/Final%20Report%20Dungeness%202dmodel%20report%20May7_07.pdf Lai, Y.G. and Bountry, J.A. (2007). "Numerical modeling study of levee setback alternatives for lower Dungeness River, Washington"] | * [http://www.usbr.gov/tsc/sediment/model/srh2d/Downloads/Final%20Report%20Dungeness%202dmodel%20report%20May7_07.pdf Lai, Y.G. and Bountry, J.A. (2007). "Numerical modeling study of levee setback alternatives for lower Dungeness River, Washington"]--> | ||
* [http://smsdocs.aquaveo.com.s3.amazonaws.com/srh_model_applications_and_progress_report_on_bank_erosion_and_turbidity_current_models.pdf Lai, Y.G. and Greimann, B.P. (2011). "SRH Model Applications and Progress Report on Bank Erosion and Turbidity Current Models"] | * [http://smsdocs.aquaveo.com.s3.amazonaws.com/srh_model_applications_and_progress_report_on_bank_erosion_and_turbidity_current_models.pdf Lai, Y.G. and Greimann, B.P. (2011). "SRH Model Applications and Progress Report on Bank Erosion and Turbidity Current Models"] | ||
Line 169: | Line 196: | ||
=== In the News === | === In the News === | ||
* [ | * [https://www.wired.com/2007/07/dam/ News article "Computer Modeling Smoothes a Dam Hard Job"] | ||
* [ | * [https://www.wired.com/2007/08/dam-gallery/ Photos related to Wired.com article "Computer Modeling Smoothes a Dam Hard Job"] | ||
== External Links – SRH-W == | == External Links – SRH-W == | ||
=== SRH-W Version 1.1 === | === SRH-W Version 1.1 === | ||
* [ | * [https://www.usbr.gov/tsc/techreferences/computer%20software/models/srh2d/downloads/SRH-W%20v1.1%20User%20Manual%20June2007.pdf SRH-W version 1.1 User Manual] | ||
* [ | <!--vlink broken * [https://www.usbr.gov/tsc/techreferences/computer%20software/models/srh2d/SRH-W1.1-Distribution-Package-2007.zip SRH-W version 1.1 Distribution Package] – Includes Software, Manual, and Tutorials --> | ||
=== Papers / Presentations === | === Papers / Presentations === |
Latest revision as of 22:07, 1 May 2023
SRH-2D | |
---|---|
SRH-2D Screenshot | |
Model Info | |
Model type | Two-dimensional (2D) hydraulic, sediment, temperature, and vegetation model for river systems |
Developer |
Yong Lai |
Web site | Bureau of Reclamation |
Tutorials |
General Section
Models Section
|
SRH-2D, or the Sedimentation and River Hydraulics – Two-Dimensional model, is a two-dimensional (2D) hydraulic, sediment, temperature, and vegetation model for river systems under development at the Bureau of Reclamation. It was known as SRH-W in version 1, and the name was changed to the current SRH-2D from version 2 onward.
The SRH-2D model is included in the Community Version of SMS. however, it is limited to a single simulation. For multiple simulation management, full post processing capabilities, and a complete suite of tools, the interface can be added to a paid edition of SMS.
Features
SRH-2D is a 2D model that is particularly useful for problems where 2D effects are important. Examples include flows with in-stream structures, through bends, with perched rivers, with multiple channel systems, and with complex floodplains. A 2D model may also be needed for local flow velocities, eddy patterns and flow recirculation, lateral variations, flow spills over banks and levees, and flow diversion and bifurcation.
- Both steady or unsteady flows may be simulated.
- All flow regimes (i.e. subcritical, transcritical, and supercritical flows) may be simulated simultaneously without the need of a special treatment.
- Solution domain may include a combination of main channels, side channels, floodplains, and overland flow areas.
- Solved variables include water surface elevation, water depth and depth averaged velocity, magnitude, and vector. Output information includes above variables, plus Froude number, and bed shear stress.
- SRH-2D can also be used to simulate sediment transport conditions, including erosion and deposition of cohesive and non-cohesive sediments. Simulations can include bed load, suspended load, and total load using a variety of sediment transport and erosion equations.
- Limitation: SRH-2D does not simulate bank erosion and channel migration processes other than what would be computed based only on the sediment transport relationships.
- Limitation: SRH-2D sediment transport does not work with 1D hydraulic structures (e.g. culverts and 1D weir over topping), nor does the current version work with bridge pressure flow internal boundary conditions.
- The additional solved variables for a sediment transport simulation include sediment transport by size fraction, bed elevation, bed elevation change (erosion and deposition), sediment concentration, changes to bed gradation, and D50 particle size.
The Bureau of Reclamation does not provide technical support for SRH-2D.
Graphical Interface
SRH-2D uses a custom interface to specify boundary conditions, model parameters, model control, and material parameters. The interface includes the following:
- SRH-2D Simulation
- SRH-2D Menu
- SRH-2D Model Control
- SRH-2D Material Properties – SRH-2D is able to work with a number of material zones. Materials may be created from a Materials coverage, an SRH-2D coverage, or directly from the mesh.
- SRH-2D Sediment Material Properties – Similar to Material Properties, a number of sediment material zones are required to describe the bed material thickness, particle sizes, and bulk densities for multiple sediment layers.
- SRH-2D Coverages
- SRH-2D input and output files – When wanting to execute an SRH-2D model, export the model native files using the Export SRH-2D Files or Save, Export, and Launch SRH-2D commands in the SRH-2D menu. The native files include:
- SRHHYDRO – Contains key information about the simulation while acting as a directory to other files for SRH-2D to use.
- SRHGEOM – Tells SRH-2D where each element is located and the characteristics of that element.
- SRHMAT – Gives each element a material type.
- SRHSEDMAT – Gives each element a sediment material type.
- SRHMPOINTS – Tells SRH-2D that there are monitor points to watch and where those points are located.
In the past, this model has been utilized through the Generic Model graphical interface. The SRH-2D version 2.0 Distribution included SRH-2D template files for both SMS 8.0 and SMS 10.0. These are no longer needed with the custom interface.
Steps to Create an SRH-2D Model
To create an SRH-2D model, the following general steps should be followed:
- Gather data pertinent to the project and location. This should include bathymetry data, roughness data (Manning's n value), coordinate system corresponding to the data, and flow data.
- Specify a coordinate system. This is done in the Display Projections dialog accessed through the Display menu.
- Add bathymetry data. This may come as survey data, Lidar data, or Raster DEM data to name a few.
- Check the triangulation or raw data display. It is important to make sure that SMS is reading the data the same way that it was measured. Turning on contours will allow the user to view what SMS sees and make adjustments as needed. Contours may be turned on using the Display Options command. Optionally, the user may use the tools available in SMS for refinement of the data.
- Create coverages. A simple SRH-2D project would likely include a Mesh Generator coverage which contains mesh type and bathymetry data, an SRH-2D Boundary Conditions coverage which holds boundary arcs and flow data, an SRH-2D Materials coverage which maps material types defined for a region to each cell/element defined in the mesh, and an SRH-2D Monitor coverage which specifies locations where results are gathered. Coverages may be created by right-clicking on the map data folder from the data tree and selecting New Coverage or by duplicating an existing coverage. The coverage is assigned a type upon creation which can be changed at any time. Changing a coverage type may necessitate other modifications to a simulation. For an SRH-2D model, select the SRH-2D coverages as they relate to the data that will correspond to that type.
- For sediment transport, an SRH-2D Sediment Materials coverage will be necessary.
- For sediment transport, set BC Type parameters for the Boundary Conditions coverage.
- Outline the workspace with arcs. Here the user is defining regions of the model location that will have unique features. For example, locations of more water interaction will need more detail which equates to more nodes; locations with different roughness values will need to be separated for material type assignments. Create polygons for areas of similar characteristics. Keep in mind that SMS has a variety of tools available to adjust the arcs to meet the modeling needs of the project.
- Build polygons in the Mesh Generator, Materials, and Sediment Transport coverages.
- Assign attributes to the polygons. Direct SMS to what materials and what mesh type is to be built over that polygon. Specify how SMS should assign elevation data by selecting the bathymetry source.
- Additional polygon attributes must be set for sediment transport in the Sediment Materials coverage.
- In the Boundary Conditions coverage, assign attributes to the arcs by giving the arcs boundary conditions. For interior arcs, the only option is a monitor line. For exterior arcs the user may choose from a variety of inlet conditions, whether subcritical or supercritical, as well as outlet or water surface elevation options. For no-flow boundaries, the options of a wall or symmetry are available.
- Prepare to build the mesh. Review the information given to SMS to ensure that the field data matches what is represented virtually for the region. After a review of the inputs for the model, the mesh is ready to be built. The mesh is built by converting the Map coverage to a 2D Mesh.
- Build the mesh from the Mesh Generator coverage. If bathymetry data was assigned to the polygons on the Mesh Generator coverage, the mesh will contain the correct data for the project. Otherwise, the bathymetry data must be interpolated to the mesh.
- Create an SRH-2D simulation.
- Add components to the simulation. Components usually include all coverages and the 2D mesh.
- Set the model run parameters in the SRH-2D Model Control dialog.
- Run the model. First, be sure that the project is saved. Next, right-click on the SRH-2D simulation and select Export SRH-2D. Finally, right-click on the simulation and select Launch SRH-2D.
- Analysis of results and post processing.
Sediment Transport
Sediment transport modeling in SRH-2D follows the same steps outlined above. However, sediment transport modeling requires the following additional data:
- Bed surface and subsurface samples to determine particle gradation which could vary through the model domain.
- Sediment transport data of suspended and bed loads over a range of flows.
- Samples of floodplain materials.
- Results of erodibility testing on floodplain materials.
To perform a sediment transport simulation, follow these general steps:
- Set up and run an SRH-2D hydraulic model. The following should be kept in mind:
- Due to the computational expense of sediment transport modeling, if the ultimate goal of a modeling project is to conduct a sediment transport study, the number of elements should be limited to less than 40,000. In fact, we recommend that the number of elements be kept to less than 30,000. When converting an existing hydraulic model into a sediment transport model the same limits apply, meaning the mesh may require coarsening. This can be accomplished by changing the parameters of the Mesh Generation coverage and generating a coarser mesh. Hydraulic results of the coarse mesh should be reviewed to ensure they still retain comparable and adequate results.
- A sediment transport simulation must have a representative condition at the beginning of the simulation to avoid unrealistic sediment transport, erosion, and deposition as the model transitions from an initial condition to a convergence condition.
- Specify sediment transport properties in a Boundary Condition coverage. Typically this would be created by copying the Boundary Condition coverage and changing the Run Type of the new coverage to "Mobile". This is accessed by right-clicking on the coverage and selecting the BC Types dialog. This allows for the specification of the sediment parameters.
- Specify sediment data on the boundary conditions which include inflow. This includes Inlet_Q and Inlet_SC boundary types.
- Additional monitor lines (in the Boundary Condition coverage) and monitor points (in the Monitor coverage) may be desired in the sediment transport simulation to monitor the sediment transport rates. In a sediment transport run two output files are create for each point. One focuses on detailed hydraulic results and the other contains detailed sediment transport results. In the "quickwin" version of SRH-2D, bed elevation plots are shown during the sediment transport simulation for the first two monitoring points (ID 1 and ID 2).
- A sediment transport simulation also requires the creation of a Sediment Materials coverage. Sediment transport material properties include sediment gradations, bulk densities, and sediment layer thicknesses for each sediment material type. SRH-2D guidance indicates that a minimum of two layers is recommended even if the bed material is vertically uniform. The top layer interacts with the sediment in the water column and is part of the active layer used for computing sediment transport. A zero-thickness layer can also be specified to set the area as non-erodible.”
- At the simulation level, the model control parameters for a sediment transport simulation are identical to those of a hydraulic simulation including run times, timestep, restart file location, etc. Note that sediment transport simulations take much longer than hydraulic simulations. A 12-hour sediment transport simulation can take several hours of computer time, so simulation times are often limited to hours or days. A one-year sediment transport simulation can take multiple weeks of computer time.
Troubleshooting
Consult the SRH-2D Manual for troubleshooting problems in executing the model.
Some common errors in running SRH-2D include:
- The quality of the 2D mesh and the placement of the boundary conditions arcs will significantly impact how SRH-2D processes the project. Errors frequently occur when a mesh is poorly generated (uneven patches, narrow triangles, extreme size differences between elements, etc.) or boundary arcs are misplaced.
- When importing SRH-2D projects created in earlier versions of SMS, the project may be slow to open in a current version of SMS. This is common with SRH-2D projects created in SMS 11.2. After the project has been opened in a newer version of SMS, it needs to be saved to prevent the project from having a slow open.
- Running SRH-2D requires both the pre-SRH-2D executable and the SRH-2D executable. If SMS cannot find these executables, the model run will encounter an error. The File Locations tab in the Preferences dialog can be used to specify the path to these executables if SMS cannot find them.
- While SRH-2D can handle a large variety of geometries, using a 2D mesh with a large number of elements, poor transition between elements, and other mesh quality issues can cause SRH-2D to become unstable. It is advised that the mesh quality be reviewed before running SRH-2D. When using sediment transport options in a simulation, it is advised that the 2D mesh contain fewer than 50,000 elements.
- Starting in SMS version 12.3.4, SRH-2D solutions are cell-centered. A post-processor was added to interpolate the solution to the mesh nodes after running SRH-2D.
- Starting in SMS version 13.0, monitor lines are no longer set on the Boundary Conditions coverage. Because of this change, when importing SRH-2D projects create in earlier version several items will happen:
- SMS will separate the monitoring lines from other arcs in the boundary conditions coverage. The monitoring lines will be placed in their own Monitor coverage.
- The former Boundary Condition coverage will be made into two coverages, one for the monitor lines and the other for the existing boundary conditions. The names of these coverages will have either "Boundary Conditions" or "Monitor" attached to designate which coverage contains the boundary condition arcs and which contains the monitor lines.
- If a Monitor Points coverage existed in the project, this coverage will be changed to a Monitor coverage and the coverage name will have "Monitor" attached.
- If the imported project had both monitor points and monitor lines, their will now be two Monitor coverages in the project.
- By default, the Monitor coverage with the monitor lines will be included in any simulations.
- The monitor points and monitor lines from the previous project will need to be merged into one Monitor coverage and linked to the simulation in order for both points and lines to be included in the simulation run.
- The default boundary condition for arcs on an SRH-2D boundary condition coverage is now the "Wall" type. The boundary condition arcs should be reviewed to ensure no parameters were lost in the conversion.
- Note that iTIME is in seconds. If the restart times are in hours, then iTIME values will half the number of seconds for those hours. The TSO file will identify the times corresponding to the restart files based on their ID number. If not needed it can be safely deleted.
- There is sometimes a disparity between the depth+elevation and the water surface elevation because of how PostSRH-2D processes the data. SRH-2D calculates the values at the centroid, and then PostSRH-2D goes back through and averages to the nodes, which can cause the depth to appear a little bumpy. This generally does not cause an issue for most projects.
See SRH-2D Error Messages for more information.
Releases
Version 1
SRH-W, or Sedimentation and River Hydraulics – Watershed, is a two-dimensional (2D) hydraulic model for river systems and watersheds developed at the Bureau of Reclamation. SRH-W was originally developed for Reclamation internal use for various projects, and version 1.1 was released for public use.
SRH-W v1.1 is used for hydraulic flow simulation in rivers and runoff from watersheds, but without the sediment capability. It solves the 2D dynamic wave equations (the standard depth-averaged St. Venant equations) that are mainly used for river simulation. In addition, the diffusive wave solver is used for watershed runoff simulation and river simulation.
Version 1.1 is comparable to many existing models such as RMA-2 (US Army Corps of Engineers, 1996) and MIKE 21 (DHI software, 1996) in its river simulation capability. For watershed applications, SRH-W v1.1 is a distributed model for event based runoff simulation and has capabilities similar to CASC2D (Julien, et al, 1995).
Version 2
In Version 2, SRH-W was renamed to SRH-2D. This is the version that is associated with the users document that is distributed by the United States Bureau of Reclamation.
Version 2 solves the 2D dynamic wave equations, i.e., the depth-averaged St. Venant equations. Its modeling capability is comparable to some existing 2D models but SRH-2D claims a few boasting features. First, SRH-2D uses a flexible mesh that may contain arbitrarily shaped cells. In practice, the hybrid mesh of quadrilateral and triangular cells is recommended though purely quadrilateral or triangular elements may be used. A hybrid mesh may achieve the best compromise between solution accuracy and computing demand. Second, SRH-2D adopts very robust and stable numerical schemes with a seamless wetting-drying algorithm. The resultant outcome is that few tuning parameters are needed to obtain the final solution. SRH-2D was evolved from SRH-W which had the additional capability of watershed runoff modeling. Many features are improved from SRH-W.
Major Features of Version 2
Major SRH-2D capabilities are listed below
- 2D depth-averaged dynamic wave equations (the standard St. Venant equations) are solved with the finite-volume numerical method
- Steady state (with constant discharge) or unsteady flows (with flow hydrograph) may be simulated
- An implicit scheme is used for time integration to achieve solution robustness and efficiency
- An unstructured arbitrarily-shaped mesh is used which includes the structured quadrilateral mesh, the purely triangular mesh, or a combination of the two. Cartesian or raster mesh may also be used. In most applications, a combination of quadrilateral and triangular meshes is the best in terms of efficiency and accuracy
- All flow regimes, i.e., subcritical, transcritical, and supercritical flows, may be simulated simultaneously without the need for special treatments
- Robust and seamless wetting-drying algorithm; an
- Solved variables include water surface elevation, water depth, and depth averaged velocity. Output variables include the above, plus Froude number, and bed shear stress.
SRH-2D is a 2D model, and it is particularly useful for problems where 2D effects are important. Examples include flows with in-stream structures, through bends, with perched rivers, with side channel and agricultural returns, and with braided channel systems. A 2D model may also be needed if one is interested in local flow velocities, eddy patterns, flow recirculation, lateral velocity variation, and flow over banks and levees.
Version 3
SRH-2D version 3 is essentially Version 2 with Sediment Transport capability added. This version is currently distributed with the SMS package. Additional solved variables include sediment concentration, erosion and deposition, bed elevation, sediment transport rates, bed material D50 size, and bed material gradations.
SRH-2D Files
See the article SRH-2D Files
External Links – SRH-2D Version 3.2
SRH-2D version 3 uses the same user manual as SRH-2D version 2 with some additional papers on sediment transport.
- SRH-2D version 3.2.0 Read Me
- SRH-2D version 3.2 Distribution – Includes software, manual, template file for integration into SMS interface, and tutorials
External Links – SRH-2D Version 2.0
- SRH-2D version 2.0 Theory and User Manual
- SRH-2D version 2.0 Distribution – Includes software, manual, template file for integration into SMS interface, and tutorials
- SRH-2D version 2.2 Distribution – Includes software, manual, Template file for integration into SMS interface, and tutorials
Papers / Presentations
- SRH-2D Theory Paper
- SRH-2D Training Presentation
- 2006 FISC Paper on Savage Rapids Dam Removal Project "Comparison of Numerical Hydraulic Models Applied To The Removal of Savage Rapids Dam Near Grants Pass, Oregon"
- List of journal articles using SMS by Prof. Greg Pasternack, UC Davis
- Lai, Y. G. (2020). A Two-Dimensional Depth-Averaged Sediment Transport Mobile-Bed Model with Polygonal Meshes. Water, 12(4), 1032.
Project Reports
- Lai, Y.G., Holburn, E.R., and Bauer, T.R. (2006)."Analysis of sediment transport following removal of the Sandy River Delta Dam."
- Lai, Y.G. and Greimann, B.P. (2011). "SRH Model Applications and Progress Report on Bank Erosion and Turbidity Current Models"
Additional References
- Launder, B. E., & Spalding, D. B. (1974). "The numerical computation of turbulent flows." Computer Methods in Applied Mechanics and Engineering, 3(2), 269-289.
- Seminara, G., Solari, L., & Parker, G. (2002). "Bed load at low Shields stress on arbitrarily sloping beds: Failure of the Bagnold hypothesis." Water Resources Research, 38(11).
- Greimann, B., Lai, Y., & Huang, J. (2008). "Two-dimensional total sediment load model equations." Journal of Hydraulic Engineering, 134(8), 1142-1146.
- Phillips, B. C., & Sutherland, A. J. (1989). "Spatial lag effects in bed load sediment transport." Journal of Hydraulic Research, 27(1), 115-133.
- Van Rijn, L. C. (1984). "Sediment transport, part I: Bed load transport." Journal of Hydraulic Engineering, 110(10), 1431-1456.
- Van Rijn, L. C. (1984). "Sediment transport, part II: Suspended load transport." Journal of Hydraulic Engineering, 110(10), 1613-1641.
In the News
- News article "Computer Modeling Smoothes a Dam Hard Job"
- Photos related to Wired.com article "Computer Modeling Smoothes a Dam Hard Job"
External Links – SRH-W
SRH-W Version 1.1
Papers / Presentations
Related Topics
SMS – Surface-water Modeling System | ||
---|---|---|
Modules: | 1D Grid • Cartesian Grid • Curvilinear Grid • GIS • Map • Mesh • Particle • Quadtree • Raster • Scatter • UGrid | |
General Models: | 3D Structure • FVCOM • Generic • PTM | |
Coastal Models: | ADCIRC • BOUSS-2D • CGWAVE • CMS-Flow • CMS-Wave • GenCade • STWAVE • WAM | |
Riverine/Estuarine Models: | AdH • HEC-RAS • HYDRO AS-2D • RMA2 • RMA4 • SRH-2D • TUFLOW • TUFLOW FV | |
Aquaveo • SMS Tutorials • SMS Workflows |