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.

Exact analytical solutions for the quasistationary electromagnetic fields in the presence of conducting objects require the field solutions both internal and external to the conductors. Such solutions are limited for certain canonically shaped objects but are useful in testing the accuracy of various approximate models and numerical methods developed to solve complex problems related to real world conducting objects and in calibrating instruments designed to measure various field quantities. Theoretical investigations of quasistationary electromagnetic fields also aid in improving the understanding of the physical phenomena of electromagnetic induction. This thesis presents rigorous analytical expressions derived as benchmark solutions for the quasistationary field quantities both inside and outside, Joule losses and the electromagnetic forces acting upon a conducting spheroid placed in the proximity of a non-uniform field produced by current-carrying turns. These expressions are used to generate numerous numerical results of specified accuracy and selected results are presented in a normalized form for extended ranges of the spheroid axial ratio, the ratio of the depth of penetration to the semi-minor axis and the position of the inducing turns relative to the spheroids. They are intended to constitute reference data to be employed for comprehensive comparisons of results from approximate numerical methods or from boundary impedance models used for real world conductors. Approximate boundary conditions such as the simpler perfect electric conductor model or the Leontovich surface impedance boundary condition model can be used to obtain approximate solutions by only analyzing the field external to the conducting object. The range of validity of these impedance boundary condition models for the analysis of axisymmetric eddy-current problems is thoroughly investigated. While the simpler PEC model can be employed only when the electromagnetic depth of penetration is much smaller than the smallest local radius of curvature, the results obtained using the surface impedance boundary condition model for conducting prolate and oblate spheroids of various axial ratios are in good agreement with the exact results for skin depths of about 1/5 of the semi-minor axis when calculating electromagnetic forces and for skin depths less than 1/20 of the semi-minor axis when calculating Joule losses.

Exact analytical solutions for the quasistationary electromagnetic fields in the presence of conducting objects require the field solutions both internal and external to the conductors. Such solutions are limited for certain canonically shaped objects but are useful in testing the accuracy of various approximate models and numerical methods developed to solve complex problems related to real world conducting objects and in calibrating instruments designed to measure various field quantities. Theoretical investigations of quasistationary electromagnetic fields also aid in improving the understanding of the physical phenomena of electromagnetic induction. This thesis presents rigorous analytical expressions derived as benchmark solutions for the quasistationary field quantities both inside and outside, Joule losses and the electromagnetic forces acting upon a conducting spheroid placed in the proximity of a non-uniform field produced by current-carrying turns. These expressions are used to generate numerous numerical results of specified accuracy and selected results are presented in a normalized form for extended ranges of the spheroid axial ratio, the ratio of the depth of penetration to the semi-minor axis and the position of the inducing turns relative to the spheroids. They are intended to constitute reference data to be employed for comprehensive comparisons of results from approximate numerical methods or from boundary impedance models used for real world conductors. Approximate boundary conditions such as the simpler perfect electric conductor model or the Leontovich surface impedance boundary condition model can be used to obtain approximate solutions by only analyzing the field external to the conducting object. The range of validity of these impedance boundary condition models for the analysis of axisymmetric eddy-current problems is thoroughly investigated. While the simpler PEC model can be employed only when the electromagnetic depth of penetration is much smaller than the smallest local radius of curvature, the results obtained using the surface impedance boundary condition model for conducting prolate and oblate spheroids of various axial ratios are in good agreement with the exact results for skin depths of about 1/5 of the semi-minor axis when calculating electromagnetic forces and for skin depths less than 1/20 of the semi-minor axis when calculating Joule losses.

La réduction des nuisances sonores induites par les aéronefs civils autour des grands aéroports est devenue un enjeu sociétal important. Pour réduire le bruit propulsif de soufflante, devenu prépondérant au cours des dernières années avec l'avènement de turboréacteurs à double flux, les constructeurs sont amenés à généraliser l'utilisation de matériaux absorbants acoustiques (également appelés « liners »). Ce sujet de thèse s'inscrit dans le cadre de l'amélioration des outils CAA relativement à la prévision numérique de ces matériaux absorbants. Cette modélisation soulève plusieurs interrogations, liées à divers aspects tels que le type d'écoulement mis en jeu (inhomogénéités, couche limite, etc.), les niveaux acoustiques en présence (effets de non linéarité), les effets de diffraction induits par les ruptures d'impédance, etc. Cette étude consiste donc à valider et à améliorer la condition limite d'impédance, implémentée dans le solveur CAA sAbrinA.v0, développé par l'Onera. Des développements théoriques sont d'abord consacrés à la modélisation de l'impédance dans le domaine temporel, et conduisent à une discussion sur la généralisation de cette modélisation. Le travail consiste ensuite à simuler plusieurs cas tests canoniques de l'absorption du bruit par un revêtement acoustique, lesquels sont validés par confrontation avec d'autres résultats analytiques et/ou expérimentaux. Ce travail fournit de nouvelles connaissances sur la façon dont les matériaux absorbants acoustiques peuvent être modélisés et simulés de manière précise dans le cadre d'une approche CAA dans le domaine temporel. / The reduction of acoustic emission induced by civil aircraft around major airports has become an important societal issue. To reduce the fan noise, induced by the engines, which has become preponderant over the past years with the advent of turbofan engines, manufacturers are led to generalize the employment of acoustic absorbing materials (or acoustic “liners”). The present thesis is related to the numerical prediction of such absorbing materials, in the context of time-domain CAA (Computational AeroAcoustics) methods. Such modeling raises several key questions, which are related to various aspects such as the type of flow involved (boundary layers effects, etc.), the sound levels considered (non-linear phenomena), the diffraction effects induced by ruptures of impedance, etc. The present study then consists in validating and improving the time-domain impedance boundary condition implemented in Onera’s structured CAA solver (named sAbrinA.v0). Theoretical developments are first devoted to the modeling of impedance in the time-domain, and lead to a discussion on the generalization of this modeling. The work then consists in CAA-simulating several canonical tests of noise absorption by acoustic liners. Outputs are compared against experimental and/or analytical results, delivering new insight in the way noise absorption materials can be accurately modeled and simulated using time-domain CAA-approaches.

