HEC-RAS Applications Guide-V1-20221124 - 133856
HEC-RAS Applications Guide-V1-20221124 - 133856
– 1
HEC-RAS Applications Guide – HEC-RAS Applications Guide
US Army Corps
of Engineers
Hydrologic Engineering Center
HEC-RAS
River Analysis System
Applications Guide
Version 6.0
May 2021
Approved for Public Release. Distribution Unlimited
CPD-70 | REPORT DOCUMENTATION PAGE | OMB No. 0704-0188
Public reporting burden for this collection of information is estimated to average 1 hour per
response, including the time for reviewing instructions, searching existing data sources, gathering
and maintaining the date needed, and completing and reviewing the collection of information.
Send comments regarding this burden estimate or any other aspect of this collection of
information, including suggestions for reducing this burden, to Washington Headquarters Services,
Directorate for Information Operations and Reports, 1215 Jefferson Davis Highway, Suite 1204,
Arlington, VA 22202-4302, and to the Office of Management and Budget, Paperwork Reduction
Project (0704-0188), Washington, DC 20503.
6. AUTHOR(S)
Gary W. Brunner, John C. Warner, Brent C. Wolfe, Steven S. Piper, and
Landon Marston
9. SPONSORING / 10.
MONITORING SPONSO
AGENCY NAME(S) RING /
AND ADDRESS(ES) MONITO
RING
AGENCY
REPORT
NUMBER
11.
SUPPLEMENTARY
NOTES
HEC-RAS is an integrated system of software, designed for interactive use in a multi-tasking, multi-
user network environment. The system is comprised of a graphical user interface (GUI), separate
hydraulic analysis components, data storage and management capabilities, graphics and reporting
facilities.
The HEC-RAS system contains four one-dimensional hydraulic analysis components for: (1) steady
flow water surface profile computations; (2) unsteady flow simulation; (3) movable boundary
sediment transport computations; and (4) temperature and water quality constituent transport
modeling. A key element is that all four components use a common geometric data representation
and common geometric and hydraulic computation routines. In addition to the four hydraulic
analysis components, the system contains several hydraulic design features that can be invoked
once the basic water surface profiles are computed.
16. PRICE
CODE
1 http://www.hec.usace.army.mil/
The United States Government, US Army Corps of Engineers, Hydrologic Engineering Center ("HEC") grants to the
user the rights to install the HEC River Analysis System (HEC-RAS) "the Software" (either from a disk copy obtained
from HEC, a distributor or another user or by downloading it from a network) and to use, copy and/or distribute
copies of the Software to other users, subject to the following Terms and Conditions for Use:
All copies of the Software received or reproduced by or for user pursuant to the authority of this Terms and
Conditions for Use will be and remain the property of HEC.
User may reproduce and distribute the Software provided that the recipient agrees to the Terms and Conditions for
Use noted herein.
HEC is solely responsible for the content of the Software. The Software may not be modified, abridged, decompiled,
disassembled, unobfuscated or reverse engineered. The user is solely responsible for the content, interactions, and
effects of any and all amendments, if present, whether they be extension modules, language resource bundles,
scripts or any other amendment.
The name "HEC-RAS" must not be used to endorse or promote products derived from the Software. Products
derived from the Software may not be called "HEC- RAS " nor may any part of the "HEC- RAS " name appear within
the name of derived products.
No part of this Terms and Conditions for Use may be modified, deleted or obliterated from the Software.
No part of the Software may be exported or re-exported in contravention of U.S. export laws or regulations.
Waiver of Warranty:
THE UNITED STATES GOVERNMENT AND ITS AGENCIES, OFFICIALS, REPRESENTATIVES, AND EMPLOYEES,
INCLUDING ITS CONTRACTORS AND SUPPLIERS PROVIDE HEC-WAT \"AS IS,\" WITHOUT ANY WARRANTY OR
CONDITION, EXPRESS, IMPLIED OR STATUTORY, AND SPECIFICALLY DISCLAIM ANY IMPLIED WARRANTIES OF TITLE,
MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND NON-INFRINGEMENT. Depending on state law, the
foregoing disclaimer may not apply to you, and you may also have other legal rights that vary from state to state.
Limitation of Liability:
IN NO EVENT SHALL THE UNITED STATES GOVERNMENT AND ITS AGENCIES, OFFICIALS, REPRESENTATIVES, AND
EMPLOYEES, INCLUDING ITS CONTRACTORS AND SUPPLIERS, BE LIABLE FOR LOST PROFITS OR ANY SPECIAL,
INCIDENTAL OR CONSEQUENTIAL DAMAGES ARISING OUT OF OR IN CONNECTION WITH USE OF HEC-WAT
REGARDLESS OF CAUSE, INCLUDING NEGLIGENCE.
THE UNITED STATES GOVERNMENT'S LIABILITY, AND THE LIABILITY OF ITS AGENCIES, OFFICIALS,
REPRESENTATIVES, AND EMPLOYEES, INCLUDING ITS CONTRACTORS AND SUPPLIERS, TO YOU OR ANY THIRD
PARTIES IN ANY CIRCUMSTANCE IS LIMITED TO THE REPLACEMENT OF CERTIFIED COPIES OF HEC-RAS WITH
IDENTIFIED ERRORS CORRECTED. Depending on state law, the above limitation or exclusion may not apply to you.
Indemnity:
As a voluntary user of HEC- RAS you agree to indemnify and hold the United States Government, and its agencies,
officials, representatives, and employees, including its contractors and suppliers, harmless from any claim or
demand, including reasonable attorneys' fees, made by any third party due to or arising out of your use of HEC- RAS
or breach of this Agreement or your violation of any law or the rights of a third party.
Assent:
By using this program you voluntarily accept these terms and conditions. If you do not agree to these terms and
conditions, uninstall the program and return any program materials to HEC (if you downloaded the program and do
not have disk media, please delete all copies, and cease using the program.)
2 Foreword
The U.S. Army Corps of Engineers' River Analysis System (HEC-RAS) is software that allows you to perform one-
dimensional steady flow hydraulics; one and two-dimensional unsteady flow river hydraulics calculations; quasi
Unsteady and full unsteady flow sediment transport-mobile bed modeling; water temperature analysis; and
generalized water quality modeling (nutrient fate and transport).
The first version of HEC-RAS (version 1.0) was released in July of 1995. Since that time there have been several
major releases of this software package, including versions: 1.1; 1.2; 2.0; 2.1; 2.2; 3.0, 3.1, 3.1.1, 3.1.2, 3.1.3, 4.0, 4.1,
5.0 and now version 6.0 in 2021.
The HEC-RAS software was developed at the Hydrologic Engineering Center (HEC), which is a division of the
Institute for Water Resources (IWR), U.S. Army Corps of Engineers.
The software was designed by Mr. Gary W. Brunner, leader of the HEC-RAS development team. The user interface
and graphics were programmed by Mr. Mark R. Jensen, Alex Kennedy, Anton Rotter-Sieren, Cameron Ackerman, and
Stanford Gibson. The steady flow water surface profiles computational module and the majority of the one-
dimensional unsteady flow computations modules was programmed by Mr. Steven S. Piper. The One-dimensional
unsteady flow matrix solution algorithm was developed by Dr. Robert L. Barkau (Author of UNET and HEC-UNET).
The two-dimensional unsteady flow modeling capabilities were developed by Gary W. Brunner, Ben Chacon
(Resource Management Consultants, RMA), Steve S. Piper, Mark R. Jensen, Alex J. Kennedy, and Alex Sanchez.
The sediment transport interface module was programmed by Mr. Stanford A. Gibson. The quasi unsteady flow
computational sediment transport capabilities were developed by Stanford A. Gibson and Steven S. Piper. The
Unsteady flow sediment transport modules were developed by Stanford A. Gibson, Steven S. Piper, and Ben Chacon
(RMA). Special thanks to Mr. Tony Thomas (Author of HEC-6 and HEC-6T) for his assistance in developing the quasi-
unsteady flow sediment transport routines used in HEC-RAS. The two-dimensional sediment transport modules
were developed by Alex Sanchez and Stanford Gibson. The Debris flow capabilities in HEC-RAS (1D and 2D) were
developed by Stanford Gibson and Alex Sanchez. Most of the sediment output was designed by Stanford Gibson
and Alex Sanchez and programmed by Anton Rotter-Sieren.
The new 2D plotting library and plots (Breach Plot, Hydrographs, and DSS viewer) were developed by Mark R.
Jensen, Anton Rotter-Sieren, and Ryan Miles (RMA).
The new 3D visualization tool was developed by Anton Rotter-Sieren and Alex Kennedy.
The water quality computational modules were designed and developed by Mr. Mark R. Jensen, Dr. Cindy Lowney
and Zhonglong Zhang (ERDC-RDE-EL-MS).
The spatial data and mapping tools (RAS Mapper) were developed by Mark R. Jensen, Cameron T. Ackerman, Alex J.
Kennedy, and Anton Rotter-Sieren. Special thanks to Mr. Will Breikreutz for his assistance in developing the RAS Tile
server.
The interface for channel design/modifications was designed and developed by Mr. Cameron T. Ackerman and Mr.
Mark R. Jensen. The stable channel design functions were programmed by Mr. Chris R. Goodell.
The routines that import HEC-2 and UNET data were developed by Ms. Joan Klipsch. The routines for modeling ice
cover and wide river ice jams were developed by Mr. Steven F. Daly of the Cold Regions Research and Engineering
Laboratory (CRREL).
Many other HEC staff members have made contributions in the development of this software, including: Vern R.
Bonner, Richard Hayes, John Peters, Al Montalvo, and Michael Gee. Mr. Matt Fleming was the Chief of the H&H
Foreword – 6
HEC-RAS Applications Guide – HEC-RAS Applications Guide
Division, and Mr. Chris Dunn was the director during the development of this version of the software.
This manual was written by John C. Warner, Gary W. Brunner, Brent C. Wolfe, Steven S. Piper, and Landon Marston.
Foreword – 7
HEC-RAS Applications Guide – HEC-RAS Applications Guide
3 Introduction
Welcome to the U.S. Army Corps of Engineers River Analysis System (HEC-RAS) developed by the Hydrologic
Engineering Center. This software allows you to perform one-dimensional steady flow, one and two-dimensional
unsteady flow calculations, sediment transport/mobile bed computations, and water temperature/water quality
modeling.
The HEC-RAS modeling system was developed as a part of the Hydrologic Engineering Center's "Next
Generation" (NexGen) of hydrologic engineering software. The NexGen project encompasses several aspects of
hydrologic engineering, including: rainfall-runoff analysis (HEC-HMS); river hydraulics (HEC-RAS); reservoir system
simulation (HEC-ResSim); flood damage analysis (HEC-FDA and HEC-FIA); and real-time river forecasting for
reservoir operations (CWMS).
Documentation Description
User's Manual This manual is a guide to using HEC-RAS. The manual provides an introduction and
overview of the modeling system, installation instructions, how to get started, a
simple example, entering and editing geometric data, detailed descriptions of each of
the major modeling components, and how to view graphical and tabular output.
2D User's Manual This document describes how to use the 2D modeling capabilities that are included in
this version of the software. It also describes how to use RAS Mapper in support of 2D
modeling (mesh generation) and inundation mapping for models containing 2D flow
areas.
HEC-RAS Mapper This document describes how to use HEC-RAS Mapper to do the following: establish a
horizontal coordinate system; develop an HEC-RAS terrain model; layout the
geometric data model and extract terrain data; visualize results in the form of maps,
plots, and tables.
Sediment Transport This manual describes how to perform sediment transport modeling. The document
User's Manual describes 1D quasi unsteady; 1D unsteady flow, and 2D sediment transport modeling.
Additionally sediment impact analysis (SIAM) and bank stability using BSTEM is also
described.
Hydraulic Reference This manual describes the theory and data requirements for the hydraulic
Manual calculations performed by HEC-RAS. Equations are presented along with the
assumptions used in their derivation. Discussions are provided on how to estimate
model parameters, as well as guidelines on various modeling approaches.
Introduction – 8
HEC-RAS Applications Guide – HEC-RAS Applications Guide
Applications Guide This document contains a series of examples that demonstrate various aspects of
HEC-RAS. Each example consists of a problem statement, data requirements, and
general outline of solution steps, displays of key input and output screens, and
discussions of important modeling aspects.
Introduction – 9
HEC-RAS Applications Guide – HEC-RAS Applications Guide
subcritical and supercritical flow regime analyses are presented to show inconsistencies that developed,
and to provide guidance when to perform a mixed flow analysis.
• Example 10, Stream Junction, demonstrates the analysis of a river system that contains a junction. This
example illustrates a flow combining of two subcritical streams, and both the energy and momentum
methods are used for two separate analyses.
• Example 11, Bridge Scour, presents the determination of a bridge scour analysis. The user should be
familiar with the procedures for modeling bridges before reviewing this example. The scour equations and
procedures are based upon the methods outlined in Hydraulic Engineering Circular No. 18 (FHWA 1995).
• Example 12, Inline Weir and Gated Spillway, demonstrates the analysis of a river reach that contains an
inline weir and a gated spillway. Procedures for entering the data to provide flexibility for the flow analysis
are provided.
• Example 13, Bogue Chitto - Single Bridge (WSPRO), performs an analysis of a river reach that contains a
bridge crossing. The example is similar to Example 2, however, all of the water surface profiles are low flow
and are computed using the WSPRO (FHWA, 1990) routines that have been adapted to the HEC-RAS
methodology of cross section locations around and through a bridge.
• Example 14, Ice-Covered River, is an example of how to model an ice covered river as well as a river ice-
jam.
• Example 15, Split Flow Junction With Lateral Weir and Spillway, is an example of how to perform a split
flow optimization with the steady flow analysis portion of the software. This example has a split of flow at a
junction, as well as a lateral weir.
• Example 16, Channel Modification. This example demonstrates how to use the channel modification
feature within the HEC-RAS Geometric Data Editor. Channel modifications are performed, and existing and
modified conditions geometry and output are compared.
• Example 17, Unsteady Flow Application. This example demonstrates how to perform an unsteady flow
analysis with HEC-RAS. Discussions include: entering storage area information; hydraulic connections;
unsteady flow data (boundary conditions and initial conditions); performing the computations; and
reviewing the unsteady flow results.
• Example 18, Advanced Inline Structure Modeling. This example demonstrates all of the inline structure
capabilities within a single structure. This includes: overflow weirs; gates; culverts; user defined rating
curves; and a flow time series outlet.
• Example 19, Hydrologic Unsteady Flow Routing – Modified Puls. This example demonstrates how to
intermix hydrologic routing in the middle of an unsteady flow model. Hydrologic routing can be used for very
steep reaches where the full unsteady flow equations may not be applicable or stable.
• Example 20, Hager's Lateral Weir Equation. This example demonstrates how to use Hager's Lateral Weir
equation for a Lateral Structure.
• Example 21, Overflow Gates. This example demonstrates how to use overflow gates in HEC-RAS.
• Example 22, Groundwater Interflow. This example demonstrates how to use the groundwater interflow
options in HEC-RAS, within a river reach or with storage areas.
• Example 23, Urbane Modeling. This example demonstrates how to model pressurized pipe systems within
the HEC-RAS unsteady flow modeling capabilities.
• Example 24, Manning's n Calibration. This example discusses how to calibrate an unsteady flow model for
base Manning's n values, as well as how to use the HEC-RAS automated Manning's n value calibration
option.
• Appendix A contains a list of references.
Introduction – 10
HEC-RAS Applications Guide – HEC-RAS Applications Guide
The cross section data were entered in the Cross Section Data Editor, which is activated by selecting the Cross
Section icon on the Geometric Data Editor (as outlined in Chapter 6 of the User's Manual). Most of the 12 cross
sections contain at least 50 pairs of X-Y coordinates, so the cross section data will not be shown here for brevity. The
distances between the cross sections are as shown in the bottom figure below (The reach lengths for cross section
12 can be seen by using the scroll bars in the window.). This summary table can be viewed by selecting Tables and
then Reach Lengths on the Geometric Data Editor.
From the geometric data, it can be seen that most of the cross sections are spaced approximately 500 feet apart.
The change in elevation from cross section 12 to cross-section 1 is approximately 56 feet along the river reach of
5700 feet. This yields a slope of approximately 0.01 ft/ft, which can be considered as a fairly steep slope. The
remaining geometric data consists of Manning's n values of 0.10, 0.04, and 0.10 in the left overbank (LOB), main
channel, and right overbank (ROB), respectively. Also, the coefficients of contraction and expansion are 0.10 and
0.30, respectively. After all the geometric data was entered, it was saved as the file "Base Geometry Data."
Next, the Reach Boundary Conditions button located at the top of the Steady Flow Data Editor was selected. The
reach was analyzed for subcritical flow with a downstream normal depth boundary condition of S = 0.01 ft/ft. This
value was estimated as the average slope of the channel near the downstream boundary. For a subcritical flow
analysis, boundary conditions must be set at the downstream end(s) of the river system. After all of the flow data
was entered, it was saved as the file "100 Year Profile."
From this profile, it can be seen that the water surface appears to approach or is equal to the critical depth at
several locations. For example, from section 12 through 8, the water surface appears to coincide with the critical
depth. This implies that the program may have had some difficulty in determining a subcritical flow value in this
region, or perhaps the actual value of the flow depth is in the supercritical flow regime. To investigate this further, a
closer review of the output needs to be performed. This can be accomplished by reviewing the output at each of the
cross sections in either graphical or tabular form, and by viewing the summary of Errors, Warnings and Notes.
First, a review of the output at each cross section will be performed. From the main program window, select View,
Detailed Output Tables, Type, and then Cross Section. Selection of cross section 12 should result in the display as
shown in the figure below. At the bottom of the table is a box that displays any errors, warnings, or notes that are
specific to that cross section. For this example, there are several warning messages at cross section 12. The first
warning is that the velocity head has changed by more than 0.5 feet and that this may indicate the need for
additional cross sections. To explain this message, it is important to remember that for a subcritical flow analysis,
the program starts at the downstream end of the reach and works upstream. After the program computed the water
surface elevation for the 11th cross section, it moved to the 12th cross section. When the program computed the
water surface elevation for the 12th cross section, the difference in the velocity head from the 11th to the 12th cross
section was greater than 0.5 feet. This implies that there was a significant change in the average velocity from
section 11 to section 12. This change in velocity could be reflecting the fact that the shape of the cross section is
changing dramatically and causing the flow area to be contracting or expanding, or that a significant change in
slope occurred. In order to model this change more effectively, additional cross sections should be supplied in the
region of the contraction or expansion. This will allow the program to better calculate the energy losses in this
region and compute a more accurate water surface profile.
The second warning at cross section 12 states that the energy loss was greater than 1.0 feet between the current
cross section (#12) and the previous cross section (#11). This warning also indicates the possible need for additional
cross sections. This is due to the fact that the rate of energy loss is usually not linear. However, the program uses, as
a default, an average conveyance equation to determine the energy losses. Therefore, if the cross sections are too
far apart, an appropriate energy loss will not be determined between the two cross sections. (The user may select
alternate methods to compute the average friction slope. Further discussion of user specified friction loss
formulation is discussed in Chapter 4 of the Hydraulic Reference Manual.)
A review of other cross sections reveals the same and additional warnings. To review the errors, notes, and
warnings for all of the cross sections, select Summary Errors, Warnings, and Notes from the View menu on the
main program window. A portion of the summary table is shown in the figure below.
The additional warnings and notes that are listed in the summary table are described as follows.
Warning - The energy equation could not be balanced within the specified number of iterations. The program used
critical depth for the water surface and continued on with the calculations.
This warning implies that during the computation of the upstream water surface elevation, the program could not
compute enough energy losses to provide for a subcritical flow depth at the upstream cross section. Therefore, the
program defaulted to critical depth and continued on with the analysis.
Warning - The energy equation could not be balanced within the specified number of iterations. The program used
critical depth for the water surface and continued on with the calculations.
After the flow depth was calculated for the cross section, the program determined that the flow was occurring in
more than one portion of the cross section. For example, this warning occurred at river station # 10 and the plot of
this cross section is shown in the figure below. From the figure, it can be seen that at approximately an X-coordinate
of 800, there exists a large vertical land mass. During this output analysis, it must be determined whether or not the
water can actually be flowing on both sides of the land mass at this flow rate. Since the main channel is on the right
side of the central land mass, could the water be flowing on the left side or should all of the flow be contained to the
right side of the land mass? By default, the program will consider that the water can flow on both sides of the land
mass. If this is not correct, then the modeler needs to take additional action.
Additional action can be one of two procedures. First, if the existing scenario is not feasible, then the water on the
left side may be considered as an ineffective flow area, where the water is accounted for volumetrically but it is not
considered in the conveyance determination until a maximum elevation is reached. Secondly, if all of the flow
should be occurring only on the right side of the land mass, then the land mass could be considered as a levee. By
defining the central vertical land mass as a levee, the program will not permit a flow onto the left side of the levee
until the flow depth overtops the levee.
For further discussion on ineffective flow areas and levees, refer to Chapter 6 of the User's Manual and Chapter 3 of
the Hydraulic Reference Manual.
• Warning - During the standard step iterations, when the assumed water surface was set equal to critical
depth, the calculated water surface came back below critical depth. This indicates that there is not a valid
subcritical answer. The program defaulted to critical depth. This warning is issued when a subcritical flow
analysis is being performed but the program could not determine a subcritical flow depth at the specified
cross section. As the program is attempting to determine the upstream depth, it is using an iterative
technique to solve the energy equation. During the iterations, the program tried critical depth as a possible
solution, which resulted in a flow depth less than critical. Since this is not possible in a subcritical analysis,
the program defaulted to using critical depth at this cross section and continued on with the analysis. This
error is often associated with too long of a reach length between cross sections or misrepresentation of the
effective flow area of the cross section.
• Warning - The parabolic search method failed to converge on critical depth. The program will try the cross
section slice/secant method to find critical depth. This message appears if the program was required to
calculate the critical depth and had difficulty in determining the critical depth at the cross section. The
program has two methods for determining critical depth: a parabolic method and a secant method. The
parabolic method is the default method (this can be changed by the user) because this method is faster and
most cross sections have only one minimum energy point. However, for cross sections with large, flat over
banks, there can exist more than one minimum energy point. For further discussion, refer to the section
Critical Depth Determination in Chapter 2 of the Hydraulic Reference Manual.
• Note - Multiple critical depths were found at this location. The critical depth with the lowest, valid, water
surface was used. This note appears when the program was required to determine the critical depth and
accompanies the use of the secant method in the determination of the critical depth (as described in the
previous warning message). This note prompts the user to examine closer the critical depth that was
determined to ensure that the program supplied a valid answer. For further discussion, refer to the section
Critical Depth Determination in Chapter 2 of the Hydraulic Reference Manual.
• Warning - The conveyance ratio (upstream conveyance divided by downstream conveyance) is less than 0.7
or greater than 1.4. This may indicate the need for additional cross sections. The conveyance of the cross
section, K, is defined by:
If the n values for two subsequent cross sections are approximately the same, it can be seen that the ratio of the
two conveyances is primarily a function of the cross sectional area. If this ratio differs by more than 30%, then this
warning will be issued. This warning implies that the cross sectional areas are changing dramatically between the
two sections and additional cross sections should be supplied for the program to be able to more accurately
compute the water surface elevation.
In summary, these warnings and notes are intended to inform the user that potential problems may exist at the
specified cross sections. It is important to note that the user does not have to eliminate all the warning messages.
However, it is up to the user to determine whether or not these warnings require additional action for the analysis.
The interpolation was adequate for the right overbank and the main channel. However, the interpolation in the left
overbank failed to connect the two existing high ground areas. These two high ground areas could be representing
a levee or some natural existing feature. Therefore, Del Interp was selected to delete the interpolation. (This only
deleted the interpolation between cross sections 11 and 10.) Then, the two high points and the low points of the
high ground areas were connected with user supplied master cords. This was accomplished by selecting the Master
Cord button and connecting the points where the master cords should be located. Finally, a maximum distance of
150 feet was entered between cross sections and Interpolate was selected. The final interpolation appeared as is
shown in the following figure.
The modeler should now go through all of the interpolated cross sections and determine that the interpolation
procedure adequately produced cross sections that depict the actual geometry. When completed, the geometric
data was saved as the new file name "Base Geometry + Interpolated." This allowed the original data to be unaltered
To investigate this further, the results will be viewed in tabular form. From the main program window, View, Profile
Summary Tables, Std. Tables, and then Standard Table 1 were selected. This table for the mixed flow analysis is
shown in the Standard Table 1 (shown below). The table columns show the default settings of river, reach, river
station, total flow, minimum channel elevation, water surface elevation, etc. The meanings of the headings are
described in a box at the bottom of the table. By selecting a cell in any column, the definition of the heading will
appear in the box for that column.
From the Standard Table 1, the water surface elevations and critical water surface elevations can be compared.
The values at river station 11.2* show that the flow is supercritical at this cross section since the water surface is at
an elevation of 1811.29 ft and the critical water surface elevation is 1811.46 ft. Additionally, it can be seen that the
flow at river station 11.0 is subcritical. (Note: the asterisks "*" denote that the cross sections were interpolated.) By
selecting the Cross Section type table (as performed above), and toggling to river station 11.0, a note appears at
the bottom of the table indicating that a hydraulic jump occurred between this cross section and the previous
upstream cross section. These results are showing that the flow is both subcritical and supercritical in this reach.
The user can continue this process of reviewing the warnings, notes, profile plot, profile tables, and cross section
tables to determine if additional cross sections are required.
The reach was initially defined with 14 cross sections beginning at river mile 5.00 as the downstream river station
and river mile 5.99 as the upstream river station. The cross sections with an asterisk "*" were added by
interpolation for the purposes of this example. When the bridge was added, it was placed at river mile 5.40 to place
it at the appropriate location. On the river schematic, some of the cross section labels may not appear due to
overlapping of the labels. If this occurs, the labels can be seen by zooming in on the location of the closely spaced
cross sections.
Reach Lengths. The distances between the cross sections are entered as the downstream reach lengths in the
Cross Section Data Editor. To view the summary of the reach lengths, the table as shown in the figure below can be
activated by selecting Tables and then Reach Lengths from the Geometric Data Edito. The reach lengths were
obtained by measuring the distances on the USGS atlas. To determine the main channel distances, it was initially
assumed that during the peak event, the major active portion of the flow will follow the course of the main channel.
If, after the analysis, it is determined that the major portion of the active flow is not following the main channel
course, then the main channel flow distances will need to be adjusted. In other words, if the major portion of the
active flow is "cutting across" the meanders of the main channel, then these reach lengths will need to be
reevaluated.
The reach lengths determine the placement of the cross sections. The placement of the cross sections relative to
the location of the bridge is crucial for accurate prediction of expansion and contraction losses. The bridge routine
utilizes four cross sections to determine the energy losses through the bridge. (Additionally the program will
interpret two cross sections inside of the bridge by superimposing the bridge data onto both the immediate
downstream and upstream cross sections from the bridge.) The following is a brief summary for the initial
estimation of the placement of the four cross sections. The modeler should review the discussion in Chapter 6 of the
User's Manual and Chapter 5 of the Hydraulic Reference Manual for further detail.
First Cross Section. Ideally, the first cross section should be located sufficiently downstream from the bridge so that
the flow is not affected by the structure (ie, the flow has fully expanded). This distance should generally be
determined by field investigation during high flows and will vary depending on the degree of constriction, the shape
of the constriction, the magnitude of the flow, and the velocity of the flow. In order to provide better guidance to
determine the location of the fully expanded cross section, a study was performed by the Hydrologic Engineering
Center [HEC-1995]. This study focused on determining the expansion reach length, the contraction reach length,
and the expansion and contraction energy loss coefficients.
For this example, cross section number 5.29 was initially considered as the cross section of fully expanded flow. This
cross section was determined by field investigations as the approximate location of fully expanded flow during the
high flow event. After the pressure/weir flow analysis was performed, the location of this cross section was
evaluated using the procedures as outlined in the recent HEC study [HEC-1995]. The procedures required flow
parameters at the initially chosen location to evaluate the location of the cross section. These procedures will be
described after the pressure/weir flow analysis is performed near the end of this example.
Second Cross Section. The second cross section used by the program to determine the energy losses through the
bridge is located a short distance downstream of the structure. This section should be very close to the bridge, and
reflect the effective flow area on the downstream side of the bridge. For this example, a roadway embankment
sloped gradually from the roadway decking on both sides of the roadway. Cross section 5.39 was located at the toe
of the roadway embankment and was used to represent the effective flow area on the downstream side of the
bridge opening. The program will superimpose the bridge geometry onto this cross section to develop a cross
section inside the bridge at the downstream end.
Third Cross Section. The third cross section is located a short distance upstream from the bridge and should reflect
the length required for the abrupt acceleration and contraction of the flow that occurs in the immediate area of the
opening. As for the previous cross section, this cross section should also exhibit the effective flow areas on the
upstream side of the bridge. For this example, cross section 5.41 was located at the toe of the roadway
embankment on the upstream side of the bridge. Similar to the previous cross section, the program will
superimpose the bridge geometry onto this cross section to develop a cross section inside the bridge at the
upstream end.
Fourth Cross Section. The fourth cross section is located upstream from the bridge where the flow lines are parallel
and the cross section exhibits fully effective flow. For this example, cross section 5.44 was initially used as this
section where the flow lines were parallel. After the pressure/weir flow analysis, the location of this cross section
was evaluated using the procedures as outlined in the HEC study [HEC-1995]. This evaluation will be presented in
the discussion near the end of this example.
Manning's n Values. The Manning's n values were obtained from the field data displayed on the USGS atlas. For
some of the cross sections, the Manning's n values changed along the width of the overbank areas and the
horizontal variation in n values option was selected, such as for cross section 5.99. This option was performed from
the Cross Section Data Editor by selecting Options and Horizontal Variation in n Values. This caused a new
column to appear under the Cross Section X-Y Coordinates heading (as shown in the first figure). For cross section
5.99, the n values changed at the X-coordinates of 518 (in the left overbank), 866 (the main channel left bank
station), and 948 (the main channel right bank station). This can be seen by scrolling down in the coordinates
window of the cross section editor. The overbank areas have densely wooded areas, which created the necessity for
the variation in n values. The final data are shown in the cross section plot in the figure below.
Levees. As can be seen in the plot of cross section 5.99 in the figure above, there exists a large area to the left of the
main channel that is lower in elevation than the invert of the main channel. During the analysis, the program will
consider the water to be able to go anywhere in the cross section. The modeler must determine whether or not the
lower area to the left of the main channel can initially convey flow. If the area cannot convey flow until the main
channel fills up and then overtops, then the levee option should be used. For this example, a left levee was
established at the left main channel bank station for river station 5.99. This prevents water from being placed to the
left of the levee until the elevation of the levee is reached. The elevation selected for this levee was the elevation of
the left side of the main channel.
To insert the levee, Options and then Levees were selected from the Cross Section Data Editor. This resulted in
the display shown in the figure below. The values for this example were station 866 and elevation 214.8 for a left
levee. As can be seen in the figure above, the levee is displayed as a small square located at the LOB station.
Additionally, a note appears identifying the selection of a levee for the specific cross section at the bottom of the
Cross Section Data Editor. This note can be seen in the box at the bottom of the top figure. Levee options were
selected for other cross sections in addition to cross section 5.99. In each case, the modeler needs to view each
cross-section and determine whether the levee option needs to be utilized.
Contraction/Expansion Coefficients. The contraction and expansion coefficients are used by the program to
determine the transition energy losses between two adjacent cross sections. From the data provided by the recent
HEC study [HEC-1995], gradual transition contraction and expansion coefficients are 0.1 and 0.3, and typical bridge
contraction and expansion coefficients are 0.3 and 0.5, respectively. For situations near bridges where abrupt
changes are occurring, the coefficients may take larger values of 0.5 and 0.8 for contractions and expansions,
respectively. A listing of the selected values for this river reach can be viewed by selecting Tables and then
Contraction/Expansion Coefficients (Steady Flow) from the Geometric Data Editor. This table is shown in the
figure below and displays the values selected for the river cross sections. Typical gradual transition values were
selected for stations away from the bridge. However, near the bridge section, the coefficients were increased to 0.3
and 0.5 to represent greater energy losses. For additional discussion concerning contraction and expansion
coefficients at bridges, refer to Chapter 5 of the Hydraulic Reference Manual.
This completed the input for the cross section geometric data. Next, the bridge geometry data was entered as
outlined in the proceeding section.
Bridge Deck and Roadway Geometry. From the Bridge/Culvert Data Editor, the Deck/Roadway icon was
selected and this activated the Deck/Roadway Data Editor, as shown in the following figure.
The first input at the top of the editor is the distance from the upstream side of the bridge deck to the cross section
immediately upstream from the bridge (cross section 5.41). This distance was determined to be 30 feet from the
USGS atlas. In the next field, the bridge deck width of 40 feet was entered. Finally, a weir flow coefficient of 2.6 was
selected for the analysis. (Additional discussion of the weir flow coefficient will be presented in the calibration
section.)
The central section of the Deck/Roadway Editor is comprised of columns for input of the station, high cord
elevation, and low cord elevation for both the upstream and downstream sides of the bridge deck. The data are
entered from left to right in cross section stationing and the area between the high and low cord is the bridge
structure. The stationing of the upstream side of the deck was based on the stationing of the cross section located
immediately upstream. Likewise, the stationing of the downstream side of the deck was based on the stationing of
the cross section placed immediately downstream.
If both the upstream and downstream data are identical, the user needs only to input the upstream data and then
select Copy Up to Down to enter the downstream data.
As a final note, the low cord elevations that are concurrent with the ground elevation were entered as a value lower
than the ground elevation. The program will automatically clip off and remove the deck/roadway area below the
ground. For example, at station 0, a low cord elevation of 200 feet was entered. However, the actual ground
elevation at this point is approximately 216 feet. Therefore, the program will automatically remove the area of the
roadway below the ground. Additionally, the last station was entered as a value of 2000 feet. This stationing
ensured that the roadway and decking extended into the limits of the cross section geometry. As described
previously, the program will clip off the area beyond the limits of the cross section geometry.
The US and DS Embankment SS (upstream and downstream embankment side slope) values were entered as 2
(horizontal to 1 vertical). These values are used for graphical representation on the profile plot and for the WSPRO
low flow method. The user is referred to Example 13 - Bogue Chitto Single Bridge (WSPRO) for a discussion on the
use of this parameter with the WSPRO method. The WSPRO method is not employed for this example.
At the bottom of the Deck/Roadway Data Editor, there are three additional fields for data entry. The first is the
Max Allowable Submergence. This input is a ratio of downstream water depth to upstream energy, as measured
above the minimum weir elevation. When the ratio is exceeded, the program will no longer consider the bridge deck
to act as a weir and will switch the computation mode to the energy (standard step) method. For this example, the
default value of 0.95 (95 %) was selected, however this value may be changed by the user.
The second field at the bottom of the editor is the Min Weir Flow Elevation. This is the elevation that determines
when weir flow will start to occur over the bridge. If this field is left blank (as for this example), the program will
default to use the lowest high cord value on the upstream side of the bridge. Finally, the last field at the bottom of
the editor is the selection of the Weir Crest Shape. This selection will determine the reduction of the weir flow
coefficient due to submergence. For this example, a broad crested weir shape was selected. Upon entering all of the
above data, the OK button was selected to exit the Deck/Roadway Data Editor.
Bridge Pier Geometry. From the Bridge/Culvert Data Editor, select the Pier icon. This will result in the display
shown in the figure below. The modeler should not include the piers as part of the ground or bridge deck/roadway
because pier-loss equations use the separate bridge pier data during the computations.
The program will establish the first pier as pier number 1. As shown in the previous figure, the upstream and
downstream stations were entered for the centerline of the first pier. The upstream and downstream stations were
based on the geometry of the cross sections located immediately upstream (cross section 5.41) and immediately
downstream (cross section 5.39) of the bridge. The user needs to be cautious placing the pier centerline stations
because the X-coordinates for the upstream and downstream cross section stationing may be different. This is to
ensure that the piers "line up" to form the correct geometry. For this example, the pier centerline stations are 470,
490, 510, 530, 550, 570, 590, 610, and 630 for the nine piers. Each pier was set to start at an elevation of 200 feet (this
elevation is below the ground level and the excess will be removed by the program) and end at an elevation of 216
feet (this elevation is inside the bridge decking and the excess was removed by the program). Additionally, each pier
had a continuous width of 1.25 feet. After entering the data, the OK button was selected and the schematic of the
bridge with the piers was displayed on the Bridge/Culvert Data Editor as shown on the figure below. (Note: The
figure in the text displays the ineffective flow areas that will be added in the next section.)
The cross sections are developed by superimposing the bridge data on the cross sections immediately upstream
(5.41) and immediately downstream (5.39) of the bridge. The top cross section reflects the geometry immediately
inside the bridge on the upstream side and the bottom cross section reflects the geometry immediately inside the
bridge on the downstream side.
While viewing the bridge, the modeler can select to view just the upstream, just the downstream, or both of the
cross section views. This is performed by selecting View and then the required option. Additionally, from the View
menu, the user should select Highlight Weir, Opening Lid and Ground as well as Highlight Piers. These options
enable the modeler to view what the program will consider as the weir length, bridge opening, and pier locations.
Any errors in the data may appear as inconsistent images with these options. Also, the zoom-in option will allow the
user to examine data details.
As a final note for the bridge geometry, a bitmap image of the bridge was added to the geometry file (denoted by a
red square on the river system schematic in the first figure on this page). The user can view this image by selecting
the View Picture icon on the Geometric Data Editor.
For this example, ineffective flow areas were included on both the upstream cross section (5.41) and the
downstream cross section (5.39). To determine an initial estimate for the stationing of the ineffective flow areas, a
1:1 ratio of the distance from the bridge to the cross section was used. For this example, section 5.41 is located 30
feet upstream of the bridge. Therefore, the left and right ineffective flow areas were set to start at 30 feet to the left
and right of the bridge opening. Similarly, cross section 5.39 is located 30 feet downstream from the bridge and the
ineffective flow areas at this cross section were set at 30 feet to the left and right of the bridge opening.
To determine the initial elevation of the ineffective flow areas for the upstream cross section, a value slightly lower
than the lowest high cord elevation was used. This ineffective flow elevation was chosen so that when the water
surface becomes greater than this ineffective elevation, the flow would most likely be weir flow and would be
considered as effective flow. At the downstream cross section, the elevation of the ineffective flow area was set to
be slightly lower than the low cord elevation. This elevation was chosen so that when weir flow occurs over the
bridge, the water level downstream may be lower than the high cord, but yet it will contribute to the active flow
area. (Additional discussion of the selection of these elevations is described in the calibration section of this
example.)
To enter the ineffective flow areas, from the Geometric Data Editor select the Cross Section icon. Toggle to cross
section 5.41 and select Options and then Ineffective Flow Areas. This will result in the display shown in the figure
below. The default option (normal) is to enter the areas as a left station and elevation and/or a right station and
elevation. For this example, both the left and right ineffective flow areas were used.
The left and right ineffective flow stations were entered as 420 and 677 feet, respectively. These values are 30 feet to
the left and right of the bridge opening, as discussed previously. The elevation was then entered as 216.7 feet, a
value slightly lower than the high cord elevation. These entries imply that all the water to the left of the left station
and to the right of the right station will be considered as ineffective flow until the water level exceeds the elevation
of 216.7 feet.
Similarly, ineffective flow areas were set at river station 5.39 with a left station at 420 and a right station at 677, both
at an elevation of 215.0 feet. The OK button was selected and the ineffective flow areas appeared as green triangles,
as shown previously on the graph above. Additionally, the ineffective flow areas will appear on the plots of the cross
sections. Finally, a note will appear in the box at the bottom of the Cross Section Data Editor that states an
ineffective flow exists for each cross section for which this option was selected.
exceed the highest point of the low cord on the upstream side of the bridge). High flow occurs when the water
surface encounters the highest point of the low cord on the upstream side of the bridge. Finally, combination flow
occurs when both low flow or pressure flow occur simultaneously with flow over the bridge. The modeler needs to
select appropriate methods for both the low flow and for the high flow methods. For the combination flow, the
program will use the methods selected for both of the flows.
From the Geometric Data Editor, select the Bridge/Culvert icon and then the Bridge Modeling Approach button.
This will activate the Bridge Modeling Approach Editor as shown in the figure below. For this example, there is
only 1 bridge opening located at this river station and therefore the bridge number was 1. The following sections
describe the additional parameters of the bridge modeling editor. The modeler is referred to Chapter 6 of the User's
Manual and Chapter 5 of the Hydraulic Reference Manual for additional discussion on the bridge modeling
approach editor.
Low Flow Methods. Low flow exists when the flow through the bridge is open channel flow. The program has the
capability of analyzing low flow with four methods:
• Energy Equation (Standard Step)
• Momentum Balance
• Yarnell Equation (Class A only)
• WSPRO Method (Class A only)
These can be seen in the following figure:
The Energy Equation (Standard Step) method considers the bridge as just being part of the natural channel and
requires Manning's n values for the friction losses through the bridge and coefficients of contraction and expansion.
The Momentum Balance method performs a momentum balance through the bridge area and requires the selection
of a drag coefficient, Cd. This coefficient is used to estimate the force due to the water moving around the piers, the
separation of flow, and the resulting downstream wake. The Yarnell Equation is an empirical equation based on lab
experiments. Finally, the WSPRO method is an energy based method developed by the USGS for the Federal
Highway Administration.
At this time, the modeler needs to select which methods the program should compute and which method the
program should use. The modeler can select to have the program compute particular methods or all of the
methods. Then, the modeler needs to select which method the program will use as a final solution. Alternatively,
the modeler can select the computation of several or all of the methods, and then have the program use the
method with the greatest energy loss for the final solution. This will allow the modeler to view the results of all the
methods and compare the results of the different techniques.
For this example, the Energy, Momentum, and Yarnell methods were selected to be computed. For the momentum
method, a drag coefficient Cd = 2.00 was entered for the square nose piers and for the Yarnell method, a value of K =
1.25 was entered. Finally, the method that resulted in the greatest energy loss was selected to be used for the
solution. (The user is referred to Example 13 for an application of the WSPRO method.)
High Flow Methods. High flows occur when the water surface elevation upstream of the bridge is greater than the
highest point on the low cord of the upstream side of the bridge. Referring to the most recent figure, the two
alternatives for the program to compute the water surface elevations during the high flows are: Energy Only
(Standard Step) or Pressure and/or Weir Flow. The Energy Only (Standard Step) method regards the flow as open
channel flow and considers the bridge as an obstruction to the flow. Typically, most bridges during high flows may
act primarily as just an obstruction to the flow and the energy method may be most applicable.
As a second method for the analysis of high flows, the program can consider the flow to be causing Pressure Flow
and/or Weir Flow. For pressure flow, there are two possible scenarios. The first is when only the upstream side of the
bridge deck is in contact with the water. For this scenario, the submerged inlet coefficient, Cd, was set to be 0.34.
(This value was arrived at during the calibration, which is described later in this example.) The second scenario for
pressure flow is when the bridge constriction is flowing completely full. For this situation, the submerged inlet and
outlet coefficient was set to 0.80.
The program will begin to calculate either type of pressure flow when the computed low flow energy grade line is
greater than the highest point of the upstream low cord. Alternatively, the user can set the elevation at which
pressure flow will begin to be checked, instead of the highest low cord value. This value can be entered as the last
input to the Bridge Modeling Approach Editor (see figure above). For this example, this field was left blank which
implies that the program used the highest value of the low cord (the default). As an additional option, the user can
select to have the program begin to calculate the pressure flow by using the value of the water surface instead of
the value of the energy grade line. This is accomplished from within the Bridge/Culvert Data Editor by selecting
Options and then Pressure flow criteria. This will result in the display shown in the figure below. For this example,
the option to use the upstream energy grade line was chosen.
Finally, for the high flow analysis, Weir Flow occurs when the upstream energy grade line elevation (as a default
setting) exceeds the lowest point of the upstream high cord. The weir flow data was entered previously in the Deck/
Roadway Data Editor. At this point, all of the bridge data have been entered. The user should exit the geometry
data editors and save the geometry data. For this example, the geometry data was saved as the file "Beaver Cr. +
Bridge - P/W."
the reach so no flow change locations were used. Additionally, the three profile names were changed from the
default values of "PF#1," etc., to "25 yr," "100 yr," and May '74 flood," respectively. These names will be used to
represent the flow profiles when viewing the output.
To enter the boundary conditions, the Reach Boundary Conditions button was selected and this resulted in the
display shown in the following figure. For this example, a subcritical analysis was performed. Therefore, a
downstream boundary condition was required for each flow value. The mouse arrow was placed over the
downstream field and then the box was selected (highlighted). Then, 1 of the 4 boundary conditions was selected
and this caused the type of boundary condition that was chosen to appear in the downstream end of the reach.
For this example, Known W. S. was selected. This caused the input editor as shown below to appear. For each of
the flows, the known downstream water surface elevations of 209.5, 210.5, and 211.8 feet were entered for the flows
1, 2, and 3, respectively. These values were obtained from observed data on the USGS Atlas.
For the purposes of the analysis, if the downstream boundary conditions are not known, then the modeler should
use an estimated boundary condition. However, this may introduce errors in the region of this estimated value.
Therefore, the modeler needs to have an adequate number of cross sections downstream from the main area of
interest so that the boundary conditions do not effect the area of interest. Multiple runs should be performed to
observe the effect of changing the boundary conditions on the output of the main area of interest. For a detailed
explanation of the types of boundary conditions, refer to Chapter 7 of the User's Manual and Chapter 3 of the
Hydraulic Reference Manual. After entering the boundary condition data, the OK button was selected to exit the
editor. This completed the necessary input for the flow data and the steady flow data was then saved as "Beaver Cr.
- 3 Flows."
Select Options and then ensure that there is a check mark "" in front of Check data before execution. This will
cause the program to check all of the input data to ensure that all pertinent information was entered. Next, the
options were saved as a plan entitled "Press/Weir Method," with a Short ID entered as "Press/Weir." Finally,
COMPUTE was selected at the bottom of the window.
graph, displaying the water surface elevations and critical depth lines for all three profiles. (Note: the variables that
are displayed can be changed by selecting Options and then Variables.)
From the graph below, it can be seen that all three of the flow profiles are occurring in the subcritical flow regime.
This ensures that for the low flow analysis, Class A low flow (subcritical flow) was occurring through the bridge. Low
flow occurred for the first (5000 cfs) and for the second (10000 cfs) flow profiles. For the high flow, the method of
analysis was chosen to be pressure/weir flow. Pressure and weir flow occurred during the third flow profile (14000
cfs). One way to determine the type of flow that occurred is by viewing the bridge only output table. This table is
presented in the table below and was activated from the main program window by selecting View, Profile Table,
Std. Tables, and then Bridge Only.
For this example there is only one bridge located at river station 5.40, as listed in the table. Pressure flow
calculations were set to begin when the energy grade line elevation of the upstream section (5.41) was greater than
the highest elevation of the upstream low cord (215.7 ft). The first column in the following table shows the energy
grade line elevation of the upstream section (EG US) and the second column shows the elevation when pressure
flow was set to begin. A comparison of these two columns shows that pressure flow occurred for the third profile.
Additionally, it can be seen that weir flow occurred for the third profile, since there is a weir flow value for the third
profile. The following sections detail the output for the first two profiles and then for the third flow profile.
First and Second Flow Profiles. The first (5000 cfs) and second (10000 cfs) flow profiles were both computed using
the low flow methods of: Energy, Momentum, and Yarnell. From the main program window, select View, Profile
Table, Standard Tables, and then Bridge Comparison. This will provide a comparison table for the different
energy loss methods and is shown in the figure below.
In the figure below, the three rows display the results for each of the three flow profiles, in ascending order. The
river station is set at 5.4 (the only bridge location for this reach). The fourth column shows the water surface
elevation immediately upstream of the bridge. The sixth, seventh and eighth columns show the results of the low
flow methods that were chosen to be computed: Energy, Momentum, and Yarnell Methods, respectively. The
program compared the results and used the value with the greatest energy loss. For the first and second profiles,
the energy method calculated the greatest energy losses and the program used the results of 213.31 and 215.67
feet, respectively.
Third Flow Profile. The third flow profile was computed for a flow of 14000 cfs. As can be seen in the table above,
approximately 3060 cfs of the total flow was weir flow. The remaining flow, approximately 10940 cfs, was pressure
flow through the bridge opening. As can be seen in the table above, the energy grade line necessary for pressure-
only flow was 221.66 feet. Since this value is greater than the upstream high cord, weir flow also developed.
Therefore, the program used the pressure/weir energy value as the solution to the high flow method, namely 217.68
feet.
calculated using the same low flow procedures. For the third profile, the energy method was used to calculate the
energy losses through the bridge for this high flow and resulted in an energy gradeline elevation of 217.37 feet.
where:
= expansion reach length, ft
= expansion ratio
= main channel Froude number at the cross section immediately downstream of the bridge (cross
section 5.39 for this example)
= main channel Froude number at the cross section of the fully expanded flow (initially cross section
5.29 for this example)
= average length of obstruction caused by the two bridge approaches, ft
= total discharge,
(Note: The subscripts used in the first equation at the top of this page and all subsequent equations reflect the river
station numbering for this example.)
From the field data, the average length of the obstruction is approximately 740 feet and the total discharge, Q, is
14000 cfs for the high flow event. From the initial analysis, the values of the Froude number at cross-section 5.39
was 0.37 and at cross section 5.29 was 0.30. (Both of these values were the same for pressure/weir and energy
methods and are displayed on Standard Table 1.) Substituting the values into the equation above yielded that the
expansion reach length, Le, was approximately 778 feet. This equation has a standard error of 96 feet, which yields
an expansion reach length range from 682 to 874 feet to define the 68% confidence band. The distance used for the
expansion reach length (the distance from cross section 5.39 to cross-section 5.29) was set to be 500 feet in the
main channel, which is less than the recommended range from the equation. The modeler now has the option to
adjust this length so that it is within the calculated range. Then, after a new analysis, the new Froude numbers
should be used to calculate a new expansion reach length. If the geometry is not changing rapidly in this region,
then only 1 or 2 iterations should be necessary to obtain a constant expansion reach length value.
For this example, the location of the fully expanded cross section was changed to reflect the new expansion reach
length (778 feet) by interpolating a new cross section to be located 778 feet downstream from cross section 5.39. To
accomplish this, the following steps were taken. First, the pressure/weir flow plan ("Pressure/Weir Method") was
activated. Then, in the Geometric Data Editor, a cross section was interpolated between cross sections 5.29 and
5.21*. This interpolated cross section (5.24*) was set to be at a distance of 278 feet downstream from cross section
5.29. Since cross section 5.29 was already located 500 feet downstream from cross section 5.39, the location of
cross section 5.24* was then 500 + 278 = 778 feet downstream from cross section 5.39. The cross sections 5.29, 5.27*,
and 5.21* were then deleted.
Finally, this new geometric data was saved as a file called "Bvr.Cr.+Bridge - P/W: New Le, Lc." On the Steady Flow
Analysis Window, a Short ID was entered as "P/w+NewLeLc" and then the new geometry file and the original
steady flow data file were saved as a plan entitled "Press/Weir Method: New Le, Lc." The Geometric Data Editor
was then reactivated and the Bridge Modeling Approach Data Editor was selected. The high flow method was
chosen to be the energy method and then the geometry data was saved as "Bvr.Cr.+Bridge - Energy: New Le, Lc."
The Steady Flow Analysis Window was activated, a Short ID was entered as "En.+New LeLc," and the new energy
geometry file and the steady flow data file were then saved as a new plan entitled "Energy Method : New Le, Lc."
This procedure created two new plans, with each plan containing the necessary interpolated geometry and
appropriate high flow calculation methods.
Each of the two new plans were then executed and the resulting flow parameter values were reentered into the
above equation. (Note: The Froude number for river station 5.29 was replaced by the Froude number at river station
5.24*.) The final mean value of the expansion reach length was then determined to be 750 feet, with a range of ±96
feet to define the 68% confidence band. The adjusted reach length value of 778 feet is within the confidence band
and no additional iterations were computed. Finally, the expansion ratio (ER) as described in the above equation
should not exceed 4:1 and should not be less than 0.5:1. For this example, the final expansion ratio was ER = (778) /
(740) = 1.05, which is within the acceptable range. In the above procedure, the modeler is directed to Chapter 4 of
the Hydraulic Reference Manual for additional information on cross section interpolation and to Chapter 5 of the
User's Manual for further discussion on working with projects.
Upon reviewing the above procedures, the modeler can open either of these new plans and the corresponding
geometry and flow data files will be activated. For this example, the results of the water surface profiles for the new
plans are approximately equal to the results obtained from the original geometry for both of the pressure/weir flow
and energy methods, respectively. However, the modeler should apply the above procedures to evaluate the
location of the expansion reach length for each specific application.
Finally, during the procedures as outlined above, if the location of the expansion reach length produces a distance
sufficiently far downstream from the bridge, then the modeler may be required to include additional cross sections
within this reach length to accurately predict the energy losses. This may be accomplished by inserting cross
sections and providing the appropriate ineffective flow areas at each cross section according to their location with
respect to the bridge opening.
The contraction reach length, Lc, is defined as the distance from the cross section located immediately upstream of
the bridge (5.41) to the cross section that is located where the flow lines are parallel and the cross section exhibits
fully effective flow (5.44). To evaluate this reach length, the regression results (shown in the equation below) from
the US Army Corps of Engineers (USACE) study \[HEC-1995\] was used.
The equation is:
where:
= contraction reach length, ft
= contraction ratio
= main channel Froude number at the cross section immediately downstream of the bridge (cross
section 5.39 for this example)
= main channel Froude number at the cross section of fully expanded flow (initially cross section 5.29
for this example)
= discharge conveyed in the two overbanks at cross section 5.44, cfs
= total discharge,
= Manning n value for the overbanks at section 5.44
= Manning n value for the main channel at section 5.44
From the field data and the results of the initial analysis, the Froude numbers at sections 5.39 and 5.29 were 0.37
and 0.30, respectively, the total over bank flow at cross section 5.44 was approximately 9780 cfs (an average of 9880
for the pressure/weir flow and 9685 for the energy method), the total flow was 14000 cfs, the weighted n value for
both of the overbanks was 0.069, the n value for the main channel was 0.04, and the average length of the
obstruction was 740 feet. Substitution of these values into this equation yielded the contracted reach length of 478
feet. This equation has a standard error of 31 feet which results in a contraction reach length range from 447 to 509
feet to define the 68% confidence band. For this example, the distance from cross section 5.44 to cross section 5.41
was set at 170 feet along the main channel. Since this value was outside of the confidence range, the reach length
was adjusted to reflect the new contraction reach length.
The adjustment of the geometry for the new contraction reach length was performed similarly to the adjustment
procedure for the expansion reach length. A cross section (5.49*) was interpolated between river stations 5.525*
and 5.44 that was set to be 308 feet upstream of river station 5.44. Then, the reach length from river station 5.49* to
river station 5.41 was 308 + 170 = 478 feet, the required contraction reach length. This new river station can be
viewed by opening the plan "Press/Weir Method: New Le, Lc" for the pressure/weir method or the plan "Energy
Method: New Le, Lc" for the energy method analysis.
Finally, after the subsequent analysis with the new contraction reach length, the new flow parameters were entered
into the second equation and yielded a contraction reach length of 499 feet, with a range from 468 to 530 feet to
define the 68% confidence band. The adjusted contraction reach length of 478 feet is within this range and no
additional iterations were necessary. (Note: The Froude number for river station 5.29 was replaced by the Froude
number at river station 5.24*.)
As a final criteria, the contraction ratio (CR) should not exceed 2.5:1 nor should it be less than 0.3:1. For this
example, the final contraction ratio was CR = (499) / (740) = 0.67, which is within the acceptable range. The following
table shows a relationship of the values computed for the expansion reach lengths and the contraction reach
lengths during the iterations as described above. Additionally, the table shows the values as determined by the
USGS and the traditional USACE methods.
As can be seen from the data in the table above, the USGS method will typically provide a minimum criteria and the
traditional USACE method will provide a maximum length. The recent study [HEC-1995] was performed to provide
better guidance for the evaluation of the reach lengths and these values fall within the range as determined by the
previous two methods.
where:
=expansion coefficient
= hydraulic depth (flow area/top width) for the overbank cross section at 5.29)
= hydraulic depth in the main channel at cross section 5.29
From the data for the analysis, the hydraulic depth for the overbank at cross section 5.29 was 4.26 feet (an average
of 5.31 and 3.20 feet for the LOB and ROB, respectively, with the values being consistent for both the pressure/weir
flow and energy method) and the hydraulic depth of the main channel was 7.20 feet. Substitution of the values for
the variables yielded an expansion coefficient of 0.34. This is the median value and the range of ± 0.2 defines the
95% confidence band for the third equation. Therefore, the modeler should use the value of 0.34 as an initial value
and vary the coefficient by ±0.2. For this example, a value of 0.5 was initially used for the expansion coefficient in
the vicinity of the bridge. During the iterations for the contraction and expansion reach lengths, this coefficient was
reevaluated for each iteration. For the final value, the hydraulic depth in the overbank at cross section 5.24* was
4.11 feet (an average of 5.06 and 3.16 feet for the LOB and ROB, respectively, with the values being consistent for
both the pressure/weir flow and energy method) and the hydraulic depth of the main channel was 8.40 feet. This
yielded an expansion coefficient of 0.27, ±0.2. For the example, a final value of 0.5 was calibrated in order to match
the observed data.
In summary, the above recommendations for the expansion reach length, the contraction reach length, and the
expansion and contraction coefficients represent an improvement in the general methodology behind the
prediction of these values. The modeler is recommended to apply these new criteria as a more substantial method
for estimating the transition reach lengths. As a final note, after the initial analysis, the expansion and the
contraction reach lengths as well as the expansion and contraction coefficients should be evaluated
simultaneously. Then, adjustments should be made to the reach length and coefficient values before a subsequent
analysis is performed. Finally, the new data should be used to reevaluate all of the reach lengths and coefficients.
This procedure will ensure that the modeler is always using the current flow data for the analysis.
Finally, for the calibration, the ineffective flow areas at cross section 5.41 and 5.39 were set to balance the flow
going over the weir. Initially, the flow in the LOB and ROB at cross section 5.41 was drastically larger than the flows
in the LOB and ROB at the downstream section 5.39. This was due to the fact that the ineffective flow area
elevations at cross section 5.39 were previously set to a value in between the low cord and the high cord. When the
flow came over the weir, the depth downstream was less than the ineffective flow elevation, so the program initially
considered the overbank area as ineffective flow. This is not a realistic answer. When the flow goes over the weir it
contributes to the overbank flow at the downstream cross section and then this downstream area should be
considered as effective flow. Therefore, the ineffective flow elevations were lowered to allow the weir flow that
entered cross section 5.39 to become effective.
These calibrations were performed because the actual flow depths for the event were known. If these actual data
were not known, then the adjustment of the Manning's n values, pressure flow coefficients, and weir flow rate could
not be conducted to the refinement in the previous discussion. However, the balancing of the weir flow to the flow
in the LOB and ROB at river stations 5.39 and 5.41 could be performed. The calibration to the model accounted for a
more accurate determination of the water surface profile to the observed data.
plan and the "Pressure/Weir Flow: New Le Lc" plan. As shown in the figures, both methods produced the same
water surface profiles for the first two flows (both of the low flows). This is as would be expected because both
methods analyzed the two low flows using the same criteria. For the high flow, the water surface profile for the
energy method and the pressure/weir flow method varied slightly upstream of the bridge. The Zoom feature under
the Options menu can be used to obtain a closer view of the profiles.
The table at the bottom of this page shows a tabular comparison of the calculated water surface elevations for both
the pressure/weir flow method and the energy method to the observed data. The observed values with a ** denote
that the value may be in question due to there being only a few observed values in the vicinity of the location on the
atlas or because those that were provided were not in the active flow area.
In comparison of the calculated values to the observed values, both of the modeling approaches were able to
predict the actual water surface elevations within a reasonable tolerance. The largest errors occurred where the
observed water surface values are in question. Additionally, the observed values in the table are an average water
surface elevation over the area of effective flow where the cross sections are located. The one dimensional model
can only predict one resulting water surface; therefore, the fluctuations across the cross section will not occur in the
model as they did during the actual event.
In comparison of the pressure/weir method to the energy method, the greatest difference occurs at the bridge
structure. The water surface elevations for the pressure/weir flow method inside the bridge are estimated using the
upstream and downstream flow depths.
Pressure/Weir Energy
5.6 Summary
This example demonstrated the use of HEC-RAS to analyze a river reach that contains a single bridge crossing. The
geometric data consisting of the cross sections and bridge geometry were entered for the reach along Beaver Creek,
as obtained from the USGS Atlas No. HA-601. Three flow values were used for the analysis, with the largest flow
coinciding with the flood event in May 1974. The first plan consisted of the geometry data (with the high flow
method selected as press/weir) and the flow data. A second plan was created with the selection of the energy
method for the high flow analysis. Review of results for these plans reflected the necessity for adjustments to the
expansion and contraction reach lengths.
After the adjustments were made, two new plans were created, one for the pressure/weir and one for the energy
method for the high flow analysis. The results of these two new plans were then compared to the observed water
surface elevations. From the comparison, the pressure/weir method resulted with the closest values to the
observed water surface elevations.
Additionally, on the right side of the data editor, the left and right stations of the main channel were entered. For
this cross section, the left side of the main channel is defined to start at station 972 and end at the right station of
1027 feet. Finally, the contraction and expansion coefficients were entered as 0.3 and 0.5, respectively. These
coefficients are used by the program to determine the energy losses due to the flow contracting or expanding as it
travels from one cross section to the next. Typical values of these coefficients for gradual transitions are 0.1 and 0.3
for contractions and expansions, respectively. However, at locations where there are sudden changes in the cross
section geometry (i.e., flow into or out of a culvert or bridge opening), the coefficients may take larger values. The
selection of these coefficients is discussed in detail in Chapter 3 of the Hydraulic Reference Manual. For this cross
section (being the section immediately upstream of the culvert opening), the coefficients were initially selected as
0.3 and 0.5 for the contraction and expansion, respectively. After the flow analysis, ranges for these values were
determined by using the methods outlined in the HEC-1995 research document. These ranges were compared to
the selected values and will be discussed near the end of this example.
The cross-section information for the other river stations were entered in a similar fashion as for river station
20.238. Finally, the ineffective flow areas of the cross sections were entered. This option allows the user to define
areas of the cross section that will contain water but the water is not flowing in the downstream direction. This
option is typically used at cross sections in the vicinity of a culvert or bridge. For this example, the ineffective flow
option was used at river station 20.238 (located immediately upstream of the culvert) and at river station 20.227
(located immediately downstream of the culvert).
From the Cross Section Data Editor, select River Station 20.238, Options, and then Ineffective Flow Areas. This
will display the Ineffective Flow Editor shown in the top figure below. River station 20.238 was surveyed at a
location 5 feet upstream from the culvert. Typically, the stationing of the ineffective flow areas are set on a 1:1 ratio
to the distance from the opening. When the culvert data are entered, the centerline stations of the two culverts will
be 996 and 1004 feet and each culvert will be 6 feet in diameter. Therefore, the left edge of the opening is at station
993 and the right edge is at station 1007. Using the 1:1 ratio, the left ineffective flow station was set to be equal to 5
feet left of the left opening. Similarly, the right ineffective flow station was set to be equal to 5 feet right of the right
side of the openings. These values were entered as 988 and 1012 feet for the left and right stations, respectively.
Finally, the elevation of the ineffective flow area was set to be equal to 33.7 feet, a value slightly lower than the high
cord on the upstream side of the roadway. Similarly, ineffective flow areas were set at cross section 20.227 at
stations of 991 and 1009 (since this cross section was only 2 feet downstream of the culvert outlet) and at an
elevation of 33.3 feet. The location of these ineffective flow areas will be discussed further during the analysis of the
output. Typically, the culvert information may be entered first and then the modeler can enter the location of the
ineffective flow areas more readily with the location of the culverts known.
The placement of the cross sections relative to the location of the culvert is crucial for accurate prediction of
expansion and contraction losses. The culvert routine (as does the bridge routine) utilizes four cross sections
specifically located on both sides of the structure to determine the energy losses through the culvert. (Additionally
the program will interpret two cross sections inside of the culvert by superimposing the culvert and roadway data
onto both the immediate downstream and immediate upstream cross sections from the culvert.) The following is a
brief summary for determining the locations of the four cross sections. This procedure is identical to the procedure
used for determining the cross section locations for a bridge analysis. The modeler should review the discussion in
Chapter 6 of the User's Manual and Chapter 6 of the Hydraulic Reference Manual for further discussion.
First Cross Section. Ideally, the first cross section should be located sufficiently downstream from the culvert so
that the flow is not affected by the structure (i.e., the flow has fully expanded). This distance, referred to as the
expansion length (Le), should be determined by: field investigation during high flows; the procedure outlined in a
recent study by the USACE \[HEC-1995\]; or other acceptable procedure. For this example, the criteria developed by
USACE \[HEC-1995\] research document was utilized to determine the expansion reach length. To utilize this
method, an initial length was estimated from values obtained in tables that are presented in the document and
provided in Appendix B of the Hydraulic Reference Manual. Then, after the flow analysis was completed, the
location was evaluated based on equations developed from the research. (This evaluation will be discussed near
the end of this example.)
First, the following criteria were required to determine the location of the first cross section:
where :
= Manning's n value for the overbank at cross section 20.251
= Manning's n value for the main channel at cross section 20.251
= culvert opening width, ft (m)
= total floodplain width, ft (m)
= slope, ft/mile
= average length of the side obstruction, ft
Substitution of the field data yields the results as shown above. With these values, the expansion ratio (ER) was
determined to range from 0.8 to 2.0 from Table B-1 in Appendix B of the Hydraulic Reference Manual. The
expansion ratio (ER), is the length of expansion (Le) divided by the average length of obstruction (Lobs). For this
example, an average value of 1.4 was initially used for the expansion ratio. Therefore, the expansion reach length
will be the expansion ratio times the average length of obstruction:
From the initial values of the cross section locations, the expansion reach length is the distance from cross section
20.227 to 20.189. This distance is initially set at 200 feet. From the above analysis, it was determined that the
distance should be approximately 100 feet. Therefore, an additional cross section was placed 100 feet downstream
from cross section 20.227.
To produce this cross section, field data should be utilized. If this data is not available, then the program has the
ability to interpolate a cross section. From the Geometric Data Editor, Tools, XS Interpolation, and then Between
2 XS's were selected. "Culvrt Reach" was selected as the river reach (the only reach in this example) and river
station 20.227 was entered as the upper river station (this will default to river station 20.189 as being the lower river
station). The maximum distance between the interpolated cross sections was set to be 100 feet and then the
interpolation was performed. This resulted in the display shown in the figure below. For additional information on
cross section interpolation, refer to Chapter 4 of the Hydraulic Reference Manual and Chapter 6 of the User's
Manual. The interpolation window was closed and the river schematic displayed the new interpolated cross section
at river station 20.208* (as shown in the first figure on this page). The number 20.208 was the default setting since
the distance chosen (100 ft) was equal to one half the previous reach length (200 ft).
After interpolating the new river station, from the Geometric Data Editor Tables and then Reach Lengths were
selected. This resulted in the table shown in the preceding figure except that the interpolated cross-section 20.208*
appeared and the distances from cross section 20.227 to 20.208* were 100 feet and from 20.208* to 20.189 were also
100 feet for the LOB, channel, and ROB. The program now considered cross section 20.208* as being the location of
the fully expanded cross section.
Second Cross Section. The second cross section used by the program to analyze the energy losses through the
culvert is located within a few feet downstream of the structure. This section should be close to the culvert (within a
few feet) and reflect the effective flow area on the downstream side of the culvert. Therefore, any ineffective flow
areas outside of the flow expanding out of the culvert, should not be used for conveyance calculations. For this
example, cross section 20.227 was located two feet downstream from the culvert. The ineffective flow areas were
developed previously using this distance. Finally, after the culvert and roadway geometry were entered, the
program superimposed the geometry onto this cross section to develop a cross section inside the culvert at the
downstream end.
Third Cross Section. The third cross section is located within a few feet upstream from the culvert and should
reflect the length required for the abrupt acceleration and contraction of the flow that occurs in the immediate area
of the opening. Similar to the second cross section, this cross section should also block the ineffective flow areas on
the upstream side of the culvert. For this example, cross section 20.238 was located five feet upstream of the
culvert. Similar to the second cross section, the program will superimpose the culvert geometry onto the third cross
section to develop a cross section inside the culvert at the upstream end.
Fourth Cross Section. The fourth cross section is located upstream from the culvert where the flow lines are
parallel and the cross section exhibits fully effective flow. The distance between the third and fourth cross section,
referred to as the contraction reach length, can be determined by: 1) field investigation during high flows, 2) the
procedure outlined in a recent study by the USACE \[HEC-1995\]; or 3) other acceptable procedure. For this
example, the criteria developed by USACE \[HEC-1995\] research document was utilized to determine the
contraction reach length. To utilize this method, an initial length was estimated from values obtained in Table B.2
in Appendix B of the Hydraulic Reference Manual. To use this method, the following criteria were necessary:
where the variables are as described previously. Substitution of the values yields the results as shown above. With
these values, the contraction ratio (CR) was determined to range from 0.8 to 1.5 from Table B-2. The contraction
ratio is the length of contraction (Lc) divided by the average length of obstruction. For this example, an average
value of 1.15 was used for the contraction ratio. Therefore, the contraction reach length will be the contraction ratio
times the average length of obstruction:
From the initial values of the cross section locations, the contraction reach length is the distance from cross section
20.251 to 20.238. This distance was initially set at 70 feet. From the above analysis, it was determined that the
distance should be approximately 80 feet. Because these values are so close together, the initial value of 70 feet will
be maintained. Finally, after the flow analysis was performed, the location of this cross section was evaluated and
will be discussed near the end of this example.
The central portion of the editor consists of fields to enter the stations and elevations of the deck/roadway. The
values for this example are as shown in the figure. If the upstream and downstream decking is identical, then the
user needs to only enter the upstream information and then select Copy Up to Down. (Note: For culverts, only the
high cord information is required. The program will automatically block out the area between the high cord and the
ground. For a bridge analysis, the low cord information is required to define the bridge opening. For a culvert
analysis, the culvert data will define the openings below the high cord for the locations of the culverts.)
The next two fields are the US and DS Embankment Side Slopes. These values were entered as 2 (horizontal to 1
vertical). For a culvert analysis, these values are only used for the profile plot.
The bottom of the editor consists of three additional fields. The first field is the Maximum Allowable submergence
ratio. This is the ratio of downstream flow depth to upstream energy, as measured from the minimum high cord of
the deck. When this ratio is exceeded for a bridge analysis, the program will switch from the weir flow equation to
the energy method to determine the upstream flow depth. For a culvert analysis, this ratio is not used because the
program cannot perform a backwater analysis through a culvert flowing full. Therefore, the weir analysis method
will always be used when overflow occurs.
The second field is the Minimum Weir Flow Elevation. This is the elevation that the program uses to determine
when weir flow will begin. If this field is left blank, the program will use the lowest value of the high cord on the
upstream side of the deck. Alternatively, the user can enter a value for the program to start checking for the
possibility of weir flow. For this example, an elevation of 33.7 feet was used. This is the elevation of the roadway
above the culvert openings on the upstream side of the culvert. (Note: This is also the minimum elevation of the
high cord and therefore, this field could have been left blank.)
Finally, the last field requires the selection of the weir crest shape: broad crested or ogee shaped. This selection is
used for the type for submergence correction. For this example, a broad crested weir shape submergence
correction was used. With all of the data entered, the OK button was selected to exit the Deck/Roadway Data
Editor.
Culvert Geometric Data. To enter the culvert geometric data, from the Bridge Culvert Data Editor, the Culvert
icon was selected. This activated the Culvert Data Editor as shown in the figure below. Each of the fields for the
editor are described in the following sections.
Culvert ID - By default, the identifier for the first culvert will be set to "Culvert #1." A culvert type is defined by the
shape, diameter (or rise and span), chart number, scale, length, n value, loss coefficients, upstream invert, and
downstream invert. If all of these parameters are the same for each culvert, then the modeler will only have one
culvert type. Then the modeler can enter up to 25 identical barrels for this culvert type, with each barrel occurring
at a different location (defined by the upstream centerline and downstream centerline). If any of the culvert
parameters change, then the modeler must define each culvert that is different as a separate type (to a maximum of
10 culvert types at the same river station), with each type containing up to 25 identical barrels. For this example, the
culvert consisted of only 1 culvert type (since all of the parameters were the same for each barrel) but it contained
two identical barrels (with each placed at separate upstream and downstream centerline locations).
Solution Criteria - The user has the option to select to use the result for inlet control or outlet control as the final
answer for the upstream energy grade line value. The default method is to use the highest of the two values, as was
selected for this example.
Shape - The culvert shape is chosen from the eight available shapes: circular, box, elliptical, arch, pipe arch, semi-
circle, low arch, or high arch. For this example, the culvert barrels were circular shape. To select the shape, press the
down arrow on the side of the shape field and highlight the desired shape.
Diameter, Rise, or Rise and Span - Depending on the shape of the culvert, the modeler must enter the inside
dimensions of the culvert. For a circular pipe, only the diameter is required. For other shapes, the rise is defined as
the inside vertical measurement and the span is the inside horizontal measurement. (Note: For box culverts with
chamfered corners, refer to the discussion in Chapter 6 of the Hydraulic Reference Manual.) For this example, the
circular culvert was set to have a diameter of 6 feet.
Chart # - Each culvert type and shape is defined by a Federal Highway Administration Chart Number. Depress the
down arrow next to this field to select the appropriate chart number. Once a culvert shape has been selected, only
the corresponding chart numbers available for that culvert shape will appear in the selections. Descriptions for the
chart numbers appear in Chapter 6 in the Hydraulic Reference Manual. For this example, the culvert chart was
selected as "1 - Concrete Pipe Culvert."
Scale - This field is used to select the Federal Highway Administration scale that corresponds to the selected chart
number and culvert inlet shape. Only the scale numbers which are available for the selected chart number will
appear for selection. Descriptions for the scale numbers appear in Chapter 6 in the Hydraulic Reference Manual.
The scale for this example was "1 - Square edge entrance with head wall."
Distance to Upstream XS - This is the distance from the inlet of the culvert to the upstream cross section (20.238). For
this example, this was a distance of 5 feet. On the Deck/Roadway Data Editor, a measure of 10 feet was entered for
the distance from the upstream side of the deck/roadway to the upstream cross section. Therefore, the culvert
entrance is located midway between the upstream side of the roadway and cross section 20.238.
Culvert Length - This field is the measure of the culvert (in feet or meters) along the centerline of the barrel. The
length of the culvert for this example was 50 feet. The program will add this 50 feet to the 5 foot distance from cross
section 20.238 (to the culvert entrance) and obtain 55 feet. The reach length from cross section 20.238 to 20.227 is
57 feet, which leaves 2 feet from the exit of the culvert to the downstream cross section.
Entrance Loss Coefficient - The value of the entrance loss coefficient will be multiplied by the velocity head at the
inside upstream end of the culvert to obtain the energy loss as the flow enters the culvert. Typical values for the
entrance loss coefficient can be obtained from Tables 6.3 and 6.4 in the Hydraulic Reference Manual. The entrance
loss coefficient for the concrete pipe in this example was set at 0.5.
Exit Loss Coefficient - To determine the amount of energy lost by the water as it exits the culvert, the exit loss
coefficient will be multiplied by the difference of the velocity heads from just inside the culvert at the downstream
end to the cross section located immediately downstream of the culvert exit. In general, for a sudden expansion,
the exit loss coefficient should be set equal to 1. However, this value may range from 0.3 to 1.0. For this example,
the exit loss coefficient was set to be equal to 1.0.
Manning's n for Top - This field is used to enter the Manning's n value of the top and sides of the culvert lining and is
used to determine the friction losses through the culvert barrel. Suggested n values for culvert linings are available
in many textbooks and also may be obtained from Table 6.1 in the Hydraulic Reference Manual. Roughness
coefficients should be adjusted according to individual judgment of the culvert condition. For this example, a
Manning's n value of 0.013 was used for the concrete culvert.
Manning's n Value for Bottom – This field is used to enter the Manning's n value of the bottom of the culvert. For
most culverts, this field will be the same as the Manning's n value for the top. However, if the culvert has a natural
bottom, or something has been placed in the bottom for fish passage, the n value may vary.
Depth to use Bottom n – This field is used to enter a depth inside of the culvert that the bottom n value is applied to.
If the bottom and top n value are the same, a value of zero should be entered.
Depth Blocked – This field is used to fill in a portion of the culvert. The user enters a depth, and everything below
that depth is blocked out.
Upstream and Downstream Invert Elevation - These two fields are used to enter the elevations of the inverts. For a
particular culvert type, all of the identical barrels will have the same upstream invert elevation and downstream
invert elevation. For this example, the upstream invert was set at an elevation of 25.1 feet and the downstream
invert was 25.0 feet.
Barrel Centerline Stations - This table is used to enter the name of each barrel, and the stationing (X-coordinates) of
the centerline of the culvert barrels. The upstream centerline is based upon the X-coordinates of the upstream cross
section (20.238) and the downstream centerline is based upon the X-coordinates of the downstream cross section
(20.227). This example employs two culvert barrels, with the centerline of the barrels occurring at stations 996 and
1004 feet, as measured on both cross sections. For this example, the X-coordinate geometry of both cross section
20.238 and cross section 20.227 are referenced from the same left station starting point. Therefore, the upstream
and downstream centerline stations are the same value and this will align the culvert in the correct configuration as
being parallel to the channel. The modeler must be cautious to ensure that the centerline stationing of the culvert
ends align the culvert in the correct position.
# Barrels - This field will automatically display the number of barrels entered by the user (determined by the
number of centerline stations entered). Up to 25 identical barrels can be entered for each culvert type, and this
example consisted of 2 identical barrels.
This completed the necessary geometric data for the analysis. The OK button at the bottom of the Culvert Data
Editor was selected and this displayed the culvert as shown in the following figure. The geometry data editors were
then closed and the geometry was saved as "Multiple Pipe Geometry."
As shown in the bottom figure, a known water surface elevation was then entered for each of the flow profiles that
were computed. For this example, the known water surface elevations of 29.8, 31.2, and 31.9 were entered for the
flows of 250, 400, and 600 cfs, respectively. Once the data were entered, the OK button was selected to exit this
window. This completed the necessary steady flow data for the analysis and the data were saved as "Multiple Pipe
Flow Data."
Expansion Reach Length. Cross section 20.208* was the interpolated section that was assumed to be at the
location where the flow became fully expanded. To evaluate the location of this cross section, the relationship
shown as Equation 3-1 was used. Equation 3-1 is applicable when the width of the floodplain and the discharge is
less than those of the regression data. The equation is:
where
= expansion ratio
= expansion reach length, ft
= average length of side obstruction, ft
= main channel Froude number at the cross section immediately downstream of the culvert (cross
section 20.227 for this example)
= main channel Froude number at the cross section of fully expanded flow (cross section
for this example)
= total discharge, ft3/s
(Note: The subscripts used in Equation 3-1 and all subsequent equations reflect the river station numbering for this
example.) From the analysis, the Froude numbers at cross sections 20.227 and 20.208*, for the flow of 600 cfs, are
0.32 and 0.14, respectively. Substituting these values into Equation 3-1 yields an expansion ratio of 1.51. This value
falls within the range of 0.8 – 2.0 as determined previously from the table values. With this new expansion ratio, the
expansion reach length is:
Additionally, the expansion ratio has a standard error of 0.26. Using the range of the ER from 1.25 (= 1.51 - 0.26) to
1.77 (= 1.51 + 0.26) yields an Le range from 88 to 124 feet to define the 68% confidence band for Equation 3-1. The
actual distance from cross section 20.227 to cross-section 20.208* was set at 100 feet. Therefore, the existing
expansion reach length seems appropriate. If the existing length had been significantly outside of the range of the
calculated expansion reach length, then a second iteration for the placement of the fully expanded cross section
and an additional analysis may be warranted. As a final check, the expansion ratio should not exceed 4:1 and should
not be less than 0.5:1.
Contraction Reach Length. Cross section 20.251 is located where the flow lines are parallel to the main channel. To
evaluate this location, Equation 3-2 from the research document \[HEC-1995\] was utilized. This equation is used
when the floodplain scale and dischargers are significantly different than those used in the regression analysis and
is:
where
= contraction ratio
= contraction reach length
= discharge conveyed by the two overbanks at cross section 20.251, cfs
= Manning's n value for the over banks at cross section 20.251
= Manning's n value for the main channel at cross section 20.251
From the analysis of the 600 cfs profile, the flow in the two overbanks at cross section 20.251 is 69.33 cfs and the
Manning's n values for the overbanks and main channel at cross section 20.251 are 0.10 and 0.04, respectively.
Substituting these values into Equation 3-2 yields a contraction ratio of 0.39. The standard error for this equation is
0.19 which yields a range of the CR from 0.20 to 0.58. From the results cited in the research document, a minimum
contraction ratio is 0.3:1 and a maximum ratio is 2.5:1. The calculated average value of 0.39 is very close to the
minimum value. If this value is used, the contraction reach length would be:
This is the median value for the range of 14 to 41 feet (using CR = 0.20 and 0.58, respectively). The actual distance
used for the contraction reach length (the length from river station 20.251 to river station 20.238) was 70 feet. For
this example, the contraction reach length was maintained at the 70 feet value. However, an additional analysis was
performed with the 27 feet contraction reach length and no appreciable difference in the water surface was
observed to reflect the necessity for the change of the contraction reach length to 27 feet.
Expansion Coefficient. The expansion coefficient can be obtained from Equation 3-3:
where
= coefficient of expansion
= hydraulic depth (flow area divided by the top width) for the overbank at cross section 20.208*
= hydraulic depth for the main channel at cross section 20.208*
From the analysis of the 600 cfs profile, the hydraulic depths for the overbanks and the main channel are 0.58 and
5.66 feet, respectively. Substitution of the values into Equation 3-3 yielded an expansion coefficient of 0.14. This is
the median value and the range of ± 0.2 defines the 95% confidence band for Equation 3-3. For this example, a value
of 0.5 was used. The value of the expansion coefficient is generally larger than the value used for the contraction
coefficient so the value of 0.5 will remain as the selected value. The regression equation (3-3) was developed for
bridges with overbank areas larger than the current example. Therefore, the data for this current example may not
be within the range of data used to develop the regression equation. The modeler can perform a sensitivity of this
coefficient by changing this coefficient and performing subsequent analyses. For this example, a value of 0.3 was
used during a subsequent analysis and no appreciable difference was observed in the resulting water surface.
Contraction Coefficient. From the research document \[HEC-1995\], the contraction coefficient is obtained by first
determining the relationship:
where
= culvert opening width, ft (m)
= total floodplain width, ft (m)
From Table B-3 of Appendix B in the Hydraulic Reference Manual, the recommended contraction coefficient range
is 0.3 - 0.5. The value selected for this example was the minimum value of 0.3, which reflects a value in between a
typical contraction and an abrupt contraction.
To investigate the first flow profile, from the main window select View, Cross Section Table, Type, and then
Culvert. Select the river reach of "Culvrt Reach," river station 20.237, profile 1, and culvert #1. This will display the
table shown in the figure below. The left column of the table shows a total flow rate of 250 cfs through the culvert.
Since the culvert has two identical barrels, this yields a flow of 125 cfs through each barrel. The table also shows
that the normal depth (3.56 ft) was greater than the critical depth (3.02 ft), which corresponds to subcritical flow
occurring through the culvert. At the bottom of the left column, the data shows that the culvert did not flow full for
any length of the culvert. Additional values such as the velocity in the culvert at the upstream and downstream ends
are displayed.
To determine the control of flow through the culvert (i.e., inlet or outlet), the values of the upstream energy grade
line necessary for inlet control (E.G. IC) and outlet control (E.G. OC) are shown in the left column of the above figure.
For the specified flow of 250 cfs, the upstream energy grade line for inlet control was 29.50 feet and the upstream
energy grade line for outlet control was 30.57 feet. The program will select the higher of these two values to
determine which type of control will occur (since "Highest Upstream EG" was selected on the Culvert Data Editor).
For this example, outlet control occurred and the energy gradeline used by the program is listed as the energy
grade line upstream of 30.57 feet. Finally, by using the values in the figure above and following the decision flow
chart shown as Figure 6.9 in the Hydraulic Reference Manual, the modeler can determine the procedure used by
the program to determine the water surface profile. For this flow of 250 cfs, the program used the Direct Step
Method to calculate the water surface profile.
For an analysis of the second flow (400 cfs), a similar procedure can be followed. For this example, the second flow
resulted in full flow along the entire length of the barrels of the culvert. The upstream water surface profile was
determined by using the FHWA full flow equations.
For the third flow (600 cfs), the inlet and outlet were submerged and weir flow occurred over the roadway. Select
the Culvert type Cross Section Table (as performed for the figure above) and select the third flow profile. This will
display the table shown in figure below. For this profile, the flow through the culvert was 525.75 cfs, 262.87 cfs
through each identical barrel. The energy grade line upstream was calculated to be 34.34 feet, which corresponds
to the outlet control energy grade line as shown.
The weir flow at river station 20.237 resulted with a value of 600 - 525.75 = 74.25 cfs. This occurred from an X-
coordinate of 945.10 to 1048.77, a distance of 103.67 feet. The main channel bank stations for cross sections 20.227
and 20.238 are at 972 and 1027. Therefore, the weir flow from 945.1 to 972 should balance with the flow in the LOB
at river stations 20.227 and 20.238. Additionally, the weir flow from 1027 to 1048.77 should balance with the flow in
the ROB at river stations 20.227 and 20.238. At river station 20.237, the amount of weir flow from 945.1 to 972 can
be approximated as: the specific length of weir divided by the total weir length times the weir flow. This was
calculated as:
Similarly, the amount of weir flow from 1027 to 1048.77 is approximately 15.4 cfs. Therefore, at river stations 20.227
and 20.238, the flow in the LOB should be approximately 19.26 cfs and the flow in the ROB should be approximately
15.4 cfs.
To determine the amount of flow in the overbanks at cross section 20.227, from the Cross Section Table window,
select Type and then Cross Section. Toggle to river station 20.227 for the third flow profile. The values for the flow
in the LOB and ROB are zero at this cross section. By toggling to river station 20.238, it was observed that the flows
in the LOB and ROB were 36.2 and 34.4 cfs, respectively. Therefore, these flow values in the LOB and ROB for both
river stations need to be adjusted to balance with the amount of weir flow.
The following discussion is provided as an example procedure that can be utilized to balance the weir flow with the
overbank flow. The modeler should compare the magnitude of the weir flow to the total flow rate to determine if
the procedure is practical for the specific situation.
To adjust the LOB and ROB flow values, first the situation at river station 20.227 was analyzed. Since there was not
any flow in the overbanks, the ineffective flow elevation was lowered from 33.3 to 32.0. This will allow for the flow
coming over the weir to become active at this downstream cross section. The elevation of 32.0 feet was chosen
because it is slightly lower than the calculated water surface (32.01 feet) at river station 20.227. (The user must be
cautious not to lower this elevation to a point where the ineffective flow will impact the second flow profile. Each
flow profile must be analyzed separately.)
As a second step to balance the weir and overbank flows, at river station 20.238 the Manning's n values in the
overbanks were raised from 0.1 to 0.4. Additionally, at river station 20.227, the Manning's n values for the overbanks
were decreased from 0.10 to 0.06. Since the n value is inversely proportional to the flow rate, the increase in n value
at river station 20.238 will cause a decrease in the flow rate in the overbank areas. Similarly, the decrease in the n
value at river station 20.227 will cause an increase in the flow rate in the overbank areas.
After these adjustments were made, the geometry file was saved as "Adjusted Ineffective + n Values." Then this
geometry file and the steady flow data file were saved as a plan entitled "Sp. Cr. Culverts - Adj. Weir Flow." The user
can activate this plan to review the remaining discussion of the output.
After the adjusted plan was executed, the weir flow at river station 20.237 was determined to be 74.20 cfs, as shown
in the Culvert Table of the figure below. This weir flow occurred from X-coordinates of 945.08 to 1048.78, a distance
of 103.65 feet. As calculated previously, the approximate portion of this weir flow that occurred from 945.11 to 972
is:
Similarly, the portion of the weir flow from 1027 to 1048.78 was approximately 15.6 cfs. These flow values were then
compared to the flow values in the LOB and ROB at river stations 20.227 and 20.238.
At river station 20.227, the flow in the LOB was 6.93 and the flow in the ROB was 5.60 cfs. These values are
approximately equal to the portions of the weir flow values as determined above. Therefore, these flow rates were
considered as being balanced with the weir flow. If the flows in the overbanks had been higher than the portions of
the weir flow, then the Manning's n values in the overbanks at river station 20.227 would have been increased until
a balance was achieved.
At river station 20.238, the flow in the left and right overbanks were 9.91 and 9.38 cfs, respectively. These flow rates
were considered to be reasonably in balance with the portions of the weir flow. If the flow rates were not in balance,
then the n values would have been adjusted further until a balance was achieved.
For both river stations 20.227 and 20.238, it should be noted that the flow rates do not exactly match the portions of
the weir flow as calculated above. An exact match is not warranted because the weir flow portions were
approximate and the weir flow that contributes to the left and right overbank is only a minor portion of the total
flow rate. If observed high water marks were available, the modeler could make adjustments to the data to more
accurately predict the actual water surface elevations.
For additional detailed analysis of the flow, the modeler should review the energy losses associated with the
contraction and expansion of the flow in the channel and the entrance and exit losses for the culvert to evaluate the
selected energy loss coefficients.
Finally, a three dimensional view of the water surface profiles is displayed in the figure below. This was activated
from the main program window by selecting View and then X-Y-Z Perspective Plots. The figure is available to aid the
user to view the calculated water surface profiles. The water surface image represents the hydraulic grade line at
the respective cross section locations.
6.4.4 Summary
For this example, a culvert was analyzed with three flows. The culvert was composed of two identical circular
barrels. During the review of the output, it was determined that initially the flow in the left and right overbanks at
the upstream and downstream cross sections from the culvert did not match the weir flow that was occurring. In
order to balance the weir flow with the overbank flows, the Manning's n values and the ineffective flow areas were
adjusted at the cross sections that bound the culvert.
The next example (Example 4) utilizes this data set and adds another culvert type to the geometric data. This
creates a multiple culvert analysis, each with multiple identical barrels.
The entering of data and the analysis of these culverts was identical to the procedures used for the single culvert
type as performed for Example 3. Additionally, the data used for Example 3 was modified for this example to include
a second culvert type. Therefore, Example 4 is presented as a continuation of Example 3 and the modeler should
review the data and the procedures as performed in Example 3 before reviewing Example 4.
The user is referred to Example 3 for the basic procedures of culvert analyses and to Chapter 6 of the Hydraulic
Reference Manual for a detailed discussion on modeling culverts. To review the procedures performed for this
example, from the main HEC-RAS window, select File and then Open Project. Go to the directory in which you have
installed the HEC-RAS example data sets. From the "Applications Guide\Example 4 – Multiple Culverts"
subdirectory, select the project labeled "Multiple Culverts - Example 4." This will open the project and activate the
following files:
Plan: "Spring Creek Multiple Culverts"
Geometry: "Multiple Culvert Geometry"
Flow: "Multiple Culvert Flow Data"
Deck/Roadway Data. The deck and roadway data editor is activated from the main program window by selecting
Edit, Geometric Data, the Bridge/Culvert icon, and then the Deck/Roadway icon. This will display the Deck/
Roadway Data Editor as shown in Example 3. The values for this current example were exactly the same as those
entered for Example 3.
Culvert Geometric Data. To enter the culvert geometric data, from the Bridge/Culvert Data Editor select the
Culvert icon. This will activate the Culvert Data Editor and display the data for culvert #1 as shown in the figure
below. Culvert #1 is defined exactly as described for Example 3: a circular culvert of diameter 6 feet with two
identical barrels and the other parameters as shown in the figure below. In the upper right corner of the editor, the
Rename button was selected and the new name "Circular" was entered. This will help the modeler during the
review of the output for this example.
The output analysis of Example 3 showed that during the flow of 600 cfs, the flow overtopped the roadway and
created weir flow. For this current example, we will consider the possibility that the flow must not overtop the
roadway and that the modeler desires to install additional culverts to alleviate this problem. Alternatively, to
reduce the upstream water surface depth, the modeler has extensive options available such as increasing the
diameter of the existing culverts or adding other identical barrels to the culvert. For this example, however, the
option of installing additional culverts was pursued.
From the Culvert Data Editor (as shown in the figure above), the Add button in the upper left corner was selected.
This cleared the entry fields and created a new input window as shown in the figure below. The identification name
"Culvert #2" appeared in the Culvert ID field in the upper right corner of the editor. The Rename button, located
immediately below the Culvert ID, was selected and a new name "Box" was entered. The additional fields of the
Culvert Data Editor will be described as follows.
Solution Criteria - The user has the option to select to use the result for inlet control or outlet control as the final
answer for the upstream energy grade line value. The default method is to use the highest of the two values, as was
selected for this example.
Shape - The culvert shape is chosen from the eight available shapes: circular, box, elliptical, arch, pipe arch, semi-
circle, low or high arch. For this second culvert, the culvert shape was selected as box. To select the shape, press the
down arrow on the side of the shape field and highlight the desired shape.
Diameter, Rise, or Rise and Span - Depending on the shape of the culvert, the modeler must enter the inside
dimensions of the culvert shape. For a box shape, the rise (vertical distance) and span (horizontal distance) must be
entered. A value of 5 feet for the span and 3 feet for the rise were entered. (Note: For box culverts with chamfered
Chart # - Each culvert type and shape is defined by a Federal Highway Administration Chart Number. Depress the
down arrow next to this field to select the appropriate chart number. Once a culvert shape has been selected, only
the corresponding chart numbers available for that culvert shape will appear in the selections. Descriptions for the
chart numbers appear in Chapter 6 on Table 6.5 of the Hydraulic Reference Manual. For this example, the culvert
chart was selected as "10 - 90 degree headwall; Chamfered or beveled inlet edges."
Scale - This field is used to select the Federal Highway Administration scale that corresponds to the selected chart
number and culvert shape. Only the scale numbers which are available for the selected chart number will appear
for selection. Descriptions for the scale numbers appear in Chapter 6 on Table 6.5 of the Hydraulic Reference
Manual. The scale for this example was "2 - Inlet edges beveled ½ in / ft at 45 degrees (1:1)."
Distance to Upstream XS - This is the distance from the inlet of the culvert to the upstream cross section (20.238). For
this example, this was a distance of 5 feet. On the Deck/Roadway Data Editor, a measure of 10 feet was entered for
the distance from the upstream side of the deck/roadway to the upstream cross section. Therefore, the culvert
entrance is located midway between the upstream side of the roadway and cross section 20.238.
Length - This field is the measure of the culvert (in feet or meters) along the centerline of the barrel. The length of
the culvert for this example is 50 feet.
The program will add this 50 feet to the 5 foot distance from cross section 20.238 (to the culvert entrance) and
obtain 55 feet. The reach length from cross section 20.238 to 20.227 is 57 feet, which leaves 2 feet from the exit of
the culvert to the downstream cross section.
Entrance Loss Coefficient - The value of the entrance loss coefficient will be multiplied by the velocity head at the
inside upstream end of the culvert to obtain the energy loss as the flow enters the culvert. Typical values for the
entrance loss coefficient can be obtained from Tables 6.3 and 6.4 in the Hydraulic Reference Manual. The entrance
loss coefficient for the concrete box culvert was set at 0.2.
Exit Loss Coefficient - To determine the amount of energy lost by the water as it exits the culvert, the exit loss
coefficient will be multiplied by the difference of the velocity heads from just inside the culvert at the downstream
end to the cross section located immediately downstream from the culvert exit. In general, for a sudden expansion,
the exit loss coefficient should be set equal to 1. However, this value may range from 0.3 to 1.0. For this example,
the exit loss coefficient was set to be equal to 1.0.
Manning's n for Top - This field is used to enter the Manning's n value of the top and sides of the culvert lining and is
used to determine the friction losses through the culvert barrel. Suggested n values for culvert linings are available
in many textbooks and also may be obtained from Table 6.1 in the Hydraulic Reference Manual. Roughness
coefficients should be adjusted according to individual judgment of the culvert condition. For this example, a
Manning's n value of 0.013 was used for the concrete culverts.
Manning's n Value for Bottom – This field is used to enter the Manning's n value of the bottom of the culvert. For
most culverts, this field will be the same as the Manning's n value for the top. However, if the culvert has a natural
bottom, or something has been placed in the bottom for fish passage, the n value may vary.
Depth to use Bottom n – This field is used to enter a depth inside of the culvert that the bottom n value is applied to.
If the bottom and top n value are the same, a value of zero should be entered.
Depth Blocked – This field is used to fill in a portion of the culvert. The user enters a depth, and everything below
that depth is blocked out.
Upstream and Downstream Invert Elevation - These two fields are used to enter the elevations of the inverts. For a
particular culvert type, all of the identical barrels will have the same upstream invert elevation and downstream
invert elevation. For this example, the upstream invert was set at an elevation of 28.1 feet and the downstream
invert was 28.0 feet.
Barrel Centerline Stations - This table is used to enter the stationing (X-coordinates) of the centerlines of the culvert
barrels. The upstream centerlines are based upon the X-coordinates of the upstream cross section (20.238) and the
downstream centerlines are based upon the X-coordinates of the downstream cross section (20.227). This example
employs two culvert barrels, with the centerlines of the barrels occurring at stations 988.5 and 1011.5 feet, as
measured on both cross sections. For this example, the X-coordinate geometry of both cross section 20.238 and
cross section 20.227 are referenced from the same left station starting point. Therefore, the upstream and
downstream centerline stations are the same value and this will align the culvert in the correct configuration as
being parallel to the channel. The modeler must be cautious to ensure that the centerline stationing of the culvert
ends align the culvert in the correct position.
# Barrels - This field will automatically display the number of barrels entered by the user (determined by the
number of centerline stations entered). Up to 25 identical barrels can be entered for each culvert type, and this box
culvert consisted of 2 identical barrels.
This completed the necessary geometric data for the culverts. The OK button was selected and the final
configuration of the culverts is as shown in the figure below. (Note: The ineffective flow areas will be adjusted
subsequently.) The modeler should now close the editor and save the data.
The stationing of the ineffective flow areas at cross section 20.238 can be determined by using a 1:1 relationship to
the distance of cross section 20.238 from the culvert entrance. For this example, cross section 20.238 is located 5
feet upstream from the culvert entrance. Therefore, the left ineffective flow station was set at 5 feet to the left of the
left edge of the culvert openings. The left edge of the culvert opening is at station 986. This places the left ineffective
flow station at station 986 - 5 = 981 feet. Similarly, the right ineffective flow station was set to be 5 feet to the right of
the right edge of the culvert opening. This equates to station 1014 + 5 = 1019. The elevation of the ineffective flow
areas will remain as originally set for Example 3, at an elevation of 33.7 feet.
The ineffective flow areas at cross section 20.227 are determined in a similar manner as used for cross section
20.238; however, cross section 20.227 is placed at a distance of 2 feet downstream of the culvert exit. Therefore, the
ineffective flow stationing at cross section 20.227 is 984 and 1016 for the left and right stations, respectively. The
elevations for the ineffective flow areas at cross section 20.227 were set at 33.6, a value slightly lower than the
lowest downstream high cord elevation.
This completed the required geometric data and the information was saved as the geometry file "Multiple Culvert
Geometry."
actual distances of the expansion reach length (from cross section 20.227 to cross section 20.189) and the
contraction reach length (from cross section 20.251 to cross section 20.238). It was then determined that the
contraction reach length was adequate for the analysis; however, the estimated expansion reach length was less
than the distance from cross section 20.227 to cross section 20.189. Therefore, an additional cross section, 20.208*,
was interpolated and included as the location where the flow would fully expand. After the analysis for Example 3,
the expansion and contraction reach lengths were computed using regression equations and compared to the
predetermined values. For this current example, the regression equations will again be used to evaluate the
expansion and contraction reach lengths.
Expansion Reach Length. Cross section 20.208* was the interpolated section that was assumed to be at the
location where the flow became fully expanded. To evaluate the location of this cross section, the relationship
shown below was used. It is applicable when the width of the floodplain and the discharge is less than those of the
regression data. The equation is:
where
= expansion ratio
= expansion reach length, ft
= average length of side obstruction, ft
= main channel Froude number at the cross section immediately downstream of the culvert (cross
section 20.227 for this example)
= main channel Froude number at the cross section of fully expanded flow (cross section 20.208* for
this example)
= total discharge, ft3/s
(Note: The subscripts used and all subsequent equations reflect the river station numbering for this example.) From
the analysis, the Froude numbers at cross sections 20.227 and 20.208* for the flow of 600 cfs are 0.18 and 0.14,
respectively. Substituting these values into this equation yields an expansion ratio of 1.06. The standard error for
this equation is 0.26, which yields a range of the ER from 0.80 to 1.32 to define the 68% confidence band. For this
current example, the average length of obstruction is determined to be approximately equal to 60 feet. With this
length of obstruction, the expansion reach length is:
This is the median value of the range from 50 to 80 feet (using ER = 0.80 and 1.32, respectively). The actual distance
from cross section 20.227 to cross section 20.208* was set at 100 feet, which is only slightly greater than the
maximum value of the range as determined from the equation above. Therefore, the existing reach length was not
adjusted. Also, since the regression equation was based on studies conducted for low flow through bridges and
since the flow for this example is less than the flow rates used to develop the equation, the expansion reach length
was not adjusted from the existing value. Finally, the resulting expansion ratio should not exceed 4:1 nor should it
be less than 0.5:1.
Contraction Reach Length. Cross section 20.251 is located where the flow lines are parallel to the main channel. To
evaluate this location, the equation shown below will be utilized. This equation is used when the floodplain scale
and dischargers are significantly different than those used in the regression analysis and is:
where
= contraction ratio
= contraction reach length
= discharge conveyed by the two overbanks at cross section 20.251, cfs
= total discharge, cfs
= Manning's n value for the overbanks at section 20.251
= Manning's n value for the main channel at cross section 20.251
From the analysis of the 600 cfs profile, the flow in the two overbanks at cross section 20.251 was 34.54 cfs and the
Manning's n values for the overbanks and main channel at cross section 20.251 are 0.10 and 0.04, respectively.
Substituting these values into the equation above yielded a contraction ratio of 0.68. The standard error for this
equation is 0.19, which yields a CR range from 0.49 to 0.87 to define the 68% confidence band. This range is less
than the range obtained from Table B-2 of 0.8 - 1.5. Using the contraction ratio of 0.8 and an average obstruction
length of 60 feet, the contraction reach length is:
Expansion Coefficient. The expansion coefficient can be obtained from the equation below:
where
= coefficient of expansion
= hydraulic depth (flow area divided by the top width) for the overbank at cross section 20.208*
= hydraulic depth for the main channel at cross section 20.208*
From the analysis of the 600 cfs profile, the hydraulic depths for the overbanks and the main channel at river station
20.208* are 0.58 and 5.66 feet, respectively. Substitution of the values into this equation yields an expansion
coefficient of 0.07. This is the median value and the range of ± 0.2 defines the 95% confidence band for the equation
above. For this example, a value of 0.5 was used. The value of the expansion coefficient is generally larger than the
value used for the contraction coefficient so the value of 0.5 will remain as the selected value. The modeler can
perform a sensitivity of this coefficient by changing this coefficient and performing subsequent analyses. For this
example, the change in velocity head from cross section 20.227 to 20.208* was only 0.06 feet. Since this value is
small, a change in the expansion coefficient will only reflect a minor change in the resulting water surface elevation.
Contraction Coefficient. From the research document \[HEC-1995\], the contraction coefficient is obtained by first
determining the relationship:
where
= culvert opening width, ft (m)
= total floodplain width, ft (m)
From Table B-3 in Appendix B of the Hydraulic Reference Manual, the recommended contraction coefficient range is
0.3 - 0.5. The value selected for this example was the minimum value of 0.3.
To analyze this further, the modeler can view the profile output table to determine the depths of flow at the
upstream side of the culverts. From the main program window, select View, Profile Table, Std. Tables, and then
Culvert Only. This will display the table as shown in the figure below.
In the this figure, the first two rows are for the flow through the box and circular culverts for the flow of 250 cfs. Then
the second set of rows are for the flow of 400 cfs, and finally the last two rows are for the flow of 600 cfs. For both
the box and the circular culverts, the upstream inside top elevation was set at 31.1 feet. By comparing this value to
the upstream water surface elevation (W. S. US.), it is determined that both of the culverts were submerged at the
entrance for the second and third flow profiles because the upstream water surface elevations were 31.83 and 33.22
feet, respectively.
Additional analysis of the table reveals that outlet control was the type of flow occurring through the culvert for all 3
flow profiles, because the outlet energy grade line was greater than the inlet control energy grade line. The flow
through each culvert, for each flow profile, is shown in the column with the heading "Culv Q." For example, for the
first flow profile (250 cfs), there was 74.73 cfs flowing through the box culverts (37.36 cfs through each barrel) and
175.27 cfs flowing through the circular culverts (87.64 cfs through each barrel). This totals 74.73 + 175.27 = 250 cfs
and accounts for the total flow rate for the first flow profile. Finally, there are no values in the weir flow column of
the table, which signifies that weir flow did not occur for these flow rates. This was an initial goal to develop a
culvert system that produced no weir flow. The modeler can now determine the available freeboard on the
upstream side of the roadway embankment and adjust the sizes and shapes of the culvert as deemed necessary.
For a more detailed analysis of each culvert, select View, Cross Section Table, Type, and then Culvert. Toggle to
profile 3 (600 cfs) and select the "Box" as the Culvert ID. This will display the table as shown in the figure below.
This table displays additional information for the selected culvert and flow profile such as the length of the culvert
flowing full and the energy losses. The modeler can compare these losses with values obtained for other entrance
and exit loss coefficients. Finally, as displayed on the previous table, it can be seen that there is no weir flow
occurring at this river station. The modeler can then toggle to the circular culvert at this river station by selecting
the Culvert ID field. If additional culverts were located at other river stations, the data for these culverts can be
viewed by selecting the appropriate river station.
Finally, the location of the ineffective flow areas will be reviewed. From the main program window, select View and
then Cross Sections. Toggle to river station 20.238 and this will display the cross section as shown in the figure
below. In the figure, it can be seen that the water surface elevations did not exceed the elevations set for the
ineffective flow areas. Additionally, toggle to river station 20.227 and it can be seen that the water surface elevation
is lower than the elevations set for the ineffective flow areas at this cross section. Since weir flow did not occur, the
ineffective flow areas set for cross sections 20.238 and 20.227 should reasonably reflect the actual flow conditions.
7.4 Summary
This example demonstrated the use of HEC-RAS to analyze a river reach that contained a multiple culvert opening.
The multiple culverts consisted of two culvert types (circular and box), each with two identical barrels. After the
flow analysis was completed, the expansion and contraction reach lengths were evaluated and the ineffective flow
areas were reviewed. Finally, the various output features (flow profiles, cross section type tables, and profile tables)
were displayed to show the features available for a review of the output.
For the current example, a flow relief bridge was added at river station 5.40. To add this bridge, first the bridge deck
and roadway data were entered by selecting the Deck/Roadway icon. This activated the Deck/Roadway Data
Editor as shown the bottom figure. The additional high and low cord information for the second bridge was then
added to provide for an opening for the second bridge. This was accomplished by scrolling down to the end of the
existing data and then entering the stations and elevations for the high and low cords of the second bridge. The
data appears in the bottom figure, showing the bridge opening from X-coordinates of 960 to 1240, with the low and
high cords of 215.7 and 216.93 feet, respectively. (These values are the same as used for the first main bridge.) After
the data were entered, the editor was closed. It should be noted that the distance, width, weir coefficient, and other
parameters of the Deck/Roadway Data Editor will be the same for the second bridge as for the first bridge.
This completed the input required for the bridge geometry and the OK button was selected at the bottom of the
editor. Next, the culvert geometry was entered.
Finally, the user should highlight the piers and bridge openings from the View menu of the Bridge/Culvert Data
Editor. Then, zoom in to view each of the bridge openings and the culvert barrels. This will enable the user to view
any inconsistencies that may have developed during the data entry. This completed the data entry for the culverts.
Next, the openings were defined for the multiple flow analysis.
First, the order of the openings are established from left to right as looking in the downstream direction. For this
example, there were three openings: a culvert group, the main bridge opening, and then the relief bridge opening.
To enter the data, the field under "Opening Type" and adjacent to number 1 was selected. Then, the Culvert Group
icon was selected, because the first opening type was a culvert. This placed the description "Culvert Group" in the
first row, under the heading "Opening Type."
Next, the upstream left and right stagnation limits were entered for the culvert group. These left and right stations
will be used to determine where the flow separated between the culvert group and the bridge. For this example, the
left and right stagnation limits were established as 0 and 390. Additionally, the downstream left and right
stagnation limits were also entered as 0 and 390. The downstream stations define the X-coordinate limits where the
flow from consecutive openings will rejoin.
Next, the field under "Opening Type" and adjacent to row number 2 was selected. Then, the Bridge icon was
selected and the description "Bridge" appeared as the second opening type. The left and right stagnation limits
were entered for the bridge as 310 and 880 for both the upstream and downstream side. Finally, the third opening
was selected as "Bridge" and the data entered as shown in the figure above.
As can be seen in the figure above, the right stagnation limit for the culvert group was set at 390 and the left
stagnation limit for the main bridge was set at 310. This created an overlap area from 310 to 390. Since this overlap
exists, the program will determine the actual location of the stagnation point (the flow separation point) between
the culvert and the main bridge. By entering the data in this fashion, the stagnation point was then allowed to vary
for each flow profile, within the limits from 310 to 390. Conversely, if the right stagnation limit of an opening
coincides exactly with the left stagnation limit of the next opening, then the user has defined a specific stagnation
point and this point will be fixed for all flow profiles. For a conveyance type opening, in the current version of the
program, a fixed stagnation point must be used on both sides of the opening. For a further discussion on multiple
opening analyses, refer to Chapter 6 of the User's Manual and Chapter 7 of the Hydraulic Reference Manual.
Once all of the data were entered into the Multiple Opening Analysis Data Editor, the OK button was selected.
Then the locations of the stagnation limits can be viewed in the display of the Bridge/Culvert Data Editor as shown
in . The limits selected for the stagnation points for each opening are shown above the two cross section plots.
For guidance on selecting the stagnation limits, there are two main objectives. First, there is a physical limit. This
implies that there might exist a physical attribute of the openings on the cross section that can be used to
determine the stagnation limits. For example, the left stagnation limit for the culvert was set at the left edge of the
cross section. Additionally, the right stagnation limit for the culvert cannot be placed in the main bridge opening.
This would imply that the flow in front of the bridge would go over to the culvert and this is not a practical.
Therefore, a physical right stagnation limit for the culvert would be the left side of the main bridge opening.
Secondly, for guidance on selecting the stagnation limits, there exists practical limits. For example, the flow in the
cross section in between the culvert and the main bridge must separate at some point to travel into either opening.
Since the bridge opening is larger than the culvert opening, more of the flow will probably go towards the bridge
opening. Therefore, the right stagnation limit for the culvert should not be set all the way over to the left edge of the
main bridge opening. The right limit for the culvert should be located at a practical limit before the bridge opening.
This practical limit can be determined by analyzing the amount of conveyance in the area between the culvert and
the bridge, and then developing a reasonable estimate of where the stagnation limit should be placed.
Finally, it is recommended that the user allow the stagnation points to migrate (where possible) rather than enter
specific fixed stagnation points. This is important when evaluating several flows during the same run. Conversely, if
the stagnation limits are allowed to migrate over a large distance, the program may experience difficulty in
converging to a solution. When this occurs, the overlaps of the stagnation limits should be reduced.
below.
Four ineffective flow area blocks were established at this river station. These were located: to the left of the culvert
group, in between the culvert group and the main bridge, in between the main bridge and the relief bridge, and
finally to the right of the relief bridge. It should be noted that during a multiple opening analysis, the ineffective flow
areas need to be described using the block method, not the normal method. The block method allows the user to
enter ineffective flow areas in between the openings. Finally, ineffective flow blocks were defined at river station
5.39 (located downstream of the multiple opening). These ineffective blocked areas are shown as connected green
lines on the Bridge/Culvert Data Editor (seen in Bridge Geometry in Example 5).
The ineffective flow areas could have been entered along with the cross section geometry (X and Y - coordinates,
etc.); however, for this example, the locations for the culvert group and relief bridge were established first. Then,
the locations of the ineffective flow areas were more readily ascertained from the locations of the culverts and
bridges.
where
= main bridge opening width, ft
= floodplain width contributing to flow through the main bridge, ft
= slope, ft/mi
= Manning's n value of the overbanks
= manning's n value of the overbanks
= average obstruction reach length for the floodplain width contributing flow through the bridge, ft
The main bridge opening width, b, is obtained from the bridge geometric data and was determined to be 200 feet.
The floodplain width for this scenario, B, will be the width of flow that contributes to the main bridge opening and is
the distance from the left stagnation point to the right stagnation point of the main bridge opening. Since these
stagnation points were entered as floating values, an approximate location was assumed for this initial
determination of the expansion reach length and the value was estimated as 540 feet. Additionally, the average
length of the side obstructions, was estimated at 170 feet. This value was determined by only considering the
floodplain width between the left and right stagnation points for the main bridge. Finally, the b/B ratio and the
slope of the river reach are used to determine an initial estimate of the expansion ratio from the "Ranges of
Expansion Ratios," Table B.1 in Appendix B of the *Hydraulic Reference Manual{*}. From the table, the ER was found
to range from 1.3 - 2.0. An average value of 1.7 was used and this resulted with an expansion reach length, , of:
For this example, the expansion reach length is the main channel distance from cross section 5.39 to cross section
5.24* which equals 778 ft. Therefore, an additional cross section was added at a distance of 290 feet downstream
from cross section 5.39 based on the estimated expansion reach length as determined above.
To insert the additional cross section, field data should be used. If this is not available, then the user can utilize the
interpolation routines. Then, the interpolated cross section should be compared with the existing geometry and
topographic maps. For this example, the interpolation method was utilized to obtain a cross section at a distance of
290 feet downstream from river station 5.39. To perform this interpolation, it would be necessary to interpolate
between river stations 5.39 and 5.24*. However, the program will not allow for interpolation between an existing
river station and a previously interpolated section. Therefore, the interpolation would be required from river station
5.39 to river station 5.13. Since this was a very long river reach, an alternative approach was employed by obtaining
data from the USGS Atlas No. HA-601 (This atlas provided the data for Example 2.). From the atlas, the data for river
station number 5.29 was entered into the existing geometry data file. Then, an interpolation was performed
between river stations 5.39 and 5.29.
The interpolation procedure was performed by opening the Geometric Data Editor, then selecting the Cross
Section icon, Options, then Add a new Cross Section. River station 5.29 was entered as the new location for the
cross section and the data for river station 5.29 were entered, with the downstream reach lengths being from river
station 5.29 to 5.13*. Since river station 5.29 was added, the existing reach lengths for river stations 5.39 were
adjusted to 320, 500, and 580 feet for the LOB, main channel, and ROB, respectively. Then, from the Geometric
Data Editor, the following was selected: Tools, XS Interpolation, Between 2 XS's. The upper river station was set
to be 5.39 and this caused the lower river station to be 5.29. A maximum distance of 100 feet was entered and the
interpolation was performed. This yielded 4 new river stations between 5.39 and 5.29, each 100 feet apart. The
interpolated river stations of 5.37*, 5.35*, and 5.31* were deleted and the program adjusted the reach lengths
accordingly. This produced the river station 5.33*, located 300 feet downstream from river station 5.39. For further
discussion on the interpolation procedures, refer to Chapter 6 of the User's Manual and Chapter 4 of the Hydraulic
Reference Manual. The goal was to obtain a river station 290 feet downstream from 5.39, and this distance of 300
feet was determined to be appropriate for an initial location.
For the current example, the contraction reach length is the distance from river station 5.49* to river station 5.41, a
main channel distance of 478 feet. Therefore, a new cross section, located a distance of 190 feet upstream from
river station 5.41 was inserted. The data for this new cross section was obtained from the USGS atlas, namely river
mile 5.44. This river station is located 170 feet upstream from river station 5.41 and was considered appropriate for
an initial estimate of the contraction reach length. Finally, the reach length values for river station 5.61 were
adjusted to account for the inclusion of the new cross section.
Editor, and then selecting Tables, and Reach Lengths. This table is shown in the figure below. Similarly, the
coefficients can be viewed by selecting Tables and then Coefficients.
range of data used to develop the regression equations. The alternate equations, provided by the document, also
produced inconsistent results. The regression equations were developed based on single bridge opening data sets.
Therefore, they may not apply for a multiple opening bridge. The modeler should always use engineering judgment
to the results obtained from these equations.
For this example, the reach lengths estimated at the beginning of the example will be used for the final analysis.
These reach lengths were based on average expansion and contraction ratios for the main bridge opening. It was
considered that these ratios provided a reasonable basis for estimating the expansion and contraction reach
lengths. It should be noted that these reach lengths were significantly less than the reach lengths used in Example
2, with just the main bridge opening. With the multiple openings, the flow does not have to contract and expand as
much as it would for a single opening.
As shown in the figure above, the rows in the table are divided into three groups, one for each profile. The first three
rows are for the first flow of 5000 cfs, the second group is for the flow of 10000 cfs, and the third group is for the flow
of 14000 cfs. The second column in the table displays the river stations and the type of the multiple opening, in the
order as they were entered in the Multiple Opening Analysis Data Editor. The third column displays the total flow
rate through each of the opening types. For example, during the first flow profile (5000 cfs), there was 106 cfs
flowing through the culvert group, 4448 cfs through the main bridge, and finally 446 cfs through the relief bridge.
The sum of these values equals 5000 cfs.
The fifth column displays the calculated upstream energy gradeline elevation for each opening. During the multiple
opening analysis, the program performs an iterative procedure to balance the upstream energy for all the openings.
To do this, the program starts with an initial flow distribution through each opening and then calculates the
upstream energy for each opening. If the energy values are within a specified tolerance (0.05 ft), then the solution is
final. If the energy values for each opening are not within the specified tolerance, then a new flow distribution is
estimated and the procedure is repeated, up to a maximum of 30 iterations. To determine the ranges for the flow
distribution, the stagnation limits are used. For a further discussion on the multiple opening solution scheme, the
user is referred to Chapter 6 of the User's Manual and Chapter 7 of the Hydraulic Reference Manual.
In reviewing the calculated energy gradeline values, it can be seen that for the three flow profiles, the energy
gradeline elevations were balanced within the specified default tolerance of 0.05 feet. With these energy gradeline
values, the program then calculates an upstream water surface elevation for each opening by subtracting the
velocity head from the energy gradeline. These water surface elevations are displayed in the sixth column.
Finally, the last two columns of the multiple opening profile table display the calculated left and right stagnation
points. These stagnation points are the flow distribution limits that were determined by the program in order to
balance the upstream energy gradelines. For the first flow profile, it can be seen that the flow distribution was
established from upstream X-coordinates of 73.6 to 310 for the culvert, from 310 to 700 for the main bridge, and
from 700 to 1287.02 for the relief bridge.
For the culvert, the program used a left stagnation point of 73.6 feet. This value is the X-coordinate where the
culvert water surface elevation of 213.22 intersected the left side of the cross section. The right stagnation point for
the culvert was determined to be 310 feet. Therefore, in order to balance the upstream energy gradelines, the
program determined that the flow rate through the culvert required an upstream conveyance that encompassed
cross section station 73.6 to 310.
Similarly, for the first flow profile, the program determined that the flow through the main bridge opening would be
from an upstream X-coordinate of 310 to 700 and the flow rate through the relief bridge would be from upstream X-
coordinates of 700 to 1287.02. The value of 1287.02 is the X-coordinate where the relief bridge water surface
elevation of 213.26 coincided with the right edge of the cross section.
Upon reviewing the left and right stagnation points for the third flow profile, it can be seen that the stagnation point
between the culvert and the main bridge is a value of 310 feet as opposed to the previous value of 390. Therefore, in
order to balance the upstream energy grade line elevations for the three openings, a greater portion of the total
flow was required to travel through the main bridge opening during this higher flow event. Similarly, the left culvert
stagnation point and the right relief bridge stagnation point have changed to reflect that the higher water surface
elevations are approaching the limits of the cross section width.
Finally, the Multiple Opening Profile Table provides additional information such as the total cross sectional flow
area and the top width of the effective flow. The descriptions for the column headings will appear in the dialog box
at the bottom of the table when an entry in the specific column is selected.
As discussed previously, the energy gradeline elevations shown in the Multiple Opening Profile Table are the
energy values for each of the openings, upstream of the openings. However, since the program is a one dimensional
model, the program must determine only one energy grade line value to use at the upstream cross section. To
determine this energy value, the program uses a flow weighting method to determine the average energy at the
upstream cross section. Then, this average energy is used to calculate one water surface elevation for the entire
upstream cross section. To review the energy gradeline and water surface elevations that were used as the final
answer for the upstream cross section, from the main program window select View, Cross Section Table, Type,
and then Cross Section. Toggle to river station 5.41 (the upstream cross section for this example) and select profile
3. This will display the table as shown in the figure below.
For the third flow profile, the program calculated that the upstream energy gradeline elevations for the culvert,
main bridge, and relief bridge were 216.55, 216.58, and 216.57 feet, respectively, as shown in the Multiple Opening
Profile Table. However, as described above, the program can only use one energy grade line elevation at the
upstream cross section. The flow weighted upstream energy gradeline elevation used for cross section 5.41 is
shown in the figure below to be 216.58. Additionally, the program had calculated the upstream water surface
elevations for each of the three openings as 216.55, 216.14, and 216.55, respectively, as shown in the Multiple
Opening Profile Table. However, the program can only use one water surface elevation at the upstream cross
section. To obtain this one water surface elevation, the program used the average energy value, subtracted the
average velocity head for the entire cross section of 0.38, and determined that the upstream water surface elevation
was 216.20, as shown in the figure below. This water surface elevation is calculated from the average energy of the
upstream section. The actual water surface elevations at the bridge are more likely reflected by the values as shown
in the Multiple Opening Profile Table.
As a final component of the analysis, the flow distribution for the three openings was compared to the flow in the
LOB, main channel, and ROB at the river stations upstream and downstream of the openings. From the Multiple
Opening Profile Table (see the figure above), the flow for the third profile through the culvert, main bridge, and
relief bridge were 566.98, 9950.34, and 3482.68 cfs, respectively. Since the culvert group is located in the LOB, the
main bridge is located over the main channel, and the relief bridge is located in the ROB, these flow values should
approximate the flows in the LOB, main channel, and ROB at river stations 5.41 and 5.39. As shown in the figure
below, the flows for the third profile in the LOB, main channel, and ROB were 676.13, 10182.13, and 3141.74 cfs,
respectively. Similarly, the flows at river station 5.39 were 705.59, 9931.67, and 3362.74 cfs. These flow distributions
are similar and reflect the transition of the flow rates across the cross sections.
8.8 Summary
The geometry of Example 2 was modified to prevent weir flow from occurring over the main bridge decking. This
was accomplished by defining a culvert group and a relief bridge on the left and right side of the main bridge,
respectively. These additional openings lowered the water surface upstream of the bridge so that weir flow did not
occur. The figure below shows a comparison of the third water surface profile for the multiple opening and for the
original river reach with just the main bridge. To develop the figure, select View and then Water Surface Profiles
from the main program window. Then select Options and Plans and choose both of the plans. Finally, Options and
Profiles were chosen and only the third profile was selected to be plotted. This resulted in the display shown in the
figure below, which clearly shows that the modified conditions decreased the upstream water surface profile so
After the analysis was performed, the user can review the flow parameters for the culvert and the relief bridge and
adjust the size of these openings to develop the most practical alternative to prevent the weir flow from occurring.
To compare the results in tabular format, from the main program window select View and then Profile Table.
Standard Table 1 was then selected from the Std. Tables menu. Finally, both of the plans and only the third profile
were selected from the Options menu and this resulted in the table as shown in the figure below.
By comparing the values in the table, the user can obtain required information to assist in the determination of any
changes that may be necessary. For example, as can be seen in the figure below, the change in the water surface
elevation at river station 5.41 was 1.23 feet (from 217.43 to 216.20 feet). The modeler can use this information to
determine if any additional decrease is necessary.
program must first determine a natural profile with no encroachments. This base profile is typically computed using
the one percent chance discharge. The computed profile will define the floodplain, as shown in the figure below.
Then, by using one of the 5 encroachment methods, the floodplain will be divided intotwo zones: the floodway
fringe and the floodway. The floodway fringe is the area blocked by the encroachment. The floodway is the
remaining portion of the floodplain in which the one-percent chance event must flow without raising the water
surface more than the target amount.
For this example, the following procedure was employed to perform the encroachment analysis:
• Determine the 100-year flood profile
• Method 5 optimization procedure
• Method 4 - with 3 target depths
• Method 4 - with 1 target depth
• Method 1 - final delineation of floodway
• Review of floodway delineation sketch
To perform the floodplain encroachment analysis for this example, the first step was to develop a model of the river
reach that would compute the 100-year flood profile. This model must be developed and calibrated to the fullest
extent possible because it defines the base flood elevations and all subsequent calculations will be based upon this
profile. This was accomplished in Example 2.
After the base profile was computed, the Method 5 procedure was chosen as an initial attempt to calculate the
encroachments. Method 5 will typically calculate reasonable encroachment stations for "well behaved" streams.
That is, for streams that exhibit minor changes in cross section geometry and have small losses due to bridges and
culverts. If the river reach has abrupt changes in geometry or orientation, contains a flow controlling structure, or
the encroachments encounter the main channel bank stations, then Method 5 may produce erratic results at these
locations.
If Method 5 produces inconsistent results, then Method 4 may be utilized. Method 4 is frequently used for a
floodplain encroachment analysis. The initial approach is to use this method with several target water surface
increases.
Finally, after the results are obtained from the Method 4 analysis, Method 1 was used to further refine the
encroachment stations. The modeler should sketch the floodway on a topographic map to visually inspect the
floodway and allow for smooth transitions. The computed floodway is considered preliminary, in that the
regulating community must approve and adopt.
Each of these steps, as performed for this example, are discussed in detail in the following sections.
station information; and method and target values. The following sections describe each of the data entry items.
For a further discussion of the data entry procedure, the user is referred to Chapter 9 of the User's Manual.
The global information is applied to all of the reaches and cross sections that will be selected for the analysis. The
first input for the global information is the Equal Conveyance Reduction box. If this box is selected, then the
program will encroach by simultaneously removing an equal amount of conveyance on both sides of the main
channel. As the amount of conveyance is removed, if one of the encroachments reaches the main channel bank
station (or the offset), then the program will continue to encroach on the other side until the target values are
obtained or until the encroachment on the other side reaches the main channel bank station (or the offset). If the
equal conveyance box is not selected, then the program will encroach by maintaining a loss of conveyance in
proportion to the distribution of the natural overbank conveyance. The equal conveyance reduction applies to
Methods 3, 4, and 5; and, for this example, the option was selected.
The next global information items are the Left bank offset and the Right bank offset. These offsets limit the
distance of the encroachments. Without an offset, encroachments can go up to the bank stations, eliminating the
entire overbank. For this example, the offsets were set to be 10 feet for both the left and right bank. Therefore, the
limit of encroachment was 10 feet to the left of the left channel bank station and 10 feet to the right of the right
channel bank station.
Reach and River Station Information. The next items to select in the encroachments editor are the River and
Reach that will be analyzed. For this example, there is only one river and reach: the "Kentwood" reach on the river
"Beaver Creek." Next, the Starting River Station and Ending River Station were entered. The analysis was set to
begin at river station 5.99 and end at river station 5.00, the entire river reach. Finally, the Profile was selected. For a
Method 5 analysis, only 2 profiles are necessary. The first profile will be used to determine the base water surface
profile and cannot be selected in the encroachments data editor. The second profile is selected to be used for the
Method 5 analysis.
Method and Target Values. Method 5 was selected for this analysis. When a method is selected, the data entry
fields required for that method will appear immediately under the method field. For Method 5, the fields Target WS
change (ft) and Target EG change (ft) appeared. A target water surface change of 1.0 foot and a target energy
change of 1.2 feet were entered. Typically, the energy target will be chosen to be slightly greater than the water
surface target. The energy target will act as an upper limit during the iterations to prevent the encroachment from
getting very large. The next step is to select the range of the river reach that will be analyzed, for the target values
that were entered. For this example, the Set Selected Range button was selected and this applied Method 5 and
the chosen targets to the selected range of river stations 5.99 to 5.00 for the Beaver Creek reach. By selecting this
button, the fields in the table were filled with the corresponding method and values. In the table, since Method 5
was selected, the heading Value 1 corresponds to the target water surface change and Value 2 corresponds to the
target energy change. When other methods are selected, the value 1 and value 2 columns will represent the specific
data input items for the method. (Note: Some methods only require 1 value.) As a final note, the user can now edit
the data table and change any of the methods and value items for any specific river station. For this example, the
table was not edited at this time.
The OK button was then selected at the bottom of the encroachment editor. This prompted the Steady Flow
Analysis Window to appear. Then the geometry file "Existing Conditions" and the steady flow file "Base + 1 ft
Target Depth" were saved as the plan "Method 5 Encroachment." Finally, the Short ID was entered as "M5" and the
COMPUTE button was selected for the subcritical flow analysis.
different flow depth at this river station. Often, the resulting depth of flow will be greater than the target
increase in water surface elevation. Therefore, the first run with Method 4 is typically applied with several
target values, usually smaller than the maximum increase.
6. The program then moves to the next river station and steps 2 through 5 are repeated. This process continues
until the last river station is evaluated.
For this review of the first run of Method 4, the main concern is with the change in water surface elevation. The
encroachment process and resulting water surface elevations were determined using the 6 steps as described
previously at the beginning of this section. For river station 5.00, the change in water surface was 1.0 foot for all of
the target values because the downstream boundary conditions were set to be 1 foot higher than the base profile.
At river station 5.13, the changes in water surface were 0.94, 0.97, and 1.01 feet for the target values of 0.8, 0.9, and
1.0 feet. This shows that the actual resulting water surface is generally greater than the target values. Further review
of the table shows the various resulting water surface elevations for the target depths.
To continue the encroachment analysis process, a Method 4 analysis was again calculated, however this was
performed with only 1 target depth for each river station. The target depths that were used in the subsequent
analysis were the target depths that resulted in a water surface change as close to 1 foot without exceeding 1 foot.
For example, at river station 5.13, the resulting change in water surface that is as close to 1 foot as possible without
exceeding 1 foot is 0.97 feet. This value was obtained from a target increase of 0.9 feet from the third profile.
Therefore, for the next trial, a target value of 0.9 was used at river station 5.13. Likewise, a target of 0.80 was used at
river station 5.24* because this target yielded a change in water surface of 0.94 feet. The following table shows the
values that were used for the subsequent analysis.
It should be noted that for the bridge river station (5.40), the user can only enter one target value into the
encroachment data editor. Therefore, only one target value appears in the table above. Additionally, since the
pressure/weir option is being used as the high flow analysis method, the program will use the encroachments
calculated at river station 5.39 as the encroachments for river stations 5.40 and 5.41. This will be discussed further
in the proceeding section and additional discussion is provided in Chapter 9 of the Hydraulic Reference Manual.
A new plan was created with the geometry data file "Existing Conditions" and the steady flow data file "Base + 1 ft
Target Depth." These files were saved as the plan "Method 4 Encroachment - Trial 2." The Short ID was entered as
"M4 - Trial 2" and the COMPUTE button was then selected.
results in the table, the change of water surface for river stations 5.00 and 5.065* are both 1.00 feet. For river station
5.13, the change in water surface is 1.01 feet. This value was obtained with a target value of 0.9 feet. Therefore, to
decrease the resulting water surface elevations at this river station, the target value for this river station must be
decreased and a subsequent Method 4 analysis performed. This iterative process of changing the target values and
reviewing the output can determine the floodway that will result in a change of water surface less than 1 foot.
When performing a subcritical flow analysis, the user should begin this iterative process at the downstream cross
section and work upstream. Additionally, the user should not attempt to adjust a large quantity of target values at
the same time. For this example, the target values for the river stations below the bridge (5.00 through 5.29) were
adjusted first. Then the bridge section was analyzed and finally the upstream river stations were adjusted. This
procedure will allow the user to focus on specific river sections and adjust these target values before moving onto
the further upstream river stations.
As an additional review of the encroachment 1 table, it can be seen that the left and right encroachment stations
that were calculated for river station 5.39 were 440.00 and 1144.54. These encroachments were then used for river
stations 5.40 and 5.41. This is a default method for the program since the pressure/weir method was used for the
bridge high flow analysis. If the energy method had been used, then the program would allow for separate
encroachment stations throughout the bridge.
Finally, the left encroachment station of 440 at river station 5.39 is 10 feet to the left of the main channel bank
station. Therefore, the program encroached up to the left offset on this side of the main channel. (This also
occurred at river stations 5.00, 5.065*, and 5.13.) If the offsets had not been used, then the left encroachment would
have continued up to the main channel bank station. When this occurs, additional wetted perimeter will be added
to the main channel. This will cause the conveyance of the main channel to decrease and the total loss of
conveyance at the cross section may be greater than if the encroachment did not encounter the main channel.
Therefore, this may reduce the amount of encroachment on the right side of the channel since an additional loss of
conveyance had already occurred.
To determine the percent reduction of conveyance removed from each side, the Encroachment Table 2 can be
viewed. A portion of this table is shown as the following figure. From the table, it can be seen that for river station
5.00, an approximately equal amount of conveyance was removed from each side of the main channel (15.43 and
15.53 percent for the left and right sides, respectively). This occurred even though the encroachment encountered
the left offset. However, for river station 5.065*, the encroachment also encountered the left offset but the percent
reduction of left and right conveyance is not equal. The percent of conveyance removed at river station 5.065* is
10.79 and 18.89 for the left and right sides, respectively. This implies that after the left encroachment reached the
left offset, the program needed to remove an additional amount of conveyance from the left. Therefore, the
program removed that amount from the right side in addition to the amount required to be removed from the right
side. This caused the percent of conveyance removed to be unequal on both sides of the main channel.
The final target water surface values used and the resulting change in water surface are listed in the table above.
These values were obtained from the Encroachment Table in trial 1 as discussed previously.
In the vicinity of the bridge, a target value of 1.1 foot at river station 5.39 (just downstream of the bridge) was found
to result in the most practical encroachments through the bridge. If a higher target value was used, the program
would encroach further towards the main channel and a slightly higher water surface value would be obtained at
river stations 5.39, 5.40, and 5.41. If the encroachments at the bridge were moved closer to the main channel, then
this would further increase the upstream water surface. The increase in upstream water surface is in effect before
the upstream encroachments are calculated. Therefore, the upstream encroachments are limited from the start
because the upstream water surface elevation is already greater than the natural profile.
As the target value at river station 5.39 was increased beyond 1.1, the program could only provide minor
encroachment distances at river station 5.49*. This created an erratic transition in the floodway at river station
5.49*. As a compensation between the bridge encroachment constriction and initial upstream rise, the target value
of 1.1 foot at river station 5.39 was determined to be most practical.
Additionally, at river station 5.40 (the bridge) and 5.41, the change in water surface was only 0.11 feet. This water
surface elevation only increased slightly due to two factors. First, the rise in water surface upstream of the bridge
(river station 5.41) is controlled by the bridge structure itself, due to the occurrence of pressure and weir flow. The
increase in water surface at section 5.39 did not cause enough submergence on the weir to increase the upstream
headwater. Secondly, the conveyance reduction at 5.40 and 5.41 removed only the weir flow that was occurring.
This reduction of weir flow was not sufficient enough to cause the water surface to rise dramatically. This implies
that the pressure flow was dominate through the bridge opening. The modeler should check the bridge solution to
determine the value of the pressure flow. For this example, the pressure flow through the bridge was 12179.63 cfs, a
major portion of the total flow of 14000 cfs.
It should be emphasized that the Method 4 iterations were continued until a practical floodway was developed.
After the Method 4 procedure was completed, the results were used in a Method 1 analysis which is discussed in the
following section.
the user can fine tune and adjust the encroachments as deemed necessary by adjusting the left and right
encroachment stations in the Encroachment Data Editor. For this example, no further adjustments were made.
To review the output in graphical form, from the main program window select View and then X-Y-Z Perspective
Plots. This will result in the display as shown in the following figure. As can be seen below, the encroachment
stations appear to follow a smooth transition throughout the river reach. However, the user must be aware of the
fact that this plot is based upon the X-coordinates as entered by the user. If the X-coordinates for the cross sections
are not all established from the same left baseline, then the plot may not be accurately portraying the correct
configuration of the floodway. The modeler should sketch the resulting encroachments and floodway on a
topographic map to view the correct alignment of the floodway. At this point, further refinement for the locations of
the encroachments stations should be made.
In addition to the 3-D plot, the user can also view the individual cross section plots to see the location of the
encroachment stations. By using the information from the encroachment tables, the cross section plots, the 3D
plot, and the user developed topographic plot, the encroachment stations should be evaluated for the required
constraints and the transitions of the floodway.
9.8 Summary
To perform the floodway analysis for this example, a Method 5 procedure was first attempted. This procedure can
yield reasonable results for a smooth transitioning river reach. However, for this example, the bridge structure in
the river reach caused difficulties in the Method 5 analysis and the program did not yield reliable results upstream
of the bridge.
To continue the analysis, a Method 4 procedure was employed. First, 3 target depths were used to obtain a first cut
at the encroachment stations. Then, the Method 4 procedure was employed iteratively with one target value for
each river station to further define the encroachments. After the encroachments appeared to be well established,
the results were imported to Method 1 for a final check on the encroachments and to perform any additional
smoothing of the floodway transitions.
When performing a floodway analysis, the general approach is to attempt to encroach on both sides of the water
course without increasing the water surface elevation by some predefined amount. The user should also be aware
of other constraints such as velocity limits and equal conveyance reduction requirements, which may be
constraining the floodway delineation. Additionally, the floodway must be consistent with local development plans
and provide reasonable hydraulic transitions throughout the study reach. These transitions must be determined by
plotting the encroachment stations onto a topographic map. The user should not rely on the 3-D plot provided with
the program due to the constraints of the plot as described previously.
Project Files
As the file types are created, the program will number the extensions in consecutive order starting at 01. The run
and output file extension numbers will always coincide with the extension number of the plan that the program
used to create the files. Finally, the project also specifies the default variables. These variables include the system
of units (English or SI) and the contraction and expansion coefficients. These variables can be changed by the user
through the program interface.
The river reach is defined by 84 river stations. The user can view the geometric data for each river station by
selecting the Cross Section icon from the Geometric Data Editor. The data for each river station is comprised of : a
description; X and Y coordinates; downstream reach lengths, Manning's n values; main channel bank stations; and
contraction and expansion coefficients.
Additionally, a levee was placed on the left side river stations 200 and 300 and ineffective flow areas were
established at river stations: 956 and 1002; 1208 and 1229; 1318 and 1383; 2484 and 2540. Each of these four
groupings are for the four bridge locations in the river reach. Finally, the expansion and contraction reach lengths
for each bridge were estimated using the expansion and contraction ratios obtained from Table B.1 and B.2 in
Appendix B of the Hydraulic Reference Manual.
The focus of this project is to replace the existing bridge structure located at river station 2512. To view the existing
bridge, from the Geometric Data Editor select the Bridge/Culvert icon and toggle to river station 2512. This will
display the Seminary Street Bridge at river station 2512 as shown in the figure below.
The bridge deck and roadway information were entered by selecting the Deck/Roadway icon on the left side of the
editor. Additionally, the Bridge Modeling Approach icon was selected and the appropriate information provided
(These procedures were also applied to develop the other 3 bridges in the river reach).
After all of the geometric data were entered, the geometry file was saved. This was performed from the Geometric
Data Editor by selecting File and then Save Geometry Data As. The title was entered as "Existing Geometry" and the
OK button was selected.
view the profile plot of the river reach, from the main program window select View and then Water Surface
Profiles. This will display the profile as shown in the figure below. The first part of the heading at the top of the
profile is the name of the project: "Napa Cr. Bridge Project - Example 7." The second part of the heading is the name
of the plan: "Existing Plan Data." Finally, the date that the output file was created also appears in the heading.
The profile of the existing conditions analysis shows that the upstream bridge at river station 2512 was overtopped
during the 100-year flow event. It was proposed that a new bridge be installed to replace the existing bridge. The
new bridge would be designed so that the 100-year flow event did not impact the bridge decking. To evaluate this
proposed bridge replacement, the existing conditions geometry file was changed and a new plan was created. This
procedure is outlined in the following sections.
Since the bridge geometry changed at river station 2512, the expansion and contraction reach lengths were
adjusted in the vicinity of the bridge. Then, the ineffective flow areas were reestablished at river stations 2484 and
2540. This concluded the changes to the geometric data, and the geometry file was saved by selecting Save
Geometry Data from the File menu of the Geometric Data Editor.
The figure above shows the profile for both plans. The heading displays the information as discussed previously.
The legend shows that there are two water surface profiles plotted in the figure. The first profile is labeled "WS 100
yr - Existing" and is a solid line. The label "WS 100 yr" refers to the label of the flow data for the first water surface
profile, while the label "Existing" refers to the Plan Short ID that was entered in the Steady Flow Analysis Window.
Therefore, this water surface is for the Existing Condition plan. Similarly, the label "WS 100 yr - Proposed" is the 100
yr water surface profile for the "Proposed" plan data.
Since both plans only had one flow entered, only one water-surface profile can be plotted for each plan. If other
flow profiles had been computed, then the user could also select to plot those profiles.
It should be noted that since the profiles were plotted from the proposed conditions plan, the geometry that is
displayed in the figure above is from the proposed conditions plan. In other words, the bridge at river station 2512
displays the geometry (elevation) of the proposed bridge. If this procedure had been performed from within the
existing conditions plan, then the geometry of the profile would exhibit the existing conditions geometry. In either
case, the water surface profiles for each plan are plotted as calculated for that plan. Finally, the figure above clearly
shows the decrease in the upstream water surface due to the proposed bridge allowing the flow to pass completely
under the bridge.
Additionally, the program will revert back to only displaying the current plan information whenever a new plan or
project is opened.
In reference to the preceding figure, the headings and legend are as described previously. The cross section was
selected from within the proposed geometry plan and therefore the bridge shown is the proposed bridge. The first
water surface is for the existing conditions plan and the second water surface is for the proposed conditions plan.
As an additional note, if a river station only appears in one of the plans, then the table will only display the data for
that plan. This occurs for river station 2549* and 2471*, which are only used in the proposed geometry plan. These
river stations can be viewed by using the down arrows on the right side of the table.
As for the standard table, the river reach and river stations are shown in the first two columns and the Short IDs are
used to identify the plans in the third column. This table shows that for the bridge at river station 2512, weir flow
occurred for the exiting conditions bridge but did not occur for the proposed bridge geometry plan. The data for the
other three bridges in the table are the same for both plans since there was only a change in the geometry at the
upstream bridge.
stations on the plot are aligned according to the configuration as drawn on the River System Schematic.
10.6 Summary
For Example 7, the concept of multiple plan analysis within a single project was discussed. There are several
advantages to using a multiple plan analysis. One of these advantages is the use of a single flow data file for
multiple geometry simulations. This reduces the need for identical flow files. Another advantage pertains to the
analysis of the output. With multiple plans, the user can select to have the data from any number of plans displayed
collectively. This allows for a proficient comparison of the plans.
For this example, the multiple plans were developed by altering the geometry of the river reach. Conversely, a
multiple plan analysis could be composed of plans that relate changes of user-selected coefficients such as energy
loss coefficients, Manning's n values, and entrance and exit loss coefficients. By understanding the concept of
multiple plan analyses, the user can employ a more efficient procedure for analyzing a project.
The next piece of information required was the Length Across Junction. These are the distances from the
downstream river station of Upper Spruce to the upstream river stations of Middle Spruce and Bear Run. In general,
the cross sections that bound a junction should be placed as close to the junction as possible. This will allow for a
more accurate calculation of the energy losses across the junction. These values were entered as 80 and 70 feet, for
the distances to Middle Spruce and Bear Run, respectively.
The last item in the junction editor is the computation mode. Either the Energy or the Momentum method must be
selected. The energy method (the default method) uses a standard step procedure to determine the water surface
across the junction. The momentum method takes into account the angle of the tributaries to evaluate the forces
associated with the tributary flows. For this example, the flow velocities were low and the influence of the tributary
angle was considered insignificant. Therefore, the energy method was selected for the analysis. For a further
discussion on stream junctions, the user is referred to example 10 and to chapter 4 of the Hydraulic Reference
Manual.
After the data were entered for the Tusseyville Junction, the Apply Data button was selected. The down arrow
adjacent to the Junction Name was depressed to activate the second junction at Coburn. At this junction, the
description "Confluence of Bear Run and Middle Spruce" was entered. Next, a length of 70 feet was entered from
Bear Run to Lower Spruce and 85 feet for the distance from Middle Spruce to Lower Spruce. Again, the energy
method was selected and the Apply Data button was chosen before closing the junction editor.
To enter the flow data, a flow value must be entered at the upstream end of each reach. The program will consider
the flow rate to be constant throughout the reach unless a change in flow location is entered. For this example, the
flow will be constant throughout each reach. The three profiles will be for flow values of 300, 800, and 1000 cfs.
These values were entered as the flow rates for Upper Spruce and Lower Spruce.
For the flow rates through Middle Spruce and Bear Run, the user must estimate the amount of flow for each reach.
Then, after the analysis, the user must compare the energy values at the upstream ends of Middle Spruce and Bear
Run. If the energy values differ by a significant amount, then the flow rates through the two reaches must be
redistributed and a second analysis performed. This process will continue until the upstream energies are within a
reasonable tolerance. This procedure implies that the upstream cross sections of Middle Spruce and Bear Run are
located close to the junction. Therefore, the energy value at these two locations should be approximately equal.
For this first attempt at a flow distribution, the values of 170 and 130 cfs were entered for the first profile for Middle
Spruce and Bear Run, respectively. Similarly, flow values of 450 and 350 were entered for the second profile and 560
and 440 for the third profile. After the analysis, the upstream energies for each profile were compared to determine
if the flow distribution was appropriate. This will be discussed in a subsequent section.
(Note: Since there was only 1 geometry file and only 1 flow file, this was not necessary.) Then, File and Save Plan As
were selected and the title "Looped Plan" was entered. The OK button was selected and the plan title appeared
near the top of the steady flow window. Finally, the Flow Regime was selected as subcritical and the COMPUTE
button was selected.
11.7 Summary
As a summary for this example, a river system that contained a loop was analyzed. The flow rates for the branches
of the loop were initially estimated and, after an initial analysis, the upstream energy values for each branch were
compared. Since the initial energy values were not within a reasonable tolerance, the flow rates through each
branch were redistributed and a subsequent analysis performed. This procedure was continued until the upstream
energy values for the two branches were within a reasonable tolerance. By performing the flow distribution and
energy comparison in this manner, it was necessary that the cross sections around the junction were spaced close
together.
To perform the interpolation, the reach of Putah Creek and the starting and ending river stations of 12 and 1 were
selected as shown in the figure above. Then, a maximum distance of 50 feet was entered and the Interpolate XS's
button was selected. This created the interpolated cross sections along the entire river reach. After the
interpolation, each interpolated cross section was reviewed to determine the adequacy of the interpolation. For
additional discussion on cross section interpolation, the user is referred to chapter 6 of the User's Manual and to
chapter 4 of the Hydraulic Reference Manual.
As a final review of the cross section locations, from the Geometric Data Editor select Tables and then Reach
Lengths. This will activate the Reach Lengths Table as shown in the figure below. As can be seen in the table, the
final channel reach lengths are 50 feet for each cross section, except through the bridge.
The central portion of the editor consists of the station and elevation data for the low and high cords of the deck/
roadway. These values were entered as shown. Finally, the bottom portion of the editor consists of the data entry
for the weir flow calculations. For this example, weir flow did not occur and this data will not be emphasized. For
additional discussion on the deck/roadway data, the user is referred to example 2 for bridges and to example 3 for
culverts. The OK button was then selected to exit the editor.
Pier Data. To enter the pier data, the Pier icon was selected from the Bridge/Culvert Data Editor. This activated
the Pier Data Editor. For this example, the bridge geometry consisted of only one pier located at a centerline
station of 150 and a width of 8 feet. This data was entered and then the editor was closed. This concluded the input
of the bridge geometry and the bridge appeared in the Bridge/Culvert Data Editor as shown in the figure below.
As a final note for the bridge geometry, there were no ineffective flow areas defined for the analysis. This was due to
the fact that the flow for the analysis remained in the main channel and the bridge geometry did not create any
appreciable ineffective flow areas. This will become more apparent during the review of the output.
Bridge Modeling Approach. The final component of the geometric data is the entering of the coefficients for the
bridge analysis. This was performed by selecting the Bridge Modeling Approach icon on the Bridge/Culvert Data
Editor. This activated the Bridge Modeling Approach Editor as shown in the figure below.
The first selection is the low flow computation methods. For this example, the energy and momentum methods
were selected. Then, the Highest Energy Answer field was selected. This will inform the program to use the greater
answer of the energy and the momentum methods for the final solution of the low flow analysis. The next selection
was the method for the high flow analysis. For this simulation, the energy method was selected. (A subsequent
analysis was performed with the pressure/weir method and will be discussed later.)
This concluded the geometric input for the analysis and the data was then saved as "Base Geometry Data - Energy."
Next, the steady flow data were entered for the simulation.
When the program calculated the water surface profile for river station 9, it determined that there was a valid
answer for both the subcritical and supercritical flow profiles. The program then compared the specific force of
both of these flow regimes and determined that the subcritical flow had a greater specific force. This implies that
the flow at river station 9 was subcritical and that a hydraulic jump developed upstream of this river station. This
can be seen in figure above to have occurred between the main channel stations of 710 and 760 (river stations 9 and
9.5*). This hydraulic jump occurred because the cross section geometry in the vicinity of the bridge caused a
constriction of the flow and a backwater was created upstream from the bridge. This backwater created a
subcritical profile immediately upstream of the bridge and a hydraulic jump was necessary for the flow to transition
from supercritical to subcritical.
To view the calculated values of the water surface elevations, critical depth, and energy gradeline, the Standard
Table 1 profile table was activated from the main program window by selecting View, Profile Summary Tables, and
then Std. Table 1. This table is shown in the figure above, and shows a water surface elevation of 1280.84 and a
critical water surface elevation of 1281.58 at river station 9.5*. This water surface elevation is less than the critical
depth and implies a supercritical flow regime. At river station 9, the water surface elevation of 1282.19 is greater
than the critical water surface of 1281.03 and implies a subcritical flow regime. The modeler can use this table to
determine all of the computed values that are displayed on the profile plot as shown in the figure preceding the one
above.
Referring back to the figure before, as the calculations continued through the bridge, Class B low flow was found to
occur since the flow did not encounter the low cord of the bridge and the flow passed through critical depth under
the bridge. For Class B low flow, the program will set the water surface at critical depth at either the upstream
inside or downstream inside cross section of the bridge. The program will calculate the specific force for critical
depth at both of these sections and set the flow at critical depth at the section that is the most constricted and has
the greater specific force. For this example, the flow was set at critical depth at the bridge upstream inside cross
section, as shown in the second most recent figure. Additionally, if the specific force of both cross-sections are
approximately equal, then the program will use the location entered by the user. The location can be selected from
within the Bridge/Culvert Data Editor by selecting Options and then Momentum Class B defaults.
Finally, a supercritical profile continued downstream of the bridge to the last downstream cross section. The
downstream boundary condition had been set at critical depth. However, the program determined that the
supercritical flow solution at the downstream end had a greater specific force than the boundary condition and
used the supercritical flow answer.
12.5.2 Water Surface Profiles for Subcritical and Supercritical Flow Analyses
To perform the analysis, the mixed flow regime had been selected. If the user had selected a subcritical or
supercritical flow regime for the analysis, the output would have reflected various warnings and notes intended to
alert the user of possible inconsistencies with the results. For example, if the user had selected a subcritical flow
regime, the water surface profile would have appeared as shown in the figure below.
As can be seen in the figure above, the water surface coincided with the critical depth line for the majority of the
river stations. A review of the Summary of Errors, Warnings and Notes would reveal the repetition of the warning:
"During the standard step iterations, when the assumed water surface was set equal to critical depth, the
calculated water surface came back below critical depth. This indicates that there is not a valid subcritical answer.
The program defaulted to critical depth." This warning is issued when the user has requested a subcritical flow
analysis but the program could not determine a subcritical flow depth at the specified cross section. Since a
subcritical solution was not possible, the program used critical depth at this location and continued on with the
calculations. This warning may be associated with too long of reach lengths between cross sections or the fact that
the flow analysis should be performed in the supercritical or mixed flow regimes.
If the user had selected to perform a supercritical flow analysis, the water surface profile would have appeared as
shown in the figure below. As can be seen in the figure, there is an inconsistent drop in the energy gradeline
immediately upstream of the bridge. A review of the warnings at the upstream inside bridge cross section revealed
that the energy equation could not be balanced within the specified number of iterations and the program
defaulted to critical depth at this location. The user should perform the computations in the mixed flow regime to
determine if a subcritical flow profile exists in the river reach. The program can only provide for both subcritical and
supercritical flow answers when the mixed flow regime is selected.
The analysis of the river reach in the subcritical and supercritical flow regimes are not provided as plans for this
example. They were computed and presented to show what would develop if these flow regimes had been selected.
However, the modeler can readily select these other flow regimes and execute the program to observe the output.
For an additional discussion for the descriptions of the warnings, errors, and notes, the user is referred to example 1
and to chapter 10 of the User's Manual.
The first two columns in the table show the reach and river station for the bridge location. For this example, the
bridge is located at river station 7. The third column shows the energy gradeline elevation that was used as the final
answer for the analysis, 1283.00 ft. The fourth column displays the water surface elevation (1282.21 ft) that
corresponds to the energy gradeline in column three. The fifth column displays the bridge method that was
selected as the final answer. Finally, the sixth and seventh columns show the results of the energy and momentum
methods since these two methods were selected to be computed. As can be seen in the table, the energy low flow
method produced a result of 1283.07 ft for the energy gradeline. The momentum method produced a result of
1283.00 for the upstream energy gradeline, and was selected as the final bridge solution. The reason the
momentum method was selected, instead of the highest upstream energy answer, is that the bridge is in a class B
flow regime (passing from subcritical, through critical, and going supercritical) inside of the bridge. For this type of
flow regime, HEC-RAS uses the Momentum equation as the best answer, due to the rapidly varying flow (i.e. the
energy method does not perform as well in rapidly varying flow regimes – where flow passes through critical depth,
and may go through a hydraulic jump). The only time the momentum method is not used for Class B low flow is if it
fails to converge on an answer, then the energy method will be used..
This completed the review for the energy method analysis. As a final review of the water surface profile as shown in
the first figure on this page, it can be seen that the water surface upstream of the bridge did not encounter the
bridge decking for the energy method analysis and therefore was calculated as a low flow profile. However, the
water surface elevation is very close to the bridge decking and, due to the turbulent wave action of the flow, may
jump to pressure flow during the event. Therefore, an additional analysis was performed with the pressure/weir
flow option selected for the high method. This is discussed further in the next section.
For the pressure/weir analysis, the editor allows three input coefficients. The first coefficient is for the flow situation
when only the upstream side of the bridge decking is submerged. For this example, the field was left blank so that
the table values would be used by the program. The second coefficient is for the flow situation when both the
upstream and downstream sides of the bridge decking are submerged. This coefficient was left at a value of 0.8, the
default value. The final field is for the modeler to enter an elevation value for the program to use to determine when
to begin the pressure flow calculations. Pressure flow calculations will begin when the low flow energy gradeline
value is greater than the value entered in this field. If the field is left blank (as for this example), then the program
will use the highest low cord of the bridge decking on the upstream side. Once all of the data had been entered, the
editor was closed by selecting the OK button.
Next, the ineffective flow areas were added at river station 8 (the river station located immediately upstream of the
bridge). A left ineffective flow area was set to begin 10 feet to the left of the bridge opening and a right ineffective
flow area was set to begin 10 feet to the right of the right side of the bridge opening. Both of these ineffective flow
areas were set at an elevation of 1285 feet, the high cord elevation of the bridge.
Then, the geometry file was saved as "Base Geometry Data - Press/Weir." Next, this geometry file and the steady
flow data file were saved as a plan. This was performed by first activating the Steady Flow Analysis Window from
the main program window. The steady flow window is shown in the figure below. Then, the geometry and steady
flow file were selected by depressing the down arrows on the right side of the window.
Next, the Short ID "Press/Weir" was entered in the upper right corner of the editor and the flow regime "Mixed" was
selected. Finally, File and Save Plan As were selected and the title "Putah Creek Bridge - Press/Weir" was entered.
This created a new plan with the pressure/weir geometry file and the steady flow data file. The plan name then
appeared at the top of the steady flow analysis window (as well as on the main program window).
After the new plan was created, the COMPUTE button was depressed to activate the calculations of the water
surface profile. The user can activate this plan by selecting File, Open Plan, and then the plan "Putah Creek Bridge -
Press/Weir" to review the remaining discussion of this example.
The flow through the bridge structure is again Class B flow because the flow passed through critical depth within
the bridge. Finally, the downstream profile is supercritical to the downstream cross section. As before, it was
necessary to perform the analysis in the mixed flow regime in order to obtain a water surface profile in both the
subcritical and supercritical flow regimes.
Expansion Reach Length. To estimate the expansion reach length from Table B.1, the following information was
required:
where:
b=bridge opening width, ft
B=floodplain width, ft
S=bed slope, ft/mi
nob =Manning's n value of the overbank
nc=Manning's n value of the main channel
Lobs=average length of obstruction, ft
From Table B.1, the expansion ratio (ER) was determined to be in the range from 1.2 - 1.5. Using an average value of
1.3 yields an expansion reach length (Le) of:
From the geometric data, the distance from river station 6 to 5.5* was set at 50 feet. This value is approximately
equal to the expansion reach length as determined above and the location of the river stations were not adjusted.
Contraction Reach Length. To estimate the contraction reach length, the contraction ratio from Table B.2 was
determined to be in the range from 0.8 - 1.4. Using an average value of 1.1 yielded a contraction reach length (Lc) of:
From the geometric data, the contraction reach length was set to be 50 feet. For this example, the flow in the main
channel did not exhibit a large degree of contraction or expansion losses. This allowed for shorter expansion and
contraction reach lengths, as determined above. As stated previously, the table values were used as a general
guide, because the data of the example were not within the range of the data used to develop the regression
equations.
The tenth column (Prs O EG) displays the energy gradeline necessary for only pressure flow to be occurring through
the bridge. For this example, this value is 1284.80 feet. The program then compared this pressure only energy
gradeline to the energy gradeline that was used for the low flow analysis (1283.07). The greater of the two values
was then used as the final answer. For this example, the pressure only energy gradeline was the greater value and
the program then used the pressure only method as the final solution.
The ninth column (Prs/Wr EG) shows the calculated energy gradeline value for the situation when both pressure and
weir flow would occur. This situation did not develop for this example. For a further discussion for bridge analyses,
the user is referred to example 2 and to chapter 6 of the User's Manual and chapter 5 of the Hydraulic Reference
Manual.
The bottom of the table displays the errors, warnings, and notes for the river station. At this bridge location, the
notes shown in the figure indicated that the sluice gate equation was used for pressure flow. This equation was
used because the water surface elevation at the river station immediately below the bridge was less than the lowest
value of the low cord for the bridge.
12.8 Summary
As a summary for this example, a profile plot for both the energy and pressure flow analysis is shown in the figure
below. Once the plot was activated, Options and Plans were chosen and both plans were selected to be displayed.
The heading at the top of the figure displays the Short ID's that were used for each plan : 1) Energy and 2) Press/
Weir. The legend displays two lines for the critical depth and two lines for the water surface. Both of the critical
depth lines will coincide since the flow rate was the same for both plans. For the water surface profiles, the lines are
labeled "WS 100yr - Energy" and "WS 100yr - Pressure" for the energy and pressure methods respectively.
(Note: For this example, each plan only had one profile.) Similarly, the legend description for the solid water surface
profile signifies that this line is for the first profile of the energy plan.
A review of this figure shows a significant difference in the calculated water surface profiles for the two plans.
During the energy analysis, the water surface does not encounter the bridge decking and this lead to a low flow
profile. However, the pressure analysis determined that the water surface came in contact with the upstream side
of the bridge. With both of these analyses, the modeler must use engineering judgment to determine which profile
is actually occurring.
Realistically, for this example, as the flow rate begins to increase to the value of 3200 cfs, the flow will most likely be
a low flow profile as calculated by the energy method. At the flow rate of 3200 cfs, the flow may also initially be
occurring as the low flow profile. However, the water surface immediately upstream of the bridge has risen due to
the constriction of the cross sections in the vicinity of the bridge. If the flow is sediment laden and as debris begins
to accumulate in the vicinity of the bridge opening, the water surface may begin to fluctuate due to the turbulent
nature of the flow. These water surface fluctuations may become great enough to cause the water surface to come
in contact with the upstream low cord of the bridge. When this occurs, the flow may "jump" to become pressure
flow through the bridge opening.
To create the river system schematic, the River Reach icon was selected and a line was drawn in the downstream
direction. Then, the program requested the name of the river and reach. The titles "Spring Creek" and "Upper
Reach" were entered respectively. Then the River Reach icon was selected again and the "Lower Reach" was
sketched. The program requested titles for this river reach and then prompted for a title of the stream junction. The
name "Pottsville" was entered for the stream junction title. Finally, the River Reach icon was selected a third time
and the Spruce Creek tributary was drawn and labeled. This created the river system schematic as shown in the
figure above (without the river stations). The cross section data for the river stations were entered next and are
described subsequently.
Finally, a third criterion is concerned with the fact that the program is a one- dimensional model. Therefore, the
cross sections should be located in regions where the flow direction is perpendicular to the section. For example, as
shown in the figure above, the cross section number 6 should be placed adequately downstream from the junction
so that the flow is predominately in the downstream direction. This cross section should not be placed in the
junction where the turbulent mixing of the flow is occurring. However, as discussed previously, the cross sections
should be placed close to the stream junction to accurately evaluate the energy losses across the junction and to
provide a reasonable result for the momentum method, which equates the two upstream water surface elevations.
Observations of water surface elevations at high flows can assist the modeler in determining the appropriate
locations for the cross sections.
The left side of the editor is where the station and elevation data for the cross section coordinates were entered.
The right side of the editor displays the entered values for the downstream reach lengths, Manning's n values, main
channel bank stations, and the contraction and expansion coefficients. The downstream reach lengths for the last
cross section of each river reach should be set at zero or left blank. Therefore, for this example, the downstream
reach lengths at the river stations of 10.106, 10.000, and 0.013 for Upper Reach, Lower Reach, and Spruce Creek,
respectively, were set at 0.
Finally, any levees, ineffective flow areas, blocked obstructions, etc. would be entered at this time. For this example,
the flow and cross sectional data did not provide for the use of any of these options.
To enter the stream junction data, first one of the stream junctions was selected by depressing the down arrow
adjacent to the Junction Name box. For this example, there was only one junction, which had been named
"Pottsville" and this junction name automatically appeared. Next, a description for the junction was entered. The
entire description can be viewed by selecting the "..." button. The "..." button can then be re-selected to exit the
description display.
The next item of information required is the stream lengths across the junction. The table at the bottom left side of
the editor will automatically display the names of the river reaches at the junction for which the user must enter the
reach lengths. For this example, a reach length of 80 feet was entered as the distance across the junction from
Upper Reach to Lower Reach. The program will then use this distance of 80 feet as the length from river station
10.106 (the downstream river station of Upper Reach) to river station 10.091 (the upstream river station of Lower
Reach). Similarly, a distance of 70 feet was entered as the length across the junction from the upstream river station
on Lower Reach (10.091) to the downstream river station on Spruce Creek (0.013).
The last item for the Junction Data Editor is the Steady Flow Computation Mode. The user must select either the
energy or the momentum method for the steady flow computational procedure. For this example, the energy
method was selected. After the discussion of the output for the energy method analysis, the momentum method
will be selected and the results from the two methods will be compared. A discussion on the computational
procedures for each method will be addressed during the review of the output for each method.
The selection of the energy method for the junction analysis completed the geometric input for this example. The
junction editor was then closed by selecting the OK button. Finally, the geometric data was then saved by selecting
File and then Save Geometric Data As from the Geometric Data Editor. The title "Base Geometry - Energy
Junction" was entered for the name of the file. The next procedure was to enter the steady flow data.
For this example, a subcritical analysis was performed and therefore a downstream boundary condition was
required to be entered for the Lower Reach of Spring Creek. To enter this boundary condition, the downstream field
for Lower Reach was selected and then the boundary condition Normal Depth was chosen. A slope of 0.001 was
then entered for this boundary condition. For a further discussion on boundary conditions, the user is referred to
chapter 7 of the User's Manual and to chapter 3 of the Hydraulic Reference Manual.
After all of the steady flow data had been entered, the steady flow data file was then saved by selecting File and
then Save Steady Flow Data As from the Steady Flow Data Editor. The title "10 Year Profile" was then entered and
the OK button selected.
At the top of the Steady Flow Analysis Window, a Short ID was entered as "Energy." The next step was to select
the appropriate flow regime for the analysis. For this example, the Subcritical regime was selected. Then, the
geometry file "Base Geometry - Energy Junction" and the steady flow file "10 Year Profile" were selected by
depressing the down arrows on the right side of the window. (Note: At this point in the example, there was only one
geometry and one flow file. Therefore, this step was not necessary.) To save these files as a plan, select File and
then Save Plan As. The title "Junction - Energy" was then entered as the plan title and the OK button was selected.
This associated the geometry and the steady flow file as a plan and the name of the plan then appeared on the
Steady Flow Analysis Window (as well as on the main program window). Finally, the COMPUTE button was
selected to perform the analysis.
At this time, the user can select which reaches to be displayed by selecting Options and then Reaches. In this
manner, the water surface profile for the combination of Spruce Creek and Lower Reach can be viewed. The user
can also select to have all three reaches displayed simultaneously.
To determine the water surface profile across the stream junction, the program used the energy based solution
routine for this plan. Since this was a subcritical flow analysis, the program started the calculations at the
downstream end of the Lower Reach and computed the water surface profile up to the upstream river station of
Lower Reach (10.091). Then the program performed standard step calculations separately across the stream
junction to each of the downstream river stations of Upper Reach and Spruce Creek. During each of the separate
calculations, the friction losses and the contraction and expansion losses were calculated to balance the energy
from river station 10.091 to each of the river stations 10.106 and 0.013 of Upper Reach and Spruce Creek,
respectively. By performing the calculations in this manner, the downstream water surface elevations of Upper
Reach and Spruce Creek do not necessarily have to coincide.
Similarly, the user can toggle the table to display the information for Spruce Creek and follow the same procedure
to determine the energy losses across the junction from the upstream river station of Lower Reach to the
downstream river station of Spruce Creek. (Note: If the table does not display the information for all of the river
reaches, the user can select the appropriate reaches for display under the Options menu.)
The components of Junction Name, Description, and Reach Lengths were kept the same as described for the
energy method analysis. However, the Steady Flow Computation Mode was selected as Momentum for the
subsequent analysis. For the momentum method, the user must enter the angles at which the reaches are entering
or leaving the junction. For a flow combining situation (such as this example), the user must enter the angles of the
inflow reaches as measured from a line perpendicular to the upstream cross section of the outflow reach. (Note:
Chapter 4 in the Hydraulic Reference Manual will assist the user to visualize the flow angles.)
For this example, the tributary angle from Upper Reach to Lower Reach was entered as 0 degrees. This implies that
the flow will continue in a straight line. The flow angle from Spruce Creek to Lower Reach was then entered as 45
degrees, as measured from the survey data.
Finally, the user must select whether the program should include the weight and/or the friction terms in the
momentum equation. For this example, the terms were included by selecting the appropriate boxes in the lower
right corner of the window. The user should refer to chapter 4 of the Hydraulic Reference Manual for a further
discussion of the momentum equations.
After these changes were made, the Apply Data button was selected and the junction editor was closed. Since this
geometry file was changed, it was then saved as a new file. This was performed by selecting File and then Save
Geometry As. The title "Base Geometry - Momentum Junction" was entered and the OK button selected. Then the
Geometric Data Editor was closed.
From the main program window, Run and then Steady Flow Analysis were then selected. This activated the
Steady Flow Analysis Window as shown in the following figure.
In the upper right corner of the window, a Short ID was entered as "Momentum." Then, the geometry file "Base
Geometry - Momentum Junction" and the flow file "10 Year Profile" were selected by depressing the down arrows
on the right side of the window. Next, the Flow Regime was selected as Subcritical. Then, File and Save Plan As
were selected and the title "Junction - Momentum" was entered. This plan title then appeared at the top of the
window (as well as on the main program window). Finally, the COMPUTE button was selected to perform the
analysis.
To determine the water surface profiles, the program started at the downstream end of the Lower Reach (since this
was a subcritical flow analysis) and used the standard step procedure up to the upstream end of Lower Reach. At
the stream junction, the program then used the momentum method to balance the forces across the junction. The
momentum method (for this example) will solve simultaneously the forces in the X-direction for the flow at river
stations 10.091 (the upstream river station of Lower River), 10.106 (the downstream river station of Upper Reach),
and 0.013 (the downstream river station of Spruce Creek). The X-direction is determined as the direction of the flow
out of the junction (the direction of flow at river station 10.091).
Since the water surface elevations at the downstream ends of Upper Reach and Spruce Creek are two unknown
values and are solved simultaneously, the program will assume that these water surface elevations are equal to
each other. Therefore, it is necessary that the cross sections be placed close to the stream junction in order to
minimize the error associated with this assumption. For a more detailed discussion on the momentum method
analysis, the user is referred to chapter 4 of the Hydraulic Reference Manual.
previously, the first two columns in the table are in fixed format and display the river reach and the river stationing.
(Note: For this table, all 3 river reaches were selected to be displayed.) The columns in the non-fixed area of the
table display the energy gradeline, the water surface elevation, the velocity head, the friction losses, and the
contraction/expansion losses. The remaining portion of the table can be viewed by depressing the left and right
arrows at the bottom of the window.
As shown in the figure below, the water surface elevations for Lower Reach are exactly the same values in the
previous run because nothing has changed in the downstream reach. However, for this plan, the water surface
elevations across the stream junction were determined using the momentum method. As discussed previously, to
perform the momentum calculations at the stream junction, the program will equate the two upstream water
surface elevations. This is seen to have occurred as the water surface elevations at river station 10.106 (the
downstream station of Upper Reach) and river station 0.013 (the downstream river station of Spruce Creek) are
both equal to 75.50 feet. Then, the program determined separately the remaining water surface profiles for Upper
Reach and Spruce Creek by using the standard step procedure in the upstream direction.
The figure below shows the water surface elevations for Upper Reach and Lower Reach for both methods of
junction analysis. At the top of the figure, the heading lists the project name "Stream Junction - Example 10" and
By visually comparing the two results, it can be seen that the calculated water surface profile for Lower Reach was
the same for both plans. However, the water surface profile for Upper Reach was different. This is due to the fact
that the starting downstream water surface elevation at river station 10.106 of Upper Reach was different, for the
two methods of computation across the junction. If the standard step procedure was allowed to continue for a
longer distance upstream on Upper Reach, both of the water surface profiles would eventually converge.
At this time, the user can select to view the combination of Spruce Creek and Lower Reach and observe a similar
result. The profile would show that the momentum method produced a higher resulting water surface profile on
Spruce Creek than the energy method. However, the difference of the water surface profiles for Spruce Creek is not
as significant as that for Upper Reach.
The energy method uses friction and the coefficients of contraction and expansion in determining the energy losses
across the junction. The user can adjust these coefficients to account for any abrupt transitions that occur in the
cross sectional area. Additionally, the user could adjust the Manning's n value at the stream junction to account for
additional internal energy losses associated with the junction. By adjusting these parameters, the calculated water
surface profile can be calibrated to actual measured water surface elevations.
For the momentum method, the program determines the water surface elevations across the stream junction by
taking into account the forces associated with the flow. To balance the forces across the junction, the program only
uses the forces in the X-direction. This direction is determined as being perpendicular to the outflow direction (for a
flow combining situation). Therefore, the tributary angles, as entered by the user, are crucial for the accurate
calculation of the forces. Additionally, as for the energy method, the Manning's n values can be adjusted to account
for additional friction losses associated with the stream junction.
13.8 Summary
Both the energy and the momentum methods were used to determine the water surface profiles across the stream
junction. For this example, the resulting water surface profiles differed upstream of the junction for the two
calculation procedures.
The energy method used the standard step procedure with the coefficients of expansion and contraction and the
Manning's n value to account for energy losses across the junction. With the energy method, the program
calculated separately the resulting upstream water surface elevations. The momentum method equated the forces
in one dimension and used the tributary angles to balance the forces at the stream junction. Additionally, the
Manning's n values are used to account for friction losses. The momentum method should be used when the
tributary flow angles play an important role in influencing the water surface around the junction. However, in order
to solve the momentum calculations, the upstream water surface elevations were equated. To reduce the error
associated with this assumption, the cross sections should be located close to the stream junction.
The momentum method is an attempt at a more theoretical analysis of the stream junction. However, the user
should be aware of the limitations of this one-dimensional analysis. To determine the most applicable method for
the analysis, the user should compare the results to observed data and calibrate the model as deemed appropriate.
Along this reach of Pine Creek, a bridge was entered at river station 10.36. The bridge data were entered by
selecting the Brdg/Culv icon on the left side of the Geometric Data Editor. This activated the Bridge/Culvert Data
Editor as shown in the figure below. Then, the bridge information for the deck/roadway, piers, sloping abutments,
and bridge modeling approach were entered by selecting the appropriate icons on the left side of the Bridge/
Culvert Data Editor. The bridge opening between the sloping abutments is approximately 600 feet wide and the
bridge is supported by six piers, each with a width of 5 feet. The high and low cord values for the bridge deck are 22
and 18 feet, respectively. The user can select the appropriate icons to review the bridge data.
After all of the geometric data had been entered, the editors were closed and the geometry was saved as "Base
Geometry." Next, the steady flow data were entered.
To perform the bridge scour calculations, the program requires detailed values of the depth and velocity within the
cross sections located just upstream from the bridge (cross section 10.37 for this example) and at the approach
section (cross section 10.48). Therefore, the modeler is required to set the flow distribution option for these two
cross sections. For this example, the flow distributions were selected for the entire river reach. As shown in the
figure above, the left and right overbanks were divided into 5 subsections each, and the main channel was divided
into 20 subsections. This will allow the program to produce detailed results of the distribution of depth and velocity
at the cross sections.
The number of subsections is dependent upon such factors as the cross section geometry, the bridge opening
width, and the number of piers. The modeler should perform the hydraulic calculations with different numbers of
subsections to evaluate the impact on the bridge scour results. It is recommended to use fewer subsections,
however, an adequate number of subsections is required to determine the hydraulic properties. For this example,
the bridge scour calculations were also performed using 20 subsections for the main channel, and no appreciable
changes were observed in the scour results. For a further discussion on the flow distribution option, the modeler is
referred to chapter 7 of the User's Manual and to chapter 4 of the Hydraulic Reference Manual.
Finally, the flow distribution editor was closed and the data were saved as a plan entitled "Scour Plan 1." The
COMPUTE button was then selected to execute the analysis.
At this point, the modeler should review the output from the hydraulic analysis and calibrate the model. It is
important to obtain a good working model of the river system before attempting to perform a bridge scour analysis.
For this example, the hydraulic analysis included the evaluation of the expansion and contraction reach lengths
according to the procedures as outlined in the Hydraulic Reference Manual. Finally, after a working model has
been developed, the user should evaluate the long-term aggradation or degradation for the river reach and
incorporate this analysis into the working model.
The remaining portion of the editor is divided into three areas: input data tabs, a graphic, and a results window.
There are three tabs, one for each of the three types of scour computations: contraction, pier, and abutment. The
graphic displays the bridge cross section (inside upstream). When the Compute button is selected, the scour results
will be displayed graphically on the cross section and in tabular format in the results window. The following
sections describe the parameters for each of the three data tabs.
Finally, the down arrows adjacent to Equation were selected and the default option was chosen. This informed the
program to use either the clear water or the live bed scour equation as determined from equation 10-1 in the
Hydraulic Reference Manual.
To perform the contraction scour computations, the Compute button at the top of the editor was selected. When
the calculations were completed, the results appeared in tabular form in the lower right corner of the editor and in
graphical form on the bridge cross-section plot, as shown in the first figure above.
As a review of the results for the contraction scour, the critical velocity (Vc) for the LOB was determined to be 2.63
ft/s, from equation 10-1. This value is greater than the velocity at the approach section (V1 = 2.00) in the LOB;
therefore the clear water scour equation was used for the LOB, as listed in the summary table. Comparatively, the
live-bed scour equation was used for the main channel because the critical velocity (Vc = 2.99) was less than the
approach section velocity (V1 = 4.43), in the main channel. Finally, the contraction scour depth (Ys) was determined
to be 2.07 and 6.67 feet for the LOB and main channel, respectively. As a final note, there was no contraction scour
in the ROB because the right abutment extended into the main channel. These contraction scour depths are also
shown on the graphic display of the bridge cross section in the first figure.
For this example, the Maximum V1 Y1 option was selected to inform the program to use the maximum value of the
depth and velocity values, as opposed to the values upstream from each pier. Then, Method was selected as the
"CSU equation." Next, the Pier # option was selected as "Apply to All Piers" to inform the program that the data will
be used for all of the piers. (The user has the option of entering the data for each individual pier.)
Next, the Shape of the piers was selected as "Round nose" which set the K1 value to be 1.00. Then, the D50 was
entered as 2.01 mm. The angle was set to be 0 degrees which set the K2 value to be 1.00. Next, the bed condition
was selected as "Clear-Water Scour" (this set K3 = 1.1) and the D95 was entered as 2.44 mm.
This completed the required user input and then the Compute button was selected. The results were then
displayed graphically and in the summary table and showed that the pier scour depth (Ys) was 10.61 feet, as shown
in the figure above. (Note: When the compute button was selected, the program automatically computed all 3 scour
depths: contraction, pier, and abutment.)
This completed the data entry for the abutment scour and the Compute button was selected. The results were then
displayed on the graphic and in the summary table, as shown in the figure below. The results show that the HIRE
equation was used for both the left and right abutment and the magnitude of the scour was 10.92 and 15.49,
respectively. Additionally, the summary table displayed the values of the Froude numbers used for the calculation.
Contraction Scour
Pier Scour
Abutment Scour
Channel: 17.28
The first three portions of the table display the results of the contraction, pier, and abutment sour, as discussed
previously. The final portion of the table displays the combined scour depths. For this example, the pier and
contraction scour was 12.67 feet (= 10.61 + 2.06) for the left bank and 17.28 feet (= 10.61 + 6.67) for the main
channel. Additionally, the total left abutment and contraction scour was 12.98 feet (= 10.92 + 2.06) and the right
abutment and contraction scour was 21.55 feet (= 14.90 + 6.67). The contraction scour for the right abutment was
the contraction scour for the main channel because the right abutment extended into the main channel.
Finally, the total scour is displayed graphically, as shown in the figure below. (Note: The graphic has been zoomed
in to see more detail.) As shown in the legend, the long dashed line represents the contraction scour and the short
dashed line portrays the total scour. This graphic was obtained by selecting Copy Plot to Clipboard from the File
menu.
14.5 Summary
To perform the bridge scour computations, the user must first develop a model of the river system to determine the
hydraulic parameters. Then, the program will automatically incorporate the hydraulic results into the bridge scour
editor. The user can adjust any of the values that the program has selected. For each particular scour computation,
the modeler is required to enter only a minimal amount of additional data. The results for the scour analysis are
then presented in tabular form and graphically. Finally, the user can select Detailed Report from the Bridge Scour
Editor to obtain a table displaying a full listing of all the input data used and the results of the analysis.
To enter the data for the weir, the Weir/Embankment icon was selected from the left side of the Inline Structure
Data Editor. This activated the Inline Structure Station Elevation Data Editor as shown in the figure below. This
editor is similar to the deck/roadway editor used for bridges and culverts.
The top row of the editor consists of three data entry fields. For the first field, the user must enter the Distance from
the upstream cross section (41.76) to the upstream side of the weir. For this example, this distance was 20 feet.
Next, the Width of the weir was entered as 50 feet. This is a total distance of 70 feet. From the figure titled, "Reach
Length Table for Nittany River', it can be seen that the distance from river station 41.76 to 41.74 is 90 feet. Therefore,
the distance from the downstream end of the weir to cross-section 41.74 is 20 feet. The last field in the top row of
the editor is the Weir Coefficient. This will be discussed shortly.
The central portion of the editor consists of a table in which the user must enter the station and elevation data for
the weir. For this example, the emergency spillway is located on the left at an elevation of 9.5 feet from Station 61 to
190; however, the entire top-of-dam is defined. The first station of the weir was entered as 0 and the last station is at
1000 feet. With these weir stations and elevation, the program will block out the entire area below the weir crest. In
this manner, the data entry is similar as that for a culvert.
Additionally, it should be noted that the weir station values of 0 and 1000 occur beyond the limits of the cross
section data. As for the bridge and culvert routines, the program will automatically "clip off" the excess area so that
the weir coincides with the cross section geometry.
The next fields are for the upstream and downstream embankment side slopes. For this example, the slope of two
was entered for the US and DS Embankment S.S. fields.
At the bottom of the editor are several other required variables. The Min Weir Flow Elevation was left blank which
implies that the lowest elevation of the weir will be used to determine when weir flow begins to occur. Finally, the
shape of the weir was entered as "Ogee" for the Submergence criteria. When ogee was selected, the editor
expanded to allow for two more fields of entry. These fields are the Spillway Approach Height and Design Energy
Head. The approach height was entered as 24 feet and the design head was 3 feet for the ogee shape. To determine
the weir coefficient with these design parameters, the Cd button was selected and the program calculated a
coefficient of 3.95, as shown in the figure below.
The "Yes" button was selected and then the coefficient appeared at the top of the Inline Structure Station
Elevation Data Editor in the Weir Coefficient field. If the weir shape had been selected as "Broad Crested", the user
is required to enter the value of the weir coefficient. This completed the data entry for the weir. Next, the data for
the gates were entered.
When the Gate Editor was activated, the first gate was labeled as "Gate #1." For this example, the Rename button
was selected and the label "Left Group" was entered. Next, the Height, Width, and Invert for the gates of the Left
Group were entered as 10, 30, and 0 feet, respectively. On the lower left side of the editor, the centerline stations for
the five gates in the "Left Group" were entered as shown in the figure above. As these values were entered, the
counter field # Openings increased to represent the total number of gates for the group (5 for this example).
The remaining portion of the editor is divided into two sections, one for the Gate Flow data and one for Weir Flow
data. The gate data is used when the water surface upstream of the gate is greater than 1.25 times the gate opening
(as measured from the gate invert). At this water surface elevation, the gate is controlling the flow rate. The weir
data is the shape of the weir under the gate and is used when the upstream water surface is less than or equal to the
gate opening. At this water surface elevation, the weir under the gate is controlling the flow through the gate
opening (i.e., the water is not in contact with the gate).
For the gate data, the first field is the Gate Type. By selecting the down arrow, the type "Radial" was chosen. When
the gate type was selected, the Radial Discharge Coefficient, Trunnion Exponent, Opening Exponent, Head Exponent,
and Trunnion Height values were automatically set to 0.80, 0.16, 0.72, 0.62, and 10.0 respectively. The orifice
coefficient of 0.8 was entered for fully submerged flow conditions. Additionally, the head reference for the gate
equation defaults to the "Sill (invert)" of the gate opening. This is appropriate for gates that release water to a
channel, rather than to the free atmosphere.
For the Weir Flow data, the Weir Shape was selected as "Ogee". This caused the editor to add the data fields for
Spillway Approach Height and Design Energy Head. The distances of 14 and 3 feet were then entered for each of
these fields, respectively. Finally, the Cd button was selected and a window appeared similar to the figure, "Ogee
Weir Shape Coefficient", with a coefficient of 3.91. The "Yes" button was selected and the weir coefficient appeared
in the weir data area at the bottom of the Gate Editor.
This completed the data entry for the gates in the "Left Group." Next, the Add button at the top of the Gate Editor
was selected and this added another gate group. The group was renamed to "Center Group" and the data for 5 new
gates were entered exactly as for the Left Group, except for the centerline stations. Finally, a third gate group was
added and renamed "Right Group," with the data entry as for the two previous groups with new centerline stations.
The OK button was selected at the bottom of the Gate Editor and the gates appeared on the Inline Structure Data
Editor as shown in the figure "Inline Structure Data Editor" above. (Note: The ineffective flow areas will be entered
subsequently.) At this point, the user should zoom in on the gate openings to ensure that they do not overlap and
appear as intended. The Inline Structure Data Editor was then closed.
where:
= gate openings area, ft
= floodplain width, ft
= channel slope, ft/mi
= Manning's n value of the overbank at river station 41.78
= Manning's n value of the main channel at river station 41.78
From Table B.1, the expansion ratio (ER) is approximately 2.0. Using an average length of obstruction (Lobs) of
approximately 100 feet yields an expansion reach length (Le) of:
The expansion reach length is measured downstream from river station 41.74. Therefore, a cross section (41.70) was
placed 200 feet downstream from river station 41.74. River station 41.70 represents the cross section where the flow
is fully expanded.
For the contraction reach length, the contraction ratio (CR) was obtained as 1.0 from Table B.2. With this value, the
contraction reach length (Lc) is:
The contraction reach length is measured upstream from river station 41.76. Therefore, a cross section (41.78) was
placed 100 feet upstream from river station 41.76. River station 41.78 represents the cross section where the flow
lines are parallel.
As a final note, the values for b, B, and Lobs were approximated for the flow rate of 75,000 cfs. As the flow rate
changes, the length of expansion and contraction would also change. For this example, the values that were
determined for this flow rate were held constant for all of the flow rates.
This concluded the entry for all of the geometric data. At this point, the geometric data was saved as the file "Gate
Geometry with 3 Gate Groups." Next, the flow data was entered.
At the top portion of the editor, the River "Nittany River," Reach "Weir Reach" and the River Station "41.75" were
selected. The Description is the same as was entered in the Inline Structure Data Editor (see third figure in the
Geometric Data Section). The # Gate Groups field shows that there are 3 gate groups at this river station. The table in
the central portion of the editor has 3 rows for this example, one row for each of the gate groups. The first column
lists the descriptions for the gate groups, as they were named in the Gate Editor (see the Geometric Data Section).
The second column displays the number of gate openings for each gate group (5 for each gate group for this
example). The third column displays the maximum gate height for each gate group (10 feet for each gate group for
this example).
The remaining portion of the editor consists of entry fields for the number of gates opened and the opening heights
of the gates for each flow profile. For this example, for profile 2, the Left and Right Gate Groups were set to have 2
gates open, each with an opening height of 3 feet. The Center Group was set to have 5 gates open, each at 5 feet.
Since the gates were entered as three groups, the flexibility existed to have each gate group opened at different
heights. If all 15 gates had been entered as only one group, then all of the gates that were open would have to have
been set at the same opening height.
The user can toggle across the table to view the number of gates open and the gate opening heights for all of the
profiles. During the analysis of the output, the various gate settings will be discussed. This concluded the data entry
for this example. At this point, the OK button at the bottom of the editor was selected and the flow data was saved
as "7 Flow Profiles."
First, the Short ID was entered as "3 groups." The geometry file was then selected as "Gate Geometry with 3 Gate
Groups" and the flow file was "8 Flow Profiles." Next, the Flow Regime was selected as "Subcritical." Then, File and
Save Plan As were chosen and the information was saved as the plan "Gated Spillway." This plan name then
appeared on the Steady Flow Analysis Window, as well as on the main program window. Finally, the COMPUTE
button was selected to perform the analysis.
The Center Group had been set (as shown in Steady Flow Data section) to have 5 gates open at a height of 5 feet.
The figure below displays this information in the right column. With a gate invert of -10 and a gate opening height of
5, the top of the gate opening was at -5 feet. The left side of the figure shows that the water surface at river station
41.76 was at an elevation of 4.56 feet. Therefore, the water surface did not come into contact with the top of the
gate opening and weir flow through the gate openings occurred. The weir data that was use to calculate the
upstream energy grade line was the data that had been entered in the Gate Editor (see the Geometry Data section),
not the weir data as entered in the Inline Structure Station Elevation Data Editor (see the Geometry Data
section). Additionally, the Gate Area field shows a gate flow area of 136.81 ft2. With a gate opening height of 5 feet
and a width of 30 feet, the total gate opening area is 150 ft2. This shows that the gate area was not flowing full.
To review the data for the second profile, the down arrow for the Profile field was depressed and "PF#2" was
selected. The Inline Structure Type Detailed Output Table for profile 2 is shown in the figure below. For this
profile for the Center Group, 5 gates were set to an opening height of 5 feet, as before (as shown in the right column
of the table).
With an invert at 0 and an opening height of 5, the top of the gate openings were at an elevation of 5.00 feet. The
upstream water surface elevation is shown to be at 4.05 feet, which is above the top of the gate openings.
Additionally, the Gate Area field shows a flow area of 121.45 ft2 (for each gate opening). This equals the gate
opening area (4.05 ft high x 30 feet wide) in which water was flowing. The depth of flow (as measured from the
invert of the gate) was 4.05 feet, which is less than the gate opening height. Therefore, the flow was still being
computed as weir flow.
Additionally, for the second profile, the Left and Right Gate Groups were each set to have 2 gates open, at a height
of 3 feet. By depressing the arrow for the Gate ID field, the Left (or Right) group can be selected. For this group, gate-
controlled flow occurred through the gate openings because the top of the gate openings were at an elevation of 3
feet. The field Gate Q Total shows that 2926.11 cfs was the total flow through the gate openings. Since the Left and
Right groups had the same gate settings, the total flow in the cross section was 2 times the flow in the Left Group
plus the flow in the Center Group = 2(2925.04) + 4149.93 = 10000 cfs, the total flow rate for the second profile.
Similarly, for the remaining profiles, the flows through the gate openings can be reviewed by pressing the down
arrow next to the Profile: box. The Weir Type Detailed Output Table for profile 7 is shown in the figure below. For
this profile, the weir flow data is shown on the left side of the table, based on the data as entered in the Inline
Structure Station Elevation Data Editor (see 2022-02-08_13-40-16_Geometric Data (Ex12)(see page 172)). The total
weir flow is shown to be 36229.06 cfs. Because the Left, Center, and Right Gate Groups were entered with the same
gate settings, the total flow in the cross section is 3(12923.65) + 36229.06 = 75,000 cfs, the total flow for the profile.
The table also displays other weir data such as the left and right station, average depth, and submergence.
15.5 Summary
This example computed 7 flow profiles for the reach of Nittany River, which included an inline weir and spillway.
The gates for the inline structure were divided into 3 groups, with 5 gates in each group. This provided for flexibility
when setting the number of gates opened and the gate opening heights for each profile because the opening
heights must be the same for all of the gates opened in each gate group.
By reviewing the water surface profiles and the inline structure tables, the user can determine the type of flow
through the gate openings and determine if adjustments to the gate settings are required to provide for a selected
water surface elevation. The Inline Structure Output table provides detailed output for each gate group, for any
profile. The Profile Summary Output Table - Inline Structure provides upstream energy and water surface
elevations along with the total weir and gate flow. For the maximum discharge profile, the entire weir profile is
overflowing, a condition that may not be structurally sound.
The editing commands of Add Points to a Reach and Move Object were used to curve the river system schematic.
This was performed for aesthetic purposes only, and does not effect the hydraulic computations.
Manning's n values, and main channel bank stations were obtained from the field data displayed on the USGS atlas.
A summary table of the reach lengths, Manning's n values, or contraction and expansion coefficients for the cross
sections can be seen by selecting the relevant listing under the Tables menu on the Geometric Data Editor. The
summary table of the reach lengths for the cross sections is shown in the bottom, figure below.
First Cross Section. Ideally, the first cross section should be located sufficiently downstream from the bridge so that
the flow is not affected by the structure (i.e., the flow has fully expanded). This distance (the expansion reach
length) should generally be determined by field investigation during high flows and will vary depending on the
degree of constriction, the shape of the constriction, the magnitude of the flow, and the velocity of the flow. In
order to provide better guidance to determine the location of the fully expanded cross section, a study was
performed by the Hydrologic Engineering Center \[HEC-1995\]. This study focused on determining the expansion
reach length, the contraction reach length, and the expansion and contraction energy loss coefficients. The results
of the study are summarized in Appendix B of the Hydraulic Reference Manual. To determine an initial estimate of
the expansion reach length, the following information was required:
where:
= bridge opening width
= total floodplain width
= slope
= Manning's n value for the overbank
Manning's n value for the main channel
= average length of the side obstructions
With the above information, an initial estimate of the expansion reach length (Le) was obtained by using the table
values in Appendix B of the Hydraulic Reference Manual. From Table B.1, an average value of the expansion ratio
(ER) was 1.2. Using this value, the expansion reach length is:
From the reach lengths table (above), the distance from cross section 52.36 to 52.00 is 1956 feet (= 426 + 1530).
Therefore, cross section 52.00 was initially considered as the cross section of fully expanded flow. Since cross
section 52.29 is in the zone of flow expansion, ineffective flow areas were placed in the overbank areas. This will be
discussed further during the determination of the ineffective flow areas. After the steady flow analysis was
performed, the expansion reach length was evaluated using regression equations that are presented in the HEC
study \[HEC-1995\]. The equations require flow parameters at the initially chosen location and cannot be evaluated
until after the steady flow analysis is performed. These procedures will be described near the end of this example.
Second Cross Section. The second cross section used by the program to analyze the energy losses through the
bridge is located a short distance downstream of the structure. This section should be very close to the bridge and
reflect the effective flow area on the downstream side of the bridge. For this example, a roadway embankment
sloped gradually from the roadway decking on both sides of the roadway. Cross section 52.36 was located at the toe
of the roadway embankment and was used to represent the effective flow area on the downstream side of the
bridge opening. The program will superimpose the bridge geometry onto this cross section to develop a cross
section inside the bridge at the downstream end.
Third Cross Section. The third cross section is located a short distance upstream from the bridge and should reflect
the length required for the abrupt acceleration and contraction of the flow that occurs in the immediate area of the
opening. As for the previous cross section, this cross section should also exhibit the effective flow areas on the
upstream side of the bridge. For this example, cross section 52.38 was located at the toe of the roadway
embankment on the upstream side of the bridge. Similar to the previous cross section, the program will
superimpose the bridge geometry onto this cross section to develop a cross section inside the bridge at the
upstream end.
Fourth Cross Section. The fourth cross section is located upstream from the bridge where the flow lines are parallel
and the cross section exhibits fully effective flow. To determine an initial estimate of the contraction reach length
, a value of the contraction ratio (CR) was obtained as 0.8 from Table B.2, Appendix B of the Hydraulic
Reference Manual. This yields a contraction reach length of:
From the reach lengths table (see recent figure above), the distance from cross section 52.46 to 52.38 is 380 feet. For
this example, cross section 52.46 was initially used as the section where the flow lines were parallel. After the steady
flow analysis, the location of this cross section was evaluated using the procedures as outlined in the recent HEC
study [HEC-1995]. This evaluation will be presented in the discussion near the end of this example.
Additionally, the HEC study [HEC-1995] provided guidance to determine the expansion and contraction coefficients.
The contraction and expansion coefficients are used by the program to determine the transition energy losses
between two adjacent cross sections. From the data provided by the recent HEC study \[HEC-1995\], gradual
transition contraction and expansion coefficients are 0.1 and 0.3, and typical bridge contraction and expansion
coefficients are 0.3 and 0.5, respectively. For situations near bridges where abrupt changes are occurring, the
coefficients may take larger values of 0.5 and 0.8 for contractions and expansions, respectively. A listing of the
selected values for the river reach can be viewed by selecting Tables and then Coefficients from the Geometric
Data Editor. This table is shown in the figure below and displays the values selected for the river cross sections.
Typical gradual transition values were selected for stations away from the bridge. However, near the bridge section,
the coefficients were increased to 0.3 and 0.5 to represent greater energy losses. For additional discussion
concerning contraction and expansion coefficients at bridges, refer to Chapter 5 of the Hydraulic Reference
Manual. After the steady flow analysis, the selection of the coefficients were evaluated and will be discussed near
the end of this example.
This completed the input for the cross section geometric data. Next, the bridge geometry data was entered as
outlined in the proceeding section.
for the bridge deck were entered. The Deck/Roadway icon was selected and the Deck/Roadway Data Editor
appeared as shown in figure below.
The first input at the top of the editor is the distance from the upstream side of the bridge deck to the cross section
immediately upstream from the bridge (cross section 52.38). This distance was determined to be 26 feet from the
USGS atlas. In the next field, the bridge deck width of 31 feet was entered. This is a total distance of 31 + 26 = 57 feet.
From the reach lengths table (see Cross Section Placement figure), the distance from cross section 52.38 to 52.36 is
84 feet. This leaves 84 - 57 = 27 feet from the downstream side of the bridge deck to cross-section 52.36. Finally, at
the top of the Deck/Roadway Data Editor, a weir flow coefficient of 2.6 was selected for the analysis.
The central section of the Deck/Roadway Editor is comprised of columns for input of the station, high cord
elevation, and low cord elevation for both the upstream and downstream sides of the bridge deck. The data were
entered from left to right in cross section stationing and the area between the high and low cord comprised the
bridge structure. The stationing of the upstream side of the deck was based on the stationing of the cross section
located immediately upstream. Likewise, the stationing of the downstream side of the deck was based on the
stationing of the cross section placed immediately downstream. If both the upstream and downstream data are
identical, the user needs only to input the upstream data and then select Copy Up to Down to enter the
downstream data.
As a final note, the low cord elevations that are concurrent with the ground elevation were entered as a value lower
than the ground elevation. The program will automatically clip off and remove the deck/roadway area below the
ground. For example, at station 0, a low cord elevation of 0 feet was entered. However, the actual ground elevation
at this point is approximately 340 feet. Therefore, the program will automatically remove the area of the roadway
below the ground. Additionally, the last station was entered as a value of 5000 feet. This stationing ensured that the
roadway and decking extended into the limits of the cross section geometry. As described previously, the program
will clip off the area beyond the limits of the cross section geometry.
Below the deck station and elevation data, there are two fields to enter the upstream and downstream
embankment side slopes. These values are used by the WSPRO bridge analysis method and for graphical
representation on the profile plot. For this example, values of 2.0 (horizontal to 1 vertical) were entered for both the
US Embankment SS and DS Embankment SS, as measured from the USGS Atlas.
At the bottom of the Deck\Roadway Data Editor, there are three additional fields for data entry. The first is the Max
Allowable Submergence. This input is a ratio of downstream water surface to upstream energy, as measured
above the minimum weir elevation. When the ratio is exceeded, the program will no longer consider the bridge deck
to act as a weir and will switch the computation mode to the energy (standard step) method. For this example, the
default value of 0.95 (95 %) was selected.
The second field at the bottom of the editor is the Min Weir Flow Elevation. This is the elevation that determines
when weir flow will start to occur over the bridge. If this field is left blank (as for this example), the program will
default to use the lowest high cord value on the upstream side of the bridge. Finally, the last field at the bottom of
the editor is the selection of the weir flow submergence method. This method will determine the reduction of the
weir flow coefficient due to submergence. For this example, a broad crested weir flow submergence method was
selected. Upon entering all of the above data, the OK button was selected to exit the Deck/Roadway Data Editor.
Bridge Pier Geometry. To enter the pier data, the Pier icon was selected from the Bridge/Culvert Data Editor.
This resulted in the display shown in the figure below. The modeler should not include the piers as part of the
ground or bridge deck/roadway because pier-loss equations use the separate bridge pier data during the
computations.
The program will establish the first pier as pier number 1. As shown in the figure above, the upstream and
downstream stations were entered for the centerline of the first pier. The upstream and downstream stations were
based on the geometry of the cross sections located immediately upstream (cross section 52.38) and immediately
downstream (cross section 52.36) of the bridge. The user needs to be cautious placing the pier centerline stations
because the X-coordinates for the upstream and downstream cross section stationing may be different. This is to
ensure that the piers "line up" to form the correct geometry. For this example, the 17 piers are placed 24 feet on
center, with centerline stationing from 2466 to 2850. The piers are 1 foot wide for their entire height. The starting
elevation is set below the ground, and the ending elevation is inside the bridge decking. With these elevations, the
program will "clip off" the excess pier heights. After entering the data, the OK button was selected.
Sloping Abutments. The next icon on the left side of the Bridge/Culvert Data Editor is Sloping Abutment. This
icon was selected to enter the abutment data, and is shown in the figure below for the first (left) abutment. The user
can view the data for the second (right) abutment by selecting the down arrow on the Sloping Abutment Data
Editor in the HEC-RAS program. After the data were entered, the OK button was selected to exit the editor and the
schematic of the bridge/deck, piers, and abutments was displayed as is shown in the figure proceeding the figure
below. (Note: The figure in the text displays the ineffective flow areas that will be added in the next section.)
The cross sections shown in the figure below are developed by superimposing the bridge data on the cross sections
immediately upstream (52.38) and immediately downstream (52.36) of the bridge. The top cross section in the
figure below reflects the geometry immediately inside the bridge on the upstream side and the bottom cross
section reflects the geometry immediately inside the bridge on the downstream side.
While viewing the bridge, the modeler can select to view just the upstream, just the downstream, or both of the
cross section views. This is performed by selecting View and then the required option. Additionally, from the View
menu, the user should select Highlight Weir, Opening Lid and Ground as well as Highlight Piers. These options
enable the modeler to view what the program will consider as the weir length, bridge opening, and pier locations.
Any errors in the data will appear as inconsistent images with these options. Also, the zoom-in option will allow the
user to examine data details. Finally, the Apply Data button was selected to accept the data. (Note: To save the
data to the hard disk, the user must select Save Geometry Data under the File menu on the Geometric Data
Editor).
For this example, ineffective flow areas were included on both the upstream cross section (52.38) and the
downstream cross section (52.36), as well as on cross section 52.29 (the cross section in the expansion reach). To
determine an initial estimate for the stationing of the ineffective flow areas on either side of the bridge, a 1:1 ratio of
the distance from the bridge to the cross section was used. For this example, cross section 52.38 is located 26 feet
upstream of the bridge. Therefore, the left and right ineffective flow areas were set to start at 26 feet to the left and
right of the bridge opening. The left side of the bridge opening is at station 2444, which leads to a left ineffective
station of 2444 - 26 = 2418. The right side of the bridge opening is at station 2864, which leads to a right ineffective
flow station of 2864 + 26 = 2890. Similarly, cross section 52.36 is located 27 feet downstream from the bridge and the
ineffective flow areas at this cross section were set at 27 feet to the left and right of the bridge opening, at stations
of 2417 (= 2444 - 27) and 2891 (= 2864 + 27), respectively.
The initial elevation of the ineffective flow areas for the upstream cross section was chosen as 340.8, a value slightly
lower than the lowest high cord elevation (341.0). This ineffective flow elevation was chosen so that when the water
surface becomes greater than this ineffective elevation, the flow would most likely be weir flow and would be
considered as effective flow.
At the downstream cross section, the elevation of the ineffective flow area was set to be 339.0, slightly lower than
the low cord elevation (340.2). This elevation was chosen so that when weir flow occurs over the bridge, the water
level downstream may be lower than the low cord, but yet it will contribute to the active flow area.
The ineffective flow areas were entered by first selecting the Cross Section icon from the Geometric Data Editor.
Then, Options and Ineffective Flow Areas were chosen for cross section 52.38 and then for cross section 52.36. The
ineffective stations and elevations as previously determined were entered, as shown in the figure below for cross
section 52.38.
Since the "Normal" type of ineffective flow was chosen, these entries imply that all the water to the left of the left
station and to the right of the right station will be considered as ineffective flow until the water level exceeds the
elevation of 340.8 feet.
Cross section 52.29 is in the area of expanding flow and therefore only a portion of the cross sectional flow will be
effective. A "Normal" type of ineffective flow areas were set at stations 1030 and 3800, both at an elevation of 338.
These stations were determined by taking into account the distance of cross section 52.29 from the bridge (426
feet), as compared to the total distance of the expansion reach length (1956 feet). With these distances, cross
section 52.29 is located approximately 25% (= 426 / 1956 * 100%) down the expansion reach length. As the flow exits
the bridge opening, it will initially expand at a greater rate than it will farther downstream. An initial estimate of the
degree of expansion was assumed to be approximately 60% at cross section 52.29. Therefore, 60% of the left
overbank width was subtracted from the left main channel bank station to determine the left ineffective flow
station. Similarly, 60% of the right overbank width was added to the right main channel bank station to determine
the right ineffective flow station. The elevation of 338 was chosen to be greater than the water surface elevations.
The OK button was selected and then Apply Data was chosen to accept the data. The ineffective flow areas then
appeared as triangles, as shown in the previous figure above for the bridge cross-sections. Additionally, the
ineffective flow areas will appear on the cross sections plots. Finally, a note will appear in the box at the bottom of
the Cross Section Data Editor that states an ineffective flow exists for each cross section that this option was
selected for, as shown in the second to top figure on this page. The ineffective flow areas could have been entered
previously along with the discussion of the cross section data (X-Y coordinates, reach lengths, etc.). However, it is
often easier to enter the bridge data first and then the ineffective flow data to easily determine where the
ineffective flow stations should be located. After the steady flow analysis was performed, the ineffective flow areas
were evaluated to determine if they were placed appropriately. This will be discussed near the end of this example.
Chapter 6 of the User's Manual and Chapter 5 of the Hydraulic Reference Manual for additional discussion on the
bridge modeling approach editor.
Low Flow Methods. The modeler needs to select which low flow methods the program should compute and which
method the program should use. The modeler can select to have the program compute particular methods or all of
the methods. Then, the modeler needs to select which method the program will use as a final solution.
Alternatively, the modeler can select the computation of several or all of the methods, and then have the program
use the method with the greatest energy loss for the final solution.
For this example, all four low flow methods were selected to be computed. For the momentum method, a drag
coefficient Cd = 2.00 was entered for the square nose piers. For the Yarnell method, a value of K = 1.25 was entered.
Finally, the WSPRO method was selected to be used for the solution. The modeler is referred to Example 2 - Beaver
Creek for an example application that further discusses the low flow methods of energy, momentum, and Yarnell.
For the WSPRO method, the user is required to enter additional data. This was performed by selecting the WSPRO
Variables button from the Bridge Modeling Approach editor. This activated the WSPRO Data Editor, as shown in
the figure below. Most of the data entered on this editor is used to determine the "Coefficient of Discharge" for the
WSPRO method. The following discussion outlines the WSPRO data that was entered for this analysis. The user is
referred to Chapter 6 of the User's Manual and to Chapter 5 of the Hydraulic Reference Manual for a more
detailed discussion of WSPRO method.
At the top of the editor, as shown in the figure below, the value of 341 feet was entered for the left and right
elevations of the top of the embankments. The elevations for the toe of the abutments were entered as 323.6 and
330.5, for the left and right abutments, respectively. Next, the abutment type was selecting as type 3: sloping
abutments and sloping embankments. The average slope of the abutments was entered as 1 (horizontal to 1
vertical) and the top width of the embankment was entered as 31 feet. All of these values were obtained from the
USGS Atlas.
The centroid stationing of projected bridge opening at the approach cross section was entered as a value of 2531 feet.
To obtain this value, the bridge opening width is projected onto the approach cross section (52.46). This projected
width is situated onto the approach section such that the main channel is in the same relative position to the bridge
opening width on the approach section as at the bride opening. The width of the projecting bridge opening should
embrace the flow, which could pass through the actual bridge opening without contraction. Then, the centroid of
the bridge opening width on the approach section was then determined to be at a stationing of 2531.
The bottom portion of the editor is divided into four sections. The first section is to describe the Wing Walls on the
abutments. For this example, there were no wing walls present. If wing walls were present, the user can select
either "Angular" or "Rounded," by depressing the down arrow adjacent to this field. Then, the appropriate data
would be entered for the angle, length and radius fields, pertaining to the type of wing wall selected.
The second section is to describe the Guide Banks. For this example, there were no guide banks. If guide banks
were present, the options are "Straight" or "Elliptical" to describe the type and then the user would enter the
appropriate dimensions of length, offset, and skew.
The third section at the bottom of the editor is Optional Contraction and Expansion Losses. The WSPRO method
includes friction losses throughout the bridge region for the calculations. However, the only contraction or
expansion loss that is included is an expansion loss from the exit section (cross section 52.29) to the cross section
immediately downstream from the bridge (52.36). The user can include additional losses at the approach section
(52.46), at the guide banks, at the upstream outside contracted section (52.38), at the inside bridge upstream
section (52.37 BU), and at the inside bridge downstream section (52.37 BD). For this example, there were no
additional energy losses selected.
The final section of the editor is used to select options for the piers and friction slope method. If the piers are
continuous, then the user should select Piers are continuous for the width of the bridge. This selection will determine
a weighting factor in the determination of the "Coefficient of Q" discharge factor. Finally, the WSPRO method
traditionally utilizes the geometric mean equation to determine the representative friction slope between two cross
sections. If the Use Geometric Mean as Friction Slope Method is selected, then the program will use this method. If it
is not selected, then the program will use the friction slope averaging method that has been selected by the user.
(This is selected in the Steady Flow Analysis window by choosing Options, Friction Slope Methods, and then the
appropriate choice. The default method is Average Conveyance. The user is referred to Chapter 7 of the User's
Manual and to Chapter 2 of the Hydraulic Reference Manual for more information on friction slope methods.)
This completed the entry of the data for the WSPRO method as the low flow method. The OK button at the bottom
of the editor was selected to exit the editor. Next, the high flow method was selected on the Bridge Modeling
Approach Editor.
High Flow Methods. High flows occur when the water surface elevation upstream of the bridge is greater than the
highest point on the low cord of the upstream side of the bridge. Referring to the seond most recent figure above,
the two alternatives for the program to compute the water surface elevations during the high flows are: Energy Only
(Standard Step) or Pressure and/or Weir Flow. The Energy Only (Standard Step) method regards the flow as open
channel flow and considers the bridge as an obstruction to the flow. The Pressure Flow computations are divided
into two scenarios. The first is when only the upstream side of the bridge deck is in contact with the water and the
second is when the bridge constriction is flowing completely full. The program will begin to calculate either type of
pressure flow when the computed low flow energy grade line is greater than the highest point of the upstream low
cord. Weir Flow occurs when the upstream energy grade line elevation (as a default setting) exceeds the lowest
point of the upstream high cord.
For this example, the high flow method was selected as energy. From the data on the USGS Atlas, the observed
water surface elevations did not encounter the bridge deck, so therefore a high flow method of analysis was not
anticipated. The user is referred to "Example 2 - Beaver Creek" for a bridge analysis that applied a high flow
method.
At this point, all of the bridge data have been entered. The OK button was selected to close the Bridge Modeling
Approach Editor and the Apply Data button was depressed on the Bridge/Culvert Data Editor before it was
closed. Then, the geometry data was saved as the file "Bridge Crossing near Johnston Station" by selecting File and
then Save Geometry Data As from the Geometric Data Editor.
As a final geometric component, a picture of the bridge has been included as a separate file. The red square on the
river system schematic indicates that a picture has been added for that cross section. This picture was scanned
from the USGS Atlas and saved as the file "chit.bmp." To view this picture, the user can select the View Picture icon
from the left side of the Geometric Data Editor. When river station 52.37 is selected, the picture should appear. The
user is referred to Chapter 6 of the User's Manual for a discussion on adding and viewing pictures to the river
system schematic.
surface of 326.0 feet. These values were obtained from the USGS Atlas. The OK button was selected to exit the
boundary conditions editor.
For the purposes of the analysis, if the downstream boundary conditions are not known, then the modeler should
use an estimated boundary condition. However, this may introduce errors in the region of this estimated value.
Therefore, the modeler needs to have an adequate number of cross sections downstream from the main area of
interest so that the boundary conditions do not effect the area of interest. Multiple runs should be performed to
observe the effect of changing the boundary conditions on the output of the main area of interest. For a detailed
explanation of the types of boundary conditions, refer to Chapter 7 of the User's Manual and Chapter 3 of the
Hydraulic Reference Manual.
A final optional component of the steady flow data was to enter the observed water surface elevation data from the
USGS Atlas. This was performed by selecting Options and then Observed WS ... from the Steady Flow Data Editor.
This activated the Observed Water Surface Data Editor, as shown in the figure below. The River was selected as
"Bogue Chitto" and the Reach was "Johnston Sta." The observed values were entered for each of the two profiles,
adjacent to the river stations that observed data was available for, as shown in the figure below. The OK button was
then selected to exit the editor.
This completed the entry for the steady flow data. The Apply Data button was selected on the Steady Flow Data
Editor and the flow data was then saved as "50 and 100 year flows."
Another type of profile table is the Bridge Only table. This table is activated from the Std. Tables menu and is
shown in the figure below, displaying the results for both profiles. The first row is for the first profile and the second
row is for the second profile. The results in the table show the calculated values of the energy grade line at the cross
section immediately upstream from the bridge (EG US) are 336.86 and 338.18 feet for the first and second profiles,
respectively. The column heading Min El Prs shows a value of 340.20, which is the highest elevation of the low cord
on the upstream side. If the energy grade line value had exceeded this elevation, then pressure flow would have
been calculated. The Min Top Rd column shows the minimum elevation of the high cord on the upstream side. If the
energy grade line had exceeded this value, then weir flow would have been calculated. Since the energy grade line
for either profile did not exceed the low cord elevation, pressure flow, nor weir flow developed. This is also
apparent since there are no values in the pressure only water surface (Prs O WS) or the weir flow (Q Weir) columns.
The user is referred to "Example 2 - Beaver Creek" for an example where pressure and weir flow occurred.
To determine the method that was used for the final energy grade line value at the bridge, the Bridge Comparison
table was activated from the Std. Tables menu, and is shown in the figure below. The two rows display the results
for each of the two flow profiles, in ascending order. The river station in the second column is set at 52.37 (the only
bridge location for this reach). The sixth through ninth columns show the results of the low flow methods that were
chosen to be computed : Energy, Momentum, Yarnell and WSPRO methods, respectively. The values of the energy
grade line that the program used is shown in the third column. As shown in Column five, the program used the
result from the WSPRO method, as had been selected in the Bridge Modeling Approach Editor. Finally, the fourth
column shows the water surface elevation that coincides with the energy elevation used by the program. (Note:
There were no results in the last two columns for the high flow methods since only low flow occurred through the
bridge opening.)
The results of the WSPRO method were only slightly greater than the results of the energy method. Both the
momentum and Yarnell methods returned valid results since the water surface did not encounter the low cord of
the deck. The user can adjust the input parameters of the WSPRO method to determine the significance of the
values on the output, such as slope of the abutments and the inclusion of wing walls or guide banks.
The final profile table that will be discussed is the Six XS Bridge table, which is also selected from the Std. Tables
menu. This table is shown in the figure below, for profile 2. This table is specifically designed to show the data for
the six cross sections that are used in the hydraulic computations through the bridge. The water surface and energy
grade line values are shown in the first two columns of data, and are used to determine the water surface profile
through the bridge. The C&E Loss column of the table shows the contraction or expansion losses. As shown in the
figure, the only river station (within the six bridge cross sections) that had an expansion or contraction loss was
river station 52.36 (to river station 52.29). This is because the WSPRO method, by default, only includes an
expansion loss between the exit section and the cross section immediately downstream from the bridge. The user
can include additional contraction/expansion losses by selecting the losses on the WSPRO Data Editor. For this
example, no additional energy losses were selected, so no additional contraction or expansion losses appeared on
the table.
At this time, the user can select to view the data for the first profile as well as viewing the cross section tables for the
cross sections. This completed the review of the output for the cross section tables. Next, the locations of the cross
sections in relation to the bridge were evaluated, by using the information presented in the various tables that were
previously discussed.
where:
= expansion reach length, ft
= expansion ratio
= main channel Froude number at the cross section immediately downstream of the bridge (cross
section 52.36 for this example)
= main channel Froude number at the cross section of fully expanded flow (initially cross section 52.00
for this example)
= = average length of obstruction caused by the two bridge approaches, ft
= total discharge, ft3/s
(Note: The subscripts used in this equation and all subsequent equations reflect the river station numbering for this
example.)
From the field data, the average length of the obstruction is approximately 1770 feet and the total discharge, Q, is
31500 cfs for the high flow event. From the initial analysis, the value of the Froude number at cross-section 52.36
was 0.40 and at cross section 52.00 was 0.27 (as shown on Standard Table 1). Substituting the values into the
equation above yielded that the expansion reach length, Le, was approximately 1874 feet. This equation has a
standard error of 96 feet, which yields an expansion reach length range from 1778 to 1970 feet to define the 68%
confidence band. The distance used for the expansion reach length (the distance from cross section 52.36 to cross
section 52.00) was set to be 1956 feet in the main channel, which is approximately equal to the calculated distance
of 1874 feet. Therefore, the expansion reach length was determined to be set appropriately.
If the reach lengths were not acceptable, then the modeler has the option to adjust the length so that it is within the
calculated range. Then, after a new analysis, the new Froude numbers should be used to calculate a new expansion
reach length. If the geometry is not changing rapidly in this region, then only 1 or 2 iterations should be necessary to
obtain a constant expansion reach length value.
where
= contract reach length, ft
= contraction ratio
= main channel Froude number at the cross section immediately downstream of the bridge (cross
section 52.36 for this example)
= main channel Froude number at the cross section of fully expanded flow (cross section 52.00 for this
example)
= discharge conveyed in the two overbanks at cross section 52.46, cfs
= total discharge, ft3/s
= Manning n value for the overbanks at section 52.46
= Manning n value for the main channel at section 52.46
From the field data and the results of the initial analysis, the Froude numbers at sections 52.36 and 52.00 were 0.40
and 0.27, respectively, the total over bank flow at cross section 52.46 was approximately 25773 cfs (15697 + 10076),
the total flow was 31500 cfs, the weighted n value for both of the overbanks was 0.13, the n value for the main
channel was 0.05 (from the Cross Section Type table), and the average length of the obstruction was 1770 feet.
Substitution of these values into the equation above yielded the contracted reach length of 685 feet. This equation
has a standard error of 31 feet which results in a contraction reach length range from 716 to 654 feet to define the
68% confidence band. For this example, the distance from cross section 52.38 to cross section 52.46 was set at 380
feet along the main channel. An additional analysis was performed with the contraction reach length set at 685 feet
and no appreciable difference was observed in the results. Therefore, the contraction reach length of 380 feet
considered as appropriate for this example.
of the reach lengths and coefficients. This procedure will ensure that the modeler is always using the current flow
data for the analysis.
16.6 Summary
This example demonstrated the use of the WSPRO routines in HEC-RAS to calculate water surface profiles along a
river reach that contained a bridge crossing. The WSPRO routines are used for low flow analysis, and have been
adapted to the HEC-RAS methodology of cross section locations around and through a bridge. The user has the
options to include additional energy losses, as well as account for the type of abutments and the presence of wing
walls and guide banks in the WSPRO routines. By adjusting these input parameters, the modeler can determine the
impact on the calculated coefficient of discharge value and the resulting water surface profiles.
Ice Jam Chan column was set to n for these two river stations. It is important to remember that for every ice jam,
there must be an n set in the Ice Jam Chan column immediately upstream and downstream of the jam.
The nominal Manning's n value was set for the channel ice cover by entering the Manning n value in the Chan ice
Mann n column using the Set Values icon. However, in this case, it is desired that HEC-RAS estimate the Manning's
n of the jam using the procedure outlined in the Hydraulic Reference Manual. It is necessary to indicate to HEC-
RAS that this simulation should estimate the Manning's n of the ice jam and not use the value entered in the Chan
ice Mann n column. This is done in the Fixed Mann n column by changing the default value of y to n. The default
value of y was entered for the upstream and downstream river stations of the Ice Jam Reach, because these
stations are considered to have a known roughness as well as a known thickness.
The other parameters, which describe the ice jam material properties, such as the friction angle, porosity, and
stress ratio, were all left at their default values. In most cases, it will not be necessary to modify these values from
their default values. After all of the ice cover data were entered, the geometry file was saved using the Geometric
Data editor by selecting File and the Save Geometry Data As. The title entered was "Thames River Ice Jam Data."
17.4 Comparison of Open Water, Ice Cover, and Ice Jam Results
To compare the output from the three plans, the user can view the results graphically and in tabular format. In the
interest of brevity, not all the graphical plots that are available will be illustrated here.
To enter the data for the weir, the Weir/Embankment icon was selected from the left side of the Lateral Weir Data
Editor. This activated the Lateral Weir/Embankment Editor as shown in the figure below. This editor is somewhat
similar to the deck/roadway editor used for bridges and culverts.
The right hand side of the editor has two columns for the station and elevation data. These data points define the
top of the lateral weir and they are entered in the upstream to the downstream direction. Similar to cross section
data, the stationing for the data can be based on any horizontal datum. In this example, the weir starts at an
elevation of 90 feet at an arbitrary horizontal station of 100 feet. Fifty feet downstream, the weir drops to an
elevation of 78 feet. After forty feet, it returns to an elevation of 90 feet at station 190. The weir continues at an
elevation of 90 feet until the final station of 200.
The left side of the editor has several fields. The first field is the Weir Width. This is width of the weir in the direction
of flow over the weir (perpendicular to the flow in the river cross section). It is not used for hydraulic computations,
but is used by the GUI for graphical purposes. In this example, it has been set to 10 feet.
The next field is labeled Weir Computations. This field allows the user to select from two different equations to be
used in computing the lateral weir flow: the Standard Weir Equation, and Hager's Lateral Weir equation.
The next lower field is Weir Flow Reference. By left clicking on the down arrow, the modeler can toggle between
"Energy Grade" and "Water Surface." Based on this switch, the program will compute the flow over the weir (and
through any gates), using either the energy grade or the water surface when calculating the head (i.e. the depth of
flow) on the weir. In this example, since the lateral weir is located in the overbank, the energy method was selected.
The next field is the Weir Coefficient (Cd). It is a coefficient that is used in the computation of weir flow. It has been
entered as 3.0.
The next field is Weir Crest Shape, and the user can select one of the following four shapes: Broad Crested; Ogee;
sharp crested; or zero height weir. For this example, broad crested has been selected. If the Ogee had been selected,
two additional fields and a Cd button would have popped up allowing the user to enter a Spillway Approach Height
and a Design Energy Head. For the Ogee weir, the user can either enter a weir coefficient or have the program
compute a weir coefficient by clicking on the Cd button.
The last field labeled HW Distance to Upstream XS, is used to locate the lateral structure along the stream spatially.
The upstream cross section is the first normal cross section immediately upstream of the lateral weir based on river
stationing. (Between two normal river cross sections, there can be one or more lateral weirs.) In this example, the
upstream cross section is 1188 and the upstream end of the lateral weir is 10 feet from this cross section. So the HW
Distance has been entered as 10 feet. The program will compute the distances to downstream cross sections based
on where the lateral weir is located. In this example, the HW Position of the lateral weir has been specified as the
right overbank (see above). This means that the program will use the right overbank reach lengths. The right
overbank length (between cross sections 1188 and 1108) is 90 feet. Therefore, the lateral weir intersects the
downstream cross section (river station 1108) twenty feet from the downstream end of the lateral weir.
It should be noted that the lateral weir cannot cross more than 498 cross sections. If this happens, the user should
break the lateral weir into two or more weirs that are physically adjacent.
The remaining portion of the editor is divided into two sections, one for the gate flow data and one for weir flow
data. The gate flow data are used when the water surface upstream of the gate is greater than 1.25 times the gate
opening (as measured from the gate invert). At this water surface elevation, the gate is in contact with the water
and is controlling the flow rate. The weir data are used when the upstream water surface is less than or equal to the
gate opening. At this water surface elevation, the weir under the gate is controlling the flow through the gate
opening (i.e., the water is not in contact with the gate). In between these two elevations, the flow is in a transition
zone.
For the gate flow data, the first field is the Gate Type. By selecting the down arrow, the type "Sluice" was chosen.
When the gate type was selected, the Trunnion Exponent, Opening Exponent, and Head Exponent values were
automatically set to 0.0, 1.0, 0.5 respectively. Since a sluice gate has been selected, the Trunnion Height field is not
visible. The Discharge Coefficient was entered as 0.6. This coefficient is used when the gate is experiencing sluice
flow (downstream end of the gate is not submerged). The orifice coefficient of 0.8 was entered for full flow
conditions (this is used when the tailwater elevation on the gate causes it to be submerged).
For the weir flow data, the Shape was selected as "Broad Crested." The weir coefficient was left as 3. This
information appears in the weir data area at the bottom of the Gate Editor.
This completed the data entry for the gates. The OK button was selected at the bottom of the Gate Editor and the
gates appeared on the Lateral Structure Data Editor as shown in the third most recent figure. At this point, the
user should zoom in on the gate openings to ensure that they do not overlap and appear as intended. The Lateral
Structure Data Editor was then closed.
At the top portion of the editor, the River "Spruce Creek," Reach "Upper River" and the River Station "1150" were
selected. The Description is the same as was entered in the Lateral Spillway Data Editor (shown on previous page).
The # Gate Groups field shows that there is only 1 gate group at this river station. The table in the central portion of
the editor has a single row for the one and only gate group. The first column lists the description for the gate group,
as it was named in the Gate Editor. The second column displays the number of gate openings for the gate group (3).
The third column displays the maximum gate height for the gate group (5 feet for this example).
The remaining portion of the editor consists of entry fields for the number of gates opened and the opening heights
of the gates for each flow profile. For the Lateral weir, all of the gates in a gate group must be operated together. For
this example, for profile 1, the three gates were opened 3 feet. To close the gates, the user could set the number of
gates opened to zero, or, alternately, the gate height could have been set to zero. However, if the user enters a non-
zero (i.e. positive) number for the number of gates and the gate height, then the program will open all of the gates
in the gate group to that height (the field for the number of gates is set up for the Inline weir where it is possible to
only open part of the gates in a gate group). This means that if one gate is opened four feet, then all of the gates will
be opened four feet. If the user wanted to only open one gate, then that gate would have to be defined as a separate
gate group. A Lateral Weir can have a maximum of ten gate groups. So ten gates (or groups of gates) could be
operated independently. If more than ten gate groups are needed, the lateral weir could be defined as two or more
parts. For instance, a thousand foot long lateral weir could be defined as two 500 foot long lateral weirs. Each "weir"
would have to be entered separately with its own river station and data. By choosing the appropriate Distance to
upstream cross section, the two weirs could be placed immediately adjacent to each other. This would allow for
twenty different gate groups to be operated independently.
The user can toggle across the table to view the number of gates open and the gate opening heights for all of the
profiles. During the analysis of the output, the various gate settings will be discussed. This concludes the data entry
for this example. At this point, the OK button at the bottom of the editor was selected and the flow data was saved
as "3 Flow Profiles."
Next, optimization for the junction was turned on. This was done by clicking on the Junction tab, which changes the
Optimization Editor as shown in the figure below. Optimization at the Meadows junction can then be toggled like it
was for the lateral weir. Flow optimization can only be performed at junctions that have more than one
downstream reach. For this reason, the Editor does not display the Pottsville junction. If the optimization flag at the
Meadows junction had been left off, the program would not optimize the flow split. Rather, the program would
maintain the same flow ratio in each reach.
For instance, assume that for the first profile, the junction optimization is off and that the upstream lateral weir has
an outflow of 500 cfs. This means that there would be 1,000 cfs flowing into the Meadows junction instead of 1,500
cfs. The program would then proportionately reduce the flow in each reach. Middle River would be 200 cfs (instead
of 300 cfs) and Bryon Creek would be 800 cfs (instead of 1,200 cfs).
Now that the optimization flags have been turned on, the files for this plan can be named and saved. First, the
geometry file was selected as "Lateral Weir with Full Looped Network" and the flow file was "3 Flow Profiles." Next,
the Flow Regime was selected as "Subcritical." Then, File and Save Plan As were chosen and the information was
saved as the plan "Split Flow." This plan name then appeared on the Steady Flow Analysis Window, as well as on
the main program window. Finally, the COMPUTE button was clicked to perform the analysis.
In order to perform the backwater computations (and get the water surfaces and energies in the vicinity of the
weir), the program must know the amount of flow in the river. To do this, it assumes the amount of flow that will be
diverted by the lateral weir. This is based on the previous iteration (for this example, for the final iteration, a flow of
2612 cfs was assumed to be diverted). After the backwater calculations have been performed, the program can
compute a flow for the lateral weir based on water surfaces and energies (in this case, 2594 cfs). If these two
numbers are within the default tolerance of 2% (2% of the assumed 2612 cfs or 52 cfs), then the lateral weir is
assumed to have converged. The default tolerances can be changed from the Steady Flow Editor by choosing
Options and then selecting Set Calculation Tolerances (this editor is not shown).
4000 cfs flow at the start of the Middle River reach. Hence, the program adds 500 cfs flow at this river location.
for the lateral weir flow can be entered. The program will then assume that this much flow is diverted when it
performs the first iteration of the backwater calculations.
18.6.1 Summary
This example computed 3 flow profiles for the Bryon Creek and Spruce Rivers. It included a lateral weir with gates
and a fully looped river network between the Meadows and Pottsville junction.
By reviewing the lateral weir and junction tables and the water surface profiles and tables, the user can determine
the amount of flow being diverted by the lateral weir and the amount of flow being split into the two reaches
downstream of the Meadows junction. The Lateral Weir/Spillway Output table provides detailed output for the
lateral weir, for a given profile.
The Junctions table provides flow split output for one or more profiles. The user can also adjust the initial starting
conditions to improve the run performance of the computer computations.
The channel modification began with the River set to "Critical Creek" and the Reach set to "Upper Reach." The
Starting River Station was set by default to 12 and the Ending River Station was changed to a value of 1, the
downstream end of this particular reach. The Project cut from upper RS at slope option was entered as 0.01. This
option projects a cut using the specified slope from the invert elevation of the Starting River Station.
Once the correct river stations were set, the first cut was established. The Center Cuts (y/n) column in the Set
Range of Values table was entered as "y." This set the centerline for the first cut at the centerline of the main
channel. Next, "100" was entered in the Bottom Width column to set the width of the first cut to 100 feet along the
centerline. The Invert Elevation column was left blank for the first cut, which defaults the program to the existing
invert elevation of the Starting River Station. Both side slopes for the banks were entered as 2 (2 horizontal to 1
vertical). The Cut n/K column for the new Manning's n-value was entered as 0.025.
The second cut was also done on the existing centerline of the main channel. A Bottom Width of 400 feet and an
Invert Elevation of 1810 feet were entered. The entry initiated the second cut at an elevation of 1810 feet, at the
Starting River Station. The cut will then project with the specified slope from this elevation. The slopes were again
entered with a value of 2 and Manning's n-value was entered as 0.03. After the data was entered, the Apply Cuts to
Selected Range button was pressed to compute the data shown in the lower table of the figure above.
An option available to the user is Cut cross section until cut daylights. For this particular example this option was
selected. As the program performs the cutting of the trapezoidal channel, the left and right banks of the channel will
initiate at the invert elevation and cut through the ground until they reach open air, then the cutting will stop. If this
option is turned off, the left and right banks of the trapezoid will be projected to infinity, continually cutting any
ground that lies above them.
This procedure was repeated for the modified channel geometry. On the Steady Flow Analysis Window the
Geometry File was selected as "Modified Geometry" and the Steady Flow File remained "100 Year Profile." The
Flow Regime was changed to "Supercritical." The title was entered as "Modified Conditions Run" and the Short
I.D. was entered as "Modified." The file was saved and then the Compute button was selected.
After pressing the OK button on the Plan Selection Window, the figure below was displayed. The figure shows the
two plans for the existing channel geometry and the modified channel geometry. It can be seen from the profile that
the existing channel geometry had a flow that was mixed between the subcritical and supercritical regime. The
modified channel geometry altered the flow to be exclusively supercritical for this reach.
19.4 Summary
The geometry of Example 1 was modified to prevent the flow from the 100- year event from overflowing the
channel. This was accomplished by modifying the existing channel conditions to include two cuts down the
centerline of the channel. By reviewing the water surface profiles and tables, the user can determine the benefits of
a specific channel modification for a given flow.
In addition to these features, a culvert is located at River Station 3.395 on the Northwest Reach and a bridge is
located at River Station 1.895 on the South Reach. The cross section geometry is prismatic with little change in
elevation with distance except for the North Reach. The level topography provides an excellent environment to
examine unsteady flow for such aspects as flow reversal.
This example assumes that the user has sufficient experience constructing the following river system features and
all attributes associated with each feature.
• Reaches
• Cross Section Geometry
• Junctions
• Bridges and Culverts
Consult with previous examples for information pertaining to the above mentioned features. The following
description commences with the creation of storage areas.
The remaining three storage areas were created in the same manner as Pyramid Lake without any connections to
reaches. The figure below displays the completed storage areas in the Geometric Data Editor before the addition
of the hydraulic connections. The hydraulic connections for each storage area are discussed in impending sections.
An additional option for calculating storage volume is the Elevation versus volume curve method, shown in the
figure below for Eagle Lake. For this method the initial elevation was entered as 0 feet with a volume of 0 acre-feet.
Data was then entered in increments of 2 feet for the cumulative volume of storage produced at each increment. It
is important to note that the user is not entering the surface area of storage at each increment but the cumulative
volume of storage at each increment.
Another lateral structure was added from the southeast reach to the Eagle lake storage area. This lateral structure
was added in the same manner as the previous one, except it had a culvert in addition to the overflow weir.
the From location, and Eagle lake as the To location. Next the Weir/Embankment icon was selected. A 500 ft long
broad crested weir was entered at an elevation of 16.5 feet.
An opening in the embankment was placed along the centerline of the weir by selecting the Culvert icon on the
Storage Area Connection Data Editor. This selection displayed the Culvert Data Editor shown on the figure
below. A circular culvert with a diameter of 5 feet and an invert elevation of 11.5 feet was entered for this hydraulic
connection. Additional parameters entered are shown in the figure below.
One additional storage area connection (Upper to Lower lake) was constructed in the same manner to connect
those two storage areas.
manually entered with a baseflow of 100 cfs and a floodwave that peaked at 5000 cfs.
The user may also select to read hydrograph data from a DSS file. To do this the user presses the button Select DSS
file and Path. When this button is pressed a DSS file and pathname selection screen will appear as shown in the
figure below. The user first selects the desired DSS file by using the browser button at the top. Once a DSS file is
selected, a list of all of the DSS pathnames within that file will be displayed in the table. The user can use the
pathname filters to reduce the number of pathnames shown in the table. When the desired DSS pathname is found
the user closes the window and the filename and pathname will be recorded in the Flow Hydrograph Window.
Flow data was entered for each reach so the program could perform a steady-flow backwater run to compute the
corresponding stages at each cross section. As shown on the figure above, the North Reach has an initial flow of 100
cfs, corresponding to the baseflow of the hydrograph entered as the upstream boundary condition. The other
various reaches then split the 100 cfs of flow down to the South Reach where the flow combines to 80 cfs. In
addition to flow data, initial elevation for each storage area is needed. Initial elevations were set to 11, 10, 12, and
10 feet for Pyramid Lake, Eagle Lake, Lower Angora, and Upper Angora, respectively. The data was then saved in the
Unsteady Flow Data Editor and the window was closed.
can be reduced by only selecting the most essential cross sections for output.
The Geometry Pre-processor box under "Programs to Run" was selected in the Unsteady Flow Analysis Window.
The Unsteady Flow Simulation and Post-Processor boxes were unselected (no checkmarks) and the Compute
button was depressed. After the geometry pre-processor finished processing the data, Hydraulic Property Plots
was selected from the View menu on the main program menu. The figure below displays a typical cross section plot
of hydraulic properties for this example. Cross sections are processed into tables of elevations versus hydraulic
properties of areas, conveyances, and storage. Viewing these plots for anomalies is a good diagnostic tool to search
for cross sections with irregular geometry.
The geometry pre-processor evaluates hydraulic structures such as bridges and culverts and relates the structures
as a function of tailwater, flow, and headwater. From the Type menu in the Hydraulic Properties Tables, Internal
Boundaries was selected. Next, the South Reach was selected displaying the figure below, the family of rating
curves for the bridge located at RS 1.895. On the plot the free flow rating function describes the flow if tailwater
submergence does not occur, such as free flow over a weir. Above the free-flow rating function is a family of
submerged flow rating curves, one for each tailwater elevation.
As shown in the figure above, there is a transitional area to left of the free flow curve between 10 and 12 feet of
elevation where the lines of constant tailwater extend vertically. This occurs because at this water elevation flow
begins to contact the upstream side of the bridge, causing backwater. This is a transition zone where free surface
flow changes to orifice flow. This type of flow is unpredictable because the flow is changing from free surface flow
to a "sluice gate type" of pressure flow and then possibly to a full flowing orifice flow.
First, the "Cross Section" option was selected from the Type menu and the reach was set to West, displaying the
figure above. As seen from the figure, the flow was initially negative, denoting that at the beginning of the
simulation time flow was moving away from Pyramid Lake. As the floodwave progressed the flow changed to
positive, implying a reversal in the direction of flow towards Pyramid Lake. After the floodwave passed, the
direction of flow returned to flowing away from Pyramid Lake. The stage and flow can also be viewed as tabular
output by selecting the Table tab located on the Stage and Flow Hydrograph Window. The output will be
displayed in the increment set on the hydrograph output interval located in the Unsteady Flow Analysis Window.
The geometric data set was constructed with relatively no slope to emphasize the ability of RAS to model unsteady
flow, including flow reversals. As seen in the previous figure for the Southwest Reach, the flow reverses direction
during the peak of the floodwave. This flow reversal occurs because water is diverted to the West Reach and to the
lateral weir on the Southwest Reach. These diversions decrease the flow in the Southwest Reach compared to the
flow in the Southeast Reach. The discrepancy in flows between the two reaches causes a significant difference in
water surface elevations at the Lower Junction. The difference in water surface elevation forces water to move
upstream on the Southwest Reach during the floodwave.
Next, the Lateral Structure option was selected from the Type menu and the reach was set to Southwest,
displaying the figure below. As shown in the figure below, the lateral weir is affected by the flow reversal on the
Southwest Reach. Additionally, the tailwater stage for the lateral weir rises above the headwater stage twice during
simulation, causing flow to enter the reach from the lateral weir.
In addition to viewing output directly from the Stage and Flow Hydrograph Window the user can open the
Geometric Data Editor and select to view output by clicking on the desired feature. From the Geometric Data
Editor the storage area "Pyramid" was selected by clicking on it with the mouse. Plot Stage and Flow Hydrograph
was selected, displaying the figure below.
As shown in the figure below, the stage of the storage area steadily decreases because flow is leaving the storage
area until the floodwave passes through, causing the stage to increase. The stage of the storage area then begins to
decrease again after the floodwave.
Finally, from the Geometric Data Editor the storage are connection "Pyramid to Eagle" was selected, displaying
the figure below. Initially, there is no flow in this connection. This occurs because the stage is below the low invert
elevation of the culvert at 11.5 feet. When the stage increases to over 11.5 feet the flow rate steadily increases, as
shown on the figure below.
20.5 Summary
The concept of unsteady flow analysis within a network of channels and storage areas was discussed. This example
demonstrates the ability of HEC-RAS to route a hydrograph through a network of channels containing junctions,
culverts, bridges, storage areas, lateral weirs, and hydraulic connections.
nsteady flow analysis can be extremely difficult compared to steady flow analysis because input parameters can
cause instabilities in calculations. Instabilities can cause the program to fail to converge on a solution. It is highly
recommended that the user have experience with unsteady flow modeling. Refer to Chapter 8 of the User's Manual
for additional information on unsteady flow analysis.
structures connected to a storage area see "Unsteady Flow Application – Example 17.(see page 240)"
Since there are no flap gates on the culverts, No Flap Gates was selected from the All Culverts dropdown menu in
the center of the Inline Structure Data Editor. "Flaps prevent Negative Flow" can be selected in cases where flap
gates prevent the inline structure's tailwater from propagating upstream through the culverts. (e.g. Immediately
downstream of the inline structure the river empties into a bay whose tidal fluctuations are impeded at the culverts
by flap gates.) "Flaps prevent Positive Flow" prevents the culverts from passing flow from upstream to downstream;
however, flow can pass from downstream to upstream.
21.2 Weir/Embankment
To enter the data for the weir, the Weir/Embankment icon was selected from the left side of the Inline Structure
Data Editor. This activated the Inline Structure Weir Station Elevation Editor as shown in the figure below. This
editor is similar to the deck/roadway editor used for bridges and culverts. Further insight and instruction on the
Weir/Embankment function can be found in "Inline Structure – Example 12(see page 172)" and Chapter 6 of the User's
Manual.
When the Inline Gate Editor was activated, the first gate was automatically labeled as "Gate #1." The Height, Width,
and Invert for the gate were entered as 5, 10, and 210 feet, respectively. On the lower left side of the editor, the
centerline station for the gate was entered as shown in the figure above. The # Openings automatically updates to
coincide with the number of centerline stations entered. In this case there is only one opening.
The remaining portion of the editor is divided into two sections, one for the Gate Flow data and one for Weir Flow
data. For the gate data, the first field is the Gate Type. For this example the gate type was left at the default type
"Sluice". The Sluice Discharge Coefficient was entered as 0.6, while the default value of 0.8 was used for the Orifice
Coefficient. The orifice coefficient is only used for fully submerged flow conditions. Additionally, the head reference
for the gate equation defaults to the "Sill (invert)" of the gate opening.
For the Weir Flow data, the Weir Shape was kept at the default shape of "Broad Crested" and the Weir Coefficient
was kept at the default value of three.
The OK button was selected at the bottom of the Gate Editor and the gate appeared on the Inline Structure Data
Editor as shown in the second figure on this page. (Note: The ineffective flow areas will be entered after entering all
the geometric data for the inline structure components.) At this point, the user should zoom in on the gate opening
to ensure that they do not overlap with other features and appear as intended. The Inline Structure Data Editor
was then closed.
21.3 Culverts
To enter the data for the culverts, the Culvert icon was selected from the Inline Structure Data Editor (second
figure on this page). This activated the Culvert Data Editor as shown in the figure below. The Culvert Data Editor
for the inline structure is identical to the one used under the Bridge/Culvert icon found in the Geometric Data
Editor. In this example two culverts were used, one circular culvert and one box culvert. For more information on
entering culvert data refer to "Single Culvert (Multiple Identical Barrels) Example 3(see page 47)" and "Multiple
Once the Outlet RC icon was selected, a window appeared that provided a table to enter the rating curve
information (below). The radio button Outlet flow based on water surface upstream was selected to signify that
the relationship between the diversion flows from the structure has a direct correlation to the upstream water
surface elevation. If the diversion flow was instead based on the upstream flow rates then the Outlet flow based on
flow upstream radio button would be selected, causing the heading of the first column of the table to change from
US WS Elev to US Flow.
The outlet rating curve table was filled in with the values seen in the figure above. It should be noted that if the
calculated water surface elevation is outside the upper or lower bounds of the rating curve the program uses the
two nearest points on the rating curve to linearly extrapolate a diversion flow value to correspond to the calculated
water surface elevation. Once the values are entered in the table, the resulting rating curve can be plotted and
graphically inspected by clicking Plot Curve. The outlet's rating curve is plotted in a separate window as seen in the
figure below. After filling in the table and viewing the plot, click the OK button in the lower right hand corner of the
Inline Structure Outlet Rating Curve window.
Note
The Outlet Rating Curve option only supports single value rating curves and non tailwater influenced outlets. If
you have a tailwater influenced outlet that requires a family of curves, use the "User Define Curves" option
within the Gate editor.
The time series outlet is time dependent; therefore, the operation of the outlet is event specific and cannot be
entered in as geometric data like the other inline features. Instead, a flow hydrograph is created at the inline
structure cross section as an internal boundary condition in the Unsteady Flow Editor. The name given to the time
series outlet (i.e. "HydroPower") is not seen on the Unsteady Flow Editor; however, the title does appear on the
various output tables once the model is computed. How to enter the flow hydrograph in the Unsteady Flow Editor
will be explained further in the next section.
appropriate units.
window was closed and the data was stored by clicking OK.
Next, the corresponding cell for the other internal boundary condition (also at RS 5.26 IS) was highlighted under the
Boundary Condition column. The Flow Hydrograph was chosen to represent the time series outlet titled
"HydroPower". The Flow Hydrograph data entry window that appears has the same inputs as the T.S. Gate
Openings boundary condition just completed. The same parameters were entered for both with the exception of
the flow values under the Flow column under the Hydrograph Data table. Since the hydrograph represents the flow
of a hydroelectric facility, there is some daily variability in the flows (as seen in the figure below). In this example,
there is zero flow during the hydro facility's non-peak hours (8 p.m. – 6 a.m.), 150 cfs flow from 6 a.m. – 3 p.m., and
during peak energy demand hours (3 p.m. – 8 p.m.) there is 220 cfs diverted through the hydropower facility. Once
all the data was entered, the window was closed and the data was stored by clicking OK.
As seen in the figure above, the Enter Initial flow distribution radio button was selected and an initial flow of 3000
was entered for the most upstream river station (5.99). In addition to flow data, the initial elevation for the storage
area is needed. The initial elevation of the storage area was set to 200. The data was then saved (File -> Save
Unsteady Flow Data As…) in the Unsteady Flow Data Editor and the window was closed.
of the weir and gate structures. The Inline Structure Output table can be seen in the figure below.
21.8 Summary
This example computed an unsteady flow event through a complex inline structure on Beaver Creek, which
included an inline weir, culverts, gated spillway, outlet rating curve, and a time series outlet. The example
demonstrated how to make use of the various inline features noted above.
In addition, a brief overview was given on the output profiles and tables that can be utilized in evaluating the
operation of the inline structure. The output tables and profiles described were the Inline Structure Output table,
Profile Summary Output Table - Inline Structure, and Profile Plot.
As seen in the figure above, the Enter Initial flow distribution radio button was selected and an initial flow of
1075.53 was entered for the most upstream river station (138154.4). The data was then saved (File -> Save
Unsteady Flow Data As) in the Unsteady Flow Data Editor and the window was closed.
After all the necessary data was entered, the Compute button was selected to perform the steady flow analysis.
The cross sections for which Modified Puls routing was to be performed were selected by clicking the button Add
Region and then choosing the upstream and downstream bounds of the river and reach which the routing was to
be carried out (see figure below). Once the desired region was selected, OK was clicked. The region selected was
automatically populated in the Modified Puls Hydrologic Routing table as seen in the figure above.
In order to populate the last column (RC's Imported), rating curves needed to be imported for the selected region.
Below the table there is an area titled Import Rating Curves for Selected Regions where a list of steady flow
output plans can be selected from a drop down menu. For this example, "SteadyNoDam" was selected and the
number of profiles within the plan (30) was listed next to the drop down menu (see the first figure on this page).
Next, the blank cell under the column heading RC's Import was selected. Then, Import Rating Curves (RC's) from
Steady Flow Output was clicked and rating curves from the steady flow output were imported for all 116 cross
sections within the designated region. The imported rating curves can be visually inspected by clicking Plot in the
lower left hand corner of the window. It should be noted that the steady flow models should have been run
previous to this step and that the same geometry file must be used.
Before clicking OK and exiting the window, the box corresponding to Tail Water Check was checked on. When
turned on, this option will monitor the downstream interface of any hydrologic routing reach to verify that the next
cross section downstream has a lower computed water surface elevation than the last section of the hydrologic
routing reach. If it does have a lower computed water surface elevation then the results are accepted and displayed.
However, if the downstream water surface (tailwater) is higher than the water surface of the last cross of the
hydrologic routing reach, then the tailwater is forced into the downstream portion of the hydrologic routing reach.
When run, the model will continue to overwrite the hydrologic routing cross sections' computed water surface
elevations until a cross section is reached where the hydrologic routing computed water surface elevation is higher
than the tailwater elevation.
Below are the water surface profiles from a similar geometry file as used in this example where the Tail Water
Check was not used (top figure below) and where it was (bottom figure below). There is an inline structure in the
middle of the reach that is causing a backwater effect upstream of the structure. As seen in the figure below, this
backwater effect is only illustrated in the portion of the reach that is modeled as unsteady flow, while the portion
modeled with Modified Puls does not consider the backwater effects of the inline structure, which causes the huge
stage difference seen in the top figure below. By turning on the Tail Water Check, this unrealistic water surface
profile is corrected by propagating the backwater effect upstream into the reach region that was modeled with
The Tail Water Check allows for downstream conditions to influence the water surface, volume, and flow rate in the
hydrologic routing reach. Although the Tail Water Check function may give a more realistic result, it may also
increase instability. It should also be noted that there must be at least one or two normal cross sections between
the end of the reach and the Modified Puls region; if there isn't then the tailwater check will be ignored. For this
specific example, the Tailwater check option was not used, because there were no downstream structures that will
cause a backwater upstream in this example.
The geometry pre-processor is used to speed up the unsteady flow calculations by processing the geometric data
into a series of hydraulic property tables and rating curves. It is highly recommended that the user run the geometry
pre-processor and examine hydraulic output for anomalies before running the unsteady flow simulator and the
post-processor. The unsteady flow simulation is what actually creates the stage and flow hydrographs. The
program is a three step process that requires the tables from the geometry pre-processer, the boundary conditions,
and flow data from the interface to operate.
The post processor is used to compute detailed hydraulic information for a set of user specified time lines during
the unsteady flow simulation period. By running the post processor, the user will have all of the available plots and
tables for unsteady flow that HEC-RAS normally produces for steady flow.
The Geometry Preprocessor, Unsteady Flow Simulation, and Post Processor were all checked on.
Transition from Modified Puls routing (upstream) to full unsteady flow (downstream).
22.7 Summary
The Hydrologic Unsteady Routing option is very useful when encountering portions of the model that are very steep
and full unsteady flow routing is either unstable or not possible at all. To perform the analysis the user must enter
the geometric data for the river reach. Within the Geometric Data Editor the hydrologic unsteady routing
parameters are entered. A steady flow model must be created using the exact same geometry file to be used for the
unsteady flow modeling and should include a set of profiles that includes all possible flows within the unsteady
flow model. The rating curves from the steady flow output are imported into the Modified Puls Hydrologic Routing
Editor so that stages and flows can be interpolated when it is evaluated in the Unsteady Flow Analysis. Although
unsteady modeling is more accurate than Modified Puls routing, there is typically very little variance between the
results.
The left side of the editor has several fields and is where Hager's Weir Equation Parameters are entered. The first
field is the Weir Width. This is width of the weir in the direction of flow over the weir (perpendicular to the flow in the
river cross section). With the exception of calculating the broad crested weir coefficients with Hager's Equation, the
weir width is not used for hydraulic computations, but is used by the GUI for graphical purposes. In this example, it
has been set to 10 feet.
The next field is labeled Weir Computations. This field allows the user to select from two different equations to be
used in computing the lateral weir flow: the Standard Weir Equation and Hager's Lateral Weir equation. Hager's Eqn
was chosen.
There are four entries for the Hager's Weir Equation Parameters field: Default Weir Coefficient (Cd), Weir Average
Height, Average Bed Slope, and Weir Angle in Degrees. The default weir coefficient (Cd) is used for the first iteration
of calculations when computing Hager's equation. Although the weir value initially entered for Hager's equation
will not be exactly the same as the final value used after the equation goes through the iterative process, it should
represent a reasonable approximation. Typically, the default weir coefficient is lower than the prescribed
coefficients for weirs of the same type. In this example, 2.2 is used for a sharp crested weir whereas the Hydraulic
Reference Manual prescribes a weir coefficient of 3.1 – 3.3 for a sharp crested weir (Table 8-1 Typical Overflow Weir
Coefficients). The coefficient is lowered because the weir equation assumes that the water flows perpendicular to
the weir face, which is not the case for lateral weirs. Therefore, the coefficient is lower than normal to accurately
account for the lower flows that will occur due to the variance in flow direction.
The Weir Average Height is the average height of the weir above the ground. The average height of the weir was 10
feet in this example.
The Average Bed Slope and Weir Angle in Degrees are both optional fields. If the user does not enter the average
slope of the stream bed, then the program will compute the slope by estimating an average bed elevation for each
cross section and then compute the slope from the average bed elevations. The Weir Angle in Degrees value is
assumed to be zero (meaning the weir is parallel to the stream) if nothing is entered. If the weir is not parallel to the
stream but instead angled towards the center of the river, an angle is required in the field (see the figure below to
The next field is Weir Crest Shape, and the user can select one of the following four shapes: Broad Crested; Ogee;
Sharp Crested; or Zero Height weir. For this example, sharp crested has been selected. If the Ogee had been
selected, two additional fields and a Cd button would have popped up allowing the user to enter a Spillway
Approach Height and a Design Energy Head. For the Ogee weir, the user can either enter a weir coefficient or have
the program compute a weir coefficient by clicking on the Cd button.
The next category is Weir Stationing Reference. The field HW Distance to Upstream XS was set to zero feet, signifying
that the lateral weir starts immediately at the upstream cross section. The field titled "TW flow goes" was left at the
default setting and should have no influence on the results because the tailwater was set to go "Out of the system"
on the Lateral Structure Editor.
The right side of the editor has two columns for the station and elevation data. These data points define the top of
the lateral weir and are entered in the upstream to the downstream direction. Similar to cross section data, the
stationing for the data can be based on any horizontal datum. In this example, the weir starts at an elevation of
11.13 feet at an arbitrary horizontal station of zero and concludes at an elevation of 10.94 feet at station 100.
The bottom left side of the editor has two buttons labeled "HW Connections" and "TW Connections", which are used
to locate the lateral structure along the stream spatially. The upstream cross section is the first normal cross
section immediately upstream of the lateral weir based on river stationing. (Between two normal river cross
sections, there can be one or more lateral weirs.) In this example, the upstream cross section is 600* and the
upstream end of the lateral weir is zero feet from this cross section. So the HW Distance has been entered as zero
feet. The program will compute the distances to downstream cross sections based on where the lateral weir is
located. In this example, the HW Position of the lateral weir has been specified as next to the right bank station (see
above). This means that the program will use the main channel reach lengths of 50 feet. Therefore, the lateral weir,
which is 100 feet long, intersects the downstream river station 500 and extends no farther. It should be noted that
the lateral weir cannot cross more than 98 cross sections. If this happens, the user should break the lateral weir into
two or more weirs that are physically adjacent. However, breaking the lateral weir into two or more weirs might
cause the program to take longer to converge during the iterative process.
compare the weir coefficient values calculated by Hager's equation for each time step.
23.5 Summary
This example used Hager's Lateral Weir Equation to compute flow profiles for a designed rectangular channel and
sharp crested lateral weir. It was seen how the Hager's Lateral Weir Equation adjust the weir coefficient for each
profile to account for difference in hydraulic conditions. By reviewing the lateral weir tables and the water surface
profiles and tables, the user can determine the amount of flow being diverted by the lateral weir and evaluate other
hydraulic properties that are influenced by the presence of the weir.
station 1000 as the upstream cross section and zero as the downstream cross section.
To enter the data for the weir, the Weir/Embankment icon was selected from the left side of the Inline Structure
Data Editor. This activated the Inline Structure Station Elevation Data Editor as shown in the figure below. This
editor is similar to the deck/roadway editor used for bridges and culverts.
The top row of the editor consists of three data entry fields. For the first field, the user must enter the Distance from
the upstream cross section (RS 300) to the upstream side of the weir. For this example, this distance was 10 feet.
Next, the Width of the weir was entered as 80 feet. This is a total distance of 90 feet. The distance from river station
300 to 200 is 100 feet. Therefore, the distance from the downstream end of the weir to cross-section 200 is 10 feet.
The last field in the top row of the editor is the Weir Coefficient and was entered as 2.6, coinciding with the Broad
Crested weir shape chosen at the bottom of the editor. The central portion of the editor consists of a table in which
the user must enter the station and elevation data for the weir. The first station of the weir was entered as zero and
the last station is at 210 feet – both have an elevation of 95 feet. With these weir stations and elevation, the program
will block out the entire area below the weir crest.
For Gate Group #1, the Height, Width, and Invert of the two openings were entered as 6, 6, and 89 feet, respectively.
On the lower left side of the editor, the centerline stations for the two gates in group one were entered as 110 and
140. As these values were entered, the counter field # Openings increased to represent the total number of gates for
the group (two for this example). The Weir Shape was chosen as "Broad Crested" from the dropdown menu and its
coefficient was entered as 2.6. Unlike the closed top overflow gate, the open air overflow gate does not require the
user to enter values for sluice gate flow or submerged orifice flow because it only experiences weir flow.
The information for gate group two was entered by first selecting it from the Gate Group dropdown menu in the top
left corner of the Inline Gate Editor. The Height, Width, and Invert of the closed top overflow gate (chosen from the
Gate type dropdown) were entered as 5, 4, and 87 feet, respectively. The centerline station is 125.
The remaining portion of the editor is divided into two sections, one for the Gate Flow data and one for Weir Flow
data. The gate data is used when the water surface upstream of the gate is greater than 1.25 times the gate opening
(as measured from the gate invert). At this water surface elevation, the gate is controlling the flow rate. The weir
data is the shape of the weir under the gate and is used when the upstream water surface is less than or equal to the
gate opening. At this water surface elevation, the weir under the gate is controlling the flow through the gate
opening (i.e., the water is not in contact with the gate).
Since the first field - Gate Type - has already been chosen as Overflow (closed top) the fields below it should have
already been visible. The Sluice Discharge Coefficient was entered as 0.6. The Orifice Coefficient was kept at its
default value of 0.8 for fully submerged flow conditions. Additionally, the head reference for the gate equation
defaults to the "Sill (invert)" of the gate opening. This is appropriate for gates that release water to a channel, rather
than to the free atmosphere.
For the Weir Flow data, the Weir Shape was selected as "Broad Crested". The weir coefficient was entered as 2.6. The
OK button was selected at the bottom of the Gate Editor and the gates appeared on the Inline Structure Data
Editor . (Note: The ineffective flow areas will be entered subsequently.) At this point, the user should zoom in on the
gate openings to ensure that they do not overlap and appear as intended. The Inline Structure Data Editor was
then closed.
To enter data for the structure, the Weir/Embankment icon was selected from the Lateral Structure Data Editor.
The data was entered as shown in the figure below on the Lateral Weir Embankment Window. The 500 foot weir is
broad crested with a width of 15 feet, placed one foot from the upstream river station 900. "Water surface" was
chosen as the weir flow reference method. Returning to the Lateral Structure Data Editor, the HW Position was
selected as "Right overbank" from the pull-down menu and the Tailwater Connection Type was chosen as Storage
Area. Lastly, a connection was made to the storage area by selecting the Set SA button. A window then popped up,
allowing the user to select a storage area to connect to. On the pull-down menu "Storage area: SA 1" was selected.
The gate type was selected as "Overflow (open air)". The Height, Width, and Invert of the gate opening were entered
as 7.5, 10, and 91 feet, respectively. On the lower left side of the editor, the centerline station for the gate was
entered as 250. The Weir Shape was chosen as "Sharp Crested" from the dropdown menu and the "User entered
coefficient" was selected as the weir method. The coefficient was entered as 3.2.
the beginning of the simulation time window. The hydrograph was then manually entered into the table.
The elevation controlled gate used for the storage area connector was entered in the Elevation Controlled Gates
editor (below) that appeared once Elev Controlled Gates was selected as the boundary condition type. The gates
operation is "Based on difference in stage," which was chosen from the dropdown of the Reference field. Once
selected, the rest of the fields appeared and were filled in as shown below. The First Reference and Second
Reference were selected as "Storage area: SA 1" and "Storage area: SA 2," respectively, by clicking the Set SA button
under each field and choosing the appropriate storage area from the dropdown. The Stage difference at which gate
begins to open/close was set to 2/0. In laymen's terms this says that when the water surface in SA 1 is two feet above
that of SA 2 the gate begins to open; when the water surface of the two storage areas are equal, the gate begins to
close back.
As shown above, the Enter Initial flow distribution radio button was selected and an initial flow of 100 was
entered for the most upstream river station (1000). The initial elevation of both storage areas was set to 90 feet. The
data was then saved (File -> Save Unsteady Flow Data As) in the Unsteady Flow Data Editor and the window was
closed.
24.6 Summary
This unsteady flow example consisted of a relatively complex system which included a storage area connection
(connecting two storage areas), lateral weir, and inline structure – all utilized the overflow gate feature. While all of
the gates modeled were overflow gates, one of the inline structure's gates was closed top while the remaining gates
were open air. The modeler was instructed on how to setup the gates and where to review the output results.
By reviewing the water surface profiles and the output tables, the user can determine the type of flow through the
gate openings and determine if adjustments to the gate settings are required to provide for a selected water surface
elevation. The Detailed Output table provides detailed output for each gate group, for any profile. The Profile
Summary Output Table provides upstream energy and water surface elevations along with the total weir and gate
flow.
a minimum elevation of 70 feet was entered for Upper Pond, while an area of 60 acres with a minimum elevation of
68 feet was entered for Lower Pond. This method computes the product of a constant area and change in water
surface elevation to calculate storage volume.
a floodwave that peaked at 3500 cfs. A flow hydrograph was entered at River Station 0.4 of Butte Creek in the same
manner.
RS locations were populated in the Boundary Conditions Table. The corresponding cell next to the river station
number was highlighted under the Boundary Condition column and the Groundwater Interflow button was
selected.
Upon selecting Groundwater Interflow as the boundary condition a window opens up titled "Groundwater
Interflow" (see image below). At the top of the window, the user selects the river stations that the interflow will be
between. The upstream river station is the one picked from the Unsteady Flow Data Editor earlier, while the
downstream value is selected from a dropdown menu. For this example, the river stations corresponding to RS 10.3
and RS 9.7 were RS 9.9 and RS 9.4, respectively. Next, the radio button Enter Table was selected. The Data time
interval dropdown was kept at the default of 1 Hour. Under Select/Enter the Data's Starting Time Reference the Use
Simulation Time radio button was selected, automatically filling in the Date and Time fields. The groundwater
The unsteady flow simulation is what actually creates the stage and flow hydrographs. The program is a three step
process that requires the tables from the geometry pre-processer, the boundary conditions, and flow data from the
interface to operate.
The post processor is used to compute detailed hydraulic information for a set of user specified time lines during
the unsteady flow simulation period. By running the post processor, the user will have all of the available plots and
tables for unsteady flow that HEC-RAS normally produces for steady flow. The Geometry Preprocessor, Unsteady
Flow Simulation, and Post Processor were all checked on. More information can be found in Example 17 Unsteady
Flow Application.
desired to compute flow distribution information for all the cross sections. A flow distribution was calculated at Fall
River Lower Reach RS 9.6 by selecting these parameters from the corresponding fields under the Set Specific
Locations Subsection Distribution. Next, the left overbank (LOB dropdown) was divided into two subsections, the
channel into 10 subsections, and the right overbank (ROB dropdown) into two subsections. Once all the parameters
were entered, the Set Selected Range was clicked, which populated the blank field below. OK was selected and all
25.5 Summary
This unsteady flow example applies the groundwater interflow boundary condition to multiple reach and storage
area locations. The program uses Darcy's equation to calculate the interchange between the surface waters and the
groundwater. The user is shown how to input all necessary parameters and where to find the appropriate output
data.
The schematic gives a visual representation of the geometric elements entered by the user.
been entered and applied. Unsteady flow models, such as this one, require that
the Add a Preissmann Slot be checked on. The Preissman Slot option can be turned on from each lidded cross
section's editor (see figure above) or from a table that shows the locations of all lidded cross sections. This
comprehensive table can be viewed by selecting Priessmann's Slot on Lidded XS's from the Tables menu on the
Geometric Data Editor. The table is show in the figure below. For more information on Preissmann Slot Theory, the
user is encouraged to review Chapter 16 of the User's Manual and Chapter 4 of the Hydraulic Reference Manual.
bend's loss value was distributed amongst the bend's surrounding cross sections to represent the fact that bend
losses typically occur over a range of pipe and not at one point.
dragging the mouse for each desired point for the storage area. The final point for the storage area was selected by
double clicking the mouse. Upon double clicking the mouse, a dialog box appears requesting a storage area name.
A name was entered for each storage area created. Points on a storage area may be added or removed by using the
Add Points to a Reach or SA and Remove Points to a Reach or SA options found in the Edit menu.
The storage areas that the connector transferred water to and from were selected by pressing the Set SA/2D
buttons next to the From and To location fields. The Structure Type was chosen as Weir from the drop down menu
for all of the storage area connections. Next the Weir/Embankment icon was selected. The weir width, coefficient,
and weir shape were all entered. In addition the weir station and elevations were entered into the table. Both the
Connection Data Editor and the Storage Area Connection Weir Data can be seen below.
Claiborne's Reach 2 has culverts along the structure (as do most of the other lateral structures in the system). To
enter the data for the culverts, the Culvert icon was selected from the Lateral Structure Data Editor (see second
most previous figure). This activated the Culvert Editor and the data was entered as shown in the figure below.
Three five by five box culverts were entered at centerline stations 200, 300, and 1000. The culverts were identical,
each having a length of 200 feet and flared wingwalls. No flap gates were used at this lateral structure. For more
information on entering culverts the user is encourage to review examples 3 or 4 and Chapter 6 of the User's
Manual.
The Breach Progression tab was chosen next and a one-to-one linear relationship was established between the
breach time and the breach completion by clicking Set to Linear. Alternatively, Sine Wave could have been
selected, giving the breach completion/time relationship a sine correlation.
Lastly, the optional Breach Repair tab was selected from the Levee (Lateral Structure) Breach Data Editor. This
option allows the user to signify that the levee will be repaired during the course of the simulation, thereby altering
the lateral structure's geometry and unsteady flow calculations. The three fields that must be filled in by the user if
selected are: Number of hours after full breach to start repair, Total repair time (hours), and Final filled in elevation.
For this example, the fields' values were 48 hours, 24 hours, and 10 feet, respectively.
The top row of the editor consists of three data entry fields.
For the first field, the user must enter the Distance from the upstream cross section (RS 21200) to the upstream side
of the weir. For this example, this distance was 20 feet. Next, the Width of the weir was entered as 40 feet. This is a
total distance of 60 feet. The distance from river station 21200 to 21120 is 80 feet. Therefore, the distance from the
downstream end of the weir to cross-section 21120 is 20 feet. The last field in the top row of the editor is the Weir
Coefficient and was entered as 2.6, coinciding with the Broad Crested weir shape chosen at the bottom of the editor.
The central portion of the editor consists of a table in which the user must enter the station and elevation data for
the weir. The first station of the weir was entered as zero and the last station is at 50 feet – both have an elevation of
20 feet. With these weir stations and elevation, the program will block out the entire area below the weir crest.
The inline structure has one gate. To enter the data for the gates, the Gate icon was selected from the Inline
Structure Data Editor (shown below) which activated the Gate Editor. The Gate type (or methodology) was left at
the default of "Sluice" from the pull down menu. The Orifice Coefficient, Weir Shape, Head Reference, and Weir
Coefficient were also kept at their default values (0.8, Sill (Invert), Broad Crested, and 3, respectively). The Sluice
Discharge Coefficient was set to 0.6. The gate was set with a centerline station of 25 and an invert of -4 feet. It was 2
feet in height and width. The OK button was selected at the bottom of the Gate Editor and the gates appeared on
the Inline Structure Data Editor as shown in the figure below. The Inline Structure Data Editor was then closed.
section corresponding to the row in which you cursor lies in the table.
radio button Use Simulation Time was selected causing the initial propagation of the hydrograph to coincide with
the beginning of the simulation time window. The hydrograph was then manually entered into the table.
Levee (Lateral Structure) Breach from the Unsteady Flow Analysis window.
theory, which is the method that allows pressurized pipes to be modeled with unsteady flow in HEC-RAS.
was activated. The rule caused all subsequent flow rates from the pump station to hold steady at 100 cfs.
26.7 Summary
Urban modeling is unique in that it often incorporates several HEC-RAS elements that many river systems don't
utilize. This example gave note to many of those features - pressurized pipe flow, how to model pipe connections,
drop inlets, minor losses, levee breaches, interconnected storage areas (surface overland flow system) and pump
stations – by modeling a portion of a New Orleans Parish system. For many of these features, the user was shown
how to access them from the Geometric Data Editor. For instance, the user was shown how to turn a typical cross
section into a pipe by using the Lid feature which is accessed through the Cross Section Editor. After entering all of
the geometric data, the user was shown how to enter the unsteady flow data, as well as the plan information. Once
the model was run, some of the relevant output data is reviewed.
To enter weir data for the structure, the Weir/Embankment icon was selected from the Lateral Structure Data
Editor. The data were entered as shown in the figure below on the Lateral Weir Embankment Window. The editor
requires an entry for the following fields: Weir Width, Weir Computations, Weir flow reference, Weir Coefficient, Weir
Crest Shape, HW Distance to Upstream XS, and TW flow goes. Further description of these fields can be found in
Chapter 6 of the User's Manual.
Returning to the Lateral Structure Data Editor, the HW Position was selected as either "Left overbank" or "Right
overbank" from the pull-down menu depending on where the lateral structure resided. The Tailwater Connection
Type was chosen as Storage Area, meaning that the lateral structure connected the river to a storage area, or Out of
System, which means that the water leaving the river via the lateral structure is no longer modeled upon leaving.
For the Storage Area type tailwater connections, an individual storage area was selected to correspond to the
lateral structure by clicking the Set SA/2DFA button. A window then popped up allowing the user to select a storage
area to connect to from a pull-down menu.
Next the weir data was entered. To enter the data for the weir, the Weir/Embankment icon was selected from the
left side of the Inline Structure Data Editor. This activated the Inline Structure Weir Station Elevation Data
Editor as shown in the figure below. This editor is similar to the deck/roadway editor used for bridges and culverts.
The top row of the editor consists of three data entry fields. For the first field, the user must enter the Distance from
the upstream cross section to the upstream side of the weir. Next, the Width of the weir was entered (both inline
weirs has a width of 30 feet). The last field in the top row of the editor is the Weir Coefficient and was entered as 2.6
for both inline weirs, coinciding with the Broad Crested weir shape chosen at the bottom of the editor.
The central portion of the editor consists of a table in which the user must enter the station and elevation data for
the weir. With these weir stations and elevation, the program will block out the entire area below the weir crest.
To enter the data for the gates, the Gate icon was selected from the Inline Structure Data Editor (above). This
activated the Gate Editor as shown in the figure below. Each inline structure had four gates, with each gate
belonging to its own gate group. A separate gate group for each gate was created by clicking the New Gate Group
button (looks like a blank sheet of paper at the top of the editor). The Gate type (or methodology) was changed to
Overflow (open air) from the pull down menu for all gates. The Weir Shape was selected as Broad Crested Weir and
the Weir Coefficient was entered as 3. Each gate's geometric properties were entered on the left side of the editor
(see figure below). The OK button was selected at the bottom of the Gate Editor and the gates appeared on the
Inline Structure Data Editor as shown in the figure below. The Inline Structure Data Editor was then closed.
The flow hydrograph associated with the upper reach of the Mississippi River (RS 110.4) differed from the Ohio River
upstream boundary condition in that it pulled its flow information from a DSS file. A DSS file was connected to the
RS by pressing the button Select DSS file and Path. When this button was pressed a DSS file and pathname
selection screen appeared as shown in the figure below. The desired DSS file was found and selected by using the
browser button at the top. Once the DSS file was selected, a list of all of the DSS pathnames within that file was
displayed in the table. The pathname filters were used to reduce the number of pathnames shown in the table.
When the desired DSS pathname was found the window was closed and the filename and pathname will be
recorded in the Flow Hydrograph Window.
Lateral inflow hydrographs were added throughout the river reaches as a means of representing flow from an un-
modeled tributary. This is often done when there is un-modeled tributary flow that has a noticeable impact on the
modeled system but the impact of the modeled floodwave on the tributary is not considered to be within the scope
of the model's objective. Therefore, the effects of the flow are captured without having to model unnecessary parts
of the system. For all five of the lateral inflow hydrographs used in this example, the hydrograph was read in from
DSS. This was done in the same manner as was the flow hydrographs for the upstream BC (explained above).
As seen in the figure above, the Enter Initial flow distribution radio button was selected and initial flows were
entered at the upstream river station of each of the three reaches. The initial water surface elevation of all twenty-
seven storage areas was set as well. The data was then saved (File -> Save Unsteady Flow Data As) in the
Unsteady Flow Data Editor and the window was closed.
window). When the DSS file was connected, the table below was automatically populated with the DSS file
pathnames, which contained the observed values.
Finally, the appropriate DSS file and pathname was connected to each individual RS in the top table by highlighting
the desired row and then selecting the corresponding gage data (pathname) from the bottom table. The
information from the pathname selected in the bottom table should then be populated in the highlighted row in the
top table. In this example, Part A, Part B, and Part C columns of the DSS file tell the user which gage location and
data type (stage, flow, etc.) the user is selecting. It is important that the user is careful to match the DSS file with the
correct location and required data type. However, matching column Part D, which has the time that the data
begins, is not pertinent. It is not required that the date listed in Part D for the row that the user selected fall within
the simulation time of the model. As long as the gage has collected data during the simulation time frame, HEC-RAS
will automatically find the corresponding gage times. The pathname data can be visually inspected by clicking the
Plot Selected Pathname button in the lower left hand corner of the window.
If it were the case that the stage values were consistently to low for RS 43.7 of the Upper Mississippi River, the
channel (highlighted in green) and overbank Manning's n values near RS 43.7 could be selected and adjusted
upward by a factor of 1.05 by clicking the Multiply Factor from the Selected Area Edit Options above the table. This
would increase all of the highlighted Manning's n values by 5%. This change would increase the stage near and
upstream of RS 43.7 when the plan was reprocessed. Alternatively, the Add Constant or Set Values could have been
selected to change the Manning's n values.
As mentioned earlier, the automated calibration procedure computes adjustments to existing flow roughness
factors to follow observed stages. An initial set of flow roughness factors must already be set in the geometry or
plan and the observed stages must already be in the unsteady flow data before implementing this calibration
procedure. In this example, the initial flow vs. roughness factors are all being set to a value of 1.0 (no change to the
base n values), since the automated calibration process will be computing flow vs roughness factors for us. After the
initial set of flow versus roughness factors were set up for all of the reaches to be calibrated, the Automated
Roughness Calibration tool was accessed from the Options menu of the Unsteady Flow Analysis Editor. The
Unsteady Roughness Calibration Editor appears as in the figure below. There are four distinct regions within the
editor: Calibration Parameters, Calibration Regions, Forcing Internal Observed Flows, and Forcing Internal Observed
Stages. The last two are not thoroughly discussed in this example because they were not used. The user is
encouraged to look to Chapter 16 of the User's Manual for more information on the Unsteady Roughness Calibration
Editor.
First, the fields within the Calibration Parameters region were entered. The Optimization Method was set to Global.
The other possible selection from the dropdown was the sequential method. The sequential method can potentially
lead to better results; however, it takes much longer to run because it goes through many more iterations.
Computationally, the primary distinction between the two methods is that the global method calibrates all
Manning's n values simultaneously, whereas the sequential method calibrates the upstream river stations first and
then works its way down to the lower river stations. Also, for sequential method calibration, the bottom two regions
– Forcing Internal Observed Flows and Forcing Internal Observed Stages – must be used; however, they are not used
with the global method.
Next, the remaining fields within the Calibration Parameters were filled in. The Error Evaluation Method was set to
Average Error for this example. The Maximum number of iterations and Maximum change in factor per iteration were
set to 10 and 0.05, respectively. The iteration tolerance is the acceptable convergence error, which basically tells the
program what amount of error the user is comfortable with and when to stop running. This value was set to 0.1 for
this example. In layman's terms - for the entries of this example - these parameters tell the program that it has up to
ten attempts to get an error under 0.1 while adjusting the factors by no more than 0.05 (5%) on any iteration and
keeping its estimates within a range of 0.5 and 1.7. The limit placed on the amount changed per iteration is used to
keep the program from overshooting the best answer, while the range ensures that the final answer isn't something
unreasonable. An optional time window was also used for this example. The start date was five days after the
simulation time and ended at the same time as the simulation (25 Sep 2004 – 05 Mar 2005). This function can be
used if the user doesn't want to evaluate a certain portion of the flow hydrograph.
The Calibration Regions portion of the window was completed next. This area was already populated with the
roughness factor river segments that were created earlier. Under the Calibration Gage column, a gage was paired
with the river segment that it corresponded to by selecting the cell and then choosing the appropriate gage from a
drop down menu that appeared. Lastly, each box under the Calibrate column heading was checked on.
The OK button was clicked and the Unsteady Flow Analysis was run. After the model ran, the new calibrated
Manning's n values can be found by reopening the Unsteady Roughness Calibration Editor and clicking the Plot
Flow-Roughness Calibration Results button in the lower left hand corner of the window. A DSS plot will appear
(see figure below) that displays pre- and post-calibration roughness factors. The first tab plots the data, while the
second tab (as seen in the figure below) displays the results in a table. If the user is satisfied with the calibrated
roughness factors, they can copy and paste them into the Flow Roughness Factor Editor so that the new values will
be used in the next run.
27.6 Summary
This example demonstrated how to set up a gaged unsteady flow model that was calibrated with HEC-RAS's
Automated Roughness Calibration feature. The user was first shown how to set up the geometric data, unsteady flow
data, and the unsteady flow analysis. The primary focus however was demonstrating how to use the Flow
Roughness Factors Editor and the Automated Roughness Calibration feature. Once the model was calibrated, the
results were found in a DSS plot from the Automated Roughness Calibration Editor. The output was then analyzed by
comparing the observed, calibrated, and originally simulated stage hydrographs. The importance of calibrating a
model should be made clear from this example. What also should be conveyed is that the Automated Roughness
Calibration feature makes calibration much quicker and simpler.
28 References
Barkau, Robert L., 1992. UNET, One-Dimensional Unsteady Flow Through a Full Network of Open Channels,
Computer Program, St. Louis, MO.
Federal Highway Administration, 1995. Evaluating Scour at Bridges, Hydraulic Engineering Circular No. 18, U.S.
Department of Transportation, Washington D.C.
Federal Highway Administration, 1990. User's Manual for WSPRO - A computer model for water surface profile
computations, Publication No. FHWA-IP-89-027, 177 p.
Federal Highway Administration, 1985. Hydraulic Design of Highway Culverts, Hydraulic Design Series No. 5, U.S.
Department of Transportation, (September) Washington D.C.
Federal Emergency Management Agency, 1985. "Flood Insurance Study Guidelines and Specifications for Study
Contractors," FEMA 37, Washington, D.C.
Hydrologic Engineering Center, 1995. Flow Transitions in Bridge Backwater Analysis, RD-42, U.S. Army Corps of
Engineers, Davis, CA.
Hydrologic Engineering Center, 1996. HEC-RAS River Analysis System, Hydraulic Reference Manual, U.S. Army Corps
of Engineers, Davis, CA.
Hydrologic Engineering Center, 1996. HEC-RAS River Analysis System, User's Manual, U.S. Army Corps of Engineers,
Davis, CA.
Hydrologic Engineering Center, 1995. UNET, A One-Dimensional Unsteady Flow Through a Full Network of Open
Channels, User's Manual, U.S. Army Corps of Engineers, Davis, CA.
King, H.W. and E.F. Brater, 1963. Handbook of Hydraulics, Fifth Edition, McGraw Hill Book Company, New York.
Shearman, J.O., 1990. User's Manual for WSPRO - A computer model for water surface profile computations, Federal
Highway Administration, Publication No. FHWA-IP-89-027, 177 p.
United States Geological Survey, 1979. Backwater at Bridges and Densely Wooded Flood Plains, Beaver Creek Near
Kentwood, Louisiana, Hydrologic Atlas No. HA-601.
Yarnell, D.L., 1934. Bridge Piers as Channel Obstructions, Technical Bulletin 442, U.S. Department of Agriculture,
Washington, D.C.
References – 378