Computational Sciences for Energy Research (13CSER052)

Computing seismic waves with minimal pollution error


Funded by:

Shell/FOM/NWO Computational Science in Energy Research Program 13CSER052


Poorvi Shukla


Jaap van der Vegt




We aim to reduce the pollution error in the computation of seismic elastic waves in heterogeneous media that may contain cracks, holes, thin layers and other local structures. To this purpose a novel space-time Trefftz discontinuous Galerkin (DG) finite element method will be developed. In this numerical technique the higher frequency part of seismic waves is represented element-wise using a Green's function approach, whereas global wave interactions are accounted for by the discontinuous Galerkin method. This approach will significantly reduce the pollution error, which currently limits time-dependent seismic wave simulations to relatively low frequencies. The discontinuous basis functions in the DG method are also extremely useful for the accurate representation of interfaces with different material properties and to capture important features that are smaller than the dominant seismic wave length, but still affect seismic waves, using hp-adaptive mesh refinement