Dans cette thèse, nous nous intéressons à l’amélioration du schéma de Yee pour traiter demanière plus efficace et pertinente les problèmes industriels auxquels nous sommes confrontés à l’heureactuelle. Pour cela, nous cherchons avant tout à diminuer les erreurs numériques de dispersion et àaméliorer les modélisations des géométries courbes ainsi que des réseaux de câbles. Pour répondre àces besoins, une solution basée sur un schéma Galerkin discontinu pourrait être envisagée. Toutefois,l’utilisation d’une telle technique sur la totalité du volume de calcul est relativement coûteuse. De plus,la prise en compte de structures filaires sur un tel schéma n’est pas encore opérationnelle. C’est pourquoi,dans l’optique d’avoir un outil industriel, et après une étude bibliographique, nous nous sommes plutôtorientés sur l’étude d’un schéma éléments finis (FEM) sur maillage cartésien qui possède toutes lesbonnes propriétés du schéma de Yee. Notamment, à l’ordre d’approximation spatiale égal à 0 ce schémaFEM est exactement le schéma de Yee, et, pour des ordres supérieurs, il permet de réduire fortementl’erreur de dispersion numérique de ce dernier. Dans le travail de cette thèse, pour ce schéma, nous avons notamment donné un critère de stabilité théorique, étudié sa convergence théorique et fait une analyse de l’erreur de dispersion. Pour tenircompte des possibilités d’ordre d’approximation spatiale variable par direction, nous avons mis en placeune stratégie d’affectation des ordres suivant le maillage donné. Ceci nous a permis d’obtenir un pas detemps optimal pour une précision souhaitée tout en réduisant les coûts de calcul. Après avoir porté ceschéma sur des machines de production, différents problèmes de CEM, antennes, IEM ou foudre ont ététraités afin de montrer les avantages et le potentiel de celui-ci. En conclusion de ces expérimentationsnumériques, il s’avère que la méthode est limitée par le manque de précision pour prendre en comptedes géométries courbes. Afin d’améliorer cela, nous avons proposé une hybridation entre ce schéma et leschéma GD que l’on peut étendre aux autres schémas comme les méthodes différences finies (FDTD) etvolumes finis (FVTD). Nous avons montré que la technique d’hybridation proposée conserve l’énergie etest stable sous une condition que nous avons évaluée de manière théorique. Des exemples de validationont ensuite été montrés. Enfin, pour tenir compte des réseaux de câbles, un modèle de fils minces d’ordred’approximation spatiale élevé a été proposé. Malheureusement, celui-ci ne peut pas couvrir l’ensembledes cas industriels et pour remédier à cela, nous avons proposé une hybridation de notre approche avecune équation de ligne de transmission. L’intérêt de cette hybridation a été montré sur un certain nombred’exemples, que nous n’aurions pas pu traiter par un modèle de structure filaire simple. / In this thesis, we study the improvement of the Yee’s scheme to treat efficiently and in arelevant way the industrial issues we are facing nowadays. For that, we first of all try to reduce thenumerical errors of dispersion and then to improve the modeling of the curved surfaces and of theharness networks. To answer these needs, a solution based on a Galerkin Discontinuous (GD) methodhas been first considered. However, the use of such method on the entire modeling volume is quite costly ;moreover the wires are not taken into account in this method. That is the reason why, with the objectiveof an industrial tool and after a large bibliographic research, we headed for the study of finite elementsscheme (FEM) on a Cartesian mesh which has all the good properties of the Yee’s scheme. Especially,this scheme is exactly the Yee’s scheme when the spatial order of approximation is set to zero. Forthe higher orders, this new scheme allows to greatly reduce the numerical error of dispersion. In theframe of this thesis and for this scheme, we give a theoretical criterion of stability, study its theoreticalconvergence and we perform an analysis of the error of dispersion. To take into account the possibilityof the variable spatial orders of approximation in each direction, we put in place a strategy of orderaffectation according to the given mesh. This strategy allows to obtain an optimal time step for a givenselected precision while reducing the cost of the calculations. Once this new scheme has been adaptedto large industrial computing means, different EMC, antennas, NEMP or lightning problems are treatedto demonstrate the advantages and the potential of this scheme. As a conclusion of these numericalsimulations we demonstrate that this method is limited by a lack of precision when taking into accountcurved geometries. To improve the treatment of the curved surfaces, we propose an hybridization between this scheme andthe GD scheme. This hybridization can also be applied to other methods such as Finite Differences(FDTD) or Finite Volumes (FVTD). We demonstrate that the technique of hybridization proposed,allows to conserve the energy and is stable under a condition that we study theoretically. Some examplesare presented for validation. Finally and to take into account the cables, a thin wire model with a highorder of spatial approximation is proposed. Unfortunately, this model does not allow to cover all theindustrial cases. To solve this issue we propose an hybridization with a transmission line method. Theadvantage of this hybridization is demonstrated thanks to different cases which would not have beenfeasible with a more simple thin wire method.

