Les schémas de calcul déterministes permettent une modélisation à moindre coût du comportement de la population de neutrons en réacteur, mais sont traditionnellement construits sur des approximations (décomposition réseau/cœur, homogénéisation spatiale et énergétique…). La thèse revient sur une partie de ces sources d’erreur, de façon à rapprocher la méthode déterministe d’un schéma de référence. L’objectif est de profiter des architectures informatiques modernes (HPC) pour résoudre le problème neutronique à l’échelle du cœur 3D, tout en préservant l’opérateur de transport et une partie des hétérogénéités de la géométrie. Ce travail est réalisé au sein du solveur cœur Sn Minaret de la plateforme de calcul Apollo3® pour des réacteurs à neutrons rapides.Une méthode de décomposition de domaine en espace, est retenue. L'idée consiste à décomposer un problème de grande dimension en sous-problèmes "indépendants" de taille réduite. La convergence vers la solution globale est assurée par échange de flux angulaires entre sous-domaines au cours d'un processus itératif. En favorisant un recours massif au parallélisme, les méthodes de décomposition de domaine contribuent à lever les contraintes en mémoire et temps de calcul. La mise en place d'un parallélisme hybride, couplant les technologies MPI et OpenMP, est en particulier propice au passage sur supercalculateur. Une méthode d'accélération de type Coarse Mesh Rebalance  est ajoutée pour pallier à la pénalité de convergence constatée sur la méthode de décomposition de domaine. Le potentiel du nouveau schéma est finalement mis en évidence sur un coeur CFV 3D, construit en préservant l'hétérogénéité des assemblages absorbants. / Deterministic calculation schemes are devised to numerically solve the neutron transport equation in nuclear reactors. Dealing with core-sized problems is very challenging for computers, so much that the dedicated core codes have no choice but to allow simplifying assumptions (assembly- then core-scale steps…). The PhD work aims to correct some of these ‘standard’ approximations, in order to get closer of reference calculations: thanks to important increases in calculation capacities (HPC), nowadays one can solve 3D core-sized problems, using both high mesh refinement and the transport operator. Developments were performed inside the Sn core solver Minaret, from the new CEA neutronics platform Apollo3® for fast neutrons reactors of the CFV-kind.This work focuses on a Domain Decomposition Method in space. The fundamental idea involves splitting a core-sized problem into smaller and 'independent' subproblems. Angular flux is exchanged between adjacent subdomains. In doing so, all combined subproblems converge to the global solution at the outcome of an iterative process. Domain decomposition is well-suited to massive parallelism, allowing much more ambitious computations in terms of both memory requirements and calculation time. An hybrid MPI/OpenMP parallelism is chosen to match the supercomputers architecture. A Coarse Mesh Rebalance accelration technique is added to balance the convergence penalty observed using Domain Decomposition. The potential of the new calculation scheme is demonstrated on a 3D core of the CFV-kind, using an heterogeneous description of the absorbent rods.

