Harmonische Analyse der Strömung: Unterschied zwischen den Versionen
Aus BAWiki
imported>Lang Guenther (→Motivation: Text und Bild ergänzt) |
(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>).) |
||
(14 dazwischenliegende Versionen von 5 Benutzern werden nicht angezeigt) | |||
Zeile 4: | Zeile 4: | ||
=Motivation= | =Motivation= | ||
Die Strömung kann für jeden Ort in einem Ästuar in einen durch die Gezeiten bedingten (periodischen) Anteil und den Rest (aperiodischer Rest) zerlegt werden. | Die Strömung kann für jeden Ort in einem [[Ästuar]] in einen durch die Gezeiten bedingten (periodischen) Anteil und den Rest (aperiodischer Rest) zerlegt werden. | ||
Der periodische Anteil wird dabei von den Partialtiden, also den astronomischen Tiden und den Flachwassertiden dominiert. | Der periodische Anteil wird dabei von den Partialtiden, also den astronomischen Tiden und den Flachwassertiden dominiert. | ||
[[Datei:Tidal_current_ellipse_definition. | [[Datei:Tidal_current_ellipse_definition.png|200px|thumb|right|Definition Tide-Ellipse]] | ||
Für jede Partialtide ergibt sich eine Tide-Ellipse, welcher von der Spitze des Strömungsvektors einmal in jeder Periode (der Partialtide) durchlaufen wird (siehe nebenstehendes Bild). Charakteristische Parameter sind | Für jede [[Partialtide]] ergibt sich eine [[Tide]]-Ellipse, welcher von der Spitze des Strömungsvektors einmal in jeder Periode (der [[Partialtide]]) durchlaufen wird | ||
(siehe nebenstehendes Bild). Charakteristische Parameter sind der Wert der großen sowie der kleinen Halbachse der [[Tide]]-Ellipse, | |||
deren Neigunggswinkel gegen z. B. Ost, sowie deren Exzentrizität. | |||
Je nach Exzentrizität degeneriert die Ellipse entweder zu einem Kreis (Exzentrizität +1.0 oder -1.0) oder einer Geraden (Exzentrizität 0.0). | |||
Mit diesen Parametern der Partialtiden können die tidebedingten Anteile an der lokalen Strömungsgeschwindigkeit quantifiziert werden. | |||
Der aperiodische Rest ergibt sich überwiegend aus den Einflüssen von [[Abfluss]], Meteorologie sowie lokalen topografischen Besonderheiten. | |||
=Grundlagen= | =Grundlagen= | ||
Grundlegende Informationen zu dem Thema ''Harmonische Analyse'' werden in [[Harmonische Analyse des Wasserstands]] dargelegt. Die bei der Berechnung der Parameter von Tide-Ellipsen in dem Programm [[NCANALYSE]] eingesetzten Methoden werden in [[ | Grundlegende Informationen zu dem Thema ''Harmonische Analyse'' werden in [[Harmonische Analyse des Wasserstands]] dargelegt. | ||
Die bei der Berechnung der Parameter von [[Tide]]-Ellipsen in dem Programm [[NCANALYSE]] eingesetzten Methoden werden in | |||
[[Medium:zhigang_xu_2000_ellipse_parameters_conversion.pdf|zhigang_xu_2000_ellipse_parameters_conversion.pdf]] | |||
ausführlich beschrieben. | |||
Der Ablauf der Berechnung ist wie folgt (die in Klammern angegebenen Bezeichnungen beziehen sich auf die Gleichungsnummern in Zhigang Xu (2000)): | |||
# Berechne Amplitude und Phase für die u-Komponente (x- oder Ost-Komponente) der Strömung (Gl. 1). | |||
# Berechne Amplitude und Phase für die v-Komponente (y- oder Nord-Komponente) der Strömung (Gl 2). | |||
# Ermittle die äquivalenten komplexen Formulierungen für u und v (Gl. 9 und Gl. 10). | |||
# Leite daraus die entsprechenden, ebenfalls komplexen Formulierungen für die (mathematisch) positiv sowie negativ rotierenden kreisförmigen Komponenten ab (Gl. 7 und Gl. 8). Aus deren Überlagerung folgt die [[Tide]]-Ellipse. | |||
# Berechne aus den zuletzt genannten (komplexen) Formulierungen jeweils Betrag und Phase der positiv und negativ rotierenden Komponente. | |||
# Damit können große und kleine Halbachse (Gl. 14 und Gl. 19), deren Orientierung (Gl. 15 und Gl. 20) sowie die Exzentrizität (Gl. 21) berechnet werden. | |||
Die vorstehenden Berechnungen müssen für jede [[Partialtide]] durchgeführt werden. | |||
=Definition der Kennwerte= | =Definition der Kennwerte= | ||
==Residuelle Strömungsgeschwindigkeit | Für die Analyse wird die Least Square Fit Methode (LSQF) verwendet. | ||
Residuelle Strömungsgeschwindigkeit (aus Partialtidenanalyse) innerhalb des Analysezeitraums. | In dem Programm [[NCANALYSE]] wird hierfür u. a. auf Methoden des Softwarepaketes LAPACK zurückgegriffen. | ||
Alle Größen werden ausschließlich für dauerhaft überflutete Gebiete berechnet. | |||
Für jede [[Partialtide]] werden die nachfolgend aufgezählten Größen ermittelt. | |||
==Residuelle Strömungsgeschwindigkeit== | |||
Residuelle Strömungsgeschwindigkeit (aus Partialtidenanalyse) innerhalb des Analysezeitraums. | |||
Diese Größe beschreibt den nicht durch Partialtiden erklärbaren Anteil an der Strömungsgeschwindigkeit. | |||
Dieser aperiodische Rest hängt einerseits von [[Abfluss]], Meteorologie und lokalen topografischen Besonderheiten ab, ist | |||
zusätzlich aber auch von der Anzahl der in der Analyse benutzten Partialtiden abhängig. | |||
<!-- [[Beispielgrafik: Residuelle Strömungsgeschwindigkeit (aus Partialtidenanalyse)]]. <br /> --> | <!-- [[Beispielgrafik: Residuelle Strömungsgeschwindigkeit (aus Partialtidenanalyse)]]. <br /> --> | ||
==Amplitude der Strömungsgeschwindigkeit== | ==Amplitude und Phase der Partialtiden== | ||
Die nachfolgenden Größen bilden die Basis zur späteren Berechnung verschiedener Parameter von [[Tide]]-Ellipsen. | |||
Zusätzlich werden für alle vier Größen auch Standard-Fehler berechnet, die aber hier nicht mit aufgeführt werden. | |||
Diese Größen sind i.d.R. für eine Visualisierung ungeeignet. | |||
===Amplitude der u-Komponente=== | |||
Amplitude der u-Komponente (x- oder Ost-Komponente) einer [[Partialtide]]. | |||
<!-- eine Beispielgrafik für diese Größe ist nicht sinnvoll --> | |||
===Amplitude der v-Komponente=== | |||
Amplitude der v-Komponente (y- oder Nord-Komponente) einer [[Partialtide]]. | |||
<!-- eine Beispielgrafik für diese Größe ist nicht sinnvoll --> | |||
===Phase der u-Komponente=== | |||
Phase (in Grad) der u-Komponente (x- oder Ost-Komponente) einer [[Partialtide]]. | |||
<!-- eine Beispielgrafik für diese Größe ist nicht sinnvoll --> | |||
===Phase der v-Komponente=== | |||
Phase (in Grad) der v-Komponente (y- oder Nord-Komponente) einer [[Partialtide]]. | |||
<!-- eine Beispielgrafik für diese Größe ist nicht sinnvoll --> | |||
==Parameter von Tide-Ellipsen== | |||
Aus den zuvor aufgeführten vier Größen (Amplituden, Phasen) können alle nachfolgenden Parameter für [[Tide]]-Ellipsen, wie eingangs skizziert, berechnet werden. | |||
===Große Halbachse der Tide-Ellipse=== | |||
Große Halbachse einer [[Partialtide]]. Bezeichnet die maximale Strömungsgeschwindigkeit (der [[Partialtide]]) innerhalb einer Periode der [[Partialtide]]. | |||
<!-- [[Beispielgrafik: Große Halbachse der Strömungsgeschwindigkeit]]. <br /> --> | |||
===Kleine Halbachse der Tide-Ellipse=== | |||
Kleine Halbachse einer [[Partialtide]]. Bezeichnet die minimale Strömungsgeschwindigkeit (der [[Partialtide]]) innerhalb einer Periode der [[Partialtide]]. | |||
<!-- [[Beispielgrafik: Kleine Halbachse der Strömungsgeschwindigkeit]]. <br /> --> | |||
===Richtung der großen Halbachse der Tide-Ellipse=== | |||
Bezeichnet die (mathematisch positive) Richtung (-180.0 bis 180.0 Grad) der großen Halbachse einer [[Partialtide]]. | |||
Beschreibt die Richtung, in welcher die maximale Strömungsgeschwindigkeit der [[Partialtide]] auftritt. | |||
<!-- eine Beispielgrafik für diese Größe ist nicht sinnvoll --> | |||
===Richtung der kleinen Halbachse der Tide-Ellipse=== | |||
Bezeichnet die (mathematisch positive) Richtung (-180.0 bis 180.0 Grad) der kleinen Halbachse einer [[Partialtide]]. | |||
Beschreibt die Richtung, in welcher die minimale Strömungsgeschwindigkeit der [[Partialtide]] auftritt. | |||
<!-- eine Beispielgrafik für diese Größe ist nicht sinnvoll --> | |||
===Exzentrizität der Tide-Ellipse=== | |||
Exzentrizität, also das Verhältnis von kleiner Halbachse zu großer Halbachse einer [[Partialtide]]. Sonderfälle: | |||
* 0.0 : geradlinige Strömung, z. B. durch topografische Führung; | |||
* 1.0 : kreisförmige Strömung mit mathematisch positivem Umlaufsinn (entgegen dem Uhrzeiger); | |||
* -1.0 : kreisförmige Strömung mit mathematisch negativem Umlaufsinn (mit dem Uhrzeiger). | |||
<!-- [[Beispielgrafik: Exzentrizität der Strömungsgeschwindigkeit]]. <br /> --> | |||
===Phase der maximalen Strömungsgeschwindigkeit=== | |||
Phase (in Grad) der maximalen Strömungsgeschwindigkeit einer [[Partialtide]]. Der Wertebereich ist auf -180 bis 180 Grad beschränkt. | |||
<!-- [[Beispielgrafik: Phase der maximalen Strömungsgeschwindigkeit]]. <br /> --> | |||
===Phase der minimalen Strömungsgeschwindigkeit=== | |||
Phase (in Grad) der minimalen Strömungsgeschwindigkeit einer [[Partialtide]]. Der Wertebereich ist auf -180 bis 180 Grad beschränkt. | |||
<!-- [[Beispielgrafik: Phase der minimalen Strömungsgeschwindigkeit]]. <br /> --> | |||
==Flow Duration Asymmetry (Schiefe FDA)== | |||
Mit Hilfe dieses Parameters werden Unterschiede in den Strömungsgeschwindigkeiten bei Kenterung des Flutstroms sowie des Ebbestroms bewertet. Positive Werte weisen auf kleinere Strömungsgeschwindigkeiten bei Flutstromkenterung als bei Ebbestromkenterung hin. Man beachte, dass eine residuelle Strömung in den verwendeten Formulierungen unberücksichtigt bleibt. | |||
==Phase | ===Schiefe FDA=== | ||
Berechnung der Schiefe FDA (Summe der Beiträge aller Partialtiden) nach [https://www.researchgate.net/publication/256682174_The_contribution_to_tidal_asymmetry_by_different_combinations_of_tidal_constituents Song (2011), Gl. 16] aus Amplitude, Phase und Frequenz aller Partialtiden.<br /> | |||
<!-- [[Beispielgrafiken: Schiefe FDA aus allen Partialtiden]]. <br /> --> | |||
''Beispielgrafik noch nicht vorhanden''. <br /> | |||
Analyseprogramme: [[NCANALYSE]]. | |||
===Schiefe FDA aller Paare und Tripletts=== | |||
Berechnung der Schiefe aller Paare und Tripletts geeigneter Partialtiden nach [https://www.researchgate.net/publication/256682174_The_contribution_to_tidal_asymmetry_by_different_combinations_of_tidal_constituents Song (2011), Gl. 20 und Gl. 21] aus Amplitude, Phase und Frequenz der beteiligten Partialtiden.<br /> | |||
<!-- [[Beispielgrafiken: Schiefe FDA aus allen Partialtiden]]. <br /> --> | |||
''Beispielgrafik noch nicht vorhanden''. <br /> | |||
Analyseprogramme: [[NCANALYSE]]. | |||
== | ==Flow Velocity Asymmetry (Schiefe FVA)== | ||
Mit Hilfe dieses Parameters werden Unterschiede in den maximalen Strömungsgeschwindigkeiten bei [[Flutstrom]] sowie [[Ebbestrom]] bewertet. Positive Werte weisen auf größere Strömungsgeschwindigkeiten bei [[Flutstrom]] als bei [[Ebbestrom]] hin. Man beachte, dass die residuelle Strömung in den verwendeten Formulierungen berücksichtigt wird. | |||
== | ===Schiefe FVA=== | ||
Berechnung der Schiefe FVA (Summe der Beiträge aller Partialtiden sowie des Reststroms) nach [https://www.researchgate.net/publication/299342942_Tidal_asymmetry_in_a_funnel-shaped_estuary_with_mixed_semidiurnal_tides Gong (2016), Gl. 24] aus Amplitude und Phase aller Partialtiden sowie einem [[Reststrom]].<br /> | |||
<!-- [[Beispielgrafiken: Schiefe FVA aus allen Partialtiden]]. <br /> --> | |||
''Beispielgrafik noch nicht vorhanden''. <br /> | |||
Analyseprogramme: [[NCANALYSE]]. | |||
== | ===Schiefe FVA aller Paare und Tripletts=== | ||
Berechnung der Schiefe aller Paare und Tripletts geeigneter Partialtiden nach [https://www.researchgate.net/publication/299342942_Tidal_asymmetry_in_a_funnel-shaped_estuary_with_mixed_semidiurnal_tides Gong (2016), Gl. 24] aus Amplitude und Phase aller Partialtiden sowie einem [[Reststrom]].<br /> | |||
<!-- [[Beispielgrafiken: Schiefe aus allen Partialtiden]]. <br /> --> | |||
''Beispielgrafik noch nicht vorhanden''. <br /> | |||
Analyseprogramme: [[NCANALYSE]]. | |||
== | ===Schiefe FVA aus Reststrom (Residuum)=== | ||
Berechnung der Schiefe aus dem [[Reststrom]] (Residuum) nach [https://www.researchgate.net/publication/299342942_Tidal_asymmetry_in_a_funnel-shaped_estuary_with_mixed_semidiurnal_tides Gong (2016), Gl. 24] aus Amplitude und Phase aller Partialtiden sowie einem [[Reststrom]].<br /> | |||
<!-- [[Beispielgrafiken: Schiefe aus allen Partialtiden]]. <br /> --> | |||
''Beispielgrafik noch nicht vorhanden''. <br /> | |||
Analyseprogramme: [[NCANALYSE]]. | |||
---- | ---- |
Aktuelle Version vom 21. Oktober 2022, 10:17 Uhr
Motivation
Die Strömung kann für jeden Ort in einem Ästuar in einen durch die Gezeiten bedingten (periodischen) Anteil und den Rest (aperiodischer Rest) zerlegt werden. Der periodische Anteil wird dabei von den Partialtiden, also den astronomischen Tiden und den Flachwassertiden dominiert.
Für jede Partialtide ergibt sich eine Tide-Ellipse, welcher von der Spitze des Strömungsvektors einmal in jeder Periode (der Partialtide) durchlaufen wird (siehe nebenstehendes Bild). Charakteristische Parameter sind der Wert der großen sowie der kleinen Halbachse der Tide-Ellipse, deren Neigunggswinkel gegen z. B. Ost, sowie deren Exzentrizität. Je nach Exzentrizität degeneriert die Ellipse entweder zu einem Kreis (Exzentrizität +1.0 oder -1.0) oder einer Geraden (Exzentrizität 0.0). Mit diesen Parametern der Partialtiden können die tidebedingten Anteile an der lokalen Strömungsgeschwindigkeit quantifiziert werden. Der aperiodische Rest ergibt sich überwiegend aus den Einflüssen von Abfluss, Meteorologie sowie lokalen topografischen Besonderheiten.
Grundlagen
Grundlegende Informationen zu dem Thema Harmonische Analyse werden in Harmonische Analyse des Wasserstands dargelegt. Die bei der Berechnung der Parameter von Tide-Ellipsen in dem Programm NCANALYSE eingesetzten Methoden werden in zhigang_xu_2000_ellipse_parameters_conversion.pdf ausführlich beschrieben.
Der Ablauf der Berechnung ist wie folgt (die in Klammern angegebenen Bezeichnungen beziehen sich auf die Gleichungsnummern in Zhigang Xu (2000)):
- Berechne Amplitude und Phase für die u-Komponente (x- oder Ost-Komponente) der Strömung (Gl. 1).
- Berechne Amplitude und Phase für die v-Komponente (y- oder Nord-Komponente) der Strömung (Gl 2).
- Ermittle die äquivalenten komplexen Formulierungen für u und v (Gl. 9 und Gl. 10).
- Leite daraus die entsprechenden, ebenfalls komplexen Formulierungen für die (mathematisch) positiv sowie negativ rotierenden kreisförmigen Komponenten ab (Gl. 7 und Gl. 8). Aus deren Überlagerung folgt die Tide-Ellipse.
- Berechne aus den zuletzt genannten (komplexen) Formulierungen jeweils Betrag und Phase der positiv und negativ rotierenden Komponente.
- Damit können große und kleine Halbachse (Gl. 14 und Gl. 19), deren Orientierung (Gl. 15 und Gl. 20) sowie die Exzentrizität (Gl. 21) berechnet werden.
Die vorstehenden Berechnungen müssen für jede Partialtide durchgeführt werden.
Definition der Kennwerte
Für die Analyse wird die Least Square Fit Methode (LSQF) verwendet. In dem Programm NCANALYSE wird hierfür u. a. auf Methoden des Softwarepaketes LAPACK zurückgegriffen.
Alle Größen werden ausschließlich für dauerhaft überflutete Gebiete berechnet.
Für jede Partialtide werden die nachfolgend aufgezählten Größen ermittelt.
Residuelle Strömungsgeschwindigkeit
Residuelle Strömungsgeschwindigkeit (aus Partialtidenanalyse) innerhalb des Analysezeitraums. Diese Größe beschreibt den nicht durch Partialtiden erklärbaren Anteil an der Strömungsgeschwindigkeit. Dieser aperiodische Rest hängt einerseits von Abfluss, Meteorologie und lokalen topografischen Besonderheiten ab, ist zusätzlich aber auch von der Anzahl der in der Analyse benutzten Partialtiden abhängig.
Amplitude und Phase der Partialtiden
Die nachfolgenden Größen bilden die Basis zur späteren Berechnung verschiedener Parameter von Tide-Ellipsen. Zusätzlich werden für alle vier Größen auch Standard-Fehler berechnet, die aber hier nicht mit aufgeführt werden. Diese Größen sind i.d.R. für eine Visualisierung ungeeignet.
Amplitude der u-Komponente
Amplitude der u-Komponente (x- oder Ost-Komponente) einer Partialtide.
Amplitude der v-Komponente
Amplitude der v-Komponente (y- oder Nord-Komponente) einer Partialtide.
Phase der u-Komponente
Phase (in Grad) der u-Komponente (x- oder Ost-Komponente) einer Partialtide.
Phase der v-Komponente
Phase (in Grad) der v-Komponente (y- oder Nord-Komponente) einer Partialtide.
Parameter von Tide-Ellipsen
Aus den zuvor aufgeführten vier Größen (Amplituden, Phasen) können alle nachfolgenden Parameter für Tide-Ellipsen, wie eingangs skizziert, berechnet werden.
Große Halbachse der Tide-Ellipse
Große Halbachse einer Partialtide. Bezeichnet die maximale Strömungsgeschwindigkeit (der Partialtide) innerhalb einer Periode der Partialtide.
Kleine Halbachse der Tide-Ellipse
Kleine Halbachse einer Partialtide. Bezeichnet die minimale Strömungsgeschwindigkeit (der Partialtide) innerhalb einer Periode der Partialtide.
Richtung der großen Halbachse der Tide-Ellipse
Bezeichnet die (mathematisch positive) Richtung (-180.0 bis 180.0 Grad) der großen Halbachse einer Partialtide. Beschreibt die Richtung, in welcher die maximale Strömungsgeschwindigkeit der Partialtide auftritt.
Richtung der kleinen Halbachse der Tide-Ellipse
Bezeichnet die (mathematisch positive) Richtung (-180.0 bis 180.0 Grad) der kleinen Halbachse einer Partialtide. Beschreibt die Richtung, in welcher die minimale Strömungsgeschwindigkeit der Partialtide auftritt.
Exzentrizität der Tide-Ellipse
Exzentrizität, also das Verhältnis von kleiner Halbachse zu großer Halbachse einer Partialtide. Sonderfälle:
- 0.0 : geradlinige Strömung, z. B. durch topografische Führung;
- 1.0 : kreisförmige Strömung mit mathematisch positivem Umlaufsinn (entgegen dem Uhrzeiger);
- -1.0 : kreisförmige Strömung mit mathematisch negativem Umlaufsinn (mit dem Uhrzeiger).
Phase der maximalen Strömungsgeschwindigkeit
Phase (in Grad) der maximalen Strömungsgeschwindigkeit einer Partialtide. Der Wertebereich ist auf -180 bis 180 Grad beschränkt.
Phase der minimalen Strömungsgeschwindigkeit
Phase (in Grad) der minimalen Strömungsgeschwindigkeit einer Partialtide. Der Wertebereich ist auf -180 bis 180 Grad beschränkt.
Flow Duration Asymmetry (Schiefe FDA)
Mit Hilfe dieses Parameters werden Unterschiede in den Strömungsgeschwindigkeiten bei Kenterung des Flutstroms sowie des Ebbestroms bewertet. Positive Werte weisen auf kleinere Strömungsgeschwindigkeiten bei Flutstromkenterung als bei Ebbestromkenterung hin. Man beachte, dass eine residuelle Strömung in den verwendeten Formulierungen unberücksichtigt bleibt.
Schiefe FDA
Berechnung der Schiefe FDA (Summe der Beiträge aller Partialtiden) nach Song (2011), Gl. 16 aus Amplitude, Phase und Frequenz aller Partialtiden.
Beispielgrafik noch nicht vorhanden.
Analyseprogramme: NCANALYSE.
Schiefe FDA aller Paare und Tripletts
Berechnung der Schiefe aller Paare und Tripletts geeigneter Partialtiden nach Song (2011), Gl. 20 und Gl. 21 aus Amplitude, Phase und Frequenz der beteiligten Partialtiden.
Beispielgrafik noch nicht vorhanden.
Analyseprogramme: NCANALYSE.
Flow Velocity Asymmetry (Schiefe FVA)
Mit Hilfe dieses Parameters werden Unterschiede in den maximalen Strömungsgeschwindigkeiten bei Flutstrom sowie Ebbestrom bewertet. Positive Werte weisen auf größere Strömungsgeschwindigkeiten bei Flutstrom als bei Ebbestrom hin. Man beachte, dass die residuelle Strömung in den verwendeten Formulierungen berücksichtigt wird.
Schiefe FVA
Berechnung der Schiefe FVA (Summe der Beiträge aller Partialtiden sowie des Reststroms) nach Gong (2016), Gl. 24 aus Amplitude und Phase aller Partialtiden sowie einem Reststrom.
Beispielgrafik noch nicht vorhanden.
Analyseprogramme: NCANALYSE.
Schiefe FVA aller Paare und Tripletts
Berechnung der Schiefe aller Paare und Tripletts geeigneter Partialtiden nach Gong (2016), Gl. 24 aus Amplitude und Phase aller Partialtiden sowie einem Reststrom.
Beispielgrafik noch nicht vorhanden.
Analyseprogramme: NCANALYSE.
Schiefe FVA aus Reststrom (Residuum)
Berechnung der Schiefe aus dem Reststrom (Residuum) nach Gong (2016), Gl. 24 aus Amplitude und Phase aller Partialtiden sowie einem Reststrom.
Beispielgrafik noch nicht vorhanden.
Analyseprogramme: NCANALYSE.
zurück zu Analyse der Berechnungsergebnisse