Dirichlet-Randbedingung

Dieser Artikel wurde auf der Qualitätssicherungsseite des Portals Mathematik eingetragen. {{#ifeq: | ja | | Dies geschieht, um die Qualität der Artikel aus dem Themengebiet Mathematik auf ein akzeptables Niveau zu bringen.}}

Bitte hilf mit, die Mängel dieses Artikels zu beseitigen, und beteilige dich bitte an der Diskussion! ({{#ifeq:Portal:Mathematik/Qualitätssicherung|Dirichlet-Randbedingung|{{#switch:JA|JA=|NEIN=}}}}[[[:Vorlage:Neuer Abschnitt/URL]] Artikel eintragen])

Vorlage:Hinweisbaustein Die Randbedingung erster Art<ref>{{#invoke:Vorlage:Literatur|f}}</ref> ist ein Begriff aus der Mathematik und gehört zum Teilbereich der Theorie der Differentialgleichungen. Sie ist eine spezielle Randbedingung eines Randwertproblems, bei dem die Werte der gesuchten Funktion am Rand ihres Definitionsbereichs vorgegeben sind. Der Begriff findet sowohl bei gewöhnlichen als auch bei partiellen Differentialgleichungen Anwendung. Insbesondere im Bereich der partiellen Differentialgleichungen und dort in der Theorie der elliptischen Differentialgleichungen wird die Randbedingung erster Art auch Dirichlet-Randbedingung<ref name=":0">{{#invoke:Vorlage:Literatur|f}}</ref> genannt. Der Name geht auf den Mathematiker Peter Gustav Lejeune Dirichlet zurück. Genauso geht auch der Name des Dirichlet-Problems – der Prototyp eines partiellen Randwertproblems mit Dirichlet-Randbedingung – auf ihn zurück.<ref name=":0" />

Im Unterschied dazu legt die Randbedingung zweiter Art – auch Neumann-Randbedingung genannt – fest, welchen Wert die Ableitung der gesuchten Funktion am Rand annimmt. Daneben gibt es weitere Randbedingungen wie die gemischte und die schiefe Randbedingung, die beide Ansätze miteinander kombinieren.

Gewöhnliche Differentialgleichung

Definition der Randbedingung erster Art

Im Falle einer gewöhnlichen Differentialgleichung wird der Definitionsbereich der Lösung in der Regel als ein abgeschlossenes Intervall festgelegt. Damit besteht der Rand aus den beiden Intervall-Endpunkten. Aufgrund der Struktur gewöhnlicher Differentialgleichungen sind Randbedingungen erster Art nur für Gleichungen zweiter oder höherer Ordnung sinnvoll. Ein Randwertproblem für eine Differentialgleichung zweiter Ordnung hat die Form:

<math>\begin{cases} f(x,y(x),y'(x),y(x))=0, \quad x\in (a,b) \\y(a) = \alpha,\quad y(b) = \beta\end{cases}</math>

Hierbei ist <math>f</math> eine vorgeschriebene Funktion, <math>\alpha</math> und <math>\beta</math> sind vorgegebene Zahlen für die Funktionswerte der Lösung an den Intervall-Endpunkten.<ref>{{#invoke:Vorlage:Literatur|f}}</ref>

Beispiel

Gegeben sei das Randwertproblem

<math>\begin{cases} y = -y \\y(0) = 0,\quad y(\pi) = 0\end{cases}</math>

auf dem Intervall <math>[0,\pi]</math>.

Mit der Theorie der linearen gewöhnlichen Differentialgleichungen mit konstanten Koeffizienten erhält man als klassische Lösung der Differentialgleichung:

mit zwei frei wählbaren reellen Konstanten <math>C</math> und <math>D</math>. Mit den Randbedingungen können diese Konstanten fixiert werden. Dabei erhält man ein lineares Gleichungssystem in den Unbekannten <math>C</math> und <math>D</math>:

<math>C = 0</math>,

<math>-C = 0</math>.

Bemerkenswerterweise ist dieses System nicht eindeutig lösbar, aber es ist für beliebiges reelles <math>D</math> eine Lösung gegeben durch

<math>y(x) = D\sin x</math>.

Existenz und Eindeutigkeit

Der folgende Satz wird für homogene (<math>\alpha=\beta=0</math>) Daten formuliert. Dies ist jedoch keine Einschränkung, denn durch eine Transformation <math>\tilde{u}(x)=u(x)-r(x)</math> mit

<math>r(x)=\frac{(b-x)\alpha+(x-a)\beta}{b-a}</math>

kann ein inhomogenes Problem stets in ein homogenes Problem überführt werden.

Gegeben sei das Randwertproblem

<math>\begin{cases}-u(x)=f(x,u(x),u'(x)),\quad x\in (a,b) \\u(a)=u(b)=0\end{cases}</math>

Dabei sei <math>f\colon [a,b]\times \mathbb{R}^2 \to \mathbb R</math> eine stetige Funktion. Außerdem erfülle sie eine Lipschitz-Bedingung, das heißt, es gebe Zahlen <math>L,K> 0 </math>, so dass für alle <math>x\in[a,b]</math> und für alle <math>s,t,s',t'\in \mathbb{R}</math> die Ungleichung

erfüllt sei. Weiterhin gelte

Sei <math>w</math> eine Lösung von

Die Lösung <math>w</math> verschwinde für <math>x=a</math> und <math>\alpha(L,K)</math> sei die erste eindeutige Zahl, so dass <math>w'(x)=0</math> für <math>x=a+\alpha(L,K)</math>. Dann hat das zugrunde liegende Problem genau eine Lösung, falls

<math>b-a<2\alpha(L,K).</math><ref>{{#invoke:Vorlage:Literatur|f}}</ref>

Gilt hingegen <math>b-a \geq 2\alpha(L,K)</math>, so muss keine Lösung existieren oder sie muss nicht eindeutig sein. Weiterhin gilt

<math>\alpha(L,K)=\begin{cases}

\frac{2}{\sqrt{4K-L^2}}\arccos \frac{L}{2\sqrt{K}},& 4K-L^2>0 \\ \frac{2}{L^2-4K}\operatorname{arcosh} \frac{L}{2\sqrt{K}}, & 4K-L^2 <0;L,K>0 \\ \frac{2}{L}, & 4K-L^2=0,L>0 \\ +\infty, & \text{sonst} \end{cases}</math>

Ist die rechte Seite <math>f</math> der Differentialgleichung jedoch nur stetig und beschränkt, dann garantiert der Satz von Scorza Dragoni die Existenz einer Lösung.

Partielle Differentialgleichungen

Definition der Dirichlet-Randbedingung

Bei einer partiellen Differentialgleichung ist die alleinige Angabe von Dirichlet-Randbedingungen nur für elliptische Gleichungen auf einem beschränkten Gebiet <math>\Omega\subset\R^n</math> sinnvoll, da die anderen Typen auch Vorgaben der Anfangswerte benötigen.<ref>{{#invoke:Vorlage:Literatur|f}}</ref> Dabei werden Dirichlet-Randbedingungen auf dem Rand des Gebietes <math>\partial\Omega</math> vorgeschrieben. Wir definieren hier das Dirichletproblem für quasilineare partielle Differentialgleichungen

<math>u\in C^2(\Omega)\cap C^0(\overline\Omega)</math>

<math>\sum_{i,j=1}^n a_{ij}(x,u,\nabla u)\frac{\partial^2}{\partial x_i\partial x_j} u = f(x,u,\nabla u)</math>

<math>u(x) = g(x),\qquad x\in\partial\Omega\, .</math>

Hierbei stellt die Funktion <math>g \colon \partial\Omega\rightarrow\R</math> die vorgeschriebenen Funktionswerte der Lösung auf dem Rand dar. Allein die Frage nach der Lösbarkeit eines solchen Problems ist komplex und ist teilweise noch ungelöst.

Beispiel

Wir betrachten in diesem Beispiel auf dem Gebiet <math>\Omega = (0,\pi)^n = \{x=(x_1, \dots, x_n)\in\R^n\,:\, 0< x_i< \pi, \quad i=1,\dots,n\},</math> das folgende Randwertproblem:

<math>u\in C^2(\Omega)\cap C^0(\overline \Omega)</math>

<math>\Delta u(x) = -nu(x),\qquad x\in \Omega</math>

<math>u(x) = 0,\qquad \qquad \quad \; x\in\partial \Omega.</math>

Hierbei bezeichnet <math>\Delta</math> den Laplace-Operator. Zunächst stellen wir fest, dass <math>u\equiv 0</math> eine Lösung des Problems ist. Wir wollen noch weitere Lösungen finden. Wir nehmen nun <math>u(x)\neq 0</math> für <math>x\in\Omega</math> an und machen den folgenden Produktansatz

Für die Funktionen <math>v_k</math> leiten wir gewöhnliche Differentialgleichungen mit entsprechenden Dirichlet-Randbedingungen her. Es folgt

<math>\begin{align}

\Delta u &= \left (\frac{\partial^2}{\partial x_1^2}+ \dots +\frac{\partial^2}{\partial x_n^2}\right )u(x) \\ &= v_1(x_1)v_2(x_2)\cdot \dots \cdot v_n(x_n) + \dots + v_1(x_1)\cdot \dots \cdot v_{n-1}(x_{n-1})v_n(x_n) \\ &= u(x) \sum_{k=1}^n \frac{v_k(x_k)}{v_k(x_k)}. \end{align}</math>

Wenn nun die <math>v_k</math> dem Randwertproblem

genügen, dann ist die oben definierte Funktion <math>u</math> eine Lösung des Dirichlet-Randwertproblems für die partielle Differentialgleichung. Mit dem Beispiel für gewöhnliche Differentialgleichungen erhalten wir

und somit

als Lösung unseres Problems partieller Differentialgleichungen zu Dirichlet-Randbedingungen. Offen bleibt die Frage, ob es noch weitere Lösungen gibt.

Einzelnachweise