Difference between revisions of "Direct Rainfall (Rain on Grid) Modelling Guidance"

From Tuflow
Jump to navigation Jump to search
 
(30 intermediate revisions by 4 users not shown)
Line 1: Line 1:
 
= Direct Rainfall =
 
= Direct Rainfall =
 
Please see <u>[https://www.tuflow.com/library/webinars/#feb2021_direct_rainfall Direct Rainfall (Rain on Grid) Webinar]</u>.<br>
 
Please see <u>[https://www.tuflow.com/library/webinars/#feb2021_direct_rainfall Direct Rainfall (Rain on Grid) Webinar]</u>.<br>
<br>
+
 
 +
= TUFLOW Executable =
 +
* The single precision version of TUFLOW is recommended for direct rainfall models using HPC.
 +
* The double precision version of TUFLOW is recommended for direct rainfall models using TUFLOW Classic. 
 +
 
 +
= Rainfall Boundary Conditions=
 +
The following links provide boundary condition guidance for applying direct rainfall:
 +
* [[TUFLOW Rainfall Control File Examples]]
 +
* [[TUFLOW NetCDF Rainfall Format]]
  
 
= Rainfall Losses and Soil Infiltration=  
 
= Rainfall Losses and Soil Infiltration=  
Line 16: Line 24:
 
</ol>
 
</ol>
  
Similarly, a variety of soil infiltration options are supported. The available options include, Initial / Continuing Infiltration, the Horton Infiltration method and the Green-Ampt Infiltration method. The following link provides some future discussion on the Green-Ampt method.  
+
Similarly, a variety of soil infiltration options are supported. The available options include, Initial / Continuing Infiltration, the Horton Infiltration method and the Green-Ampt Infiltration method. The following link provides further discussion on the Green-Ampt method.  
 
* <u>[[Green-Ampt Infiltration Parameters | Green-Ampt Infiltration Parameters]]</u>
 
* <u>[[Green-Ampt Infiltration Parameters | Green-Ampt Infiltration Parameters]]</u>
 
  
 
= Common Questions Answered (FAQ)=
 
= Common Questions Answered (FAQ)=
Line 41: Line 48:
  
 
== What is the best approach for modelling buildings in rain on grid model? ==
 
== What is the best approach for modelling buildings in rain on grid model? ==
There are numerous different industry standard ways to represent buildings in a direct rainfall model. Australian Rainfall and Runoff Guideline, Project 15 (REPRESENTATION OF BUILDINGS IN 2D NUMERICAL FLOOD MODELS) discusses some of teh available options:
+
There are numerous different industry standard ways to represent buildings in a direct rainfall model. <u>[https://downloads.tuflow.com/_archive/Australian_Rainfall_Runoff_Project_15_Subproject_report_buildings_final.pdf Australian Rainfall and Runoff Guideline, Project 15 (Representation of Buildings in 2D Numerical Flood Models]</u> discusses some of the available options. Common TUFLOW modelling approaches are summarised below:
# Using depth varying Manning's n, apply a low value (n = 0.02) at shallow depth ( d < 0.03m) and a high value (>0.1  to the building area:
+
# Using depth-varying Manning's ''n'' over the area of the building footprint, with the application of a lower value (''n'' = 0.02) at shallow depths (d < 0.03m) and a higher value (''n'' > 0.3) at a more significant depths (d > 0.1m):
#* This is a very common and quick to setup approach that significantly slows down water passing through the buildings.
+
#* This is a very common and easy to implement option. The low Manning's ''n'' value aims to mimic the quick runoff response associated with drainage from the roof. The higher Manning's ''n'' value aims to represent the losses associated with deeper floodwater impacting the side of the building.
# Raise elevation with topographic updates:
+
# Raise the building footprint elevation using TUFLOW's topographic update features (eg. <font color="blue"><tt>Read GIS Z Shape</tt></font>):
#* Raising elevation eliminates water passing through the buildings, however as rainfall is assigned to all active cells, it will be falling on top of the buildings and from there to the ground.
+
#* Raising the model topography creates an obstruction to flow. It prevents floodwater from passing through buildings (as is the case with the Manning's ''n'' approach)
#* With bigger rainfall, this might create a waterfall like output around the buildings.
+
#* The application of rainfall on top of the building can however produce some undesirable results that require further consideration. Those being, water falling from the building roof to the ground can require a reduced model timestep to maintain model stability which can slow the simulation speed. Also, depending on the Map Cutoff Depth assumptions, water may be present in the results on the building rooftops. This may not be desired for mapping purposes. There are a variety of options available to resolve these issues:
# Exclude buildings from rainfall polygon:
+
:::* Retain the model design described above, with subsequent post-processing of the mapped results before reporting. Either, delete the model result where there is overlap with the building footprint, or overlay the building footprint polygon objects over/above the result dataset in the GIS Map Layout (hiding the flood model result within the building footprints).
#* This removes the rainfall from the model that supposed to fall on the buildings and underestimates the amount of rainfall entering the model.
+
:::* Exclude buildings from the rainfall polygon: This removes the rainfall from the model that would otherwise fall on the buildings. This approach will under-estimate the amount of rainfall entering the model. If the collective building footprint area is negligible in comparison to the entire model, this approach may be acceptable.
#* If the area of the buildings is negligible in comparison to the full model, it might be acceptable to leave the rainfall out for simplicity.
+
:::* Exclude buildings from the rainfall polygon and represent the rainfall that would be falling on the building using <font color="blue"><tt>Read GIS SA RF</tt></font> inflow boundaries. To do this, digitise a 2d_sa_rf polygon for each building (with a buffer of one of more 2D cells) where the building footprint has been excluded from the direct rainfall region. The 2d_sa_rf input will convert the input rainfall hyetograph to flow, deposited initially on the lowest 2D cell, then for subsequent timesteps distributed over all wet cells, within the 2d_sa_rf regions (ie. on the ground surrounding the building). Refer to TUFLOW <u>[[TUFLOW_Example_Models#Boundary_Condition_Options | Example model EG03_005.tcf]]</u> for a demonstration of this inflow boundary condition option.
# Exclude buildings from rainfall polygon and include it with 2d_sa_rf:
+
:::* Exclude buildings from the rainfall polygon and represent the rainfall that would fall on the building using <font color="blue"><tt>Read GIS SA RF PITS</tt></font> inflow boundaries. This approach is similar to the previous method, although instead of directing the inflow to the ground surrounding the building, it is directed into the sub-surface drainage (underground pipe) system. To implement this approach every 2d_sa_rf polygon must encompass at least one 1D pit. If multiple 1D pits are within a single 2d_sa_rf region, the flow from the polygon is split equally between the 1D pits. The pits can be snapped to a 1D node, connected via the pit-search distance or via an x-connector, to allow captured by the pit to enter a 1D representation of the sub-surface drainage system. Refer to TUFLOW <u>[[TUFLOW_Example_Models#Multiple_Domain_Model_Design | Example model EG15_000.tcf]]</u> for a demonstration of this inflow boundary condition option.
#* The cut out rainfall from the buildings can be included in the model again with <font color="blue"><tt>Read GIS SA RF</tt></font> command. Digitise a small 2d_sa_rf polygon for each building on the ground where runoff from the building is expected and assign attributes that convert rainfall hyetograph into flow and input it into the model as if it was 2d_sa. Check TUFLOW <u>[[TUFLOW_Example_Models#Boundary_Condition_Options | example model EG03_005.tcf]]</u> on how to setup 2d_sa_rf.
+
 
#* Most rainfall that falls on buildings collects on the roof and through gutters and downpipes finds its way directly into the sub-surface drainage system. In such case, <font color="blue"><tt>Read GIS SA RF PITS</tt></font> command can be used to feed the water directly into the pits. Every 2d_sa_rf polygon needs to have at least one 1D pit within to automatically deposit the flow. If multiple pits are present, the flow from the polygon splits equally.
+
== What is the recommended cell wet/dry depth for direct rainfall models? ==
 +
For models with SGS turned off, a reduced cell wet/dry depth of 0.0002m (0.0007ft) is recommended due to the substantial amount of shallow sheet flow. <br>
 +
If SGS is turned on, a reduced cell wet/dry depth is not necessary because the cell wet/dry calculation conducted on the SGS storage curve captures a greater change in depth. In the case, the default of 0.002m (0.007ft) may be sufficient; however, it should be reviewed and selected according to the magnitude of flooding depths.
 +
 
 
<br>
 
<br>
  

Latest revision as of 13:16, 23 August 2024

Direct Rainfall

Please see Direct Rainfall (Rain on Grid) Webinar.

TUFLOW Executable

  • The single precision version of TUFLOW is recommended for direct rainfall models using HPC.
  • The double precision version of TUFLOW is recommended for direct rainfall models using TUFLOW Classic.

Rainfall Boundary Conditions

The following links provide boundary condition guidance for applying direct rainfall:

Rainfall Losses and Soil Infiltration

TUFLOW supports a range of rainfall loss and soil infiltration options. These are split into two broad categories, defined by their respective calculation methods:

  • Rainfall Excess Losses
  • Soil infiltration Losses

The fundamental difference in calculation approach is described in the below Common Questions Answered (FAQ) section.

Rainfall Excess losses can be implement in a variety of ways, including

  1. Loss values specified within the TCF Read Materials File command, linking to TGC Read GIS Mat inputs;
  2. Loss values specified within the GIS layers associated with the TBC Read GIS SA RF command; or
  3. Global Rainfall Initial Loss and Global Rainfall Continuing Loss TBC commands.

Similarly, a variety of soil infiltration options are supported. The available options include, Initial / Continuing Infiltration, the Horton Infiltration method and the Green-Ampt Infiltration method. The following link provides further discussion on the Green-Ampt method.

Common Questions Answered (FAQ)

What is the difference between rainfall excess and soil infiltration?

Rainfall Excess Approach - Read Materials File and Read GIS SA RF continuing loss specification:

  • This initial and continuing loss approach is a simplistic calculation method comparable to the loss methods included in traditional hydrology models (e.g. RORB, URBS, WBNM etc).
  • The calculation approach is as follows:
    • The user defines the rainfall hyetograph (time vs depth (mm)) boundary condition inputs.
    • The rainfall value is reduced by the loss value (i.e. rainfall excess) before the boundary condition input is applied to the 2D cells.
    • The GIS SA RF takes the rainfall excess (units = m) and multiples the value by the Area attribute (units = m2) in the GIS object to convert rainfall depth to a volume (m3).
    • TUFLOW applies the calculated rainfall excess flow to the lowest cell in each GIS polygon during the first timestep that wetting occurs. Every timestep thereafter the inflow is distributed over the wet cells within the polygon.

Soil Infiltration Approach - TUFLOW Soils File (.tsoilf):

  • This approach is a more realistic representation of the actual physics associated with water infiltration into the soil.
  • The calculation approach is as follows:
    • The user defines the rainfall hyetograph (time vs depth (mm)) boundary condition inputs.
    • The total rainfall value is applied directly to every 2D cell within the 2d_rf polygon.
    • When a 2D cell is wet the soil infiltration function subtracts the appropriate loss volume of water from it. Computationally this is referred to as a “Sink” term.

What is the best approach for modelling buildings in rain on grid model?

There are numerous different industry standard ways to represent buildings in a direct rainfall model. Australian Rainfall and Runoff Guideline, Project 15 (Representation of Buildings in 2D Numerical Flood Models discusses some of the available options. Common TUFLOW modelling approaches are summarised below:

  1. Using depth-varying Manning's n over the area of the building footprint, with the application of a lower value (n = 0.02) at shallow depths (d < 0.03m) and a higher value (n > 0.3) at a more significant depths (d > 0.1m):
    • This is a very common and easy to implement option. The low Manning's n value aims to mimic the quick runoff response associated with drainage from the roof. The higher Manning's n value aims to represent the losses associated with deeper floodwater impacting the side of the building.
  2. Raise the building footprint elevation using TUFLOW's topographic update features (eg. Read GIS Z Shape):
    • Raising the model topography creates an obstruction to flow. It prevents floodwater from passing through buildings (as is the case with the Manning's n approach)
    • The application of rainfall on top of the building can however produce some undesirable results that require further consideration. Those being, water falling from the building roof to the ground can require a reduced model timestep to maintain model stability which can slow the simulation speed. Also, depending on the Map Cutoff Depth assumptions, water may be present in the results on the building rooftops. This may not be desired for mapping purposes. There are a variety of options available to resolve these issues:
  • Retain the model design described above, with subsequent post-processing of the mapped results before reporting. Either, delete the model result where there is overlap with the building footprint, or overlay the building footprint polygon objects over/above the result dataset in the GIS Map Layout (hiding the flood model result within the building footprints).
  • Exclude buildings from the rainfall polygon: This removes the rainfall from the model that would otherwise fall on the buildings. This approach will under-estimate the amount of rainfall entering the model. If the collective building footprint area is negligible in comparison to the entire model, this approach may be acceptable.
  • Exclude buildings from the rainfall polygon and represent the rainfall that would be falling on the building using Read GIS SA RF inflow boundaries. To do this, digitise a 2d_sa_rf polygon for each building (with a buffer of one of more 2D cells) where the building footprint has been excluded from the direct rainfall region. The 2d_sa_rf input will convert the input rainfall hyetograph to flow, deposited initially on the lowest 2D cell, then for subsequent timesteps distributed over all wet cells, within the 2d_sa_rf regions (ie. on the ground surrounding the building). Refer to TUFLOW Example model EG03_005.tcf for a demonstration of this inflow boundary condition option.
  • Exclude buildings from the rainfall polygon and represent the rainfall that would fall on the building using Read GIS SA RF PITS inflow boundaries. This approach is similar to the previous method, although instead of directing the inflow to the ground surrounding the building, it is directed into the sub-surface drainage (underground pipe) system. To implement this approach every 2d_sa_rf polygon must encompass at least one 1D pit. If multiple 1D pits are within a single 2d_sa_rf region, the flow from the polygon is split equally between the 1D pits. The pits can be snapped to a 1D node, connected via the pit-search distance or via an x-connector, to allow captured by the pit to enter a 1D representation of the sub-surface drainage system. Refer to TUFLOW Example model EG15_000.tcf for a demonstration of this inflow boundary condition option.

What is the recommended cell wet/dry depth for direct rainfall models?

For models with SGS turned off, a reduced cell wet/dry depth of 0.0002m (0.0007ft) is recommended due to the substantial amount of shallow sheet flow.
If SGS is turned on, a reduced cell wet/dry depth is not necessary because the cell wet/dry calculation conducted on the SGS storage curve captures a greater change in depth. In the case, the default of 0.002m (0.007ft) may be sufficient; however, it should be reviewed and selected according to the magnitude of flooding depths.


Up
Go-up.png Back to TUFLOW Modelling Guidance