FMA Challenge 3: Difference between revisions

Content deleted Content added
 
(47 intermediate revisions by 4 users not shown)
Line 2:
This demo comprises the third and final model of the FMA Challenge series and assumes that you're an intermediate to advanced TUFLOW user. <br>
 
This challenge investigates a typical non-urban stream from the California Central Valley. The stream originates near the Sierra Nevada Foothills and conveys runoff to the West into the Central Valley Flood Management System. The streams are natural, with large amounts of conveyance in their upper reaches. As the streams progress downstream, their conveyance capacity gradually reduces. At some point nearing the valley floor, typical levees have been constructed in an effor to containmaintain normalflooding floodsto within the the main flow conveyance systemsarea.
 
From this challenge, it is expected you will develop skills in:<br>
*Understanding of the influence of levees on flood behaviour;
*Understanding of the influence of channel roughness on flood behaviour;
*Understanding of the influence of infiltration on flood behaviour;
*Nested 1D/2D models; and
*Understanding the Green and Ampt infiltration method and USDR soil types.<br>
 
DataDData for this model is provided viain ZIPa compressedvariety fileof availabledifferent [http://www.tuflow.com/Tuflow%20Tutorial%20Models.aspxGIS compatible onformats. Download the TUFLOWdataset Downloadsthat Page]matches underthe 'DemoGIS Models'.<br>software you are using:
*[https://www.tuflow.com/Download/TUFLOW/Demo_Models/FMA_Challenge_Model_3_QGIS.zip QGIS Data Download]
*[https://www.tuflow.com/Download/TUFLOW/Demo_Models/FMA_Challenge_Model_3_Mapinfo.zip MapInfo Data Download]
*[https://www.tuflow.com/Download/TUFLOW/Demo_Models/FMA_Challenge_Model_3_ArcGIS.zip ArcGIS Data Download]
 
=Relevant Tutorials=
Other tutorials that are relevant to this challenge that may help refresh your memory are as follows:<br>
*1D-2D Linking - <u>[[Tutorial Module02Tutorial_M03|Tutorial Module 23]] </u>
*Running1D EventsNested and ScenariosChannel - <u>[[Tutorial Module10Tutorial_M11|Tutorial Module 1011]]</u>
*Running Scenarios - <u>[[Tutorial_M08|Tutorial Module 8]]</u>
*Running Events - <u>[[Tutorial_M09|Tutorial Module 09]]</u>
<br>
 
Line 26 ⟶ 32:
Cell size selection is managed through the use of variables.
 
Within the river bed, 1D sections spacing wereare generally between 500ft and 1000ft,. The exception exceptis at structures where theymore frequent sections are morenecessary frequentto adequately model changes in flow width and velocity in these areas.<br>
 
 
{| align="center" class="wikitable" width="50%"
Line 34 ⟶ 39:
! style="background-color:#005581; font-weight:bold; color:white;"| Active 2D Cells
|-
| 100ft || 338, 020
|-
| 200ft || 84, 818
|}
 
Linking between 1D/2D was along the top of the levees.
 
All outflow from the model was assumed to be only via the main creek (ie. there was no water level boundary applied to the 2D overbank domain).
 
The number of cross sections, grids, or mesh elements used were as follows:<br>
 
==Topography and Bed Resistance==
==Computational Domain Assembly and Execution==
The provided DEM via GMG format was too coarse for hydraulic modeling as key hydraulic features (such as levee crests) were not adequately represented. The first step was therefore to process the 20Gb of LiDAR data into a higher resolution DEM. A 10ft DEM was created, a small part of which is illustrated below.
 
Line 75:
|}
 
To check these Manning's 'n' assumptions, a sensitivity testing was completed on the in-bank roughness and is detailed further in the following sections.
 
==Boundary Conditions==
The 1d/2d linking was along the top of the levee. The elevations along these levees were extracted from the 10ft DEM, these and other significant features were included in the TUFLOW model as 3D GIS breakline layers, ensuring the hydraulic control is represented in the grid regardless of cell size.
 
===Open Boundaries===
The flood extent from the 100ft 2D grid model is shown below.
A single hydrograph, derived from a previous hydrologic analysis was applied as a QT boundary via a 1d_bc.
All outflow from the model was assumed to be only via the main creekchannel (ie. there was no water level boundary applied to the 2D overbank domain).
 
===1D/2D Linking===
The 1d/2d1D linkingriver channel model was linked to the 2D model using HX cells along the top of either natural or build levees on either side of the leveechannel. These locations defined the location of spill into the floodplain and are the key control on how much flow can interact between the floodplain and river conveyance area. The elevations along these levees were extracted from the 10ft DEM, these and other significant features were included in the TUFLOW model as 3D GIS breakline layers, ensuring the hydraulic control is represented in the grid regardless of cell size.
 
=Results - No Infiltration=
The flood extent from the 100ft 2D grid model with no infiltration and Manning's 'n' within the channel of 0.2 is shown below.
 
[[File:FMA3_3.jpg|600px]]
 
==Levee Overtopping Assessment==
===Sensitivity Creek Manning's n Test (Scenario 100ft n0.1)
To assess the timing and location of where levee banks were overtopped, the TGC was setup with the evacuation route command:<br> <font color="blue"><tt>Read GIS Z Shape Route </tt></font> <font color="red"><tt>==</tt></font> shp\2d_zshr_T3_levees_001_L.shp | shp\2d_zsh_T3_levees_001_P.shp.<br>
As discussed in the Manning’s n table above, the main creek n value of 0.20 is considered very high, especially in the lower reaches of the study area. A sensitivity analysis was carried out by lowering all the Manning’s n values in the main channel (modeled as 1D cross-sections) to 0.10.
Via review of the RC Map Output Data Type and the _RCP output point layer we can identify sections of the levee where a breach has occured (refer below).<br>
[[File:FMA3_6.jpg|600px]]
 
===Sensitivity Creek Manning's n Test (Scenario 100ft n0.1)==
As discussed in the Manning’s n table above, the main creek n value of 0.20 is considered very high, especially in the lower reaches of the study area. A sensitivity analysis was carried out by lowering all the Manning’s n values in the main channel (modeledmodelled as 1D cross-sections) to 0.10.
The image below shows the flood depths and extent, and for the outflowtwo isnon-infiltration includedsimulations inat the100ft ftpresolution download.(left with Ofchannel particular'n' interest= is0.2, thatright reducingwith thechannel 'n' value to= 0.1). has aSome significant effectdeviations on the arrival time of thein flood watersextent atcan thebe model outlet (much earlier)observed, andparticularly reduces the volumesouth of waterthe flowingmain ontochannel. the floodplainThis bycan aroundbe 20%largely dueattributed to the higher conveyance of the creek. Also of interest is that for the n=0.2 scenario, some overbank floodwaters return to the main creek near the model outlet causing a delayed second rise in the outlet flow hydrographs as illustrated in the chart further below. This effect does not occurchannel for the 'n'=0.1 scenario,case. with all overbank floodwaters remaining on the floodplain.<br>
[[File:FMA3_3.jpg|600px]][[File:FMA3_4.jpg|600px]]<br>
 
The effect of channel roughness can also be seen in hydrograph outflows from the model downstream boundary (refer figure below).
[[File:FMA3_4.jpg|600px]]
Of particular interest is that reducing the n value to 0.1 has a significant effect on the arrival time of the flood waters at the model outlet (much earlier), and reduces the volume of water flowing onto the floodplain by around 20% due to the higher conveyance of the creek. Also of interest is that for the n=0.2 scenario, some overbank floodwaters return to the main creek near the model outlet causing a delayed second rise in the outlet flow hydrographs as illustrated in the chart further below. This effect does not occur for the n=0.1 scenario, with all overbank floodwaters remaining on the floodplain.<br>
 
Please note the scenarios provided in the legend below are as follows:<br>
*100ft = 100ft cell resolution, no infiltration and channel roughness of 'n' = 0.2<br>
*200ft = 200ft cell resolution, no infiltration and channel roughness of 'n' = 0.2<br>
*100ft (n0.1) = 100ft cell resolution, no infiltration and channel roughness of 'n' = 0.1<br>
*100ft GA = 100ft cell resolution, Green and Ampt infiltration with channel roughness of 'n' = 0.2<br>
*Inflow = Hydrograph input at northeastern QT boundary.<br>
 
[[File:FMA3_5.jpg|600px]]
 
=Green and Ampt Infiltration Setup=
'''Chart of flow out of the model for the different scenarios simulated. Note that “100ft GA” is the infiltration scenario modeled for Challenge 3+.'''
 
TheProvided soils layers provided were roughly classified into two soil types for the purposes of demonstrating the Green-Ampt infiltration feature in TUFLOW. The soils are shown in the image below where:
 
 
 
The sections of levee that were overtopped were post-processed using GIS for the 100ft with n=0.2 scenario as illustrated by the magenta lines in the image below. The GIS layer showing the overtopped sections is provided as part of the ftp download.
 
[[File:FMA3_6.jpg|600px]]
 
=Solution to Challenge 3+=
Where unstated, the same parameters were used as per Challenge 3.
 
==Computational Domain Assembly and Execution==
The 100ft and 200ft models were used from Challenge 3 with the Creek Manning’s n value set to 0.2.
 
The soils layers provided were roughly classified into two soil types for the purposes of demonstrating the Green-Ampt infiltration feature in TUFLOW. The soils are shown in the image below where:
<ol>
<li>Red indicates a Sandy Loam.
Line 131 ⟶ 140:
[[File:FMA3_7.jpg|600px]]
 
The 100ft and 200ft simulations were re-run with the Green-Ampt infiltration switched on. The resulting peak flood depths map for the 100ft simulation is shown in the image below (right) whilst the 100ft without infiltration is provided on the right. As can be seen, the extent of inundation is significantly reduced with floodwaters infiltrating before reaching the model extents or returning to the main channel.
 
The mass balance reporting from TUFLOW indicates over half (58%) of the water infiltrates into the ground during the 171 hour simulation. For the flow out of the model see “100ft GA” in the chart of flow out of the model presented in Challenge 3.
 
The mass balance reporting from TUFLOW indicates over half (58%) of the water infiltrates into the ground during the 171 hour simulation. ForThis thereduction flowin outvolume ofcan thealso modelbe seeobserved “100ft GA” invia the chart“100ft ofGA” flow out of the model presented inhydrograph Challengeshow 3above.
[[File:FMA3_8.jpg|600px]]
 
[[File:FMA3_3.jpg|600px]][[File:FMA3_8.jpg|600px]]
==Challenges==
The soils GIS layers were difficult to work with. There is a lot of detail, but no clear way to correlate the different soils with Green-Ampt parameters (porosity, hydraulic conductivity, suction, etc).
 
=Conclusion=
Line 144 ⟶ 150:
In this challenge, we explored typical non-urban stream of the California Central Valley, with scenarios of various infiltration and flood levees adopted. From this, we gained a better understanding of the influence of flood levees on surface water behaviour, understanding of the Green Ampt infiltration method and USDR soil types, and a better understanding of nested 1D/2D models.
 
Congratulations on finishing Challengethe FMA Demo Model 3!. Please click here to return the [[Main_Page| Wiki Main Page]].