Skip to Main content Skip to Navigation

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

Abstract : 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.
Document type :
Complete list of metadatas

Cited literature [69 references]  Display  Hide  Download
Contributor : Abes Star :  Contact
Submitted on : Thursday, October 3, 2019 - 1:57:07 PM
Last modification on : Wednesday, October 14, 2020 - 3:46:34 AM


Version validated by the jury (STAR)


  • HAL Id : tel-02304642, version 1



Jérôme Tomezyk. Numerical resolution of some Helmholtz-type problems with impedance boundary condition or PML. Numerical Analysis [math.NA]. Université de Valenciennes et du Hainaut-Cambresis, 2019. English. ⟨NNT : 2019VALE0017⟩. ⟨tel-02304642⟩



Record views


Files downloads