Au cours des dernières décennies, les plasmas créés par une décharge micro-ondes ont de plusen plus attiré l’attention de la communauté scientifique aérospatiale sur le sujet du contrôled’écoulements. En effet, il a été démontré expérimentalement que le dépôt d’énergie obtenu parle plasma peut modifier les propriétés aérodynamiques de l’écoulement autour d’un objet telleque la trainée de frottement. Or, la conception et l’optimisation de ces actionneurs plasma entant que technique de contrôle d’écoulements nécessitent une compréhension approfondie de laphysique sous-jacente que les seules expériences sont incapables de fournir.Dans ce contexte, nous nous intéressons à la modélisation numérique de l’interaction desondes électromagnétiques avec un plasma et le gaz afin de mieux comprendre la nature desdécharges micro-ondes et leur applicabilité. La modélisation de ces phénomènes présente desdifficultés importantes en raison du couplage multi-physique et donc de la multitude des échellesspatiales et temporelles qui apparaissent. Ce travail de thèse traite des questions de physiqueet de mathématiques appliquées soulevées par la modélisation numérique de ces plasmas.La première partie du travail se focalise sur les questions de validité du modèle physique duclaquage micro-onde fondé sur l’approximation de champ effectif local. En raison des gradientsde densité du plasma très élevés, la validité du concept de champ effectif local peut être misen doute. Pour cela, un modèle fluide du second ordre est développé en incluant une equationd’énergie électronique non-locale. Cette modélisation permet de décrire de façon plus précisele dépôt d’énergie par plasma induisant la formation d’ondes de choc dans le gaz. Une analysedimensionnelle du système d’équations fluide permet de caractériser la non-localité en espace dubilan d’énergie électronique en fonction du champ électrique réduit et de la fréquence de l’onderéduite. Une discussion est également menée sur d’autres approximations des coefficients detransport électronique. Dans une deuxième partie, la construction et l’analyse d’une méthode multi-échelles derésolution numérique du problème de propagation des ondes électromagnétiques dans le plasmasont réalisées. Il s’agit du couplage entre les équations de Maxwell dans le domaine temporel avecune équation de quantité de mouvement pour les électrons. L’approche s’appuie sur la méthodede décomposition de domaine de type Schwartz, basée sur une formulation variationnelle duschéma de Yee et utilisant deux niveaux de grilles Cartésiennes emboitées. Une grille locale,appelée patch, est utilisée pour calculer de manière itérative la solution dans la région du plasmaoù une meilleure précision est requise. La méthode proposée permet le raffinement local etdynamique du maillage spatial tout en conservant l’énergie du système. Une analyse théorique dela convergence de l’algorithme pour les résolutions temporelles explicite et implicite est égalementréalisée. Dans la dernière partie, des simulations numériques sur le claquage micro-ondes et la formation de structures filamentaires de plasma sont conduites. Les effets de différents types d’approximations sur le modèle physique du plasma sont analysés. Puis, ces expériences numériques démontre la précision et l’efficacité, en terme de temps de calcul, de la méthode multi-échelleproposée. Enfin, on étudie les effets de chauffage du gaz sur la formation et l’entretien de structures filamentaires dans l’air à pression atmosphérique. Pour cela, le modèle micro-onde-plasma développé est couplé avec les équations de Navier-Stokes instationnaires pour les écoulements compressibles. Les simulations montrent des caractéristiques intéressantes de la dynamique deces structures plasma pendant le processus de chauffage du gaz, qui sont en accord étroit avec les données expérimentales. / In recent decades, microwave discharge plasmas have attracted increasing attention of aerospace scientific community to the subject of aerodynamic flow control because of their capability of sub- stantially modifying the properties of the flow around bodies by effective energy deposition. The design and optimization of these plasma actuators as flow control technique require a compre- hensive understanding of the complex physics involved that the sole experiments are incapable to provide.In this context, we have interest in the numerical modeling of the mutual interaction of elec- tromagnetic waves with plasma and gas in order to better understand the nature of microwave discharges and their applicability. A challenging problem arises when modeling such phenomena because of the coupling of different physics and therefore the multiplicity of spatial and tempo- ral scales involved. A solution is provided by this thesis work which addresses both physics and applied mathematics questions related to microwave plasma modeling.The first part of this doctorate deals with validity matters of the physical model of microwave breakdown based on the local effective field concept. Because of large plasma density gradients, the local effective field approximation is questionable and thus a second-order plasma fluid model is developed, where the latter approximation is replaced by the local mean energy approximation. This modeling approach enables to take into account the non-locality in space of the electron energy balance that provides a more accurate description of the energy deposition by microwave plasma leading to the shock waves formation into the gas. A dimensionless analysis of the plasma fluid system is performed in order to theoretically characterize the non-locality of the introduced electron energy equation as function of the reduced electric field and wave frequency. It also discusses other approximations related to the choice and method of calculation of electron transport coefficients.Concerning the mathematical aspects, the thesis work focuses on the design and the analysis of a multiscale method for numerically solving the problem of electromagnetic wave propagation in microwave plasma. The system of interest consists of time-dependent Maxwell’s equations coupled with a momentum transfer equation for electrons. The developed approach consists of a Schwartz type domain decomposition method based on a variational formulation of the standard Yee’s scheme and using two levels of nested Cartesian grids. A local patch of finite elements is used to calculate in an iterative manner the solution in the plasma region where a better precision is required. The proposed technique enables a conservative local and dynamic refinement of the spatial mesh. The convergence behavior of the iterative resolution algorithm both in an explicit and implicit time-stepping formulation is then analyzed.In the last part of the doctorate, a series of numerical simulations of microwave breakdown and the filamentary plasma array formation in air are performed. They allow to study in detail the consequences of the different types of physical approximations adopted in the plasma fluid model. Then, these numerical experiments demonstrate the accuracy and the computational efficiency of the proposed patch correction method for the problem of interest. Lastly, a numerically investigation of the effects of gas heating on the formation and sustaining of the filamentary plasma array in atmospheric-pressure air is carried out. For doing this, the developed microwave-plasma model is coupled with unsteady Navier-Stokes equations for compressible flows. The simulations provide interesting features of the plasma array dynamics during the process of gas heating, in close agreement with experimental data.

