High-Order Hybridized Discontinuous Galerkin (HDG) Method for Wave Propagation Simulation in Complex Geophysical Media – Elastic, Acoustic and Hydro-Acoustic – an Unifying Framework to Couple Continuous Spectral Element and Discontinuous Galerkin Methods.

Thursday, 18 December 2014
Terrana Sébastien1,2, Jean-Pierre Vilotte1, Laurent Guillot2 and Christian Mariotti2, (1)Institut de Physique du Globe de Paris, Paris, France, (2)CEA Commissariat à l'Energie Atomique DAM, DIF, Arpajon Cedex, France
Today seismological observation systems combine broadband seismic receivers, hydrophones and micro-barometers antenna that provide complementary observations of source-radiated waves in heterogeneous and complex geophysical media. Exploiting these observations requires accurate and multi-physics – elastic, hydro-acoustic, infrasonic - wave simulation methods.

A popular approach is the Spectral Element Method (SEM) (Chaljub et al, 2006) which is high-order accurate (low dispersion error), very flexible to parallelization and computationally attractive due to efficient sum factorization technique and diagonal mass matrix. However SEMs suffer from lack of flexibility in handling complex geometry and multi-physics wave propagation. High-order Discontinuous Galerkin Methods (DGMs), i.e. Dumbser et al (2006), Etienne et al. (2010), Wilcox et al (2010), are recent alternatives that can handle complex geometry, space-and-time adaptativity, and allow efficient multi-physics wave coupling at interfaces. However, DGMs are more memory demanding and less computationally attractive than SEMs, especially when explicit time stepping is used.

We propose a new class of higher-order Hybridized Discontinuous Galerkin Spectral Elements (HDGSEM) methods for spatial discretization of wave equations, following the unifying framework for hybridization of Cockburn et al (2009) and Nguyen et al (2011), which allows for a single implementation of conforming and non-conforming SEMs. When used with energy conserving explicit time integration schemes, HDGSEM is flexible to handle complex geometry, computationally attractive and has significantly less degrees of freedom than classical DGMs, i.e., the only coupled unknowns are the single-valued numerical traces of the velocity field on the element's faces. The formulation can be extended to model fractional energy loss at interfaces between elastic, acoustic and hydro-acoustic media.

Accuracy and performance of the HDGSEM are illustrated and assessed agaings classical SEMs and DGMs through two-dimensional geophysical examples, as well as the coupling between hybridized SEMs and DGMs in the same computational domain. Finally on-going extension of this HDG method for three-dimensional wave propagation phenomena will be outlined in conclusion.