Dirichlet-Randbedingung

aus Wikipedia, der freien Enzyklopädie
Wechseln zu: Navigation, Suche

Als Dirichlet-Randbedingung (nach Peter Gustav Lejeune Dirichlet) bezeichnet man im Zusammenhang mit Differentialgleichungen (genauer: Randwertproblemen) Werte, die auf dem jeweiligen Rand des Definitionsbereichs von der Funktion angenommen werden sollen. Weitere Randbedingungen sind beispielsweise Neumann-Randbedingungen oder schiefe Randbedingungen.

Gewöhnliche Differentialgleichung[Bearbeiten]

Das Dirichletproblem[Bearbeiten]

Im Falle einer gewöhnlichen Differentialgleichung ist der Definitionsbereich der Funktion ein abgeschlossenes Intervall. Folglich besteht der Rand des Definitionsbereiches nur aus dem rechten und dem linken Intervall-Ende. Aufgrund der Freiheit in gewöhnlichen Differentialgleichungen sind Dirichlet-Randbedingungen nur für Gleichungen von zweiter oder höherer Ordnung sinnvoll. In diesem Fall sieht ein Dirichletproblem, d.h. eine Differentialgleichung mit Dirichlet-Randbedingung, folgendermaßen aus:

\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}

Hierbei ist f eine vorgeschriebene Funktion, \alpha und \beta sind vorgeschriebene reelle Zahlen für die Funktionswerte einer Lösung an den Intervallenden. Schließlich suchen wir eine (klassische) Lösung y \in C^2(a,b)\cap C^0[a,b] aus der angegebenen Regularitätsklasse.

Beispiel[Bearbeiten]

Wir wählen als unser Intervall [0,\pi] und betrachten das folgende Dirichletproblem:

\begin{cases} y'' = -y \\y(0) = 0,\quad y(\pi) = 0\end{cases}

Mit der Theorie der linearen gewöhnlichen Differentialgleichungen mit konstanten Koeffizienten erhalten wir zunächst als allgemeine (klassische) Lösung der Differentialgleichung:

 y (x) = C\cos x + D\sin x

mit zwei frei wählbaren reellen Konstanten C und D. Wir benutzen die Randbedingungen, um diese Konstanten zu fixieren. Dabei erhalten wir ein lineares Gleichungssystem in den Unbekannten C und D:

C = 0,
-C = 0.

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

y(x) = D\sin x.

Existenz und Eindeutigkeit[Bearbeiten]

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

r(x)=\frac{(b-x)\alpha+(x-a)\beta}{b-a}

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

Gegeben sei die Aufgabe

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

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

|f(x,s,s')-f(x,t,t')|\leq L|s-t|+K|s'-t'|

erfüllt sei. Weiterhin gelte

L\frac{(b-a)^2}{8}+K\frac{b-a}{2}<1.

Sei w eine Lösung von

 w''(x)+Lw'(x)+Kw(x)=0,\quad x\in(a,b).

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

b-a<2\alpha(L,K).

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

\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}{\rm 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}

Einen Beweis dieses Satzes findet man in Bailey, Shampine, Waltman. Nonlinear two-point boundary value problems. Academic Press, 1968.

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

Partielle Differentialgleichungen[Bearbeiten]

Das Dirichletproblem[Bearbeiten]

Bei einer partiellen Differentialgleichung ist die alleinige Angabe von Neumann-Randbedingungen nur für elliptische Gleichungen auf einem beschränkten Gebiet \Omega\subset\R^n sinnvoll, da die anderen Typen auch Vorgaben der Anfangswerte benötigen. Dabei werden Dirichlet-Randbedingungen auf dem Rand des Gebietes \partial\Omega vorgeschrieben. Wir definieren hier das Dirichletproblem für eine quasilineare partielle Differentialgleichung:

u=u(x)\in C^2(\Omega)\cap C^0(\overline\Omega)
\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)
u(x) = g(x),\qquad x\in\partial\Omega

Hierbei stellt die Funktion g=g(x):\partial\Omega\rightarrow\R die vorgeschriebenen Funktionswerte unserer Lösung dar. Allein die Frage nach der Lösbarkeit eines solchen Problems ist schon sehr anspruchsvoll und steht im Mittelpunkt der aktuellen Forschung. Es ist auch sehr schwierig, eine allgemeingültige Lösungsmethode anzugeben.

Beispiel[Bearbeiten]

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

u=u(x)\in C^2(\Omega)\cap C^0(\overline \Omega)
\Delta u(x) = -nu(x),\qquad x\in \Omega
u(x) = 0,\qquad x\in\partial \Omega.

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

 u(x) = v_1(x_1)\cdot ...\cdot v_n(x_n) = \prod_{k=1}^n v_k(x_k).

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

 \Delta u = \left (\frac{\partial^2}{\partial x_1^2}+...+\frac{\partial^2}{\partial x_n^2}\right )u(x)
   = v_1''(x_1)v_2(x_2)\cdot...\cdot v_n(x_n) + ... + v_1(x_1)\cdot ...\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)}.

Wenn nun die v_k dem Randwertproblem

v_k \in C^2(0,\pi)\cap C^0[0,\pi]
v_k'' \,=\, -v_k
v_k(0) = 0,\quad v_k(\pi) = 0

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

v_k(x) \,=\, D_k \sin(x_k)

und somit

u(x) = D\prod_{k=1}^n\sin(x_k)

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

Literatur[Bearbeiten]

  • D. Gilbarg, N.S. Trudinger: Partial Differential Equations of Second Order, Springer-Verlag, Berlin 1998, ISBN 3-540-41160-7.