Professional Documents
Culture Documents
FishXing - Gradually Varied Flow
FishXing - Gradually Varied Flow
Overview
The hydraulic calculations from the culvert outlet to the culvert headwater immediately
upstream of the inlet are based on the conservation of energy and mass. This is described in
the basic energy balance equation:
Where:
y = Depth (ft, m)
V = Velocity (ft/s, m/s)
Z = Change in elevation (ft,m)
TW = Tailwater
HW = Headwater
g = Accceleration due to gravity (ft/s2)
In most cases the approach velocity (VHW) is low and the approach velocity head is neglected.
Similarly the exit velocity can be neglected in the energy equation if the upstream and
downstream channels are similar, reducing the energy equation to:
FishXing assumes the headwater (approach) velocity and tailwater velocity are zero and does
not include bend losses, junction losses, or grate losses. For determining water depths within
the culvert, FishXing solves the appropriate form of the energy equation using a step method.
Solutions are obtained at regular intervals throughout the culvert.
Below is an outline of the generalized procedure used in FishXing for determining the water
surface profile and water velocities within a culvert at a specific flow:
1.
2.
Determine Boundary Conditions for solving the Gradually Varied Flow equations.
o
3.
Determine starting location and depth for the downstream and/or upstream
boundary.
Solve the Gradually Varied Flow or Full Flow equations to obtain a water surface
profile through the culvert.
o
For steep slopes, check for a Hydraulic Jump within the culvert.
4.
5.
Calculate average cross sectional water velocities within the culvert and
the Contraction Velocity within the Inlet Zone.
6.
http://www.fsl.orst.edu/geowater/FX3/help/4_Calculations/Culvert_Hydraulic_Calculations_Overvie
w.htm
Horizontal (H) if So = 0
Critical (C) if Yo = yc
For culverts with slopes > 0, the slope classifications change as flows change. It is not
uncommon for a culvert to switch from being Mild Slope to Steep Slope as flows increase.
Type 2 curve: Depth is between ycand yn, flow can be either subcritical or
supercritical
Type 3 curve: Depth is less than both ycand yn, flow is supercritical.
Since normal depth (yn) is undefined for Horizontal and Adverse slopes, they only experience
Type 2 and Type 3 curves.
Note: Since FishXing does not model slope breaks within a culvert, therefore type 3 curves
do not occur within the model.
http://www.fsl.orst.edu/geowater/FX3/help/4_Calculations/Classification_of_Water_Surface_Profiles.
htm
To develop a water surface profile (WSP) through the culvert using Gradually Varied
Flow calculations one or more boundary condition must be defined. A boundary condition is a
section of the channel where the depth of flow is known at a given flow rate. For culverts,
these boundary conditions occur at or near the inlet and outlet.
For backwater calculations (WSP calculated from downstream to upstream) a downstream
boundary condition is necessary. When performing frontwater calculations WSP calculated
going downstream) an upstream boundary condition is required. Determining the type of
boundary condition and calculation (backwater verses frontwater) needed is a function of
the hydraulic
slope and tailwater depth. In FishXing the tailwater depth is defined by one
Steep Slopes
GVF Boundary Conditions = Critical Depth at inlet.
Frontwater GVF Calculations begin at critical depth at inlet and end at outlet
When the Tailwater Depth (yTW) is greater than Critical Depth (yc ) the culvert is controlled by
the downstream water surface. Backwater GVF calculations begin at the outlet and proceed
upstream to the inlet.
When Tailwater Depth (yTW) is less than the Free Surface Depth (yfs) the flow is passing
through Critical Depth (yc ) and entering a zone of Rapidly Varied Flow. RVF is approximated
linearly by drawing the flow profile from critical depth at a distance of 4y c from the outlet
through yfs at the outlet. When yTW is negative, an outlet plunge exists.
When Tailwater Depth (yTW) is less than the Critical Depth (yc ) and greater than Free Surface
Depth (yfs) the flow is passing through critical and entering a zone of Rapidly Varied Flow. The
water surface within the RVF zone is approximated by drawing the flow profile from the
tailwater elevation at the outlet to critical depth at a distance of 4y c upstream from the outlet.
Steep Slopes
For Steep slope culverts, a downstream boundary condition is only needed if a backwater
calculation is required.
In this case the Tailwater Depth (yTW) is below critical depth, so no backwater calculation is
required.
For hydraulically Steep slopes, if the Tailwater Depth (yTW) is greater than the Critical Depth
(yc ) then the downstream boundary condition is yTW and a backwater calculation is performed.
In addition, for all Steep sloping culverts a frontwater calculation is also performed. Determine
the extent the flow profile is influenced by the tailwater depth requires comparing the
frontwater and backwater generated profiles and identifying the location of a
hydraulic
See: Outlet
http://www.fsl.orst.edu/geowater/FX3/help/4_Calculations/Boundary_Conditions.htm
2.
One Dimensional (can only calculate average cross sectional water velocity)
Where:
So = Bottom slope, positive in the downward direction
Sf = Friction slope, positive in the downward direction
y = Water depth, measured from culvert bottom to water surface
x = Longitudinal distance, measured along the culvert bottom
Fr = Froude number
The Friction slope is approximated from Mannings Equation:
Where:
Sf = Friction slope, positive in the downward direction
n = Mannings roughness coefficient
V = Average cross section velocity
= Constant equal to 1.49 for English units and 1.00 for SI units.
R = Hydraulic radius, (Wetted Area / Wetted Perimeter)
Numerical Integration of the GVF Equation:
The GVF differential equation does not have an analytical solution. Therefore, FishXing uses
numerical integration to generate a water surface profile. Numerical integration is a technique
of dividing the channel, or culvert, into numerous short reaches and then performing the
computations from one end of the reach to the other.
FishXing primarily uses the Standard Step Method of numerical integration. The following
form of the equation is used:
Where:
E = Change in specific energy from one end of the reach to the other
Sfave = Average friction slope across the reach
x = Longitudinal distance from one end of the reach to the other
y = Depth of water
Q = Flow rate
g = Gravitational acceleration
A = Wetted cross sectional area
Since the friction slope and wetted area are functions of depth, solving for depth at a given
distance (x) requires an iterative solution. FishXing uses a bisection method to find the
solution.
Other Considerations
If only a portion of the culvert becomes pressurized, FishXing will switch between the GVF
and full flow equations.
http://www.fsl.orst.edu/geowater/FX3/help/4_Calculations/Gradually_Varied_Flow.htm
Full Flow
The velocity within a full flowing culvert is calculated as:
Where:
V = Average velocity in the culvert barrel [ft/s or m/s]
Q = Flow rate through culvert [ft3/s or m3/2]
A = Full cross sectional area of the flow [ft2 or m2]
When flowing full, the friction loss within the culvert barrel, Hf, over a given length of culvert, L,
is calculated based on the friction slope, Sf:
To calculate the friction slope, Mannings equation is rearranged into the following form:
Where:
If the hydraulic grade line intersects the top of the culvert, FishXing switches between full flow
and open channel flow calculations.
http://www.fsl.orst.edu/geowater/FX3/help/4_Calculations/_Full_Flow.htm
http://www.fsl.orst.edu/geowater/FX3/help/4_Calculations/_Rapidly_Varied_Flow.htm
Hydraulic Jumps
Whenever the flow profile changes from supercritical to subcritical, hydraulic jumps will occur.
A hydraulic jump represents a significant head loss that manifests in available energy for scour
and creation of turbulence. Hydraulic jumps are one of the three occurrences of Rapidly Varied
Flow that FishXing approximates. Hydraulic jumps are generally an undesirable condition for
fish passage and erosion control.
1.
Compute the upstream supercritical water surface profile by solving the Gradually
Varied Flow equations from the inlet depth equal to critical depth. Calculations
proceeding in the downstream direction are called "frontwater calculations".
2.
Starting at the downstream boundary condition at the outlet, compute the subcritical
water surface profile in the upstream direction. Calculations proceeding in the
upstream direction are called "backwater calculations".
3.
At each node, compute the momentum (specific force) for associated with the two
depths (supercritical and subcritical).
Where:
M = Momentum or specific force
Q = Flow rate
g =Acceleration due to gravity
A = Cross sectional area
zbar = is the distance from the water surface to the centroid of the cross sectional
area of flow.
4. Beginning at the outlet and proceeding towards the inlet, compare the momentum
associated with the corresponding supercritical and subcritical depths. When the
momentum associated with the subcritical profile becomes less than the momentum
associated with the supercritical profile, a hydraulic jump is assumed to have
occurred between the two nodes.
http://www.fsl.orst.edu/geowater/FX3/help/4_Calculations/_Hydraulic_Jumps.htm
Inlet Zone The region within the culvert inlet where entering flows contract and
then expand rapidly. FishXing reports the inlet zone as the first two nodes (points)
down the culvert. For culverts with a width or diameter of 9 feet or less the inlet zone
is 2 feet in length, while for culverts larger than 9 feet the inlet zone is 3 feet in
length (Belhke, 1992). At this location within the culvert, FishXing calculates
the Inlet Contraction Velocity.
Barrel - All the nodes between the first two nodes and last two nodes in the culvert
Outlet Zone - The distance from the outlet to a distance equal to 4yc (four times
critical depth) within the culvert. This zone consists of the last two nodes of the
culvert. The area of rapidly varying flow associated with a full or partially free
surface outlet will be contained within the Outlet Zone.
The zones also define where the inlet, barrel, and outlet Velocity Reduction Factors are
applied.
For this application, the Barrel Velocity (VB) is the average cross-sectional flow in the barrel
approximately one culvert diameter downstream from the inlet and represents the area
immediately downstream of the inlet Contraction and Expansion zone.
See also: Inlet Head Loss Coefficient, Inlet Contraction Velocity, Velocity Reduction Factors
http://www.fsl.orst.edu/geowater/FX3/help/4_Calculations/_Inlet,_Barrel_and_Outlet_Zones.htm
Headwater Calculations
The depth of water above the culvert inlet bottom is known as the Headwater Depth. This
depth represents the amount of energy available to convey water through the culvert.
Headwater depths are a function of the entrance shape, along with the depth and velocity
immediately inside the culvert.
Headwater depths are determined by summing the energy losses associated with entrance
shape, exit expansion and friction of the culvert. This is described in the basic energy balance
equation:
Where:
Y = Depth (ft, m)
V = Velocity (ft/s, m/s)
Z = Change in elevation (ft,m)
TW = Tailwater
HW = Headwater
g = Accceleration due to gravity (ft/s2)
In most cases the approach velocity (V) is low and the approach velocity head is neglected.
Similarly the exit velocity can be neglected in the energy equation if the upstream and
downstream channels are similar, reducing the Headwater calculation to:
where yHW is the sum of all losses and represents the difference in water surface elevation at
the Outlet (headwater) and outlet (tailwater) FishXing does not include bend losses, junction
losses, or grate losses.
Entrance Loss depends on the geometry of the inlet. This loss is expressed as the velocity
head immediately inside the culvert reduced by the entrance loss coefficient, Ke.
Where:
HL = Head Loss (ft)
Ke = Entrance Loss Coefficient
V = Velocity in the culvert barrel (ft/s)
g = Acceleration due to gravity
FishXing calculates the headwater depth using the following Equation:
FishXing calculates total headloss as the sum of the Entrance Loss, Exit Loss, and Friction
Loss.
Note: High exit velocities can cause downstream scour and present a barrier to fish passage
so Exit Velocity should be considered in the design of culverts.
http://www.fsl.orst.edu/geowater/FX3/help/4_Calculations/Headwater_Calculations.htm
Illustration of inlet zone hydraulics, adapted from Behlke et, al. 1992.
The hydraulic environment within the inlet zone of a culvert is typically classified as rapidly
varying flow. Two distinct sections exist within the inlet zone: a contraction zone where
streamlines contract as the water enters the culvert and an expansion zone that exists
immediately downstream of the contraction. This contraction occurs as water from the
upstream channel enters into the culvert inlet. If the culvert is narrower than the channel the
abrupt change in width causes flows to be constricted.
Head loss occurs as water passes through the inlet zone. The magnitude of the head-loss is
dependent on the velocity within the culvert barrel and inlet geometry, which determines the
appropriate entrance loss coefficient. The head-loss coefficient is a function of transition
efficiency as the water enters the culvert from the upstream channel.
As the water enters the culvert, velocities increase until they reach the point of maximum
contraction. In this zone there is a drop in the hydraulic grade line, which is evident by a drop
in water surface. At the same time there is an increase in kinetic energy that is equal to the
drop in the hydraulic grade line minus frictional losses. Just downstream of the contraction is
an area of rapidly expanding flow. A large portion of the entrance headloss occurs within the
expansion
The contraction zone at a culvert inlet can hinder fish passage due to the increased velocities
and steeper water slope. FishXing approximates the maximum contraction velocity by
assuming the calculated entrance headloss through the inlet zone is first converted entirely to
kinetic energy, and that this gained kinetic energy is then lost in the expansion.
FishXing uses the following equation to calculate the contraction velocity, and reports it as the
velocity immediately inside the culvert in the Water Surface Profile Results:
Where:
Vcntr = Contraction velocity
Ke = Entrance loss coefficient
VB = Average cross sectional water velocity within the culvert barrel, located at the
point where the velocity contraction occurs.
For this application the Barrel Velocity (VB) is the average cross-sectional velocity at the
location of the contraction, calculated assuming gradually varied flow conditions. For
calculations, the location of the contraction is defined as 2 feet (0.61 m) from the inlet for
culverts 9 feet wide (2.75 m) or smaller and 3 feet (0.91 m) from the inlet for culverts larger
than 9 feet (2.75 m) (Behlke, 1992). FishXing assumes the approach velocity is zero just
upstream of the inlet.
See also:
Inlet Head Loss Coefficient, Continuity Equation
http://www.fsl.orst.edu/geowater/FX3/help/4_Calculations/_Inlet_Contraction_Velocity.htm
Tailwater Calculations
FishXing allows the user to define the tailwater elevation using one of three methods:
1.
2.
User Defined Tailwater Rating Curve: Inputting a rating table that relates tailwater
elevation to flow.
3.
FishXing calculates the tailwater elevation - flow relationship for the cross section assuming
uniform flow, and using Mannings Equation. Calculations begin with the tailwater set at the
lowest elevation in the cross section and a corresponding flow equal to zero. Flow through the
tailwater cross section is then calculated at regular intervals in tailwater elevation. Calculations
end once the water surface encounters one of the cross section ends (whichever side is
lower).
Input a rating curve that includes the entire range of flows you plan to model (User
Defined Tailwater Rating Curve Method)
2.
Enter a channel cross section that will contain all of the flows you plan to model
(Channel Cross Section Method)
http://www.fsl.orst.edu/geowater/FX3/help/4_Calculations/_Tailwater_Calculations.htm
bottom culverts that are embedded. In these situations, FishXing assumes the bottom
segment to be horizontal.
For culverts, FishXing uses the following form of the Horton-Einstein equation:
Where:
ncomposite = Mannings roughness coefficient for multiple materials
Pside = Perimeter of side material
nculvert = Mannings roughness coefficient of culvert material
Pbottom = Perimeter of bottom material
nbottom = Mannings roughness coefficient of bottom material
Where:
ncomposite = Mannings roughness coefficient for multiple materials
P = Perimeter of material
i = subsection of crossing
If a Darcy-Weisbach friction factor or Chezy coefficient is used instead of Mannings n,
FishXing converts the alternative roughness coefficient into an equivalent Mannings n to
calculate a composite roughness.
http://www.fsl.orst.edu/geowater/FX3/help/4_Calculations/_Composite_Roughness_Calculations.htm
Plunging Flow
Water leaving a perched culvert has a velocity with both horizontal and vertical components.
By neglecting all frictional losses, a simplified projectile equation can be applied to the exiting
water. Using this equation, FishXing calculates the path of the falling water from exiting the
culvert outlet until it impacts the downstream outlet pool. FishXing assumes the location in
which the plunging flow impacts the outlet pool is where the fish will attempt to leap from.
The height the water plunges, H, and the horizontal distance the water travels from the outlet,
L, is described by the following standard projectile equations:
Where:
H = Vertical plunge distance.
L = Horizontal plunge distance
Vout = Exiting water velocity at the culvert outlet [ft/s or m/s]
g = Gravitational acceleration [32.2 ft/s2 or 9.8 m/s2]
t = Time for the exiting water to fall from the culvert outlet to the pool [s],
= Angle at which the water exits the culvert outlet [degrees or radians]
Since H is known, FishXing solves for the time, t, and the horizontal plunge distance, L.
Exit Angle
The exit angle is the angle the exiting water surface makes with the horizontal. It is calculated
using the last two nodes within the culvert, which is defined as the Outlet Zone. The water
surface slope at the outlet is determined by the applicable Free Surface Outlet conditions.
See Also: Defining Normal Depth, Froude Number and Flow States, Open Channel
Flow, Perched Outlet
http://www.fsl.orst.edu/geowater/FX3/help/4_Calculations/Plunging_Flow.htm
Steep Culverts
At hydraulically steep slopes (yn < yc), if the tailwater is below critical depth (y tw < yc)
supercritical flow will occur throughout the culvert. In this case, the controlling water depth is
at the upstream end of the culvert. If the tailwater depth is below the culvert outlet bottom
(ytw < 0) the culvert is considered perched. Under supercritical flow conditions, FishXing
assumes the water depth at the perched culvert outlet is not influenced by the free overfall.
Non-steep Culverts
For culverts with a hydraulic slope that is not steep (instead having either a mild, critical,
horizontal, or adverse slope) and having a perched culvert outlet (y tw < 0) , the water surface
profile will be on a draw-down curve as it approaches the free overfall. From where the profile
crosses through critical depth to the free outfall the flow enters into a rapidly varied flow
regime. Rouse (1936) examined the case of an overfall at the end of a mild slope and found
that the true critical section is located at the overfall crest. The depth of this crest is a function
of the wetted area at the computed critical depth, yc:
The exact location of the calculated critical depth is indeterminate, as it will move upstream
with increasing flow and downstream with increasing boundary roughness and surface slope.
For relatively smooth rectangular channels, Rouse found this distance approximately equal to
four times the critical depth (4yc) upstream from the crest.
For non-steep culvert with a free overfall FishXing approximates the water surface profile near
the outlet by assuming the water depth is equal to critical depth a distance of 4y c upstream of
the outlet. FishXing then assumes the free surface depth located at the outlet (Y fs) is a
function of 0.71Ac. The water surface is then assumed to be linear between y c and yfs.
http://www.fsl.orst.edu/geowater/FX3/help/4_Calculations/_Free_Surface_Outlet.htm