PARTRACE: Unterschied zwischen den Versionen
Aus BAWiki
imported>BAWiki 3 KKeine Bearbeitungszusammenfassung |
(Arbeitsgruppe SIM) |
||
| (7 dazwischenliegende Versionen von 4 Benutzern werden nicht angezeigt) | |||
| Zeile 2: | Zeile 2: | ||
|name_en=PARTRACE | |name_en=PARTRACE | ||
|name=PARTRACE | |name=PARTRACE | ||
|version= | |version=Oktober 2011 | ||
|version_beschr= | |version_beschr=September 2022 | ||
|stichworte=Analyse von Berechnungsergebnissen (2D)<br /> | |stichworte=Analyse von Berechnungsergebnissen (2D)<br /> | ||
Stochastisches Partikelmodell<br /> | Stochastisches Partikelmodell<br /> | ||
| Zeile 9: | Zeile 9: | ||
turbulente Dispersion<br /> | turbulente Dispersion<br /> | ||
Sinkgeschwindigkeit<br /> | Sinkgeschwindigkeit<br /> | ||
Driftgeschwindigkeit<br /> | |||
Sedimenttransport | Sedimenttransport | ||
|kurzbeschreibung= | |kurzbeschreibung= | ||
| Zeile 16: | Zeile 17: | ||
* Vertikale Konvergenz, bzw. Divergenz der Stromlinien aufgrund der Neigung von Boden und Wasserspiegel | * Vertikale Konvergenz, bzw. Divergenz der Stromlinien aufgrund der Neigung von Boden und Wasserspiegel | ||
* Sinkgeschwindigkeit | * Sinkgeschwindigkeit | ||
* Stochastische Diffusion aufgrund von Turbulenz und Dispersion infolge Tiefenmittelung des Geschwindigkeitsfeldes | * Stochastische Diffusion aufgrund von Turbulenz und Dispersion infolge Tiefenmittelung des Geschwindigkeitsfeldes | ||
* Konstante Driftgeschwindigkeit an der Wasseroberfläche | |||
Am Ende der Simulation wird die zeitliche Abfolge der Ortskoordinaten der Teilchen in die Datei insel.dat geschrieben und kann danach z.B. mit dem Visulisierungsprogramm HVIEW2D als Partikelbahnen im Strömungsgebiet dargestellt werden. | Am Ende der Simulation wird die zeitliche Abfolge der Ortskoordinaten der Teilchen in die Datei insel.dat geschrieben und kann danach z.B. mit dem Visulisierungsprogramm [[HVIEW2D]] als Partikelbahnen im Strömungsgebiet dargestellt werden. | ||
|eingabedateien= | |eingabedateien= | ||
# allgemeine '''Eingabedaten''' (Dateityp [[PARTRACE.DAT|partrace.dat]]) | # allgemeine '''Eingabedaten''' (Dateityp [[PARTRACE.DAT|partrace.dat]]) | ||
# '''Gitternetz''' mit Modelltopographie (Dateityp [[GITTER05.DAT und GITTER05.BIN|gitter05.dat/bin oder untrim_grid.dat]]) | # '''Gitternetz''' mit Modelltopographie (Dateityp [[GITTER05.DAT und GITTER05.BIN|gitter05.dat/bin]] oder [[UNTRIM_GRID.DAT|untrim_grid.dat]]) | ||
# '''Kantenverzeichnis''' des Dreiecksgittersnetzes (Datei des Typs [[FKVZ.BIN|fkvz.bin]]) | # '''Kantenverzeichnis''' des Dreiecksgittersnetzes (Datei des Typs [[FKVZ.BIN|fkvz.bin]]) | ||
# '''Nachbarelementverzeichnis''' des Dreiecksgittersnetzes (Datei des Typs [[FKEZ|fkez.bin]]) | # '''Nachbarelementverzeichnis''' des Dreiecksgittersnetzes (Datei des Typs [[FKEZ.BIN|fkez.bin]]) | ||
# '''synoptische Berechnungsergebnisse''' (Dateien des Typs [[DIRZ.BIN.R|dirz.bin.r]], [[DIRZ.BIN.I|dirz.bin.i]] und [[DIRZ.BIN|dirz.bin]]) | # '''synoptische Berechnungsergebnisse''' (Dateien des Typs [[DIRZ.BIN.R|dirz.bin.r]], [[DIRZ.BIN.I|dirz.bin.i]] und [[DIRZ.BIN|dirz.bin]]) | ||
# (optional) '''Startposition von Partikeln''', z.B. aus einem vorhergehenden Programmlauf (Datei ftn41 auf HP, bzw. fort.41 auf SGI) | # (optional) '''Startposition von Partikeln''', z.B. aus einem vorhergehenden Programmlauf (Datei ftn41 auf HP, bzw. fort.41 auf SGI) | ||
| Zeile 37: | Zeile 39: | ||
Um den Einfluss des tiefengemittelten Geschwindigkeitsfeldes zu erfassen, werden deren Komponenten U und V an den Knoten der Dreieckszelle, in der sich ein Teilchen gerade befindet, jeweils zu Zeiten vor und nach der aktuellen Simulationszeit verfügbar gemacht. Diese Werte werden linear auf die aktuelle Simulationszeit sowie bilinear auf den Teilchenort interpoliert. Die Bewegung des Partikels in den drei Raumrichtungen wird durch drei gewöhnliche gekoppelte Differentialgleichungen in der Zeit beschrieben, in denen die interpolierten Strömungsgeschwindigkeiten, sowie Terme, die die übrigen genannten Effekte charakterisieren, auftreten. Das Differentialgleichungssystem wird mit einem Standard-Runge-Kuttaverfahren (siehe Numerical Recipes) schrittweise gelöst. Daraus ergibt sich schließlich die Bahn des Partikels im Strömungsgebiet. | Um den Einfluss des tiefengemittelten Geschwindigkeitsfeldes zu erfassen, werden deren Komponenten U und V an den Knoten der Dreieckszelle, in der sich ein Teilchen gerade befindet, jeweils zu Zeiten vor und nach der aktuellen Simulationszeit verfügbar gemacht. Diese Werte werden linear auf die aktuelle Simulationszeit sowie bilinear auf den Teilchenort interpoliert. Die Bewegung des Partikels in den drei Raumrichtungen wird durch drei gewöhnliche gekoppelte Differentialgleichungen in der Zeit beschrieben, in denen die interpolierten Strömungsgeschwindigkeiten, sowie Terme, die die übrigen genannten Effekte charakterisieren, auftreten. Das Differentialgleichungssystem wird mit einem Standard-Runge-Kuttaverfahren (siehe Numerical Recipes) schrittweise gelöst. Daraus ergibt sich schließlich die Bahn des Partikels im Strömungsgebiet. | ||
Nähere Informationen zu PARTRACE finden sich unter [[Mathematisches Verfahren PARTRACE|Das Mathematische Verfahren PARTRACE]]. Eine umfassende Dokumentation seiner physikalischen Grundlagen und seiner Anwendungsmöglichkeiten sind der Programmbeschreibung des Partikelverfahrens PARTRACE zu entnehmen. | Nähere Informationen zu PARTRACE finden sich unter der Seite "[[Mathematisches Verfahren PARTRACE|Das Mathematische Verfahren PARTRACE]]". Eine umfassende Dokumentation seiner physikalischen Grundlagen und seiner Anwendungsmöglichkeiten sind der Programmbeschreibung des Partikelverfahrens PARTRACE zu entnehmen. | ||
|preprozessor=DATACONVERT, TELEMAC2D, TR2DIDA, TM2DIDA, UNTRIM | |preprozessor=[[DATACONVERT]], [[TELEMAC-2D|TELEMAC2D]], [[TR2DIDA]], [[TM2DIDA]], [[UNTRIM]] | ||
|postprozessor=HVIEW2D, xmgr | |postprozessor=[[HVIEW2D]], davit, xmgr | ||
|programmiersprache=Fortran, teilw. C | |programmiersprache=Fortran 95, teilw. C | ||
|zus_software= - | |zus_software= - | ||
|kontakt_original= | |kontakt_original=F. Bergemann | ||
|kontakt_pflege=[mailto: | |kontakt_pflege=[mailto:sim.proghome@baw.de Arbeitsgruppe SIM] | ||
|dokumentation=siehe $PROGHOME/examples/partrace/ | |dokumentation=siehe $PROGHOME/examples/partrace/ | ||
Download der Programmbeschreibung als PDF-Datei [http://www.baw.de/downloads/wasserbau/mathematische_verfahren/pdf/partrace1.pdf partrace.pdf] | Download der Programmbeschreibung als PDF-Datei [http://www.baw.de/downloads/wasserbau/mathematische_verfahren/pdf/partrace1.pdf partrace.pdf] | ||
}} | }} | ||
Aktuelle Version vom 6. September 2022, 11:00 Uhr
Basisinformationen
Programm-Name
PARTRACE
Version
Oktober 2011
Beschreibung
September 2022
Stichworte
Analyse von Berechnungsergebnissen (2D)
Stochastisches Partikelmodell
Partikelbahn
turbulente Dispersion
Sinkgeschwindigkeit
Driftgeschwindigkeit
Sedimenttransport
Kurzbeschreibung
PARTRACE ist ein Programm zur Simulation der Bewegung von Partikeln mit wählbaren Eigenschaften in einem zeitabhängigen 2D-tiefengemittelten Strömungsfeld, welches zuvor z.B. mit TRIM-2D, TELEMAC-2D oder UNTRIM (zwei-dimensionale, tiefengemittelte Ergebnisse) berechnet wurde. Es ist ein Werkzeug zur Analyse der Bewegung des Wasserkörpers und des Transportes von Sedimenten z.B. in Tidegewässern. PARTRACE benötigt die Strömungsgrößen in Form von Direktzugriffsdateien des Typs dirz.bin.r, dirz.bin.i und dirz.bin. In der Eingabedatei partrace.dat spezifiziert der Benutzer Zahl und Lage der Partikelquellen im Strömungsgebiet, sowie Zahl, Zeitraum des Aussendens von Teilchen und Eigenschaften der Partikel, wie z.B. Massendichte, Durchmesser, Sinkgeschwindigkeit, Diffusionseigenschaften, usw.. Folgende Einflüsse auf die Partikelbewegung sind berüksichtigt:
- 2D-tiefengemitteltes Geschwindigkeitsfeld
- Vertikale Konvergenz, bzw. Divergenz der Stromlinien aufgrund der Neigung von Boden und Wasserspiegel
- Sinkgeschwindigkeit
- Stochastische Diffusion aufgrund von Turbulenz und Dispersion infolge Tiefenmittelung des Geschwindigkeitsfeldes
- Konstante Driftgeschwindigkeit an der Wasseroberfläche
Am Ende der Simulation wird die zeitliche Abfolge der Ortskoordinaten der Teilchen in die Datei insel.dat geschrieben und kann danach z.B. mit dem Visulisierungsprogramm HVIEW2D als Partikelbahnen im Strömungsgebiet dargestellt werden.
Eingabe-Dateien
- allgemeine Eingabedaten (Dateityp partrace.dat)
- Gitternetz mit Modelltopographie (Dateityp gitter05.dat/bin oder untrim_grid.dat)
- Kantenverzeichnis des Dreiecksgittersnetzes (Datei des Typs fkvz.bin)
- Nachbarelementverzeichnis des Dreiecksgittersnetzes (Datei des Typs fkez.bin)
- synoptische Berechnungsergebnisse (Dateien des Typs dirz.bin.r, dirz.bin.i und dirz.bin)
- (optional) Startposition von Partikeln, z.B. aus einem vorhergehenden Programmlauf (Datei ftn41 auf HP, bzw. fort.41 auf SGI)
Ausgabe-Dateien
- Partikelbahnen (Dateityp insel.dat)
- Für die ersten vier Partikel als Ascii-Dateien:
- z-Koordinaten des Partikels, des Bodens und der freien Oberfläche, letztere am Ort des Teilchens als Funktion der Zeit (Dateien ftn21, ftn22, ftn23, ftn24 auf HP , bzw.
- fort.21, ... auf SGI)
- die drei Geschwindigkeitskomponenten des Partikels als Funktion der Zeit (Dateien ftn31, ftn32, ftn33, ftn34, auf HP, bzw. fort.31, ... auf SGI)
- Druckerdatei mit Informationen zum Programmablauf (Dateityp partrace.sdr)
- (optional) Datei mit Testausgaben (Dateityp partrace.trc)
Methode
Um den Einfluss des tiefengemittelten Geschwindigkeitsfeldes zu erfassen, werden deren Komponenten U und V an den Knoten der Dreieckszelle, in der sich ein Teilchen gerade befindet, jeweils zu Zeiten vor und nach der aktuellen Simulationszeit verfügbar gemacht. Diese Werte werden linear auf die aktuelle Simulationszeit sowie bilinear auf den Teilchenort interpoliert. Die Bewegung des Partikels in den drei Raumrichtungen wird durch drei gewöhnliche gekoppelte Differentialgleichungen in der Zeit beschrieben, in denen die interpolierten Strömungsgeschwindigkeiten, sowie Terme, die die übrigen genannten Effekte charakterisieren, auftreten. Das Differentialgleichungssystem wird mit einem Standard-Runge-Kuttaverfahren (siehe Numerical Recipes) schrittweise gelöst. Daraus ergibt sich schließlich die Bahn des Partikels im Strömungsgebiet.
Nähere Informationen zu PARTRACE finden sich unter der Seite "Das Mathematische Verfahren PARTRACE". Eine umfassende Dokumentation seiner physikalischen Grundlagen und seiner Anwendungsmöglichkeiten sind der Programmbeschreibung des Partikelverfahrens PARTRACE zu entnehmen.
Vorlauf-Programme
DATACONVERT, TELEMAC2D, TR2DIDA, TM2DIDA, UNTRIM
Nachlauf-Programme
HVIEW2D, davit, xmgr
Weitere Informationen
Programmiersprache
Fortran 95, teilw. C
zusätzliche Software
-
Originalversion
F. Bergemann
Programmpflege
Dokumentation/Literatur
siehe $PROGHOME/examples/partrace/
Download der Programmbeschreibung als PDF-Datei partrace.pdf
zurück zu Programmkennblätter