Le traitement des systèmes pleins issus de modélisation d'interactions de type particulaires est un problème lourd que les méthodes multipôles proposent d'optimiser grâce à des approximations hiérarchiques basées sur la régularité des potentiels sous-jacents. Cette thèse offre des développements sur une nouvelle méthode multipôle pour l'électromagnétisme à la fois en terme d'un traitement unifié de l'ensemble des fréquences et informatique dans le cadre du paralléélisme et d'optimisations algorithmiques telles que le partitionnement par graphes et le recouvrement des communications par des calculs.

Le travail effectué dans cette thèse a consisté à analyser différents aspects mathématiques et numériques d'une stratégie de résolution des problèmes de propagation d'onde acoustique et électromagnétique en domaine extérieur. Nous nous intéressons plus particulièrement à la méthode de couplage entre éléments finis et représentation intégrale (CEFRI) où nous analysons un algorithme de résolution itérative par analogie avec une méthode de décomposition de domaine ainsi que l'utilisation de la méthode multipôles rapide (FMM). Le système à résoudre fait intervenir des opérateurs intégraux ce qui rend crucial le recours à des méthodes rapides telles que la FMM. L'analogie avec une méthode de décomposition de domaine s'obtient par extension au problème de Maxwell des résultats établis par F. Ben Belgacem et al. pour le problème de Helmholtz posé en domaine non borné. Pour cela, nous avons montré le lien entre la méthode CEFRI et la méthode de Schwarz avec recouvrement total pour la résolution du problème de Maxwell en domaine non borné. Cette relecture de la méthode CEFRI offre également une technique de préconditionnement pour les solveurs de Krylov et nous a permis d'avoir une idée préliminaire sur la convergence de ces méthodes. Ainsi, nous nous intéressons plutôt à des méthodes itératives rapides. Pour cela, nous avons mené une analyse théorique afin de montrer la convergence superlinéaire du GMRES dans une configuration sphérique. La validation de ces aspects a été réalisée par l'enrichissement de nombreux intégrants de la librairie éléments finis Mélina++, en C++. / We are concerned with the study of different aspects of a numerical strategy for the resolution of acoustic and electromagnetic scattering problems. We focus more particu- larly on a coupling of finite element and integral representation (CEFRI) : we study an iterative algorithm by analogy with a domain decomposition method, and consider the use of the Fast Multipole Method (FMM). The system to be solved involves integral operators which requires the use of fast methods such as the FMM. The correspondence with a domain decomposition method is obtained by extending to the exterior Maxwell problem the results derived by F. Ben Belgacem et al. for the Helmholtz problem posed in unbounded domain. To this aim, we show the analogy to the Schwarz method with total overlap. This interpretation of CEFRI suggests a preconditioner for Krylov solvers and enables us to have a preliminary idea of their convergence. We derive in this context an analytical proof of a superlinear convergence of GMRES in a spherical configuration. The validation of these aspects has been achieved by the enrichment of the finite element library Mélina++ in C++.

Nous réalisons dans cette thèse l'étude d'interactions intermoléculaires d'un point de vue d'orbitales localisées. Cela concerne d'une part les interactions intermoléculaires produites par des orbitales localisées occupées comme l'interaction électrostatique, et d'autre part les interactions intermoléculaires qui engagent aussi les orbitales localisées virtuelles comme l'interaction de dispersion. Nous évaluons dans un premier temps les interactions électrostatiques produites par des distributions multipolaires en orbitales localisées, qui donnent une représentation chimique intuitive d'une densité de charges moléculaire. Nous montrons que les distributions multipolaires en orbitales localisées sont raisonnables pour décrire les interactions électrostatiques des densités de charges gelées si l'interaction multipolaire est tronquée à un ordre bien choisi, ce qui rend les orbitales localisées potentiellement intéressantes pour modéliser les interactions électrostatiques dans un champ de force. Nous utilisons ensuite les propriétés des orbitales localisées a priori dans un complexe pour définir une référence multipolaire dans le cas de l'interaction électrostatique des densités de charges relaxées. Nous évaluons ensuite la capacité de distributions multipolaires issues d'orbitales localisées relaxées pour décrire l'interaction électrostatique relaxée. Dans un second temps, localiser les orbitales occupées et virtuelles dans un cadre intermoléculaire nous permet d'une part d'attribuer des orbitales à chaque fragment d'un système fragmenté non covalent, et donc de diviser les excitations en classes et sélectionner uniquement les excitations les plus importantes à l'énergie de corrélation intermoléculaire post-Hartree-Fock. Nous proposons deux méthodes différentes que nous avons développé dans cette thèse pour sélectionner des excitations dans le cadre général de la DFT à séparation de portée couplée à l'approximation des phases aléatoires (RPA). La première méthode de sélection est basée sur un simple critère énergetique tandis que la seconde est basée sur la sélection d'une seule classe d'excitation à savoir la classe de dispersion. Enfin, nous montrons l'intérêt et les limites de ces deux méthodes de sélection pour des complexes à interactions variées.

Dans cette thèse, nous nous intéressons à la modélisation mathématique des hétérogénéités de longueurs caractéristiques beaucoup plus petites que la longueur d'ondes. La thèse consiste en deux parties. La partie théorique est dédiée à l'obtention d'un développement asymptotique raccordé: la solution est décrite à l'aide d'un développement de champ proche au voisinage de l'obstacle et par un développement de champ lointain hors de ce voisinage. Le développement de champ lointain met en jeu des solutions singulières de l'équation des ondes tandis que le champ proche lui est régi par un modèle quasi-statique. Ces deux développements sont alors raccordés dans une zone intermédiaire dite de raccord. Nous obtenons alors des estimations d'erreurs permettant de rendre rigoureux ce développement asymptotique formel. La deuxième partie est numérique. Elle décrit à la fois la méthode de Galerkine discontinue, une méthode de raffinement de maillage espace-temps et propose une discrétisation des modèles asymptotiques obtenues précédemment. Elle est illustrée par un certain nombre de tests numériques. / In this thesis, we focus our attention on the modeling of heterogeneities which are smaller than the wavelength. The document is decomposed into two parts : a theoretical one and a numerical one. In the first part, we derive a matched asymptotic expansion composed of a far-field expansion and a near-field expansion. The terms of the far-field expansion are singular solutions of the wave equation whereas the terms of the near-field expansion satisfy quasistatic problems. These expansions are matched in an intermediate region. We justify mathematically this theory by proving error estimates. In the second part, we describe the Discontinuous Galerkin method, a local time stepping method and the implementation of the matched asymptotic method. Numerical simulations illustrate these results.