SMS:SRH-2D: Difference between revisions
Line 167: | Line 167: | ||
===Additional References=== | ===Additional References=== | ||
* [https://www.researchgate.net/profile/Brian_Launder/publication/222171646_The_Numerical_Computation_of_Turbulent_Flow_Computer_Methods/links/5720899108aeaced788bc5d8.pdf Launder, B. E., & Spalding, D. B. (1974). The numerical computation of turbulent flows. ''Computer methods in applied mechanics and engineering'', 3(2), 269-289.] | * [https://www.researchgate.net/profile/Brian_Launder/publication/222171646_The_Numerical_Computation_of_Turbulent_Flow_Computer_Methods/links/5720899108aeaced788bc5d8.pdf Launder, B. E., & Spalding, D. B. (1974). The numerical computation of turbulent flows. ''Computer methods in applied mechanics and engineering'', 3(2), 269-289.] | ||
* [http://onlinelibrary.wiley.com/doi/10.1029/2001WR000681/full 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).] | * [http://onlinelibrary.wiley.com/doi/10.1029/2001WR000681/full 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).] | ||
=== In the News === | === In the News === |
Revision as of 17:51, 7 February 2017
SRH-2D | |
---|---|
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.
Currently, it is recommended that the SRH-2D model only be used with the 64-bit version of SMS. Errors may occur if using the 32-bit version.
The SRH-2D model is included in the Community Version of SMS for 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 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 multiple channel systems, and with complex floodplains. A 2D model may also be needed if one is interested in local flow velocities, eddy patterns and flow recirculation, lateral variations, flow spills over banks and levees, and flow diversion and bifurcation.
- 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 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.
- 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 condition coverage which holds boundary arcs and flow data, a materials coverage which maps material types defined for a region to each cell/element defined in the mesh, and a monitor points coverage which specifies locations where results 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 modification 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, a 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 option of a wall or symmetry is 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 o floodplain materials.
A sediment transport simulation starts by setting up and running an SRH-2D hydraulic model. 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.
Sediment transport properties are specified in a boundary condition coverage by setting the coverage BC Types and specifying the Run Type as Mobile. This allows the specification of the sediment parameters and requires sediment data on the boundary conditions which include inflow.
Additional monitor lines (in the boundary condition coverage) and monitor points (in the monitor points 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. Sediment materials can also be specified as nonerodible.
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 difficulty in executing the model.
Some common errors in running SRH-2D include:
- 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.
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 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
- 2005 US-China Workshop Paper[1]
- 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"
- FISC 2006 Short Course Presentation
- Using Bathymetric LiDAR and a 2-D Hydraulic Model to Identify Aquatic River Habitat
- List of journal articles using SMS by Prof. Greg Pasternack, UC Davis
Project Reports
- Bountry J.A. and Lai, Y.G. (2006)."Numerical modeling of flow hydraulics in support of the Savage Rapids Dam removal."
- 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 Bountry, J.A. (2006). "Numerical hydraulic modeling and assessment in support of Elwha Surface Diversion Project."
- Lai, Y.G. and Bountry, J.A. (2007). "Numerical modeling study of levee setback alternatives for lower Dungeness River, Washington"
- 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).
In the 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
- SRH-W version 1.1 User Manual
- SRH-W version 1.1 Distribution Package – Includes Software, Manual, and Tutorials
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 |