Greensche Formeln

In der Mathematik, speziell der Vektoranalysis, sind die beiden greenschen Formeln (manchmal auch greensche Identitäten, greensche Sätze oder Theoreme) spezielle Anwendungen des gaußschen Integralsatzes. Sie sind benannt nach dem Mathematiker George Green. Anwendung finden sie unter anderem in der Elektrostatik bei der Berechnung von Potentialen. Die Formeln sind nicht zu verwechseln mit dem Satz von Green, bei dem es um ebene Integrale geht.

Im Folgenden sei $U \subset \mathbb{R}^n$ kompakt mit abschnittweise glattem Rand und $\phi$ und $\psi$ seien zwei Funktionen auf , wobei $\phi$ einfach und $\psi$ zweifach stetig differenzierbar sei. $\nabla$ ist der Nabla-Operator.

Erste greensche Identität

$\int \limits _{{U}}(\phi \nabla ^{2}\psi +\nabla \phi \cdot \nabla \psi )\,{\mathrm {d}}U=\int \limits _{{\partial U}}\phi {\frac {\partial \psi }{\partial n}}{\mathrm {d}}S$ ,

wobei ${\tfrac {\partial \psi }{\partial n}}=\nabla \psi \cdot {\vec {n}}$ die Normalenableitung von $\psi$ , also die Normalkomponente des Gradienten von $\psi$ auf dem Rand $\partial U$ bezeichnet.

Diese Identität lässt sich wie folgt beweisen:

$\;\int \limits _{\partial U}\phi {\frac {\partial \psi }{\partial n}}\mathrm {d} S=\int \limits _{\partial U}(\phi \,\nabla \psi )\cdot {\vec {n}}\,\mathrm {d} S=\int \limits _{U}\nabla \cdot (\phi \,\nabla \psi )\,\mathrm {d} U=\int \limits _{U}(\phi \nabla ^{2}\psi +\nabla \phi \cdot \nabla \psi )\,\mathrm {d} U$ ,

wobei im zweiten Schritt der gaußsche Integralsatz in der Form

$\;\int \limits _{\partial U}{\vec {F}}\cdot {\vec {n}}\;\mathrm {d} S=\int \limits _{U}\nabla \cdot {\vec {F}}\;\mathrm {d} U$

benutzt wurde.

Zweite greensche Identität

$\int \limits _{{U}}(\phi \nabla ^{2}\psi -\psi \nabla ^{2}\phi )\,{\mathrm {d}}U=\int \limits _{{\partial U}}\left(\phi {\frac {\partial \psi }{\partial n}}-\psi {\frac {\partial \phi }{\partial n}}\right)\,{\mathrm {d}}S$

Die zweite greensche Identität folgt aus der ersten greenschen Identität, wobei nun vorausgesetzt wird, dass auch $\phi$ zweimal stetig differenzierbar ist:

$\int \limits _{{U}}(\phi \nabla ^{2}\psi +\nabla \phi \cdot \nabla \psi )\,{\mathrm {d}}U=\int \limits _{{\partial U}}\phi {\frac {\partial \psi }{\partial n}}\,{\mathrm {d}}S$ ,

$\int \limits _{{U}}(\psi \nabla ^{2}\phi +\nabla \psi \cdot \nabla \phi )\,{\mathrm {d}}U=\int \limits _{{\partial U}}\psi {\frac {\partial \phi }{\partial n}}\,{\mathrm {d}}S$

Subtrahiert man nun die zweite Gleichung von der ersten Gleichung, so ergibt sich die zweite greensche Identität.

Anwendungen in der Elektrostatik (3D)

Eindeutigkeitssatz

Für ein elektrostatisches Potential $\phi \!\,$ gilt die Poissongleichung $\nabla ^{2}\phi =-4\pi \rho$ wobei $\rho \!\,$ die elektrische Ladungsdichte ist (gaußsches Einheitensystem). Wenn in einem Volumen $U\!\,$ die Ladungsdichte gegeben ist, und wenn zusätzlich auf dem Rand $\partial U$ die Werte von $\phi \!\,$ gegeben sind (Dirichlet-Randbedingung), dann gilt:

Innerhalb von $U\!\,$ ist $\phi ({\vec {r}})$ eindeutig bestimmt.

Beweis: Es seien $\phi _{1}({\vec {r}})$ und $\phi _{2}({\vec {r}})$ zwei Potentiale, die dieselben Vorgaben über Ladungsdichte und Randwerte erfüllen. Für die Differenzfunktion gilt dann

$\chi ({\vec {r}}):=\phi _{1}({\vec {r}})-\phi _{2}({\vec {r}})\qquad {\begin{array}{rl}\nabla ^{2}\chi ({\vec {r}})=0&{\vec {r}}\in U\\\chi ({\vec {r}})=0&{\vec {r}}\in \partial U\end{array}}$

