Quantcast
Channel: The RAS Solution
Viewing all 200 articles
Browse latest View live

Shifting web imagery with world files for HEC-RAS

$
0
0

Written by Krey Price  |  Surface Water Solutions
Copyright © The RAS Solution 2019.  All rights reserved. 

Using RAS Mapper as a GIS tool Part 2:
Shifting Web Imagery with World Files



Have you ever tried pulling web imagery into RAS Mapper only to find that it doesn't line up with your existing shape files, terrain surfaces, or other geospatial data?  In this second of a three-part series on using RAS Mapper as a GIS tool we'll cover the modification of world files for georeferencing.

Photo by Franck V. on Unsplash

In some cases, the misalignment of web imagery may result from an incorrectly applied projection file. If needed, projection files can be downloaded for free from spatialreference.org and other online data sources. You can check the projection of your terrain file by double-clicking on it in RAS Mapper, then select the "Source Files" tab and view the metadata using the Info column. 



If you have confirmed that the correct projection file is being used and you are still having issues, there are a few additional workarounds you may wish to try. One option is to select the "Alternate HEC-RAS Raster Warping Method" under "Tools | Options" in RAS Mapper. This method applies the gDAL OGR vector reprojection. I have found this method useful for resolving discrepancies in European data sets in particular.





In some cases, the alternate method doesn't do the trick either. RAS Mapper uses the ESRI projection file format, which as I understand it does not recognize TOWGS84 parameters and in some cases improperly applies the proj.4 projection parameters. In any case, you may find yourself stuck with misaligned web imagery.


One workaround is to save static images of your web imagery and then shift the image to the correct location by making adjustments to the world file (Wikipedia has a good summary of the six lines of code comprising a world file here).


Here's a video walk-through of the process:


https://youtu.be/5REWQ8Z_08I


The "World File Calculator" spreadsheet referenced in the video can be downloaded here:


http://www.surfacewater.biz/wp-content/uploads/2018/12/Surface-Water-Solutions-World-File-Calculator.xlsx


Here are the steps covered in the video:


