Der allgegenwärtige Planewave-Ansatz

In der Physik wird der Planewave- Ansatz (was bedeutet: eine fundierte Lösungsschätzung) sehr allgegenwärtig verwendet, wenn Differentialgleichungen in verschiedenen Bereichen der Physik gelöst werden. ZB um die Dispersionsrelation von Magnonen oder Phononen in Festkörpern zu lösen, auch sehr häufig in der QM zB in semiklassischen Näherungen, und so weiter. Vor allem die Beispiele entsprechen linearen Differentialgleichungssystemen 2. Ordnung.

Meine Frage:

  • Ist die Verwendung des Planewave-Ansatzes rein mathematisch motiviert? (dh Form/Art der Differentialgleichung).

  • Wenn man weiß, dass ebene Wellen in Wirklichkeit nicht wirklich existieren, da sie eine unendliche Ausbreitung erfordern würden, wie kommt es, dass es immer noch ein theoretisch geeigneter Ansatz ist (zur Beschreibung physikalischer Phänomene)? (z. B. Phononen, die sich sicher bis ins Unendliche ausbreiten) Ich denke, diese letzte Frage könnte man auch anders formulieren: In welchen Fällen kann man sicher sein, dass ein Planewave-Ansatz (Typenlösung) absolut scheitern würde?

Antworten (3)

In vielen Fällen werden unsere Systeme durch lineare Differentialgleichungen beschrieben, und diese haben die Eigenschaft, dass jede Linearkombination von Lösungen der Differentialgleichung auch eine Lösung der Differentialgleichung ist.

Dies ist nützlich, da normalerweise jede beliebige Lösung Fourier-transformiert werden kann, um sie als Summe ebener Wellen auszudrücken. Wenn wir also ebene Wellenlösungen für unsere LDE finden, können wir sie kombinieren, um kompliziertere Lösungen zu schaffen.

Exakt. Wenn es richtig gemacht wird, verlieren wir weder physikalisch noch mathematisch etwas, weil unsere "Vermutung" lediglich auf eine vernünftige Wahl der Basis für unseren Funktionsraum hinausläuft.

John hat dies teilweise beantwortet, jedoch fehlt die grundlegende mathematische Idee:

Wenn wir uns vorstellen, dass eine Funktion Mitglied eines Vektorraums ist, dann existieren Basisvektoren. Dieses Konzept ist Ihnen sicherlich aus der Quantenmechanik bekannt.

Aber auch von dort erinnern wir uns, dass Probleme viel einfacher zu handhaben waren, wenn wir die Eigenbasis der Operatoren kennen, die wir auf Funktionen werfen. Dann ändern wir einfach die Basis der Funktionen in die Operatoren Eigenbasis und fertig.

Wie machen wir das? Wie konstruieren wir den Eigenraum für Differentialoperatoren, die Differentialgleichungen bilden? Lassen Sie uns darüber nachdenken ... Lassen Sie uns zuerst das Eigenwertproblem in abstrakter Notation formulieren

Ö ^ F = C F
und setzen Sie dann das ein, was wir verwenden Ö ^ = D / D X :
D F D X = C F
Und wir sehen sofort, dass dies eine bekannte Gleichung ist! Wofür ist die Lösung F e X P ( C X )

Daher wissen wir, dass, wenn wir uns mit linearen Differentialoperatoren befassen, eine exponentielle Basis unsere Berechnungen immer erleichtern wird. Alles, was wir tun müssen, ist, die richtigen Argumente zu verwenden C um die Reihenfolge der Gleichung anzupassen (dh addieren Sie einige ich s oder ± ).

Die Freiheit in der Wahl C denn diese noch funktionierende Methode hat dann nicht nur zur Fourier-Transformation und ihren von John beschriebenen Eigenschaften geführt, sondern auch zur Laplace-Transformation .

Durch Wählen C passenderweise machen wir natürlich auch eine 'begründete Vermutung' über die Form unserer exponentiellen Lösung. Zum Beispiel ein DEQ

T 2 F + X 2 F = 0
geben oszillierende Lösungen für C = ich ( k X ± ω T ) , sondern etwas, das einen "Startpunkt" hat und dann für echte cs zerfällt oder wächst. Das ist die Idee einer Laplace-Transformation. Sie können die Beispiele auch auf der Wikipage nachschlagen (nach unten scrollen).


Zu deiner anderen Frage:

Der Plane-Wave-Ansatz versagt sicherlich, sobald Nichtlinearitäten in die Differentialgleichung eingehen, wie diese maximal vereinfachte 1-D Navier-Stokes-Gleichung:

