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

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

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