Lösung der elastischen Wellengleichung auf einem variablen FD-Gitter Daniel Köhn und Thomas Bohlen...

23
Lösung der elastischen Lösung der elastischen Wellengleichung auf einem variablen Wellengleichung auf einem variablen FD-Gitter FD-Gitter Daniel Köhn und Thomas Bohlen Daniel Köhn und Thomas Bohlen Graz, den 24. Januar 2005 Graz, den 24. Januar 2005

Transcript of Lösung der elastischen Wellengleichung auf einem variablen FD-Gitter Daniel Köhn und Thomas Bohlen...

Page 1: Lösung der elastischen Wellengleichung auf einem variablen FD-Gitter Daniel Köhn und Thomas Bohlen Graz, den 24. Januar 2005.

Lösung der elastischen Wellengleichung auf Lösung der elastischen Wellengleichung auf einem variablen FD-Gittereinem variablen FD-Gitter

Daniel Köhn und Thomas BohlenDaniel Köhn und Thomas Bohlen

Graz, den 24. Januar 2005Graz, den 24. Januar 2005

Page 2: Lösung der elastischen Wellengleichung auf einem variablen FD-Gitter Daniel Köhn und Thomas Bohlen Graz, den 24. Januar 2005.

Einführung in die adaptive FD-Modellierung

Kantenlänge LKantenlänge L

L/2L/2

Vs1

Vs2 >> Vs1

Page 3: Lösung der elastischen Wellengleichung auf einem variablen FD-Gitter Daniel Köhn und Thomas Bohlen Graz, den 24. Januar 2005.

FD-Diskretisierung auf einem homogenen Gitter

Vs1

Beachte Gitterdispersion:dh 1/nmit1=Vs1/fmax

n = 12 (2.Ordnung) = 8 (4.Ordnung)

Vs2>> Vs1

L/2L/2

„ oversampled “

Page 4: Lösung der elastischen Wellengleichung auf einem variablen FD-Gitter Daniel Köhn und Thomas Bohlen Graz, den 24. Januar 2005.

FD-Diskretisierung auf einem adaptiven Gitter

Vs1

Beachte Gitterdispersion:dh 1/nmit1=Vs1/fmax

n = 12 (2.Ordnung) = 8 (4.Ordnung)

Vs2>> Vs1

L/2L/2

Page 5: Lösung der elastischen Wellengleichung auf einem variablen FD-Gitter Daniel Köhn und Thomas Bohlen Graz, den 24. Januar 2005.

Rechenzeitersparnis

Rechenzeit 2D: 1/25

Rechenzeit 3D: 1/125

Page 6: Lösung der elastischen Wellengleichung auf einem variablen FD-Gitter Daniel Köhn und Thomas Bohlen Graz, den 24. Januar 2005.

Implementierung des adaptiven FD-Codes

Nach Jastram (1992)

Page 7: Lösung der elastischen Wellengleichung auf einem variablen FD-Gitter Daniel Köhn und Thomas Bohlen Graz, den 24. Januar 2005.

Testproblem:homogener Raum (1x1 m), umgeben von Luft

100x100 Gitterpunktedt = 3e-6 sRechenzeit: 40000 Zeitschritte

Test: 2D-Modellierung eines homogenen Raumes

Verwendeter FD-Code: - Velocity-Stress Formulierung der elastischen Wellengleichung- Standard Staggered Grid (SSG)- FD-Operatoren 2.Ordnung

=> zu interpolierende Variablen: syy , vx

Page 8: Lösung der elastischen Wellengleichung auf einem variablen FD-Gitter Daniel Köhn und Thomas Bohlen Graz, den 24. Januar 2005.

Test: 2D-Modellierung eines homogenen Raumes

Page 9: Lösung der elastischen Wellengleichung auf einem variablen FD-Gitter Daniel Köhn und Thomas Bohlen Graz, den 24. Januar 2005.

Analyse der Instabilität

Page 10: Lösung der elastischen Wellengleichung auf einem variablen FD-Gitter Daniel Köhn und Thomas Bohlen Graz, den 24. Januar 2005.

Definition des Instabilitätszeitpunktes

Page 11: Lösung der elastischen Wellengleichung auf einem variablen FD-Gitter Daniel Köhn und Thomas Bohlen Graz, den 24. Januar 2005.

Einfluss unterschiedlicher Interpolationsverfahren auf den Instabilitätszeitpunkt

Hypothese: Addition der Interpolationsfehler führt zur Instabilität ?

Interpolationsverfahren Instabilität nach L1-Norm [Zeitschritte]

gesamte Rechenzeit [s]

homogenes Gitter 8 332.1nearest neighbor 4233 368.1trigonometrisch 16277 401.03

linear 17600 247.6kubisch 19000 249.4

