Weir formula (Water Overlay)
Revision as of 10:04, 15 April 2019 by Frank@tygron.nl (talk | contribs)
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:
- hs = max(0, max( wl, wr ) - zw)
- hd = max(0, min( wl, wr ) - zw)
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:
- hratio = hd / hs
- C = min( Csubmerged , Cfree) if hratio > 0,5
- C = Cfree otherwise
For free flow, capacity is calculated directly:
- Cfree = fw,d * cw * ww * ( hs - hd )3/2
For submerged flow, a culvert-like calculation is used:
- Csubmerged = Uloss * A * sqrt( 2 * g * ( hs - hd ) )
Finally the actual amount of water flow is calculated:
- Δw = Δt * C / cell
Where:
- hs = The height of the water column relative to the top of the weir, on the side with the highest water level.
- hd = The height of the water column relative to the top of the weir, on the side with the lowest water level.
- wl = The water level on the left side of the weir, relative to datum.
- wr = The water level on the right side of the weir, relative to datum.
- zw = The WEIR_HEIGHT attribute of the weir.
- fw,d = Dutch weir factor, set to 1.7.
- cw = The WEIR_COEFFICIENT attribute of the weir.
- ww = The WEIR_WIDTH attribute of the weir.
- Cfree = The potential rate of water flow across the weir, based on a free flow calculation.
- Csubmerged = The potential rate of water flow across the weir, based on a submerged calculation.
- C = The potential rate of water flow across the weir.
- hratio = The ratio of water heights on either side of the culvert.
- Uloss = Loss coefficient for submerged weirs (0,9).
- A = Flow area, based on the highest water column height relative to the top of the weir, and the weir width, defined by the WEIR_WIDTH attribute of the weir
- g = Acceleration factor of gravity.
- Δw = The water flow which takes place.
- Δt = Computational timestep.
- cell = Cell size.