Résolution numérique de quelques problèmes du type Helmholtz avec conditions au bord d'impédance ou des couches absorbantes (PML)

Numerical resolution of some Helmholtz-type problems with impedance boundary condition or PML

• Auteur / Autrice : Jérôme Tomezyk
• Direction : Serge Nicaise
• Type : Thèse de doctorat
• Discipline(s) : Mathématiques. Mathématiques appliquées
• Date : Soutenance le 02/07/2019
• Etablissement(s) : Valenciennes
• Ecole(s) doctorale(s) : École doctorale Sciences pour l'ingénieur (Lille)
• Partenaire(s) de recherche : Laboratoire : Laboratoire de mathématiques et leurs applications de Valenciennes (2006-2021)
• Partenaire(s) de recherche : Communauté d'Universités et Etablissements (ComUE) : Communauté d'universités et d'établissements Lille Nord de France (2009-2013)
• Jury : Président / Présidente : Emmanuel Creusé
• Jury : Examinateurs / Examinatrices : Serge Nicaise, Victor Péron, Stefan A. Sauter, Marion Darbas
• Jury : Rapporteurs / Rapporteuses : Victor Péron, Stefan A. Sauter

Mots clés

Mots clés contrôlés

Éléments finis, Méthode des

Maxwell, Équations de

Helmholtz, Équation d'

Mots clés libres

Effet de pollution

Condition au bord d'impédance

Finite element method

Helmholtz equation

Pollution effect

Maxwell's equations

Impedance boundary condition

Perfectly Matched Layer (PML)

Résumé

Dans cette thèse, nous étudions la convergence de méthode de type éléments finis pour les équations de Maxwell en régime harmonique avec condition au bord d'impédance et l'équation de Helmholtz avec une couche parfaitement absorbante(PML). On étudie en premier, la formulation régularisée de l'équation de Maxwell en régime harmonique avec condition au bord d'impédance (qui consiste à ajouter le term ∇ div à l'équation originale pour avoir un problème elliptique) et on garde la condition d'impédance comme une condition au bord essentielle. Pour des domaines à bord régulier, le caractère bien posé de cette formulation est bien connu mais cela n'est pas le cas pour des domaines polyédraux convexes. On commence alors le premier chapitre par la preuve du caractère bien posé dans le cas du polyèdre convexe, qui est basé sur le fait que l'espace variationnel est inclus dans H¹. Dans le but d'avoir des estimations explicites en le nombre d'onde k de ce problème, il est obligatoire d'avoir des résultats de stabilité explicites en ce nombre d'onde. C'est aussi proposé, pour quelques situations particulières, dans ce chapitre. Dans le second chapitre on décrit les singularités d'arêtes et de coins pour notre problème. On peut alors déduire la régularité de la solution du problème original, ainsi que de son adjoint. On a tous les ingrédients pour proposer une analyse de convergence explicite en k pour une méthode d'éléments finis avec éléments de Lagrange. Dans le troisième chapitre, on considère une méthode d'éléments finis hp non conforme pour un domaine à bord régulier. Pour obtenir des estimations explicites en k, on introduit un résultat de décomposition, qui sépare la solution du problème original (ou de son adjoint) en une partie régulière mais fortement oscillante et une partie moins régulière mais peu oscillante. Ce résultat permet de montrer des estimations explicites en k. Le dernier chapitre est dédié à l'équation de Helmholtz avec une PML. L'équation de Helmholtz dans l'espace entier est souvent utilisée pour modéliser la diffraction d'onde acoustique (en régime harmonique), avec la condition de radiation à l'infini de Sommerfeld. L'ajout d'une PML est une façon pour passer d'un domaine infini à un domaine fini, elle correspond à l'ajout d'une couche autour du domaine de calcul qui absorbe très vite toutes les ondes sortantes. On propose en premier un résultat de stabilité explicite en k. On propose alors deux schémas numériques, une méthode d'éléments finis hp et une méthode multi- échelle basée sur un sous-espace local de correction. Le résultat de stabilité est utilisé pour mettre en relation de choix des paramètres des méthodes numériques considérées avec k. Nous montrons aussi des estimations d'erreur a priori. A la fin de ces chapitres, des tests numériques sont proposés pour confirmer nos résultats théoriques.

In this thesis, we propose wavenumber explicit convergence analyses of some finite element methods for time-harmonic Maxwell's equations with impedance boundary condition and for the Helmholtz equation with Perfectly Matched Layer (PML). We first study the regularized formulation of time-harmonic Maxwell's equations with impedance boundary conditions (where we add a ∇ div-term to the original equation to have an elliptic problem) and keep the impedance boundary condition as an essential boundary condition. For a smooth domain, the wellposedness of this formulation is well-known. But the well-posedness for convex polyhedral domain has been not yet investigated. Hence, we start the first chapter with the proof of the well-posedness in this case, which is based on the fact that the variational space is embedded in H¹. In order to perform a wavenumber explicit error analysis of our problem, a wavenumber explicit stability estimate is mandatory. We then prove such an estimate for some particular configurations. In the second chapter, we describe the corner and edge singularities for such problem. Then we deduce the regularity of the solution of the original and the adjoint problem, thus we have all ingredients to propose a explicit wavenumber convergence analysis for h-FEM with Lagrange element. In the third chapter, we consider a non conforming hp-finite element approximation for domains with a smooth boundary. To perform a wavenumber explicit error analysis, we split the solution of the original problem (or its adjoint) into a regular but oscillating part and a rough component that behaves nicely for large frequencies. This result allows to prove convergence analysis for our FEM, again explicit in the wavenumber. The last chapter is dedicated to the Helmholtz equation with PML. The Helmholtz equation in full space is often used to model time harmonic acoustic scattering problems, with Sommerfeld radiation condition at infinity. Adding a PML is a way to reduce the infinite domain to a finite one. It corresponds to add an artificial absorbing layer surrounding a computational domain, in which scattered wave will decrease very quickly. We first propose a wavenumber explicit stability result for such problem. Then, we propose two numerical discretizations: an hp-FEM and a multiscale method based on local subspace correction. The stability result is used to relate the choice of the parameters in the numerical methods to the wavenumber. A priori error estimates are shown. At the end of each chapter, we perform numerical tests to confirm our theoritical results.