Critical Depth Determination

Critical depth for a cross section will be determined if any of the following conditions are satisfied:

1. The supercritical flow regime has been specified.
2. The calculation of critical depth has been requested by the user.
3. This is an external boundary cross section and critical depth must be determined to ensure the user entered boundary condition is in the correct flow regime.
4. The Froude number check for a subcritical profile indicates that critical depth needs to be determined to verify the flow regime associated with the balanced elevation.
5. The program could not balance the energy equation within the specified tolerance before reaching the maximum number of iterations.

The total energy head for a cross section is defined by:

The critical water surface elevation is the elevation for which the total energy head is a minimum (i.e., minimum specific energy for that cross section for the given flow). The critical elevation is determined with an iterative procedure whereby values of WS are assumed and corresponding values of H are determined with (1) until a minimum value for H is reached.

The HEC-RAS program has two methods for calculating critical depth: a "parabolic" method and a "secant" method. The parabolic method is computationally faster, but it is only able to locate a single minimum energy. For most cross sections there will only be one minimum on the total energy curve, therefore the parabolic method has been set as the default method (the default method can be changed from the user interface). If the parabolic method is tried and it does not converge, then the program will automatically try the secant method.

In certain situations it is possible to have more than one minimum on the total energy curve. Multiple minimums are often associated with cross sections that have breaks in the total energy curve. These breaks can occur due to very wide and flat overbanks, as well as cross sections with levees and ineffective flow areas. When the parabolic method is used on a cross section that has multiple minimums on the total energy curve, the method will converge on the first minimum that it locates. This approach can lead to incorrect estimates of critical depth. If the user thinks that the program has incorrectly located critical depth, then the secant method should be selected and the model should be re-simulated.

The "parabolic" method involves determining values of H for three values of WS that are spaced at equal ΔWS intervals. The WS corresponding to the minimum value for H, defined by a parabola passing through the three points on the H versus WS plane, is used as the basis for the next assumption of a value for WS. It is presumed that critical depth has been obtained when there is less than a 0.01 ft. (0.003 m) change in water depth from one iteration to the next and provided the energy head has not either decreased or increased by more than .01 feet (0.003 m).

The "secant" method first creates a table of water surface versus energy by slicing the cross section into 30 intervals. If the maximum height of the cross section (highest point to lowest point) is less than 1.5 times the maximum height of the main channel (from the highest main channel bank station to the invert), then the program slices the entire cross section into 30 equal intervals. If this is not the case, the program uses 25 equal intervals from the invert to the highest main channel bank station, and then 5 equal intervals from the main channel to the top of the cross section. The program then searches this table for the location of local minimums. When a point in the table is encountered such that the energy for the water surface immediately above and immediately below are greater than the energy for the given water surface, then the general location of a local minimum has been found. The program will then search for the local minimum by using the secant slope projection method. The program will iterate for the local minimum either thirty times or until the critical depth has been bounded by the critical error tolerance. After the local minimum has been determined more precisely, the program will continue searching the table to see if there are any other local minimums. The program can locate up to three local minimums in the energy curve. If more than one local minimum is found, the program sets critical depth equal to the one with the minimum energy. If this local minimum is due to a break in the energy curve caused by overtopping a levee or an ineffective flow area, then the program will select the next lowest minimum on the energy curve. If all of the local minimums are occurring at breaks in the energy curve (caused by levees and ineffective flow areas), then the program will set critical depth to the one with the lowest energy. If no local minimums are found, then the program will use the water surface elevation with the least energy. If the critical depth that is found is at the top of the cross section, then this is probably not a real critical depth. Therefore, the program will double the height of the cross section and try again. Doubling the height of the cross section is accomplished by extending vertical walls at the first and last points of the section. The height of the cross section can be doubled five times before the program will quit searching.