Eléments finis mixtes spectraux et couches absorbantes parfaitement adaptées pour la propagation d'ondes élastiques en régime transitoire

Sandrine Fauqueux

tel-00007445, version 1

Eléments finis mixtes spectraux et couches absorbantes parfaitement adaptées pour la propagation d'ondes élastiques en régime transitoire

Sandrine Fauqueux ^1, 2, 3

ENSTA ParisTech (06/02/2003), COHEN Gary (Dir.)

1 : Département Géophysique à l'IFP (Institut Français du Pétrole)
Aucune France
2 : Unité de Mathématiques Appliquées (UMA)
ENSTA ParisTech France
3 : ONDES (INRIA Rocquencourt)
INRIA – CNRS : UMR2706 – ENSTA ParisTech France

Références bibliographiques

Type de publication : Thèses
titre : Eléments finis mixtes spectraux et couches absorbantes parfaitement adaptées pour la propagation d'ondes élastiques en régime transitoire
titre en anglais : Mixed spectral finite elements method and perfectly matched layers to model the propagation of elastic waves during time
date de soutenance : 06/02/2003
résumé : Nous nous intéressons à la propagation d'ondes en milieu élastique non-borné. Nous développons une nouvelle formulation mixte H(div)-L2 du système de l'élastodynamique linéaire et lui appliquons la "méthode des éléments finis mixtes spectraux". Cette nouvelle méthode permet, par un choix judicieux d'espaces d'approximation et une condensation de masse, d'obtenir un schéma explicite de stockage réduit, en donnant la même solution que la méthode des éléments finis spectraux. Nous introduisons ensuite des couches absorbantes parfaitement adaptées pour modéliser les milieux non-bornés. Des phénomènes d'instabilité sont révélés et analysés pour certains matériaux élastiques 2D. La méthode numérique obtenue est validée et testée sur des modèles réalistes en acoustique et élastique. Une analyse par ondes planes donne des résultats de dispersion numérique et montre la supériorité des maillages adaptés aux vitesses du milieu. Enfin, une extension au couplage fluide-structure 2D est mise en place.
résumé en anglais : We consider the propagation of elastic waves in unbounded domains. A new formulation of the linear elasticity system as an H(div)-L2 system enables us to use the "mixed spectral finite element method". This new method is based on the definition of new spaces of approximation and the use of mass-lumping. It leads to an explicit scheme with reduced storage and provides the same solution as the spectral finite element method. Then, we modelize unbounded domains by using Perfectly Matched Layers. Instabilities in the PML in the case of particular 2D elastic media are pointed out and investigated. The numerical method is validated and tested in the case of acoustic and elastic realistic models. A plane wave analysis gives results about numerical dispersion and shows that meshes adapted to the physical and geometrical properties of the media are more accurate than the others. Then, an extension of the method to fluid-solid coupling is introduced for 2D seismic propagation.
domaine :
Informatique/Modélisation et simulation

Mathématiques

Planète et Univers/Sciences de la Terre/Géophysique

Physique/Physique/Géophysique

Sciences de l'environnement/Milieux et Changements globaux
organisme de délivrance : ENSTA ParisTech
Numéro national (NNT) : 2003PA090002
langue : Français
directeur de thèse : COHEN Gary
courriel du directeur : Gary.Cohen@inria.fr
composition du Jury : SVAY-LUCAS Julie (co-directeur)
mots-clés : élastodynamique linéaire – éléments spectraux – condensation de masse – couche absorbante parfaitement adaptée – analyse par ondes planes
mots-clés en anglais : elasticity, spectral element, mass-lumping, perfectly matched layers, plane wave analysis
commentaire : Directeur de thèse: M. Gary COHEN Promoteur IFP: Mme Julie SVAY-LUCAS Président: M. Yvon MADAY Rapporteurs: Mme Laurence HALPERN M. Peter MONK Examinateurs: M. Patrick JOLY M. Jean-Pierre VILOTTE

Liste des fichiers attachés à ce document :

PDF

tel-00007445.pdf

(7.3 MB)

tel-00007445, version 1
http://tel.archives-ouvertes.fr/tel-00007445
oai:tel.archives-ouvertes.fr:tel-00007445
Contributeur : Sandrine Fauqueux <>
Soumis le : Jeudi 18 Novembre 2004, 16:24:53
Dernière modification le : Jeudi 10 Novembre 2011, 15:49:26

Voir la fiche concise

Documents associés

PDF :

Exporter

Bibtex EndNote TEI RefWorks