Notre sujet de recherche porte sur la résolution du problème de transport avec coûts fixes (FCTP). Le problème de transport classique consiste à déterminer le schéma optimal de distribution dans un réseau. Le réseau est divisé en deux sous-ensembles de sommets : les origines caractérisées par une offre et les destinations caractérisées par une demande. À chaque arc reliant une origine et une destination, est associé un coût variable. L'objectif est de satisfaire toutes les demandes en minimisant la somme des coûts variables de transport. Dans le FCTP, il y a aussi des coûts fixes associés à tous les arcs, en supplément de tout ce qui décrit un problème de transport classique. Ainsi, à chaque arc, est associé un coût variable et un coût fixe qui est considéré si et seulement si l'arc est utilisé. L'objectif est désormais de minimiser la somme totale des coûts, variables et fixes. Le FCTP nous confronte donc à un modèle différent et plus complexe. La complexité de résolution est accrue pour les instances de grande taille. Dans ce mémoire, nous étudions et présentons une nouvelle méthode de résolution pour les instances de grande taille du FCTP. Il s'agit d'une méthode de décomposition lagrangienne qui utilise la relaxation lagrangienne et un algorithme de sous-gradient pour trouver une borne inférieure au problème global. Nous avons intégré à la méthode une heuristique lagrangienne, incluant une procédure de ``slope scaling'' afin d'améliorer notre algorithme de sous-gradient et le résultat final de la méthode. À l'issue de notre processus de résolution, nous trouvons, pour certaines instances de grande taille, un moyen d'améliorer la solution proposée par CPLEX pour le FCTP en donnant comme paramètre à CPLEX la solution finale de notre méthode. / Our research topic focuses on solving the fixed charge transportation problem (FCTP). The classic transportation problem is to determine the optimal distribution pattern in a network. The network is divided into two subsets of vertices : the origins characterized by a supply and the destinations characterized by a demand. Each arc connecting an origin and a destination has a variable cost associated with it. The objective is to satisfy all demands while minimizing the sum of variable transportation costs. In FCTP, there are also fixed costs associated with all arcs, in addition to all others things that describe a typical transportation problem. So, each arc is associated a variable cost and a fixed cost which is considered if and only if the arc is used. The objective is now to minimize the total sum of costs, variable and fixed. The FCTP therefore confronts us with a different and more complex model. The resolution complexity is even increased for large instances. In this thesis, we study and present a new resolution method for large FCTP instances. This is a lagrangian decomposition method which uses lagrangian relaxation and a subgradient algorithm to find a lower bound to the global problem. We have integrated into the method a lagrangian heuristic, including a “slope scaling” procedure in order to improve our sub-gradient algorithm and the final result of the method. At the end of our resolution process, we find, for some large instances, a way to improve the solution proposed by CPLEX for the FCTP by giving as parameter to CPLEX the final solution of our method.

