An hp-adaptive discontinuous Galerkin finite-element method for 3-D elastic wave modelling - Université Pierre et Marie Curie Accéder directement au contenu
Article Dans Une Revue Geophysical Journal International Année : 2010

An hp-adaptive discontinuous Galerkin finite-element method for 3-D elastic wave modelling

Résumé

We present a discontinuous Galerkin finite-element method (DG-FEM) formulation with Convolutional Perfectly Matched Layer (CPML) absorbing boundary condition for 3-D elastic seismic wave modelling. This method makes use of unstructured tetrahedral meshes locally refined according to the medium properties (h-adaptivity), and of approximation orders that can change from one element to another according to an adequate criterion (p-adaptivity). These two features allow us to significantly reduce the computational cost of the simulations. Moreover, we have designed an efficient CPML absorbing boundary condition, both in terms of absorption and computational cost, by combining approximation orders in the numerical domain. A quadratic interpolation is typically used in the medium to obtain the required accuracy, while lower approximation orders are used in the CPMLs to reduce the total computational cost and to obtain a well-balanced workload over the processors. While the efficiency of DG-FEMs have been largely demonstrated for high approximation orders, we favour the use of low approximation orders as they are more appropriate to the applications we are interested in. In particular, we address the issues of seismic modelling and seismic imaging in cases of complex geological structures that require a fine discretization of the medium. We illustrate the efficiency of our approach within the framework of the EUROSEISTEST verification and validation project, which is designed to compare high-frequency (up to 4 Hz) numerical predictions of ground motion in the Volvi basin (Greece). Through the tetrahedral meshing, we have achieved fine discretization of the basin, which appears to be a sine qua non condition for accurate computation of surface waves diffracted at the basin edges. We compare our results with predictions computed with the spectral element method (SEM), and demonstrate that our method yields the same level of accuracy with computation times of the same order of magnitude.
Fichier principal
Vignette du fichier
183-2-941.pdf (1.44 Mo) Télécharger le fichier
Origine : Fichiers éditeurs autorisés sur une archive ouverte

Dates et versions

insu-00565022 , version 1 (03-03-2021)

Identifiants

Citer

V. Etienne, E. Chaljub, J. Virieux, N. Glinsky. An hp-adaptive discontinuous Galerkin finite-element method for 3-D elastic wave modelling. Geophysical Journal International, 2010, 183, pp.941-962. ⟨10.1111/J.1365-246X.2010.04764.X⟩. ⟨insu-00565022⟩
288 Consultations
63 Téléchargements

Altmetric

Partager

Gmail Facebook X LinkedIn More