Preview of the new Finite-Volume Approach for 1D Reaches
One of the most anticipated new features soon to come in the next major version of HEC-RAS (Version 5.1) is the option of running unsteady 1D reaches with a finite volume solution scheme. This will be a fantastic addition to HEC-RAS. Gary Brunner recently gave me a brief overview of the new finite volume feature we can expect. But before you ask, there is no set release date for Version 5.1 yet. But I’m hoping we’ll see it within the next year or two.
1D Finite Volume Solution Algorithm
By Gary W. Brunner, P.E., D.WRE
Senior Technical Hydraulic Engineer
Hydrologic Engineering Center
A brand new solution algorithm has been developed for 1D modeling. A Finite-Volume solution approach, similar to what was added for 2D modeling will be available for 1D modeling in HEC-RAS version 5.1.
The current 1D Finite Difference solution scheme has the following deficiencies:
- Cannot handle starting or going dry in a cross section
- Low flow model stability issues with irregular cross section data
- Extremely rapidly rising hydrographs can be difficult to get stable
- Mixed flow regime (i.e. flow transitions) approach is approximate
- Stream junctions do not transfer momentum
The new 1D Finite Volume approach has the following positive attributes:
- Can start with cross sections completely dry, or they can go dry during a simulation (wetting/drying)
- Very stable for low flow modeling
- Can handle extremely rapidly rising hydrographs without going unstable
- Handles subcritical to supercritical flow, and hydraulic jumps better.
- Junction analysis is performed as a single 2D cell when connecting 1D reaches (continuity and momentum is conserved through the junction).
Additionally, the new 1D Finite Volume approached is solved in the same matrix as the 2D equations. Solving in the same matrix allows for faster 1D/2D model solutions and more accurate flow transfers between 1D and 2D elements. The equations are solved together and all hydraulic connections are updated together on an iteration by iteration approach, rather than separately, as in previous versions of HEC-RAS.