Analysis of a high-order space and time discontinuous Galerkin method for elastodynamic equations. Application to 3D wave propagation
- Creators
- Delcourte, Sarah
- Glinsky, Nathalie
- Others:
- Université Claude Bernard Lyon 1 (UCBL) ; Université de Lyon
- Kinetic models AppLIed for Future of Fusion Energy (KALiFFE) ; Inria Grenoble - Rhône-Alpes ; Institut National de Recherche en Informatique et en Automatique (Inria)-Institut National de Recherche en Informatique et en Automatique (Inria)-Institut Camille Jordan (ICJ) ; École Centrale de Lyon (ECL) ; Université de Lyon-Université de Lyon-Université Claude Bernard Lyon 1 (UCBL) ; Université de Lyon-Institut National des Sciences Appliquées de Lyon (INSA Lyon) ; Université de Lyon-Institut National des Sciences Appliquées (INSA)-Institut National des Sciences Appliquées (INSA)-Université Jean Monnet - Saint-Étienne (UJM)-Centre National de la Recherche Scientifique (CNRS)-École Centrale de Lyon (ECL) ; Université de Lyon-Université de Lyon-Université Claude Bernard Lyon 1 (UCBL) ; Université de Lyon-Institut National des Sciences Appliquées de Lyon (INSA Lyon) ; Université de Lyon-Institut National des Sciences Appliquées (INSA)-Institut National des Sciences Appliquées (INSA)-Université Jean Monnet - Saint-Étienne (UJM)-Centre National de la Recherche Scientifique (CNRS)-Institut National des Sciences Mathématiques et de leurs Interactions (INSMI)
- Modélisation mathématique, calcul scientifique (MMCS) ; Institut Camille Jordan (ICJ) ; École Centrale de Lyon (ECL) ; Université de Lyon-Université de Lyon-Université Claude Bernard Lyon 1 (UCBL) ; Université de Lyon-Institut National des Sciences Appliquées de Lyon (INSA Lyon) ; Université de Lyon-Institut National des Sciences Appliquées (INSA)-Institut National des Sciences Appliquées (INSA)-Université Jean Monnet - Saint-Étienne (UJM)-Centre National de la Recherche Scientifique (CNRS)-École Centrale de Lyon (ECL) ; Université de Lyon-Université de Lyon-Université Claude Bernard Lyon 1 (UCBL) ; Université de Lyon-Institut National des Sciences Appliquées de Lyon (INSA Lyon) ; Université de Lyon-Institut National des Sciences Appliquées (INSA)-Institut National des Sciences Appliquées (INSA)-Université Jean Monnet - Saint-Étienne (UJM)-Centre National de la Recherche Scientifique (CNRS)
- Numerical modeling and high performance computing for evolution problems in complex domains and heterogeneous media (NACHOS) ; Inria Sophia Antipolis - Méditerranée (CRISAM) ; Institut National de Recherche en Informatique et en Automatique (Inria)-Institut National de Recherche en Informatique et en Automatique (Inria)-Laboratoire Jean Alexandre Dieudonné (JAD) ; Université Nice Sophia Antipolis (1965 - 2019) (UNS) ; COMUE Université Côte d'Azur (2015-2019) (COMUE UCA)-COMUE Université Côte d'Azur (2015-2019) (COMUE UCA)-Centre National de la Recherche Scientifique (CNRS)-Université Côte d'Azur (UCA)-Université Nice Sophia Antipolis (1965 - 2019) (UNS) ; COMUE Université Côte d'Azur (2015-2019) (COMUE UCA)-COMUE Université Côte d'Azur (2015-2019) (COMUE UCA)-Centre National de la Recherche Scientifique (CNRS)-Université Côte d'Azur (UCA)
- Département Géotechnique, Eau et Risques (IFSTTAR/GER) ; Institut Français des Sciences et Technologies des Transports, de l'Aménagement et des Réseaux (IFSTTAR)-PRES Université Paris-Est
- Centre d'Etudes et d'Expertise sur les Risques, l'Environnement, la Mobilité et l'Aménagement - Equipe-projet MOUVGS (Cerema Equipe-projet MOUVGS) ; Centre d'Etudes et d'Expertise sur les Risques, l'Environnement, la Mobilité et l'Aménagement (Cerema)
Description
in this paper, we introduce a high-order discontinuous Galerkin method, based on centerd fluxes and a family of high-order leap-frog time schemes, for the solution of the 3D elastodynamic equations written in velocity-stress formulation. We prove that this explicit scheme is stable under a CFL type condition obtained from a discrete energy which is preserved in domains with free surface or decreasing in domains with absorbing boundary conditions. Moreover, we study the convergence of the method for both the semi-discrete and the fully discrete schemes and we illustrate the convergence results by the propagation of an eigenmode. We also propose a series of absorbing conditions which allow improving the convergence of the global scheme. Finally, several numerical applications of wave propagation, using a 3D solver, help illustrating the various properties of the method.
Abstract
Accepté pour publication
Abstract
International audience
Additional details
- URL
- https://hal.inria.fr/hal-01109424
- URN
- urn:oai:HAL:hal-01109424v1
- Origin repository
- UNICA