Difference between revisions of "Flood Modeller Tutorial Module01"

From Tuflow
Jump to navigation Jump to search
Line 127: Line 127:
 
<li> The internal 1D/2D boundaries section needs to be added. We need to update the culvert boundaries. We also need to include the HX boundaries for the open channel sections. </li>
 
<li> The internal 1D/2D boundaries section needs to be added. We need to update the culvert boundaries. We also need to include the HX boundaries for the open channel sections. </li>
  
The boundaries 1D/2D boundaries we need to add are
+
* The boundaries 1D/2D boundaries we need to add are
 +
'''Mapinfo'''
 +
:<font color="blue"><tt>Read GIS BC</tt></font> <font color="red"><tt>== </tt></font> mi\2d_bc_M04_culverts_001.MIF <font color="green"><tt>! This command reads in SX boundaries linking the 1D ESTRY culverts to the 2D domain</tt></font>
 +
:<font color="blue"><tt>Read GIS BC </tt></font> <font color="red"><tt>== </tt></font> mi\2d_bc_FMT_M01_HX_001.MIF <font color="green"><tt>! This command reads in HX boundaries linking the 1D Flood Modeller watercourse to the 2D domain</tt></font>
 +
 
 +
'''Other GIS'''
 
:<font color="blue"><tt>Read GIS BC</tt></font> <font color="red"><tt>== </tt></font> gis\2d_bc_M04_culverts_001_L.shp <font color="green"><tt>! This command reads in SX boundaries linking the 1D ESTRY culverts to the 2D domain</tt></font>
 
:<font color="blue"><tt>Read GIS BC</tt></font> <font color="red"><tt>== </tt></font> gis\2d_bc_M04_culverts_001_L.shp <font color="green"><tt>! This command reads in SX boundaries linking the 1D ESTRY culverts to the 2D domain</tt></font>
 
:<font color="blue"><tt>Read GIS BC </tt></font> <font color="red"><tt>== </tt></font> gis\2d_bc_FMT_M01_HX_001_p.shp | 2d_bc_FMT_M01_HX_001_P_L.shp <font color="green"><tt>! This command reads in HX boundaries linking the 1D Flood Modeller watercourse to the 2D domain</tt></font>
 
:<font color="blue"><tt>Read GIS BC </tt></font> <font color="red"><tt>== </tt></font> gis\2d_bc_FMT_M01_HX_001_p.shp | 2d_bc_FMT_M01_HX_001_P_L.shp <font color="green"><tt>! This command reads in HX boundaries linking the 1D Flood Modeller watercourse to the 2D domain</tt></font>

Revision as of 01:47, 27 October 2016

Introduction

In this module we will link an existing 2D TUFLOW domain to an existing Flood Modeller 1D model. The 2D TUFLOW domain will represent the floodplain in the study area, and the 1D Flood Modeller model will represent the watercourse and all online structures.

1D-2D linked models are able to utilise the individual benefits of 1D and 2D solution schemes. In this example, the 1D Flood Modeller scheme is used to represent the watercourses where the flow is essentially one-directional. A 2D scheme is suited to the representation of floodplains where a more detailed flow patterns may occur.

TUFLOW may be dynamically linked to 1D networks using the hydrodynamic solutions of ESTRY (TUFLOW 1D), Flood Modeller (previously ISIS), XP-SWMM and 12D Solutions’ Dynamic Drainage.

Existing Model Data

Introduction

This section of the tutorial provides a summary of the background Flood Modeller and TUFLOW model data that has been provided for the purposes of this tutorial. This tutorial builds upon the 2D TUFLOW domain that was constructed as part of Module 1 and Module 2 of the TUFLOW Tutorial Model. This TUFLOW model consists of:

  • Definition of Active/In-Active Areas
  • Definition of Land Use
  • 1D ESTRY culverts
  • 1D/2D boundary links to connect the 1D ESTRY culverts to the 2D TUFLOW domain.

The provided model differs from Module 2 in that the 2D boundary conditions (upstream inflows and downstream stage-discharge boundary) have been removed from the model. These will instead be represented in Flood Modeller as it is a more typical schematisation for a 1D/2D linked model. Hydraulic structures along the watercourse have also been removed and have been represented in the Flood Modeller 1D model. Amended versions of the TUFLOW control files and GIS files have been provided as part of this tutorial in preparation for linking to Flood Modeller.

