Helmholtz-Gleichung

Die Helmholtz-Gleichung (nach Hermann von Helmholtz) ist eine partielle Differentialgleichung. Sie lautet:

$\mathrm{\Delta }\phi =\lambda \cdot \phi$

in einem Gebiet $\mathrm{\Omega }$ mit geeigneten Randbedingungen auf dem Rand $\partial \mathrm{\Omega }$. Dabei ist

$\mathrm{\Delta }=\sum _{k=1}^n\frac{{\partial }^2}{\partial x_k^2}.$

der Laplace-Operator in kartesischen Koordinaten.

Die Helmholtz-Gleichung ist dementsprechend eine partielle Differentialgleichung (PDGL) zweiter Ordnung aus der Klasse der elliptischen PDGL. Sie ergibt sich auch z. B. aus der Wellengleichung nach Trennung der Variablen und Annahme harmonischer Zeitabhängigkeit.

Setzt man $\lambda =0$, so erhält man die Laplace-Gleichung.

Beispiel: Partikuläre Lösung der inhomogenen Maxwellgleichungen

Eine Anwendung aus der Physik ist z. B. die Lösung der inhomogenen Maxwellgleichungen (Maxwellgleichungen mit Strömen und Ladungen). Aus diesen folgen in Gaußschen Einheiten mit der Lorenz-Eichung

$\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{A}+\frac{1}{c^2}\frac{\partial \mathrm{\Phi }}{\partial t}=0$

die inhomogenen Wellengleichungen für das elektrische Skalarpotential

$\mathrm{\Phi }$ sowie für das magnetische Vektorpotential $\overrightarrow{A}$:

$\mathrm{\Delta }\mathrm{\Phi }(\overrightarrow{r},t)-\frac{1}{c^2}\frac{{\partial }^2\mathrm{\Phi }(\overrightarrow{r},t)}{\partial t^2}=-4\pi \varrho (\overrightarrow{r},t)$

$\mathrm{\Delta }{A}_i(\overrightarrow{r},t)-\frac{1}{c^2}\frac{{\partial }^2{A}_i(\overrightarrow{r},t)}{\partial t^2}=-\frac{4\pi }{c}{j}_i(\overrightarrow{r},t)$

(hier für die einzelnen Komponenten mit: $\overrightarrow{A}=\sum _{i=1}^3{A}_i{\hat{e}}_i$)

Exemplarisch wird nun die Lösung für $\mathrm{\Phi }$ durchgeführt, die Herleitung für $\overrightarrow{A}$ geht analog.

Die allgemeine Lösung dieser Differentialgleichungen ist die Linearkombination der allgemeinen Lösung der dazugehörigen homogenen DGL sowie einer partikulären Lösung der inhomogenen DGL:

$\mathrm{\Phi }={\mathrm{\Phi }}_{\mathrm{h}\mathrm{o}\mathrm{m}\mathrm{.}}+{\mathrm{\Phi }}_{\mathrm{p}\mathrm{a}\mathrm{r}\mathrm{t}\mathrm{.}}$

Die Lösung der homogenen DGL sind die elektromagnetischen Wellen; wir beschränken uns hier auf die Herleitung einer partikulären Lösung.

Um die Wellengleichung auf die Helmholtz-Gleichung zurückzuführen betrachten wir die Fourier-Transformation von $\mathrm{\Phi }$ und $\varrho$ bezüglich $t$:

$\mathrm{\Phi }(\overrightarrow{r},t)=\frac{1}{\sqrt{2\pi }}\int d\omega {\mathrm{\Phi }}_{\omega }(\overrightarrow{r})e^{-i\omega t}$

$\varrho (\overrightarrow{r},t)=\frac{1}{\sqrt{2\pi }}\int d\omega {\varrho }_{\omega }(\overrightarrow{r})e^{-i\omega t}$

Einsetzen in die Wellengleichung liefert:

