Efficiency of Old and New Triangular Finite Elements for Wavefield Modelling in Time

Finite elements with mass lumping allow for explicit time stepping when modelling wave propagation and can be more efficient than finite differences in complex geological settings. In 2D on quadrilaterals, spectral elements are the obvious choice. Triangles are more flexible for meshing, but the construction of polynomial elements is less straightforward. So far, elements up to degree 9 have been found. Some years ago, an accuracy criterion that is sharper and less restrictive than the customary one led to new tetrahedral elements that are considerably more efficient than those previously known. Applying the same criterion to triangular elements provides infinitely many new elements of degree 5, with the same number of nodes as the old one, and two elements of degree 6 with less nodes than the known ones. Their efficiency, measured in terms of the compute time needed to obtain a solution with a given accuracy, is determined for a homogeneous problem and compared to that of the old elements of degree 1 to 8. For moderate accuracy, elements of degree 3 are the most efficient. For high accuracy, one of the new degree-6 elements performs best.Green Open Access added to TU Delft Institutional Repository 'You share, we take care!' - Taverne project https://www.openaccess.nl/en/you-share-we-take-care Otherwise as indicated in the copyright section: the publisher is the copyright holder of this work and the author uses the Dutch legislation to make this work public.Applied Geophysics and Petrophysic