L'équation intégrale du champ électrique (EFIE) est très utilisée pour résoudre des problèmes de diffusion d'ondes électromagnétiques grâce à la méthode aux éléments de frontière (BEM). En domaine fréquentiel, les systèmes matriciels émergeant de la BEM souffrent, entre autres, de deux problèmes de mauvais conditionnement : l'augmentation du nombre d'inconnues et la diminution de la fréquence entrainent l'accroissement du nombre de conditionnement. En conséquence, les solveurs itératifs requièrent plus d'itérations pour converger vers la solution, voire ne convergent pas du tout. En domaine temporel, ces problèmes sont également présents, en plus de l'instabilité DC qui entraine une solution erronée en fin de simulation. La discrétisation en temps est obtenue grâce à une quadrature de convolution basée sur les méthodes de Runge-Kutta implicites.Dans cette thèse, diverses formulations d'équations intégrales utilisant notamment des conditions d'impédance aux frontières (IBC) sont étudiées et préconditionnées. Dans une première partie en domaine fréquentiel, l'IBC-EFIE est stabilisée pour les basses fréquences et les maillages denses grâce aux projecteurs quasi-Helmholtz et à un préconditionnement de type Calderón. Puis une nouvelle forme d'IBC est introduite, ce qui permet la construction d'un préconditionneur multiplicatif. Dans la seconde partie en domaine temporel, l'EFIE est d'abord régularisée pour le cas d'un conducteur électrique parfait (PEC), la rendant stable pour les pas de temps larges et immunisée à l'instabilité DC. Enfin, unerésolution efficace de l'IBC-EFIE est recherchée, avant de stabiliser l'équation pour les pas de temps larges et les maillages denses. / The Electric Field Integral Equation (EFIE) is widely used to solve wave scattering problems in electromagnetics using the Boundary Element Method (BEM). In frequency domain, the linear systems stemming from the BEM suffer, amongst others, from two ill-conditioning problems: the low frequency breakdown and the dense mesh breakdown. Consequently, the iterative solvers require more iterations to converge to the solution, or they do not converge at all in the worst cases. These breakdowns are also present in time domain, in addition to the DC instability which causes the solution to be completely wrong in the late time steps of the simulations. The time discretization is achieved using a convolution quadrature based on Implicit Runge-Kutta (IRK) methods, which yields a system that is solved by Marching-On-in-Time (MOT). In this thesis, several integral equations formulations, involving Impedance Boundary Conditions (IBC) for most of them, are derived and subsequently preconditioned. In a first part dedicated to the frequency domain, the IBC-EFIE is stabilized for the low frequency and dense meshes by leveraging the quasi-Helmholtz projectors and a Calderón-like preconditioning. Then, a new IBC is introduced to enable the development of a multiplicative preconditioner for the new IBC-EFIE. In the second part on time domain,the EFIE is regularized for the Perfect Electric Conductor (PEC) case, to make it stable in the large time step regime and immune to the DC instability. Finally, the solution of the time domain IBC-EFIE is investigated by developing an efficient solution scheme and by stabilizing the equation for large time steps and dense meshes.

[pt] Os efeitos das variações laterais de um terreno irregular na propagação de ondas eletromagnéticas são considerados pela representação dos campos vetoriais em termo de dois potenciais escalares Hertzianos em coordenadas esféricas. A combinação da equação parabólica para esses potenciais com uma condição de contorno de impedância para o solo, seguida por uma transformação de variáveis, define um problema de condição de contorno caracterizado por equações exibindo coeficientes que dependem da função altura do terreno e de suas derivadas parciais. A solução do problema através do esquema de Crank-Nicolson leva a um sistema esparso de equações lineares que é resolvido por um método direto. O modelo numérico resultante é aplicado a terrenos irregulares, representando configurações tridimensionais hipotéticas. / [en] The effects from lateral variations of irregular terrain on the propagation of radio waves are considered by the representation of the vector fields in terms of two scalar Hertz potentials in spherical coordinates. The combination of three-dimensional parabolic equations for these potentials with an impedance boundary condition for the ground, followed by a transformation of variables, will define a boundary-condition problem characterized by equations displaying coefficients that depend on the terrain height function and its partial derivatives. The problem solution through the Crank-Nicolson scheme will lead to a sparse system of linear equations, which will be solved by a direct method. The resulting numerical model will be applied to irregular terrain, representing hypothetical three-dimensional configurations.