Flood Modeller Tutorial Module01

From Tuflow
Revision as of 21:35, 11 November 2016 by Rachel (talk | contribs) (→‎TGC File)
Jump to navigation Jump to search

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.

As discussed in Module 1, the main creek channel is not very well represented using the 5m 2D cell size. In parts, the creek is only 5-10m wide and the 5m cell size could be considered too coarse to accurately represent the creek topography. Refer to the below figure.
Poor 2d rep.png
Using a cell size that is coarse relative to the width of the creek channel may reduce the accuracy of the conveyance in the channel and introduce a source of mass error to the model. There are two options for improving the representation of the creek channel:

  • decrease the width of the 2D cells; and/or
  • model the channel as a 1D network, dynamically linked to the 2D domain (the floodplain).

In the optional section of module 1, we looked at reducing the cell size to get a better representation of the channel. In this module we will adopt the second approach of modelling the creek as 1D elements in Flood Modeller.

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.
Setting up a 1D/2D model where the 1D channel cuts through the 2D domain is probably the most time-consuming type of a model to setup. However, the reduction in simulation time can be beneficial and make this a good approach. For this module, the complete Flood Modeller 1D model has been provided, so you can progress through the module in a relatively short period of time.

Existing Model Data

This tutorial builds upon the 2D TUFLOW domain that was constructed as part of Module 1 and Module 2 of the TUFLOW Tutorial Models.
The model developed in these tutorial modules already contains some culverts modelled as 1D elements. The culverts are modelled in ESTRY, TUFLOW's internal 1D engine. One of these culverts will be kept in ESTRY and the other will be changed to Flood Modeller.
The 2D boundary conditions (upstream inflows and downstream stage-discharge boundary) will be 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.
The existing 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.

If you need more information on these elements, we recommend that you work through Module 1 and Module 2.

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:

  • 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 further information on each file type, please refer to the Flood Modeller User Manual.
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 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 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. We need to remove the simulation time control commands (shown below) which specify timestep, start time and end time from the tcf. These will be moved to the Flood Modeller IEF file soon.
    ! 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. We need to add the Flood modeller commands, 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

    Now we will update the TUFLOW geometry Control (TGC) file.
    The following changes have been changes made to the geometry:

    • We have turned off the cells in the channel (as these are now modelled in Flood Modeller)
    • We have added the bank lines as 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 correct tab.
  14. We need to add the internal 1D/2D boundaries. We need to update the culvert boundaries and include the HX boundaries for the open channel sections.
    • The 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.