$\mathrm{\Delta }\int d\omega {\mathrm{\Phi }}_{\omega }(\overrightarrow{r})e^{-i\omega t}-\frac{1}{c^2}\frac{{\partial }^2}{\partial t^2}\int d\omega {\mathrm{\Phi }}_{\omega }(\overrightarrow{r})e^{-i\omega t}=-4\pi \int d\omega {\varrho }_{\omega }(\overrightarrow{r})e^{-i\omega t}$

$\Rightarrow \int d\omega (\mathrm{\Delta }-\frac{1}{c^2}\frac{{\partial }^2}{{\partial }^{2}t}){\mathrm{\Phi }}_{\omega }(\overrightarrow{r})e^{-i\omega t}=-4\pi \int d\omega {\varrho }_{\omega }(\overrightarrow{r})e^{-i\omega t}$

$\Rightarrow \int d\omega (\mathrm{\Delta }+\frac{{\omega }^2}{c^2}){\mathrm{\Phi }}_{\omega }(\overrightarrow{r})e^{-i\omega t}=-4\pi \int d\omega {\varrho }_{\omega }(\overrightarrow{r})e^{-i\omega t}$

Beide Integranden müssen gleich sein, da sich die $d\omega$-Integration über die gleichen Bereiche erstreckt:

$(\mathrm{\Delta }+\frac{{\omega }^2}{c^2}){\mathrm{\Phi }}_{\omega }(\overrightarrow{r})=-4\pi {\varrho }_{\omega }(\overrightarrow{r})$

Für die homogene Wellengleichung $(\varrho (\overrightarrow{r},t)=0)$ erkennen wir mit $(\mathrm{\Delta }+\frac{{\omega }^2}{c^2}){\mathrm{\Phi }}_{\omega }(\overrightarrow{r})=0$ die Helmholtz-Gleichung wieder.

Zur Lösung der inhomogenen Gleichung $(\varrho (\overrightarrow{r},t)\ne 0)$ kann eine Greensche Funktion $G(\overrightarrow{r},{\overrightarrow{r}}^{\prime })$ verwendet werden, welche die Gleichung

$(\mathrm{\Delta }+\frac{{\omega }^2}{c^2})G(\overrightarrow{r},{\overrightarrow{r}}^{\prime })=-4\pi \delta (\overrightarrow{r}-{\overrightarrow{r}}^{\prime })$

erfüllt.

Diese lautet:

$G(\overrightarrow{r},{\overrightarrow{r}}^{\prime })=\frac{\exp (\pm i\omega |\overrightarrow{r}-{\overrightarrow{r}}^{\prime }|/c)}{|\overrightarrow{r}-{\overrightarrow{r}}^{\prime }|}$

Physikalisch beschreibt diese Funktion eine Kugelwelle.

Damit erhalten wir für die gesamte Ladungsverteilung:

${\mathrm{\Phi }}_{\omega }(\overrightarrow{r})=\int d^3{r}^{\prime }{\varrho }_{\omega }({\overrightarrow{r}}^{\prime })G(\overrightarrow{r},{\overrightarrow{r}}^{\prime })=\int d^3{r}^{\prime }{\varrho }_{\omega }({\overrightarrow{r}}^{\prime })\frac{\exp (\pm i\omega |\overrightarrow{r}-{\overrightarrow{r}}^{\prime }|/c)}{|\overrightarrow{r}-{\overrightarrow{r}}^{\prime }|}$

Dieses Ergebnis setzen wir in die Fourierdarstellung von $\mathrm{\Phi }(\overrightarrow{r},t)$ ein und erhalten

$\mathrm{\Phi }(\overrightarrow{r},t)=\frac{1}{\sqrt{2\pi }}\int d\omega \int d^3{r}^{\prime }{\varrho }_{\omega }({\overrightarrow{r}}^{\prime })\frac{\exp (\pm i\omega |\overrightarrow{r}-{\overrightarrow{r}}^{\prime }|/c)}{|\overrightarrow{r}-{\overrightarrow{r}}^{\prime }|}{e}^{-i\omega t}=\frac{1}{\sqrt{2\pi }}\int d\omega \int d^3{r}^{\prime }\frac{{\varrho }_{\omega }({\overrightarrow{r}}^{\prime })}{|\overrightarrow{r}-{\overrightarrow{r}}^{\prime }|}\exp (-i\omega (\mp |\overrightarrow{r}-{\overrightarrow{r}}^{\prime }|/c+t))$

