A flow retardation parameter, f, is used for conduits that have slight or flat slopes. The definition of a slight slope as used in the model is conduits with slopes < 1/10,000. The default value of f (currently 0.05) is used to moderate the solution of the conduit flow.



The above parameter is used as follows in EXTRAN:

Equation 53:  

where:

LHS = dynamic wave equation variables at the new time step or iteration, and

RHS  = dynamic wave equation variables at the old time step or iteration.

Generally, the solution of the flow at the new time step (Qn+1) is solved as follows:

Equation 54: 

This solution is further moderated by using under-relaxation on the new values of RHS and LHS. This control makes the solution smoother and eliminates the growth of potential instabilities.

Equation 55: LHSn+1 = (1-w)·LHSn + w·LHSn+1/2

Equation 56: RHSn+1 = (1-w)·RHSn + w·RHSn+1/2

The combined pressure-bed slope term in the dynamic wave equation is solved at the new and old time step by factoring out a Q from the A in  as follows:

Equation 57 

EXTRAN uses an iterative solution to solve the gradually varied unsteady flow (St. Venant) equations. The convergence criteria for this method is related to the maximum number of iterations and the relative accuracy of nodal and conduit computations as specified in the simulation tolerances dialog.

The weighting coefficient for the non-linear momentum equation is calculated based on the conduit Froude number. Furthermore, when both ends of the conduit are surcharged the non-linear acceleration term is neglected by the program.

The finite difference approximation for the  term is calculated as follows in EXTRAN:

First Iteration:

Equation 58:  


Subsequent Iterations:

Equation 59:  


The model does not allow a flow reversal in one time step. The intervening flow is 1/1000 the previous conduit flow. This prevents oscillations in the conduit from arising from alternating positive and negative flows. This reversal check is performed on all conduits as well as orifices, pumps and weirs. The minimum flow in any conduit or diversion is set at 1.0e-10 (cms or cfs) to prevent underflow and overflow errors using MICROSOFT FORTRAN.

Control over the flow direction is also possible with internal flap gates in conduits that only allow flow either from the upstream to the downstream node or from the downstream to the upstream node.

The vertical differentiation of conduit roughness as defined by the special conduit factor dialog is a quadratic equation ranging from the peak at zero depth to normal roughness at the transition depth.

Initial transients in the model are dampened by increasing the roughness coefficient of each conduit by a factor of ten. This optional technique decreases oscillations at the beginning of a simulation. If used, the transition to normal roughness takes 9 large time steps .

The combined momentum-continuity equation uses the difference of the average cross-section area (A) at the n-1 and n+1 time steps with the denominator being the larger plus the sum of the smaller time steps . The numerator uses the average velocity at the current iteration.