T v + v X v = 0
Schon jeder Versuch, eine Fourier-Transformation zu verwenden v führt nicht zu einer algebraisch lösbaren Gleichung.
Numerische Modellierer gehen dieses Problem ee durch pseudospektrale Verfahren an, bei denen die gesamte Gleichung mit Ausnahme von Nichtlinearitäten Fourier-transformiert wird. Diese werden im Realraum berechnet und später numerisch transformiert.

Das andere von Ihnen erwähnte Problem der Unendlichkeit kann auch als Skalenproblem angesehen werden: Mit einer Box-Größe L auf dem Sie Ihr Problem lösen wollen, erhalten Sie immer Diskretisierungseffekte der Wellen, die in dieser Box herumwackeln. Wenn jedoch die Wellenzahlen k << L dann ist deine Wellendichte hoch genug, um noch ein physikalisch sinnvolles Ergebnis erzielen zu können. Die Korrektheit steigt dann, je besser obige Bedingung erfüllt ist.

Vielen Dank für diese ausführliche Antwort, Sie haben alles angesprochen, wonach ich gefragt hatte. Ein bisschen, das ich nicht verstanden habe, war die Schlussfolgerung des ersten Teils Ihrer Antwort, nämlich: "Die Freiheit bei der Wahl von c, damit diese Methode immer noch funktioniert, hat dann nicht nur zur Fourier-Transformation und ihren von John beschriebenen Eigenschaften geführt, sondern auch zur Laplace-Transformation." Könnten Sie vielleicht etwas näher erläutern, warum Sie sagen, dass es zu einer Fourier- und einer Laplace-Transformation geführt hat? vielen Dank.
Erledigt. Sag mir, wenn du noch etwas wissen willst.
Danke für die Bearbeitung. Sollte die Bedingung im letzten Absatz nicht eher sein λ << L ?

Es gibt einige Situationen, in denen der Ebene-Wellen-Ansatz nützlich ist. Es ist eine Lösung für viele Wellengleichungen. Ebene Wellen sind ebenfalls bekannt und wir haben mathematische Techniken, um sie zu handhaben, wie zum Beispiel Fourier-Reihen.

Allerdings kann man mit dem Planewave-Ansatz in manchen Situationen ganz schön daneben liegen. Insbesondere wenn die Welle mit einer Struktur interagiert, die Merkmale in der Größenordnung der Wellenlänge aufweist, können Sie ein Problem haben. Ich werde den Fall von Beugungsgittern erörtern, wo dieses Thema ausführlich diskutiert wurde. Einiges davon kann auf jedes Problem zutreffen, an dem Sie interessiert sind.

Betrachten Sie den Fall eines auf ein Beugungsgitter einfallenden elektromagnetischen Feldes mit einer Periode ähnlich der des einfallenden Feldes und einer Tiefe ähnlich seiner Wellenlänge. Das Feld in der Nähe dieser Struktur hat oft eine ebene Wellenausdehnung, die sich nicht gut verhält. Angenommen, die Ausdehnung der ebenen Welle hat Amplituden a J für J = 1 . Die Folge von a J ist konvergent, wenn lim J = a für einige endlich a . Wenn Sie eine endliche Menge dieser Terme nehmen und sie addieren, dann haben Sie eine Reihe A N = J = 1 N a J . Diese Reihe ist konvergent, falls lim N = A für einige endlich A . Im Allgemeinen kann die Folge konvergieren, während die Reihe divergiert oder umgekehrt. Oder sie konvergieren zu unterschiedlichen Werten. Sehen

JP Hugonin, R. Petit und M. Cadilhac, "Expansion ebener Wellen zur Beschreibung des von einem Gitter gebeugten Feldes", JOSA, Vol. 3, No. 71, Heft 5, S. 593-598 (1981)

für weitere Diskussionen.

In Fällen, in denen ebene Wellenmethoden zur Berechnung von interessierenden Größen nicht funktionieren, können andere Methoden nützlich sein. Beispielsweise das Chandezon-Verfahren zum Berechnen des Feldes, das durch ein Gitter mit einem glatten Profil ungeachtet der Tiefe des Gitters gebeugt wird. Das Chandezon-Verfahren beinhaltet das Umschreiben der Maxwell-Gleichungen in einem Koordinatensystem, in dem das Gitter flach ist, und das anschließende Erweitern der Welle in diesem Koordinatensystem. Für eine Diskussion darüber, wie die Chandezon-Methode implementiert wird, siehe

Lifeng Li, Jean Chandezon, Gérard Granet und Jean-Pierre Plumey, "Rigorose und effiziente Gitteranalysemethode leicht gemacht für Optikingenieure", Applied Optics, Vol. 38, Ausgabe 2, S. 304–313 (1999).

Für einen Vergleich zwischen dem Chandezon-Verfahren und dem Expansionsverfahren mit ebenen Wellen siehe

http://www.ip.univ-bpclermont.fr/Personnel/Kofi.Edee/rayleih_method_C.pdf .