The second-order formulation of the wave equation is often used for spectral-element discretizations. For some applications, however, a first-order formulation may be desirable. It can, in theory, provide much better accuracy in terms of numerical dispersion if the consistent mass matrix is used and the degree of the polynomial basis functions is odd. However, we find in the 1-D case that the eigenvector errors for elements of degree higher than one are larger for the first-order than for the second-order formulation. These errors measure the unwanted cross talk between the different eigenmodes. Since they are absent for the lowest degree, that linear element may perform better in the first-order formulation if the consistent mass matrix is inverted. The latter may be avoided by using one or two defect-correction iterations. Numerical experiments on triangles confirm the superior accuracy of the first-order formulation. However, with a delta-function point source, a large amount of numerical noise is generated. Although this can be avoided by a smoother source representation, its higher cost and the increased susceptibility to numerical noise make the second-order formulation more attractive.
Original languageEnglish
Title of host publication78th EAGE Conference and Exhibition 2016
Subtitle of host publicationVienna, Austria
Number of pages5
Publication statusPublished - 2016
Event78th EAGE Conference and Exhibition 2016 - Messe Wien, Exhibition and Congress Center, Vienna, Austria
Duration: 30 May 20162 Jun 2016
Conference number: 78


Conference78th EAGE Conference and Exhibition 2016
Abbreviated titleEAGE 2016
Internet address

ID: 4700064