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.@en

We analyse the dispersion properties of two types of explicit finite element methods for modelling acoustic and elastic wave propagation on tetrahedral meshes, namely mass-lumped finite element methods and symmetric interior penalty discontinuous Galerkin methods, both combined with a suitable Lax–Wendroff time integration scheme. The dispersion properties are obtained semi-analytically using standard Fourier analysis. Based on the dispersion analysis, we give an indication of which method is the most efficient for a given accuracy, how many elements per wavelength are required for a given accuracy, and how sensitive the accuracy of the method is to poorly shaped elements.@en

We present a new accuracy condition for the construction of continuous mass-lumped elements. This condition is less restrictive than the one currently used and enabled us to construct new mass-lumped tetrahedral elements of degrees 2 to 4. The new degree-2 and degree-3 tetrahedral elements require 15 and 32 nodes per element, respectively, while currently, these elements require 23 and 50 nodes, respectively. The new degree-4 elements require 60, 61, or 65 nodes per element. Tetrahedral elements of this degree had not been found until now. We prove that our accuracy condition results in a mass-lumped finite element method that converges with optimal order in the $L^2$-norm and energy-norm. A dispersion analysis and several numerical tests confirm that our elements maintain the optimal order of accuracy and show that the new mass-lumped tetrahedral elements are more efficient than the current ones. @en