1. Adjust transparency as needed so that you can see a known point in both the web imagery and in the survey or LiDAR data.
2. Using the measure tool in RAS Mapper, click on a known point from the web imagery (with incorrect location), then double-click on the the same known point based on the correct location.
3. Select "Copy coordinates to clipboard".
4. Zoom to preferred extents for results and save view (click here for further details on saved views and static imagery).
5. Right-click on web imagery and select "export layer". Note: RAS Mapper currently has two export options: tif and jpg. The tif format will embed the metadata and georeferencing information into the encoded text of a geotif file without generating a separate world file. A tfw file can still be created to override the embedded georeferencing, but I suggest using the jpg format so that the world file is automatically generated. Save the image (I suggest creating an "aerial imagery" subfolder within the current project's directory structure to store static images).
6. Open the world file calculator spreadsheet.
7. Paste the coordinates from RAS Mapper into the blue cell.
8. Open the newly created jgw file (jpg format world file) in a text editor and copy the values.
9. Paste the values into the "Original World File" column of the spreadsheet.
10. Copy the "New World File" values over the original values in the jgw file, save, and close.  
11. Under Map Layers in RAS Mapper, right-click and select "Add Existing Layer". Be sure to drag down file types to show all image files and browse to the newly created file.
12. Adjust transparency as needed to confirm common points are now co-located.

The process can be repeated to fine-tune the positioning. Because a static image won't be pyramided like web imagery, you will unfortunately need to complete these steps at every desired zoom level and view extent, the but the process of saving static images allows you to control the appearance of your figures without having to worry about losing your internet connection while presenting or viewing results.


Note that these steps assume the units and north alignment are the same (or close enough) between the data sets and that the static image can be simply shifted without scaling or rotation factors, which we'll cover in more detail in our next post.
Please let us know any feedback or suggestions for improvements or additional efficiencies in adjusting web imagery.



HEC-RAS 1D-2D Course in Atlanta

$
0
0
Hey RAS Users-

This is a "Save the Date" notice for an upcoming 1D-2D HEC-RAS 3-day training course I'll be giving in Atlanta on May 7-9, 2019.  We'll be covering everything from setting up a 2D model, a 1D-2D combined model, RAS Mapper pre- and post-processing, viewing output, river restoration, to dam/levee breaching.  The class will consist of a mix of informative lectures focusing on practical aspects of 1D-2D HEC-RAS modeling and interactive exercises.


Come join me for a great opportunity to learn the tricks of 1D-2D modeling in HEC-RAS and to network with fellow HEC-RAS Users.

If you are interested, please follow the following link and express your interest.  We will send more details, including the venue and registration instructions as soon as they are available.






Chris Goodell
Kleinschmidt Associates
The RAS Solution



Georeferencing hard copy or pdf maps

$
0
0
Written by Krey Price  |  Surface Water Solutions
Copyright © The RAS Solution 2019.  All rights reserved. 

Using RAS Mapper as a GIS tool Part 3:
Georeferencing hard copy or pdf maps





Have you ever received your background data in hard copy or pdf format and wanted to view it in its georeferenced location? For this final topic in our three-part series on using RAS Mapper as a GIS tool, we'll cover the steps required to georeference a map using world files. [For additional background see the previous posts Part 1: Terrain modification and Part 2: Web imagery with world files.] 








Before you begin any georeferencing process, always check that you are using the correct projection (and confirm the desired projection with the client or end user of the data.) If you don't know which projection to use, you may need to check with the provider of your LiDAR data or other geospatial data. If your terrain file has been provided in geotif format, another option is to add the terrain file to RAS Mapper without assigning a projection, then double click on the name of the terrain, select the "Source Files" tab, and click on the "Info" button to view the GDAL metadata. The projection should show up under the "PROJCS" tag:


If you don't have the relevant projection file available, most projections (with the exception of local project grids) are available online for free download. Projections are typically catalogued according to European Petroleum Survey Group (EPSG) codes maintained by the International Association of Oil and Gas Producers. [When you think about it, it makes sense that the oil and gas industry would have a substantial interest in pinpointing global locations for exploratory wells and other critical geospatial data!]
Online repositories for projection files such as spatialreference.org include the EPSG code with each spatial reference system:


http://www.spatialreference.org


Projection files can also be copied from the prj file associated with any shape file that matches the target coordinate system. In any case, once I've got the correct prj files, I like to place it in a separate subfolder named "projection" that I create under my project directory, and then assign it a file name such as "GDA 94 MGA Zone 55 Projection.prj" that clearly distinguishes it is a projection file (as opposed to a HEC-RAS project prj file!) I also like to make sure that the projection folder includes only a single prj file so that anyone who picks up the project in the future knows the intended projection. [RAS Mapper sometimes loses the spatial reference system association and that can avoid confusion when it is re-assigned.]
Once you've confirmed that you are using the correct projection (typically by turning on web imagery and checking that it lines up with your terrain, shape files, or other geospatial data for your project) you can georeference any hard copy, scans, image files, or pdf files within the specified projection.
The georeferencing steps are covered on our new YouTube channel in this video walk-through beginning at 6:37:

https://youtu.be/5REWQ8Z_08I
Here are the steps covered in the video: 


1. For hard copies, scan your image to a raster format (commonly jpg or tif). For pdf files, use "save as" in Adobe Acrobat and select jpg (resolution can be adjusted under settings). For other file types, you may wish to use the Snipping tool that ships with Windows and save the screenshot as a jpg (best if it's on a high-resolution screen as you'll be stuck with the screen resolution and not the original source file's resolution).
2. Create a new file in Word, Notepad, Wordpad, or other text editor or word processor (or copy an existing world file and skip to Step #5).
3. Enter six lines of numerical text. Use dummy values from "original world file" column in the World File Calculator spreadsheet (available for free downloaded here) or take an initial guess at values based on the assumed coordinates of the upper left pixel.
4. Save file as plain text format (*.txt) using the same file name as the image file (select "default settings" if prompted).
5. Change file extension in Windows Explorer to *.tfw, *.jgw, *.bpw, etc. as appropriate to match your image file format (click yes if prompted with unusable file warning).




6. Under Map Layers in RAS Mapper, right-click and select "Add Existing Layer". Be sure to drag down file types to show all image files. Browse to file and select OK.
7. Using the measure tool in RAS Mapper, measure the distance between two points that are a known distance apart (preferably from a scale bar or coordinate tick marks on the map).
8. Open the World File Calculator spreadsheet (or create your own - there's really not that much to it!) and enter the measured and actual distances under Line #1 (cells D3 and E3) and again in Line #4 (cells D6 and E6).
9. Copy values from "New World File" column, open the world file in a text editor, and paste the updated values over the old values.
10. Right click on the image file and select "remove layer".
11. Re-add the image file and adjust transparency as needed.
12. Using the measure tool in RAS Mapper, click on a known point from the original image (with incorrect location), then double-click on the same known point based on the correct location.
13. Select "Copy coordinates to clipboard".
14. Paste the coordinates into Cell C10 in the spreadsheet (Note: values will be replaced in the cell range C10:F12 but only C10:D11 are used in the spreadsheet).
15. Copy "New World File" column over the previous values in the world file.
15. Remove and re-add the image file under map layers in RAS Mapper.
16. Adjust transparency as needed to confirm common points are now co-located. Measure the scale bar to confirm correct scaling.






Hopefully the first try will get you close enough for your purposes, but the process can be repeated to fine-tune the results as needed. You can also use these steps to convert geospatial data to or from a local project grid for which projection details are unavailable.


Keep in mind that these adjustments only provide a visual approximation and should be treated cautiously if the results are to be used for permitting or construction purposes. In that case, you'd want to get a surveyor on board with the proper tools to apply all of the required warping factors and address any other discrepancies. 




The full set of parameters in a projection file can get quite complex; projecting a curved surface to a flat plane comes with mathematical difficulties that have plagued mapmakers for centuries. Here's a great video explaining why all world maps are wrong.


www.surfacewater.biz/projection/




With these difficulties in mind, it is no surprise that even if you get two points to line up perfectly, a third point can still be shifted, particularly for points located a significant distance away from the alignment points. For manual georeferencing, I suggest using alignment points that are as far apart as possible for your selected zoom extent to avoid major discrepancies.
Note on rotation factors: The above steps assume that north is straight up in your hard copy or pdf map based on the applied projection (or at least close enough to straight up to suit your purposes). If not, rotation terms can be added to the world file, but in my experience this can lead to some confusion, as the pixels shift on the fly at different zoom levels.




Although rotation terms can be added to the world file by measuring the angular difference between lines drawn between two known points, my preference is to make rotation adjustments graphically before diving into any world file adjustments. Raster images can be rotated in a number of programs (Photoshop, Paintshop, Word, PowerPoint, etc.) to achieve a proper north alignment. [The "Z axis rotation" can be used in Microsoft products to apply specific sub-degree rotation factors.] The rotation terms can then remain as zero values in the world file.
We hope you have found this useful for your work. Please let us know if you have any comments or suggestions for improving these processes! Thanks for tuning in, and as always, let us know your recommendations for upcoming blog topics.



Collecting ideas for future posts: One subject we are working on for a future post is hacks for removing flow from HEC-RAS models. Internal boundary conditions are very useful for adding flow anywhere in your model, but removing flow can be a tedious process requiring a bit of creativity. Some new features coming in Version 5.1 may simplify the process, but in the meantime I've seen some clever workarounds being applied that we would like to share with this forum; please contact me if you have done this successfully so we can perhaps feature your method in an upcoming blog post.




And speaking of Version 5.1, right at the top of my wish list for HEC-RAS updates is for project files to be called *.ras files, *.hrp (HEC-RAS Project) files or any unique extension that doesn't happen to coincide with thousands of other files on my computer (in this case, of course the ESRI-format *.prj projection files). I understand the Corps may have been first in this case (calling their project files prj's well before ESRI did) but the rest of the geospatial industry is not about to change their format, so maybe HEC-RAS project files can get a unique identity in the next version - or at least a search tool within RAS that can recognize and distinguish HEC-RAS-format project files. Let us know if you agree and maybe we can turn this into a grass roots, crowd-based request!






HEC-RAS Version 5.0.7 is now available!

$
0
0
Fellow RAS Users:  Please take note that a new version of HEC-RAS, Version 5.0.7, is now available to download.  This is purely a bug fix release, so there are no new (advertised) features.  You may download the new version here:





Please read the following announcement from Gary Brunner and the list of bug fixes described below.
  
             "We have released a new version of HEC-RAS called version 5.0.7 today.  This version is a bug fix only version (i.e. No new features).  Unfortunately there were some bugs in version 5.0.6 that were significant enough to warrant a new release of the software.  Enclosed are the HEC-RAS 5.0.7 release notes describing the bugs that were fixed.

                We are sorry for any inconvenience this may have caused you while using HEC-RAS.  We hope you enjoy the new version of the software.  As always, please keep us informed on any bugs you find, or your ideas for new features.  Here is the link to the HEC-RAS 5.0.7 download area of our webpage:


Warning:  A few new software checks have been added for Lateral Structures and SA/2D Hydraulic Connections.  When you run version 5.0.7 of HEC-RAS these data checks may stop your data set from running, until you address the issues listed in the error messages."


The following is a list of bugs that were found in HEC-RAS Version 5.0.6 and fixed for Version 5.0.7:


  1. Manning’s n Regions Override Bug. If you used the Manning’s n Region option in the RAS Geometry editor (Not RAS Mapper), there may be an issue with your final Manning’s n values not being correct within that Region. Unfortunately, in previous versions of HEC-RAS Geometry editor (5.0.6 and earlier), if you draw a Manning’s n Region polygon in the counter clockwise direction, a negative area was computed for the polygon (the shapefile standard is to store it in the clockwise direction). Because the polygon had a negative area, the polygon was actually not being used to override any of the Manning’s n values.  This problem has been fixed in version 5.0.7. Users can draw these polygons in any direction.
  2. Uniform Lateral Inflow Bug. When using a uniform lateral inflow along a 1D river reach, HEC-RAS was not correctly putting all of the flow into the river from the uniform lateral inflow hydrograph when the lateral inflow extended across a bridge/culvert or an inline structure. The bug was caused due to HEC-RAS not using the length across the bridge/culvert or inline structure to compute the total lateral structure length, which is then used in computing the flow per foot along the lateral structure. The flow is then put in between the cross sections based on the flow per foot computed, times the length between the cross sections. However, no flow was being put in between the cross sections upstream and downstream of a bridge/culvert or an inline structure. Therefore that portion of the flow was not being put into the river. This mistake was introduced in version 5.0.4 and has been in the program since this version. This mistake has been fixed for HEC-RAS version 5.0.7.  If you used a work around for this problem in the previous versions, remember to change the model back to the normal way of using uniform lateral inflows with no flow modifiers.
  3. Floodway Inundation Boundary Mapping Issue: Inundation Boundary polygons created by RAS Mapper for floodway analysis were not properly created if cross-section station-elevation data started with non-zero stationing. This problem did not occur for stationing data that starts at zero.
  4. RAS Mapper not Saving all Breakline Options Issue: The program was not saving all options when the user activates the option to “Enforce a 1 Cell Protection Radius” for some data sets. This was a problem as you continue to edit/modify your mesh, as those breaklines may not function as you intended them to.
  5. External SA/2D Hydraulic Structures Issue. The issue would happen for external SA/2D Hydraulic Connections (i.e. connecting a storage area to a 2D area, or a 2D area to another 2D area) when the length of the station/elevation data did not match the length of the structure centerline for cases when a geo-referenced centerline was used for the structure. The program connects to the 2D cells only based on the geospatial centerline and its length. If the station/elevation data has a length that is longer than this line, then the computational engine has a length of the structure that is not connected to any actual cell, and this causes the program to make incorrect flow calculations if that part of the structure is wet. The program now checks for this and does not allow a user to have a centerline and station/elevation data that are not the same length.
  6. Lateral Structure Connection to 2D Area Issue. There was a problem where the GUI was sending an extra 2D face point at either the beginning or end of a lateral structure connection. If a user tried to start a lateral structure with exactly the same coordinates as a facepoint, sometimes a round off of the numbers cause the interface to think it needed an additional facepoint in front of or at the end of the lateral structure connection to the 2D area. This was causing unsteady flow computations to go unstable with no message as to what the problem was.
  7. 1D Velocity and DV Mapping in RAS Mapper Issue. In version 5.0.6, if you have cross sections that do not start at a zero stationing, RAS Mapper was not correctly mapping velocity and Depth x Velocity in the 1D model region.
  8. Alternative Overbank Deposition Method Issue. The new overbank deposition algorithm in 5.0.6 was not connected to the Exner 5 and Active Layer methods (they used the old method by default). Only the Exner 7 method was being used (so turning it off did not revert to the old method). The new method (which limits overbank bed change to deposition for a whole computation increment, rather than a mixing time step) and the old method are now available for all three mixing methods.
  9. Max Water Surface Profile Issue. If a user had model “Warmup” turned on, the software was incorrectly including the initial condition at the beginning of the model warmup, as possible results for the maximum water surface elevations during the run. The maximum water surface elevations should not include any of the warmup time solutions. This has been fixed in 5.0.7. Also, if no model warmup is turned on, then the maximum water surface elevations should potentially include water surfaces at time zero (from the initial conditions backwater computations). This is how this works now in version 5.0.7.
  10. 1D Unsteady Flow Encroachments. The program was not allowing the user to put an encroachment at the cross section just upstream of a culvert. If the user put an encroachment upstream of the culvert, a bad water surface was being computed at this location (It was getting a water surface like 4000 ft in the air).
  11. RAS Interface is not Catching Bad Data for Internal SA/2D Hyd Conn. A 2D flow area with an internal SA/2D hydraulic structure must start and end at a face point. Enforcing it as a breakline will snap it to the face points. The program will no longer let you run a model where the internal structures do not start and stop at a face point.  WARNING: when users first try a simulation in version 5.0.7 and there are internal structures that do not start and stop at face points, then the simulation will not run and the program will tell you the locations in which you need to fix this problem.
  12. RAS Unsteady Flow Sediment Transport Computations. The software was not running any unsteady flow sediment transport models. This was due to a mistake added just before release with the HDF file reads.
  13. Bridge/Culvert Editor HTAB Curves. When you open HEC-RAS and a geometry file, the bridge/culvert curves get loaded into memory. From the Bridge/Culvert editor there is a button called "HTAB Curves" that allows you to plot these curves (and edit them). However, if you change the HTAB parameters that went into making the curves, rerun the model, and then click the HTAB Curves button, you would still see the original curves when the project was loaded, and not the new updated curves from the run.
  14. Plotting Observed Daily Instantaneous Data Problem. HEC-RAS was not plotting daily instantaneous time series data correctly for the hydrograph plots and profile plots. The data normally has meta data that says when the data was measured, for example “0600” on the day. But HEC-RAS was always plotting it at 2400 on the day.
  15. SA/2D Hydraulic Connection Culvert Issue. We received some data sets where users were drawing the centerline for a hydraulic connection from left to right in a specific direction, but then they were telling the software that the Headwater and Tailwater were the opposite of that. Additionally, they were drawing culverts from downstream to upstream. This was causing several issue.  WARNING!!! The software now checks for this issue (above) and will not run if encountered. A message will appear telling the user that the direction they drew the structure centerline is inconstant with how they have headwater and tailwater defined, or they have drawn their culverts, gates, rating curves, or flow time series from downstream to upstream. HEC-RAS requires the user to draw the structure centerline from left to right while looking downstream, and for the headwater and tailwater connections to reflect that (HEC-RAS will automatically connect it that way, but a user could change it). Additionally, all centerlines drawn for hydraulic outlets (culverts, gates, rating curves, time series of flows) must be drawn across the structure from upstream to downstream.
  16. Storage Area Initial Elevations Problem. Unfortunately, a bug was added to version 5.0.6 that stopped it from correctly initializing the elevations of storage areas whose initial elevations were left blank (undefined) in the Unsteady Data Flow editor. In 5.0.6, the program did not initialize the storage area elevations at all. SA's that are hydraulically connected to the river with a lateral structure should be set to the same water surface as the river across the lateral structure. This was not happening. All storage areas were being left dry at the start of the simulation.
  17. Storage Area Elevation Points Filter in RAS Mapper. There was an occasional issue in RAS Mapper that when a storage area elevation relationship was extracted from the terrain, occasionally it would end up with thousands of points. The points filter routine had a precision issue that was causing this.
  18. Maximum Depth for Storage Areas. If a user had storage areas that were not connected to anything else (i.e. they only had a flow hydrograph attached to it). HEC-RAS was not correctly reporting the maximum water surface elevation. The hydrograph plots were correct, but the maximum inundation plot in RAS mapper was wrong.
  19. Dredge Material Unit Weight Issue. Dredged cohesive material used the cohesionless unit weight to convert from mass to volume. Dredged material now uses the correct unit weight.
  20. Import of Vertical Manning’s n values. When importing an existing HEC-RAS model into another existing model, if the mode to be imported had vertical Manning’s n roughness specified at cross sections, those values were not being imported.

HEC-RAS 1D-2D Course in Atlanta, Georgia

Those Nuisance Cell Errors and Wedge Depressions

$
0
0

Written by Chris Goodell | Kleinschmidt Associates
Copyright © The RAS Solution.  All rights reserved

If you've been doing some 2D modeling, you most likely have run into what I call nuisance cell errors.  You see them while your model is running and they can really slow down your simulation as the RAS chugs through iterations to try to come to a solution.

The errors here are in the second to last column and are in units of feet.  So as you can see, they are quite small, and normally I wouldn't be too concerned with them.  However, if the errors persist, they can really slow down your simulation time.  Sometimes adding several hours to the overall time, in my experience.  So while the errors are small, there is a need to get rid of them, at least the ones that persist for many timesteps.  

A common cause of nuisance cell errors is what I call "wedge depressions".  In the following figure, notice the wedge of water that sits below a well defined ridge line.  The red arrow points to it.  


As you can see by the blue hydraulic connectivity lines that there is hydraulic connectivity between Cell 1 and Cell 2.  But water fills up cells from the lowest spot, and in Cell 1, this is below the ridge.  So you get a narrow wedge of water (wedge depression) that builds up causing some local fragmentation and a resulting error.  Notice the water just abruptly stops at the cell face between Cell 1 and Cell 2.  Wedge depression doesn't always cause an error (by error, I mean RAS cannot converge on a solution within the error tolerance in 20 iteration tries).  But it did in Cell 2.  Notice Cell 4 has the same wedge depression, but it was not producing errors in this simulation (ahhh, the mysteries of HEC-RAS!).  

Fortunately, the fix for this is easy.  You draw a breakline along the ridge line so that a cell face resides on the ridge line itself.  This removes the low spot (and thus the wedge depression) in Cell 2.  the next figure shows the placement of the breakline (but before enforcement).  In this case, to preserve the spacing of the cells on the north side of the line, I specified a cell center spacing of 55 feet on the breakline (actually, I started with 40 feet, then 50 feet, but these produced the dreaded "red dots".  55 feet worked.)


After a bit more tweaking of the breakline alignment, here is the result with the breakline enforced and the wedge depression eliminated.  Notice without the wedge depression, flow is able to move more naturally towards the upper left from what used to be Cell 1.  Oh yeah...and no more nuisance cell errors!  


Yes...there's still some more fragmentation that could be addressed, but this is early on in the simulation and my objective is to ultimately present the max ws maps.  So my concern at this point in the simulation is to have decent results with fast simulation times.  I'm not concerned with little bits of fragmentation that don't throw errors.  

HEC-RAS Pub and Grub in Portland-Save the Date!

$
0
0

For those of you who will be in the Portland Oregon area, I invite you to join me for this special event Kleinschmidt is hosting at the Lucky Lab on Hawthorne.  Krey Price will be in town from Australia and we along with Sean Welch of BPA will be giving fun and engaging presentations on HEC-RAS.  This will be an informal, casual gathering with lots of time for Q&A.  Hope to see you there!  No RSVP needed.  Just show up, get a pint, and grab a chair.  



McBreach Probabilistic Dam Breach Modeling Software

$
0
0
Fellow HEC-RAS modelers. Today is National Dam Safety Day.  Dam Breach modeling is a key component to a well-rounded and robust dam safety program.  If you're doing dam breach modeling, you need to check out McBreach. Probabilistic hydraulic modeling is where we're heading with dam breach analysis. Why not get a head start with McBreach. This is free software and works seamlessly with HEC-RAS.  https://spark.adobe.com/page/RXu9m2OZJFOuc/


2D HEC-RAS Training - Atlantic City September 24-26, 2019

$
0
0
I hope you will be able to join me in Atlantic City for our second hosted HEC-RAS course on 2D and combined 1D/2D Modeling.  The one we did back in May in Atlanta Georgia had a great turnout and was a lot of fun.  Don't miss out on this great opportunity to learn 2D HEC-RAS, network with other modeling professionals, and to have a great time doing it!

Don't delay, this one will fill up fast!

https://www.eventbrite.com/e/1d-2d-modeling-with-hec-ras-professional-development-course-tickets-67259295363


Go here for more information and to sign up: 

HEC-RAS 5.0.4 is released!

$
0
0

Message from HEC on the new release of HEC-RAS version 5.0.4.  Make sure to also check on the new Verification/Validation and the Benchmarking documents. 

***Note.  If you originally downloaded 5.0.4 before May 4 at 4:00 pm PDT then please redownload and install from the links below.  There were some minor errors found in the original release.

Download HEC-RAS Version 5.0.4 here:

Download RD-52  HEC-RAS Verification and Validation Tests document here: 


Download RD-51 Benchmarking of the HEC-RAS Two-Dimensional Hydraulic Modeling Capabilities here:

There are several other new useful documents related to HEC-RAS here:
http://www.hec.usace.army.mil/software/hec-ras/documentation.aspx


From HEC:
We have released the final version of HEC-RAS 5.0.4 today.  It is now available for download from our webpage at the link below.  I have enclosed the 5.0.4 release notes for your review.  I have also included a separate document that describes how to use the new features contained in this version. Additionally  I have also included a document on how to spec out a computer for doing 2D modeling with HEC-RAS.

            We hope you enjoy the new version of the software.  As always, please keep us informed on any bugs you find, or your ideas for new features.  Here is the link to the HEC-RAS 5.0.4 web page:


            The following is a list of the major new features you can read about in the Release Notes and the Supplemental User’s Manual:

  1. New RAS Mapper Pre-Processor Tools
  2. Variable Time Step Capabilities (improved accuracy and faster run times)
  3. Increased Computational Speed for 2D Computations (2 to 5 times faster run times)
  4. Enhancements to SA/2D Hydraulic Connections
  5. New 2D Polygon Mesh Refinement Tool
  6. Internal Boundary Condition Lines for 2D Areas
  7. New Velocity Term for 2D Flow Area Boundary Conditions
  8. 64–bit Version of RAS Mapper and 64-bit 2D Pre-Processor (RAS now can handle much larger 2D models)
  9. Improved Terrain Generation Utility
  10. Improved Manning’s n Layer Development
  11. Sediment – New Options in the User Defined Grain Class Editor
  12. Sediment – New Transport Function Calibration Options

WARNING: Some of the bug fixes “may” cause the computed results from 5.0.3 to 5.0.4 to change.  Please read the Release Notes carefully for information about these changes.


There are two new verification and validation documents for HEC-RAS now available on our web page in the HEC-RAS area.  Here is a short description of each:


  1. RD-52  HEC-RAS Verification and Validation Tests (Enclosed in this email)

            This is the first ever full blown Verification and Validation document for HEC-RAS.  In this report we compare HEC-RAS to: analytical and text book data sets; lab studies, and real world flood events with good observed data.  This report covers 1D steady flow; 1D unsteady flow, and 2D unsteady flow computations.  The report verifies that HEC-RAS correctly solves the equations used in each of these computational frameworks, and it verifies that HEC-RAS can be used for a very wide range of riverine modeling conditions.  It can aslo be download here:



  1. RD-51 Benchmarking of the HEC-RAS Two-Dimensional Hydraulic Modeling Capabilities

This research document summarizes how the Hydrologic Engineering Center's (HEC) River Analysis System (HEC-RAS) software performed in the two-dimensional (2D) modeling benchmark tests developed by the United Kingdom's (UK) Joint Defra (Department for Environmental Food and Rural Affairs) Environment Agency.  This is an updated version of this report that we previously developed.  We re-ran all of the test, plus one new test, with version 5.0.4 of HEC-RAS.  The results are even better than before, and the HEC-RAS computational speed is much faster than previously document.  This report can be downloaded from our web page at the following link:



2D HEC-RAS Class in Boise October 28-30

$
0
0
HEC-RAS Modelers!  Don't miss this opportunity for 2D HEC-RAS training in beautiful Boise Idaho.  I'll be teaching this class and hope to see you there!

There are limited seats for this class so don't delay in signing up.

This will be a three-day course with a mix of lectures and workshop exercises.  You will learn how to set up, run, and troubleshoot 2D and combined 1D-2D HEC-RAS models with an emphasis on river restoration projects.  Five separate workshops will give you practical experience setting up and running the software.



1D/2D Modeling with HEC-RAS
Professional Development Course

Purpose and Background
This intensive, workshop-oriented, three day seminar will prepare the engineer and water resource professional to use the HEC-RAS computer program for modeling two-dimensional (2D) and combined one-dimensional/two-dimensional (1D/2D) unsteady flow applications. Led by Kleinschmidt Associates, participants will learn how to approach and construct a 2D model for unsteady flow conditions, and to effectively view and analyze results. 

The seminar includes lectures on 2D flow theory, RAS Mapper, an introduction to the new capabilities and features of HEC-RAS, post processing and analysis of results, and procedures for creating a stable and calibrated 2D model.  Workshops focus on giving students hands on experience with building and pre-processing the computational mesh, performing offline and inline 2D projects, and using the built-in feature RAS Mapper to spatially analyze results. The seminar can also be customized to fit an organization or groups interest in areas such as dam breach analysis, river restoration, or rain-on-grid.


HEC-RAS incorporates various aspects of 2D hydraulic modeling, including floodplain hydraulics, dam and levee breaches, rain-on-grid applications, and the interaction between 1D and 2D systems.  Version 5.0.7 of HEC-RAS includes features in 2D open channel hydraulic analysis such as:
  • Two-dimensional flow analysis using the full St. Venant or diffusion wave equations in 2D.
  • Ability to perform a combination of 1D and 2D flow analysis in the same model.
  • The use of unstructured or structured computational meshes for the 2D flow areas.
  • Dam and Levee breaching in 1D and 2D areas.
  • Rain-on-grid modeling.
  • Full pre- and post-processing of geometry in 1D and 2D. 
Seminar Benefits/Learning Outcomes
  • Use the HEC-RAS (River Analysis System) computer program to model 2D and combined 1D/2D unsteady flow hydraulics.
  • Use RAS Mapper to pre-process both 1D and 2D geometric elements. 
  • Understand 2D flow theory and the differences between 1D and 2D modeling.
  • Gain hands-on HEC-RAS experience by participating in practical computer workshops.
  • Understand how to develop a stable and calibrated 1D/2D flow model.  
  • Know how to post-process and analyze 1D and 2D results. 
  • Obtain valuable insights in methods for minimizing computation errors and instabilities for 2D unsteady hydraulic models.  
  • Learn from real world projects and applications.
Who Should Attend
Consulting engineers, water resource planners, engineers employed by local, state, or federal government agencies. Participants should have some experience in floodplain hydrology and hydraulics, and some experience in HEC-RAS steady and unsteady flow computer modeling. They should also be able to follow simple computer instructions. 


DAY 1
  • Introduction to HEC-RAS 1D/2D Modeling
  • Building a 1D/2D Model
  • RAS Mapper
  • Computer Workshop on Creating a 2D Model
  • 1D/2D Flow Data
DAY 2
  • Computer Workshop on Channel and Floodplain Modeling
  • 2D Theory
  • Computer Workshop on Dam Breach
  • Performing Computations
  • 2D Output
DAY 3
  • Computer Workshop on Levee Breach Modeling
  • 2D Advanced Topics
  • 2D Modeling Stream Restoration
  • Computer Workshop on River Restoration

HEC-RAS in Three Dimensions

$
0
0
Now this is the kind of post that gets you on The RAS Solution!  Cutting edge creative solutions to complex problems.  This is great stuff.  Thanks Michael for all the effort you put in to this.

Written by Michael Link | AECOM
Copyright © The RAS Solution.  All rights reserved.  

In today’s post we’ll cover two methodologies that leverage HEC RAS results to create impressionable and multi-dimensional flood visualizations. If these methodologies are implemented industry-wide, it is likely that the public will become more engaged with our work and our clients wowed. Products of the two methods are shown below.

Figure 1: 3DV-Flo Output of Muncie, Indiana (left) and 3D Printed Flood of Puerto Rico (right)

Introduction
The mitigation of flood risk is traditionally accomplished through regulation, correct science, and public buy-in. The federal government writes the regulation, we make the floodplain maps, and homeowners begrudgingly purchase flood insurance. How can we as engineers change the public’s sentiment towards our work and how can we further reduce infrastructure damage? One way we can do this is by transforming HEC RAS results into a format that is visually stimulating and easily digestible for nearly anyone. 

FIRM maps, WSEL/Depth rasters, and inundation boundaries serve the purpose of communicating information but do not necessarily communicate flood risk in a simple and powerful manner. This post is targeted at transforming RAS results into just that, a simple and powerful resource for capturing the attention of the public. If flood risk is internalized by communities, we can expect infrastructure loss to decrease and funding for our line of work to increase. Let’s dive in!

3DV-Flo (3D Visualization of Flooding)
3DV-Flo is a methodology that takes in RAS results and generates a 3D flood visualization in Google Earth. Here are a few examples of 3DV-Flo output.

Figure 2: Austin, Texas flooding visualized with 3DV-Flo


Figure 3: Google Street View of 3DV-Flo output


Figure 4: Muncie, Indiana flooding visualized with 3DV-Flo

The following YouTube playlist shows examples of 3DV-Flo output (the first video in the playlist is embedded below).  If you would like to explore the output yourself in Google Earth, the non-regulatory KML’s can be downloaded here.



3DV-Flo allows the information stored in WSEL and depth rasters to be viewed three-dimensionally in Google Earth. This output is an improvement upon RAS Mapper generated KML’s in that the results are no longer clamped to the ground. The 3DV-Flo method relies on three inputs (an inundation boundary, depth raster, WSEL raster), two free software (HEC-RAS, Google Earth), and one proprietary software (ESRI ArcGIS*). The 3DV-Flo toolbox has been tested on ArcMap for Desktop versions 10.3.1 and 10.6.1. Compatibility of this toolbox cannot be guaranteed with ArcGIS Pro or other Desktop versions. Additionally, there are two model types within the 3DV-toolbox. The first uses tools from the Advanced ArcGIS License. The advanced model is slightly faster and provides the functionality of adding breaklines to your resulting mesh. Breaklines are helpful for modelling levees, elevated roads, and dams. The second model uses tools solely from the Basic ArcGIS license, but you must have the 3D Analyst and Spatial Analyst extensions activated as described in this ESRI post. If this tool has a red ‘X’ next to it in ArcGIS follow the steps discussed in this video.

Assuming you have a georeferenced 1D or 2D RAS model and access to ArcGIS, let us proceed with the tutorial. You can follow along with the video below or you can scan the detailed written instructions.
*3DV-Flo methods are currently dependent on ESRI ArcGIS. These methods can surely be recreated in QGIS or in a standalone script given enough time and ingenuity. Any individual with the time and desire to make this a fully open-source method can reach out to Michael Link for further guidance. 


3DV-Flo (3D Visualization of Flooding) Steps:


                  1.  Download 3DV-Flo Files

a.       Download 3DV-Flo zipped files from Github here

                  2.  Within HEC RAS

a.       Open geo-referenced HEC RAS project
b.       Open RAS Mapper
c.       Import best available terrain if not already there
d.       Right-click results and click ‘Add new results layer’ for the max inundation boundary, depth raster, and WSEL raster
e.       Compute/update layers
f.        Export break lines as shapefile if applicable

                  3.  Within ArcGIS

a.       Open blank MXD
b.       Import model output (inundation boundary, depth raster, WSEL raster, and breaklines) to MXD
c.       Import symbology shapefile. Adjust symbology to be blue with no border and displayed with a transparency of 40%
d.       Open the 3DV-Flo toolbox by navigating there in ArcCatalog
e.       Open the 3DV-Flo_WSEL tool
                                                               i.      The 3DV-Flo_Depthtool is used to reference flood depths to the ground rather than to mean sea level. This tool can be used when the vertical datums between RAS and Google Earth differ.
                                                             ii.      To use this tool right click the KML in Google Earth, click properties, click altitude, and switch the altitude to be ‘relative to ground’
f.        Populate tool parameters and specify where the Google Earth KMZ is to be saved
g.       If breaklines are not applicable, then bring in the ‘Arbitrary_Breakline’ shapefile from the unzipped folder
h.       Run 3DV-Flo tool

                 4A - Within Google Earth for Desktop

a.       Download software here if not installed on your machine
b.       Open Google Earth
c.       Click File>Import>and then navigate to the KMZ you created in ArcGIS

                4B - Within Google Earth for Chrome

a.       Open Google Earth in chrome browser
                                                               i.      https://Earth.Google.com/web/
b.       Click three horizontal lines in the top left of window
c.       Click my places
d.       Click import KML file
e.       Click open file
f.        Navigate to KML file of interest and import

There are two known 3DV-Flo ‘bugs’. The first of which deals with poor terrains spatially restricting results. The second deals with 3DV-Flo output that appears unrealistically high in Google Earth. The first bug can be fixed by updating the underlying terrain. The second bug can be fixed by relating 3DV-Flo output relative to the ground rather than to mean sea level. These bugs and their respective workarounds are discussed in this video.

The vision for 3DV-Flo in the future includes 1. widespread adoption of the method by industry, 2. conversion of the entire US regulatory floodplain into 3DV-Flo format, and 3. the combination of 3DV-Flo with the National Water Model forecasts to help cities prepare for incoming floods. Advancing goal #1 is in part achieved through the readership of this post. Advancing goal #2 is theoretically possible by splitting the NFHL polygon layer nationwide by all XS’s with regulatory elevations. The resulting polygons could be transformed into 3DV-Flo output. Lastly, goal #3 has been prototyped in this video and this post. Any assistance to advance these goals is welcomed.

3D Printed Floods
According to the Wohlers Report 2019, the 2020 value of 3D Printing and Additive Manufacturing will be $15.8 billion! Can we get a piece of that??? One way that we can participate in this trend is to overlay RAS results onto 3D printed terrains. Here is a YouTube playlist link to some examples of just that for a pluvial flooding model of Puerto Rico (the first video in the playlist is embedded below).



The intended effect of this visualization is to intuitively convey how floods develop and where they pose the largest threat. 3D printed terrains are often presented at a large scale due to the coarse nature of the underlying digital elevation model. This visualization is unlikely to be used as a structure by structure evaluation of risk. It is more likely to be used in a city hall meeting or in the classroom. The Austin, Texas Watershed Concepts Group from AECOM was the first (to my knowledge) to prototype a 3D Printed Flood. Details on their FEMA funded San Marcos, Texas project can be found here. That same group generated a second 3D printed flood for Puerto Rico to give locals a high-level understanding of, go figure, watershed concepts. The projector and mount costed roughly $500. The San Marcos 3D print costed $1400 and the Puerto Rican model $250. The difference in price was related to print infill density and size.

If you are interested in making one of these models you can follow along with the video tutorial below and/or you can scan the detailed written instructions.


3D Printed Flood Steps:

                  1.  Develop 3D Printed Terrain

a.       Check 3D model repository (Thingiverse or Google searching 3D printed terrain of “Location”) to see if your location already has an STL 3D print file
b.       If there is no pre-existing STL file or you would like to make your own, you have a few options. 
                                                               i.      Create STL from the Terrain2STL (Ideal due to time saved and nonideal due to fixed rectangular shape of STL)
                                                             ii.      Create terrain from Touch Terrain (Ideal due to time saved and ability to specify bounding coordinates)
                                                           iii.      Create STL by converting Lidar to greyscale image (process described here) and load into Blender to convert the greyscale into an STL file (process described here).
                                                           iv.      More options detailed here.
c.       Modify STL file as needed in Meshmixer
                                                               i.      The most important functions in Meshmixer are:
1.       Edit > Plane Cut – This function allows you to discard unnecessary pieces of your model and to break your model down into smaller chunks to be printed by small 3D printers
2.       Edit > Transform – This function allows you to stretch your model solely in the Z direction. For more hash marks click the up arrow.
3.       Analysis > Units/Dimensions – This function allows you to rescale your model proportionally in all directions
d.       If you do not have a 3D printer at home or in the office, shop around to see where you can get the best deal. For prototyping, I found that 3D Hubs was the cheapest and easiest to use. Their banana reference was amusing and useful for catching extremely large or small prints. For the final 1.5-foot-long print of Puerto Rico, I used a local printing service so that I could guarantee the quality and interact with a human. Here is a list of the top 10 online 3D printing services in 2019.

                  2.  Develop Movie of Flooding in HEC RAS

a.       Download ShareX or comparable screen capture software
b.       Create or load a pluvial or fluvial model in HEC RAS
c.       Open RAS Mapper
d.       Load web imagery, terrain, or basemap to be displayed in movie
e.       Open ShareX and create screen capture video of unsteady RAS simulation


3.  Project Movie of Flooding onto 3D Printed Terrain
a.       Download VLC media player or comparable video software
b.       Open the flooding movie in VLC
c.       Connect the computer to projector via HDMI cable
d.       Setup a projector and table mount
e.       Pause flooding movie at peak flooding and calibrate the movie extent to 3D print extent by
                                                               i.      Raising or lowering the Pixar lamp stand
                                                             ii.      Adjust the aspect ratio of video
f.        Equipment setup is detailed further in this video
g.       Project flooding onto terrain with basemap


In Conclusion
As illustrated in this post, the impact of our RAS modeling can be greatly bolstered with spectacular Google Earth imagery and mesmerizing 3D prints. By presenting our clients and communities with these flood visualization resources, our work can be used outside of the narrow confines of designating whether someone should or should not buy flood insurance. 

Acknowledgements
The inspiration for the 3DV-Flo methodology would not have been possible without the catalysts listed below.
1.       a Scottish scientist’s Google Earth clamped-to-ground flood simulation
2.       a breakline/polygon GIS tool from Ryan Dalton
3.       a sea-level rise Google Earth tutorialfrom David Sadoff
4.       and a Google Earth flood simulationfrom Mariusz Krukay

Much gratitude goes out to Muhammad Ashraf and Yuxiang Kang for technical guidance and edits, John Wade for wise GIS counsel, Yacoub Raheem, Clint Kimball, and Justin Baker for patient 2D HEC-RAS Training, and Chris Wright for overall mentorship.


Contact Details:
For more information on this topic, please message Michael Link through LinkedIn - https://www.linkedin.com/in/michael-link-b88329118/

Michael Link works for AECOM in Austin, Texas.  He received a B.S. in Ecological Engineering from Oregon State University and an M.S. in Environmental and Water Resources Engineering from The University of Texas at Austin.  Michael is an EIT and certified floodplain manager.  Outside of work,  enjoys combining hydraulics, spatial analysis, and data science in novel ways.


McBreach Probabilistic Dam Breach Modeling - Free Webinar

$
0
0
ICEWaRM will be hosting this free webinar where I'll be discussing the software program McBreach and its application to probabilistic dam breach modeling.  Come join in!  Please note the the webinar start time for your particular time zone.  Starts at 10:30 AM local Sydney Australia time on Wednesday, Nov 6th.  That is 3:30 PM on Tuesday, Nov 5th if you live on the west coast of the USA.  Find your local time.  After the webinar, get your free copy of McBreach and the User's Manual here.




Manage risk by leveraging the power of HEC-RAS through Monte Carlo simulations

Introducing McBreach! Faced with aging infrastructure, limited resources, dynamic weather events, and the uncertainty of climate change, today’s dam owners must manage risk to protect the lives of the public and our built environment. The free McBreach program, quantifies uncertainty and informs better decision making before an event.


“McBreach takes dam breach modelling and analysis to the next level. By complimenting the traditional deterministic dam breach methods, McBreach’s probabilistic approach provides a more meaningful and more genuine way to communicate the potential effects from a dam or levee breach.” says Chris Goodell, Principal Consultant for Hydraulics & Hydrology, Kleinschmidt.

Very fine 2D Modeling

$
0
0
Written by Mike Hross and Chris Goodell | Kleinschmidt Associates
Copyright © The RAS Solution.  All rights reserved. 


I'm a big fan of testing HEC-RAS to the limits.  What can it do?  What can't it do?  These are questions we should all ask at the beginning of a project.  Can HEC-RAS actually prove to be a useful tool in answering hydraulic-related questions for our specific problem?  Shortly after coming to work at Kleinschmidt I was posed with the question, "Can we use HEC-RAS to model a nature-like fishway?  And more specifically, can we extract data from HEC-RAS to inform us of velocities and shear stresses between the many boulders and in the many pools that make up such a complicated feature?"  The answer was YES!  Check this out:

Here's an example of making HEC-RAS work at a very small cell level, in order to get detailed results for design considerations.







The fishway is constructed using 11 concrete weirs and is approximately 300 feet wide by 400 feet long. The boulders are approximately 5 feet in diameter, with 2.25 feet projecting above the crest of the weirs. The gaps between the boulders range from 3 to 5 feet. In total, there are over 700 boulders in the fishway. 

For the simulation shown, the drop from normal pond to the tailwater downstream of the last weir is approximately 6.4 feet or 0.58 foot per pool. The fishway is passing 890 cfs, 237 cfs of which is spilling over the top most weir from the headpond and the remainder is being supplied by 5 overflow gates. The 5 gates have varying levels of submergence on their tailwater sides, depending upon each one’s location, and they are each passing from 108 to 137 cfs. The velocity varies throughout the fishway, but is generally less than 6 feet/second. The nominal cell size in the mesh is 3 feet but decreases to 0.75 foot on the tops of the weirs to capture hydraulics between the boulder gaps. The largest cells have sizes of 96 feet and are located primarily in the headpond away from the fishway. 

The model is run using the full momentum equation, with an eddy viscosity mixing coefficient of 0.44. The computational timestep is 0.1 second and the Courant Numbers max out at around 0.6, with the highest values being in the cells on the downstream sides of the gates where water is flowing into the fishway. The timestep is very small, but the model does not need to simulate a very long time since it is a quasi-steady model, meaning there are constant boundary conditions (a constant inflow hydrograph at the upstream end of the model and a stage hydrograph at the downstream end of the model).  The model achieves steady state in less than 1 hour of simulated time.  Model run-times averaged about 3 hours.  



HEC-RAS Pub & Grub 2! HEC is coming to town!

$
0
0
HEC-RAS...Portland...Food...Beer.  They seem to just go perfectly together.  That is why we're hosting...

HEC-RAS Pub & Grub 2

You will want to save this date!  On Feb 20, 2020

HEC's very own Gary Brunner, Cameron Ackerman, and Mark Jensen

will be joining me in Portland Oregon to speak at the next offering of the HEC-RAS Pub & Grub. The HEC-RAS developers themselves!  They're going to present all of the latest features coming out in the soon-to-be-released HEC-RAS Version 5.1 and plans for future versions.  

Bring your HEC-RAS questions and come join us talking HEC-RAS at this awesome venue over some food and drink.  

And as a bonus, the Craft Brew Alliance's pH Experiment will be on hand to share some of the latest innovative brews they have in the works.    And this HEC-RAS Pub & Grub will have a very special HEC-RAS-related musical feature you won't want to miss.  Yes, that's a teaser. 😊  

For now, just mark your calendar and keep an eye out for more details to come.  If you want more information or are interested in helping to sponsor the HEC-RAS Pub & Grub, please email me Chris.Goodell@KleinschmidtGroup.com or GinaRenee Autry at GinaRenee.Autrey@KleinschmidtGroup.com



Galileo and Missoula Flood Modeling

HEC-RAS Pub & Grub - HEC-RAS 5.1 and Beyond

$
0
0
Don't miss this great opportunity to listen to and interact with the HEC-RAS Developers in a fun and casual atmosphere!  Gary Brunner, Cameron Ackerman, and Mark Jensen will all be speaking at this Pub & Grub on Feb 20, 2020 from 5:30 pm to 8:30 pm in Portland Oregon.  They'll be covering all of the latest features being added to the upcoming Version 5.1 and beyond!  Show up, get some food, grab a pint and enjoy!

This is a free event, but we need you to sign up on our Eventbrite page so we can get a good head count.  Space is limited.  Hope to see you there!


Thanks to our Platinum sponsors Wolf Water Resources, Mead & Hunt, and Galileo and our Gold sponsors David Evans & Associates, OTAK, and Geosyntec.  Your generous contributions have made it possible to bring our guest speakers and provide the food.  A special thanks to the Craft Brew Alliance and pH Experiment for providing the awesome venue.  




1D/2D Offline Connection Issues

$
0
0
Written by Chris Goodell | Kleinschmidt Associates
Copyright © The RAS Solution.  All rights reserved. 

I’ve mentioned this before, but I am thoroughly impressed with the robustness of the finite volume solution scheme used in 2D areas in HEC-RAS.  As long as your Courant Numbers are in a good spot, you rarely get errors and instabilities in 2D areas.  HOWEVER, the transition from 1D to 2D and back is another story.  In fact, the majority of errors I get in 1D/2D models occur there (either at the cross section to 2D area interface for inline connections, or in the cells adjacent to the lateral structure for offline connections).  In this post, I’m going to talk about the errors next to lateral structures at 1D/2D boundaries.

Take this for example.  The model runs well with friendly blue bars everywhere.  But I just can’t let go of the 5.154-ft error.  It only happens once and obviously doesn’t cause the model to crash.  But it’s there, and a 5-ft error is a little more than I can stomach. 

My usual go-to output for cross section errors, the water surface profile plot, gives me no clues as to what is wrong.  In RAS Mapper, I have what appears to be a flooded situation around the time where the error occurs, but reasonable looking results.  I turned on the Update per Screen option in the Velocity Layer Properties window, so that I could optimize the velocity scale for this view.  Here I noticed that the maximum velocity is 11.5 feet per second (fps) (which is about 3.5 meters per second).  This is way too high for this river, and in fact while hovering around the centerline of the river, I get maximum velocities of around 2.8 to 2.9 fps.  But since the Update per Screen option is turned on for velocity here, that indicates that there is a hotspot velocity of 11.5 fps somewhere in this view.  I just have to find it. 

On closer inspection, I can see some lighter shades of blue and green (i.e. higher velocities) in the floodplain adjacent to the lateral structures (as highlighted in the white circles).  Knowing that cells adjacent to lateral structures are typically going to be the culprits in 1D/2D models, I zoomed in to the cells around the lateral structures to get a closer look.  I started with the circled area to the right, since that was closest to the cross section that generated the numerical error of 5.154 ft.    



While zoomed in, it was hard to see at first, with the terrain turned on, so I turned off the terrain and noticed this little sliver of high velocity right at the boundary of the 1D river and the 2D area.


Turning the velocity off and the terrain back on, you can see how the 2D area dips into the main channel just slightly, and in fact it overlaps the lateral structure as well.  This reveals a very typical problem that you can run into with 1D/2D offline connections with lateral structures.  For a given timestep, RAS will compute a volume of water going over (or through in the case of gates and culverts) the lateral structure into a given cell.  But if the receiving cell has a small sliver of low lying area, like in this example, that relatively small amount of volume could be enough to significantly raise the water surface in the cell.  Perhaps enough so that it is even higher than the water surface adjacent to it in the 1D reach.  This would then send water back the other direction the next time step.  Besides the numerical shock of a sudden and large rise in stage, the oscillating effect of sending water back and forth can set up errors that persist, grow, and lead to an instability.

Now if you’re using the weir equation on the lateral structure, you could change your weir submergence decay exponent from the default value of 1 to 3 (1 is the most accurate, 3 is the most stable).  This has a dampening effect on the oscillations and errors you get from this situation.  Read more about the weir submergence decay exponent in the HEC-RAS User's Manual on page 8-41.  You might also reduce the weir coefficient.  This will reduce the volume of water that transfers from the river to the cell for a given timestep.  If there is no elevated feature represented by the lateral structure (e.g. it is not a levee), you would want to use a very low weir coefficient on the order of 0.1 to 0.5, as discussed in the 2D Modeling User's Manual on page 3-50.  However, in this case, it might be better to just use the 2D equations over the lateral structure instead of the weir equation.  Using the “Normal 2D Equation Domain” (this is just a funny way of saying “Use 2D Equations”) is a relatively new feature available in lateral structures.  But if your lateral structure does not represent an elevated terrain feature (e.g. a weir or levee), then this might be the better option to use. 



However, in recognizing that the 2D area perimeter and the lateral structure are poorly located, I will fix this first to see if that’s all that is needed to solve the problem. 
First, I’ll pull in the 2D area to just beyond the high ground.  This can easily be done in RAS Mapper on the 2D Area Perimeter Layer while in edit mode.  Next, I’ll pull over the lateral structure so that it resides ON the high ground. 


In the current version of HEC-RAS (as I type this post), Version 5.0.7, you cannot edit lateral structures in RAS Mapper.  So you have to do it in the Geometry Editor for now, using the Edit…Move Points/Objects command.  After moving the lateral structure, I adjusted the location of the 2D perimeter again, so it is just inside the lateral structure.  Here you can see a much better placement of both the 2D perimeter and lateral structure. 



And don’t forget, since I moved the lateral structure placement, I have to re-extract the terrain onto it.  Fortunately, HEC-RAS gives us a short cut to do this with the Terrain Profile button. 


After re-running, the error is gone and the results look much better.  Notice the peak velocity in the scale is back to a normal value of 2.6 fps.  And the 5.154-ft error is gone!





Online Class - Using HEC-RAS to Modeling Bridges, Culverts, and Floodplains

$
0
0
Hey RAS Users-

This is a basic level steady flow HEC-RAS class.  That's right, we start from the very beginning with HEC-RAS-steady flow theory, building 1D geometry, steady flow boundary conditions, and even introductions to RAS Mapper and 2D modeling.  I'll be hosting weekly live sessions for this class to review the course material and workshops and for Q&A.

And it's an online course-that means you can take it from anywhere in the world-just need a computer with a wifi connection.

If you're relatively new the HEC-RAS, or just need a refresher, come and join me for this excellent intro class on HEC-RAS modeling.  Don't delay, there are limited seats. 









 
 
 

The RAS Solution has moved to a address

Viewing all 200 articles
Browse latest View live