Lors d'un hypothétique accident grave de réacteur à eau sous pression, un mélange de matériaux fondus, appelé corium, issu de la fusion du cœur peut se relocaliser dans le puits de cuve constitué par un radier en béton. Les codes d'évaluation réacteur pour simuler la phénoménologie de l'interaction corium-béton sont basés sur une description à grande échelle des échanges qui soulève de nombreuses questions, tant sur la prise en compte des phénomènes multi-échelles mis en jeu que sur la structure adoptée de la couche limite au voisinage du front d'ablation. Dans ce contexte, l'objectif principal de ce travail consiste à aborder le problème de la structure de la couche limite par simulation numérique directe. Ce travail s'inscrit dans le cadre plus général d'une description et d'une modélisation multi-échelle des échanges, c'est-à-dire de l'échelle locale associée au voisinage du front d'ablation jusqu'à l'échelle du code d'évaluation réacteur. Une telle description multi-échelle des échanges soulève le problème de la description locale de l'écoulement multiphasique multiconstituant mais aussi le problème du changement d'échelle et en particulier le passage de l'échelle locale à l'échelle de description supérieure dite macroscopique associée aux mouvements convectifs dans le bain de corium. Parmi les difficultés associées au changement d'échelle, nous nous intéressons à la problématique de la construction de conditions aux limites effectives ou lois de parois pour les modèles macroscopiques. Devant la complexité du problème multiphasique multiconstituant posé au voisinage du front, cette contribution a été abordée sur un problème modèle. Des conditions aux limites dites effectives ont été construites dans le cadre d'une méthode de décomposition de domaine puis testées pour un problème d'écoulement laminaire de convection naturelle sur parois rugueuses. Mˆeme si le problème traité reste encore éloigné des applications visées, cette contribution offre de nombreuses perspectives et constitue une première étape d'une modélisation multiéchelle des échanges pour la problématique de l'interaction corium-béton. Dans le cas plus complexe des écoulements multiphasiques multiconstituants et devant les difficultés expérimentales associées, le développement de lois de parois pour les outils existants aux échelles de description supérieures nécessite, au préalable, de disposer d'un outil de simulation numérique directe de l'écoulement au voisinage du front d'ablation. L'outil développé dans ce travail correspond à un modèle de Cahn-Hilliard/Navier-Stokes pour un mélange diphasique (liquide-gaz) compositionnel (corium-béton fondu) s'appuyant sur une description du système selon trois paramètres d'ordre associés respectivement aux fractions volumiques du gaz et aux deux espèces miscibles de la phase liquide ainsi que sur une décomposition de l'énergie libre selon une contribution diphasique et compositionnelle. Les équations de transport sont dérivées dans le cadre de la thermodynamique des processus irréversibles et résolues sur la base d'une application éléments finis de la plate-forme PELICANS. Plusieurs expériences numériques illustrent la validité et les potentialités d'application de cet outil sur des problèmes diphasiques et/ou compositionnels. Enfin, à partir de l'outil développé, nous abordons par simulation numérique directe une étude de la structure de la couche limite au voisinage du front d'ablation pour des bétons siliceux et silico-calcaire. / In the late phases of some scenario of hypothetical severe accident in Pressurized Water Reactors, a molten mixture of core and vessel structures, called corium, comes to interact with the concrete basemat. The safety numerical tools are lumped parameter codes. They are based on a large averaged description of heat and mass transfers which raises some uncertainties about the multi-scale description of the exchanges but also about the adopted boundary layer structure in the vicinity of the ablation front. In this context, the aim of this work is to tackle the problem of the boundary layer structure by means of direct numerical simulation. This work joins within the more general framework of a multi-scale description and a multi-scale modeling, namely from the local scale associated with the vicinity of the ablation front to the scale associated with the lumped parameter codes. Such a multi-scale description raises not only the problem of the local description of the multiphase multicomponent flow but also the problem of the upscaling between the local- and the macro-scale which is associated with the convective structures within the pool of corium. Here, we are particularly interested in the building of effective boundary conditions or wall laws for macro-scale models. The difficulty of the multiphase multicomponent problem at the local scale leads us to consider a relatively simplified problem. Effective boundary conditions are built in the frame of a domain decomposition method and numerical experiments are performed for a natural convection problem in a stamp shaped cavity to assess the validity of the proposed wall laws. Even if the treated problem is still far from the target applications, this contribution can be viewed as a first step of a multi-scale modeling of the exchanges for the molten core concrete issue. In the more complicated case of multiphase multicomponent flows, it is necessary to have a direct numerical simulation tool of the flow at the local scale to build wall laws for macro-scale models. Here, the developed tool corresponds to a Cahn-Hilliard/Navier-Stokes model for a two-phase compositional system. It relies on a description of the system by three volume fractions and on a free energy composed by a two-phase part and a compositional part. The governing equations are derived in the frame of the thermodynamic of irreversible processes. They are solved on the basis of a finite element application of the object-oriented software component library PELICANS. Several numerical experiments illustrate the validity and the potentialities of application of this tool on two-phase compositional problems. Finally, using the developed tool, we tackle by means of direct numerical simulation the problem boundary layer structure in the vicinity of the ablation front for limestone-sand and siliceous concretes.

