Grid dispersion and stability criteria of some common finite‐difference and finite‐element methods for acoustic and elastic wave propagation

2007 ◽  
Author(s):  
Jonas D. De Basabe ◽  
Mrinal K. Sen
Keyword(s):  
Finite Element ◽  
Wave Propagation ◽  
Finite Difference ◽  
Elastic Wave ◽  
Finite Element Methods ◽  
Elastic Wave Propagation ◽  
Stability Criteria ◽  
Grid Dispersion

Grid dispersion and stability criteria of some common finite-element methods for acoustic and elastic wave equations

Geophysics ◽  
2007 ◽  
Vol 72 (6) ◽  
pp. T81-T95 ◽  
Author(s):  
Jonás D. De Basabe ◽  
Mrinal K. Sen
Keyword(s):  
Finite Element ◽  
Wave Propagation ◽  
Elastic Wave ◽  
Finite Element Methods ◽  
Elastic Wave Propagation ◽  
Staggered Grid ◽  
Element Approximation ◽  
Stability Criteria ◽  
Grid Dispersion ◽  
Order Four

Purely numerical methods based on finite-element approximation of the acoustic or elastic wave equation are becoming increasingly popular for the generation of synthetic seismograms. We present formulas for the grid dispersion and stability criteria for some popular finite-element methods (FEM) for wave propagation, namely, classical and spectral FEM. We develop an approach based on a generalized eigenvalue formulation to analyze the dispersive behavior of these FEMs for acoustic or elastic wave propagation that overcomes difficulties caused by irregular node spacing within the element and the use of high-order polynomials, as is the case for spectral FEM. Analysis reveals that for spectral FEM of order four or greater, dispersion is less than 0.2% at four to five nodes per wavelength, and dispersion is not angle dependent. New results can be compared with grid-dispersion results of some classical finite-difference methods (FDM) used for acoustic or elastic wave propagation. Analysis reveals that FDM and classical FEM require a larger sampling ratio than a spectral FEM to obtain results with the same degree of accuracy. The staggered-grid FDM is an efficient scheme, but the dispersion is angle dependent with larger values along the grid axes. On the other hand, spectral FEM of order four or greater is isotropic with small dispersion, making it attractive for simulations with long propagation times.

Sign in / Sign up

Export Citation Format

Share Document