Weir formula (Water Overlay): Difference between revisions

From Tygron Support wiki
Jump to navigation Jump to search
No edit summary
No edit summary
 
(17 intermediate revisions by the same user not shown)
Line 9: Line 9:
:<math>z^{*}_{w,t} =  
:<math>z^{*}_{w,t} =  
\begin{cases}
\begin{cases}
z_{th}, & \text{if} & \|{w_{u,t} - \tau_w}\| > \mu \\
z_{th}, & \text{if} & \|{w_{u,t} - \tau_w}\| > \mu \text{ and } \tau_w > -10000 \\  
z^{*}_{w,t-1}, & \text{if} & T_{wm} > T \\
z^{*}_{w,t-1}, & \text{if} & T_{wm} > T \\
z_{w,t}, & \text{otherwise}  
z_{w,t}, & \text{otherwise}  
Line 42: Line 42:




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:
Based on the relative water heights, the weir uses either a submerged flow <math>Q_s</math> or a free flow <math>Q_f</math> formula, based on the following ratio:
: <math>r_h = \frac{h_d}{h_s}</math>
: <math>r_h = \frac{h_d}{h_s}</math>
:<math>
:<math>
Line 61: Line 61:
: <math>A = b \cdot (h_s-h_d)</math>
: <math>A = b \cdot (h_s-h_d)</math>


Finally the actual amount of water flow is calculated:
Finally the actual amount of water level change is calculated:
: <math>\Delta w = \frac{\Delta t \cdot Q}{\Delta x}</math>
: <math>\Delta w = \frac{\Delta t \cdot Q}{\Delta x \cdot \Delta x}</math>


Where:
Where:
Line 78: Line 78:
* <math>z_{min,w}</math> = The minimum allowed height for weir w.
* <math>z_{min,w}</math> = The minimum allowed height for weir w.
* <math>z_{max,w}</math> = The maximum allowed height for weir w.
* <math>z_{max,w}</math> = The maximum allowed height for weir w.
* <math>\mu</math> = The [[Weir move step m model attribute (Water Overlay)|WEIR_MOVE_STEP_M]] of the water overlay, applicable to all weirs.
* <math>\mu</math> = The [[Weir move step m (Water Overlay)|WEIR_MOVE_STEP_M]] of the water overlay, applicable to all weirs.
* <math>\rho</math> = The [[Weir move range m model attribute (Water Overlay)|WEIR_MOVE_RANGE_M]] of the water overlay, applicable to all weirs.
* <math>\rho</math> = The [[Weir move range m (Water Overlay)|WEIR_MOVE_RANGE_M]] of the water overlay, applicable to all weirs.
* <math>T</math> = Current simulated time in seconds.
* <math>T</math> = Current simulated time in seconds.
* <math>t_{wm}</math> = The [[Weir move interval model attribute (Water Overlay)|Weir move interval]] of the water overlay, applicable to all weirs.
* <math>t_{wm}</math> = The [[Weir move interval s (Water Overlay)|WEIR_MOVE_INTERVAL_S]] of the water overlay, applicable to all weirs.