Mit $t^{\prime }:=t\mp |\overrightarrow{r}-{\overrightarrow{r}}^{\prime }|/c$ folgt:

$\mathrm{\Phi }(\overrightarrow{r},t)=\frac{1}{\sqrt{2\pi }}\int d\omega \int d^3{r}^{\prime }\frac{{\varrho }_{\omega }({\overrightarrow{r}}^{\prime })}{|\overrightarrow{r}-{\overrightarrow{r}}^{\prime }|}\exp (-i\omega t^{\prime })=\frac{1}{\sqrt{2\pi }}\int d^3{r}^{\prime }\frac{1}{|\overrightarrow{r}-{\overrightarrow{r}}^{\prime }|}\int d\omega {\varrho }_{\omega }({\overrightarrow{r}}^{\prime })e^{-i\omega t^{\prime }}$

$\Rightarrow \mathrm{\Phi }(\overrightarrow{r},t)=\int d^3{r}^{\prime }\frac{\varrho ({\overrightarrow{r}}^{\prime },t^{\prime })}{|\overrightarrow{r}-{\overrightarrow{r}}^{\prime }|}$

Dies ist die gesuchte partikuläre Lösung der inhomogenen Gleichung. Für $A_i$ folgt analog:

$A_i(\overrightarrow{r},t)=\frac{1}{c}\int d^3{r}^{\prime }\frac{j_i({\overrightarrow{r}}^{\prime },t^{\prime })}{|\overrightarrow{r}-{\overrightarrow{r}}^{\prime }|}$

$\Rightarrow \overrightarrow{A}(\overrightarrow{r},t)=\frac{1}{c}\int d^3{r}^{\prime }\frac{\overrightarrow{j}({\overrightarrow{r}}^{\prime },t^{\prime })}{|\overrightarrow{r}-{\overrightarrow{r}}^{\prime }|}$

Die physikalische Bedeutung ist, dass das zur Zeit t am Ort $\overrightarrow{r}$ beobachtete Potential von Ladungen bzw. Strömen zur Zeit t' am Ort ${\overrightarrow{r}}^{\prime }$ verursacht wurde.

Diskussion: Retardierte und avancierte Lösung

Noch steht das Vorzeichen im Argument $t\pm |\overrightarrow{r}-{\overrightarrow{r}}^{\prime }|/c$ nicht fest. Physikalisch scheint aber plausibel, dass die zeitliche Änderung einer Ladungsverteilung bei ${\overrightarrow{r}}^{\prime }$ erst zu einem späteren Zeitpunkt bei $\overrightarrow{r}$ beobachtet werden kann, da sich elektromagnetische Wellen mit der (konstanten) Lichtgeschwindigkeit c ausbreiten. Daher wählen wir das Minuszeichen als physikalisch praktikable Lösung:

$\mathrm{\Phi }(\overrightarrow{r},t{)}_{\mathrm{r}\mathrm{e}\mathrm{t}\mathrm{.}}=\int d^3{r}^{\prime }\frac{\varrho ({\overrightarrow{r}}^{\prime },t-|\overrightarrow{r}-{\overrightarrow{r}}^{\prime }|/c)}{|\overrightarrow{r}-{\overrightarrow{r}}^{\prime }|}$

Man nennt das Potential bei Wahl des Minuszeichens auch retardiertes Potential. Wählt man das Pluszeichen, so spricht man vom avancierten Potential.

Siehe auch

• Bessel-Strahl

Weblinks

• Helmholtzgleichung bei Wolfram MathWorld (engl.)