kubische Hermite Polynome 19789 355.1kubische Splines 20011 420.3

Page 12: Lösung der elastischen Wellengleichung auf einem variablen FD-Gitter Daniel Köhn und Thomas Bohlen Graz, den 24. Januar 2005.

Modellierung des Interpolationsfehlers durch multiplikativen Gauss’schen Noise

Syy(j,i) = Syy(j,i) (1+) vx(j,i) = vx(j,i) (1+ wobei, [-0.02, 0.02]

Page 13: Lösung der elastischen Wellengleichung auf einem variablen FD-Gitter Daniel Köhn und Thomas Bohlen Graz, den 24. Januar 2005.

Modellierung des Interpolationsfehlers durch multiplikativen Gauss’schen Noise

=> Instabilität ist nicht allein auf Interpolationsfehler zurückzuführen.

Entwicklung der L1-Norm (vx, multiplikativer Noise)

Page 14: Lösung der elastischen Wellengleichung auf einem variablen FD-Gitter Daniel Köhn und Thomas Bohlen Graz, den 24. Januar 2005.

Analyse des Wellenzahl-Spektrums im Übergangsbereich zwischen groben und feinem

Gitter

Page 15: Lösung der elastischen Wellengleichung auf einem variablen FD-Gitter Daniel Köhn und Thomas Bohlen Graz, den 24. Januar 2005.

Wellenzahl-Spektrum und L1-Norm als Funktion der Zeit

Page 16: Lösung der elastischen Wellengleichung auf einem variablen FD-Gitter Daniel Köhn und Thomas Bohlen Graz, den 24. Januar 2005.

Der Einfluß des groben Gitters

Nyquist-Wellenzahl feines Gitter

Nyquist-Wellenzahl grobes Gitter

Quellsignal

„Coarse Grid Nyquist Peaks“

Page 17: Lösung der elastischen Wellengleichung auf einem variablen FD-Gitter Daniel Köhn und Thomas Bohlen Graz, den 24. Januar 2005.

Stabilisierung des adaptiven FD-Codes

Page 18: Lösung der elastischen Wellengleichung auf einem variablen FD-Gitter Daniel Köhn und Thomas Bohlen Graz, den 24. Januar 2005.

Response des 1D-Butterworth-Filters zur Unterdückung der Response des 1D-Butterworth-Filters zur Unterdückung der “Coarse Grid Nyquist Peaks”“Coarse Grid Nyquist Peaks”

Nyquist-Wellenzahl feines Gitter

Nyquist-Wellenzahl grobes Gitter

Page 19: Lösung der elastischen Wellengleichung auf einem variablen FD-Gitter Daniel Köhn und Thomas Bohlen Graz, den 24. Januar 2005.

Unterdrückung des Noise durch Tiefpass-FilterungUnterdrückung des Noise durch Tiefpass-Filterung

Page 20: Lösung der elastischen Wellengleichung auf einem variablen FD-Gitter Daniel Köhn und Thomas Bohlen Graz, den 24. Januar 2005.

Unterdrückung des Noise durch Tiefpass-FilterungUnterdrückung des Noise durch Tiefpass-Filterung

Page 21: Lösung der elastischen Wellengleichung auf einem variablen FD-Gitter Daniel Köhn und Thomas Bohlen Graz, den 24. Januar 2005.

Rechenzeitersparnis: TestproblemRechenzeitersparnis: Testproblem

Homogenes Gitter ... 433.09 sAdaptives Gitter ... 293.77 s

Relative Rechenzeit adaptiv/homogen ... 67 %Theoretisch maximal möglich ... 63 %

Page 22: Lösung der elastischen Wellengleichung auf einem variablen FD-Gitter Daniel Köhn und Thomas Bohlen Graz, den 24. Januar 2005.

Zusammenfassung

Die Amplitude des Interpolationsfehlers (multiplikativer Noise) wächst während der Rechnung stetig an.

Das Auftreten der „Coarse Grid Nyquist Peaks“ führt aufgrund der Verletzung des Nyquist-Kriteriums des groben Gitters zur Entstehung einer numerischen Instabilität.

Diese wächst exponentiell mit der Zeit an.

Eine Tiefpassfilterung des k-Bereichs unterhalb der Nyquist-Wellenzahl des groben Gitters führt zu einer Stabilisierung des adaptiven FD-Codes.

Page 23: Lösung der elastischen Wellengleichung auf einem variablen FD-Gitter Daniel Köhn und Thomas Bohlen Graz, den 24. Januar 2005.

Ausblick

Implementierug einer Filterung im Orts-Bereich (Parallelisierung).

Genauigkeitstest anhand geologischer Beispiele:

- Probleme mit hohen Vp/Vs-Verhältnissen - Auflösung kleinskaliger Strukturen

3D-Parallelisierung