Weir formula (Water Overlay): Difference between revisions

From Tygron Support wiki
Jump to navigation Jump to search
No edit summary
No edit summary
 
(54 intermediate revisions by the same user not shown)
Line 1: Line 1:
Flow across [[Weir (Water Overlay)|weir]]s is calculated differently for free flow and submerged flow.
Flow across [[Weir (Water Overlay)|weir]]s 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.  


The height of the water at each end of the weir, relative to the weir, is calculated:
First the upstream water level is calculated as followed:
: <math>h_s = max(0, max( w_l, w_r ) - z_{w,t})</math>
: <math>h_d = max(0, min( w_l, w_r ) - z_{w,t})</math>


Optionally, when a Weir drop threshold is configured, the weir height <math>z_w,t</math> at time <math>t</math> can update depending on the water level and the remaining weir drop time:
:<math>w_{u,t} =  max( w_{l,t}, w_{r,t})</math>
:<math>z_{w,t} =
 
Next, the adjusted weir height is determined. It either originates from the supplied weir height(s) <math>z_{w,t}</math> 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. 
 
:<math>z^{*}_{w,t} =  
\begin{cases}
\begin{cases}
z_b, & \mbox{if }t_{wd,t} > 0 \\
z_{th}, & \text{if} & \|{w_{u,t} - \tau_w}\| > \mu \text{ and } \tau_w > -10000 \\  
z_w, & \mbox{otherwise }  
z^{*}_{w,t-1}, & \text{if} & T_{wm} > T \\
z_{w,t}, & \text{otherwise}  
\end{cases}
\end{cases}
</math>
</math>
:<math>t_{wd,t} =  
 
The height adjustment range values, defined by the min and max, are determined next:
 
:<math>z_{min,w} = max ( z_{w,t} - \rho, min (z_{b,l}, z_{b,r} ) )</math>
 
:<math>z_{max,w} =  
\begin{cases}
\begin{cases}
t_{wd}, & \mbox{if } max( w_l, w_r) > w_{dt} \mbox{ and } t_{wd,t} <= 0 \\
z_{w,t}, & \text{if} &  w_{u,t} < z^{*}_{w,t-1} \\
t_{wd, t-1} - \Delta t, & \mbox{otherwise}
z_{w,t} + \rho, & \text{otherwise}
\end{cases}
</math>
 
Finally, the adjusted weir height is calculated and stored, as well as the moment in time at which the weir can be adjusted again.
:<math>z_{th} =
\begin{cases}
min(z_{max,w}, max( z_{min,w}, z^{*}_{w,t-1} + \mu )), & \text{if} & w_{u,t} < \tau_w \\
min(z_{max,w}, max( z_{min,w}, z^{*}_{w,t-1} - \mu )), & \text{otherwise}
\end{cases}
\end{cases}
</math>
</math>
:<math> T_{wm} = T + t_{wm}</math>
:<math>z_{w,t} = z^{*}_{w,t}</math>
Now knowing the height of the weir, the height of the water at each end of the weir, relative to the weir, is calculated:
: <math>h_s = max(0, max( w_{l,t}, w_{r,t} ) - z_{w,t})</math>
: <math>h_d = max(0, min( w_{l,t}, w_{r,t} ) - z_{w,t})</math>




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 31: Line 53:


For free flow, it is calculated directly:
For free flow, it is calculated directly:
: <math>Q_f = f_{w} \cdot c_w \cdot b \cdot ( h_s - h_d )^{3/2}</math>
: <math>Q_f = f_{dw} \cdot C_w \cdot b \cdot ( h_s - h_d )^{3/2}</math>


For submerged flow, the following calculation is used:
For submerged flow, the following calculation is used:
: <math>Q_s = U_{loss} \cdot A \cdot \sqrt{ 2 \cdot g \cdot ( h_s - h_d ) }</math>
: <math>Q_s = f_{loss} \cdot A \cdot \sqrt{ 2 \cdot g \cdot ( h_s - h_d ) }</math>


Finally the actual amount of water flow is calculated:
with:
: <math>\Delta w = \frac{\Delta t \cdot Q}{\Delta x}</math>
: <math>A = b \cdot (h_s-h_d)</math>
 
