Computational Sciences for Energy Research (15CSER049)

Port-Hamiltonian Discontinuous Galerkin Methods for seismic wave computations 

Organization:

Funded by:	Shell/FOM/NWO Computational Science in Energy Research Program 15CSER049
PhD:
Supervisors:	prof.dr.ir. J.J.W. van der Vegt (Jaap) Professor in Mathematics of Computational Mechanics +31534895628  j.j.w.vandervegt@utwente.nl Building: Zilverling 3034 Personal page prof.dr. H.J. Zwart (Hans) Professor in Physical Systems and Control/Chair Mathematics of Systems Theory +31534893464  h.j.zwart@utwente.nl Building: Zilverling 3042 Personal page
Collaboration:	-

Description:

We aim to increase the accuracy of seismic wave computations in heterogeneous elastic media that may contain cracks, holes, thin layers and other local structures by developing discontinuous Galerkin (DG) discretizations that preserve the mathematical structure of the (an)isotropic elastic wave equations. To this purpose we combine a DG finite element method with port-Hamiltonian system theory. This mathematical systems approach provides a new framework to derive DG discretizations that have small dispersion and dissipation errors in seismic wave propagation on general unstructured meshes. 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.

Output:

Pictures: