Breach flow formula (Water Overlay): Difference between revisions

From Tygron Support wiki
Jump to navigation Jump to search
No edit summary
No edit summary
 
(26 intermediate revisions by the same user not shown)
Line 1: Line 1:
Flow through [[Breach (Water Overlay)|breach]]es is calculated based on the [[Weir formula (Water Overlay)|weir formula]], including the consideration between free flow and submerged flow situations.
Flow through [[Breach (Water Overlay)|breach]]es is calculated based on the [[Weir formula (Water Overlay)|weir formula]], including the consideration between free flow and submerged flow situations.
The weir formula is supplied with the breach width (as stand-in for the weir width), which is always calculated with the [[Breach growth formula (Water Overlay)|breach growth formula]] before calculating the weir flow. The downstream water level can optionally be measured with an additionally configured [[Breach level area_(Water_Overlay)|Breach Level Area]], or an automatically placed  
The weir formula is supplied with the breach width (as stand-in for the weir width), which is always calculated with the [[Breach growth formula (Water Overlay)|breach growth formula]] before calculating the weir flow. The downstream water level can optionally be measured with an additionally configured [[Breach level area_(Water_Overlay)|Breach Level Area]], or an automatically placed at a distance based on the [[Breach measurement distance m (Water Overlay)|configured measurement distance]] and the [[Breach angle (Water Overlay)|orientation]] of the breach.


For the entry area, the area where the water originates from, the following also applies:
For the entry area, the area where the water originates from, the following also applies:
Line 11: Line 11:
Each timestep, the external water level is changed based on the amount of water flowing in or out.
Each timestep, the external water level is changed based on the amount of water flowing in or out.


<math>\Delta V = min(h_e \cdot A_s , Q_w \cdot \Delta t)</math>
:<math>h_e = w_{e,t} - z_b</math>
In case <math>A_s == A_b</math>:
<math>w_{e,t+1} = w_{e,t} - \frac{\Delta V}{A_b}</math>
Otherwise, the external water area is represented as a trapezoid prism and the new water level has to be calculated differently:
:<math>l = \sqrt{A_s}</math>


:<math>w_s = l</math>
:<math>\Delta V = min(h_e \cdot A_s , Q_{w,t} \cdot \Delta t)</math>


:<math>w_b = \frac{A_b}{l}</math>
In case <math>A_s = A_b</math>:


:<math>a = \frac{w_s-w_b}{4 \cdot h_0}</math>
:<math>w_{e,t+1} = w_{e,t} - \frac{\Delta V}{A_b}</math>


:<math>A_{t,t} = 2 \cdot (a \cdot h_t^2 + 0.5 \cdot w_b)</math>
Where:
:<math>\Delta A = \frac{\Delta V}{l}</math>
: <math>z_{b}</math> = Datum height of the bottom of the area water flows from.
: <math>w_{e,t}</math> = The [[Surface water level formula (Water Overlay)|water level]] of the entry area at time <math>t</math>; In case of an external area: [[External water level (Water Overlay)|EXTERNAL_WATER_LEVEL]] of the [[Breach (Water Overlay)|breach]].
: <math>Q_{w,t}</math> =  The calculated water flow which takes place at time <math>t</math>, based on the [[Weir formula (Water Overlay)|weir formula]].
: <math>\Delta t</math> = The [[Timestep formula (Water Overlay)|timestep]] in seconds.
: <math>\Delta V</math> =  The limited volume of water transported, based on the remaining water.
: <math>A_{s}</math> = The size of in the entry area; In case of an external area: [[External area (Water Overlay)|EXTERNAL_AREA]] of the breach.
: <math>A_{b}</math> = The (optional) size of the bottom of the entry area; In case of an external area: the second value provided for the [[External area (Water Overlay)|EXTERNAL_AREA]] of the breach. In case only one value is provided, <math>A_{b} = A_{s}</math>
==Trapezoidal prism==
[[File:external_water_body_trapezoid.png|right|External water body as a trapezoid]]
The [[External area (Water Overlay)|external water area]] can also be represented as a trapezoidal prism. The new water level then has to be calculated differently.
 
The length and width of the surface area are currently considered equal. Therefore:
:<math>L = \sqrt{A_s}</math>
 
:<math>W_s = L</math>
 
The width of the bottom area is then calculated as:
:<math>W_b = \frac{A_b}{L}</math>
 
