Efficient Quadrature Rules for Computing the Stiffness Matrices of Mass-Lumped Tetrahedral Elements for Linear Wave Problems

Journal Article (2019)
Author(s)

S. Geevers (University of Twente)

W. Mulder (TU Delft - Applied Geophysics and Petrophysics, Shell Global Solutions International B.V.)

J. J.W. van der Vegt (University of Twente)

Research Group
Applied Geophysics and Petrophysics
Copyright
© 2019 S. Geevers, W.A. Mulder, J.J.W. van der Vegt
DOI related publication
https://doi.org/10.1137/18M1198557
More Info
expand_more
Publication Year
2019
Language
English
Copyright
© 2019 S. Geevers, W.A. Mulder, J.J.W. van der Vegt
Research Group
Applied Geophysics and Petrophysics
Issue number
2
Volume number
41
Pages (from-to)
A1041–A1065
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.

Abstract

We present new and efficient quadrature rules for computing the stiffness matrices of mass-lumped tetrahedral elements for wave propagation modeling. These quadrature rules allow for a more efficient implementation of the mass-lumped finite element method and can handle materials that are heterogeneous within the element without loss of the convergence rate. The quadrature rules are designed for the specific function spaces of recently developed mass-lumped tetrahedra, which consist of standard polynomial function spaces enriched with higher-degree bubble functions. For the degree-2 mass-lumped tetrahedron, the most efficient quadrature rule seems to be an existing 14-point quadrature rule, but for tetrahedra of degrees 3 and 4, we construct new quadrature rules that require fewer integration points than those currently available in the literature. Several numerical examples confirm that this approach is more efficient than computing the stiffness matrix exactly and that an optimal order of convergence is maintained, even when material properties vary within the element.