PARTRACE: Difference between revisions
From BAWiki
imported>BAWiki 3 mNo edit summary |
Günther Lang (talk | contribs) (working group SIM) |
||
| (4 intermediate revisions by 4 users not shown) | |||
| Line 2: | Line 2: | ||
|name_de=PARTRACE | |name_de=PARTRACE | ||
|name=PARTRACE | |name=PARTRACE | ||
|version= | |version=October 2011 | ||
|version_descr= | |version_descr=September 2022 | ||
|catchwords=analysis of numerically calculated results (2D)<br /> | |catchwords=analysis of numerically calculated results (2D)<br /> | ||
stochastic particle model<br /> | stochastic particle model<br /> | ||
| Line 9: | Line 9: | ||
turbulent dispersion<br /> | turbulent dispersion<br /> | ||
sinking velocity<br /> | sinking velocity<br /> | ||
drift velocity<br /> | |||
sediment transport | sediment transport | ||
|shortdescription= | |shortdescription= | ||
PARTRACE simulates the movement of particles with user-defined properties in a time-dependent 2D-depth-averaged flowfield, which has been computed e.g. with [[TRIM-2D]], [[TELEMAC-2D]] | PARTRACE simulates the movement of particles with user-defined properties in a time-dependent 2D-depth-averaged flowfield, which has been computed e.g. with [[TRIM-2D]], [[TELEMAC-2D]] or [[UNTRIM]] (two-dimensional, depth averaged results only). It is a tool for analyzing the motion of the water-body and the transport of sediments in estuarine flows. PARTRACE needs the flowfield data as direct access files of the type [[DIRZ.BIN.R|dirz.bin.r]], [[DIRZ.BIN.I|dirz.bin.i]] and [[DIRZ.BIN|dirz.bin]]. In the input file [[PARTRACE.DAT|partrace.dat]] the user specifies the number of particle-emitting sources, their location as well as number of particles to be tracked, seeding time intervall, and particle properties such as mass density, diameter, settling velocity, diffusion properties, ect.. The following influences on the particle motion are taken into account: | ||
* 2D-depth-averaged velocity field | * 2D-depth-averaged velocity field | ||
| Line 17: | Line 18: | ||
* settling velocity | * settling velocity | ||
* random diffusion due to turbulence and dispersion due to depth-averaging of the velocity field | * random diffusion due to turbulence and dispersion due to depth-averaging of the velocity field | ||
* constant particle drift velocity at the water surface | |||
To the end of the simulation the history of particle locations is written to file insel.dat and can be visualized e.g. with HVIEW2D as pathlines in the domain. | To the end of the simulation the history of particle locations is written to file insel.dat and can be visualized e.g. with HVIEW2D as pathlines in the domain. | ||
|inputfiles= | |inputfiles= | ||
# general '''input data''' (filetype [[PARTRACE.DAT|partrace.dat]]) | # general '''input data''' (filetype [[PARTRACE.DAT|partrace.dat]]) | ||
# '''grid''' for 2D/3D-data (filetype [[GITTER05.DAT | # '''grid''' for 2D/3D-data (filetype [[GITTER05.DAT and GITTER05.BIN|gitter05.dat/bin]] or filetype [[UNTRIM_GRID.DAT|untrim_grid.dat]]) | ||
# '''egde list''' of triangular grid (file type [[FKVZ.BIN|fkvz.bin]]) | # '''egde list''' of triangular grid (file type [[FKVZ.BIN|fkvz.bin]]) | ||
# '''neighbor element list''' of triangular grid (file type [[FKVZ.BIN|fkez.bin]]) | # '''neighbor element list''' of triangular grid (file type [[FKVZ.BIN|fkez.bin]]) | ||
# '''synoptic flow field''' results (filetype [[DIRZ.BIN.R|dirz.bin.r]], | # '''synoptic flow field''' results (filetype [[DIRZ.BIN.R|dirz.bin.r]], [[DIRZ.BIN.I|dirz.bin.i]], and [[DIRZ.BIN|dirz.bin]]) | ||
# (optional) initial location of particles, e.g. from a previous run (file name ftn41 on HP, or fort.41 on SGI) | # (optional) initial location of particles, e.g. from a previous run (file name ftn41 on HP, or fort.41 on SGI) | ||
|outputfiles= | |outputfiles= | ||
| Line 36: | Line 38: | ||
|methodology= | |methodology= | ||
To capture the effect of the depth-averaged velocity field, its components U and V at the 3 nodes of the triangle where the particle is currently located in, at time instances before and after the current simulation time, are made available. These data are linearly interpolated to the current simulation time and bilinearly interpolated onto the location of the particle. The motion of the particle in the 3 spatial directions is described by 3 coupled ordinary differential equations (ODEs) in time. They contain the interpolated velocity components and terms that characterize the other physical effects listed above. The system of ODEs is solved stepwise using a standard-Runge-Kutta method (see Numerical Recipes). This finally yields the path of a particle in the flowfield. | To capture the effect of the depth-averaged velocity field, its components U and V at the 3 nodes of the triangle where the particle is currently located in, at time instances before and after the current simulation time, are made available. These data are linearly interpolated to the current simulation time and bilinearly interpolated onto the location of the particle. The motion of the particle in the 3 spatial directions is described by 3 coupled ordinary differential equations (ODEs) in time. They contain the interpolated velocity components and terms that characterize the other physical effects listed above. The system of ODEs is solved stepwise using a standard-Runge-Kutta method (see Numerical Recipes). This finally yields the path of a particle in the flowfield. | ||
More information about PARTRACE is given in [[Mathematical Model PARTRACE|The Mathematical Model PARTRACE]]. A thorough documentation of the physical modelling used its features can be found in the Programmbeschreibung des Partikelverfahrens PARTRACE (currently in german only, sorry!). | More information about PARTRACE is given in [[Mathematical Model PARTRACE|The Mathematical Model PARTRACE]]. A thorough documentation of the physical modelling used its features can be found in the Programmbeschreibung des Partikelverfahrens PARTRACE (currently in german only, sorry!). | ||
|preprocessor=[[DATACONVERT]], [[ | |preprocessor=[[DATACONVERT]], [[TELEMAC-2D]], [[TR2DIDA]], [[TM2DIDA]], [[UNTRIM]] | ||
|postprocessor=[[HVIEW2D]], xmgr | |postprocessor=[[HVIEW2D]], davit, xmgr | ||
|language= | |language=Fortran 95, partly C | ||
|add_software=- | |add_software=- | ||
|contact_original= | |contact_original=F. Bergemann | ||
|contact_maintenance=[mailto: | |contact_maintenance=[mailto:sim.proghome@baw.de Working group SIM] | ||
|documentation= | |documentation= | ||
please refer to $PROGHOME/examples/partrace/ | |||
Download the program description (in german language only) as PDF-file [http://www.baw.de/downloads/wasserbau/mathematische_verfahren/pdf/partrace1.pdf partrace.pdf] | Download the program description (in german language only) as PDF-file [http://www.baw.de/downloads/wasserbau/mathematische_verfahren/pdf/partrace1.pdf partrace.pdf] | ||
}} | }} | ||
Latest revision as of 11:02, 6 September 2022
Basic Information
Name of Program
PARTRACE
Version-Date
October 2011
Description-Date
September 2022
Catchwords
analysis of numerically calculated results (2D)
stochastic particle model
particle path
turbulent dispersion
sinking velocity
drift velocity
sediment transport
Short Description of Functionality
PARTRACE simulates the movement of particles with user-defined properties in a time-dependent 2D-depth-averaged flowfield, which has been computed e.g. with TRIM-2D, TELEMAC-2D or UNTRIM (two-dimensional, depth averaged results only). It is a tool for analyzing the motion of the water-body and the transport of sediments in estuarine flows. PARTRACE needs the flowfield data as direct access files of the type dirz.bin.r, dirz.bin.i and dirz.bin. In the input file partrace.dat the user specifies the number of particle-emitting sources, their location as well as number of particles to be tracked, seeding time intervall, and particle properties such as mass density, diameter, settling velocity, diffusion properties, ect.. The following influences on the particle motion are taken into account:
- 2D-depth-averaged velocity field
- vertical convergence or divergence of streamlines due to inclination of the seabed and of the water-level surface
- settling velocity
- random diffusion due to turbulence and dispersion due to depth-averaging of the velocity field
- constant particle drift velocity at the water surface
To the end of the simulation the history of particle locations is written to file insel.dat and can be visualized e.g. with HVIEW2D as pathlines in the domain.
Input-Files
- general input data (filetype partrace.dat)
- grid for 2D/3D-data (filetype gitter05.dat/bin or filetype untrim_grid.dat)
- egde list of triangular grid (file type fkvz.bin)
- neighbor element list of triangular grid (file type fkez.bin)
- synoptic flow field results (filetype dirz.bin.r, dirz.bin.i, and dirz.bin)
- (optional) initial location of particles, e.g. from a previous run (file name ftn41 on HP, or fort.41 on SGI)
Output-Files
- particle paths (file insel.dat)
- for the first 4 particles, ascii-files containing:
- z-coordinates of particle, of the seabed and of the water-level surface, the ladder at the paticle location, respectively, as function of time (file names ftn21, ftn22,
- ftn23, ftn24 on HP, or fort.21, ... on SGI)
- the 3 velocity components of the particle as a function of time (file names ftn31, ftn32, ftn33, ftn34 on HP, or fort.31, ... on SGI)
- printer file (filetype partrace.sdr)
- (optional) trace of program execution (filetype partrace.trc)
Methodology
To capture the effect of the depth-averaged velocity field, its components U and V at the 3 nodes of the triangle where the particle is currently located in, at time instances before and after the current simulation time, are made available. These data are linearly interpolated to the current simulation time and bilinearly interpolated onto the location of the particle. The motion of the particle in the 3 spatial directions is described by 3 coupled ordinary differential equations (ODEs) in time. They contain the interpolated velocity components and terms that characterize the other physical effects listed above. The system of ODEs is solved stepwise using a standard-Runge-Kutta method (see Numerical Recipes). This finally yields the path of a particle in the flowfield. More information about PARTRACE is given in The Mathematical Model PARTRACE. A thorough documentation of the physical modelling used its features can be found in the Programmbeschreibung des Partikelverfahrens PARTRACE (currently in german only, sorry!).
Program(s) to run before this Program
DATACONVERT, TELEMAC-2D, TR2DIDA, TM2DIDA, UNTRIM
Program(s) to run after this Program
HVIEW2D, davit, xmgr
Additional Information
Language
Fortran 95, partly C
Additional software
-
Original Version
F. Bergemann
Maintenance
Documentation/Literature
please refer to $PROGHOME/examples/partrace/
Download the program description (in german language only) as PDF-file partrace.pdf
back to Program Descriptions