Flood Modeller

A complete Flood Modeller 1D model of the watercourse has been provided for the purposes of this tutorial. The model data has been arranged in the following folder structure. For further information on each file type, please refer to the Flood Modeller User Manual:

  • DAT: Contains the Flood Modeller 1D network file.
  • IED: Contains event data of the model.
  • IEF: A blank folder in which to store Flood Modeller Event Files used to simulate the model.
  • RES: A blank folder in which to write the Flood Modeller result files.

For this tutorial, the floodplain of the study area will be modelled entirely in TUFLOW. The cross-sections in the Flood Modeller 1D model have been trimmed to the top of bank to ensure there is no double-counting of storage within the floodplain. Upstream and downstream boundary data is contained separately within the IED files. These will be referenced and linked to the Flood Modeller 1D model via the IEF 1D simulation file.

GIS and Model Inputs

The steps necessary to modify each of the GIS inputs are demonstrated in MapInfo, ArcGIS and QuantumGIS (QGIS). At each stage please select your GIS package to view relevant instructions.

Define the External 1D Networks

This part of the module creates the GIS layers that specify the geographic location of the Flood Modeller nodes that are to be connected to the 2D domain. Follow the instructions below for your preferred GIS package.

Define the 1D/2D Boundary Links

This part of the module creates the 1D/2D boundaries to link the Flood Modeller 1D component to the TUFLOW 2D domain. Follow the instructions below for your preferred GIS package.

Define Bank Elevations

This part of the module defines the bank elevations of the watercourse which are the elevations of the 1D/2D boundary links created in the previous section. Follow the instructions below for your preferred GIS package.

Deactivate 2D cells

This part of the module describes the steps to deactivate the 2D cells where the 1D model is replacing the 2D solution. Follow the instructions below for your preferred GIS package.

Modify Simulation Control Files

Now that we have made all of the necessary changes to the GIS layers, we need to update the TUFLOW control files and Flood Modeller Simulation Files to create a linked model.

TUFLOW Control Files

