Full Equations Utilities (FEQUTL) Model for the Approximation of Hydraulic Characteristics of Open Channels and Control Structures During Unsteady Flow

# 5.6 EMBANKQ Command

Update available for vertical scale factor, vertical shift, and horizontal scale factor

Purpose: The Hulsing (1967) procedure for flow over embankment-shaped weirs is implemented as outlined in section 4.3. Flow over a variety of weirs can be computed if the appropriate tables are provided.

Notes: A warning message has been added to EMBANKQ to alert the user to potentially unrealistic results for flow over the embankment when the height of the embankment above the approach section becomes small relative to the head on the embankment. If the ratio of piezometric head to the weir height becomes much greater than 4 or 5, the coefficients used for flow over the weir are inaccurate. The depth at which the weir-flow coefficients become inaccurate is not known. However, if the head to height ratio goes above 7, the weir-flow coefficients are certainly inaccurate. The message is only a warning, and the flow will be computed in FEQUTL for whatever ratio is present. The user must decide if the result is reasonable. The flow at each point along the crest of the embankment is limited to be no more than the critical flow rate in the approach section. Therefore, when weir subsections with small heights are present in the stream system, the critical flow in the approach section is applied. However, no warning is issued when this is done. The best method to apply for weir subsections with small heights is to divide the embankment into subsections and apply CHANRAT (section 5.3) for those sections with crest elevations very close to the approach section elevation.
LINE 1
Variable: TABLE, TYPE, HLCRIT, HLMAX
Format: 7X, I5, A5, A5
Example: TABLE #= 9900
Explanation:
 TABLE is table number for the table to be computed in FEQUTL. TYPE is the table type and can be either 6 or 13 with 13 applied by default if TYPE is omitted. HLCRIT is the critical ratio between approach head and embankment width that distinguishes between high-head and low-head flow. Ratios at or greater than HLCRIT correspond to high-head flow and ratios less than HLCRIT correspond to low-head flow. The default value is 0.15 based on Hulsing (1967) and Kindsvater (1964). HLMAX is the maximum ratio of approach head to embankment width that can be applied without generating a warning message in EMBANKQ. This warning message alerts the user to the possibility that the weir coefficients may no longer be valid at head ratios greater than HLMAX. The approximate upper limit of the head-to-width ratio in the experiments defining the weir coefficients is represented in HLMAX. HLMAX defaults to 0.32 based on Hulsing (1967) and Kindsvater (1964).
LINE 2
Variable: TAB(1)
Format: 7X, I5
Example: PLCWTB = 9994
Explanation:
 TAB(1) is the table number for the function table listing the low-head weir coefficient for a paved surface. The table number listed in this example is the default value, and this table may be retrieved electronically as described in section 5.1.
LINE 3
Variable: TAB(2)
Format: 7X, I5
Example: GLCWTB = 9995
Explanation:
 TAB(2) is the table number for the function table listing the low-head weir coefficient for a graveled surface. The table number listed in this example is the default value, and this table may be retrieved electronically as described in section 5.1.
LINE 4
Variable: TAB(3)
Format: 7X, I5
Example: PHCWTB = 9996
Explanation:
 TAB(3) is the table number for the function table listing the high-head weir coefficient for a paved surface. The table number listed in this example is the default value, and this table may be retrieved electronically as described in section 5.1.
LINE 5
Variable: TAB(4)
Format: 7X, I5
Example: GHCWTB = 9997
Explanation:
 TAB(4) is the table number for the function table listing the high-head weir coefficient for a graveled surface. The table number listed in this example is the default value, and this table may be retrieved electronically as described in section 5.1.
LINE 6
Variable: TAB(5)
Format: 7X, I5
Example: PSUBTB = 9998
Explanation:
 TAB(5) is the table number for the function table listing the submergence correction factor for a paved surface. The table number listed in this example is the default value, and this table may be retrieved electronically as described in section 5.1.
LINE 7
Variable: TAB(6)
Format: 7X, I5
Example: GSUBTB = 9999
Explanation:
 TAB(6) is the table number for the function table listing the submergence correction factor for the graveled surface. The table number listed in this example is the default value, and this table may be retrieved electronically as described in section 5.1. If PSUBTB and GSUBTB are given as zero, then no submergence computations are done and the resulting function table will be type 2. If they are nonzero, then submergence computations are done and the resulting function table will be type 13.
LINE 8
Variable: LABEL
Format: 6X, A80
Example: LABEL = FLOW OVER BUTTERFIELD ROAD AT DRY CREEK
Explanation:
 LABEL is a user-supplied label that will be output as part of the heading for the table produced in FEQUTL computations.
LINE 9
Format: A80
Example: OFFSET CREST WIDTH APPROACH SURFACE
Explanation:
LINE 10
Variables: OFF(I), CREST(I), WIDTH(I), APPROC(I), SURF(I)
Format: 4F10.0, 1X, A8
Explanation:
 The values given on this line define, in a subsection-by-subsection approach, the geometry of the crest of the embankment-shaped weir applied in computation of the flow over the embankment. OFF is the horizontal offset. CREST is the crest elevation. WIDTH is the width of the crest in the direction of flow. APPROC is the elevation of the approach channel. If the approach velocity should be ignored in the computations, then an approach elevation that is much lower than the crest elevation should be specified so that the velocity head computed will be small. SURF is the type of surface, and the designation applies to the line segment beginning at the offset for which the surface type is given. For example, if the first line of the weir-crest specification gives the type of surface as GRAVEL and the second line of the weir-crest specification gives the type of surface as PAVED, then the weir crest between the first and second offset is taken to be of a roughness similar to a graveled roadway. Two closely spaced points should be placed on the crest at the change in roughness to represent the transition between the two surface conditions. The end of the weir specification is indicated by giving END for the type of surface. Thus, if the third line of the input in this example had END for the type of surface, then the weir crest between the second and third offset would have a roughness similar to a paved roadway.

These geometric characteristics are input for the boundaries of the subsections. The crest elevation, crest width, and approach-channel elevation are assumed to change linearly between the given points on the crest, as illustrated in figure 18. To make input easier, the values for width, approach-channel elevation, and surface type propagate downward into fields left blank. Thus, if the elevation of the approach channel is constant and the width is constant, only the first line of the specification must contain the elevation of the approach channel and the width of the weir crest. The surface value also propagates so that for constant surface type, only the first line must have PAVED or GRAVEL. END is specified for the final line to indicate the end of geometric input.

LINE 11
Format: A80
Example: UPSTREAM HEADS TO USE IN COMPUTING THE TABLE
Explanation:

The following two lines (11a and 11b) are optional. However, if one is applied, the other also must be applied. If they are omitted, NFRAC defaults to 21 and POWER defaults to 2.

LINE 11a
Variable: CHAR5, NFRAC
Format: A5, 1X, I5
Example: NFRAC = 11
Explanation:
 CHAR5 must be NFRAC. NFRAC is the number of partial free drops used in computing the table. Complete details on partial free drops are given in the discussion of 2-D function tables in section 11.2 of the documentation report for FEQ (Franz and Melching, 1997).
LINE 11b
Variable: CHAR6, POWER
Format: A6, 1X, F10.0
Example: POWER = 2.
Explanation:
 CHAR6 must be POWER =. POWER is the power applied to distribute the proportion of free drop from a value of 0 to 1. The proportion of free drop is given by for i = 1, ..., NFRAC.
LINE 12 (Repeated as necessary to input all the upstream heads used in computing the 2-D table.)
Variable: HUVEC(*)
Format: F10.0
Explanation: