Poisson- Gleichung und Greensche Funktion

Elektrodynamikvorlesung von Prof. Dr. E. Schöll, PhD

Der Artikel Poisson- Gleichung und Greensche Funktion basiert auf der Vorlesungsmitschrift von Franz- Josef Schmitt des 1.Kapitels (Abschnitt 3) der Elektrodynamikvorlesung von Prof. Dr. E. Schöll, PhD.

   |}}

Poisson- Gleichung und Greensche Funktion

$\overline{E}\left(\overline{r}\right)=-\mathrm{\nabla }\mathrm{\Phi }(\overline{r})$ in $\mathrm{\nabla }\cdot \overline{E}\left(\overline{r}\right)=\frac{\rho \left(\overline{r}\right)}{{\epsilon }_0}$

liefert:

$\mathrm{\Delta }\mathrm{\Phi }(\overline{r})=-\frac{\rho \left(\overline{r}\right)}{{\epsilon }_0}$

Dies ist nicht anderes als die berühmte Poisson- Gleichung

Eine partielle DGL zur Berechnung des elektrischen Potenzials für eine vorgegebene Ladungsverteilung.

Die Eindeutigkeit kommt aus den Randbedingungen:

Entweder: 1)

$\mathrm{\Phi }(\overline{r})\to 0$

hinreichend rasch für

$r\to \mathrm{\infty }$

oder 2)

$\mathrm{\Phi }(\overline{r})$

sei gegeben auf Flächen im Endlichen, Beispielsweise Leiteroberflächen

Lösung zu 1):

$\mathrm{\Phi }(\overline{r})=\frac{1}{4\pi {\epsilon }_0}\underset{R^3}{\overset{}{\int }}{d}^{3}r\stackrel{´}{}\frac{\rho \left(\overline{r}\stackrel{´}{}\right)}{|\overline{r}-\overline{r}\stackrel{´}{}|}$

für hinreichend rasch abfallendes

$\rho \left(\overline{r}\stackrel{´}{}\right)$

Einsetzen in Poisson- Gleichung:

$\mathrm{\Delta }\mathrm{\Phi }(\overline{r})=\frac{1}{4\pi {\epsilon }_0}{\mathrm{\Delta }}_r\underset{R^3}{\overset{}{\int }}{d}^{3}r\stackrel{´}{}\frac{\rho \left(\overline{r}\stackrel{´}{}\right)}{|\overline{r}-\overline{r}\stackrel{´}{}|}=\frac{1}{4\pi {\epsilon }_0}\underset{R^3}{\overset{}{\int }}{d}^{3}r\stackrel{´}{}{\mathrm{\Delta }}_r\frac{\rho \left(\overline{r}\stackrel{´}{}\right)}{|\overline{r}-\overline{r}\stackrel{´}{}|}$,

falls Integration und Differenziation vertauschbar, also über verschiedene Koordinaten ausgeführt wird.

Man definiere für ein festes

$\overline{r}\stackrel{´}{}$,

dass

$\begin{array}{rl}& \overline{s}:=\overline{r}-\overline{r}\stackrel{´}{}\\ & {\mathrm{\nabla }}_r={\mathrm{\nabla }}_s\\ \end{array}$

Also:

$\begin{array}{rl}& {\mathrm{\Delta }}_r\frac{1}{|\overline{r}-\overline{r}\stackrel{´}{}|}={\mathrm{\nabla }}_S\left({\mathrm{\nabla }}_S\frac{1}{s}\right)=-{\mathrm{\nabla }}_S\frac{1}{s^2}\frac{\overline{s}}{s}=-\frac{1}{s^3}{\mathrm{\nabla }}_S\overline{s}-\overline{s}{\mathrm{\nabla }}_S\frac{1}{s^3}\\ & {\mathrm{\nabla }}_S\overline{s}=3\\ & \Rightarrow {\mathrm{\Delta }}_r\frac{1}{|\overline{r}-\overline{r}\stackrel{´}{}|}=-\frac{1}{s^3}{\mathrm{\nabla }}_S\overline{s}-\overline{s}{\mathrm{\nabla }}_S\frac{1}{s^3}=-\frac{3}{s^3}+\frac{1}{s^3}=0\\ \end{array}$

Dies ist aber ein Widerspruch zu

$\mathrm{\Delta }\mathrm{\Phi }(\overline{r})=-\frac{\rho \left(\overline{r}\right)}{{\epsilon }_0}$

Grund ist, dass die Vertauschung von

${\mathrm{\Delta }}_r$ und $\underset{R^3}{\overset{}{\int }}{d}^{3}r\stackrel{´}{}$

sowie auch die obige Umformung nicht erlaubt ist für

$\overline{r}=\overline{r}\stackrel{´}{}$,

also s=0  (Singularität!!)

Stattdessen für beliebige V:

Nun kann man

$\underset{\partial V}{\oint }d\overline{f}\cdot {\mathrm{\nabla }}_r$ mit $\underset{R^3}{\overset{}{\int }}{d}^{3}r\stackrel{´}{}$

vertauschen. Dies ist erlaubt, falls der Integrand von

$\underset{R^3}{\overset{}{\int }}{d}^{3}r\stackrel{´}{}$

nach der Vertauschung stetig ist!:

$\begin{array}{rl}& \underset{V}{\overset{}{\int }}{d}^{3}r\mathrm{\Delta }\mathrm{\Phi }(\overline{r})=\frac{1}{4\pi {\epsilon }_0}\underset{R^3}{\overset{}{\int }}{d}^{3}r\stackrel{´}{}\rho \left(\overline{r}\stackrel{´}{}\right)\underset{\partial V}{\oint }d\overline{f}\cdot {\mathrm{\nabla }}_r\frac{1}{|\overline{r}-\overline{r}\stackrel{´}{}|}\\ & {\mathrm{\nabla }}_r\frac{1}{|\overline{r}-\overline{r}\stackrel{´}{}|}=-\frac{(\overline{r}-\overline{r}\stackrel{´}{})}{|\overline{r}-\overline{r}\stackrel{´}{}{|}^3}\\ \end{array}$

Somit:

$\begin{array}{rl}& \underset{V}{\overset{}{\int }}{d}^{3}r\mathrm{\Delta }\mathrm{\Phi }(\overline{r})=\frac{1}{4\pi {\epsilon }_0}\underset{R^3}{\overset{}{\int }}{d}^{3}r\stackrel{´}{}\rho \left(\overline{r}\stackrel{´}{}\right)\underset{\partial V}{\oint }d\overline{f}\cdot {\mathrm{\nabla }}_r\frac{1}{|\overline{r}-\overline{r}\stackrel{´}{}|}=-\frac{1}{4\pi {\epsilon }_0}\underset{R^3}{\overset{}{\int }}{d}^{3}r\stackrel{´}{}\rho \left(\overline{r}\stackrel{´}{}\right)\underset{\partial V}{\oint }d\overline{f}\frac{\overline{r}-\overline{r}\stackrel{´}{}}{|\overline{r}-\overline{r}\stackrel{´}{}{|}^3}\\ & \underset{\partial V}{\oint }d\overline{f}\frac{\overline{r}-\overline{r}\stackrel{´}{}}{|\overline{r}-\overline{r}\stackrel{´}{}{|}^3}=\underset{}{\overset{}{\int }}d\mathrm{\Omega }\\ \end{array}$

aber:

$\underset{\partial V}{\oint }d\overline{f}\frac{\overline{r}-\overline{r}\stackrel{´}{}}{|\overline{r}-\overline{r}\stackrel{´}{}{|}^3}=\underset{}{\overset{}{\int }}d\mathrm{\Omega }=4\pi$,

falls

$\overline{r}\stackrel{´}{}\in V$

$\underset{\partial V}{\oint }d\overline{f}\frac{\overline{r}-\overline{r}\stackrel{´}{}}{|\overline{r}-\overline{r}\stackrel{´}{}{|}^3}=\underset{}{\overset{}{\int }}d\mathrm{\Omega }=0$ falls $\overline{r}\stackrel{´}{}\notin V$

Somit:

$\underset{V}{\overset{}{\int }}{d}^{3}r\mathrm{\Delta }\mathrm{\Phi }(\overline{r})=\frac{1}{4\pi {\epsilon }_0}\underset{R^3}{\overset{}{\int }}{d}^{3}r\stackrel{´}{}\rho \left(\overline{r}\stackrel{´}{}\right)\underset{\partial V}{\oint }d\overline{f}\cdot {\mathrm{\nabla }}_r\frac{1}{|\overline{r}-\overline{r}\stackrel{´}{}|}=-\frac{1}{{\epsilon }_0}\underset{V}{\overset{}{\int }}{d}^{3}r\stackrel{´}{}\rho \left(\overline{r}\stackrel{´}{}\right)$

