xpswmm/xpstorm Resource Center xps

Versions Compared

Key

  • This line was added.
  • This line was removed.
  • Formatting was changed.

The EXTRAN Model uses the momentum equation in the links and a special lumped continuity equation for the nodes. Thus, momentum is conserved in the links and continuity in the nodes. The model utilizes the relation Q = A·V frequently in the solution and the derivation of the differential equations used in the solution.

...

The basic unsteady flow continuity equation with lateral inflow and with cross-sectional area (A) and flow (Q) as dependent variables is (Yen, 1986; Lai, 1986):

Equation 1: 

Image Modified 

where, q is the lateral inflow. In EXTRAN this term is zero and the inflows enter the network at the nodes. (The notation for the variables used in this section are listed in the table below:

VariableDescription
As Surface Area of Node
C# Courant Number for the Conduit
D Current Conduit Depth
V Average Conduit Velocity
A Conduit Cross-sectional Area
T Conduit Flow Width
g Gravitational Acceleration
n Manning’s Roughness
R Hydraulic Radius
H Hydraulic Head (z + h)
h Water Depth in Conduit
Q Conduit Flow
Qf Conduit Design Flow
Hf Distance between Junction Invert and Junction Crown
Sc Expansion/Contraction Loss Slope
Se Entrance/Exit loss in Junctions
Sf Friction (Energy) Slope
w Under-relaxation Parameter
L Distance along the Conduit
Dt Time step
z Invert Elevation
y Water Depth in Conduit


The conduit momentum equation may be written in several forms depending on the choice of dependent variables. Using dependent variables flow [Q] and hydraulic head [H] the momentum equation is written (Lai, 1986):

Equation 2: 

Image Modified

Equation 2 is the form of the momentum equation used in EXTRAN. The seven terms in the momentum equation are, respectively:

  1. Local inertia: Image Modified
  2. Convective inertia Image Modified 
  3. Pressure slope Image Modified
  4. Entrance/exit loss (Se
  5. Contraction/expansion loss (Sc)
  6. Friction slope (Sf)
  7. Bed slope (So)

In the EXTRAN Model the bottom slope (So, or Image Modified) is incorporated into the gradient of H, with Image Modified defined as Image Modified. The revised momentum equation using these definitions of bottom slope can be written as:

Equation 3

Image Modified

This momentum equation must be modified to suit the requirements of EXTRAN. The purpose of these modifications to the continuity and momentum equations used by EXTRAN are:

  1. Eliminate the need for the Image Modified term in the continuity equation. EXTRAN uses a mean or center Q in conduits and the term Image Modified is not strictly valid.
  2. Eliminate one equation from the solution by creating a combined continuity-momentum equation.
  3. Linearize Image Modified by expanding the partial differential.

The possible substitutions are numerous and found by substituting A·V for Q or Q/V for A in the Image Modifiedterm in Equation 3:

Equation 4: Image Modified

Equation 5: Image Modified

and linearizing the equations,

Equation 6Image Modified

Equation 7Image Modified

Equation 7 is the conservative form of the convective inertia term. After substituting equation 7 into equation 3, the equivalent form of the momentum equation follows:

Equation 8Image Modified

The continuity equation 1 may be manipulated to replace the third term of equation 8 as follows:

Equation 9Image Modified

Substituting equation 9 into equation 8 to eliminate Image Modified leads to the following equation:

Equation 10Image Modified

The friction slope, Sf, is defined by Manning’s equation and is calculated based on the center hydraulic radius, R and alternatively the upstream and downstream hydraulic radius, Rup and Rdn, respectively, as follows:

Equation 11Image Modified

where k = (n/1.49)2 for U.S. customary units and n2 for metric units. By using the absolute value of the velocity term, Sf is a directional quantity and ensures that the frictional force always opposes the flow.

Equation 12Image Modified

Equation 11 is used whenever Rup and Rdn vary by greater than 300 percent, and equation 12 otherwise. This is a very important term in the momentum equation because the primary balance to the pressure and bed slope is the friction slope. This slope is calculated two ways to better handle the change from a steep water surface slope to a flatter water surface slope.

Incorporating the formulation for Sf (equation 11) into equation 10 yields:

Equation 13Image Modified

This is the form of the momentum equation used by EXTRAN and it has the dependent variables Q, A, V, and H. Equation 13 is still a partial differential equation and must be translated into a finite difference form for its solution using EXTRAN. The finite difference formulation is discussed in the next section.

Additional terms are incorporated into the dynamic wave equation for the entrance/exit loss from junctions and the expansion/contraction loss in conduits. The expansion/contraction conduit slope (Sc) is given by (Fread, 1977):

Equation 14Image Modified

where:

Kec is the expansion contraction coefficient, and the expansion/contraction losses are:

Equation 15

Image Modified

Expansion and contraction losses between conduits in the program are calculated using the expression Image Modified where A is the mean of the current conduit and either the conduit immediately upstream or downstream. The squared velocity difference is the difference between the mean velocities of the three conduits.

The entrance/exit conduit slope is given by:

Equation 16Image Modified

where:

Kee is the entrance-exit loss coefficient, and the expansion/contraction losses are:

Equation 17Image Modified

The entrance and exit head loss is actually calculated in the conduit by adding an additional term to the dynamic flow equation. The velocity at both the upstream and downstream end of the conduit is used to calculate the entrance and exit loss depending on the flow direction.

Equations 14 and 16 are added to the combined continuity-momentum equation to account for these additional losses in head. The losses are also used when the kinematic wave equation is used for supercritical flow. These losses are computed but not listed in the following discussions of the finite difference approximations of the partial differential equations.

...

Stylesheet
stylepadding: 20px; border-bottom: none !important;
classindent-list
Table of Contents
maxLevel2
minLevel2
indent20px
printablefalse

On this section:

Children Display
alltrue
depth1
pageEXTRAN Theory - Dynamic Wave Solution