Dieuwertje Alblas - MIA
Riccardo Bardin - MACS
Vincent Bosboom - MACS
Nicoló Botteghi - MIA
Xiaoyu Cheng - MACS
Giacomo Cristinelli - MIA
Sven Dummer - MIA
Sagy Ephrati - MMS
Arnout Franken - MMS
Elena Giamatteo - MACS
Leonardo del Grande - MIA
source: http://www.malinc.se/math/trigonometry/geocentrismen.php - Heeringa - MIA
Lucas Jansen Klomp - MIA
Muhammad Hamza Khalid - MACS
Nishant Kumar - MACS
Kaifang Liu - MACS
Xiangyi Meng - MACS
Nida Mir - MIA / MDI-TNW
Hongliang Mu - MAST
Kevin Redosado - 3MS
Julian Suk - MIA
Hannah van Susteren - MIA
Elina Thibeua-Sutre - MIA
Alexander Wierzba - MAST
Jens de Vries - MAST
Fengna Yan - MACS
Weihao Yan - MIA

Elastic wave propagation modelling in complex geological models

Organization:

Funded by: Shell

PhD: Sjoerd Geevers

Supervisor: Jaap van der Vegt

Collaboration

Description:

Numerical techniques for simulating seismic waves is an important aspect of predicting earth quakes, but also plays a central role in exploration geophysics. The propagation of these waves is modeled by the linear elastic wave equations, for which several numerical methods have already been developed, including finite difference methods, finite element methods and spectral methods. However, for many practical applications there is still a need for more efficient and flexible methods. Some of the challenges in practical applications are dealing with large domains, complex boundaries and material interfaces and detailed local structures. Large domains result in high computational costs and long running times, while complex interfaces require a flexible mesh. Small local structures require local refinement of the mesh and time step size. A Discontinuous Galerkin (DG) finite element method, combined with a suitable time integration method, should meet these challenges. Discontinuous Galerkin methods can deal with unstructured nonconforming-meshes and allow local refinement of the approximation order (p) as well as mesh size (h), known as hp-refinement. This makes the method suitable for dealing with complex domains, while this gets troublesome for standard finite difference methods. Since the method is based on discontinuous functions with support on only a single element, the mass matrix is block-diagonal and therefore solving the mass matrix equations for explicit time integration methods is significantly faster than with conforming finite element methods. The DG method can also deal with arbitrary high order polynomial basis functions on any standard element shape, while this is more challenging for mass-lumped finite element methods. In this research I will develop and analyze a suitable DG method for the anisotropic elastic wave equations and construct a suitable time integration method for dealing with detailed local structures. The code for this method will be implemented in our software package hpGEM and should be able to run parallel, such that it can run on large supercomputers


Publications
Loading...