Space-Time Trefftz - Discontinuous Galerkin Approximation for Elasto-Acoustics

Abstract : Wave reflection imaging for complex media can be effectively done by using advanced numerical methods (see [19]). In the context of the collaborative research program Depth Imaging Partnership (DIP) between Inria and Total, team-project Magique-3D and Prospective Lab of Houston develop high-order numerical schemes based mostly on discontinuous finite element approximation of wave fields. This technique, known as Discontinuous Galerkin (DG) method, is preferred because it takes into account geometrical and physical features of environment, and it is well-adapted for parallel computation [4, 19]. Recently it has been implemented for coupled elasto-acoustic problems, which led to the development of new propagators in time and frequency domains [6, 34]. However, when comparing to the conventional methods based on continuous approximation, the number of degrees of freedom required by DG method to achieve the same accuracy is significantly higher. To avoid this difficulty, Hybridizable Discontinuous Galerkin (HDG) methods have been developed and their integration into DIP is under way for both acoustic and elastic domains, with possibility of numerical coupling (see [22] and references therein). Another idea to explore consists in using Trefftz approximation space, whose elements are themselves discrete local solutions of the Acoustic System (AS) and Elastodynamic System (ES) [18, 32]. By its construction, Trefftz method reduces degrees of freedom, since it requires computing the surface integrals only to build variational formulation. Thus, we may consider the following advantages of Trefftz method compared to the standard ones: better order of convergence, flexibility in the choice of basis functions, low dispersion, incorporation of wave propagation directions in the discrete space, adaptivity and local space-time mesh refinement [18, 32]. Trefftz type methods have been widely used with time-harmonic formulations by Farhat, Tezaur, Harari, Hetmaniuk (2003 - 2006) (see [15, 31]), Gabard (2007) (see [17]), Badics (2014) (see [3]), Hiptmair, Moiola, Perugia (2011 - 2013) (see [20, 21, 28]) and others, while studies are still limited for reproducing temporal phenomena. Only few papers are interested in Maxwell equations in time [14, 23, 24, 30], but they are mostly devoted to a theoretical analysis of the method, showing the convergence and stability, and numerical tests with plane waves approximation are restricted to 1D + time dimensional case. Space-time Trefftz approximation by Lagrange multipliers for the second order formulation of the transient wave equation was explored in [5, 33]. Trefftz - Discontinuous Galerkin (Trefftz-DG) method for the first-order transient acoustic wave equations in arbitrary space dimensions extending the one-dimensional scheme of Kretzschmar et al. [23] has been introduced in recent paper of Moiola and Perugia (2017) [29]. the authors propose a complete a priori error analysis in both mesh-dependent and mesh-independent norms. In order to move on numerical simulations of geophysical phenomena and to consider more realistic application, DIP aims to develop new approximation techniques, which retain DG based methods but operates the Trefftz approach.
Complete list of metadatas

Cited literature [33 references]  Display  Hide  Download
Contributor : Elvira Shishenina <>
Submitted on : Tuesday, October 10, 2017 - 3:09:59 PM
Last modification on : Friday, June 7, 2019 - 3:16:15 PM
Long-term archiving on : Thursday, January 11, 2018 - 2:28:05 PM


Files produced by the author(s)


  • HAL Id : hal-01614126, version 1


Hélène Barucq, Henri Calandra, Julien Diaz, Elvira Shishenina. Space-Time Trefftz - Discontinuous Galerkin Approximation for Elasto-Acoustics. [Research Report] RR-9104, Inria Bordeaux Sud-Ouest; UPPA (LMA-Pau); Total E&P. 2017. ⟨hal-01614126⟩



Record views


Files downloads