Technische Universität Berlin
Institut für Mathematik
Dr. M. Ehrhardt
|
Vorlesung im Sommerersemester 2006
Theorie und Numerik hyperbolischer Erhaltungsgleichungen
Gliederung
18.04.2006
Teil I: Mathematische Theorie
§ 1. Skalare Erhaltungsgleichungen
1.1 Beispiele
- Beispiel 1: Lineare Advektionsgleichung
- Beispiel 2: Advektionsgleichung mit variablen Koeffizienten
- Beispiel 3: (viskose) Burgers-Gleichung
1.2 Schockentwicklung, schwache Lösungen
- Beispiel 1: Schockentwicklung bei Burgers-Gleichung
- Definition 1: schwache Lösung
UE 21.04.2006
1.3 Das Riemann-Problem - Stoßwellen - Verdünnungswellen
- Riemann-Problem
- Rankine-Hugoniot Sprungbedingung
- Beispiel 1: Burgers-Gleichung
- Gebietsregel
- Manipulation mit Erhaltungssätzen
- Handout: Wiederholung des Erlehrnten anhand der (viskosen) Burgers-Gleichung
25.04.2006
1.4 Die Entropiebedingung
- Entropiebedingung (1.Version)
- Entropiefunktionen
- Entropiebedingung (2.Version)
1.5 Zusammenfassung
1.6 Weitere Beispiele skalarer Gleichungen
02.05.2006
§ 2. Lineare hyperbolische Systeme
- Definition 1: hyperbolisch
- Definition 2: streng hyperbolisch
2.1 Beispiele
2.2 Charakteristische Variablen
- Abhängigkeitsbereich
- p-Charakteristiken
- Einfache Wellen
- Beispiel 1: Wellengleichung
UE 05.05.2006
2.3 Linearisierung von nichtlinearen Systemen
- Beispiel 1: 1D Euler-Gleichungen der Gasdynamik, Schallwellen
- Beispiel 2: isotherme Euler-Gleichungen
09.05.2006
2.4 Das Riemann-Problem
- Riemann-Problem
- Rankine-Hugoniot Sprungbedingung für Systeme
2.5 Der Hugoniot-Locus im Phasenraum
- Definition 1: Phasenraum, Hugoniot-Locus
- Beispiel 1:
16.05.2006
§ 3. Nichtlineare hyperbolische Systeme
3.1 Hugoniot-Locus, Schocks
- Riemann-Problem
- Definition 1: Hugoniot-Locus
- Hugoniot-Kurven
- Satz 1: lokale Existenz von Hugoniot-Kurven
- Beispiel 1: isotherme Euler-Gleichungen
- Satz 2: (lokale) eindeutige Lösbarkeit des Riemann-Problems
3.2 Die Entropiebedingung von Lax
- Entropiebedingung von Lax
- echt nichtlineares Feld
- Beispiel 1: isotherme Euler-Gleichungen
- Lineare Degeneration
- Kontaktunstetigkeit
- modifizierte Entropiebedingung
23.05.2006
3.3 Verdünnungswellen, Riemann-Invarianten
- Beispiel 1: isotherme Euler-Gleichungen
- Berechnung der Verdünnungswellen (1.Methode)
- Allgemeine Lösung des Riemann-Problems
- Beispiel 1: isotherme Euler-Gleichungen
- Satz 1: (lokale) eindeutige Lösbarkeit des Entropie-erfüllenden Riemann-Problems
- Berechnung von Kontaktunstetigkeiten und Verdünnungswellen mittels Riemann-Invarianten (2.Methode)
- Definition 1: p-Riemann-Invarianten
- Beispiel 1: isotherme Euler-Gleichungen
3.4 Das Riemann-Problem für die vollen Euler-Gleichungen
- isentrope Euler-Gleichungen
- Riemann-Invarianten
- Lösung des Riemann-Problems
30.05.2006
Teil II: Numerische Methoden
§ 4. Numerische Methoden für lineare Gleichungen
- Gitter(-punkte), Näherung, hyperbolisches Schrittweitenverhältnis
- Finite Differenzenverfahren
- Beispiel 1: explizites Euler-Verfahren
- Beispiel 2: implizites Euler-Verfahren
- Beispiel 3: Lax-Friedrichs-Verfahren
- Beispiel 4: Leapfrog-Verfahren
- Beispiel 5: Lax-Wendroff-Verfahren
- allgemeine explizite Einschrittverfahren Un+1=Hk(Un)
4.1 Fehler, Konsistenz
- Fehler(-funktion)
- Definition 1: Konvergenz
- Beispiel 1: explizites Euler-Verfahren
- Definition 2: lokaler Diskretisierungsfehler, Konsistenz(-ordnung)
- Beispiel 2: Lax-Friedrichs-Verfahren
UE 02.06.2006
4.2 Stabilität
- Rekursionsformel für den Fehler
- Definition 1: Stabilität
- Satz 1: Lax-Richtmyer'scher Äquivalenzsatz
- Beispiel 1: Lax-Friedrichs-Verfahren
- Definition 2: (numerischer) Abhängigkeitsbereich
- CFL (Courant-Friedrichs-Levy)-Bedingung
- Beispiel 2: einseitige Methode
- Beispiel 3: 3-Punkt-Schema
- Beispiel 1: Lax-Friedrichs-Verfahren (Fortsetzung)
- Upwind-Methoden
- Beispiel 2: einseitige Methode
- Beispiel 4: Upwind-Methode
06.06.2006
4.3 Berechnung unstetiger Lösungen
- Beispiel 1: ut+ux=0, Methoden 1./2. Ordnung, numerische Resultate
- Modifizierte Gleichungen (für Verfahren 1.Ordnung)
- Beispiel 2: Lax-Friedrichs-Verfahren
- künstliche Diffusion
- Beispiel 3: Upwind-Methode
- Fehler für unstetige Daten
- Modifizierte Gleichungen (für Verfahren 2.Ordnung)
- Beispiel 4: Lax-Wendroff-Verfahren
- Beispiel 5: Beam-Warming-Verfahren
- Dispersionsrelation
- Phasengeschwindigkeit, Dispersion
- Gruppengeschwindigkeit
13.06.2006
§ 5. Konservative Methoden für nichtlineare Gleichungen
5.1 Der Satz von Lax-Wendroff
- Beispiel 1: Upwind-Methode für ut+u ux=0, numerische Resultate
- Definition 1: konservativ, numerischer Fluß, konsistent
- Interpretation
- Beispiel 2: (konservative) Upwind-Methode für ut+(u2/2)x=0
- Beispiel 3: Lax-Friedrichs-Verfahren für ut+f(u)x=0
- Beispiel 4: Lax-Wendroff-Verallgemeinerungen für ut+f(u)x=0
- Beispiel 5: Zweischritt-Verfahren
- Richtmyer-Zweischritt-Lax-Wendroff-Verfahren
- MacCormack-Verfahren
- Definition 2: Totalvariation
- Definition 3: Konvergenz
- Satz 1: Satz von Lax-Wendroff
- Entropiebedingung
- Beispiel 6: Riemann-Problem für ut+(u2/2)x=0, verallg. Upwind
- transonische Verdünnungswellen
- Approximation der Entropielösung
- Diskrete Entropiebedingung
- Numerischer Entropiefluß
- Satz 2: Grenzlösung ist Entropielösung
UE 16.06.2006
5.2 Das Godunov-Verfahren
- Beispiel 1: Upwind, Lax-Friedrichs für ut+aux=0
- Beispiel 2: CIR (Courant-Isaacson-Rees)-Verfahren
- Godunov-Verfahren
- Berechnung der Entropielösung
- Satz 1: Godunov-Verfahren erfüllt diskrete Entropiebedingung
- Korollar 1: Grenzlösung der Godunov-Methode ist Entropielösung
23.06.2006
5.3 Näherungsweise Lösung des Riemann-Problems
- approximative Riemann-Löser
- 2 Ansätze für approximatives Godunov-Verfahren
- Satz 1: Konservativität und Konsistenz des approximativen Godunov-Verfahrens
- Das Verfahren von Roe
- Definition 1: Roe-Matrix
- Vor- und Nachteile des Roe-Verfahrens
- Konstruktion eines "sonic entropy fix"
- skalare Probleme
- Beispiel 1:
- Beispiel 2:
- Systeme
- Beispiel 3:
- Beispiel 4:
- Konstruktion einer Roe-Matrix
- Beispiel 5: Roe-Matrix für isotherme Euler-Gleichungen
27.06.2006
§ 6. Konvergenz skalarer Verfahren (nichtlineare Stabilität)
6.1 TV-Stabilität
- Konvergenzbegriff
- Definition 1: TV-Stabilität
- Satz 1: Kriterium für TV-Stabilität
- Satz 2: Konvergenzsatz
6.2 Stabilitätskriterien
- Satz 3: TVD-Eigenschaft von Entropielösungen
- Definition 2: TVD (total variation diminishing)-Methode
- Satz 4: Monotonieerhaltungseigenschaft von Entropielösungen
- Definition 3: monotonieerhaltend
- Satz 5: TVD-Methoden sind monotonieerhaltend
- Satz 6: Monotonieerhaltende 3-Punkt-Verfahen sind von 1.Ordnung genau
04.07.2006
- Satz 7: L1-Kontraktionseigenschaft von Entropielösungen
- Definition 4: l1-kontrahierend
- Satz 8: l1-kontrahierende Methoden sind TVD
- Beispiel 1: linksseitige Upwind-Methode
- Satz 9: Monotonieeigenschaft von Entropielösungen
- Korollar 1: Minimum-Maximum-Eigenschaft von Entropielösungen
- Definition 5: monoton
- Beispiel 2: Lax-Friedrichs-Verfahren
- Satz 10: monotone Verfahren sind l1-kontrahierend
- Satz 11: monotone Methoden besitzen höchstens Konsistenzordnung 1
- Satz 12: Grenzlösung von konsistenten, monotonen Methoden ist Entropielösung
- Zusammenfassung
UE 07.07.2006
§ 7. Hochgenaue Methoden (Konstruktion hochauflösender TVD-Methoden)
7.1 Künstliche Viskosität
- Beispiel 1: Lax-Wendroff für ut+aux=0 mit künstlicher Viskosität
- Satz 1: Lineare monotonieerhaltende Methoden haben höchstens Konsistenzordnung 1
- Korollar 1: Stabilitätsbegriffe aus Kapitel 6.2 sind für lineare Methoden äquivalent
- lösungsabhängige künstliche Viskosität
7.2 Methoden mit Flußbegrenzung (flux-limiter)
- Beispiel 2: FL: Upwind-Fluß, FH:Lax-Wendroff-Fluß für ut+aux=0
- Satz 2: TVD-Methoden in Nähe von Extremalpunkten nicht 2.Ordnung
- Auflösung 2.Ordnung
- Satz 3: TVD-Kriterium für allgemeine Methode mit Flußbegrenzung
11.07.2006
- ``Superbee limiter`` von Roe
- Flußbegrenzung nach van Leer
- Verallgemeinerung für lineare Systeme
7.3 Methoden mit Anstiegsbegrenzung (slope-limiter)
- Algorithmus: Rekonstruktion, Evolution, Projektion
- Satz 4: TVD-Kriterium für allgemeine Methode mit Anstiegsbegrenzung
- Beispiel 3: Interpretation von Lax-Wendroff für ut+aux=0
- minmod-Anstieg
- Verallgemeinerung für lineare Systeme
- Verallgemeinerung für nichtlineare Systeme
UE 14.07.2006
§ 8. Semidiskrete Methoden
- Linienmethode
- Evolutionsgleichungen für Zellenmittel
- bessere Approximation des Flusses
- Beispiel 1: stückweise lineare Approximation
- Rekonstruktion mittels Stammfunktionen
- ENO (essentially nonoscillatory)-Verfahren
- Definition 1: ENO-Verfahren p-ter Ordnung
18.07.2006
§ 9. Mehrdimensionale Probleme
- Finite Volumen-Methoden
- Semidiskrete Methoden
- Splitting-Verfahren
- Beispiel 1: skalare Advektionsgleichung in 2D
- Splitting 1.Ordnung
- Strang-Splitting
- Beispiel 2: Lineare Erhaltungsgleichung in 2D
- Satz 1:
- Satz 2:
- TVD-Methoden in 2D
- Satz 3: TVD-Methoden in 2D haben höchstens Konsistenzordnung 1
ehrhardt@math.tu-berlin.de