Greens-Funktion für die Helmholtz-Gleichung

Ich habe Probleme beim Ableiten der Greens-Funktion für die Helmholtz-Gleichung. Ich weiß zufällig, was die Antwort ist, aber ich habe Mühe, sie mit typischen Tools zur Berechnung von Greens-Funktionen tatsächlich zu berechnen.

Insbesondere löse ich diese Gleichung:

( X 2 + k 2 ) G ( X , X ' ) = δ ( X X ' ) X R 3

Ich weiß, dass die Lösung ist

G ( X , X ' ) = e ich k | X X ' | 4 π | X X ' |

aber ich würde gerne verstehen, wie man diese Antwort erhält (damit ich Greens-Funktionen für andere, kompliziertere Gleichungen finden kann). Der typische Ansatz scheint zu sein, im Fourier-Raum zu lösen und dann wieder in den realen Raum zu transformieren, aber ich habe große Probleme mit der Transformation zurück in den realen Raum. Insbesondere die Auswertung dieses Integrals ist schwierig:

G ( X , X ' ) = D 3 P e ich P ( X X ' ) P 2 + k 2

(Wo P ist meine Fourier-Raum-Variable, und k ist dasselbe k wie in der ursprünglichen Gleichung oben). Gibt es einen einfacheren Weg? Mache ich etwas falsch?

Haben Sie schon vom Cauchy-Residuensatz und dem Jordan-Lemma gehört? siehe en.wikipedia.org/wiki/Jordan%27s_lemma

Antworten (1)

Das Integral ist nicht schwer! Die Maßnahme D 3 P ist gleich | P | 2 D | P | D ( cos θ ) D ϕ und das Exponential exp { ich P ( X X ' ) } Ist exp { ich | P | | X X ' | cos θ } , also ist alles ein einfaches Integral über entkoppelte Skalare.

Das sollen wir zeigen

G ( | X | ) = R 3 D 3 P ( 2 π ) 3 e ich P X | P | 2 + M 2 = 1 4 π | X | e M | X |
(Beachten Sie, dass die OP-Gleichung nicht ganz korrekt ist: Es gibt keine ich im Exponenten auf der rechten Seite, wenn sein k 2 (Mein M 2 ) ist positiv). Wir verwenden das Maß, das ich oben angegeben habe, um die Integrale über die Winkel zu berechnen θ Und ϕ um zu sehen, dass dies gleich ist
G ( | X | ) = C 1 | X | 0 D | P | | P | Sünde | P | | X | | P | 2 + M 2
für einige konstant C . Wir müssen das Lemma von Jordan nicht verwenden, um das verbleibende Integral zu machen. Wir beobachten, dass das elementare Integral
0 e M X Sünde ( P X ) D X = P P 2 + M 2
ist eine Halbzeilen-Fourier-Transformation und das Gewünschte | P | Integral ist seine inverse Fourier-Transformation.

Wenn die OP ist k 2 negativ wird, haben wir eine Wellengleichung. Dies hat keine eindeutige Green-Funktion, und in diesem Fall gibt es Singularitäten auf der Integrationskontur des Finales | P | Integral. Wie wir die Kontur um sie herum routen, indem wir hinzufügen ± ich ϵ Na dann bestimmt ob wir tauschen M von ich | k | oder ich | k | und physikalisch bezieht sich auf ausgehende Wellen (und eine kausale Green-Funktion) oder eingehende Wellen (und eine antikausale Green-Funktion). In diesem Fall sind einige Kenntnisse über komplexe Variablentechniken nützlich, aber immer noch nicht notwendig, da wir immer noch eine Fourier-Transformation invertieren können.

Das Integral ist singulär, mit zwei einfachen Polen auf der reellen Achse. Der gebräuchlichste Lösungsweg ist die Verwendung des Residue-Theorems, aber selbst in diesem Fall ist es nicht trivial
@Alessandro Zunino: Ich habe die erste Antwort bearbeitet, um dieses Problem zu beheben.
Danke schön! Eine Folgefrage: Welche Parameter dieser allgemeinen Lösung kann ich anpassen, um die Greens-Funktion zu zwingen, bestimmte Randbedingungen zu erfüllen? Insbesondere lebt die Funktion von My Green in einem gebundenen sphärischen Bereich | X X ' | < R , und an der Domainwand habe ich D D X G ( | X X | ) G ( | X X ' | | X X ' |
@alexvas Ich würde eine Methode zur Bildlösung ausprobieren, aber ich bin mir nicht sicher, ob es einfach wäre, da die übliche grüne Coulomb-Funktion Eigenschaften hat, die nicht gelten, wenn k ungleich Null ist. Die andere Möglichkeit ist eine Reihenentwicklung in den eigentlichen Eigenfunktionen für den Definitionsbereich. Dabei handelt es sich um sphärische Bessel-Funktionen
Was meinen Sie, wenn Sie "Eigenfunktionen des Bereichs" sagen?
@alexvas Sie sagen, dass Sie innerhalb einer Domäne lösen möchten, die eine Kugel ist. Finden Sie die Eigenfunktionen φ N ( R ) und erweitern G ( R , R ' ) = N φ N ( R ) φ N ( R ' ) / λ N
@AlessandroZunino Ich weiß, dass einfache Pole auf der realen Achse erfordern, dass die Summe der eingeschlossenen Residuen gemäß Jordans Lemma mit pi multipliziert wird, aber ich folge der Logik nicht vollständig. Wie würde man das Integral über Jordans Lemma berechnen?