Les travaux de cette thèse concerne l’étude d’une méthode élément finis d’ordre spatial élevé que l’on peut assimilé à une extension du schéma de Yee. On parle alors de méthode différences finies d’ordre élevé. Après avoir donné, dans un premier chapitre, un historique non exhaustif des principales méthodes utilisées pour résoudre les équations de Maxwell dans le cadre de problèmes de CEM et montré l’ intérêt de disposer d’un solveur de type "différences finies d’ ordre élevé", nous présentons dans un deuxième chapitre le principe de la méthode. Nous donnons pour cela les caractéristiques du schéma spatial et temporel en précisant les conditions de stabilité de la méthode. En outre, dans une étude purement numérique, nous étudions la convergence du schéma. On se focalise ensuite sur la possibilité d’utiliser des ordres spatiaux variable par cellules dans chaque direction de l’espace. Des comparaisons avec le schéma de Yee et un schéma de Galerkin Discontinu particulier sont ensuite effectuées pour montrer les gains en coûts calcul et mémoire et donc l’intérêt de notre approche par rapport aux deux autres. Dans un troisième chapitre, nous nous intéressons à l’étude de modèles physiques indispensable au traitement d’un problème de CEM. Pour cela, nous nous focalisons particulièrement sur un modèle de fil mince oblique, des modèles de matériaux volumiques et minces et enfin sur la prise en compte de sol parfaitement métallique dans une agression de type onde plane. Chaque modèle est détaillé et validé par comparaison avec des solutions analytiques ou résultant de la littérature, sur des exemples canoniques. Le quatrième chapitre est dédié à une technique d’hybridation entre notre méthode et une approche Galerkin Discontinu en vue de traiter des géométries possédant des courbures. Nous donnons pour cela une stratégie d’hybridation basée sur l’échange de flux qui garantie au niveau continue la conservation d’une énergie. Nous présentons ensuite quelques exemples montrant la validité de notre approche dans une stratégie multi-domaines/multi-méthodes que nous précisons. Enfin le dernier chapitre de cette thèse concerne l’exploitation de notre méthode sur des cas industriels en comparaisons avec d’autres méthodes ou des résultats expérimentaux. / This thesis is about the study of a high spatial finite element method whichcan be assimilated at an extension of the Yee schema. In the next, this method is also called high order finite difference method. In the first chapter, we give a non exhaustive recall of the major methods used to treat EMC problems and we show the necessity to have this kind of schema to simulate efficiently some EMC configurations. In the second chapter, the principle of the numerical method is presented and a stability condition is given. A numerical study analysis of the schema convergence is also done. Next, we show the interest to have the possibility to use local spatial order by cell in each direction of the computational domain. Some canonic examples are given to show the advantages interms of CPU time and memory storage of the method by comparison with Yee’s scheme and DG approach. In the third chapter, we define and validate on several examples,some physical models as thin wire, materials and perfectly metallic ground in presence of a plane wave, to have the possibility to treat EMC problems. The fourth chapter is about a hybridization strategy between our high order FDTD method and a DG schema.We focalize our study on a hybrid method which provides an energy conservation of the continuous problem. A numerical example is given to validate the method. Finally, in the last chapter, we present some simulations on industrial problems to show the possibility of the method to treat realistic EMC problems.