Dieser Artikel oder nachfolgende Abschnitt ist nicht hinreichend mit
Belegen (beispielsweise
Einzelnachweisen) ausgestattet. Angaben ohne ausreichenden Beleg könnten demnächst entfernt werden. Bitte hilf Wikipedia, indem du die Angaben recherchierst und
gute Belege einfügst.
Der Titel dieses Artikels ist mehrdeutig. Zur gleichnamigen Geheimdienstkooperation siehe
Operation Eikonal.
Als Eikonal (Altgriechisch εἰκών eikon = Bild, Abbild) wird in der geometrischen Optik die Strecke eines Lichtstrahls zwischen Ausgangs- und Endpunkt bezeichnet; mittlerweile bezeichnet der Begriff meist das Bruns-Eikonal.
Das Bruns-Eikonal oder Brunssche Eikonal ist eine Funktion, die gemäß dem Fermatschen Prinzip den kürzesten Weg zwischen zwei durch optische Medien getrennten Punkten beschreibt. Sie wurde vom deutschen Mathematiker Heinrich Bruns 1895 veröffentlicht und in der Strahlenoptik benutzt. Der Name Eikonal stammt von Bruns, das Verfahren war aber schon William Rowan Hamilton bekannt, der es charakteristische Funktion nannte (Hamilton-Jacobi-Gleichung) und in Optik und Mechanik anwandte.
Das Bruns-Eikonal wird bei akustischen Wellen und anderen Wellenphänomenen angewendet, z. B. in der Seismologie zur Berechnung der Ausbreitung seismischer Wellen.
Nachfolgend soll die Eikonalgleichung als Hochfrequenzapproximation der akustischen Wellengleichung hergeleitet werden. In der Quantenmechanik wird ein ähnliches Verfahren verwendet, die semiklassische WKB-Näherung.
Wir gehen also von der akustischen Wellengleichung mit dem Druck
, dem Ortsvektor
, der ortsabhängigen Ausbreitungsgeschwindigkeit
und konstanter Dichte aus
![{\displaystyle \nabla ^{2}p-{\frac {1}{c^{2}}}{\frac {\partial ^{2}p}{\partial t^{2}}}=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/137049930ced8a29b4a9211f9cd9c8f9a25cf353)
Gesucht ist ein zeitlich harmonischer Hochfrequenzansatz, für den eine
frequenz- und zeitunabhängige Amplitude
und die Laufzeitfunktion
angenommen werden kann. Sie hat die Form
![{\displaystyle p({\vec {x}},t)=P({\vec {x}})e^{\mathrm {i} \omega (t-\phi ({\vec {x}}))}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/19370f5c99ddc11307170fb9c94e1766ca6f8ee6)
Zunächst berechnet man die Zeitableitungen der Wellengleichung:
![{\displaystyle {\frac {\partial p}{\partial t}}=\mathrm {i} \omega P({\vec {x}})e^{i\omega (t-\phi ({\vec {x}}))};\quad {\frac {\partial ^{2}p}{\partial t^{2}}}=-\omega ^{2}P({\vec {x}})e^{i\omega (t-\phi ({\vec {x}}))}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/59f973e072921d5fd26c05b1b5437857e02ea43c)
Nun folgen die Ortsableitungen:
![{\displaystyle \nabla p=\nabla Pe^{\mathrm {i} \omega (t-\phi ({\vec {x}}))}-\mathrm {i} \omega Pe^{\mathrm {i} \omega (t-\phi ({\vec {x}}))}\nabla \phi =(\nabla P-\mathrm {i} \omega P\nabla \phi )e^{\mathrm {i} \omega (t-\phi ({\vec {x}}))}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3c346eca6822ec5f8bbed0b1361dab2bdc2f7d65)
Wegen der vektoriellen Identität
gilt weiter:
![{\displaystyle \nabla ^{2}p=\nabla \cdot \nabla p}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ebb9f39dcf0b677f82e73c319063eada5095c4fa)
![{\displaystyle =\nabla \cdot (\nabla P-\mathrm {i} \omega P\nabla \phi )e^{\mathrm {i} \omega (t-\phi ({\vec {x}}))}+(\nabla P-\mathrm {i} \omega P\nabla \phi )\cdot \nabla e^{\mathrm {i} \omega (t-\phi ({\vec {x}}))}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/49205116c69f2cc783bc8aafe3ffdf640eb774e3)
![{\displaystyle =\left(\nabla ^{2}P-\mathrm {i} \omega \nabla P\cdot \nabla \phi -\mathrm {i} \omega P\nabla ^{2}\phi -\mathrm {i} \omega (\nabla P-\mathrm {i} \omega P\nabla \phi )\cdot \nabla \phi \right)e^{\mathrm {i} \omega (t-\phi ({\vec {x}}))}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/026ed85cff71c55230325b58fd9a7856a38248e1)
![{\displaystyle =\left(\nabla ^{2}P-2\mathrm {i} \omega \nabla P\cdot \nabla \phi -\mathrm {i} \omega P\nabla ^{2}\phi -\omega ^{2}P(\nabla \phi )^{2}\right)e^{\mathrm {i} \omega (t-\phi ({\vec {x}}))}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d377201d2f2f4e27556ac0bfc463758459f586c1)
Die beiden Ableitungen in die Wellengleichung eingesetzt ergeben nach Division durch
![{\displaystyle -\omega ^{2}P\left((\nabla \phi )^{2}-{\frac {1}{c^{2}}}\right)-\mathrm {i} \omega \left(2\nabla P\cdot \nabla \phi +P\nabla ^{2}\phi \right)+\nabla ^{2}P=0.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e751ef9c14babd3774f739b24f1b2ba70c1c98e0)
Eine Division durch
führt dann zu
![{\displaystyle \left((\nabla \phi )^{2}-{\frac {1}{c^{2}}}\right)+{\frac {\mathrm {i} }{\omega P}}\left(2\nabla P\cdot \nabla \phi +P\nabla ^{2}\phi \right)-{\frac {1}{\omega ^{2}P}}\left(\nabla ^{2}P\right)=0.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e50a681e00ce96e0536defab81edba51234591f3)
Da Real- und Imaginärteil der Gleichung unabhängig voneinander gleich null sein müssen, folgt:
![{\displaystyle \left((\nabla \phi )^{2}-{\frac {1}{c^{2}}}\right)-{\frac {1}{\omega ^{2}P}}\left(\nabla ^{2}P\right)=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bf01ccbb71167f1df9f3542672bf88a8e2acb951)
Bei der Näherung geht man davon aus, dass die Amplitude
nur schwach ortsabhängig,
also beschränkt ist. Da gleichzeitig weder die Laufzeit
noch die Amplitude
frequenzabhängig sind, ist der zweite Term für sehr hohe Frequenzen klein gegenüber dem ersten Term und die Gleichung vereinfacht sich auf:
![{\displaystyle \left(\nabla \phi \right)^{2}={\frac {1}{c^{2}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/08e44f93855c818d21ab98d99ca831ff71195787)
Die Lösung
der Eikonalgleichung ordnet jedem Punkt im Ortsraum die Laufzeit der Welle zu. Linien gleicher Laufzeit lassen sich entsprechend als Wellenfronten interpretieren.