Neumann-Randbedingung

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

Eine Neumann-Randbedingung (nach Carl Gottfried Neumann) bezeichnet im Zusammenhang mit Differentialgleichungen (genauer: Randwertproblemen) Werte, die auf dem Rand des Definitionsbereichs für die Normalableitung der Lösung vorgegeben werden. Weitere Randbedingungen sind beispielsweise Dirichlet-Randbedingungen oder schiefe Randbedingungen.

Gewöhnliche Differentialgleichung[Bearbeiten]

Das Neumannproblem[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 Intervallende. Aufgrund der Freiheit in gewöhnlichen Differentialgleichungen sind Neumann-Randbedingungen nur für Gleichungen von zweiter oder höherer Ordnung sinnvoll. In diesem Fall sieht ein Neumannproblem, d.h. eine Differentialgleichung mit Neumann-Randbedingung folgendermaßen aus:

y \in C^2(a,b)\cap C^1[a,b]
y^{\prime\prime} = f(x,y,y')
y^\prime(a) = \alpha,\quad y'(b) = \beta

Hierbei ist die rechte Seite f der Differentialgleichung eine vorgeschriebene Funktion, \alpha und \beta sind vorgeschriebene reelle Zahlen für die Werte der ersten Ableitung einer Lösung an den Intervallenden. Schließlich wird eine Lösung y aus der angegebenen Regularitätsklasse gesucht.

Beispiel für eine gewöhnliche Differentialgleichung[Bearbeiten]

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

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

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

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

mit der Ableitung

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

und 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:

D = 0,
-D = 0.

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

y(x) = C\cos x.

Partielle Differentialgleichungen[Bearbeiten]

Das Neumannproblem[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 Neumann-Randbedingungen auf dem Rand des Gebietes \partial\Omega vorgeschrieben. Es wird also die Ableitung der Lösung in Richtung der äußeren Normalen vorgeschrieben. Damit die Ableitung in Richtung der äußeren Normalen an das Gebiet sinnvoll ist, muss dabei notwendig vorausgesetzt werden, dass es sich um einen C^1-Rand handelt.

Wir definieren hier das Neumannproblem für eine quasilineare partielle Differentialgleichung:

u=u(x)\in C^2(\Omega)\cap C^1(\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)
\frac{\partial u(x)}{\partial\mathbf n} = g(x),\qquad x\in\partial\Omega

Hierbei stellt die Funktion g=g(x):\partial\Omega\rightarrow\R die vorgeschriebene Ableitung in Richtung der äußeren Normalen \mathbf n(x) an \partial\Omega von unserer Lösung dar. Allein die Frage nach der Lösbarkeit eines solchen Problemes ist schon sehr anspruchsvoll und steht im Mittelpunkt der aktuellen Forschung. Es ist auch sehr schwierig, eine allgemeingültige Lösungsmethode anzugeben.

Ermittlung notwendiger Bedingungen[Bearbeiten]

Es ist jedoch zu beachten, dass allein die Gültigkeit des gaußschen Integralsatzes eine weitere (notwendige) Bedingung an die Daten und an Lösungen unseres Neumannproblems darstellt. Wir haben hierzu lediglich den gaußschen Integralsatz auf das Vektorfeld f(x) = \nabla u(x) anzuwenden.

Wenn wir beispielsweise eine Lösung eines einfachen linearen Neumannproblems mit dem Laplace-Operator \Delta betrachten:

u=u(x)\in C^2(\Omega)\cap C^1(\overline\Omega)
\Delta u(x) = f(x),\qquad x\in\Omega
\frac{\partial u(x)}{\partial\mathbf n} = g(x),\qquad x\in\partial\Omega,

so erhalten wir unter Anwendung des gaußschen Integralsatzes die Bedingung an die Daten f und g:

\int_\Omega f(x)\,d^n x = \int_\Omega \operatorname{div}\nabla u(x) \,d^nx = \oint_{\partial\Omega}\nabla u(x)\cdot \mathbf n(x)\,d^{n-1}\sigma (x) = \oint_{\partial\Omega} g(x)\,d^{n-1}\sigma(x).

Folglich ist die Gültigkeit der Gleichung

\int_\Omega f(x)\,d^n x = \oint_{\partial\Omega} g(x)\,d^{n-1}\sigma(x)

notwendig für die Lösbarkeit dieses Neumannproblems. Bei anderen Problem ist es gegebenenfalls hilfreich geeignete andere Vektorfelder zu betrachten.

Beispiel für eine partielle Differentialgleichung[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\}, mit dem regulären Rand

\delta\Omega = \{x=(x_1,...,x_n)\in\partial\Omega\,:\, für genau ein i_0\in\{1,...,n\} gilt x_{i_0}\in \{0,\pi\}\}

das folgende Randwertproblem:

u=u(x)\in C^2(\Omega)\cap C^0(\overline \Omega)
\Delta u(x) = -nu(x),\qquad x\in \Omega
\frac{\partial u(x)}{\partial\mathbf n} = 0,\qquad x\in\delta \Omega.

Hierbei bezeichnet \Delta den Laplace-Operator. Zunächst stellen wir fest, dass u\equiv 0 eine Lösung des Problems ist. Um weitere Lösungen zu finden, können wir rein formal dem Beispiel zu Dirichlet-Randbedingungen partieller Differentialgleichungen folgen, und erhalten nach einem Produktansatz:

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

Wir müssen aber beachten, dass wir hier eigentlich nicht die Nullstellenfreiheit von u fordern können, da die Cosinusfunktion bekanntermaßen eine Nullstelle bei \tfrac \pi 2 hat. Das bedeutet, dass wir nicht wissen, ob unsere formale Lösung auch wirklich Lösung unseres Neumannproblems ist. Wenn wir dies aber einsetzen, stellen wir fest, dass wir Glück haben und unser u tatsächlich Lösung unseres Problems ist.

Verallgemeinerung für partielle Differentialgleichungen[Bearbeiten]

Häufig ist es ratsam, allgemeinere Randwertprobleme wie

u=u(x)\in C^2(\Omega)\cap C^1(\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)
\frac{\partial u(x)}{\partial\mathbf \nu} = g(x),\qquad x\in\partial\Omega

zu betrachten. In diesem Fall ist \mathbf \nu eine Richtungsableitung in eine äußere Richtung. Das heißt es gilt \mathbf \nu(x)\cdot\mathbf n(x)>0 für alle x\in\partial\Omega. Wir beachten aber, dass der Richtungsvektor \mathbf \nu(x) ein Datum des Problems ist.

Literatur[Bearbeiten]

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