Next, we will calculate the new water level based on <math>\Delta V</math>:
 
:<math> h_0 = w_{e,0} - z_b</math>
 
:<math>a = \frac{W_s-W_b}{4 \cdot h_0}</math>
 
:<math>A_{t,t} = 2 \cdot (a \cdot h_t^2 + 0.5 \cdot W_b \cdot h_t)</math>
:<math>\Delta A = \frac{\Delta V}{L}</math>


:<math>c = \frac{-\Delta A * A_{t,t}}{2}</math>
:<math>c = \frac{-\Delta A * A_{t,t}}{2}</math>


:<math>D = (\frac{w_b}{2})^2 - 4 \cdot a \cdot c</math>
:<math>D = (\frac{W_b}{2})^2 - 4 \cdot a \cdot c</math>


if <math> D < 0</math>:
if <math> D < 0</math>, determine the new water level based on a box shape:


<math>w_{e,t+1} = w_{e,t} + \frac{\Delta V}{A_b}</math>
:<math>w_{e,t+1} = w_{e,t} + \frac{\Delta V}{A_b}</math>


else:
otherwise, determine the new water level based on the determinant and flow direction:


<math>w_{-} = \frac{-\frac{w_b}{2} - \sqrt{D}}{2 \cdot a} + z_b</math>
:<math>w_{-} = \frac{-\frac{W_b}{2} - \sqrt{D}}{2 \cdot a} + z_b</math>


<math>w_{+} = \frac{-\frac{w_b}{2} + \sqrt{D}}{2 \cdot a} + z_b</math>
:<math>w_{+} = \frac{-\frac{W_b}{2} + \sqrt{D}}{2 \cdot a} + z_b</math>


<math>w_{e,t+1} =  
:<math>w_{e,t+1} =  
\begin{cases}
\begin{cases}
w_{+} & \mbox{ if } \Delta V < 0 & \mbox{ and } w_{+} \le w_{e,t}\\
w_{+} & \mbox{ if } \Delta V < 0 & \mbox{ and } w_{+} \le w_{e,t}\\
Line 50: Line 72:


Where:
Where:
: <math>Q_w</math> = The calculated water flow which takes place, based on the [[Weir formula (Water Overlay)|weir formula]].
: <math>L</math> = Length of the trapezoidal prism
: <math>Q</math> = The limited water flow which takes place, based on the remaining water.
: <math>A_{s}</math> = The size of in the entry area; In case of an external area: [[External area (Water Overlay)|EXTERNAL_AREA]] of the breach.
: <math>W_s</math> = Width of the original top (surface at <math>t = 0</math>) of the trapezoid
: <math>A_{b}</math> = The (optional) size of the bottom of the entry area; In case of an external area: the second value provided for the [[External area (Water Overlay)|EXTERNAL_AREA]] of the breach. In case only one value is provided, <math>A_{b} = A_{s}</math>: <math>z_{b}</math> = Datum height of the bottom of the area water flows from.
: <math>W_b</math> = Width of the bottom of the trapezoid
: <math>w_{e,t}</math> = The [[Surface water level formula (Water Overlay)|water level]] of the entry area at time <math>t</math>; In case of an external area: [[External water level (Water Overlay)|EXTERNAL_WATER_LEVEL]] of the [[Breach (Water Overlay)|breach]].
: <math>w_{e,t}</math> = The [[Surface water level formula (Water Overlay)|water level]] of the entry area at time <math>t</math>; In case of an external area: [[External water level (Water Overlay)|EXTERNAL_WATER_LEVEL]] of the [[Breach (Water Overlay)|breach]].
: <math>A_{s}</math> = The size of in the entry area; In case of an external area: [[External area (Water Overlay)|EXTERNAL_AREA]] of the breach.
: <math>z_{b} = </math> Datum height of the bottom of the area water flows from.
: <math>A_{b}</math> = The (optional) size of the bottom of the entry area; In case of an external area: the second value provided for the [[External area (Water Overlay)|EXTERNAL_AREA]] of the breach. In case only one value is provided, <math>A_{b} = A_{s}</math>
: <math>h_0</math> = The initial relative water height in the external area.
: <math>z_{b}</math> Datum height of the bottom of the trapezoidal prism.
: <math>a</math> = Slope of the trapezoid
: <math>A_{t,t}</math> = Area of the trapezoidal slide of the trapezoid prism at time <math>t</math>
: <math>c</math> = Constant based on new area of the trapezoid <math>A_{t,t+1}</math>
: <math>D</math> = Polynomial discriminant
: <math>\Delta V</math> =  The limited volume of water transported, based on the remaining water.
 