Mathematisch streng gilt im Distributionen- Sinn:

${\mathrm{\Delta }}_r\frac{1}{|\overline{r}-\overline{r}\stackrel{´}{}|}=-4\pi \delta (\overline{r}-\overline{r}\stackrel{´}{})$

Mit Hilfe der delta- Distribution, die auf eine Testfunktion anzuwenden ist!

Greensche Funktion der Poisson- Gleichung

$\mathrm{\Delta }\mathrm{\Phi }(\overline{r})=-\frac{\rho \left(\overline{r}\stackrel{´}{}\right)}{{\epsilon }_0}\Rightarrow$

Invertierung

$\Rightarrow \mathrm{\Phi }(\overline{r})=\hat{G}\rho \left(\overline{r}\stackrel{´}{}\right)$

Mit dem Greenschen Operator

$\hat{G}$

Eine Fourier- Transformation von

$\mathrm{\Delta }\mathrm{\Phi }(\overline{r})=-\frac{\rho \left(\overline{r}\stackrel{´}{}\right)}{{\epsilon }_0}\Rightarrow$ liefert $-k^2\tilde{\mathrm{\Phi }}=-\frac{\tilde{\rho }}{{\epsilon }_0}\Rightarrow$

Man kann schreiben:

$\begin{array}{rl}& \tilde{\mathrm{\Phi }}=\tilde{\hat{G}}\tilde{\rho }\\ & \tilde{\hat{G}}:=\frac{1}{{\epsilon }_0{k}^2}\\ \end{array}$

Die einfache Fourier- Transformierte Form von

$\Rightarrow \mathrm{\Phi }(\overline{r})=\hat{G}\rho \left(\overline{r}\stackrel{´}{}\right)$,

nur dass der Fourier- transformierte Greens- Operator angegeben werden kann.

Die Rücktransformation löst dann die Poiisson-gleichung:

$\Rightarrow \mathrm{\Phi }(\overline{r})=\underset{}{\overset{}{\int }}{d}^{3}r\stackrel{´}{}\hat{G}(\overline{r}-\overline{r}\stackrel{´}{})\rho \left(\overline{r}\stackrel{´}{}\right)$

Es gilt:

${\mathrm{\Delta }}_r\hat{G}(\overline{r}-\overline{r}\stackrel{´}{})=-\frac{1}{{\epsilon }_0}\delta (\overline{r}-\overline{r}\stackrel{´}{})$

Das heißt, die Greensfunktion ist eine Lösung der Poissongleichung für eine Punktladung q=1 an

$\overline{r}\stackrel{´}{}$

Insbesondere bei speziellen Randbedingungen

$\begin{array}{c}\lim \\ \overline{r}\to \mathrm{\infty }\\ \end{array}\mathrm{\Phi }(\overline{r})=0$

ist die Greensfunktion dann:

$G(\overline{r}-\overline{r}\stackrel{´}{})=\frac{1}{4\pi {\epsilon }_0}\frac{1}{|\overline{r}-\overline{r}\stackrel{´}{}|}$

Denn

${\mathrm{\Delta }}_{r}G=\mathrm{\Delta }\frac{1}{4\pi {\epsilon }_0}\frac{1}{|\overline{r}-\overline{r}\stackrel{´}{}|}=-\frac{1}{{\epsilon }_0}\delta (\overline{r}-\overline{r}\stackrel{´}{})$

Für eine beliebige Ladungsverteilung

$\rho$

ist also die Lösung der Poissongleichung

$\mathrm{\Phi }(\overline{r})=\underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}\frac{1}{4\pi {\epsilon }_0}\frac{\rho (\overline{r}\stackrel{´}{})}{|\overline{r}-\overline{r}\stackrel{´}{}|}{d}^{3}r\stackrel{´}{}=\underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}G(\overline{r}-\overline{r}\stackrel{´}{})\rho (\overline{r}\stackrel{´}{})d^{3}r\stackrel{´}{}$

wobei die Identität insgesamt nur für die Randbedingungen

$\begin{array}{c}\lim \\ \overline{r}\to \mathrm{\infty }\\ \end{array}\mathrm{\Phi }(\overline{r})=0$

gilt, ansonsten ist G durch die andern Randbdingungen festgelegt.