| Title: |
An hp-adaptive discontinuous Galerkin finite-element method for 3-D elastic wave modelling |
| Authors: |
Etienne, V.; Chaljub, E.; Virieux, J.; Glinsky, N. |
| Contributors: |
Géoazur (GEOAZUR 6526); Institut de Recherche pour le Développement (IRD)-Université Pierre et Marie Curie - Paris 6 (UPMC)-Université Nice Sophia Antipolis (1965 - 2019) (UNS)-Institut national des sciences de l'Univers (INSU - CNRS)-Observatoire de la Côte d'Azur; Université Côte d'Azur (UniCA)-Université Côte d'Azur (UniCA)-Centre National de la Recherche Scientifique (CNRS); Laboratoire de Géophysique Interne et Tectonophysique (LGIT); Observatoire des Sciences de l'Univers de Grenoble (OSUG); Université Joseph Fourier - Grenoble 1 (UJF)-Institut polytechnique de Grenoble - Grenoble Institute of Technology (Grenoble INP)-Institut national des sciences de l'Univers (INSU - CNRS)-Institut national de recherche en sciences et technologies pour l'environnement et l'agriculture (IRSTEA)-Université Savoie Mont Blanc (USMB Université de Savoie Université de Chambéry )-Centre National de la Recherche Scientifique (CNRS)-Université Joseph Fourier - Grenoble 1 (UJF)-Institut polytechnique de Grenoble - Grenoble Institute of Technology (Grenoble INP)-Institut national des sciences de l'Univers (INSU - CNRS)-Institut national de recherche en sciences et technologies pour l'environnement et l'agriculture (IRSTEA)-Université Savoie Mont Blanc (USMB Université de Savoie Université de Chambéry )-Centre National de la Recherche Scientifique (CNRS)-Laboratoire Central des Ponts et Chaussées (LCPC)-Institut de Recherche pour le Développement (IRD)-Centre National de la Recherche Scientifique (CNRS); Centre Inria d'Université Côte d'Azur; Institut National de Recherche en Informatique et en Automatique (Inria); ERA 6 Risque sismique (ERA 6 Risque sismique - Equipe recherche associée au LCPC); Avant création Cerema; BP, CGG-VERITAS, EXXON MOBIL, SHELL and TOTAL; CEA/CASHIMA |
| Source: |
ISSN: 0956-540X. |
| Publisher Information: |
CCSD; Oxford University Press (OUP) |
| Publication Year: |
2010 |
| Collection: |
Université Savoie Mont Blanc: HAL |
| Subject Terms: |
Surface waves and free oscillations; Computational seismology; Site effects; Wave propagation; [SDU.STU.GP]Sciences of the Universe [physics]/Earth Sciences/Geophysics [physics.geo-ph]; [PHYS.PHYS.PHYS-GEO-PH]Physics [physics]/Physics [physics]/Geophysics [physics.geo-ph]; [SDE.MCG]Environmental Sciences/Global Changes |
| Description: |
International audience ; 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. |
| Document Type: |
article in journal/newspaper |
| Language: |
English |
| Relation: |
BIBCODE: 2010GeoJI.183.941E |
| DOI: |
10.1111/J.1365-246X.2010.04764.X |
| Availability: |
https://insu.hal.science/insu-00565022; https://insu.hal.science/insu-00565022v1/document; https://insu.hal.science/insu-00565022v1/file/183-2-941.pdf; https://doi.org/10.1111/J.1365-246X.2010.04764.X |
| Rights: |
info:eu-repo/semantics/OpenAccess |
| Accession Number: |
edsbas.F0C12FE |
| Database: |
BASE |