Difference between revisions of "Sediment Transport with Particle Tracking"

From TUFLOW FV Wiki
Jump to: navigation, search
(Particle Tracking)
(Particle Tracking)
Line 22: Line 22:
 
Generate a new file called FMA2_SED_003a.fvptm.  This is where we will generate the characteristics of our particle tracking module.  We will use the same sediment characteristics that we have already applied to the sediment transport model.  Add the following commands
 
Generate a new file called FMA2_SED_003a.fvptm.  This is where we will generate the characteristics of our particle tracking module.  We will use the same sediment characteristics that we have already applied to the sediment transport model.  Add the following commands
  
<font color="green">! PTM to assess fate of STP plume - TEST POLYGON SEEDING
+
<font color="green">! PTM to assess fate of STP plume - TEST POLYGON SEEDING
! use nodestring definition from HD model to allow particles to leave model domain. </font>
+
! use nodestring definition from HD model to allow particles to leave model domain. </font>
Open boundary nodestring == 2  ! downstream tidal boundary
+
<font color="blue"><tt>Open boundary nodestring</tt></font><font color="red"><tt> ==</tt></font> <font color="green"! downstream tidal boundary</font>
  
!TIME COMMANDS
+
<font color="green">!!TIME COMMANDS
!_________________________________________________________________
+
!_________________________________________________________________</font>
lagrangian timestep == 120. ! seconds
+
<font color="blue"><tt>lagrangian timestep</tt></font><font color="red"><tt> ==</font> 120. <font color="green"! seconds</font>
eulerian timestep == 120. ! seconds
+
<font color="blue"><tt>eulerian timestep</tt></font><font color="red"><tt> ==</font> 120. <font color="green"! seconds</font>
 
+
<font color="green">! Sediment Transport COMMANDS (required if a particle group interacts with bed)</font>
! Sediment Transport COMMANDS (required if a particle group interacts with bed)
+
<font color="blue"><tt>bed roughness model</tt></font><font color="red"><tt> ==</font> ks
bed roughness model == ks
+
<font color="blue"><tt>bed roughness parameters</tt></font><font color="red"><tt> ==</font> 0.01,0.01 <font color="green"! ksc, ksw</font>
bed roughness parameters == 0.01,0.01 ! ksc, ksw
+
<font color="green">!Global COMMANDS</font>
 
+
<font color="blue"><tt>Nscalar</tt></font><font color="red"><tt> ==</font> 1
!Global COMMANDS
+
<font color="green">!PARTICLE GROUP COMMANDS </font>
Nscalar == 1
+
<font color="green">!_________________________________________________________________</font>
 
+
<font color="blue"><tt>Group</tt></font><font color="red"><tt> ==</font> fineSed
!PARTICLE GROUP COMMANDS  
+
<font color="blue"><tt>d50</tt></font><font color="red"><tt> ==</font> 0.000002
!_________________________________________________________________
+
      <font color="blue"><tt>particle density</tt></font><font color="red"><tt> ==</font> 2650
 
+
      <font color="blue"><tt>Settling model</tt></font><font color="red"><tt> ==</font> constant
Group == fineSed
+
      <font color="blue"><tt>settling parameters</tt></font><font color="red"><tt> ==</font> 0.0002 <font color="green"!(m/s)</font>
    d50 == 0.000002
+
      <font color="blue"><tt>deposition model</tt></font><font color="red"><tt> ==</font> ws0
    particle density == 2650
+
      <font color="blue"><tt>Erosion Model</tt></font><font color="red"><tt> ==</font> Mehta
    Settling model == constant
+
      <font color="blue"><tt>Erosion parameters</tt></font><font color="red"><tt> ==</font> 0.1, 0.5, 0.5 <font color="green"!Er, taucr, alpha</font>
    settling parameters == 0.0002 !(m/s)
+
<font color="blue"><tt>End Group</font>
    deposition model == ws0
+
<font color="blue"><tt>Group</tt></font><font color="red"><tt> ==</font> Gravel
    Erosion Model == Mehta
