MATHEMATISCHESINSTITUT
PROF. DR. CHRISTIANEHELZEL
DAVIDKERKMANN
23. NOVEMBER2017
16 17 18 19 Σ
NAME: MAT-NR.:
NAME: MAT-NR.:
Numerische Verfahren f¨ur hyperbolische Erhaltungsgleichungen – 6. ¨Ubungsblatt
Aufgabe 16: Beweisen Sie Satz 1.23 der Vorlesung.
Satz 1.23: Sei q eine schwache L¨osung von (1) und S eine glatte Kurve in R×R+, entlang derer q unstetig ist.
Seien f, η, F ∈ C2(R), f00 > 0, η00 > 0 und F erf¨ulle F0 =η0f0. Sei weiter q eine Entropiel¨osung im Sinne von Definition 1.22.
Dann erf¨ulltq an allen Unstetigkeitsstellen die Lax-Entropiebedingung.
Siehe Kr¨oner, Numerical Schemes for Conservation Laws.
Aufgabe 17: Schreiben Sie ein Programm zur Approximation der Burgersgleichung mit periodi- schen Anfangsdaten mittels
• Lax-Friedrichs Verfahren Qn+1i = 1
2 Qni−1+Qni+1
− ∆t
2∆x f(Qni+1)−f(Qni−1)
• Richtmyer 2-Schritt Lax-Wendroff Verfahren.
Qn+1i =Qni − ∆t
∆x
Fn+
1 2
i+12 −Fn+
1 2
i−12
mit
Qn+
1 2
i−1
2
= 1
2 Qni−1+Qni
− ∆t
2∆x f(Qni)−f(Qni−1) Fn+
1 2
i−12 =f(Qn+
1 2
i−12 )
Testen Sie Ihr Programm f¨ur verschiedene Anfangsdaten.
b.w.
Aufgabe 18: Seien u,v klassische L¨osungen von
∂tu+∂xf(u) =ε∂xxu u(x,0) =u0(x) bzw.
∂tv+∂xf(v) =ε∂xxv v(x,0) =v0(x).
Zeigen Sie, dass
ku(·, t)−v(·, t)kL1 ≤ ku0−v0kL1
gilt.
Aufgabe 19: Betrachten Sie die Buckley-Leverett Gleichung, d.h. die Erhaltungsgleichung
∂tq+∂xf(q) = 0 mit der Flussfunktion
q2
q2+a(1−q)2, 0< a <1.
Bestimmen Sie die L¨osungsstruktur des Riemann Problems q(x,0) =
1 : x <0 0 : x >0.
Lesen Sie Abschnitt 16.1 aus dem Buch von LeVeque zur Bedeutung der Buckley-Leverett Gleichung.
Abgabe am 30. November 2017 am Beginn der Vorlesung.
Abgabe der Programmieraufgaben bis zum 30. November 2017 um 14:00 an david.kerkmann@hhu.de.
Besprechung in der ¨Ubung am 8. Dezember 2017.