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.
And a zoomed in view:
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.
Good work Chris! Thanks for sharing!
ReplyDeleteIt it possible to model even smaller (narrow) fishways with smaller boulders (weirs) between little ponds (gaps between boulders?
Yes, in theory you can go smaller and smaller. However once the water depth to cell size ratio gets too high, you can get numerical dispersion issues. This is also referred to as the "Abbott Number" (named after famed IHE-Delft hydro-informatics professor Michael Abbott). Page 10-166 of the Austrailian Rainfall & Runoff Project 15 report summarizes this well. http://arr.ga.gov.au/__data/assets/pdf_file/0019/40573/ARR_Project15_TwoDimensional_Modelling_DraftReport.pdf
DeleteThanks Chris!
DeleteI'll read about it!
This is a particularly interesting, and seemingly novel, application of HEC-RAS - I am excited that it is becoming a general use quasi-CFD tool. One question - what is the thinking behind the inclusion of the eddy mixing term in a model with such fine grid resolution? Is it to compensate for the lack of three dimensionality in the flow field?
ReplyDeleteNo it is more to capture the diffusion and energy loss associated with turbulence in two dimensions.
DeleteThis is a particularly interesting, and seemingly novel, application of HEC-RAS - I am excited that it is becoming a general use quasi-CFD tool. One question - what is the thinking behind the inclusion of the eddy mixing term in a model with such fine grid resolution? Is it to compensate for the lack of three dimensionality in the flow field?
ReplyDeleteThis is a particularly interesting, and seemingly novel, application of HEC-RAS - I am excited that it is becoming a general use quasi-CFD tool. One question - what is the thinking behind the inclusion of the eddy mixing term in a model with such fine grid resolution? Is it to compensate for the lack of three dimensionality in the flow field?
ReplyDeleteHi Chris,
ReplyDeleteI would like to model a 2D flood simulation, it's possible to apply Precipitation for the Rainfall Runoff simulation? What is your assumption and suggestion to model 2D flood simulation by applying Rainfall as the Unsteady boyndary condition?
Kind Regards,
This is interesting... It could be also interesting to know something more about the eddy viscosity mixing coefficient. Why the default is blank (=zero). When is critical to use the eddy viscosity? If is not use, the simulation gives greater depth value (or velocity)? E.g. for a Dam break simulation, for the sake of safety, it could be better to not use it?
ReplyDeleteHi Chris! What do you think about a 2d model of hydraulic structures like a over crest spillway? I need to available water profile, velocities, where ocours hydraulic jump.
ReplyDeleteWhat do you recommend for this?
Thanks!
Chris, Ive never seen that document before. That was a great read that every modeler should be familiar with
ReplyDeleteChris, Ive never seen that document before but it was a greattttt read that every modeler should read through
ReplyDeleteHi Chris. How do you create a mesh with varying cell sizes, as can be seen in the above images? I'm aware of break lines (with their associated cell spacing options) but it looks like you've defined cells to be smaller across a polygon area. I figure I'm missing something...
ReplyDelete