TUFLOW 2D Hydraulic Structures: Difference between revisions
Content deleted Content added
No edit summary |
|||
Line 4:
*<u>[https://www.tuflow.com/library/webinars/#structures Webinar Link: Modelling Energy Losses at Structures]</u><br>
*<u>[https://www.tuflow.com/library/webinars/#nov2022_hydraulic_modelling_bridge Webinar Link: 1D, 2D & 3D Hydraulic Modelling of Bridges]</u>
<br>
= 2D Bridge Modelling in TUFLOW - Overview =
The TUFLOW 2D solution explicitly predicts the majority of “macro” losses due to the expansion and contraction of water through a constriction, or around a bend, provided the resolution of the grid is sufficiently fine (<u>[https://www.tuflow.com/Download/Publications/Flow%20Through%20an%20Abrupt%20Constriction%20-%202D%20Hydrodynamic%20Performance%20and%20Influence%20of%20Spatial%20Resolution,%20Barton,%202001.pdf Barton, 2001]; [https://www.tuflow.com/Download/Publications/Modelling%20of%20Bends%20and%20Hydraulic%20Structures%20in%20a%202D%20Scheme,%20Syme,%202001.pdf Syme, 2001]; [https://www.tuflow.com/Download/Technical_Memos/Modelling%20Bridge%20Piers%20in%202D%20using%20TUFLOW.pdf Ryan, 2013]</u>). Where the 2D model is not of fine enough resolution to simulate the “micro” losses (e.g. from bridge piers, vena contracta, losses in the vertical (3rd) dimension), additional form loss coefficients and/or modifications to the cells widths and flow height need to be added.
==Contraction/Expansion Losses (“Macro” Losses)==
Loss of energy is caused by the flow contraction during the expansion of water after the vena-contracta inside a bridge section and the flow expansion downstream a bridge. As discussed above, this type of "macro" losses can be explicitly resolved by the TUFLOW 2D solver, provided that a proper turbulence model and mesh size are used
<br>
[[File:FC_Velocity_Example.PNG|600px]] <br>
Line 26 ⟶ 27:
==Bridge Deck and Rail (Super Structure)==
When a bridge deck become partially or completely submerged, the deck could generate extra afflux resulting in increased water levels and flood extents upstream of the structure. The flow around the deck is highly 3-dimentional and complexed due to the different deck designs/profiles and/or the occurrence of pressure flow. In 2D SWE solver, depth-varying form loss values are often needed to reproduce the afflux caused by such structure. Due to the complexity of the flow, guidelines on how to set the form loss coefficient for the bridge deck are rare. We have carried out a joint research with QLD TMR (Queensland Department of Transport and Main Roads) regarding how to choose a proper form loss value for the bridge deck
<br>
[[File:CFD_study.png|600px]]
Line 35 ⟶ 36:
[[File:FormLoss_vs_TWT.png|600px]]
<br>
=== Bridge Design (hB/T) vs Form Loss Coefficient Table===
The peak loss coefficient value is a function of the ratio of the depth underneath the deck (hB) and the thickness of the deck (T). This table can be used to estimate the deck form loss coefficient based on the bridge design (hB/T).
<ol>
{| style="text-align: center;" class="wikitable" width="35%"
Line 49 ⟶ 52:
</ol>
*For bridge with more complicated designs (e.g. girders), higher form loss might be required due to the higher surface roughness of the bridge.
*If the hB/T ratio is less than 2 or greater than 6, use a peak form loss coefficient of 0.42 (minimum) or 0.20 (maximum), respectively.
'''NOTE''': This form loss value should not be confused with the value of 1.56 used in the pressure flow approached adopted in <u>[[1D_Bridges | TUFLOW 1D "B" and "BB" bridge]]</u>. TUFLOW 1D bridge pressure flow approach is based on the section 4.13.2 "All Girders in Contact with Flow (Case II)" of ''Guide to Bridge Technology Part 8: Hydraulic Design of Waterway Structures'' (AUSTROADS, 2019). The original hydraulic experiment conducted by <u>[https://hdl.handle.net/10217/39009 Liu et al (1957)]</u> in a laboratory flume with a pair of bridge abutments and a deck. The flow conditions were similar to orifice flow due to the high blockage ratio caused by the abutments and the deck. When modelling bridges in 2D, the contraction/expansion losses caused by the abutments would be handled explicitly by the 2D solver, so a value 1.56 can lead to duplication of the contraction/expansion losses caused by the bridge abutments.<br>
Line 55 ⟶ 60:
=TUFLOW 2D Bridge Setup=
There are two methods available to model depth varying form loss of a bridge structure:
* <u>[[TUFLOW_2D_Hydraulic_Structures#2D_Layered_Flow_Constriction|2D Layered Flow Constriction (2d_lfcsh)]]</u>
:The traditional method used to model depth-varying form loss through bridge components such as piers, decks, and rails.
*<u>[[TUFLOW_2D_Hydraulic_Structures#2D_BG_Shape |2D BG Shape (2d_bg)]]</u> (introduced in the 2023 release)
:A simplified approach developed to simplify the model input based on the findings from the joint TMR Study <u>[https://tuflow.com/media/7554/2022-bridge-deck-afflux-modelling-benchmarking-of-cfd-and-swe-codes-to-real-world-data-collecutt-et-al-hwrs.pdf (Collecutt et al, 2022)]</u>.
Both methods provide options for representing flow surcharging, the pressure flow of bridge decks and eventually submerged bridge flow at higher water levels. During the surcharging of bridge decks, higher energy losses can be specified to simulate the pressure flow.
Examples for how to configure both approaches are provided in the 2D structures section of the <u>[[TUFLOW_Example_Models#2D_Structures |TUFLOW Wiki Example Models]]</u> and <u>[[Tutorial_M04 |Tutorial Module 4]]</u> - 2D Bridges.
==2D Layered Flow Constriction (2d_lfcsh)==
Four flow constriction layers are represented in a 2d_lfcsh layer. The lower three layers represents the pier, the bridge deck and the rails. Each layer has its own attributes to specify the blockage and the form loss coefficient. The top (fourth) layer assumes the flow is unimpeded, representative of flow over the top of a bridge. Within the same shape, the invert of the bed, and thickness of each layer can vary in 3D.
The following table provides an overview for how to determine the blockage and form loss coefficient for each layer:<br>
{| style="text-align: left; margin-left: 0; " class="wikitable" width="80%"
!colspan="1" style="background-color:#005581; font-weight:bold; color:white;"| Layer
!colspan="1" style="background-color:#005581; font-weight:bold; color:white;"| Description
!colspan="1" style="background-color:#005581; font-weight:bold; color:white;"| Blockage (%)
!colspan="1" style="background-color:#005581; font-weight:bold; color:white;"| Form Loss Coefficient (FLC)
!colspan="1" style="background-color:#005581; font-weight:bold; color:white;"| Notes
|-
| 1 || Pier layer || ~5% (can be omitted if included in FLC) || Estimate using <u>[[TUFLOW_2D_Hydraulic_Structures#Pier_Losses | Pier Losses]]</u> || Represents flow obstruction from piers beneath the bridge deck
|-
| 2 || Bridge deck || 100% || Use calibration data, if available, to determine FLC. <br> If no calibration is available, estimate using <u>[[TUFLOW_2D_Hydraulic_Structures#Bridge_Design_.28hB.2FT.29_vs_Form_Loss_Coefficient_Table | hB/T vs FLC]]</u> table || Full blockage, no flow through the deck
|-
| 3 || Bridge rails || 10% – 100% || Use calibration data, if available, to determine FLC. <br>
If no calibration data is available, combined FLC for Layers 2 and 3 should be estimated using the <u>[[TUFLOW_2D_Hydraulic_Structures#Bridge_Design_.28hB.2FT.29_vs_Form_Loss_Coefficient_Table | hB/T vs FLC]]</u> table, where T = L2_Depth + (pBlockage × L3_Depth)
*(pBlockage × L3_Depth) represents the solid portion of the rails
*L2 FLC and L3 FLC should sum to the combined FLC
|Blockage and FLC depends on rail type <br> Sensitivity testing with 100% blockage is recommended due to potential for debris during flood
|-
| 4 || Above rails || 0% || 0 || Represents unimpeded overtopping flow
|}
<ol>
[[File:2d_lfcsh_attributes.png | 500px ]]
</ol>
<br>
===Blockage===
The 2d_lfcsh functions by adjusting the flow width and the form loss of 2D cell faces. The combined blockage across the 4 layers is calculated at each simulation timesteps:
Line 78 ⟶ 112:
</ol>
===Form Loss Approach===
The combined form loss coefficient is determined using one of three methods. The form loss coefficient method can be specified either individually using the 2d_lfcsh “Shape_Options” attribute or globally using the .tcf command: <br>
<tt><font color="blue">Layered FLC Default Approach</font> <font color="red">==</font> [ METHOD A | {METHOD B} | METHOD C | METHOD D]</tt> <br>
:<b>METHOD A</b>: The losses are accumulated as the water level rises through the layers. <br>
<ol>
[[File:Eq_flc_cumulate.png |450px]]
</ol>
:*Layer 1: Constant form loss (L1_FLC)
:*Layers 2 & 3: Form loss increases based on the depth of water in layer 2 & 3; peak form loss at top of Layer 3
:*Above Layer 3: Applies the full accumulated form loss continuously, even when overtopping begins (no reduction)
:Note: Simpler method but tends to overestimate losses when the structure is submerged or overtopped<br>
:<b>METHOD B</b> (default): the losses are applied pro-rata according to the depth of water in each layer. <br>
<ol>
[[File:Eq_flc_portion.png |430px]]
</ol>
:*Layer 1: Constant form loss (L1_FLC)
:*Layers 2 & 3: Form loss increases based on the depth of water in layer 2 & 3; peak form loss at top of Layer 3
:*Above Layer 3: Total form loss gradually reduces as water overtops the structure
:Note: Maintains backward compatibility but may underrepresent losses during pressurised or overtopped flows<br>
:<b>METHOD C</b> (recommended): hybrid approach combining Method A and Method B. <br>
<ol>
[[File:Eq_flc_methodC.png |520px]]
</ol>
:*Layer 1: Constant form loss (L1_FLC)
:*Layers 2 & 3: Gradual increase in form loss with water level, following Method A
:*Above Layer 3: Total form loss gradually reduces as water overtops the structure, following Method B
:Note: Recommended method; aligns closest to CFD modelling results and TUFLOW HPC behaviour.
:<b>METHOD D</b>: Allows the modeller to control the depth at which the losses start to reduce when the flow transitions between pressure flow and drowned flow.
:This approach is the same used by the 2d_bg layer (introduced in the 2023-03 release). It is recommended to use the 2d_bg layer as it has the benefit of a simplified attribute table, for easier user input.
===Form Loss Calibration Example - Iowa River Flood Study===
In this study, a combined form loss coefficient of 0.35 was used to match observed head loss during slight overtopping of a bridge. The FLC values for each layer were adjusted to achieve the correct combined form loss. The table and plot show how each layer contributes to the total form loss and highlight the differences in calculated form loss between the three methods.
{| style="text-align: left; margin-left: 0;" class="wikitable" width="60%"
!colspan="1" rowspan="2" style="background-color:#005581; font-weight:bold; color:white;" width=6%| Layer
!colspan="1" rowspan="2" style="background-color:#005581; font-weight:bold; color:white;" width=10%| Depth (m)
!colspan="1" rowspan="2" style="background-color:#005581; font-weight:bold; color:white;" width=12%| Blockage (%)
!colspan="2" style="background-color:#005581; font-weight:bold; color:white;" width=20%| Method A
!colspan="2" style="background-color:#005581; font-weight:bold; color:white;" width=20%| Method B
!colspan="2" style="background-color:#005581; font-weight:bold; color:white;" width=20%| Method C
|-
! style="background-color:#005581; font-weight:bold; color:white;"
! style="background-color:#005581; font-weight:bold; color:white;"
! style="background-color:#005581; font-weight:bold; color:white;"
! style="background-color:#005581; font-weight:bold; color:white;"| Combined FLC
! style="background-color:#005581; font-weight:bold; color:white;"| Layer FLC
! style="background-color:#005581; font-weight:bold; color:white;"| Combined FLC
|-
| 1 || 5.0 || 5 || 0.07 || 0.07 || 0.07 || 0.07 || 0.07 || 0.07
|-
| 2 || 1.5 || 100 || 0.15 || 0.22 || 1.05 || 0.30 || 0.15 || 0.22
|-
| 3 || 1.0 || 50 || 0.13 || 0.35 || 0.70 || 0.35 || 0.13 || 0.35
|}
</ol>
<br>
<ol>
[[File:
</ol>
==2D BG Shape (2d_bg)==
2D BG Shape is similar to the Layered Flow Constriction, but has several
The following table provides an overview of how to determine the blockage and form loss coefficient for each layer:<br>
{| style="text-align: left; margin-left: 0; " class="wikitable" width="80%"
!colspan="1" style="background-color:#005581; font-weight:bold; color:white;"| Layer
!colspan="1" style="background-color:#005581; font-weight:bold; color:white;"| Description
!colspan="1" style="background-color:#005581; font-weight:bold; color:white;"| Blockage (%)
!colspan="1" style="background-color:#005581; font-weight:bold; color:white;"| Form Loss Coefficient (FLC)
!colspan="1" style="background-color:#005581; font-weight:bold; color:white;"| Notes
|-
| 1 || Pier layer || ~5% (can be omitted if included in FLC) || Estimate using <u>[[TUFLOW_2D_Hydraulic_Structures#Pier_Losses | Pier Losses]]</u> || Represents flow obstruction from piers beneath the bridge deck
|-
| 2 || Bridge deck || 100% || rowspan="2" | The Super Structure (Super_S) is the bridge deck and rails layers combined. <br>
Use calibration data, if available, to determine FLC. <br>
If no calibration data is available, the Super_S FLC should be estimated using the <u>[[TUFLOW_2D_Hydraulic_Structures#Bridge_Design_.28hB.2FT.29_vs_Form_Loss_Coefficient_Table | hB/T vs FLC]]</u> table, where T = Deck_Depth + (Rail_pBlockage*Rail_Depth)
*(Rail_pBlockage*Rail_Depth) represents the solid portion of the rails
|| Full blockage, no flow through the deck
|-
| 3 || Bridge rails || 10% – 100% || Sensitivity testing with 100% blockage is recommended due to potential for debris during flood events
|-
| 4 || Above rails || 0% || 0 || Represents unimpeded overtopping flow
|}
<ol>
[[File:2d_bg_attributes.png | 700px ]]
</ol>
===Inflection Point===
Based on findings from the joint study <u>[https://tuflow.com/media/7554/2022-bridge-deck-afflux-modelling-benchmarking-of-cfd-and-swe-codes-to-real-world-data-collecutt-et-al-hwrs.pdf (Collecutt et al, 2022)]</u>, the head loss peaks when the water level is approximately 1.6*T above the bridge soffit, and decays slowly as the bridge becomes progressively drowned out. The 'SuperS_IPf' attribute (inflection point factor, default = 1.6) can be used to define the height of the inflection point. The solid portion of the rail layer is also added to the deck thickness to calculate the depth to the inflection point (D<sub>IP</sub>), i.e.:
<ol>
[[File:eq_flc_bg_infection_point.png | 520px ]]
</ol>
===Form Loss Approach===
The form loss approach is similar to the FLC approach METHOD C, with L2/L3 replaced by a single super structure layer:
<ol>
Line 136 ⟶ 227:
</ol>
===Form Loss Calibration Example - Iowa River Flood Study===
This example uses the same bridge setup described in the<u>[[TUFLOW_2D_Hydraulic_Structures#Form_Loss_Calibration_Example_-_Iowa_River_Flood_Study | 2D Layered Flow Constriction]]</u> section, with the following parameters applied:
*SuperS_FLC = 0.28
*SuperS_Ipf = 1.6,
The Depth to Inflection Point (D<sub>IP</sub>) is calculated as 3.2m above the bridge soffit.
The table and figure below show how the form loss value varies with water depth.
<ol>
{| style="text-align: center;" class="wikitable" width="32%"
Line 146 ⟶ 244:
| Pier || 5.0 || 5 || 0.07
|-
| Deck || 1.5 || 100 || rowspan="2" | 0.28
|-
| Rail || 1.0 || 50
|}
Line 155 ⟶ 253:
== 2D Bridges Line vs Polygon Layer ==
The form loss coefficient (FLC) is applied differently when using a line compared to a polygon for both 2d_lfcsh and 2d_bg inputs. The FLC is applied at cell sides (u and v faces) as this is where velocities are calculated. <br>
For larger bridges that spread across multiple cells, it is recommended to use a polygon layer, which selects all u and v faces falling within the polygon.
<br>
<
'''2D Layered Flow Constriction (2d_lfcsh)'''
{| style="text-align: left; margin-left: 0;" class="wikitable" width="80%"
!colspan="1" style="background-color:#005581; font-weight:bold; color:white; width: 7.5%;"| Geometry
!colspan="1" style="background-color:#005581; font-weight:bold; color:white; width: 7.5%;"| Line Type
!colspan="1" style="background-color:#005581; font-weight:bold; color:white; width: 11%;"| Width Attribute
!colspan="1" style="background-color:#005581; font-weight:bold; color:white; width: 22%;"| FLC Input
!colspan="1" style="background-color:#005581; font-weight:bold; color:white; width: 22%;"| FLC Applied
!colspan="1" style="background-color:#005581; font-weight:bold; color:white; width: 30%;"| Notes
|-
!rowspan="3" | Line
| Thin
| zero
| Total form loss of the bridge
| Applies the FLC to a single row of cell sides
| This approach is cell size independent. It is the easiest setup and the preferred / recommended approach when using 2d_lfcsh.
|-
| Thick
| between zero and 1.5 times the cell size
| Half of total form loss of the bridge
| FLC applied to all sides of selected cells
| A cell is selected if the polyline intersects the cell crosshair. Caution should be taken when using a "thick" line, as changes in cell size can cause it to become a "wide" line. If this occurs, the FLC attribute should be recalculated to avoid overestimating losses.
|-
| Wide
| larger than 1.5 times the cell size
| Total form loss divided by number of cell sides in the direction of flow <br>
''(number of cell sides in the direction of flow is calculated as line width divided by cell size)''
| FLC applied to all sides of selected cells
| Caution should be taken when using a "wide" line. Changes in cell size may require recalculating losses.
|-
!rowspan="1" | Polygon
| -
| -
| Total loss per unit length (meters or feet) in the direction of flow
| FLC * cell size applied to all sides of selected cells
|
|}
'''2D Bridge (2d_bg)'''
{| style="text-align: left; margin-left: 0;" class="wikitable" width="80%"
!colspan="1" style="background-color:#005581; font-weight:bold; color:white; width: 7.5%;"| Geometry
!colspan="1" style="background-color:#005581; font-weight:bold; color:white; width: 7.5%;"| Line Type
!colspan="1" style="background-color:#005581; font-weight:bold; color:white; width: 11%;"| Width Attribute
!colspan="1" style="background-color:#005581; font-weight:bold; color:white; width: 22%;"| FLC Input
!colspan="1" style="background-color:#005581; font-weight:bold; color:white; width: 22%;"| FLC Applied
!colspan="1" style="background-color:#005581; font-weight:bold; color:white; width: 30%;"| Notes
|-
!rowspan="3" | Line
| Thin
| zero
| Total form loss of the bridge
| Applies the FLC to a single row of cell sides.
| This approach is cell size independent.
|-
| Thick
| larger than zero
| Total form loss of the bridge
| FLC/2 applied to all sides of the selected cells
| This approach is cell size independent. A cell is selected if the polyline intersects the cell crosshair.
|-
| Wide
| Not supported
| –
| –
| BG polygon shapes are recommended if more than 3 rows of faces must be selected.
|-
!rowspan="1" |Polygon
| -
| ''(used to automatically distribute the total FLC to the selected faces)''
| Total form loss of the bridge
| FLC / Deck_Width * cell size applied to all sides of selected cells
| For bridges modelled using a 2d_bg polygon the relative ratio of the bridge width to the 2D cell size should be 4 or greater. For more information on this see <u>[https://downloads.tuflow.com/Other/2d_bg_R_Bridge_Configuration_Advice_202503.pdf 2d_bg_R_Bridge_Configuration_Advice.pdf]</u>.
|}
<br>
The following diagrams demonstrate how the input FLC is applied for the four geometry options for 2d_lfcsh and 2d_bg layers: <br>
[[File:2dlfcsh 2dbg combined.png|1200px]]
It is good modelling practice to check the <u>[[Check_Files_2d_lfcsh_uvpt | lfcsh_uvpt_check]]</u> and <u>[[Check Files 2d bg uvpt check | bg_uvpt_check]]</u> files to confirm the number of faces selected and the FLC values assigned. It is also strongly recommended to undertake a sensitivity analysis on the applied form losses in the model to check if it makes any difference to the results and/or double check against other methods (hand calculations, other software, CFD modelling), especially if the bridge is near an area of interest. If calibration data is available, this should be used to guide the form loss value specification.<br>
<br>
| |||