Finally the actual amount of water level change is calculated:
: <math>\Delta w = \frac{\Delta t \cdot Q}{\Delta x \cdot \Delta x}</math>


Where:
Where:
* <math>w_{l,t}</math> = The water level on the left side of the weir, relative to {{datum}}, at time <math>t</math>.
* <math>w_{l,t}</math> = The water level on the left side of the weir, relative to {{datum}}, at time <math>t</math>.
* <math>w_{r,t}</math> = The water level on the right side of the weir, relative to {{datum}}, at time <math>t</math>.
* <math>w_{r,t}</math> = The water level on the right side of the weir, relative to {{datum}}, at time <math>t</math>.
* <math>w_{u,t}</math> = The calculated upstream water level, relative to {{datum}}, at time <math>t</math>.
* <math>z_{th}</math> = The height of the weir at time t, according to the height adjustment mechanism.
* <math>\tau_w</math> = The [[Weir target level (Water Overlay)|WEIR_TARGET_LEVEL]] of the the weir.
* <math>T_{wm}</math> = Moment in time in seconds after which the weir can adjust its height again.
* <math>z_{w,t}</math> = The [[Weir height (Water Overlay)|WEIR_HEIGHT]] of the weir, at time <math>t</math>.
* <math>z_{w,t}</math> = The [[Weir height (Water Overlay)|WEIR_HEIGHT]] of the weir, at time <math>t</math>.
* <math>z^{*}_{w,t}</math> = The optionally adjusted height of the weir at time t, depending on the weir height adjustment mechanism.
* <math>z_{b,l}</math> = The water bottom elevation at left side of the weir the weir.
* <math>z_{b,r}</math> = The water bottom elevation at the right side of the weir.
* <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>\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 (Water Overlay)|WEIR_MOVE_RANGE_M]] of the water overlay, applicable to all weirs.
* <math>T</math> = Current simulated time in seconds.
* <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>w_{dt}</math> = The [[Weir drop threshold (Water Overlay)|weir drop threshold]] of the weir.
* <math>r_h</math> = The ratio of water heights on either side of the culvert.
* <math>z_{w,t}</math> = The height of the weir at time t, depending on drop mechanism.
* <math>Q</math> = The potential rate of water flow across the weir.
* <math>z_w</math> = The [[Weir height (Water Overlay)|WEIR_HEIGHT]] of the weir.
* <math>Q_f</math> = The potential rate of water flow across the weir, based on a free flow calculation.
* <math>z_b</math> = The water bottom elevation at the weir.
* <math>Q_s</math> = The potential rate of water flow across the weir, based on a submerged calculation.
* <math>t_{wd}</math> = The [[Weir drop time s (Water Overlay)|total amount of time]] a weir is dropped until the drop threshold condition is re-evaluated.
* <math>f_{dw}</math> = Dutch weir factor, set to 1.7.  
* <math>t_{wd,t}</math> = The remaining time the weir is dropped.
 
* <math>f_{w}</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>Q</math> = The potential rate of water flow across the weir.
* <math>f_{loss}</math> = Loss coefficient for submerged weirs, set to 0.9.
* <math>r_h</math> = The ratio of water heights on either side of the culvert.
* <math>Q_{free}</math> = The potential rate of water flow across the weir, based on a free flow calculation.
* <math>Q_{submerged}</math> = The potential rate of water flow across the weir, based on a submerged calculation.
* <math>U_{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===
To clarify the free flow formula in comparison with other definitions:
To clarify the free flow formula in comparison with the ISO{{ref|ISO}} standard definitions:
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>
For <math>h</math> 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.


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.
* <math>g</math> = acceleration due to gravity.
* <math>g</math> = acceleration due to gravity.
* <math>h</math> = Weir head; the height of the water column relative to the top of the weir.
* <math>h</math> = Weir head; the height of the water column relative to the top of the weir.
For <math>h</math> 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==
==Related==
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
Line 102: Line 140:
   }}
   }}
</references>
</references>
}}


{{WaterOverlay formula nav}}
{{WaterOverlay formula nav}}

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)