Quelltext der Seite Neumann-Randbedingung

Eine '''Neumann-Randbedingung''' (nach [[Carl Gottfried Neumann]]) bezeichnet im Zusammenhang mit [[Differentialgleichung]]en (genauer: [[Randwertproblem]]en) Werte, die auf dem [[Rand (Topologie)|Rand]] des Definitionsbereichs für die [[Richtungsableitung|Normalableitung]] der Lösung vorgegeben werden. Bei Neumann-Randwertproblemen werden nicht Funktionswerte, sondern Ableitungswerte vorgegeben. Weitere Randbedingungen sind beispielsweise [[Dirichlet-Randbedingung]]en (bei denen die Funktionswerte auf dem Rand vorgegeben sind) oder [[Randwertproblem#Schiefe Randbedingung|schiefe Randbedingungen]].

== Gewöhnliche Differentialgleichung ==
=== Das Neumannproblem ===
Im Falle einer [[gewöhnliche Differentialgleichung|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.&nbsp;h. eine Differentialgleichung mit Neumann-Randbedingung folgendermaßen aus:

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

Hierbei ist die rechte Seite <math>f</math> der Differentialgleichung eine vorgeschriebene Funktion, <math>\alpha</math> und <math>\beta</math> sind vorgeschriebene reelle Zahlen für die Werte der ersten Ableitung einer Lösung an den Intervallenden. Schließlich wird eine Lösung <math>y</math> aus der angegebenen [[Regularitätsklasse]] gesucht.<!-- Bei einem solchen Problem liefert der [[Satz von Scorza Dragoni]] ein Existenzresultat.-->

=== Beispiel für eine gewöhnliche Differentialgleichung ===
Wir wählen als unser Intervall <math>[0,\pi]</math> und betrachten das folgende Problem:

:<math>y \in C^2(0,\pi)\cap C^1[0,\pi]</math>
:<math>y'' = -y</math>
:<math>y'(0) = 0,\quad y'(\pi) = 0</math>

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

:<math> y (x) = C\cos x + D\sin x</math>

mit der Ableitung

:<math> y'(x) = -C\sin x + D\cos x</math>

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

:<math>D = 0,</math>
:<math>-D = 0.</math>

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

:<math>y(x) = C\cos x.</math>

== Partielle Differentialgleichungen ==
=== Das Neumannproblem ===
Bei einer [[Partielle Differentialgleichung|partiellen Differentialgleichung]] ist die alleinige Angabe von Neumann-Randbedingungen nur für elliptische Gleichungen auf einem [[Beschränktheit|beschränkten]] [[Gebiet (Mathematik)|Gebiet]] <math>\Omega\subset\R^n</math> sinnvoll, da die anderen Typen auch Vorgaben der Anfangswerte benötigen. Dabei werden Neumann-Randbedingungen auf dem Rand des Gebietes <math>\partial\Omega</math> 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 <math>C^1</math>-Rand handelt.

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

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

Hierbei stellt die Funktion <math>g\colon\partial\Omega\rightarrow\R</math> die vorgeschriebene Ableitung in Richtung der äußeren Normalen <math>\mathbf n(x)</math> an <math>\partial\Omega</math> von unserer Lösung dar. Allein die Frage nach der Lösbarkeit eines solchen Problemes ist 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ßscher Integralsatz|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 <math>f(x) = \nabla u(x)</math> anzuwenden.

Wenn wir beispielsweise eine Lösung eines einfachen linearen Neumannproblems mit dem [[Laplace-Operator]] <math>\Delta</math> betrachten:

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

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

:<math>\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).</math>

Folglich ist die Gültigkeit der Gleichung

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

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 ===
Wir betrachten in diesem Beispiel auf dem Gebiet <math>\Omega = (0,\pi)^n = \{x=(x_1,...,x_n)\in\R^n\,:\, 0< x_i< \pi, \quad i=1,\dots,n\},</math> mit dem regulären Rand

:<math>\delta\Omega = \{x=(x_1,...,x_n)\in\partial\Omega\,:\, </math> für genau ein <math>i_0\in\{1,...,n\}</math> gilt <math>x_{i_0}\in \{0,\pi\}\}</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>\frac{\partial u(x)}{\partial\mathbf n} = 0,\qquad x\in\delta \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. Um weitere Lösungen zu finden, können wir rein formal dem Beispiel zu [[Dirichlet-Randbedingung]]en partieller Differentialgleichungen folgen, und erhalten nach einem Produktansatz:

:<math>u(x) = D\prod_{k=1}^n\cos(x_k).</math>

Wir müssen aber beachten, dass wir hier eigentlich nicht die Nullstellenfreiheit von <math>u</math> fordern können, da die Cosinusfunktion bekanntermaßen eine Nullstelle bei <math>\tfrac \pi 2</math> 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 <math>u</math> tatsächlich Lösung unseres Problems ist.

== Verallgemeinerung für partielle Differentialgleichungen ==
Häufig ist es ratsam, allgemeinere Randwertprobleme wie

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

zu betrachten. In diesem Fall ist <math>\frac{\partial}{\partial \mathbf \nu}</math> eine [[Richtungsableitung]] in eine äußere Richtung. Das heißt, es gilt <math>\mathbf \nu(x)\cdot\mathbf n(x)>0</math> für alle <math>x\in\partial\Omega</math>. Wir beachten aber, dass der Richtungsvektor <math>\mathbf \nu(x)</math> 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.

[[Kategorie:Theorie der Differentialgleichungen]]