Le contexte d'étude est celui de la Compatibilité ÉlectroMagnétique (CEM). L'objectif de la CEM est, comme son nom l'indique, d'assurer la compatibilité entre une source de perturbation électromagnétique et un système électronique victime. Or, la prédiction de ces niveaux de perturbation ne peut pas s'effectuer à l'aide d'un simple calcul analytique, en raison de la géométrie qui est généralement complexe pour le système que l'on étudie, tel que le champ à l'intérieur d'un cockpit d'avion par exemple. En conséquence, nous sommes contraints d'employer des méthodes numériques, dans le but de prédire ce niveau de couplage entre les sources et les victimes. Parmi les nombreuses méthodes numériques existantes à ce jour, les méthodes Multi-Domaines (MD) sont très prisées. En effet, elles offrent la liberté aux utilisateurs de choisir la méthode numérique la plus adaptée, en fonction de la zone géométrique à calculer. Au sein de ces méthodes MD, la " Domain Decomposition Method " (DDM) présente l'avantage supplémentaire de découpler chacun de ces domaines. En conséquence, la DDM est particulièrement intéressante, vis-à-vis des méthodes concurrentes, en particulier sur l'aspect du coût numérique. Pour preuve, l'ONERA continue de développer cette méthode qui ne cesse de montrer son efficacité depuis plusieurs années, notamment pour le domaine des Surfaces Équivalentes Radar (SER) et des antennes. L'objectif de l'étude est de tirer profit des avantages de cette méthode pour des problématiques de CEM. Jusqu'à maintenant, de nombreuses applications de CEM, traitées par le code DDM, fournissaient des résultats fortement bruités. Même pour des problématiques électromagnétiques très simples, des problèmes subsistaient, sans explication convaincante. Ceci justifie cette étude. Le but de cette thèse est de pouvoir appliquer ce formalisme DDM à des problématiques de CEM. Dans cette optique, nous avons été amenés à redéfinir un certain nombre de conventions, qui interviennent au sein de la DDM. Par ailleurs, nous avons développé un modèle spécifique pour les ouvertures, qui sont des voies de couplage privilégiées par les ondes, à l'intérieur des cavités que représentent les blindages. Comme les ouvertures sont, en pratique, de petites dimensions devant la longueur d'onde, on s'est intéressé à un modèle quasi-statique. Nous proposons alors un modèle, qui a été implémenté, puis validé. Suite à ce modèle, nous avons développé une méthode originale, basée sur un calcul en deux étapes, permettant de ne plus discrétiser le support des ouvertures dans les calculs 3D.

