Weir formula (Water Overlay)

Flow across weirs is calculated differently for free flow and submerged flow.

The height of the water at each end of the weir, relative to the weir, is calculated:

h_s = max(0, max( w_l, w_r ) - z_w)

h_d = max(0, min( w_l, w_r ) - z_w)

Based on the relative water heights, the weir is judged to have either a submerged flow or a free flow, based on the following ratio:

h_ratio = h_d / h_s

C = min( C_submerged , C_free) if h_ratio > 0,5

C = C_free otherwise

For free flow, capacity is calculated directly:

C_free = f_w,d * c_w * w_w * ( h_s - h_d )^3/2

For submerged flow, the following calculation is used:

C_submerged = U_loss * A * sqrt( 2 * g * ( h_s - h_d ) )

Finally the actual amount of water flow is calculated:

Δw = Δt * C / Δx

Where:

h_s = The height of the water column relative to the top of the weir, on the side with the highest water level.
h_d = The height of the water column relative to the top of the weir, on the side with the lowest water level.
w_l = The water level on the left side of the weir, relative to datum.
w_r = The water level on the right side of the weir, relative to datum.
z_w = The WEIR_HEIGHT of the weir.
f_w,d = Dutch weir factor, set to 1.7.
c_w = The WEIR_COEFFICIENT of the weir.
w_w = The WEIR_WIDTH of the weir.
C = The potential rate of water flow across the weir.
h_ratio = The ratio of water heights on either side of the culvert.
C_free = The potential rate of water flow across the weir, based on a free flow calculation.
C_submerged = The potential rate of water flow across the weir, based on a submerged calculation.
U_loss = Loss coefficient for submerged weirs, set to 0.9.
A = Flow area, based on the highest water column height relative to the top of the weir, and WEIR_WIDTH of the weir.
g = Acceleration factor of GRAVITY, defined for the Water Overlay.
Δw = The water flow which takes place.
Δt = Computational timestep.
Δx = Cell size.

Related