* <math>h_s</math> = The height of the water column relative to the top of the weir, on the side with the highest [[Surface water level formula (Water Overlay)|water level]], at time <math>t</math>.
* <math>h_s</math> = The height of the water column relative to the top of the weir, on the side with the highest [[Surface water level formula (Water Overlay)|water level]], at time <math>t</math>.
* <math>h_d</math> = The height of the water column relative to the top of the weir, on the side with the lowest [[Surface water level formula (Water Overlay)|water level]], at time <math>t</math>.
* <math>h_d</math> = The height of the water column relative to the top of the weir, on the side with the lowest [[Surface water level formula (Water Overlay)|water level]], at time <math>t</math>.
* <math>r_h</math> = The ratio of water heights on either side of the culvert.
* <math>Q</math> = The potential rate of water flow across the weir.
* <math>Q</math> = The potential rate of water flow across the weir.
* <math>Q_f</math> = The potential rate of water flow across the weir, based on a free flow calculation.
* <math>Q_f</math> = The potential rate of water flow across the weir, based on a free flow calculation.
Line 90: Line 91:
* <math>f_{dw}</math> = Dutch weir factor, set to 1.7.  
* <math>f_{dw}</math> = Dutch weir factor, set to 1.7.  
* <math>C_w</math> = The [[Weir coefficient (Water Overlay)|WEIR_COEFFICIENT]] of the weir.
* <math>C_w</math> = The [[Weir coefficient (Water Overlay)|WEIR_COEFFICIENT]] of the weir.
* <math>C_d</math> = Corresponding Coefficient of discharge. Part of the Dutch weir factor.
* <math>b</math> = The breadth of weir crest, adjustable using the [[Weir width (Water Overlay)|WEIR_WIDTH]].
* <math>b</math> = The breadth of weir crest, adjustable using the [[Weir width (Water Overlay)|WEIR_WIDTH]].
* <math>r_h</math> = The ratio of water heights on either side of the culvert.
* <math>f_{loss}</math> = Loss coefficient for submerged weirs, set to 0.9.
* <math>f_{loss}</math> = Loss coefficient for submerged weirs, set to 0.9.
* <math>A</math> = Flow area, based on the highest water column height relative to the top of the weir, and [[Weir width (Water Overlay)|WEIR_WIDTH]] of the weir.
* <math>A</math> = Flow area, based on the highest water column height relative to the top of the weir, and [[Weir width (Water Overlay)|WEIR_WIDTH]] of the weir.
* <math>g</math> = acceleration due to gravity.
* <math>g</math> = {{gravity}}.
* <math>\Delta w</math> = The water flow which takes place.
* <math>\Delta w</math> = The water level change in meters, which takes place.
* <math>\Delta t</math> = Computational timestep.
* <math>\Delta t</math> = Computational [[Timestep formula (Water Overlay)|timestep]] in seconds.
* <math>\Delta x</math> = Cell size.
* <math>\Delta x</math> = Cell size in meters.


===Free flow coefficients===
===Free flow coefficients===
Line 104: Line 103:
When <math>h \approx 0.01</math> meters, the following holds.
When <math>h \approx 0.01</math> meters, the following holds.


:<math>f_{w} \cdot c_w \approx C_d \cdot \sqrt{g}</math>
:<math>f_{dw} \cdot C_w \approx C_d \cdot \sqrt{g}</math>


:<math>C_d=0,633 \cdot (1-\frac{0.0003}{h})^\frac{3}{2}</math>
:<math>C_d=0,633 \cdot (1-\frac{0.0003}{h})^\frac{3}{2}</math>


where:
where:
* <math>f_{w}</math> = Dutch weir factor, set to 1.7.  
* <math>f_{dw}</math> = Dutch weir factor, set to 1.7.  
* <math>C_w</math> = The [[Weir coefficient (Water Overlay)|WEIR_COEFFICIENT]] of the weir, set to a default of 1.1
* <math>C_w</math> = The [[Weir coefficient (Water Overlay)|WEIR_COEFFICIENT]] of the weir, set to a default of 1.1
* <math>C_d</math> = Calculated Coefficient of discharge.
* <math>C_d</math> = Calculated Coefficient of discharge.
Line 120: Line 119:
The following topics are related to this formula.
The following topics are related to this formula.
; Structures
; Structures
: [[Weir_(Water_Overlay)|Weir]]
: [[Weir (Water Overlay)|Weir]]
: [[Breach (Water Overlay)|Breach]]
; Models
; Models
: [[Surface model (Water Overlay)|Surface model]]
: [[Surface model (Water Overlay)|Surface model]]
; Formulas
: [[Breach flow formula (Water Overlay)|Breach flow formula]]
; Test cases
: [[Weir height test case (Water Module)|Weir height test case]]


