Neumannproblem

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

Das Neumannproblem

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 machen Neumann-Randbedingungen nur für Gleichungen von zweiter oder höherer Ordnung Sinn. 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'' = f(x,y,y')

$y'(a) = \alpha,\quad y'(b) = \beta$

Hierbei ist die rechte Seite f der Differentialgleichung eine vorgeschriebene Funktion, α und β 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

Wir wählen als unser Intervall [0,π] 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) = Ccosx + Dsinx

mit der Ableitung

y'(x) = − Csinx + Dcosx

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.

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

y(x) = Ccosx.

Partielle Differentialgleichungen

Das Neumannproblem

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 sich der Rand $\partial\Omega$ lokal als Nullstellenmenge einer einmal stetig differenzierbaren Funktion schreiben lässt. Wir sagen dann, es handelt sich um einen C¹-Rand.

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

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 Δ 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 eines Neumannproblems. Bei anderen Problem ist es gegebenenfalls hilfreich geeignete andere Vektorfelder zu betrachten.

Beispiel für eine partielle Differentialgleichung

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 Δ 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 $\frac \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

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

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