Jump to content

Unidirectional flow through 2D cells?

Recommended Posts

I have a simple but large 2d only model with direct rainfall applied to create runoff.  I have introduced embankments to protect certain areas from flooding using some zlns.  But with the direct rainfall approach (which I cannot move away from for reasons beyond the scope of this query) I obviously get ponding behind the embankments in the areas I am trying to protect and sometimes this is deeper than the flooding from the river without the embankments.  Because of the level of study I want to introduce something to allow flow to pass through the embankments as would be allowed by the downstream flood level (to all intents and purposes a flapped culvert).  But given the nature of the model that would a very onerous process, so to get to the question... can unidirectional flow be set up within the 2d domain without having to introduce any 1d components? 

Share this post

Link to post
Share on other sites

Hi Oppley,

Thank you for your query.  An unidirectional flow feature is not available in TUFLOW 2D modelling at the moment. However, you can use a 2d_sapo layer to pass flow through the embankments based on the downstream flood level. The details of the "Read GIS SA PO" option is documented in the Section of the TUFLOW manual. Basically, you need to:

  1. Create a 2d_po point object of Type “H_” at the water level location that you want to use to control the flow rate.
  2.  Add “Read GIS PO == ….” to the .tcf file if not already there. 
  3. Create a new 2d_sa layer (call it 2d_sapo) and add two new attributes, as listed in Table 7-6 of the TUFLOW manual.  
    (i)  PO_Type:   Char of length 16 
    (ii)  PO_Label:  Char of max length 40 
  4. Digitise two SA polygons, one at the area upstream the embankments (where "ponding" happens), and another one downstream the embankments. For the attributes: 
    (i)  Set the Name attribute to the name of the water level vs flow curve in the BC database. 
    (ii)  Set PO_Type to “H_”. 
    (iii) Set PO_Label to the PO Label of the relevant 2d_po “H_” point to be used to determine the flow from the h vs Q curve. 
  5. Create a .csv file to define the stage-discharge relationships. Note that for the SA polygons upstream the embankments, the flow rate must be set as negative values to extract water, while for the SA polygons downstream the embankments, use exactly same flow rate but in positive values. 
  6.  Add “Read GIS SA PO == …” to the .tbc file, and the the "Names" of the SA polygons to the BC database

Please view the _MB.csv files to check the SA in/outflow from the polygons. The extracted flow may be smaller than the inflow to the downstream, if the water upstream the embankments is "dried" by the SA polygon. In that case, you can simply use a QT inflow for the downstream SA polygon based on the actually extracted flow rate to ensure the mass balance is conserved. But this process will be iterative.

Hope this helps and please feel free to get back to us or send email queries to support@tuflow.com if you have any troubles setting up the model.




Share this post

Link to post
Share on other sites

Join the conversation

You can post now and register later. If you have an account, sign in now to post with your account.

Reply to this topic...

×   Pasted as rich text.   Paste as plain text instead

  Only 75 emoji are allowed.

×   Your link has been automatically embedded.   Display as a link instead

×   Your previous content has been restored.   Clear editor

×   You cannot paste images directly. Upload or insert images from URL.

  • Create New...