1 / 29

Session : Computational Wave Propagation: Basic Theory

Igel H., Fichtner A., Käser M., Virieux J., Seriani G., Capdeville Y., Moczo P. The finite-difference method The pseudo-spectral method. Session : Computational Wave Propagation: Basic Theory. The problem. The wave equation.

hea
Download Presentation

Session : Computational Wave Propagation: Basic Theory

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Igel H., Fichtner A., Käser M., Virieux J., Seriani G., Capdeville Y., Moczo P. • The finite-difference method • The pseudo-spectral method Session: Computational Wave Propagation: Basic Theory

  2. The problem

  3. The wave equation Displacement u, elastic parameters E(x) vary only in one direction. Density r and Lame parameters l and m. shear waves P waves with the corresponding velocities shear waves P waves

  4. … we disregard … 3-dimensionality External forces Viscoelastic behaviour Anisotropic behaviour Non-linear elasticity Etc. … and even further simplify and re-write to …

  5. … motion of a string … Displacement u, shear modulus m(x), density r vary only in one direction x.

  6. The velocity-stress formulation stress-strain Generalized Interpolation Low-order to spectral accuracy Extrapolation usually low-order

  7. Other forms of the elastic wave equation Strong form (see Moczo et al.) Displacement – stress Displacement – velocity – stress Velocity – stress Displacement Weak form (more later) Projection on test functions

  8. Numerical methods Finite differences Pseudospectral methods Finite (spectral) elements (Fichtner) Discontinuous Galerkin (Käser) Frequency domain methods (Virieux) Cell methods (Seriani) (Boundary integral methods, lattice solids, discrete wavenumber, finite volumes, etc).

  9. Finite Differences forward difference backward difference centered difference

  10. ... or even … … in space and for time … This turns out to be a useful approximation!

  11. Discrete space-time Space and time discretization Regular, 1D

  12. Our first FD wave algorithm Space discretization Regular, 1D

  13. … and the explicit extrapolation scheme … Source term Do for all times … Stop!

  14. … staggered grids … ? the only unknown l+1/2  Time l    l-1/2  m-1/2 m m+1/2 Space

  15. Results

  16. Numerical dispersion Amplitude Time

  17. Von Neumann Analysisplane wave analysis continuous discrete Conditional stability

  18. Numerical dispersion phase velocity - group velocity Theoretical prediction of phase and group velocities as a function of wavelength. In 2-D, 3-D this effect depends on direction (-> grid anisotropy)

  19. Well, let us have a closer look at the space derivatives (because there is not so much we can do about the time derivative …)

  20. The pseudospectral approach in the discrete world that means, we can … Discrete Spectrum Exact derivative

  21. PS Algorithm (displacement) Derivatives calculated by DFT (forward and inverse transform and multiplication with –ik) …. cool algorithm, but …

  22. … some more thoughts on the spectral derivative operator „-ik“ … using the convolution theorem … Differential operator in space domain The function we want to differentiate

  23. „Convolution“ d*f f 0 1 0 0 0 d 1 -2 1 0 1 0 0 0 1 -2 1 0 1 0 0 1 1 -2 1 0 1 0 0 -2 1 -2 1 0 1 0 0 1 1 -2 1 0 1 0 0 0 1 -2 1

  24. The numerical wave number as a function of true wave number (ik as a function of k) This is the derivative operater d(x) and it s discrete version (dots)

  25. Finite difference - Pseudospectral nx-point 2-4-6 point Length of (convolutional) operator Increasing accuracy Increasing floating point operations Fourier (Chebysheff) coefficients Taylor coefficients PS FD

  26. FD in real 3-D applications Cartesian Geometry Spherical coordinates Sections Axisymmetric Rheologies Viscolastic Anisotropic Poroelastic Dynamic rupture

  27. Finite Differences - Pros and Cons PROS simple theory Explicit scheme – no matrix inversion Easy parallelization Easy model generation (on regular grids) Easy adaptation to specific problems CONS Boundary conditions difficult to implement Requires large number of grid points per wavelength (in particular surface waves) Large memory requirement Inefficient for models with strong velocity variations Topography not easily implemented

  28. Pseudospectral Method - Pros and Cons PROS Beautiful! Exact space derivatives Explicit scheme – no matrix inversion Centred scheme (anisotropy) Accurate implementation of boundary condition (Chebyshev) Memory efficient (less points per wavelength w.r.t. FD) CONS Global communication scheme (inefficient parallelization) Irregular grid points (Chebyshev) –> stability problems Boundary conditions difficult for Fourier Method

  29. Final thoughts Some of the fundamental concepts of computational seismology can be understood from FD schemes (e.g., stability, dispersion) Each QUEST researcher should be able to Understand what time steps and grid distances are appropriate for a specific Earth velocity models Understand what affects the accuracy of wave simulations (np/l; operators used, wave type, complexity of velocity model, propagation distance, etc.) Know the possible traps of (community) algorithms (black boxes) Know how to benchmark simulation codes (for accuracy)

More Related