Speeding up a mass-lumped tetrahedral finite-element method for wave propagation

Conference paper (2019)

Authors

W.A. Mulder Shell Global Solutions International B.V., Applied Geophysics and Petrophysics - CEG
S. Geevers Universität Wien
Jaap van der Vegt University of Twente
Research Group
Applied Geophysics and Petrophysics (CEG) (TU Delft)
Copyright: © 2019 W.A. Mulder, S. Geevers, J. van der Vegt
DOI: https://doi.org/10.3997/2214-4609.201900659
More Info
expand_more

Published Date

2019

Language

English

Reuse Rights

Other than for strictly personal use, it is not permitted to download, forward or distribute the text or part of it, without the consent of the author(s) and/or copyright holder(s), unless the work is under an open content license such as Creative Commons.

Faculty

Civil Engineering & Geosciences

Department

Geoscience and Engineering

Research Group

Applied Geophysics and Petrophysics

Abstract

Mass-lumped finite elements on tetrahedra offer more flexibility than their counterpart on hexahedra for the simulation of seismic wave propagation, but there is no general recipe for their construction, unlike as with hexahedra. Earlier, we found new elements up to degree 4 that have significantly less nodes than previously known elements by sharpening the accuracy criterion. A similar approach applied to numerical quadrature of the stiffness matrix provides a speed improvement in the acoustic case and an additional factor 1.5 in the isotropic elastic case. We present numerical results for a homogeneous and heterogeneous isotropic elastic test problem on a sequence of successively finer meshes and for elements of degrees 1 to 4. A comparison of their accuracy and computational efficiency shows that a scheme of degree 4 has the best performance when high accuracy is desired, but the one of degree 3 is more efficient at intermediate accuracy.

Files

Apure.pdf
(.pdf | 0.595 Mb)
Tu_P02_11.pdf
(.pdf | 0.807 Mb)