Setzt man $\chi \!\,$ in der ersten greenschen Formel für $\phi \!\,$ und auch für $\psi \!\,$ ein, so folgt

$\int _{{U}}\nabla \chi \cdot \nabla \chi \,{\mathrm {d}}U=0$

Also muss der Gradient $\nabla \chi$ überall in $U\!\,$ verschwinden, somit $\chi \,$ konstant sein, und wegen seines Null-Randwerts sogar konstant gleich null sein. Also gilt $\phi _{1}({\vec {r}})=\!\,\phi _{2}({\vec {r}})$ innerhalb von $U\!\,$ .

N.B. Bei dem Beweis wird die Poissongleichung und somit die Ladungsdichte nur innerhalb von $U\!\,$ benutzt.

Abschirmung durch geschlossene Leiterfläche

$\partial U$ sei eine geschlossene Leiterfläche, so dass das elektrostatische Potential $\phi \!\,$ auf $\partial U$ einen konstanten Wert $\phi _{0}\!\,$ hat (Äquipotentialfläche). Zum Beispiel lässt sich $\phi _{0}\!\,=0$ physikalisch realisieren, indem die Leiterfläche geerdet wird. Nach dem Eindeutigkeitssatz ist der Potentialverlauf innerhalb von $U\!\,$ bereits durch die Ladungsverteilung in $U\!\,$ und durch den Randwert bestimmt. Folglich haben elektrische Ladungen im Aussenraum keinen Einfluss auf den Potentialverlauf im Innenraum.

Wenn die geschlossene Leiterfläche nicht geerdet ist, dann sind die Randwerte von $\phi \!\,$ immer noch konstant, aber mit unbekanntem Wert. Dieser Wert kann davon abhängen, welche Ladungen außerhalb von $U\!\,$ vorhanden sind. Der Beweis des Eindeutigkeitssatzes lässt sich dahingehend verallgemeinern, dass die Differenzfunktion $\chi \!\,$ noch konstant in $U\!\,$ , aber nicht mehr gleich null ist. Für die elektrische Feldstärke, die durch Ableitungen aus dem Potential gewonnen wird, spielt die Konstante keine Rolle; die elektrische Feldstärke ist also auch ohne Erdung abgeschirmt.

Symmetrie der greenschen Funktion

Die greensche Funktion mit Dirichlet-Randbedingung und mit vektoriellem Parameter ${\vec {r}}_{1}\in U$ ist definiert durch

${\begin{array}{cl}\nabla ^{2}G_{D}({\vec {r}},{\vec {r}}_{1})=\delta ({\vec {r}}-{\vec {r}}_{1})&{\vec {r}}\in U\\G_{D}({\vec {r}},{\vec {r}}_{1})=0&{\vec {r}}\in \partial U\end{array}}$

Bis auf einen Faktor $-4\pi$ entspricht das der Poissongleichung für ein Potential $\phi ({\vec {r}})$ , das von einer Punktladung am Ort $\vec{r}_1$ erzeugt wird, und das auf der geerdeten Oberfläche $\partial U$ den Randwert 0 hat. Die Existenz einer solchen Funktion ist physikalisch klar, und wegen des Eindeutigkeitssatzes ist sie eindeutig bestimmt. Obwohl die Rollen von ${\vec {r}}$ (Messpunkt) und $\vec{r}_1$ (Position der Ladung) physikalisch verschieden sind, besteht mathematisch eine Symmetrie:

$G_{D}({\vec {r}}_{2},{\vec {r}}_{1})=G_{D}({\vec {r}}_{1},{\vec {r}}_{2})$

Beweis: Setzt man in der zweiten greenschen Formel

$\phi ({\vec {r}})=G_{D}({\vec {r}},{\vec {r}}_{1})\qquad \psi ({\vec {r}})=G_{D}({\vec {r}},{\vec {r}}_{2})$

so erhält man auf der linken Seite Integrale mit Delta-Funktionen, die $G_{D}({\vec {r}}_{2},{\vec {r}}_{1})-G_{D}({\vec {r}}_{1},{\vec {r}}_{2})$ ergeben. Auf der rechten Seite verschwinden die Integranden wegen der Randwerte von $G_{D}$ .

Potential ausgedrückt durch Ladungsdichte und Randwerte

Verwendet man in der zweiten greenschen Formel als Integrationsvariable $\vec{r}_1$ und lässt man $\phi \!\,$ das elektrostatische Potential sein, so erhält man mit $\psi ({\vec {r}}_{1})=G_{D}({\vec {r}},{\vec {r}}_{1})$ und mit Hilfe der Symmetrie von $G_{D}$ den expliziten Ausdruck

$\phi ({\vec {r}})=-4\pi \int _{U}G_{D}({\vec {r}},{\vec {r}}_{1})\rho ({\vec {r}}_{1})dU_{1}+\int _{{\partial U}}\phi ({\vec {r}}_{1}){\frac {\partial G_{D}({\vec {r}},{\vec {r}}_{1})}{\partial n_{1}}}dS_{1}$

Integralgleichung für das Potential

Unter Anwendung der oben gezeigten greenschen Formeln lassen sich Ausdrücke für das elektrostatische Potential einer Ladungsverteilung herleiten. Dabei sei $\rho ({\vec r}\,')$ die Ladungsdichte am Ort ${\vec {r}}\,'$ . Mit $\phi ({\vec r})$ werde das Potenzial am Ort ${\vec {r}}$ bezeichnet. Gesucht ist die Funktion $\phi$ .

Wir setzen für $\psi ({\vec {r}}\,')={\frac {1}{|{\vec {r}}-{\vec {r}}\,'|}}$ . Es gilt dann:

,
- wobei $\Delta =\nabla ^{2}$ der Laplace-Operator ist,
- der Strich anzeigt, dass dieser Operator auf die gestrichene Variable wirkt
- und $\delta$ die Delta-Distribution ist.
Diese Identität ist also im Sinne von distributionellen Ableitungen zu verstehen.
$\Delta '\phi ({\vec {r}}\,')=-4\pi \rho ({\vec {r}}\,')$ mit der Ladungsverteilung $\rho$ am Ort ${\vec {r}}\,'$ .

Setzen wir beides in die zweite greensche Identität ein, erhalten wir auf der linken Seite:

$\int \limits _{{V}}(\phi ({\vec {r}}\,')(-4\pi \delta ({\vec {r}}-{\vec {r}}\,'))-{\frac {1}{|{\vec {r}}-{\vec {r}}\,'|}}(-4\pi \rho ({\vec {r}}\,')))\,{\mathrm {d}}V'=-4\pi \int \limits _{V}\phi ({\vec {r}}\,')\delta ({\vec {r}}-{\vec {r}}\,')\,{\mathrm {d}}V'+4\pi \int \limits _{V}{\frac {\rho ({\vec {r}}\,')}{|{\vec {r}}-{\vec {r}}\,'|}}\,{\mathrm {d}}V'$ .

Die rechte Seite der Identität ist:

$\int \limits _{F}\left(\phi ({\vec {r}}\,'){\frac {\partial }{\partial n'}}{\frac {1}{|{\vec {r}}-{\vec {r}}\,'|}}-{\frac {1}{|{\vec {r}}-{\vec {r}}\,'|}}{\frac {\partial }{\partial n'}}\phi ({\vec {r}}\,')\right){\mathrm {d}}F'$ .

Als Identität geschrieben:

$-4\pi \int \limits _{V}\phi ({\vec {r}}\,')\delta ({\vec {r}}-{\vec {r}}\,')\,{\mathrm {d}}V'+4\pi \int \limits _{V}{\frac {\rho ({\vec {r}}\,')}{|{\vec {r}}-{\vec {r}}\,'|}}\,{\mathrm {d}}V'=\int \limits _{F}\left(\phi ({\vec {r}}\,'){\frac {\partial }{\partial n'}}{\frac {1}{|{\vec {r}}-{\vec {r}}\,'|}}-{\frac {1}{|{\vec {r}}-{\vec {r}}\,'|}}{\frac {\partial }{\partial n'}}\phi ({\vec {r}}\,')\right){\mathrm {d}}F'$ .

Innerhalb des Volumens gilt an der Stelle ${\vec {r}}$ wegen der $\delta$ -Funktion

$-4\pi \int \limits _{V}\phi ({\vec {r}}\,')\delta ({\vec {r}}-{\vec {r}}\,'){\mathrm {d}}V'=-4\pi \phi ({\vec {r}})$

Damit können wir schließlich obige Identität nach dem Potential auflösen und erhalten:

$\phi ({\vec {r}})=\int \limits _{V}{\frac {\rho ({\vec {r}}\,')}{|{\vec {r}}-{\vec {r}}\,'|}}\,{\mathrm {d}}V'-{\frac {1}{4\pi }}\int \limits _{F}\left(\phi ({\vec {r}}\,'){\frac {\partial }{\partial n'}}{\frac {1}{|{\vec {r}}-{\vec {r}}\,'|}}-{\frac {1}{|{\vec {r}}-{\vec {r}}\,'|}}{\frac {\partial }{\partial n'}}\phi ({\vec {r}}\,')\right){\mathrm {d}}F'$ .

Basierend auf einem Artikel in:

Wikipedia.de