Mathematisches Verfahren UNTRIM: Unterschied zwischen den Versionen
Aus BAWiki
imported>BAWiki Glossar K (Edited by Wikibot) |
(Die LinkTitles-Erweiterung hat automatisch Links zu anderen Seiten hinzugefügt (<a target="_blank" rel="nofollow noreferrer noopener" class="external free" href="https://github.com/bovender/LinkTitles">https://github.com/bovender/LinkTitles</a>).) |
||
| Zeile 3: | Zeile 3: | ||
==Kurzbeschreibung== | ==Kurzbeschreibung== | ||
Das [[mathematische Verfahren]] [[UNTRIM]] wurde von Herrn Prof. Vincenzo Casulli (Universität Trient, Italien) entwickelt. [[UNTRIM]] ist ein semi-implizites Finite-Differenzen (-Volumen) Verfahren zur numerischen Lösung der dreidimensionalen Flachwassergleichungen sowie der dreidimensionalen Transportgleichung für Salz, Wärme, Schwebstoffe sowie suspendierte Sedimente. | Das [[mathematische Verfahren]] [[UNTRIM]] wurde von Herrn Prof. Vincenzo Casulli (Universität Trient, Italien) entwickelt. [[UNTRIM]] ist ein semi-implizites Finite-Differenzen (-Volumen) Verfahren zur numerischen Lösung der dreidimensionalen Flachwassergleichungen sowie der dreidimensionalen Transportgleichung für Salz, Wärme, [[Schwebstoffe]] sowie suspendierte Sedimente. | ||
UNTRIM arbeitet auf einem unstrukturierten orthogonalen Gitter (UOG). Hierbei wird das Lösungsgebiet von einer endlichen Anzahl konvexer Polygone überlappungsfrei überdeckt. Ein Gitter entspricht genau dann einem UOG, wenn innerhalb eines jeden Polygons ein Punkt (Zentrum) bestimmt werden kann, so dass jede Verbindungslinie zu einem Zentrum eines Nachbarpolygons die gemeinsame Seite der Polygone senkrecht schneidet. | [[UNTRIM]] arbeitet auf einem unstrukturierten orthogonalen Gitter (UOG). Hierbei wird das Lösungsgebiet von einer endlichen Anzahl konvexer Polygone überlappungsfrei überdeckt. Ein Gitter entspricht genau dann einem UOG, wenn innerhalb eines jeden Polygons ein Punkt (Zentrum) bestimmt werden kann, so dass jede Verbindungslinie zu einem Zentrum eines Nachbarpolygons die gemeinsame Seite der Polygone senkrecht schneidet. | ||
==Physikalische Prozesse== | ==Physikalische Prozesse== | ||
| Zeile 29: | Zeile 29: | ||
** horizontale turbulente Diffusivität (lokal isotrop, zeit- und ortsvariabel) | ** horizontale turbulente Diffusivität (lokal isotrop, zeit- und ortsvariabel) | ||
** turbulente Diffusivität in Vertikalrichtung unter Berücksichtigung der vertikalen Dichteschichtung | ** turbulente Diffusivität in Vertikalrichtung unter Berücksichtigung der vertikalen Dichteschichtung | ||
** Sinkgeschwindigkeit, Deposition und Erosion (bei Schwebstoffen) | ** [[Sinkgeschwindigkeit]], Deposition und [[Erosion]] (bei Schwebstoffen) | ||
** Wärmeaustausch mit der Atmosphäre und der Gewässersohle (bei Wärmetransport) | ** Wärmeaustausch mit der Atmosphäre und der [[Gewässersohle]] (bei Wärmetransport) | ||
** Quellen und Senken | ** Quellen und Senken | ||
** Senken mit unmittelbarer Wiedereinleitung an einem anderen Ort, mit der Möglichkeit zur Abwandlung z.B. der Einleitungstemperatur sowie des Einleitungssalzgehalts gegenüber den entsprechenden Entnahmewerten. | ** Senken mit unmittelbarer Wiedereinleitung an einem anderen Ort, mit der Möglichkeit zur Abwandlung z.B. der Einleitungstemperatur sowie des Einleitungssalzgehalts gegenüber den entsprechenden Entnahmewerten. | ||
| Zeile 38: | Zeile 38: | ||
* Wasserspiegelauslenkung der freien Oberfläche | * Wasserspiegelauslenkung der freien Oberfläche | ||
* Strömungsgeschwindigkeit | * Strömungsgeschwindigkeit | ||
* Tracerkonzentration (z.B. Salzgehalt, Temperatur, Schwebstoffgehalt) | * Tracerkonzentration (z.B. Salzgehalt, Temperatur, [[Schwebstoffgehalt]]) | ||
* hydrodynamischer Druck | * hydrodynamischer Druck | ||
Anmerkung: Bei zwei-dimensionaler Modellrechnung entsprechen die Ergebnisse den über die Wassertiefe gemittelten Werten. | Anmerkung: Bei zwei-dimensionaler Modellrechnung entsprechen die Ergebnisse den über die [[Wassertiefe]] gemittelten Werten. | ||
==Veröffentlichungen== | ==Veröffentlichungen== | ||
| Zeile 52: | Zeile 52: | ||
==Weitere Anwender== | ==Weitere Anwender== | ||
Verschiedene Versionen des von Herrn Prof. Casulli bereit gestellten Berechnungskerns UnTRIM werden von folgenden Institutionen verwendet: | Verschiedene Versionen des von Herrn Prof. Casulli bereit gestellten Berechnungskerns [[UNTRIM|UnTRIM]] werden von folgenden Institutionen verwendet: | ||
* [http://www.deltamodeling.com/ ''Delta Modeling Associates, Inc.'', San Francisco, California, USA] | * [http://www.deltamodeling.com/ ''Delta Modeling Associates, Inc.'', San Francisco, California, USA] | ||
* [http://www.rma.com/ ''Resource Management Associates'', Fairfield, California, USA] | * [http://www.rma.com/ ''Resource Management Associates'', Fairfield, California, USA] | ||
| Zeile 71: | Zeile 71: | ||
# Physikalisches System, | # Physikalisches System, | ||
# Modellfunktionalität | # Modellfunktionalität | ||
# konzeptionelles Modell, | # konzeptionelles [[Modell]], | ||
# algorithmische Implementierung, | # algorithmische Implementierung, | ||
# Software-Implementierung, | # Software-Implementierung, | ||
| Zeile 103: | Zeile 103: | ||
==MPI-Parallelisierung== | ==MPI-Parallelisierung== | ||
Der Berechnungskern des mathematischen Verfahrens UNTRIM wurde von [mailto:jacek.jankowski@baw.de Jacek Jankowski] mit MPI parallelisiert. Eine ausführliche Beschreibung ist in dem technischen Bericht | Der Berechnungskern des mathematischen Verfahrens [[UNTRIM]] wurde von [mailto:jacek.jankowski@baw.de Jacek Jankowski] mit MPI parallelisiert. Eine ausführliche Beschreibung ist in dem technischen Bericht | ||
(ca. 1.1 MB) [http://www.baw.de/downloads/wasserbau/mathematische_verfahren/pdf/manual_mpi_untrim.pdf MPI Version Manual] (in Englisch) enthalten. Diese Version wird derzeit hausintern vor allem in der Abteilung Wasserbau im Binnenbereich eingesetzt. | (ca. 1.1 MB) [http://www.baw.de/downloads/wasserbau/mathematische_verfahren/pdf/manual_mpi_untrim.pdf MPI Version Manual] (in Englisch) enthalten. Diese Version wird derzeit hausintern vor allem in der Abteilung Wasserbau im Binnenbereich eingesetzt. | ||
| Zeile 110: | Zeile 110: | ||
===Gitternetzgenerierung=== | ===Gitternetzgenerierung=== | ||
Für das Erzeugen des von UNTRIM benötigten unstrukturierten orthogonalen Gitters wird der Gitternetz-Generator [[JANET]] eingesetzt, eine Entwicklung der Firma [http://www.smileconsult.de/ SmileConsult]. Weitere Informationen zur Einbindung von JANET bei der BAW finden sich auf dem entsprechenden [[JANET|JANET-Programmkennblatt]]. | Für das Erzeugen des von [[UNTRIM]] benötigten unstrukturierten orthogonalen Gitters wird der Gitternetz-Generator [[JANET]] eingesetzt, eine Entwicklung der Firma [http://www.smileconsult.de/ SmileConsult]. Weitere Informationen zur Einbindung von [[JANET]] bei der BAW finden sich auf dem entsprechenden [[JANET|JANET-Programmkennblatt]]. | ||
===Simulation=== | ===Simulation=== | ||
Das mathematische Verfahren UNTRIM ist vollständig in die Softwareumgebung der BAW eingebunden. Weitere Detailinformationen finden sich auf dem [[UNTRIM|UNTRIM-Programmkennblatt]]. | Das [[mathematische Verfahren]] [[UNTRIM]] ist vollständig in die Softwareumgebung der BAW eingebunden. Weitere Detailinformationen finden sich auf dem [[UNTRIM|UNTRIM-Programmkennblatt]]. | ||
===Grafische Darstellung der Berechnungsergebnisse=== | ===Grafische Darstellung der Berechnungsergebnisse=== | ||
Für die grafische Darstellung der von UNTRIM erzeugten Ergebnisse stehen mehrere Methoden bereit. Zu den wichtigsten zählen, | Für die grafische Darstellung der von [[UNTRIM]] erzeugten Ergebnisse stehen mehrere Methoden bereit. Zu den wichtigsten zählen, | ||
* [[HVIEW2D]], für flächenhaft vorliegende Daten, | * [[HVIEW2D]], für flächenhaft vorliegende Daten, | ||
| Zeile 130: | Zeile 130: | ||
===Kopplung mit weiteren Sub-Modellen=== | ===Kopplung mit weiteren Sub-Modellen=== | ||
Das mathematische Verfahren UNTRIM ist ein seiner Einbindung bei der BAW mit den nachfolgend genannten Sub-Modellen verknüpft: | Das [[mathematische Verfahren]] [[UNTRIM]] ist ein seiner Einbindung bei der BAW mit den nachfolgend genannten Sub-Modellen verknüpft: | ||
# Spektrales Seegangsmodell ''K-MODELL'' (siehe [[K_MODEL.DAT|k_model.dat]]); | # Spektrales Seegangsmodell ''K-[[Modell|MODELL]]'' (siehe [[K_MODEL.DAT|k_model.dat]]); | ||
# Sedimentologisches Modell ''SEDIMORPH'' (siehe [[SEDIMORPH.DAT|sedimorph.dat]]). | # Sedimentologisches [[Modell]] ''SEDIMORPH'' (siehe [[SEDIMORPH.DAT|sedimorph.dat]]). | ||
Diese Sub-Modelle können in direkter Kopplung mit UNTRIM verwendet werden. | Diese Sub-Modelle können in direkter Kopplung mit [[UNTRIM]] verwendet werden. | ||
---- | ---- | ||
Aktuelle Version vom 21. Oktober 2022, 08:56 Uhr
Kurzbeschreibung
Das mathematische Verfahren UNTRIM wurde von Herrn Prof. Vincenzo Casulli (Universität Trient, Italien) entwickelt. UNTRIM ist ein semi-implizites Finite-Differenzen (-Volumen) Verfahren zur numerischen Lösung der dreidimensionalen Flachwassergleichungen sowie der dreidimensionalen Transportgleichung für Salz, Wärme, Schwebstoffe sowie suspendierte Sedimente.
UNTRIM arbeitet auf einem unstrukturierten orthogonalen Gitter (UOG). Hierbei wird das Lösungsgebiet von einer endlichen Anzahl konvexer Polygone überlappungsfrei überdeckt. Ein Gitter entspricht genau dann einem UOG, wenn innerhalb eines jeden Polygons ein Punkt (Zentrum) bestimmt werden kann, so dass jede Verbindungslinie zu einem Zentrum eines Nachbarpolygons die gemeinsame Seite der Polygone senkrecht schneidet.
Physikalische Prozesse
- reynoldsgemittelte Navier-Stokes-Gleichung (RANS)
- lokale Beschleunigung (Massenträgheit)
- advektive Beschleunigung
- Coriolisbeschleunigung
- barotroper Druckgradient
- barokliner Druckgradient
- hydrostatische oder nicht-hydrostatische Druckverteilung
- horizontale turbulente Viskosität (lokal isotrop, zeit- und ortsvariabel)
- turbulente Viskosität in Vertikalrichtung unter Berücksichtigung der vertikalen Dichteschichtung
- Bodenreibung
- Impulseintrag durch den Wind
- Quellen und Senken
- horizontale seegangsinduzierte Beschleunigung (durch Radiation Stress)
- Transport konservativer Tracer
- lokale Veränderung
- advektiver Transport durch die Strömung
- optionaler flux limiter : Minmod, van Leer oder Superbee
- horizontale turbulente Diffusivität (lokal isotrop, zeit- und ortsvariabel)
- turbulente Diffusivität in Vertikalrichtung unter Berücksichtigung der vertikalen Dichteschichtung
- Sinkgeschwindigkeit, Deposition und Erosion (bei Schwebstoffen)
- Wärmeaustausch mit der Atmosphäre und der Gewässersohle (bei Wärmetransport)
- Quellen und Senken
- Senken mit unmittelbarer Wiedereinleitung an einem anderen Ort, mit der Möglichkeit zur Abwandlung z.B. der Einleitungstemperatur sowie des Einleitungssalzgehalts gegenüber den entsprechenden Entnahmewerten.
Berechnungsergebnisse
- Wasserspiegelauslenkung der freien Oberfläche
- Strömungsgeschwindigkeit
- Tracerkonzentration (z.B. Salzgehalt, Temperatur, Schwebstoffgehalt)
- hydrodynamischer Druck
Anmerkung: Bei zwei-dimensionaler Modellrechnung entsprechen die Ergebnisse den über die Wassertiefe gemittelten Werten.
Veröffentlichungen
- Casulli, V. and Zanolli, P. (1998), A Three-Dimensional Semi-Implicit Algorithm for Environmental Flows on Unstructured Grids, Proc. of Conf. on Num. Methods for Fluid Dynamics, University of Oxford.
- Casulli, V. (1999), A Semi-Implicit Finite Difference Method for Non-Hydrostatic, Free-Surface Flows, International Journal for Numerical Methods in Fluids, 30: 425 - 440.
- Casulli, V. and R.A. Walters (2000), An unstructured grid, three-dimensional model based on the shallow water equations. International Journal for Numerical Methods in Fluids, Volume 32, Issue 3, pages 331 - 348.
- Casulli, V. and Zanolli, P. (2002), Semi-Implicit Numerical Modelling of Non-Hydrostatic Free-Surface Flows for Environmental Problems, Mathematical and Computer Modelling, 36: 1131 - 1149.
- Casulli, V. and Zanolli, P. (2004), High Resolution Methods for Multidimensional Advection-Diffusion Problems in Free-Surface Hydrodynamics, Ocean Modelling.
Weitere Anwender
Verschiedene Versionen des von Herrn Prof. Casulli bereit gestellten Berechnungskerns UnTRIM werden von folgenden Institutionen verwendet:
- Delta Modeling Associates, Inc., San Francisco, California, USA
- Resource Management Associates, Fairfield, California, USA
- Stanford University, Stanford, California, USA
- University of California, Davis, California, USA
- Virginia Institute of Marine Science, Gloucester Point, Maryland, USA
Weitere Software
Die von Herrn Prof. Casulli entwickelten und in mehreren Publikationen veröffentlichten Konzepte haben die Entwicklung folgender Modellverfahren inspiriert:
- D-Flow Flexible Mesh
- ELCIRC, SELFE
- SUNTANS
Validierungsdokument
Das Validierungsdokument enthält neben einer allgemeinen Einführung in das Simulationsverfahren detaillierte Informationen zu folgenden Themen:
- Physikalisches System,
- Modellfunktionalität
- konzeptionelles Modell,
- algorithmische Implementierung,
- Software-Implementierung,
- Validierungsstudien, und
- Literatur.
Die PDF-Version der englischen Dokument-Fassung kann frei heruntergeladen werden:
- (ca. 1.2 MB) UNTRIM standard validation document (in Englisch)
Benutzer-Schnittstellen-Beschreibung
Dieses Dokument enthält eine detaillierte Beschreibung aller dem Anwender des Simulationsverfahrens zur Verfügung stehenden Schnittstellen-Funktionen. Folgende Schnittstellen-Themen werden behandelt:
- Setzen von Daten (set-interfaces),
- Holen von Daten (get-interfaces),
- Testen des Gitters und der Genauigkeit der iterativ berechneten Lösung (check-routines),
- vom Berechnungskern gerufene externe Routinen (user-interface-routines) zum,
- Definieren der Namen der Dateien mit den Standard-Eingabedaten, zum
- Setzen des Anfangszustands, zum
- Setzen der Randwerte für jeden Berechnungszeitschritt, und zum
- Holen der Berechnungsergebnisse.
- tabellarische Aufzählung aller SET- und GET-Schnittstellen,
- Beispiele zu den Dateien mit den Standard-Eingabedaten.
Die PDF-Version der englischen Dokument-Fassung kann frei heruntergeladen werden:
- (ca. 0.65 MB) UNTRIM user interface description (in Englisch)
- (ca. 1.41 MB) UNTRIM2 (version August 2009) user interface description (in Englisch)
MPI-Parallelisierung
Der Berechnungskern des mathematischen Verfahrens UNTRIM wurde von Jacek Jankowski mit MPI parallelisiert. Eine ausführliche Beschreibung ist in dem technischen Bericht (ca. 1.1 MB) MPI Version Manual (in Englisch) enthalten. Diese Version wird derzeit hausintern vor allem in der Abteilung Wasserbau im Binnenbereich eingesetzt.
BAW-spezifische Informationen
Gitternetzgenerierung
Für das Erzeugen des von UNTRIM benötigten unstrukturierten orthogonalen Gitters wird der Gitternetz-Generator JANET eingesetzt, eine Entwicklung der Firma SmileConsult. Weitere Informationen zur Einbindung von JANET bei der BAW finden sich auf dem entsprechenden JANET-Programmkennblatt.
Simulation
Das mathematische Verfahren UNTRIM ist vollständig in die Softwareumgebung der BAW eingebunden. Weitere Detailinformationen finden sich auf dem UNTRIM-Programmkennblatt.
Grafische Darstellung der Berechnungsergebnisse
Für die grafische Darstellung der von UNTRIM erzeugten Ergebnisse stehen mehrere Methoden bereit. Zu den wichtigsten zählen,
- HVIEW2D, für flächenhaft vorliegende Daten,
- VVIEW2D und/oder LQ2PRO, für auf Quer- und Längsschnitten vorliegende Daten, sowie
- GVIEW2D, für an Einzelpositionen vorliegende Daten.
Analyse der Berechnungsergebnisse
Für eine Analyse der Berechnungsergebnisse stehen verschiedene Methoden für unterschiedlichste Fragestellungen bereit.
Kopplung mit weiteren Sub-Modellen
Das mathematische Verfahren UNTRIM ist ein seiner Einbindung bei der BAW mit den nachfolgend genannten Sub-Modellen verknüpft:
- Spektrales Seegangsmodell K-MODELL (siehe k_model.dat);
- Sedimentologisches Modell SEDIMORPH (siehe sedimorph.dat).
Diese Sub-Modelle können in direkter Kopplung mit UNTRIM verwendet werden.
zurück zu Modellverfahren für den Küstenbereich und Ästuare