Space-Time Trefftz - Discontinuous Galerkin Approximation for Elasto-Acoustics - Inria - Institut national de recherche en sciences et technologies du numérique Accéder directement au contenu
Rapport (Rapport De Recherche) Année : 2017

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

Approximation Espace-Temps de Trefftz - Galerkin Discontinue pour l’Élasto-Acoustique

Résumé

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, thenumber of degrees of freedom required by DG method to achieve the same accuracy is significantlyhigher. To avoid this difficulty, Hybridizable Discontinuous Galerkin (HDG) methods have been developedand their integration into DIP is under way for both acoustic and elastic domains, with possibilityof 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 thediscrete 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 limitedfor 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 theconvergence and stability, and numerical tests with plane waves approximation are restricted to1D + time dimensional case.Space-time Trefftz approximation by Lagrange multipliers for the second order formulation of thetransient wave equation was explored in [5, 33].Trefftz - Discontinuous Galerkin (Trefftz-DG) method for the first-order transient acoustic waveequations in arbitrary space dimensions extending the one-dimensional scheme of Kretzschmar etal. [23] has been introduced in recent paper of Moiola and Perugia (2017) [29]. the authors proposea 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 realisticapplication, DIP aims to develop new approximation techniques, which retain DG based methodsbut operates the Trefftz approach.
L’imagerie de milieux complexes par réflexion d’ondes peut se faire de manière efficace en utilisant des méthodes numériques avancées. Dans le contexte du programme de recherche collaborative DIP liant Inria et Total, l’équipe Magique-3D développe avec le Prospective Lab de Houston des schémas d’ordre élevé, essentiellement basés sur des approximations discontinues par éléments finis des champs d’ondes. Cette technique, connue sous le nom de méthode DG est privilégiée car elle permet de tenir compte des caractéristiques géometriques et physiques du milieu et qu’elle est adaptée au calcul parallèle [4, 19]. Récemment, elle a été mise en oeuvre pour des problèmes couplés d’élasto-acoustique, ce qui a donné lieu au développement des nouveax propagateurs dans le domaine temporel et le domaine harmonique [6, 34].Si les méthodes d’éléments finis discontinus ont fait leurs preuves en matière de précision et de flexibilité, elles sont critiquées pour le nombre de degrés de liberté qu’elles utilisent qui s’avère beaucoup plus élevé que les méthodes classiques basées sur des approximations continues.Pour surmonter cette difficulté, des méthodes hybrides dites "Hybridizable Discontinuous Galerkin methods" ont été développées et leur intégration dans le projet DIP est en cours, aussi bien dans le domaine acoustique qu’élastique (voir [22]).Une autre idée consiste à utiliser des méthodes de Trefftz pour lesquelles on utilise des espacesdiscrets dont les éléments sont des solutions locales des équations à résoudre.Cette approche a été souvent utilisée pour résoudre des équations d’ondes harmoniques parFarhat, Tezaur, Harari, Hetmaniuk (2003 - 2006) (voir [15, 31]), Gabard (2007) (voir [17]),Badics (2014) (voir [3]), Hiptmair, Moiola, Perugia (2011 - 2013) (voir [20, 21, 28]) et les autres,mais très peu pour reproduire des phénomènes temporels. Quelques articles récents de 2014-2017[14, 23, 24, 29, 30] s’intéressent aux équations temporelles de Maxwell mais l’essentiel des papiersse consacre à l’analyse théorique de la méthode, démontrant la convergence et stabilité.Dans le but de continuer à faire évoluer son expertise dans le domaine de la simulationnumérique de phénomènes géophysiques, le projet DIP souhaite développer de nouvelles tech-niques d’approximation qui conservent pour socle les méthodes DG mais qui exploitent de Trefftzapproche dans le but de considérer des applications de plus en plus réalistes.
Fichier principal
Vignette du fichier
RR_9104.pdf (13.63 Mo) Télécharger le fichier
Origine : Fichiers produits par l'(les) auteur(s)
Loading...

Dates et versions

hal-01614126 , version 1 (10-10-2017)

Identifiants

  • HAL Id : hal-01614126 , version 1

Citer

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⟩
304 Consultations
193 Téléchargements

Partager

Gmail Facebook X LinkedIn More