Stabilizing a Dynamic Unsteady HEC-RAS Model

Written by Chris Goodell
Copyright © RASModel.com. 2013.  All rights reserved.

One of the most frustrating aspects of unsteady HEC-RAS modeling can be the model stabilization process. You know, you’ve gone to great lengths collecting the best survey/topo data and solid hydrology. Then you’ve painstakingly spent hours…possibly days entered all of that data only to find that once you press the “Compute” button, the model crashes. The dreaded “Red Bar”! Sometimes you can get your simulation to complete without crashing, but the listed numerical errors are so high that you can’t with good conscience submit that as your final simulation. Either way, approaching an unsteady HEC-RAS model (especially a dynamic one) as a beginner with little experience and understanding of how to stabilize it can cause significant delays in your project and worse, completely blow up your budget. I’ve uploaded two HEC-RAS projects to the following Google Drive Site: https://drive.google.com/folderview?id=0B_s8OLJOgOi0YU92SGZQZm9raWM&usp=sharing



RawleyResUnstable.prj is an unsteady flow dam breach HEC-RAS project recently sent to me for help. Although the model ran to completion without crashing, it had unacceptably high errors.
RawleyResStable.prj is the fully stabilized version of the model with no numerical errors.

The following lists out the courses of action taken to stabilize the model. Feel free to download the “Unstable” and “Stable” models and try these techniques on your own. The links following some of these items will take you to more information about that particular technique.

1. Cross Section Spacing. The initial spacing was way too coarse. A visual check alone of the geometry schematic and profile plot should encourage you to investigate a finer cross section spacing.

Geometry Schematic

Profile Plot Samuels equation suggests anywhere from 15 ft to 50 ft spacing (depending on what bed slope you use). I interpolated to 50 ft for the entire reach. http://hecrasmodel.blogspot.com/2008/12/samuels-equation-for-cross-section.html

2. Dam Breach models typically have time steps on the order of a minute or less. This model was initially set with a computation interval of 10 minutes which is high even for the largest and “slowest” of dam breach models. I changed the time step to 10 seconds. The selection of 10 seconds was based on “gut” feel and lots of experience doing dam breach models. There are some methods for approximating good timesteps, notably the Courant Condition and Fread’s equation. Also, setting the cross section spacing/timestep ratio equal to a representative stream velocity will get you close. In this case, the cross section spacing is 50 ft, and I was able to pull some velocities (prior to the model crashing) at about 6 ft/s. That suggests a time step of 8.3 seconds. Close enough to 10 seconds, so we’ll stick with that.

3. I changed the downstream boundary normal depth slope from 0.01 to 0.001. Not sure what is downstream of the first cross section, but 0.01 is awfully steep and was setting up a very low depth at the downstream boundary (which was causing instabilities). http://hecrasmodel.blogspot.com/2010/01/downstream-boundary-normal-depth.html

Notice in the following profile plot of the downstream end of the reach how the water surface at the boundary cross section is below critical depth (the red dot). This creates an overestimation of the water depth at the next upstream cross section, which in turns creates some instability over the next several timesteps.

The proper way to handle this would be to find out what is downstream of your model and select a boundary condition that best represents those conditions. In absence of downstream conditions, the 0.001 slope for normal depth provides a reasonable solution.  Either way, this underscores the importance of moving your downstream boundary far away from your area of interest in your study reach.  That way the errors that do originate from your downstream boundary assumption will have diminished to negligible levels before impacting your area of interest. 

4. In the unsteady flow editor, the initial flow and the first time step flow should always be equal. In the original model, the initial flow is left blank (which is actually okay because RAS will use the first timestep flow if left blank). However, the first time step flow is very low at 0.14 cfs. Peak discharge of your inflow hydrograph is around 600 cfs. I put in a baseflow (minimum flow) of 5% of this, which is 30 cfs. Just to be safe, I put 30 cfs in the initial flow input box as well. It’s important to make sure that the 30 cfs baseflow does not have an impact on the peak of the breach outflow hydrograph. http://hecrasmodel.blogspot.com/2009/02/minimum-flow-requirements.html

5. The upstream end of your reservoir is very steep (14%). This one was easy to spot in the profile plot because of the very nice smooth profile plot followed by a sudden spike in the energy level (green dashed line) at the upstream end. Notice to that there is a corresponding supercritical solution at the upstream most cross section. These are both indications that more cross sections are needed here. http://hecrasmodel.blogspot.com/2010/01/another-reason-for-interpolated-cross.html

