Water surface profiles are computed from one cross section to the next by solving the Energy equation with an iterative procedure called the standard step method. The Energy equation is written as follows:

1) Z_2+Y_2+\frac{\alpha_2 V_2^2}{2g} = Z_1+Y_1+\frac{\alpha_1 V_1^2}{2g}+h_e
SymbolDescriptionUnits

Z_1, Z_2

elevation of the main channel inverts

Y_1, Y_2

depth of water at cross sections

V_1, V_2

average velocities (total discharge/ total flow area)

a_1, a_2

velocity weighting coefficients

g

gravitational acceleration

h_e

energy head loss

A diagram showing the terms of the energy equation is shown in the figure below.

Representation of Terms in the Energy Equation

The energy head loss (h_e) between two cross sections is comprised of friction losses and contraction or expansion losses. The equation for the energy head loss is as follows:

2) h_e=L\overline{S}_f+C \displaystyle\left( \frac{\alpha_2 V_2^2}{2g}-\frac{\alpha_1 V_1^2}{2g} \right)
SymbolDescriptionUnits

L

discharge weighted reach length

\overline{S}_f

representative friction slope between two sections

C

expansion or contraction loss coefficient

The discharge weighted reach length, L, is calculated as:

3) L=\frac{L_{lob}\ \overline{Q}_{lob}+L_{ch}\ \overline{Q}_{ch}+L_{rob}\ \overline{Q}_{rob}}{\overline{Q}_{lob}+\overline{Q}_{ch}+\overline{Q}_{rob}}
SymbolDescriptionUnits

L_{lob}, L_{ch}, L_{rob}

cross section reach lengths specified for flow in the left overbank, main channel, and right overbank, respectively

\overline{Q_{lob}} + \overline{Q_{ch}} + \overline{Q_{rob}}

arithmetic average of the flows between sections for the left overbank, main channel, and right overbank, respectively