• ggy337

    Final published version, 1.69 MB, PDF document


Finite elements can, in some cases, outperform finite-difference methods for modelling wave propagation in complex geological models with topography. In the weak form of the finiteelement method, the delta function is a natural way to represent a point source. If, instead of the usual second-order form, the first-order form of the wave equation is considered, this is no longer true. Fourier analysis for a simple case shows that the spatial operator corresponding to
the first-order form has short-wavelength null-vectors. Once excited, these modes are not seen by the spatial operator but only by the time- stepping scheme and show up as noise. A sourcewith a larger spatial extent, for instance a Gaussian or a tapered sinc, can avoid the excitation of problematic short wavelengths. A series of numerical experiments on a 2-D problem with an exact solution provides a suggestion for the best choice of parameters for these source
term distributions. The tapered sinc provided the best results and the resulting accuracy can be better than that of the second-order form. The higher operation count of the former, however, does not make it more efficient in terms of accuracy for a given computational effort, at least not for the 2-D examples considered here
Original languageEnglish
Pages (from-to)1231–1240
JournalGeophysical Journal International
Issue number2
Publication statusPublished - 2018

ID: 46748164