I added more cross sections by interpolating the steep slope at the upstream end of the reservoir (10 ft spacing).

Notice the energy spike is not completely gone, but it is much better. Refining the HTAB parameters (http://hecrasmodel.blogspot.com/2009/02/crazy-energy-grade-line.html), turning on mixed flow and bumping up n values will collectively take care of the rest of this (see numbers 8, 9, and 10 below).

7. From RS 4577 to 3791 you have your main channel defined in a small elevated side channel. The right overbank is lower in elevation than the channel.

This is problematic in RAS. Using the graphical cross section editor, I redefined the bank stations to get rid of this problem.

8. Notice now in the HTAB parameters property table, the HTAB definition begins at the elevation of the invert of the side channel that was previously defined as the main channel. There is a big gap in the computation points in the newly defined main channel. This poses problems with HEC-RAS, particularly at low discharges/stages, since RAS will have to extrapolate to obtain a solution over that range.

I readjusted the cross section HTAB parameters to reflect the new bank station definitions (I used the “Copy Invert” button for all cross sections). While in there, I maximized the resolution of the HTAB parameters by changing number of points to 100 and minimizing the increment as much as possible while still having full coverage of each cross section. http://hecrasmodel.blogspot.com/2009/02/crazy-energy-grade-line.html http://hecrasmodel.blogspot.com/2011/03/more-on-htab-parameters.html

9. In the figure below, notice all of the “red” areas in this zoomed in section of the  upstream end of the reach.

These are areas that are near critical depth. I turned on the “Mixed Flow” option in the Unsteady Flow Analysis window, which helps stabilize the near critical depths at the upstream end of the reservoir. http://hecrasmodel.blogspot.com/2011/04/mixed-flow-regime-options-lpi-method.html If you’ve been following along correcting your own copy of the model, you’ll notice that the model now finally runs to completion without crashing. There are still a few minor errors as shown in the Computation Messages.

It’s debatable whether you need to take care of these relatively small errors. Particularly for a dam breach model, where the presumed errors of the input data probably far outweigh these small numerical errors. Nevertheless, there’s a sense of pride in putting out a very robust model, completely free of numerical errors. So…let’s continue.

10. The errors in the computation window above suggest a problem around river station 6200. That happens to be the upstream end of the reservoir and the source of the error can clearly be seen with the energy grade line and critical depth turned on.

The section at the upstream end of the reservoir is very steep at about 13.5%. In fact, that’s greater than the suggested maximum slope of 10%, as stated in the HEC-RAS manuals. You can see as the pool starts to lower, the steep reach is exposed and because of the low n values, the water surface is calculated to be supercritical. Let’s assume that the geometry is correct, and the upstream end of the reservoir really is that steep. Jarrett’s Equation suggests very high n values should be used here (~0.2), based on the bed slope and hydraulic radius (after the breach when the reservoir has drained). The original model had 0.07 and 0.035 for the overbanks and main channel, respectively. I changed all the n values in this steep reach to 0.15, because I think 0.2 might be a little overkill. http://hecrasmodel.blogspot.com/2009/12/n-values-in-steep-streams.html

There…no errors! However, there is still a Warning about extrapolating above/beyond the rating curve at a bridge (R.S. 277).

11. To get rid of that warning, I increased the HTAB headwater maximum elevation at River Station 277 (the bridge) from 278 ft to 280 ft so that RAS doesn’t have to extrapolate.

There we have it. A clean solution. No errors, no warnings.

Initial Conditions Trick

Written by Chris Goodell, P.E., D. WRE | WEST Consultants
Copyright © RASModel.com. 2011. All rights reserved.
I’ve mentioned many times on this blog and in class about the importance of having your initial conditions flow match your first time step flow, when running an unsteady flow model.  As with any “rule” in RAS modeling, there are exceptions to this, but generally speaking, if your first timestep flow is “X”, then your initial conditions flow should also be “X”. 

This is generally very easy to do-if you have a single reach.  The hard part is just remembering to do it.  So, if you know you’re going to be changing your initial flow a lot, maybe while testing different hydrologic scenarios, or just trying to stabilize your model, you can save yourself some time and extra mouse clicks by leaving the initial flow cell blank.  By leaving it blank, you’re telling RAS to just use the first timestep flow provided in the flow hydrograph.  Give it a try, it works great.

This works perfectly for single reach systems.  For more complex systems, there is a catch.  For dendritic systems, this trick only applies to the upstream end of Rivers, not at junctions.  For looped systems, this trick only applies to the upstream end of Rivers that do not originate as distributaries from a junction.  Don’t worry, if you fail to enter an initial flow value where one is required, RAS will tell you about this when you try to compute. 

How to Create a Hotstart File in HEC-RAS for Dam Breach Analysis

Written by Aaron A. Lee   | WEST Consultants
Copyright © RASModel.com. 2010. All rights reserved.
While running unsteady flow simulations in HEC-RAS instabilities may occur when transitioning from the automatically created initial condition file to the first computed time step. These instabilities can be caused by mixed flow conditions, flow splits, or poorly defined initial conditions. A hotstart is another option available for defining initial conditions for the project model. This article presents one technique for setting up a hotstart run to help with initial conditions problems and to troubleshoot problem areas in your project model.

This is done by creating a new plan, using a flow file with a constant discharge as the upstream boundary, and a stage hydrograph as the downstream boundary over a 24 hour period. Typically 24 hours is long enough, but you may find that a longer hotstart period is required. The downstream boundary water surface elevation is defined artificially high, and over the simulation it is gradually reduced until it reaches the true starting depth for your project model. At this point the hotstart file will be written by HEC-RAS. In HEC-RAS lingo, the term “hotstart file” is used interchangeably with “Restart File” and “Initial Conditions file”.

STEP 1. Create a new flow file by opening the current unsteady flow file and navigate to File, Save As, and name it “hotstart”. Once saved, change the upstream boundary condition by selecting Flow Hydrograph under the Boundary Conditions tab. Change the discharge to a constant flow equal to that of the first timestep for the 24-hr period. In this example, shown in Figure 1, the entire Flow column should be modified to contain 155 for the full 24-hr simulation.



Select Stage Hydrograph as the downstream boundary condition. The beginning stage should start at an artificially high elevation - somewhere near the invert of the upstream-most cross section. The final elevation at the end of the 24-hr simulation should be equal to the starting downstream water surface elevation of the project plan. Once these values are added into the Stage column, use the Interpolate Missing Values button to add the missing elevations.

The elevation at each time interval should gradually decrease. This makes it easier for the modeler to observe problems as they occur, and where they happen in the model. For this example, the starting elevation is 820 ft and the final elevation is 637.25 ft. This is shown in Figure 2. Save the “hotstart” flow file.

STEP 2. Navigate to the Unsteady Flow Analysis window and save as a new plan named “hotstart plan”. Name the short ID as “hotstart”. This will create a new plan that will be used to define the initial conditions for the project plan. Under the Unsteady Flow Analysis window select the Ending Date and Time to 24 hours after the Starting Date and Time (or whatever time frame you want to use-just be sure it is consistent with the hotstart flow hydrograph and stage hydrograph you created in Step 1).

STEP 3. The next step is to set up the model so that the hotstart file will be written. Under the Unsteady Flow Analysis window navigate to Options, then Output Options. Figure 3 shows this window.



Check the boxes that write the initial conditions file at the Fixed Reference of the hotstart simulation ending date and time. At the end of the specified simulation time (24 hours) HEC-RAS will automatically write the initial condition file. The final step is to ensure that the hotstart plan is using the correct geometry file, and the created “hotstart” flow file. Once the plan is completed and saved, compute the Hotstart simulation.

STEP 4. The profile plot should be reviewed for problems with the hotstart model. This will indicate areas that may cause problems in your project model. Over the course of the hotstart simulation, as the water surface drops into place along your bed profile, look for hints of instabilities. If you hotstart simulation crashes, you’ll know exactly where to investigate-the intersection between the artificially high horizontal pool and the bed profile at the time of the crash. To use the hotstart file as your initial conditions, go to the Unsteady Flow editor of your project plan. Click on the “Use a Restart File” box and browse for the initial conditions file. This file will have an extension that indicates that hotstart plan number that created it, the simulation time (day-month-year) when it was created, and .rst. In this example, it should look like:

Workshop4.p03.10NOV2006.rst

*Warning-If you make a change to the geometry in your project model, you’ll have to re-run your hotstart simulation. However, once everything is set up, this is very easy to do. Simply open the hotstart plan and run it. Then open the project plan and run it.