TUFLOW 1D Channels and Hydraulic Structures: Difference between revisions
Content deleted Content added
| (78 intermediate revisions by 7 users not shown) | |||
Line 1:
=Introduction=
The objective of the following pages is to supplement the <u>[https://docs.tuflow.com/classic-hpc/manual/latest/ TUFLOW
=1D Channel Types=
1d_nwk channels represent open channels, hydraulic structures, operational structures and other flow controls. A channel is either digitised as a line or point with the relevant hydraulic properties entered into the appropriate GIS attributes, details for this can be found within the links below or the <u>[https://docs.tuflow.com/classic-hpc/manual/latest/ TUFLOW
==Open Channels==
{| class="wikitable" width="50%"
Line 23:
| [[1D_Bridges | Bridges]]|| Loss theory, irregular shaped bridges and common check files used.
|-
| [[1D_Weirs | Weirs]]||
|-
| [[ 1D_Pumps | Pumps]]|| Pump attributes, 2D-2D & 1D-2D configurations, Estry setup, depth-discharge database and 1D results file.
Line 33 ⟶ 31:
| [[1D_Manholes | Manholes]]|| Losses and storage chambers.
|-
|[[1D Syphons|Syphons]]
|1D configuration, losses.
|-
|}
Line 61 ⟶ 48:
== How can I read approach-capture flow curves for on-grade pits into TUFLOW? ==
TUFLOW currently doesn't support the use of approach-capture flow curves directly.
The pit depth-flow curve relationship can be derived as follows:
*Select a few representative slope increments for the project area and conduct a Manning's equation calculation for each slope to calculate the gutter flow rates for incrementally increasing flow depths.
*Multiply the calculated flows by the % value in the on grade curve. This turns the % capture vs. flow curve into a depth vs. flow relationship.
*Repeat for different gutter shapes
== Should I model culverts in the 1D or 2D domain? ==
Line 71 ⟶ 60:
* Benefits:
** Good stability and improved hydraulic behaviour at the culvert entrance/exit.
** The expansion loss (exit loss) can be explicitly handled in 2D provided the 2D cell resolution is sufficiently fine to model the expansion of flow downstream of the culvert.
* Challenges:
** The contraction loss (entry loss), which is related to the expansion of water after the
** Side and soffit wall friction are not modelled unless Manning's
** The vertical walls not only create extra friction but also straightens the flow in the direction of the wall. Thin breaklines can be used to represent these walls in 2D, but it is likely to cause saw-tooth effect if sub-grid sampling (SGS) is not used, i.e. extra numerical head loss, if there are too few cells between the walls.
** Flow overtopping can be represented to some extent by assuming 100% blockage of layer 2 in a 2d_lfcsh, but the flow upstream and downstream before the culvert is overtopped is hard to model.
Line 81 ⟶ 70:
* Benefits:
** A more appropriate approach where the 2D cell size is greater than around half the total culvert width.
** Contraction losses (entry losses) are handled better.
** Flow overtopping can be modelled in 2D.
* Challenges:
** Expansion losses (exit losses) are very dependent on the 2D cell resolution. A 2D cell size much larger than the culvert width will not reproduce the expansion losses very well (even with SGS) and the culvert's exit loss needs to cover this. A finer 2D cell size (several or more cells across the culvert) will reproduce the expansion losses much better and the culvert's exit loss need to be reduced to compensate it. This usually only happens for large 1D culverts with high velocities as there can be losses duplicated in the 2D on the exit side as the 2D flow expands (i.e. duplication of the exit/expansion loss). The latest release has a new feature to automatically adjust 1D culvert losses based on the 2D approach/departure velocities as what happens with a 1D-1D-1D arrangement. See this <u>[https://downloads.tuflow.com/_archive/Publications/Modelling%20of%20Bends%20and%20Hydraulic%20Structures%20in%20a%202D%20Scheme,%20Syme,%202001.pdf paper]</u> for more information.
** HX connections may cause instability, especially with a skewed culvert outlet, but can
** Given that the SX connection applies the flow going out of the 1D culvert as a source term without momentum, it is difficult to completely prevent the water from piling up. If required, wingwalls can be modelled as thin breaklines to help guide the water away. Using the SX boundary Z flag lowers other SX cells below the 1D culvert invert level and it can mitigate the water from piling up.
== What entry/exit loss and contraction coefficients should I use for 1D culverts? ==
We don’t provide hard recommendations on the exit and entry losses to use for culverts as we have found different organisations around the world, typically government, have their own guidelines for different types of inlets configurations and require these to be used, for example, the <u>[https://www.ipwea-qnt.com/products-resources/qudm/ Queensland Urban Drainage Manual]</u> (QUDM). However, it is very important to understand how losses are applied and that different 1D solvers may treat them differently. For cross-checking your results from any hydraulic modelling software, a simple calculation applying the entry and exit losses (allowing for any automatic adjustments as discussed below) to the computed head (V2/2g), plus allowing any surface roughness losses (Manning's equation) for longer culverts, is the best practice for culverts flowing in a sub-critical flow condition (i.e. downstream controlled flow).<br>
For the entrance loss values, the approach should be to use values as quoted in the literature or guidelines for the inlet shape and design unless there is evidence to use another value (e.g. comparison with reliable calibration data would indicate different energy losses).
Line 97 ⟶ 86:
Occasionally there are situations where non-standard entrance and exit loss values are needed. A good example is if the approach or departure flow is skewed to the culvert direction. In these situations there may also be a significant bend (energy) loss occurring as the water changes direction entering or leaving the structure. To account for this the modeller may need to increase the entrance and/or exit loss values.<br>
By default, TUFLOW adjusts the entrance and exit losses of 1D structures flowing sub-critical every timestep based on the approach/departure velocities for 1D-1D-1D and with a new feature in the 2020-10-AA build for 2D-1D-2D. The entrance losses are adjusted based on an empirical relationship from flume testing whilst the exit loss equation is theoretically derived as mentioned above -
TUFLOW, by default, allows the losses to reduce to effectively a zero loss coefficient (i.e. 0.0001). A zero loss occurs where the approach and departure velocity is the same as the structure velocity. For example, a clear-spanning bridge over a concrete lined channel with the water level below the bridge deck will experience no energy losses until the bridge deck is surcharged so if your software is applying unadjusted or limited energy loss coefficients there will be an unrealistic energy loss at the structure for flow below the bridge deck. For culverts, in most cases there will be some losses as it is rare that the channel is of identical shape and slope to the culvert with usually the culvert being more constrictive and therefore a higher velocity so the adjusted coefficients are nearly always non-zero. At the other extreme is flow from or into a near still body of water (e.g. a lake or the ocean). In this situation the loss coefficient(s) will not be reduced and the maximum energy loss possible should occur.
If the default adjust losses approach is used (<font color="blue"><tt>Structure Losses</tt></font> <font color="red"><tt>==</tt></font><tt> ADJUST</tt>) the recommendations are to use industry guidelines for the entrance loss coefficient based on the shape/design of the inlet (these coefficients are typically based on a near zero approach velocity), and to use 1.0 for the exit loss coefficient. This applies to 1D culverts connected to 1D channels. The adjusted entrance and exit losses can be viewed over time in the _TSL layer,
Since the 2020-10-AA release, TUFLOW has a new (beta functionality) to have the losses automatically adjusted for linked 1D culverts and other structures connected to 2D domains through SX links (<font color="blue"><tt>Structure Losses SX</tt></font> <font color="red"><tt>==</tt></font><tt> ADJUST</tt>), see Section 6.2 from the <u>[https://downloads.tuflow.com/TUFLOW/Releases/2020-10/TUFLOW%20Release%20Notes.2020-10-AA.pdf 2020-10 Release Notes]</u>
The other values to consider for modelling culverts are the inlet contraction coefficients used when the flow is upstream controlled flow should this occur. Typically the TUFLOW default values for these values should be used unless the inlet shape and design indicates otherwise.<br>
Line 152 ⟶ 141:
== How do negative "M", "N" and "R" values deactivate 1D cross-sections? ==
When negative “M” values are used to deactivate 1D cross-sections, it changes the flow area, wetted perimeter, hydraulic radius, and the conveyance, while keeping the storage width unchanged. It is like inserting a thin plate at the middle of the channel, as shown in the figure below. This approach is preferred to retain the storage in the system but reduce the conveyance (e.g., off-stream storage not contributing to conveyance or sections either side of a constriction to generate better approach/departure velocities if using the Total Area approach). <br>
'''Note:''' The negative "M" value must be a negative of a valid material ID.<br>
<br>
<ol>
[[File:Section deactivation M.png | 200px]]<br>
</ol>
When using negative (-1) “N” (Manning’s n) or “R” (relative resistance factor) to deactivate the 1D cross-sections, it reduces both the flow area and the storage width, as shown in the figure below.<br>
<br>
<ol>
[[File:Section deactivation N R.png | 200px]]
</ol>
==Calculating a 1d_bc head-discharge curve==
As of the 2023-03-AD release, it is possible for TUFLOW to automatically create a Head-Discharge curve for 1d_bc HQ boundaries. For more information see Section 3.5 of the <u>[https://docs.tuflow.com/classic-hpc/release/2023-03-AE/OneDSolv-1.html#Auto1DHQ-2 2023-03 Release Notes]</u>.
If using older releases, this can be calculated manually as the condition to be satisfied for a normal depth is the Manning’s equation. The HQ depth should be set up with elevation and the discharge calculated from Manning’s equation (including Area and Hydraulic Radius for each elevation point) and can be calculated in a few steps:<br>
<ol>
<li>Set the model up your ESTRY model using a dummy downstream level boundary (e.g. a HQ boundary with no name, an error will be reported however the required check file will have been produced, otherwise any boundary will work) and run the model in test model (-t) so that check files are produced.
<li>From the check folder, open the '''*_1d_ta_tables.csv''' file.
<li>Find the entry for the most downstream cross-section (directly upstream of the 1d_bc HQ point). The '''*_1d_ta_tables.csv''' will contain a table of hydraulic parameters. The K entry highlighted is the cross-section conveyance in m<sup>3</sup>/s and is calculated as: <br>
K = 1/n * A * R<sup>2/3</sup><br>
Where
* A is the Eff Area
* R is the Radius
<br>
[[File:Ta tables example.png | 800px ]]<br><br>
<li>Determine a representative slope (S), this could be taken from the last two cross-sections.
<li>Calculate the downstream Discharge (this is the same as calculated discharge from Mannings equation): <br>
Q = K * SQRT(S)<br>
Where:
* S is a representative slope value for the downstream reach.
* K is obtained from the '''1d_ta_tables.csv'''
<li>Use the Elevation (Head) vs Discharge pairs to populate the HQ table.
</ol>
==Why is the maximum 1D water level higher than 2D water level?==
TUFLOW calculates storage at 1D nodes (including manholes) and calculates flux at pipe mid-sections. The storages at the manholes are saved in the format of ‘elevation vs nodal area’ tables, which can be reviewed in the .eof file:<br>
<ol>
[[File:1DpitHmax 1.png]] <br>
</ol>
The change in water volume in a manhole is calculated based on the inflow from the pit above (if there is one) and the inflow/outflow from the connected pipes. The change in water level is then calculated based on the change in volume and the nodal area at that elevation.<br>
<ol>
[[File:1DpitHmax 2.png | 300px ]] <br>
</ol>
TUFLOW extends the ‘elevation vs nodal area’ table by 5 m (the last 2 numbers in the table above), which allows for the calculation of pressurised pipe flow when a manhole is drowned. If the maximum 1D water level is higher than the DEM, that means the manhole was drowned and pressurised during the peak of the flood. <br>
<ol>
[[File:1DpitHmax 3.png | 300px ]] <br>
</ol>
==How to represent the hydraulic behavior at the outlet of a pipe, including headwalls and transitions to open channels? ==
In general, pipe outlet headwalls affect the exit loss at the pipe outlet, and this can be modeled by changing the 1d_nwk “ExitC_or_WSb” attribute - exit loss coefficient for outlet controlled flow. The default exit loss coefficient is 1.0, assuming there is no headwall. For more details on these factors, refer to the <u>[https://docs.tuflow.com/classic-hpc/manual/latest/ TUFLOW Manual]</u>.
More guidance on applying entry/loss and contraction coefficients is available [[TUFLOW_1D_Channels_and_Hydraulic_Structures#What_entry.2Fexit_loss_and_contraction_coefficients_should_I_use_for_1D_culverts.3F | here]].
Once the water is in the 2D domain, headwalls can also affect the local flow field. The modelled hydraulic behavior depends on the model cell size. If the model has grid size smaller than the length of the headwall, thin breaklines or 2d_zsh tin can be added to represent the geometry in the 2D model. If the grid size is larger than the headwall, it only calculates ‘cell-averaged’ velocity including the impact of the headwall.
==Calculating velocity and bed shear stress at a pipe outlet==
The velocity and bed shear stress reported in the 2D domain at a pipe outlet depends on the model cell size. If the cell size is small enough and any outlets structures (i.e. wingwalls) are properly represented by the 2d cell/face elevations, TUFLOW can reasonably represent the flow field at the pipe outlet and estimate the bed shear stress base on the manning’s equation. However, be aware that TUFLOW solves 2D depth-averaged shallow water equation, so there is limitation of using it to estimate bed shear stress around complex structures or geometry that could experience vertical acceleration or re-circulation.
== How do I model blockage for irregular shaped culverts or channels? ==
For irregular shaped culverts and channels (1d_nwk "I" type), it is possible to model blockage by using the Skew (in degrees) attribute in the 1d_xs layer. This reduces the effective flow width of the cross section by rotating it relative to the direction of flow and can be used to approximate reduced capacity.
A skew angle of 0 means the cross section is aligned with the flow. As the angle increases, the effective width is reduced. For example, a skew of 60 degrees results results in a 50% reduction (as cos(60) = 0.5). While this does reduce the effective width of the structure, it isn’t expected to see the same reduction to conveyance as conveyance depends on the Manning’s n value, flow area and hydraulic radius.
== Do calculations differ if I specify multiple culverts using the Number_of attribute, even if the total dimensions are the same as a single culvert? ==
Yes. This can slightly affect the modelling results.
The key difference between a single wide culvert and multiple culverts with the same total width is the wetted perimeter. Multiple culverts have a slightly larger wetted perimeter due to the internal walls, which increases friction in downstream-controlled regimes. While friction losses are generally small compared to entry and exit losses, the impact on hydraulic behaviour is usually minor.
However, it’s important to model the configuration that best represents the actual structure when assessing culvert performance.
== How do I model a scruffy dome? ==
A scruffy dome is a type of dome inlet screen installed over a drop inlet, commonly used in stormwater systems. These inlets can be modelled in TUFLOW using a 1D Q type pit with a depth–discharge curve that reflects the hydraulic performance of the dome.
[https://www.ipwea-qnt.com/products-resources/qudm/ <u>QUDM</u>] recommends using the lesser of the weir flow or orifice flow equations to define inlet capacity (see Section 7.5.4). A blockage factor is also recommended to account for potential debris build-up during flood events (see Table 7.5.1 Provision for blockage at kerb inlets). This factor may vary depending on local guidelines, so regional standards should be consulted.
Once the rating curve is calculated, it can be added to the pit inlet database and linked to the Q type pit.
<br>
<ol>
[[File:How_do_I_model_a_scruffy_dome_in_TUFLOW.png|300x300px]]<br>
</ol>
== How can trash screens be represented? ==
A trash screen is a physical barrier, typically made from bars or mesh, placed at the inlet or outlet of a culvert, pipe, or channel to prevent debris from entering or exiting the structure.
Trash screens can be represented by applying additional head losses and blockage effects in accordance with <u>[https://www.ipwea-qnt.com/products-resources/qudm/ QUDM]</u> guidance. The approach depends on the location of the screen and the degree of blockage.
Possible use cases include:
* Where the screen is attached directly to the culvert inlet or outlet, the losses are applied to the culvert itself.
* Where the screen is located upstream of a culvert inlet or downstream of a culvert outlet, it can be modelled as a short zero length CF, IF, or RF channel before (upstream) or after (downstream) the culvert.
* In all cases, the net flow area should be reduced to account for the bars and any debris blockage, either by adjusting the flow area directly or by using the pBlockage attribute.
For example, the image below could be represented as a rectangular shaped screen located upstream or downstream of a culvert, modelled as a short zero length RF (rectangular culvert with fixed losses) channel using a 1d_nwk GIS layer.
<ol>
[[File:Trash_screen_schematic.png|alt=|300x300px]]
</ol>
== How can fishway baffles in a culvert be modelled? ==
<span style="color:red; font-weight:bold;">UNDER CONSTRUCTION</span>
Fishway baffles are barriers inside a culvert that slow the water and help fish swim through.
These approaches could be used to represent them inside a model:
* Short culvert with smaller opening.
** A short culvert section is added with its bottom and top levels set to match the fish baffle. The width is set to the total width of the fish passage gaps. This represents the reduced space for water to pass through.
* Two culverts in a row.
** Two culverts are connected in series. The first is the normal culvert, and the second is a very short culvert with a smaller opening sized to the fish baffle gaps.
For both options, sensitivity testing is recommended because at low flows the way the culvert is defined can change the results. These methods do not include turbulence or detailed 3D water movement around the baffles. If accurate fish passage flows are required, a CFD 3D model can be used to determine the flows and those flows can then be applied using a 1D Q type channel.
<ol>
[[File:Rectangular culvert baffled.png|alt=|300x300px]]
</ol>
<br>
{{Tips Navigation
|uplink=[[ TUFLOW_Modelling_Guidance | Back to TUFLOW Modelling Guidance]]
| |||