==Notes==
{{article end
|notes=
*The water level w<sub>e</sub> cannot become lower than the surface height s<sub>e</sub>, defined by [[External_surface_level_(Water_Overlay)|EXTERNAL_SURFACE_LEVEL]] in case of an external entry area. In case of an input area within the project, the [[terrain height (Water Overlay)|surface height]] is determined per grid cell.
*The water level w<sub>e</sub> cannot become lower than the surface height s<sub>e</sub>, defined by [[External_surface_level_(Water_Overlay)|EXTERNAL_SURFACE_LEVEL]] in case of an external entry area. In case of an input area within the project, the [[terrain height (Water Overlay)|surface height]] is determined per grid cell.
*When water is flowing from the project area to an external area, the slope of the external area remains the same, even when the water level rises above the original water level. In the case of a trapezoidal prism, is extrapolated upwards.
}}


==Related==
==Related==

Latest revision as of 13:27, 5 March 2024

Flow through breaches is calculated based on the weir formula, including the consideration between free flow and submerged flow situations. The weir formula is supplied with the breach width (as stand-in for the weir width), which is always calculated with the breach growth formula before calculating the weir flow. The downstream water level can optionally be measured with an additionally configured Breach Level Area, or an automatically placed at a distance based on the configured measurement distance and the orientation of the breach.

For the entry area, the area where the water originates from, the following also applies:

  • The water level used is defined by either an input area or the external water level.
  • In case of an external area, the surface height used is defined by external surface level. This height limits how far the water level can be lowered.
  • Optionally a second value can be supplied for the external surface level. In that case, the second value represents the area of the bottom.
  • In case of an input area, the average water level is based on the water level on the individual grid cells of the input area.

Each timestep, the external water level is changed based on the amount of water flowing in or out.

In case :

Where:

= Datum height of the bottom of the area water flows from.
= The water level of the entry area at time ; In case of an external area: EXTERNAL_WATER_LEVEL of the breach.
= The calculated water flow which takes place at time , based on the weir formula.
= The timestep in seconds.
= The limited volume of water transported, based on the remaining water.
= The size of in the entry area; In case of an external area: EXTERNAL_AREA of the breach.
= The (optional) size of the bottom of the entry area; In case of an external area: the second value provided for the EXTERNAL_AREA of the breach. In case only one value is provided,

Trapezoidal prism

External water body as a trapezoid

The external water area can also be represented as a trapezoidal prism. The new water level then has to be calculated differently.

The length and width of the surface area are currently considered equal. Therefore:

The width of the bottom area is then calculated as:

Next, we will calculate the new water level based on :

if , determine the new water level based on a box shape:

otherwise, determine the new water level based on the determinant and flow direction:

Where:

= Length of the trapezoidal prism
= The size of in the entry area; In case of an external area: EXTERNAL_AREA of the breach.
= Width of the original top (surface at ) of the trapezoid
= The (optional) size of the bottom of the entry area; In case of an external area: the second value provided for the EXTERNAL_AREA of the breach. In case only one value is provided, : = Datum height of the bottom of the area water flows from.
= Width of the bottom of the trapezoid
= The water level of the entry area at time ; In case of an external area: EXTERNAL_WATER_LEVEL of the breach.
Datum height of the bottom of the area water flows from.
= The initial relative water height in the external area.
= Slope of the trapezoid
= Area of the trapezoidal slide of the trapezoid prism at time
= Constant based on new area of the trapezoid
= Polynomial discriminant
= The limited volume of water transported, based on the remaining water.


Notes

  • The water level we cannot become lower than the surface height se, defined by EXTERNAL_SURFACE_LEVEL in case of an external entry area. In case of an input area within the project, the surface height is determined per grid cell.
  • When water is flowing from the project area to an external area, the slope of the external area remains the same, even when the water level rises above the original water level. In the case of a trapezoidal prism, is extrapolated upwards.

Related

The following topics are related to this formula.

Features
Breach
Formulas
Weir formula
Breach growth formula
Models
Surface model