Quelltext der Seite Mathematisches Verfahren UNTRIM2

[[Kategorie:Mathematische Verfahren]]
[[en: Mathematical Model UNTRIM2]]
==Kurzbeschreibung==

Das mathematische Verfahren [[UNTRIM2]] ist eng mit den Verfahren [[Mathematisches Verfahren UNTRIM|UNTRIM]] und [[Mathematisches Verfahren UNTRIM|UNTRIM2007]] verwandt. [[UNTRIM2]] 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. Als wesentliche Neuerung gegenüber den ''klassischen'' UNTRIM-Verfahren kann die Bathymetrie des Modellgebietes mit Hilfe der SubGrid-Technologie wesentlich genauer als bisher, und unabhängig von der Auflösung des Berechnungsgitters, wiedergegeben werden.

==SubGrid-Technologie==

===Berechnungsgitter und SubGrid===

* Berechnungsgitter: 
*: Es wird ein unstrukturiertes orthogonales Gitter (UOG) benutzt. 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.
* SubGrid:
*: Die Modelltopografie wird auf dem SubGrid wiedergegeben. Die Auflösung des SubGrid kann deutlich größer als die des Berechnungsgitters sein.
* Nichtlinearitäten:
*: Das Wasservolumen in einem Berechnungspolygon ist nichtlinear vom Wasserstand abhängig.
*: Die durchströmte Fläche entlang einer Kante kann ebenfalls nichtlinear vom Wasserstand abhängen.

In dem Bild [http://www.baw.de/downloads/wasserbau/mathematische_verfahren/pdf/sub_grid_bathymetry_g_lang.pdf Bathymetrie mit SubGrid] (ca. 0,5 MB) erfolgt die Diskretisierung der Modelltopografie im ''klassischen'' Gitter (links) sowie im Gitter mit SubGrid (rechts). In dem ''klassischen'', aus Dreiecken aufgebauten, unstrukturierten Gitter, variiert die Tiefe in einem Polygon sowie entlang einer Kante nicht. Wird demgegenüber in einem Berechnungsgitter die SubGrid-Technologie zur Diskretisierung der Bathymetrie genutzt, so erkennt man sehr schnell den großen Detailreichtum dieses Wattgebiets. Dies ist darauf zurückzuführen, dass bei Verwenden der SubGrid-Technologie die Tiefe in jedem Polygon und entlang einer jeden Kante veränderlich sein kann.

===Berechnungsergebnisse===

Die Ergebnisse für z. B. Wasserstand und Strömung liegen auf dem (groben) Berechnungsgitter und nicht auf dem (feineren) SubGrid vor. Typischer Weise entspricht der Wasserstand dem mittleren Wasserstand innerhalb eines Polygons, und die Strömungsgeschwindigkeit gibt für jede Kante jeweils den Mittelwert für die aktuell durchströmte Fläche wieder.

In dem Bild [http://www.baw.de/downloads/wasserbau/mathematische_verfahren/pdf/sub_grid_flow_g_lang.pdf Strömung auf Gitter mit SubGrid] (ca. 0,25 MB) wird links die Strömung für das ''klassische'', und rechts für das Gitter mit SubGrid dargestellt. Man erkennt, dass im klassischen Fall jedes Polygon nur gänzlich von Wasser bedeckt oder vollkommen trockengefallen sein kann. In dem Gitter mit SubGrid kann ein Polygon (und eine Kante) hingegen auch ''teilweise'' von Wasser bedeckt sein. Damit gelingt eine naturähnlichere Simulation der Überflutung von Wattgebieten.

===Vorteile===

# Die Topografie kann unabhängig von der Auflösung des Berechnungsgitters beschrieben werden.
# Die Genauigkeit, mit der die Topografie wiedergegeben werden kann, wird prinzipiell nur durch die Genauigkeit der gemessenen Daten begrenzt.
# Für jede Wasserspiegellage entspricht das Wasservolumen in dem Modell dem in der Natur vorhandenen Wasservolumen. Gleiches gilt für die durchströmten Querschnitte.
# Überflutete und trocken gefallene Gebiete können präzise beschrieben werden.
# Die gleiche Genauigkeit hinsichtlich der Diskretisierung der Bathymetrie wird mit SubGrid zu deutlich geringeren CPU-Kosten erzielt, als dies mit einem ''klassischen'' Berechnungsgitter, ohne SubGrid-Technologie, möglich wäre.

==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. (2008), ''A high-resolution wetting and drying algorithm for free-surface hydrodynamics'', International Journal for Numerical Methods in Fluids, Volume 60, Issue 4, pages 391 - 408. See [http://onlinelibrary.wiley.com/doi/10.1002/fld.1896/abstract abstract].
# Casulli, V. and Stelling, G. S. (2010), ''Semi-implicit sub grid modelling of three-dimensional free-surface flows'', International Journal for Numerical Methods in Fluids, online available, advance of print. See [http://onlinelibrary.wiley.com/doi/10.1002/fld.2361/abstract abstract].

==Vorträge==

* [http://www.baw.de/downloads/wasserbau/mathematische_verfahren/pdf/a_395_5_03_70150_tn12_2010_11_30_UnTRIM_SubGrid_g_lang.pdf UNTRIM<sup>2</sup> (sub grid) at BAW Hamburg], Vortragsfolien vom 11. November 2010, ca. 15 MB.

==Validierungsdokument==

Derzeit nicht verfügbar.

==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. 1.41 MB) [http://www.baw.de/downloads/wasserbau/mathematische_verfahren/pdf/Simulationsverfahren_Kueste_untrim2-ui.pdf UNTRIM<sup>2</sup> (version August 2009) user interface description] (in Englisch)

==MPI-Parallelisierung==

Derzeit nicht realisiert.

==BAW-spezifische Informationen==

===Gitternetzgenerierung===

Für das Erzeugen des von UNTRIM2 benötigten unstrukturierten orthogonalen Gitters mit SubGrid-Bathymetrie 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===

Das mathematische Verfahren UNTRIM2 ist derzeit nur unvollständig in die Softwareumgebung der BAW eingebunden. Weitere Detailinformationen finden sich auf dem [[UNTRIM2|UNTRIM2-Programmkennblatt]].

===Grafische Darstellung der Berechnungsergebnisse===

Für die grafische Darstellung der von UNTRIM2 erzeugten Ergebnisse stehen mehrere Methoden bereit. Zu den wichtigsten zählen,

* SGHVIEW2D, 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. Man beachte jedoch, dass eine Analyse derzeit nur für die auf Profilen sowie an Einzelpositionen vorliegenden Berechnungsergebnisse durchgeführt werden kann.

===Kopplung mit weiteren Sub-Modellen===

Derzeit nicht realisiert.

----
zurück zu [[Modellverfahren für den Küstenbereich und Ästuare]]

----
[[Strukturübersicht]]