%% tel-00007445, version 1
%% http://tel.archives-ouvertes.fr/tel-00007445
@phdthesis{fauqueux:tel-00007445,
    hal_id = {tel-00007445},
    url = {http://tel.archives-ouvertes.fr/tel-00007445},
    title = {{El{\'e}ments finis mixtes spectraux et couches absorbantes parfaitement adapt{\'e}es pour la propagation d'ondes {\'e}lastiques en r{\'e}gime transitoire}},
    author = {Fauqueux, Sandrine},
    abstract = {{Nous nous int{\'e}ressons {\`a} la propagation d'ondes en milieu {\'e}lastique non-born{\'e}. Nous d{\'e}veloppons une nouvelle formulation mixte H(div)-L2 du syst{\`e}me de l'{\'e}lastodynamique lin{\'e}aire et lui appliquons la "m{\'e}thode des {\'e}l{\'e}ments finis mixtes spectraux". Cette nouvelle m{\'e}thode permet, par un choix judicieux d'espaces d'approximation et une condensation de masse, d'obtenir un sch{\'e}ma explicite de stockage r{\'e}duit, en donnant la m{\^e}me solution que la m{\'e}thode des {\'e}l{\'e}ments finis spectraux. Nous introduisons ensuite des couches absorbantes parfaitement adapt{\'e}es pour mod{\'e}liser les milieux non-born{\'e}s. Des ph{\'e}nom{\`e}nes d'instabilit{\'e} sont r{\'e}v{\'e}l{\'e}s et analys{\'e}s pour certains mat{\'e}riaux {\'e}lastiques 2D. La m{\'e}thode num{\'e}rique obtenue est valid{\'e}e et test{\'e}e sur des mod{\`e}les r{\'e}alistes en acoustique et {\'e}lastique. Une analyse par ondes planes donne des r{\'e}sultats de dispersion num{\'e}rique et montre la sup{\'e}riorit{\'e} des maillages adapt{\'e}s aux vitesses du milieu. Enfin, une extension au couplage fluide-structure 2D est mise en place.}},
    keywords = {{\'e}lastodynamique lin{\'e}aire, {\'e}l{\'e}ments spectraux, condensation de masse, couche absorbante parfaitement adapt{\'e}e, analyse par ondes planes},
    language = {Fran{\c c}ais},
    affiliation = {D{\'e}partement G{\'e}ophysique {\`a} l'IFP (Institut Fran{\c c}ais du P{\'e}trole) , Unit{\'e} de Math{\'e}matiques Appliqu{\'e}es - UMA , ONDES - INRIA Rocquencourt},
    school = {ENSTA ParisTech},
    type = {THESE},
    note = {Directeur de th{\`e}se: M. Gary COHEN Promoteur IFP: Mme Julie SVAY-LUCAS Pr{\'e}sident: M. Yvon MADAY Rapporteurs: Mme Laurence HALPERN M. Peter MONK Examinateurs: M. Patrick JOLY M. Jean-Pierre VILOTTE },
    year = {2003},
    month = Feb,
    pdf = {http://tel.archives-ouvertes.fr/tel-00007445/PDF/tel-00007445.pdf},
}

%F tel-00007445, version 1
%0 Thesis
%U http://tel.archives-ouvertes.fr/tel-00007445
%T Eléments finis mixtes spectraux et couches absorbantes parfaitement adaptées pour la propagation d'ondes élastiques en régime transitoire
%A Fauqueux, Sandrine
%+ Département Géophysique à l'IFP (Institut Français du Pétrole) , Unité de Mathématiques Appliquées - UMA , ONDES - INRIA Rocquencourt
%T Mixed spectral finite elements method and perfectly matched layers to model the propagation of elastic waves during time
%8 2003-02-06
%X Nous nous intéressons à la propagation d'ondes en milieu élastique non-borné. Nous développons une nouvelle formulation mixte H(div)-L2 du système de l'élastodynamique linéaire et lui appliquons la "méthode des éléments finis mixtes spectraux". Cette nouvelle méthode permet, par un choix judicieux d'espaces d'approximation et une condensation de masse, d'obtenir un schéma explicite de stockage réduit, en donnant la même solution que la méthode des éléments finis spectraux. Nous introduisons ensuite des couches absorbantes parfaitement adaptées pour modéliser les milieux non-bornés. Des phénomènes d'instabilité sont révélés et analysés pour certains matériaux élastiques 2D. La méthode numérique obtenue est validée et testée sur des modèles réalistes en acoustique et élastique. Une analyse par ondes planes donne des résultats de dispersion numérique et montre la supériorité des maillages adaptés aux vitesses du milieu. Enfin, une extension au couplage fluide-structure 2D est mise en place.
%X We consider the propagation of elastic waves in unbounded domains. A new formulation of the linear elasticity system as an H(div)-L2 system enables us to use the "mixed spectral finite element method". This new method is based on the definition of new spaces of approximation and the use of mass-lumping. It leads to an explicit scheme with reduced storage and provides the same solution as the spectral finite element method. Then, we modelize unbounded domains by using Perfectly Matched Layers. Instabilities in the PML in the case of particular 2D elastic media are pointed out and investigated. The numerical method is validated and tested in the case of acoustic and elastic realistic models. A plane wave analysis gives results about numerical dispersion and shows that meshes adapted to the physical and geometrical properties of the media are more accurate than the others. Then, an extension of the method to fluid-solid coupling is introduced for 2D seismic propagation.
%2 INFO:INFO_MO;MATH;SDU:STU:GP;PHYS:PHYS:PHYS_GEO-PH;SDE:MCG
%I ENSTA ParisTech
%2 2003PA090002
%G Français
%Z COHEN Gary
%Z Gary.Cohen@inria.fr
%Z SVAY-LUCAS Julie (co-directeur)
%K élastodynamique linéaire, éléments spectraux, condensation de masse, couche absorbante parfaitement adaptée, analyse par ondes planes
%K elasticity, spectral element, mass-lumping, perfectly matched layers, plane wave analysis
%Z Directeur de thèse: M. Gary COHEN Promoteur IFP: Mme Julie SVAY-LUCAS Président: M. Yvon MADAY Rapporteurs: Mme Laurence HALPERN M. Peter MONK Examinateurs: M. Patrick JOLY M. Jean-Pierre VILOTTE