EDP Sciences Journals List
Advanced Search
Subscriber Authentication Point
Home All issues Volume 27 Article
Free access article

Issue ESAIM: Proc.
Volume 27, 2009
CANUM 2008
Page(s) 70 - 89
DOI http://dx.doi.org/10.1051/proc/2009020
Published online 25 June 2009

ESAIM: Proc., 2009, Vol. 27, pp. 70-89
DOI: 10.1051/proc/2009020

A high-order Discontinuous Galerkin method for the seismic wave propagation

Sarah Delcourte1, Loula Fezoui1 and Nathalie Glinsky-Olivier2

1  INRIA, NACHOS project-team, 2004 route des Lucioles, BP 93, 06902 Sophia Antipolis.
2  CERMICS and INRIA, NACHOS project-team, 2004 route des Lucioles, BP 93, 06902 Sophia Antipolis.

Sarah.Delcourte@inria.fr
Loula.Fezoui@inria.fr
Nathalie.Glinsky@inria.fr

Published online: 25 June 2009

Abstract
We are interested in the simulation of P-SV seismic wave propagation by a high-order Discontinuous Galerkin method based on centered fluxes at the interfaces combined with a leap-frog time-integration. This non-diffusive method, previously developed for the Maxwell equations [4, 9, 20], is particularly well adapted to complex topographies and fault discontinuities in the medium. We prove that the scheme is stable under a CFL type condition and that a discrete energy is preserved on an infinite domain. Convergence properties and efficiency of the method are studied through numerical simulations in two and three dimensions of space.


Résumé
Nous nous intéressons à la propagation d'ondes sismiques de types P et SV par une méthode de Galerkin Discontinue d'ordre élevé basée sur des flux centrés aux interfaces combinés à un schéma saute-mouton en temps. Cette méthode non-dissipative, précédemment développée pour les équations de Maxwell [4, 9, 20], est particulièrement bien adaptée à des milieux présentant une topographie complexe ou contenant des failles. On prouve la stabilité du schéma sous une condition de type CFL ainsi que la conservation d'une energie discrète dans un domaine infini. L'efficacité de la méthode est illustrée par des simulations numériques en deux et trois dimensions d'espace.



© EDP Sciences, ESAIM 2009


What is OpenURL?

The OpenURL standard is a protocol for transmission of metadata describing the resource that you wish to access. An OpenURL link contains article metadata and directs it to the OpenURL server of your choice. The OpenURL server can provide access to the resource and also offer complementary services (specific search engine, export of references...). The OpenURL link can be generated by different means.
  • If your librarian has set up your subscription with an OpenURL resolver, OpenURL links appear automatically on the abstract pages.
  • You can define your own OpenURL resolver with your EDPS Account. In this case your choice will be given priority over that of your library.
  • You can use an add-on for your browser (Firefox or I.E.) to display OpenURL links on a page (see http://www.openly.com/openurlref/). You should disable this module if you wish to use the OpenURL server that you or your library have defined.