TCF File

  1. Firstly save a copy of the TUFLOW file created as a part of module 2 (M02_5m_001.tcf) as FMT_M01_001.tcf.
  2. The simulation time control commands (shown below) which specify timestep, start time and end time need to be removed from the tcf. and moved to an IEF file.
    ! SIMULATION TIME CONTROL COMMANDS
    Timestep == 1.5 ! Specifies a 2D computational timestep of 1.5 seconds
    Start time == 0 ! Specifies a simulation start time of 0 hours
    End Time == 3 ! Specifies a simulation end time of 3 hours
  3. Delete the 1D Domain commands (shown below) from the tcf. These commands should be moved to the ecf. This an optional step.
    !1D DOMAIN COMMANDS
    Start 1D Domain ! Indicates the start of a block of 1D or ecf (Estry Control File) commands
    Output Folder == ..\results\M02\1d ! Specifies the location of the 1D output files
    Write Check Files == ..\check\1d\ ! Specifies the location of the 1D check files and prefixes them with the .tcf filename
    Output Interval (s) == 180 ! Writes the 1D output every 3 minutes. Note the (s) is required to indicate the value is in seconds
    Timestep == 0.75 ! Specifies a 1D computational timestep of 0.75 seconds. We will set the 1D timestep to be half the 2D timestep
    Read GIS Network == ..\model\gis\1d_nwk_M02_culverts_001_L.shp ! Reads the 1d_nwk GIS layer defining culverts within the model
    Read GIS Network == ..\model\gis\1d_nwk_M02_culverts_001_P.shp ! Reads the 1d_nwk GIS layer defining nodes used to create automatic SX 1D/2D boundaries
    End 1D Domain
  4. Flood modeller commands need to be added, these are:
  5. Mapinfo
    Read GIS X1D Nodes == ..\model\mi\1d_x1d_FMT_M01_nodes_001.MIF ! GIS layer referencing node IDs from Flood Modeller
    Read GIS X1D Network == ..\model\mi\1d_x1d_FMT_M01_nwk_001.MIF ! GIS layer representing channels to allow for the digitisation of Water Level Lines (optional)
    Read GIS X1D WLL == ..\model\mi\1d_x1d_WLL_FMT_M01_001.MIF ! GIS layer containing WLLs for visualising 1D results in 2D
    Other GIS
    Read GIS X1D Nodes == ..\model\gis\1d_x1d_FMT_M01_nodes_001_p.shp ! GIS layer referencing node IDs from Flood Modeller
    Read GIS X1D Network == ..\model\gis\1d_x1d_FMT_M01_nwk_001_L.shp ! GIS layer representing channels to allow for the digitisation of Water Level Lines (optional)
    Read GIS X1D WLL == ..\model\gis\1d_x1d_WLL_FMT_M01_001_p.shp ! GIS layer containing WLLs for visualising 1D results in 2D
  6. This concludes the changes that we need to make to the .tcf, but before we run the model we need to update the .ecf, .tgc, .tbc and bc_dbase files.
  7. TGC File

    There have been changes made to the geometry:

    • We have added a third model grid command.
    • We have added topography amendments.
  8. Open FMT_M01_001.tgc or navigate to the correct tab. Add an extra line in  ! Model Grid Commands reading:
  9. Mapinfo
    Read GIS Code BC == mi\2d_bc_FMT_M01_001_HX_001.MIF ! Deactivates the cells where the creek channel has been modelled in 1D
    Other GIS
    Read GIS Code BC == gis\2d_bc_FMT_M01_001_HX_001_p.shp | 2d_bc_FMT_M01_HX_001_L.shp  ! Deactivates the cells where the creek channel has been modelled in 1D
  10. Topography amendments should be added in a new section. These are:
  11. Mapinfo
     ! Read GIS Z Shape == mi\2d_zsh_M03_thalweg_001.MIF  ! This GIS layer is commented out / removed from the model as the creek channel has been modelled in 1D
    Read GIS Z Shape == mi\2d_zsh_M03_Rd_Crest_001.MIF ! GIS layer used to define the road crest
    Read GIS Z HX Line MAX == mi\2d_bc_FMT_M01_HX_001.MIF! Defines the bank crest levels (1D/2D boundary cell elevations). The 'MAX' option prevents any zpt elevations from being lowered
    Other GIS
     ! Read GIS Z Shape == gis\2d_zsh_M03_thalweg_001.shp  ! This GIS layer is commented out / removed from the model as the creek channel has been modelled in 1D
    Read GIS Z Shape == gis\2d_zsh_M03_Rd_Crest_001_L.shp | gis\2d_zsh_M03_Rd_Crest_001_P.shp  ! GIS layer used to define the road crest
    Read GIS Z HX Line MAX == gis\2d_bc_FMT_M01_HX_001_L.shp | gis\2d_bc_FMT_M01_HX_001_p.shp ! Defines the bank crest levels (1D/2D boundary cell elevations). The 'MAX' option prevents any zpt elevations from being lowered
  12. Save the file. The geometry file is now ready to be used.

    TBC File

  13. Open FMT_M01_001.tbc or navigate to the tab.
  14. The internal 1D/2D boundaries section needs to be added. We need to update the culvert boundaries. We also need to include the HX boundaries for the open channel sections.
    • The boundaries 1D/2D boundaries we need to add are
    Mapinfo
    Read GIS BC == mi\2d_bc_M04_culverts_001.MIF ! This command reads in SX boundaries linking the 1D ESTRY culverts to the 2D domain
    Read GIS BC == mi\2d_bc_FMT_M01_HX_001.MIF ! This command reads in HX boundaries linking the 1D Flood Modeller watercourse to the 2D domain
    Other GIS
    Read GIS BC == gis\2d_bc_M04_culverts_001_L.shp ! This command reads in SX boundaries linking the 1D ESTRY culverts to the 2D domain
    Read GIS BC == gis\2d_bc_FMT_M01_HX_001_p.shp | 2d_bc_FMT_M01_HX_001_P_L.shp ! This command reads in HX boundaries linking the 1D Flood Modeller watercourse to the 2D domain
  15. Save the file. The boundary file is now ready to be used

    Flood Modeller Simulation Files

    IEF File

    Run the Simulation

    There are a number of different methods in which to start a Flood Modeller – TUFLOW model. Please refer to the following page for a detailed discussion on each of these methods. Using your preferred method, start the model FMT_0100F_001.ief. If the simulation fails to start, please refer to the troubleshooting guidance on this page.


    Click here to return back to the main Flood Modeller - TUFLOW Tutorial page.