+
      <font color="blue"><tt>d50</tt></font><font color="red"><tt> ==</font> 0.032
    Erosion parameters == 0.1, 0.5, 0.5 !Er, taucr, alpha
+
      <font color="blue"><tt>particle density</tt></font><font color="red"><tt> ==</font> 2650
End Group
+
      <font color="blue"><tt>Settling model</tt></font><font color="red"><tt> ==</font> constant
 
+
      <font color="blue"><tt>settling parameters</tt></font><font color="red"><tt> ==</font> 0.6 <font color="green"!(m/s)</font>
Group ==  Gravel
+
      <font color="blue"><tt>Critical stress model</tt></font><font color="red"><tt> ==</font> Soulsby
  d50 == 0.032
+
      <font color="blue"><tt>Bed load model</tt></font><font color="red"><tt> ==</font> MPM_Shimizu
  particle density == 2650
+
      <font color="blue"><tt>Bed load parameters</tt></font><font color="red"><tt> ==</font> 8.0, -1 ,1.5
  Settling model == constant
+
<font color="blue"><tt>End Group</font>
  settling parameters == 0.6 !(m/s)
+
<font color="green"!_________________________________________________________________</font>
  Critical stress model == Soulsby
+
<font color="green"! This is required due to adding deposition and settling.</font>
  Bed load model == MPM_Shimizu
+
<font color="blue"><tt>Material</tt></font><font color="red"><tt> ==</font> 0
  Bed load parameters == 8.0, -1 ,1.5
+
      <font color="blue"><tt>Layer</tt></font><font color="red"><tt> ==</font> 1
End Group
+
            <font color="blue"><tt>dry density</tt></font><font color="red"><tt> == </font>1890.,1890
 
+
      <font color="blue"><tt>End layer</font>
!_________________________________________________________________
+
<font color="blue"><tt>end material</font>
 
+
<font color="green"!Upstream</font>
! This is required due to adding deposition and settling.
+
<font color="blue"><tt>seed particles</tt></font><font color="red"><tt> ==</font> point,  10796.514,8285.014
Material == 0
+
      <font color="blue"><tt>particle groups</tt></font><font color="red"><tt> ==</font> fineSed, Gravel
Layer == 1
+
      <font color="blue"><tt>group mass</tt></font><font color="red"><tt> ==</font> 100,100
dry density == 1890.,1890
+
<font color="blue"><tt>end seed</font>
End layer
+
<font color="green"! OUTPUT SETTINGS</font>
 
+
<font color="green"!_________________________________________________________________</font>
end material
+
<font color="blue"><tt>output dir</tt></font><font color="red"><tt> ==</font> ..\results\
 
+
<font color="blue"><tt>output</tt></font><font color="red"><tt> ==</font> ptm_netcdf
!Upstream
+
      <font color="blue"><tt>output groups</tt></font><font color="red"><tt> ==</font> all
seed particles == point,  10796.514,8285.014
+
      <font color="blue"><tt>output parameters</tt></font><font color="red"><tt> ==</font> age, state_age, mass, uvw, uvw_water, depth, water_depth
particle groups == fineSed, Gravel
+
      <font color="blue"><tt>output interval</tt></font><font color="red"><tt> ==</font> 300.
group mass == 100,100
+
<font color="blue"><tt>end output</font>
end seed
+
 
+
! OUTPUT SETTINGS
+
!_________________________________________________________________
+
 
+
output dir == ..\results\
+
 
+
output == ptm_netcdf
+
  output groups == all
+
  output parameters == age, state_age, mass, uvw, uvw_water, depth, water_depth
+
  output interval == 300.
+
end output
+

Revision as of 02:40, 27 April 2021

Particle Tracking

The TUFLOW Particle Tracking Module (PTM) allows the 2D or 3D simulation of discrete Lagrangian particles as they are transported by the flow field (or other drivers such as wind or waves). Particle behavior such as settling, buoyancy, decay, sedimentation and re-suspension can all be simulated. It can be run in conjunction with the sediment transport module to simulated the fate and age of sediment particles within the hydraulic model.

