Advection Dispersion Modelling: Difference between revisions

Content deleted Content added
 
(9 intermediate revisions by 4 users not shown)
Line 1:
=Introduction=
TUFLOW’s Advection Dispersion (AD) Modulefunctionality is an extension of the TUFLOW Classic/HPC engines available within the TUFLOW CATCH module. It adds to the hydrodynamic capabilities of TUFLOW Classic/HPC by simulating depth-averaged, two and one-dimensional constituent fate and transport. An example of such a constituent might include salinity. Both dissolved and particulate constituents can be simulated. TUFLOW AD takes depth and velocity fields computed by the TUFLOW Classic and HPC solvers and uses this information, together with initial and boundary conditions, to simulate the advection and dispersion of user-defined constituents.
 
TUFLOW AD is specifically oriented towards such analyses in systems including coastal waters, estuaries, rivers, floodplains and urban areas. The AD modulefunctionality is discussed in detail in the <u>[https://downloadsdocs.tuflow.com/TUFLOWclassic-hpc/Releasesmanual/2020-10latest/TUFLOW%20AD%20Manual.2020-10-AF.pdf TUFLOW AdvectionManual Dispersion- ManualChapter 9]</u>.
 
TUFLOW FV also has Advection Dispersion functionality. In many cases, it's preferable to use TUFLOW FV for AD modelling due to its functionality with flexible mesh. See the <u>[https://docs.tuflow.com/fv/manual/latest/ TUFLOW FV Manual]</u> for details.
 
=Model Development=
==Setting Up a New Model==
The steps below describe the process for setting up a TUFLOW AD model. It is assumed that the user is familiar with TUFLOW Classic/HPC and that the folder structure for TUFLOW has been setup with all required files. The user should run the TUFLOW Classic/HPC model without the AD modulefunctionality first to make sure that it is appropriately configured and stable.
 
<ol>
Line 19 ⟶ 21:
 
<li>Set up the boundary condition tables (.csv file(s)) to define the time-varying constituent concentrations at any input boundaries.
<ul><li>Set up the constituent boundary condition table(s) in the “bc_dbase” folder . For example in the below, the time-varying concentration of constituents Conc_AD1 and Conc_AD2 are set:<br>
[[File:EG17 conc 001.png]]<br></ul>
<li>Set up up the boundary condition database (.csv file)
Line 25 ⟶ 27:
[[File:EG17 AD 001.png]]<br>
<li>In the .adcf file use the "AD BC Database" command to set the location of the bc database as follows.<br><font color="blue"><tt>AD BC Database</tt></font><font color="red"><tt> == </tt></font><font color="black"><tt>..\bc dbase\my_ad_bc_dbase.csv</tt></font></ul>
<li>Setup up TUFLOW to activate the AD modulefunctionality (.tcf file)
<ul><li>In the .tcf file use the command "AD Control File" to set the location of the adcf and activate execution of the AD modulefunctionality as follows.<br><font color="blue"><tt>AD Control File</tt></font><font color="red"><tt> == </tt></font><font color="black"><tt>ad_run.adcf</tt></font></ul>
<li>Run the model
<li>Run TUFLOW as normal. The AD modulefunctionality will be calledutilised and appropriate constituent result output files written.
 
</ol>
Line 37 ⟶ 39:
= Common Questions Answered (FAQ) =
 
== How can the Advection Dispersion modulefunctionality be used to determine the time of concentration in a 1D-2D TUFLOW model? ==
The Advection Dispersion modulefunctionality can track particles and determine the time of concentration by simulating how a particle of water travels from an upstream to a downstream location. However, the AD modulefunctionality is only available for 2D domains and cannot directly operate within a 1D channel.
 
To utilise the AD modulefunctionality in this case, the 1D channel would need to be converted into a 2D domain using a Quadtree grid. This conversion involves refining the 2D cells within the channel to smaller sizes, ensuring the model accurately represents the flow behaviour. Once this modification is complete, the AD modulefunctionality can provide detailed insights into particle travel times and help demonstrate the effects of reprofiling or attenuation measures on downstream areas.
 
== Can TUFLOW generate 2D Plot (Time-Series) Output for Advection Dispersion results? ==
Line 47 ⟶ 49:
However, results can be extracted using output zones. Defining smaller output zones allows high-frequency data to be generated for areas of interest while managing file sizes efficiently. Multiple output zones can also be used to monitor widely separated locations.
 
== How can initial tracer concentrations and SGS parameters be managed in the Advection Dispersion modulefunctionality? ==
Initial tracer concentrations can be applied to dry cells, and these concentrations are mobilised as the cells become wet during a simulation. The initial water level in dry cells is set as the bed elevation plus the Cell Wet/Dry Depth. This depth determines the initial tracer volume available for advection once the cell becomes inundated.
 
When SGS (Sub-Grid Sampling) is used, the initial water volume is derived from a pre-calculated “level vs cell volume” curve. Tracer concentrations are distributed across this calculated volume. This ensures accurate representation of tracer movement, even in partially wet cells.
 
== How can the Advection Dispersion modulefunctionality be used to determine water residency time? ==
The Advection Dispersion modulefunctionality in TUFLOW can be used to calculate water residency time by modelling it as a scalar variable. This approach provides a method for tracking the duration water has spent within a specific area, such as a wetland, and is visualised in the model output as a time-based scalar field. This approach has some limitations:
* Output Capabilities: While TUFLOW supports various output formats (e.g., XMDF, DAT, NC), extracting detailed time-series data for specific constituents at individual locations may require additional post-processing. The current AD modulefunctionality does not explicitly support direct Point Output (PO) functionality for constituent data.
* Post-Processing Requirements: To obtain detailed residency time information at specific points, output zones may be required along with refined post-processing techniques. Defining output zones allows high-frequency scalar data to be captured in areas of interest, which can then be analysed to estimate residency times.
* Engine-Specific Features: Unlike the TUFLOW FV engine, which includes a particle tracking module for explicit tracking of water age, the fixed grid engine’s AD module relies on scalar-based methods to approximate residency time. For calculating water residency time with the AD modulefunctionality, properly configuring scalar outputs and planning post-processing steps are essential for accurate results.
 
== How can the Advection Dispersion modulefunctionality simplify firewater containment modelling? ==
The Advection Dispersion modulefunctionality in TUFLOW can simplify firewater containment modelling by using passive tracers to track firewater flow and concentration. Instead of running separate simulations for rainfall and firewater scenarios, the AD modulefunctionality enables a single simulation where tracers represent the firewater. This approach reduces modelling complexity while maintaining accuracy.
 
For example, a model with direct rainfall over the entire domain applies a passive tracer via 2d source area (2d_sa) polygons. The output can be set up to identify areas with tracer concentrations above a certain threshold, distinguishing firewater extents from other inundated areas. Areas outside of this represent zones with zero tracer concentration. Tracers can also include decay and settling parameters for added flexibility. This method not only simplifies the process but also ensures compliance with the UK CIRIA (Construction Industry Research and Information Association) guidance by integrating rainfall and firewater scenarios into a single simulation.
 
== What guidance is available for Non-Newtonian mixing exponents and dispersion coefficients in the Advection Dispersion modulefunctionality? ==
The Non-Newtonian Mixing Exponents (m, o, and p) were introduced in the 2023-03-AC release to improve how TUFLOW models non-Newtonian fluids. These exponents control how yield stress and density change as fluid concentration varies.
 
Line 78 ⟶ 80:
Currently, there is no guidance for dispersion coefficients when mixing pure water with non-Newtonian fluids. Suitable values should be determined based on laboratory tests or studies specific to the fluid being modelled.
 
== What are the limitations of the Advection Dispersion modulefunctionality when modelling Non-Newtonian flow through 1D elements? ==
When using the Advection Dispersion modulefunctionality for non-Newtonian flow in models that include 1D elements, the following simplifications and limitations apply:
 
Flow Calculation:
Line 90 ⟶ 92:
These limitations mean that while the AD module can be used in models with 1D elements, it does not fully capture the complexities of non-Newtonian fluid behaviour in those elements. If higher accuracy is required, 2D elements where non-Newtonian properties are more comprehensively represented may be considered.
 
== Can the Advection Dispersion modulefunctionality be used in a fixed grid TUFLOW HPC or SGS model to simulate salinity? ==
Yes, the Advection Dispersion (AD) modulefunctionality can be used with fixed grid TUFLOW HPC and SGS models to simulate salinity transport, as long as the system is relatively well mixed vertically.
 
TUFLOW HPC is a 2D shallow water equation solver, so it does not account for vertical salinity gradients. This means it is best suited for rivers, shallow lakes, and estuaries where haloclines are weak or absent. Different boundary inflows can be assigned varying time series of salinity concentrations, and the module will simulate how these mix over space and time.
Line 97 ⟶ 99:
If the system has strong density stratification or distinctly three dimensional flow behaviour, a fully 3D solver such as TUFLOW FV is recommended.
 
== Can TUFLOWthe HPCAdvection Dispersion functionality simulate pollutant runoff and transport from catchments? ==
Yes, TUFLOW HPC’sthe Advection Dispersion modulefunctionality can also be used to simulate pollutant runoff from catchments. This includes modelling the generation, transport, and fate of pollutants such as suspended sediment, nitrogen, and phosphorus.
 
Pollutants can be generated based on modelled bed shear stress, then transported through the domain using the flow field. The model can also apply decay and settling rates as part of the simulation.
 
The TUFLOW CATCH module, with which the Advection-Dispersion functionality is included, provides the capability to simulate pollutant runoff and transport from catchments. See the <u>[https://www.tuflow.com/products/tuflow-catch-module/ TUFLOW Website]</u> for further information.
For more complex water quality modelling needs, contact TUFLOW Support (support@tuflow.com) for guidance on best practices.<br>
 
<br>{{Tips Navigation
|uplink=[[ TUFLOW_Modelling_Guidance | Back to TUFLOW Modelling Guidance]]
}}