### Abstract

Original language | Undefined |
---|---|

Pages (from-to) | 112-121 |

Number of pages | 10 |

Journal | Physics of the earth and planetary interiors |

Volume | 171 |

Issue number | 1-4 |

DOIs | |

Publication status | Published - 2008 |

### Keywords

- METIS-254768
- Wave propagation
- Numerical methods
- Scattering
- Finite difference method
- Finite Element Method
- Analytical solution
- IR-79094

### Cite this

*Physics of the earth and planetary interiors*,

*171*(1-4), 112-121. https://doi.org/10.1016/j.pepi.2008.07.003

}

*Physics of the earth and planetary interiors*, vol. 171, no. 1-4, pp. 112-121. https://doi.org/10.1016/j.pepi.2008.07.003

**Comparison of finite difference and finite element methods for simulating two-dimensional scattering of elastic waves.** / Frehner, Marcel; Schmalholz, Stefan M.; Saenger, Erik H.; Steeb, Holger Karl.

Research output: Contribution to journal › Article › Academic › peer-review

TY - JOUR

T1 - Comparison of finite difference and finite element methods for simulating two-dimensional scattering of elastic waves

AU - Frehner, Marcel

AU - Schmalholz, Stefan M.

AU - Saenger, Erik H.

AU - Steeb, Holger Karl

PY - 2008

Y1 - 2008

N2 - Two-dimensional scattering of elastic waves in a medium containing a circular heterogeneity is investigated with an analytical solution and numerical wave propagation simulations. Different combinations of finite difference methods (FDM) and finite element methods (FEM) are used to numerically solve the elastodynamic wave equations. Finite difference and finite element techniques are applied to approximate both the time and space derivatives and are combined in various ways to provide different numerical algorithms for modeling elastic wave propagation. The results of the different numerical algorithms are compared for simulations of an incident plane P-wave that is scattered by a mechanically weak circular inclusion whereby the diameter of the inclusion is of the same order than the P-wave's wavelength. For this scattering problem an analytical solution is available and used as the reference solution in the comparison of the different numerical algorithms. Staircase-like spatial discretization of the inclusion's circular shape with the finite difference method using a rectangular grid provides accurate velocity and displacement fields close to the inclusion boundary only for very high spatial resolutions. Implicit time integration based on either finite differences or finite elements does not provide computational advantages compared to explicit schemes. The best numerical algorithm in terms of accuracy and computation time for the investigated scattering problem consists of a finite element method in space using an unstructured mesh combined with an explicit finite difference method in time. The computational advantages and disadvantages of the different numerical algorithms are discussed.

AB - Two-dimensional scattering of elastic waves in a medium containing a circular heterogeneity is investigated with an analytical solution and numerical wave propagation simulations. Different combinations of finite difference methods (FDM) and finite element methods (FEM) are used to numerically solve the elastodynamic wave equations. Finite difference and finite element techniques are applied to approximate both the time and space derivatives and are combined in various ways to provide different numerical algorithms for modeling elastic wave propagation. The results of the different numerical algorithms are compared for simulations of an incident plane P-wave that is scattered by a mechanically weak circular inclusion whereby the diameter of the inclusion is of the same order than the P-wave's wavelength. For this scattering problem an analytical solution is available and used as the reference solution in the comparison of the different numerical algorithms. Staircase-like spatial discretization of the inclusion's circular shape with the finite difference method using a rectangular grid provides accurate velocity and displacement fields close to the inclusion boundary only for very high spatial resolutions. Implicit time integration based on either finite differences or finite elements does not provide computational advantages compared to explicit schemes. The best numerical algorithm in terms of accuracy and computation time for the investigated scattering problem consists of a finite element method in space using an unstructured mesh combined with an explicit finite difference method in time. The computational advantages and disadvantages of the different numerical algorithms are discussed.

KW - METIS-254768

KW - Wave propagation

KW - Numerical methods

KW - Scattering

KW - Finite difference method

KW - Finite Element Method

KW - Analytical solution

KW - IR-79094

U2 - 10.1016/j.pepi.2008.07.003

DO - 10.1016/j.pepi.2008.07.003

M3 - Article

VL - 171

SP - 112

EP - 121

JO - Physics of the earth and planetary interiors

JF - Physics of the earth and planetary interiors

SN - 0031-9201

IS - 1-4

ER -