Parabolische partielle Differentialgleichung
Parabolische partielle Differentialgleichungen sind eine spezielle Klasse partieller Differentialgleichungen (PDG) zweiter oder höherer Ordnung, die bei der Beschreibung einer breiten Palette wissenschaftlicher Probleme zur Anwendung kommen. Es handelt sich dabei um sogenannte Evolutionsprobleme, in denen eine „Zeitvariable“ auftaucht und die Entwicklung in der „Zeit“ über eine Ableitung erster Ordnung beschrieben wird. Die Lösungen parabolischer Differentialgleichungen verhalten sich häufig wie die Lösungen der Wärmeleitungsgleichung, die die Wärmeleitung in Festkörpern oder die Diffusion in Flüssigkeiten und Gasen beschreibt.
Verallgemeinert man die Wärmeleitungsgleichung, erhält man die wichtige Klasse linearer parabolischer PDG zweiter Ordnung. Diese finden außer bei der Wärmeleitung zusätzlich Anwendung zum Beispiel bei der Berechnung der Ausbreitung von Schall im Meer oder der Entwicklung von Aktienoptionen (Black-Scholes-Modell). Im Folgenden werden nur parabolische Differentialgleichungen zweiter Ordnung betrachtet.
Definition im linearen Fall
Zwei Dimensionen
Die allgemeine lineare partielle Differentialgleichung 2. Ordnung mit zwei Variablen
heißt parabolisch im Punkt , wenn die Koeffizientenfunktionen der höchsten Ableitungen im Punkt die Bedingung
erfüllen. Dies bedeutet, dass die Determinante der Koeffizienten-Matrix
im Punkt den Wert 0 annimmt. Der Ursprung der Bezeichnung parabolisch kommt von der Analogie der obigen Koeffizientenbedingung zu der allgemeinen Kegelschnittgleichung
- .
Gilt bei dieser Gleichung , dann repräsentiert die Gleichung eine Parabel. Analoge Einteilungen existieren für elliptische und hyperbolische Differentialgleichungen.
n Dimensionen
Eine Verallgemeinerung auf mehrere Variablen ist die lineare partielle Differentialgleichung zweiter Ordnung
- .
In Verallgemeinerung des zweidimensionalen Falls bezeichnet man die Differentialgleichung als parabolisch im Punkt , falls die Koeffizientenmatrix positiv semidefinit und singulär ist. Dies bedeutet, dass alle Eigenwerte der Koeffizientenmatrix nichtnegativ sind und ein Eigenwert verschwindet.
Zeitabhängige Schreibweise
Im letzten Abschnitt wurde die abstrakte Klassifikation als parabolische Differentialgleichung erklärt. In vielen Anwendungen hat die singuläre Richtung der Koeffizientenmatrix die Bedeutung der Zeit . Dann ist die Lösung eine Funktion , die von der Zeit und Ortsvariablen abhängt. Da die Typeinteilung nur von den Koeffizienten der höchsten Ableitungen abhängt, kann man auch einfach nichtlineare Abhängigkeiten bei den niederen Ableitungen zulassen. Mit Koeffizientenfunktionen und einer Funktion stellt die Gleichung
eine semilineare parabolische Differentialgleichung dar, wenn die Matrix der Koeffizienten überall positiv definit ist. Die obige Form der parabolischen Differentialgleichung wird Divergenzform (in Anlehnung an den Divergenz-Operator) genannt. In nicht Divergenzform wird ein parabolischer Differentialoperator durch
notiert, wobei überall positiv definit ist.[1] Die rechte der Seite der parabolischen Differentialgleichung ist ein elliptischer Differentialoperator.
Beispiele
Wichtige Beispiele und Klassen parabolischer Differentialgleichungen sind
- die Wärmeleitungsgleichung,
- Reaktions-Diffusions-Gleichungen, bei denen die Funktion nicht vom Gradienten abhängt und die
- Konvektions-Diffusions-Gleichungen wie zum Beispiel die kompressiblen Navier-Stokes-Gleichungen oder die Fokker-Planck-Gleichung.
Wärmeleitungsgleichung
Der wichtigste Vertreter der linearen parabolischen Differentialgleichung ist die Wärmeleitungsgleichung. In einer Raumdimension lautet sie
- .
Hierbei ist die Temperatur am Ort zur Zeit , die Konstante bezeichnet die Temperaturleitfähigkeit. In mehreren Dimensionen wird die Gleichung durch
notiert. Setzt man nämlich in der Divergenz-Form der parabolischen Differentialgleichung für die Einheitsmatrix ein, dann ist der Hauptteil dieser Gleichung gerade der Laplace-Operator. Dieser ist das Standardbeispiel eines elliptischen Differentialoperators.
Anfangs- und Randwerte
Meist betrachtet man parabolische Differentialgleichungen entsprechend ihrer Struktur in „Raum“- und „Zeit“-Variablen als ein kombiniertes Anfangs- und Randwertproblem. Wird die Lösung im Innern eines räumlichen Gebiets für Zeiten gesucht, so gibt man die Anfangswerte zur Zeit durch eine Funktion
vor, die Randwerte auf dem Rand des räumlichen Gebiets werden für Zeiten Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle t>0} durch eine Funktion Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle h} (oder deren erste räumliche Ableitung)
- Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle u(x,t) = h (x,t) \quad \forall x \in \partial G , t > 0}
vorgegeben. Insgesamt erhält man das Anfangs- und Randwertproblem
- Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \begin{cases} \frac{\partial u}{\partial t}(x,t) + Lu(x,t) = f(x,t) & \text{in}\ G \times ]0,T]\\ u(x,t) = h(x,t) & \text{auf}\ \partial G \times [0,T]\\ u(x,0) = g(x) & \text{auf}\ G \times \{0\} \end{cases}\,,}
wobei Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \textstyle L = \sum_{i,j=1}^n \widetilde{a}_{i,j}(t,\vec{x}) \frac{\partial^2} {\partial x_j \partial x_i} + \sum_{i=1}^n b_i(t,\vec{x}) \frac{\partial} {\partial x_i} + c(t,\vec{x})} der räumliche Teil des parabolischen Differentialoperators ist.
Harnack-Ungleichung
Sei Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \textstyle L = \sum_{i,j=1}^n \widetilde{a}_{i,j}(t,\vec{x}) \frac{\partial^2} {\partial x_j \partial x_i} + \sum_{i=1}^n b_i(t,\vec{x}) \frac{\partial} {\partial x_i} + c(t,\vec{x})} der räumliche Teil des parabolischen Differentialoperators und Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle u} eine klassische Lösung der parabolischen Differentialgleichung
- Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \frac{\partial u}{\partial t}(x,t) + Lu(x,t) = f(x,t)\ \text{in}\ G \times ]0,T]}
mit Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle u \geq 0} in Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle G \times ]0,T]} . Sei außerdem Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle V \subset \subset G} eine echte zusammenhängende Teilmenge. Dann existiert für Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle 0 < t_1 < t_2 \leq T} eine Konstante Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle C} , so dass
- Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \sup_V u(\cdot, t_1) \leq C \inf_V u(\cdot, t_2)}
gilt. Die Konstante ist abhängig von Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle V} , Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle t_1} , Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle t_2} und den Koeffizienten von Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle L} .[2]
Maximumprinzip
Sei Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \textstyle L = \sum_{i,j=1}^n \widetilde{a}_{i,j}(t,\vec{x}) \frac{\partial^2} {\partial x_j \partial x_i} + \sum_{i=1}^n b_i(t,\vec{x}) \frac{\partial} {\partial x_i} + c(t,\vec{x})} ebenfalls wieder der räumliche Teil des parabolischen Differentialoperators, die Funktion Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle c} sei nicht negativ in Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle G \times ]0,T]} also Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle u \geq 0} und Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle u} sei eine klassische Lösung der parabolischen Differentialgleichung
- Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \frac{\partial u}{\partial t}(x,t) + Lu(x,t) = f(x,t)\ \text{in}\ G \times ]0,T]}
Außerdem sei Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle G} zusammenhängend.
- Falls
- Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \frac{\partial}{\partial t} u(x,t) + Lu(x,t) \leq 0\ \text{in}\ G \times ]0,T]}
- gilt und Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle u} ein nicht negatives Maximum über Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \overline{G} \times ]0,T]} im Punkt Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle (x_0,t_0) \in G \times ]0,T]} annimmt, dann ist Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle u} konstant in Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle G \times ]0,t_0]} .
- Analog, falls
- Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \frac{\partial}{\partial t} u(x,t) + Lu(x,t) \geq 0\ \text{in}\ G \times ]0,T]}
- gilt und Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle u} ein nicht positives Maximum über Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \overline{G} \times ]0,T]} im Punkt Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle (x_0,t_0) \in G \times ]0,T]} annimmt, dann ist Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle u} konstant in Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle G \times ]0,t_0]} .[3]
Numerische Verfahren für parabolische Anfangs-Randwert-Probleme
Wenn sich das Definitionsgebiet Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle G} der Gleichung mit der Zeit nicht verändert, stellt das parabolische Anfangs-Randwert-Problem in Zeitrichtung Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle t} ein Anfangswertproblem und in Ortsrichtung ein Randwertproblem für eine elliptische Differentialgleichung dar. Bei der numerischen Behandlung kann man diese beiden Probleme i. w. getrennt angehen. Dabei gibt es zwei Ansätze:
- Linienmethode (engl. MOL=method-of-lines): man diskretisiert zuerst im Ort, indem man Standardverfahren für elliptische Randwertprobleme einsetzt, wie das Differenzenverfahren oder Finite-Elemente-Methoden. Dadurch erhält man ein gewöhnliches Anfangswertproblem sehr großer Dimension für die Freiheitsgrade der Diskretisierung. Dieses ist aber ein steifes Anfangswertproblem und sollte mit impliziten oder linear-impliziten Verfahren, wie Rosenbrock-Wanner-Verfahren oder BDF-Verfahren gelöst werden. Der Vorteil dieses Zugangs ist, dass man für die Zeitintegration die genannten Standard-Verfahren einsetzen kann. Der Nachteil ist, dass die Ortsdiskretisierung fest ist, und daher lokale, zeitabhängige Verfeinerungen nicht möglich sind.
- Rothe-Methode: Man diskretisiert zuerst in der Zeit mit einem der gerade genannten Verfahren für steife Anfangswertprobleme. Dadurch erhält in jedem Zeitschritt ein elliptisches Randwertproblem für die aktuelle Lösung Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle u(t_k,\vec{x})} auf dem Gebiet Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle G} . Zur Lösung dieses Randwertproblems können jetzt auch z. B. Finite-Elemente-Verfahren mit adaptiver Gitteranpassung eingesetzt werden. Die Programmierung ist aber viel aufwändiger als bei der Linienmethode.
Ein einfaches numerisches Verfahren für parabolische Probleme ist das Crank-Nicolson-Verfahren. Dieses verwendet einerseits für die Ortsdiskretisierung die Finite-Differenzen-Methode mit festem Gitter und als Zeitdiskretisierung die implizite Trapez-Methode.
Literatur
- Gerhard Dziuk: Theorie und Numerik partieller Differentialgleichungen. de Gruyter, Berlin 2010, ISBN 978-3-11-014843-5, S. 183–253.
Einzelnachweise
- ↑ Lawrence C. Evans: Partial Differential Equations. Reprinted with corrections. American Mathematical Society, Providence RI 2008, ISBN 978-0-8218-0772-9, S. 350.
- ↑ Lawrence C. Evans: Partial Differential Equations. Reprinted with corrections. American Mathematical Society, Providence RI 2008, ISBN 978-0-8218-0772-9, S. 370f.
- ↑ Lawrence C. Evans: Partial Differential Equations. Reprinted with corrections. American Mathematical Society, Providence RI 2008, ISBN 978-0-8218-0772-9, S. 376f.