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

Marcel Frehner, Stefan M. Schmalholz, Erik H. Saenger, Holger Karl Steeb

Research output: Contribution to journalArticleAcademicpeer-review

31 Citations (Scopus)

Abstract

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.
Original languageUndefined
Pages (from-to)112-121
Number of pages10
JournalPhysics of the earth and planetary interiors
Volume171
Issue number1-4
DOIs
Publication statusPublished - 2008

Keywords

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

Cite this

@article{8fbdca0af76e483e87a1031e53a5d097,
title = "Comparison of finite difference and finite element methods for simulating two-dimensional scattering of elastic waves",
abstract = "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.",
keywords = "METIS-254768, Wave propagation, Numerical methods, Scattering, Finite difference method, Finite Element Method, Analytical solution, IR-79094",
author = "Marcel Frehner and Schmalholz, {Stefan M.} and Saenger, {Erik H.} and Steeb, {Holger Karl}",
year = "2008",
doi = "10.1016/j.pepi.2008.07.003",
language = "Undefined",
volume = "171",
pages = "112--121",
journal = "Physics of the earth and planetary interiors",
issn = "0031-9201",
publisher = "Elsevier",
number = "1-4",

}

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.

In: Physics of the earth and planetary interiors, Vol. 171, No. 1-4, 2008, p. 112-121.

Research output: Contribution to journalArticleAcademicpeer-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 -