Cette thèse a pour but de développer des approches déterministes et heuristiques pour résoudre certaines classes des problèmes d'optimisation en télécommunication et la mobilité d'un réseau de transport : problèmes de routage, problèmes de covoiturage, problèmes de contrôle de l'alimentation dans un réseau sans fil, problèmes d'équilibrage du spectre dans les réseaux DSL. Il s'agit des problèmes d'optimisation non convexe de très grande taille. Nos approches sont basées sur la programmation DC&DCA, méthode de décomposition proximale et la méthode d'étiquetage des graphes. Grâce aux techniques de formulation/reformulation et de pénalité exacte, nous avons établi des programmes DC équivalents en vue de leur résolution par DCA. Selon la structure de ces problèmes, on peut fournir des décompositions DC appropriées ou de bons points initiaux de DCA. Nos méthodes ont été programmées sous MATLAB, C/C++. Ils montrent la performance de nos algorithmes par rapport à des méthodes existantes. / In this thesis, we focus on developing deterministic and heuristic approaches for solving some classes of optimization problems in Telecommunication and Mobility & Transport domain: Routing problems, Car pooloing problems, Power control problems in wireless network, Optimal spectrum balancing problems in DSL networks. They are large-scale nonconvex optimization problems. Our methodologies are focus on DC programming and DCA, Proximal decomposition method and Labeling method in graph theory. They are well-known as powerful tools in optimization. The considered problems were reformulated using the DC formulation/reformulation and exact penalty techniques and the DCA was used to obtain the solution. Also, depending on the structure of considered problems, we can provide appropriate DE decompositions or good initial points for DCA. All these proposed methods have been implemented with MATLAB, C/C++ to confirm the practical aspects and enhance our research works.

Cette thèse aborde deux thématiques différentes. On s’intéresse d’abord au développement et à l’analyse de méthodes pour le sens électrique appliqué à la robotique. On considère en particulier la méthode des réflexions permettant, à l’image de la méthode de Schwarz, de résoudre des problèmes linéaires à partir de sous-problèmes plus simples. Ces deniers sont obtenus par décomposition des frontières du problème de départ. Nous en présentons des preuves de convergence et des applications. Dans le but d’implémenter un simulateur du problème direct d’électrolocation dans un robot autonome, on s’intéresse également à une méthode de bases réduites pour obtenir des algorithmes peu coûteux en temps et en place mémoire. La seconde thématique traite d’un problème inverse dans le domaine de la chimie quantique. Nous cherchons ici à déterminer les caractéristiques d’un système quantique. Celui-ci est éclairé par un champ laser connu et fixé. Dans ce cadre, les données du problème inverse sont les états avant et après éclairage. Un résultat d’existence locale est présenté, ainsi que des méthodes de résolution numériques. / This thesis tackle two different topics.We first design and analyze algorithms related to the electrical sense for applications in robotics. We consider in particular the method of reflections, which allows, like the Schwartz method, to solve linear problems using simpler sub-problems. These ones are obtained by decomposing the boundaries of the original problem. We give proofs of convergence and applications. In order to implement an electrolocation simulator of the direct problem in an autonomous robot, we build a reduced basis method devoted to electrolocation problems. In this way, we obtain algorithms which satisfy the constraints of limited memory and time resources. The second topic is an inverse problem in quantum chemistry. Here, we want to determine some features of a quantum system. To this aim, the system is ligthed by a known and fixed Laser field. In this framework, the data of the inverse problem are the states before and after the Laser lighting. A local existence result is given, together with numerical methods for the solving.

