Poisson- Gleichung und Greensche Funktion

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

65px|Kein GFDL

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.

   |}}

{{#set:Urheber=Prof. Dr. E. Schöll, PhD|Inhaltstyp=Script|Kapitel=1|Abschnitt=3}} Kategorie:Elektrodynamik __SHOWFACTBOX__

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.