In this optional exercise we will take the existing FMA2_SED_003 model and add some particle tracking to track the input of particles.

Copy and paste the existing FMA2_SED_003.fvsed file and rename as FMA2_SED_003a.fvsed.

Under the Boundary Conditions Block add the following particle tracking block to reference a particle tracking file that we will generate.

!PARTICLE TRACKING INPUT
Particle Tracking Control File == FMA2_SED_003a.fvptm	

In the output commands block, under the current specification for the XMDF file add the following block to output results to a netcdf file. This will output hydraulic parameters, sediment transport parameters as well as results from the particle tracking module.

output == netcdf
      output parameters == h,v,d, Rhow, Taub, TauC, PTM_1, PTM_2, PTM_BED_2
      output interval == 900.
end output

Save and close the file.

Generate a new file called FMA2_SED_003a.fvptm. This is where we will generate the characteristics of our particle tracking module. We will use the same sediment characteristics that we have already applied to the sediment transport model. Add the following commands

! PTM to assess fate of STP plume - TEST POLYGON SEEDING
! use nodestring definition from HD model to allow particles to leave model domain. 
Open boundary nodestring == 2  <font color="green"! downstream tidal boundary</font>
!!TIME COMMANDS
!_________________________________________________________________
lagrangian timestep ==</font> 120. <font color="green"! seconds</font>
<tt>eulerian timestep</tt><tt> == 120. <font color="green"! seconds</font>
! Sediment Transport COMMANDS (required if a particle group interacts with bed)
<tt>bed roughness model</tt><tt> == ks
<tt>bed roughness parameters</tt><tt> == 0.01,0.01	<font color="green"! ksc, ksw</font>
!Global COMMANDS
<tt>Nscalar</tt><tt> == 1
!PARTICLE GROUP COMMANDS 
!_________________________________________________________________
<tt>Group</tt><tt> == fineSed
<tt>d50</tt><tt> == 0.000002
      <tt>particle density</tt><tt> == 2650
      <tt>Settling model</tt><tt> == constant
      <tt>settling parameters</tt><tt> == 0.0002 <font color="green"!(m/s)</font>
      <tt>deposition model</tt><tt> == ws0
      <tt>Erosion Model</tt><tt> == Mehta	
      <tt>Erosion parameters</tt><tt> == 0.1, 0.5, 0.5 <font color="green"!Er, taucr, alpha</font>
<tt>End Group
<tt>Group</tt><tt> ==  Gravel
      <tt>d50</tt><tt> == 0.032
      <tt>particle density</tt><tt> == 2650
      <tt>Settling model</tt><tt> == constant
      <tt>settling parameters</tt><tt> == 0.6 <font color="green"!(m/s)</font>
      <tt>Critical stress model</tt><tt> == Soulsby
      <tt>Bed load model</tt><tt> == MPM_Shimizu
      <tt>Bed load parameters</tt><tt> == 8.0, -1 ,1.5	
<tt>End Group
<font color="green"!_________________________________________________________________</font>
<font color="green"! This is required due to adding deposition and settling.</font>
<tt>Material</tt><tt> == 0
      <tt>Layer</tt><tt> == 1
            <tt>dry density</tt><tt> == 1890.,1890	
      <tt>End layer
<tt>end material
<font color="green"!Upstream</font>
<tt>seed particles</tt><tt> == point,  10796.514,8285.014
      <tt>particle groups</tt><tt> == fineSed, Gravel
      <tt>group mass</tt><tt> == 100,100
<tt>end seed
<font color="green"! OUTPUT SETTINGS</font>
<font color="green"!_________________________________________________________________</font>
<tt>output dir</tt><tt> == ..\results\
<tt>output</tt><tt> == ptm_netcdf
      <tt>output groups</tt><tt> == all
      <tt>output parameters</tt><tt> == age, state_age, mass, uvw, uvw_water, depth, water_depth
      <tt>output interval</tt><tt> == 300.
<tt>end output