{{Article end
{{Article end

Latest revision as of 13:21, 5 March 2024

Flow across weirs is calculated differently for free flow and submerged flow. Optionally, the height of the weir can variate based on a provided height values or an automatic adjustment. Therefore, we determined the height of the weir first.

First the upstream water level is calculated as followed:

Next, the adjusted weir height is determined. It either originates from the supplied weir height(s) or it is adjusted according to the weir threshold level. When it is adjusted, a moment in time is set for that weir, during which it cannot be adjusted.

The height adjustment range values, defined by the min and max, are determined next:

Finally, the adjusted weir height is calculated and stored, as well as the moment in time at which the weir can be adjusted again.

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


Based on the relative water heights, the weir uses either a submerged flow or a free flow formula, based on the following ratio:

For free flow, it is calculated directly:

For submerged flow, the following calculation is used:

with:

Finally the actual amount of water level change is calculated:

Where:

  • = The water level on the left side of the weir, relative to datum, at time .
  • = The water level on the right side of the weir, relative to datum, at time .
  • = The calculated upstream water level, relative to datum, at time .
  • = The height of the weir at time t, according to the height adjustment mechanism.
  • = The WEIR_TARGET_LEVEL of the the weir.
  • = Moment in time in seconds after which the weir can adjust its height again.
  • = The WEIR_HEIGHT of the weir, at time .
  • = The optionally adjusted height of the weir at time t, depending on the weir height adjustment mechanism.
  • = The water bottom elevation at left side of the weir the weir.
  • = The water bottom elevation at the right side of the weir.
  • = The minimum allowed height for weir w.
  • = The maximum allowed height for weir w.
  • = The WEIR_MOVE_STEP_M of the water overlay, applicable to all weirs.
  • = The WEIR_MOVE_RANGE_M of the water overlay, applicable to all weirs.
  • = Current simulated time in seconds.
  • = The WEIR_MOVE_INTERVAL_S of the water overlay, applicable to all weirs.
  • = The height of the water column relative to the top of the weir, on the side with the highest water level, at time .
  • = The height of the water column relative to the top of the weir, on the side with the lowest water level, at time .
  • = The ratio of water heights on either side of the culvert.
  • = The potential rate of water flow across the weir.
  • = The potential rate of water flow across the weir, based on a free flow calculation.
  • = The potential rate of water flow across the weir, based on a submerged calculation.
  • = Dutch weir factor, set to 1.7.
  • = The WEIR_COEFFICIENT of the weir.
  • = The breadth of weir crest, adjustable using the WEIR_WIDTH.
  • = Loss coefficient for submerged weirs, set to 0.9.
  • = Flow area, based on the highest water column height relative to the top of the weir, and WEIR_WIDTH of the weir.
  • = the acceleration due to gravity, set to 9.80665.
  • = The water level change in meters, which takes place.
  • = Computational timestep in seconds.
  • = Cell size in meters.

Free flow coefficients

To clarify the free flow formula in comparison with the ISO[1] standard definitions: When meters, the following holds.

where:

  • = Dutch weir factor, set to 1.7.
  • = The WEIR_COEFFICIENT of the weir, set to a default of 1.1
  • = Calculated Coefficient of discharge.
  • = acceleration due to gravity.
  • = Weir head; the height of the water column relative to the top of the weir.

For generally larger than that, the flow can be calculated up to 5.5% smaller than expected when using the second pair of coefficients. However, this can be corrected with the weir coefficient if necessary.

Related

The following topics are related to this formula.

Structures
Weir
Breach
Models
Surface model
Formulas
Breach flow formula
Test cases
Weir height test case


References

  1. ISO FDIS 4360, 2020 Edition, March 3, 2020 - HYDROMETRY - OPEN CHANNEL FLOW MEASUREMENT USING TRIANGULAR PROFILE WEIRS ∙ Technical Committee: ISO/TC 113/SC 2 Flow measurement structures ∙ Found at: https://www.iso.org/standard/70915.html ∙ (last visited: 28-11-2022)