Approximation par éléments finis de problèmes d'Helmholtz pour la propagation d'ondes sismiques / Théophile Chaumont Frelet ; sous la direction de Christian Gout et de Hélène Barucq

Date :

Type : Livre / Book

Type : Thèse / Thesis

Langue / Language : anglais / English

Ondes -- Propagation

Helmholtz, Équation d'

Éléments finis, Méthode des

Analyse multiéchelle

Pollution -- Effets physiologiques

Méthode sismique-réflexion

Gout, Christian (1969-.... ; chercheur en mathématiques appliquées) (Directeur de thèse / thesis advisor)

Barucq, Hélène (19..-.... ; chercheuse en mathématiques) (Directeur de thèse / thesis advisor)

Nicaise, Serge (1960-....) (Président du jury de soutenance / praeses)

Allaire, Grégoire (1963-....) (Rapporteur de la thèse / thesis reporter)

Abgrall, Rémi (1961-....) (Rapporteur de la thèse / thesis reporter)

Melenk, Jens Markus (1967-.... ; mathématicien) (Rapporteur de la thèse / thesis reporter)

Le Guyader, Carole (1979-.... ; chercheuse en mathématiques appliquées) (Membre du jury / opponent)

Calandra, Henri (1960-.... ; ingénieur de recherche en méthodes numériques) (Membre du jury / opponent)

Institut national des sciences appliquées Rouen Normandie (Saint-Etienne-du-Rouvray ; 1985-....) (Organisme de soutenance / degree-grantor)

École doctorale sciences physiques mathématiques et de l'information pour l'ingénieur (Saint-Etienne-du-Rouvray, Seine-Maritime ; ....-2016) (Ecole doctorale associée à la thèse / doctoral school)

Laboratoire de mathématiques de l'INSA Rouen (Saint Etienne du Rouvray (Seine-Maritime) ; 1987-....) (Laboratoire associé à la thèse / thesis associated laboratory)

Laboratoire de mathématiques et de leurs applications (Pau) (Laboratoire associé à la thèse / thesis associated laboratory)

Institut national de recherche en informatique et en automatique (France). Centre de recherche Bordeaux - Sud-Ouest (Laboratoire associé à la thèse / thesis associated laboratory)

Magique 3D. Modélisation Avancée en GéophysIQUE 3D (Pau) (Equipe de recherche associée à la thèse / thesis associated research team)

Résumé / Abstract : Dans cette thèse, on s'intéresse à la propagation d'ondes en milieu fortement hétérogène modélisée par l'équation d'Helmholtz. Les méthodes numériques permettant de résoudre ce problème souffrent de dispersion numérique, en particulier à haute fréquence. Ce phénomène, appelé "effet de pollution", est largement analysé dans la littérature quand le milieu de propagation est homogène et l'utilisation de "méthodes d'ordre élevé" est souvent proposée pour minimiser ce problème. Dans ce travail, on s'intéresse à un milieu de propagation hétérogène, cas pour lequel on dispose de moins de connaissances. On propose d'adapter des méthodes éléments finis d'ordre élevé pour résoudre l'équation d'Helmholtz en milieu hétérogène, afin de réduire l'effet de pollution. Les méthodes d'ordre élevé étant généralement basées sur des maillages "larges", une stratégie multi-échelle originale est développée afin de prendre en compte des hétérogénéités de petite échelle. La convergence de la méthode est démontrée. En particulier, on montre que la méthode est robuste vis-a-vis de l'effet de pollution. D'autre part, on applique la méthode a plusieurs cas-tests numériques. On s'intéresse d'abord à des problèmes académiques, qui permettent de valider la théorie de convergence développée. On considère ensuite des cas-tests "industriels" appliqués à la Géophysique. Ces derniers nous permettent de conclure que la méthode multi-échelle proposée est plus performante que les éléments finis "classiques" et que des problèmes 3D réalistes peuvent être considérés.

Résumé / Abstract : The main objective of this work is the design of an efficient numerical strategy to solve the Helmholtz equation in highly heterogeneous media. We propose a methodology based on coarse meshes and high order polynomials together with a special quadrature scheme to take into account fine scale heterogeneities. The idea behind this choice is that high order polynomials are known to be robust with respect to the pollution effect and therefore, efficient to solve wave problems in homogeneous media. In this work, we are able to extend so-called "asymptotic error-estimate" derived for problems homogeneous media to the case of heterogeneous media. These results are of particular interest because they show that high order polynomials bring more robustness with respect to the pollution effect even if the solution is not regular, because of the fine scale heterogeneities. We propose special quadrature schemes to take int account fine scale heterogeneities. These schemes can also be seen as an approximation of the medium parameters. If we denote by h the finite-element mesh step and by e the approximation level of the medium parameters, we are able to show a convergence theorem which is explicit in terms of h, e and f, where f is the frequency. The main theoretical results are further validated through numerical experiments. 2D and 3D geophysica benchmarks have been considered. First, these experiments confirm that high-order finite-elements are more efficient to approximate the solution if they are coupled with our multiscale strategy. This is in agreement with our results about the pollution effect. Furthermore, we have carried out benchmarks in terms of computational time and memory requirements for 3D problems. We conclude that our multiscale methodology is able to greatly reduce the computational burden compared to the standard finite-element method