Le contexte d’étude est celui de la Compatibilité ÉlectroMagnétique (CEM). L’objectif de la CEM est, comme son nom l’indique, d’assurer la compatibilité entre une source de perturbation électromagnétique et un système électronique victime. Or, la prédiction de ces niveaux de perturbation ne peut pas s’effectuer à l’aide d’un simple calcul analytique, en raison de la géométrie qui est généralement complexe pour le système que l’on étudie, tel que le champ à l’intérieur d’un cockpit d’avion par exemple. En conséquence, nous sommes contraints d’employer des méthodes numériques, dans le but de prédire ce niveau de couplage entre les sources et les victimes. Parmi les nombreuses méthodes numériques existantes à ce jour, les méthodes Multi-Domaines (MD) sont très prisées. En effet, elles offrent la liberté aux utilisateurs de choisir la méthode numérique la plus adaptée, en fonction de la zone géométrique à calculer. Au sein de ces méthodes MD, la « Domain Decomposition Method » (DDM) présente l’avantage supplémentaire de découpler chacun de ces domaines. En conséquence, la DDM est particulièrement intéressante, vis-à-vis des méthodes concurrentes, en particulier sur l’aspect du coût numérique. Pour preuve, l’ONERA continue de développer cette méthode qui ne cesse de montrer son efficacité depuis plusieurs années, notamment pour le domaine des Surfaces Équivalentes Radar (SER) et des antennes. L’objectif de l’étude est de tirer profit des avantages de cette méthode pour des problématiques de CEM. Jusqu’à maintenant, de nombreuses applications de CEM, traitées par le code DDM, fournissaient des résultats fortement bruités. Même pour des problématiques électromagnétiques très simples, des problèmes subsistaient, sans explication convaincante. Ceci justifie cette étude. Le but de cette thèse est de pouvoir appliquer ce formalisme DDM à des problématiques de CEM. Dans cette optique, nous avons été amenés à redéfinir un certain nombre de conventions, qui interviennent au sein de la DDM. Par ailleurs, nous avons développé un modèle spécifique pour les ouvertures, qui sont des voies de couplage privilégiées par les ondes, à l’intérieur des cavités que représentent les blindages. Comme les ouvertures sont, en pratique, de petites dimensions devant la longueur d’onde, on s’est intéressé à un modèle quasi-statique. Nous proposons alors un modèle, qui a été implémenté, puis validé. Suite à ce modèle, nous avons développé une méthode originale, basée sur un calcul en deux étapes, permettant de ne plus discrétiser le support des ouvertures dans les calculs 3D. / The context of the study is the ElectroMagnetic Compatibility (EMC). Principal aim of the EMC is to ensure the compatibility between an electromagnetic perturbance source and an electronic device victim. Unfortunately, the perturbation levels prediction can not be made using an analytic formula, because the geometry which is generally complex for the interesting system, for example the field inside an aircraft’s cockpit. Therefore, we are contrained to use numerical methods, to be able to evaluate this coupling level between sources and victims. Among several existing numerical methods, Multi-Domains (MD) methods are very interesting. They offer to users the freedom to choose the most powerfull numerical method, in terms of the geometrical zone evaluated. With the MD methods, « Domain Decomposition Method » (DDM) has the avantage of decouplingeach of theses areas. Therefore, DDM is very interesting, compared to other methods, in particular on the numerical cost. ONERA keeps on developing this method, which has not stop showing his efficiency since several years, in particular in Radar Cross Section (RCS) and antennas. The objective of this study is to take the benefits of this method for EMC problems. Up to now, several EMC applications treated by the DDM code provided results strongly noisy. Even for with very simple electromagnetic cases, some problems remained without convincing explanations. This justifies this study. The aim of this thesis is to can be able to apply DDM formalism to EMC problems. Then, we have been induced to redefine a number of conventions which are involved in the DDM. Otherwise, we have developed a specific model for the apertures which are privilegied tracts of the coupling by the penetration of waves inside cavities (shieldings). As the apertures have in practice smaller dimensions compared to the wavelength, we have been interested to a quasistatic model which was developped, implemented and validated. Following this model, we have developed an original method, based on a two step calculation, able to do not discretize the apertures support in 3D computations.