Le but de cette thèse est d’étudier certaines équations aux dérivées partielles hyperboliques-dispersives. Une part importante est consacrée à l’analyse numérique et plus particulièrement à la convergence de schémas aux différences finies pour l’équation de Korteweg-de Vries et les systèmes abcd de Boussinesq. L’étude numérique suit les étapes classiques de consistance et de stabilité. Nous transposons au niveau discret la propriété de stabilité fort-faible des lois de conservations hyperboliques. Nous déterminons l’ordre de convergence des schémas et le quantifions en fonction de la régularité de Sobolev de la donnée initiale. Si nécessaire, nous régularisons la donnée initiale afin de toujours assurer les estimations de consistance. Une étape d’optimisation est alors nécessaire entre cette régularisation et l’ordre de convergence du schéma. Une seconde partie est consacrée à l’existence d’ondes progressives pour l’équation de Korteweg de Vries-Kuramoto-Sivashinsky. Par des méthodes classiques de systèmes dynamiques : système augmenté, fonction de Lyapunov, intégrale de Melnikov, par exemple, nous démontrons l’existence d’ondes oscillantes de petite amplitude. / The aim of this thesis is to study some hyperbolic-dispersive partial differential equations. A significant part is devoted to the numerical analysis and more precisely to the convergence of some finite difference schemes for the Korteweg-de Vries equation and abcd systems of Boussinesq. The numerical study follows the classical steps of consistency and stability. The main idea is to transpose at the discrete level the weak-strong stability property for hyperbolic conservation laws. We determine the convergence rate and we quantify it according to the Sobolev regularity of the initial datum. If necessary, we regularize the initial datum for the consistency estimates to be always valid. An optimization step is thus necessary between this regularization and the convergence rate of the scheme. A second part is devoted to the existence of traveling waves for the Korteweg-de Vries-Kuramoto-Sivashinsky equation. By classical methods of dynamical systems : extended systems, Lyapunov function, Melnikov integral, for instance, we prove the existence of oscillating small amplitude traveling waves.

L'introduction sur le marché financier du premier dérivé climatique portant sur les degrés-jours en 1997 aux Etats-Unis, a donné lieu à un nombre important de travaux sur la valorisation des instruments climatiques et sur la modélisation de la température moyenne journalière. Cependant, aucune étude conjointe de l'ensemble des approches d'évaluation et des représentations de la température suggérées dans la littérature n'avait été menée jusqu'à présent sur le plan empirique. Nous nous sommes donc proposés dans la présente thèse de calculer les prix des contrats à terme et des options d'achat climatiques à partir des méthodes en l'absence d'arbitrage, actuarielle et fondée sur la consommation. Nous nous sommes particulièrement intéressés au calcul des prix des contrats sur les degrés-jours des villes de Chicago, de Cincinnati et de New York pour lesquelles nous avions constaté des transactions fréquentes sur le Chicago Mercantile Exchange. L'analyse conjointe des cours estimés des contrats à terme climatiques à partir des différentes méthodologies d'évaluation a montré que le calibrage des modèles d'évaluation était nécessaire pour obtenir des prévisions de prix proches des cotations et plus particulièrement de la réalisation réelle de l'indice des degrés-jours à l'échéance.

Un algorithme permettant de discrétiser les équations aux dérivées partielles (EDP) tout en préservant leurs symétries de Lie est élaboré. Ceci est rendu possible grâce à l'utilisation de dérivées partielles discrètes se transformant comme les dérivées partielles continues sous l'action de groupes de Lie locaux. Dans les applications, beaucoup d'EDP sont invariantes sous l'action de transformations ponctuelles de Lie de dimension infinie qui font partie de ce que l'on désigne comme des pseudo-groupes de Lie. Afin d'étendre la méthode de discrétisation préservant les symétries à ces équations, une discrétisation des pseudo-groupes est proposée. Cette discrétisation a pour effet de transformer les symétries ponctuelles en symétries généralisées dans l'espace discret. Des schémas invariants sont ensuite créés pour un certain nombre d'EDP. Dans tous les cas, des tests numériques montrent que les schémas invariants approximent mieux leur équivalent continu que les différences finies standard. / An algorithm discretizing partial differential equations (PDEs) while preserving their Lie symmetries is provided. This is made possible by the use of discrete partial derivatives transforming as their continuous counterparts under the action of local Lie groups. In applications, many PDEs are invariant under the action of Lie point symmetries of infinite dimension designated as Lie pseudo-groups. To extend the invariant discretization method to such equations, a discretization of pseudo-groups is proposed. The pseudo-group action discretization transforms the continuous point symmetries into generalized symmetries in the discrete space. Invariant schemes are then created for a number of PDEs. In all cases, numerical tests demonstrate that invariant schemes are better approximations of their continuous equivalents than standard finite differences.

Ette thèse traite d’écoulements de fluides à seuil (ou viscoplastiques) en milieu confiné. Les difficultés analytiques et numériques sont dues à la multivaluation du tenseur des contraintes dans les zones plastiques ainsi qu’à la non-différentiabilité du problème de minimisation associé. Cette thèse s’articule en deux parties.Dans un premier temps, des simulations numériques parallèles très précises à l’aide d’algorithmes de dualité ont été effectuées. Elles ont permis de retrouver des résultats observés expérimentalement dont l’existence d’une ligne de glissement pour l’écoulement au dessus d’un obstacle et le caractère quasi-Poiseuille de la vitesse au-delà de cette ligne. Par ailleurs, la théorie de couche limite viscoplastique définie par Oldroyd (1947, à nombre de Bingham asymptotiquement grand) a été revisitée à nombre de Bingham modéré en milieu confiné. L’étude a mis en œuvre des allers-retours entre ces simulations et les expériences physiques de Luu et al. d’IRSTEA ainsi qu’une dérivation théorique. L’approximation de couche limite est vérifiée dans une certaine mesure à l’intérieur de la cavité. Une adaptation de la notion de couche limite viscoplastique est alors exhibée et permet d’étendre les scalings dérivés par Oldroyd (1947) et Balmforth et al. (J. of Fluid Mech, 2017). Ces scalings sont aussi généralisés au cas de la loi d’Herschel-Bulkley. Dans un second temps, on présente une analyse asymptotique des champs de vitesses et de contraintes pour des écoulements en faible épaisseur (ε). Un développement à l’ordre ε2 de la vitesse permet de trouver une équation de Reynolds à la même précision. Cette équation de Reynolds prolonge les résultats déjà existants dans le cadre newtonien, d’une part et dans le cadre fluide à seuil avec une surface libre, d’autre part. / This thesis is devoted to the flow of yield stress (or viscoplastic) fluids in pipes.Analytical and numerical difficulties lie in the multivaluation of the stress tensor in the plastic regions and in the non-differentiability of the associated minimization problem. This manuscript is organized following two main axes.First, very accurate numerical simulations were carried out using duality methods and parallel multifrontal solvers. Thus, experimental observations were recovered, namely the existence of a slip line for the flow over an obstacle and the Poiseuille-like behaviour of the velocity above this line. Moreover, the viscoplastic boundary layer theory defined by Oldroyd (1947 at high Bingham numbers) was revisited at moderate Bingham numbers in confined areas. This study provided an opportunity to go back and forth between these simulations and the physical measures of Luu et al. from IRSTEA and to perform a theoretical derivation. The boundary layer approximation is valid up to a certain extent in the cavity. An adaptation of the viscoplastic boundary layer definition is then given and allows to generalize the scalings shown by Oldroyd (1947) and Balmforth et al. (JFM 2017). These scalings are also generalized to the Herschel-Bulkley case. Then, an asymptotic analysis of the velocity and stress fields for thin layer (ε) flows is presented. A velocity development up to ε2 lets us find a Reynolds equation of same accuracy. This Reynolds equation extends the already existing results, on the one hand in the newtonian case and on the second hand for free surface flows.

Dans le monde économique, les contrats d'options sont très utilisés car ils permettent de se couvrir contre les aléas et les risques dus aux fluctuations des prix des actifs sous-jacents. La détermination du prix de ces contrats est d'une grande importance pour les investisseurs.Dans cette thèse, on s'intéresse aux problèmes d'évaluation des options, en particulier les options Européennes et Quanto sur un actif financier dont le prix est modélisé en multi dimensions par un modèle de diffusion-saut à volatilité stochastique avec sauts (1er cas considère la volatilité sans sauts, dans le 2ème cas les sauts sont pris en compte, finalement dans le 3ème cas, l'actif sous-jacent est sans saut et la volatilité suit un CEV modèle sans saut). Ce modèle permet de mieux prendre en compte certains phénomènes observés dans les marchés. Nous développons des méthodes numériques qui déterminent les valeurs des prix de ces options. On présentera d'abord le modèle qui s'écrit sous la forme d'un système d'équations intégro-différentielles stochastiques "EIDS", et on étudiera l'existence et l'unicité de la solution de ce modèle en fonction de ses coefficients, puis on établira le lien entre le calcul du prix de l'option et la résolution de l'équation Intégro-différentielle partielle (EIDP). Ce lien, qui est basé sur la notion des générateurs infinitésimaux, nous permet d'utiliser différentes méthodes numériques pour l'évaluation des options considérées. Nous introduisons alors l'équation variationnelle associée aux EIDP et démontrons qu'elle admet une unique solution dans un espace de Sobolev avec poids en s'inspirant des travaux de Zhang [106].Nous nous concentrons ensuite sur l'approximation numérique du prix de l'option en considérant le problème dans un domaine borné, et nous utilisons pour la résolution numérique la méthode des éléments finis de type (P1), et un schéma d'Euler-Maruyama, pour se servir, d'une part de la méthode de différences finies en temps, et d'autre part de la méthode de Monté Carlo et la méthode Quasi Monte Carlo. Pour cette dernière méthode nous avons utilisé les suites de Halton afin d'améliorer la vitesse de convergence.Nous présenterons une étude comparative des différents résultats numériques obtenus dans plusieurs cas différents afin d'étudier la performance et l'efficacité des méthodes utilisées. / In the modern economic world, the options contracts are used because they allow to hedge against the vagaries and risks refers to fluctuations in the prices of the underlying assets. The determination of the price of these contracts is of great importance for investors.We are interested in problems of options pricing, actually the European and Quanto options on a financial asset. The price of that asset is modeled by a multi-dimentional jump diffusion with stochastic volatility. Otherwise, the first model considers the volatility as a continuous process and the second model considers it as a jump process. Finally in the 3rd model, the underlying asset is without jump and volatility follows a model CEV without jump. This model allow better to take into account some phenomena observed in the markets. We develop numerical methods that determine the values of prices for these options. We first write the model as an integro-differential stochastic equations system "EIDS", of which we study existence and unicity of solutions. Then we relate the resolution of PIDE to the computation of the option value. This link, which is based on the notion of infinitesimal generators, allows us to use different numerical methods. We therefore introduce the variational equation associated with the PIDE, and drawing on the work of Zhang [106], we show that it admits a unique solution in a weights Sobolev space We focus on the numerical approximation of the price of the option, by treating the problem in a bounded domain. We use the finite elements method of type (P1), and the scheme of Euler-Maruyama, for this serve, on the one hand the finite differences method in time, and on the other hand the method of Monte Carlo and the Quasi Monte Carlo method. For this last method we use of Halton sequences to improve the speed of convergence.We present a comparative study of the different numerical results in many different cases in order to investigate the performance and effectiveness of the used methods.

L'objet de ce travail concerne l'application de la méthode de simulation des grandes échelles au mélange de deux fluides à propriétés thermodynamiques différentes.<br />L'origine des erreurs numériques liées à la discrétisation des équations de Navier-Stokes ainsi que leur interaction avec un modèle sous-maille sont étudiées pour une turbulence homogène et isotrope en auto-amortissement. Un code de calcul à haut pouvoir de solution est alors développé pour la simulation de couches de mélange bi-espèces. La réduction de l'amplitude des ondes acoustiques initiales est assurée par l'utilisation d'un champ en similitude temporelle.<br />L'importance relative des termes sous-maille issus des équations filtrées est mesurée à partir du filtrage explicite des champs de simulations numériques directes des couches de mélange temporelles $N_2/O_2$ et $H_2/O_2$.<br />L'utilisation d'une fermeture implicite autour d'un schéma numérique dissipatif est ensuite évoquée.

Cette thèse constitue un ensemble de travaux et de réflexions sur la question de la modélisation d’expériences en nano-optique utilisant la méthode des différences finies dans le domaine fréquentiel (FDFD) et la méthode des différences finies dans le domaine temporel (FDTD). D’abord, un code FDFD bidimensionnel, dédié au calcul de modes propres de guides d’ondes optiques, a été mis en œuvre et testé à travers une comparaison avec des résultats publiés. Dans une deuxième grande partie, nous étudions le piégeage optique de petites particules (de taille microscopique) à l’aide d’une antenne à nano-ouverture papillon (BNA) gravée à l’extrémité d’une sonde de microscope optique métallisée. Le confinement de lumière obtenue à la résonance de la nano-antenne permet un piégeage 3-D des nanoparticules de latex. Une étude systématique a été menée pour quantifier la puissance de la lumière incidente nécessaire pour un piégeage stable. Un bon accord entre les résultats expérimentaux et numériques a été obtenu dans le cas d’une BNA opérant dans l’eau à _ = 1064 nm pour le piégeage de particules de latex de 250 nm de rayon. En outre, les résultats numériques pour de plus petites particules sont présentés et montrent qu’une telle configuration est capable de piéger des particules avec des rayons aussi petits que 30 nm. Troisièmement, nous avons étudié le processus de piégeage optique basé sur l’amélioration du confinement, non seulement du champ électrique comme dans le cas de la BNA, mais aussi du magnétique que peut exhiber l’antenne métallique type diabolo (DA). Cette dernière a été récemment proposée car elle présente une résonance avec un fort confinement magnétique. Nous avons amélioré le design afin qu’une double résonance, électrique et magnétique, ait lieu au centre de la nano-antenne. Ce double confinement a ensuite été exploité pour exalter le gradient de champ au voisinage de l’antenne et ainsi aboutir à de meilleures efficacités de piégeage (moindre puissance). De plus, les résultats des simulations montrent que le processus de piégeage dépend fortement des dimensions des particules et que, pour des géométries particulières, un piégeage sans contact peut être réalisé. Cette structure doublement résonnante ouvre la voie à la conception d’une nouvelle génération de nano-pinces optiques à forte efficacit / This thesis is a set of work and reflections on modeling the experiments in nano-optics by using the finite difference method in the frequency domain (FDFD), and in time domain (FDTD). First, a two-dimensional code FDFD, dedicated to the calculation the eigenmodes of optical waveguides, has been implemented and tested through a comparison with results found in the literature. In a second large part, we study the optical trapping of small particles (of microscopic size) by using a bowtie nanoaperture antenna (BNA) engraved at the end of a metal-coated near-field optical microscope tip. The confinement of light obtained at the resonance of the nano-antenna allows 3-D trapping of latex nanoparticles. A systematic study was conducted to quantify the power of incident light necessary for stable trapping. Good agreement between the experimental and numerical results was obtained in the case of a BNA operating in water at _ = 1064 nm for the trapping of latex particles having a radius of 250 nm-radius. In addition, numerical results for smaller particles are presented and show that such configuration is capable of trapping particles with radii reaching 30 nm. Third, we studied the optical trapping process based on improved confinement of the electric field as in the case of the BNA, but also of the magnetic field, by using a metallic diabolo shape antenna (DA). This latter has been recently proposed because it exhibits resonance with a strong magnetic field confinement. We have improved the design in such a way that a double resonance, electric and magnetic, takes place in the center of the nano-antenna. This dual confinement was then used in order to enhance the field gradient in its vicinity and thus obtain better efficiencies of the trapping (less power). In addition, the simulation results show that the trapping process is greatly dependent of the particles size, and also show that, for specificl geometries, a trapping without contact can be achieved. This doubly resonant structure opens the way to the conception of a new generation of optical nano-tweezers with high efficiency.

Les besoins grandissants de prévisions robustes pour des systèmes complexes conduisent à introduire des modèles mathématiques considérant un nombre croissant de paramètres. Au temps s'ajoutent l'espace, l'aléa, les échelles de dynamiques, donnant lieu à des modèles stochastiques multi-échelles avec dépendance spatiale (modèles spatiaux). Cependant, l'explosion du temps de simulation de tels modèles complique leur utilisation. Leur analyse difficile a néanmoins permis, pour les modèles à une échelle, de développer des outils puissants: loi des grands nombres (LGN), théorème central limite (TCL), ..., puis d'en dériver des modèles simplifiés et algorithmes accélérés. Dans le processus de dérivation, des modèles et algorithmes dits hybrides ont vu le jour dans le cas multi-échelle, mais sans analyse rigoureuse préalable, soulevant ainsi la question d'approximation hybride dont la consistance constitue l'une des motivations principales de cette thèse.En 2012, Crudu, Debussche, Muller et Radulescu établissent des critères d'approximation hybride pour des modèles homogènes en espace de réseaux de régulation de gènes. Le but de cette thèse est de compléter leur travail et le généraliser à un cadre spatial.Nous avons développé et simplifié différents modèles, tous des processus de Markov de sauts pures à temps continu. La démarche met en avant, d'une part, des conditions d'approximations déterministes par des solutions d'équations d'évolution (type réaction-advection-diffusion), et, d'autre part, des conditions d'approximations hybrides par des processus stochastiques hybrides. Dans le cadre des réseaux de réactions biochimiques, un TCL est établi. Il correspond à une approximation hybride d'un modèle homogène simplifié à deux échelles de temps (suivant Crudu et al.). Puis, une LGN est obtenue pour un modèle spatial à deux échelles de temps. Ensuite, une approximation hybride est établie pour un modèle spatial à deux échelles de dynamique en temps et en espace. Enfin, des comportements asymptotiques en grandes populations et en temps long sont présentés pour un modèle d'épidémie de choléra, via une LGN suivie d'une borne supérieure pour les sous-ensembles compacts, dans le cadre d'un principe de grande déviation (PGD) correspondant.À l'avenir, il serait intéressant, entre autres, de varier la géométrie spatiale, de généraliser le TCL, de compléter les estimations du PGD, et d'explorer des systèmes complexes issus d'autres domaines. / The growing needs of precise predictions for complex systems lead to introducing stronger mathematical models, taking into account an increasing number of parameters added to time: space, stochasticity, scales of dynamics. Combining these parameters gives rise to spatial --or spatially inhomogeneous-- multiscale stochastic models. However, such models are difficult to study and their simulation is extremely time consuming, making their use not easy. Still, their analysis has allowed one to develop powerful tools for one scale models, among which are the law of large numbers (LLN) and the central limit theorem (CLT), and, afterward, to derive simpler models and accelrated algorithms. In that deduction process, the so-called hybrid models and algorithms have arisen in the multiscale case, but without any prior rigorous analysis. The question of hybrid approximation then shows up, and its consistency is a particularly important motivation of this PhD thesis.In 2012, criteria for hybrid approximations of some homogeneous regulation gene network models were established by Crudu, Debussche, Muller and Radulescu. The aim of this PhD thesis is to complete their work and generalize it afterward to a spatial framework.We have developed and simplified different models. They all are time continuous pure jump Markov processes. The approach points out the conditions allowing on the the one hand deterministic approximations by solutions of evolution equations of type reaction-advection-diffusion, and, on the other hand, hybrid approximations by hybrid stochastic processes. In the field of biochemical reaction networks, we establish a CLT. It corresponds to a hybrid approximation of a simplified homogeneous model (due to Crudu et al.). Then a LLN is obtained for a spatial model with two time scales. Afterward, a hybrid approximation is established, for a two time-space scales spatial model. Finally, the asymptotic behaviour in large population and long time are respectively presented for a model of cholera epidemic, through a LLN followed by the upper bound for compact sets, in the context of a corresponding large deviation principle (LDP).Interesting future works would be, among others, to study other spatial geometries, to generalize the CLT, to complete the LDP estimates, and to study complex systems from other fields.

L'apparition de CUDA en 2007 a rendu les GPU hautement programmables permettant ainsi aux applications scientifiques et techniques de profiter de leur capacité de calcul élevée. Des solutions ultra-rapides pour la résolution des transferts de chaleur par rayonnement et par conduction sur GPU sont présentées dans ce travail. Tout d'abord, la méthode MACZM pour le calcul des facteurs de transferts radiatifs directs en 3D et en milieu semi-transparent est représentée et validée. Ensuite, une implémentation efficace de la méthode à la base d'algorithmes de géométrie discrète et d'une parallélisation optimisée sur GPU dans CUDA atteignant 300 à 600 fois d'accélération, est présentée. Ceci est suivi par la formulation du NRPA, une version non-récursive de l'algorithme des revêtements pour le calcul des facteurs d'échange radiatifs totaux. La complexité du NRPA est inférieure à celle du PA et sont exécution sur GPU est jusqu'à 750 fois plus rapide que l'exécution du PA sur CPU. D'autre part, une implémentation efficace de la LOD sur GPU est présentée, consistant d'une alternance optimisée des solveurs et schémas de parallélisation et achevant une accélération GPU de 75 à 250 fois. Finalement, toutes les méthodes sont appliquées ensemble pour la résolution des transferts de chaleur en 3D dans un four de réchauffage sidérurgique de brames d'acier. Dans ce but, MACZM est appliquée avec un maillage multi-grille et le NRPA est appliqué au four en le découpant en zones, permettant d'avoir un temps de calcul très rapide une précision élevée. Ceci rend les méthodes utilisées de très grande importance pour la conception de stratégies de contrôle efficaces et précises. / The release of CUDA by NVIDIA in 2007 has tremendously increased GPU programmability, thus allowing scientific and engineering applications to take advantage of the high GPU compute capability. In this work, we present ultra-fast solutions for radiation and diffusion heat transfer on the GPU. First, the Multiple Absorption Coefficient Zonal Method (MACZM) for computing direct radiative exchange factors in 3D semi-transparent media is reviewed and validated. Then, an efficient implementation for MACZM is presented, based on discrete geometry algorithms, and an optimized GPU CUDA parallelization. The CUDA implementation achieves 300 to 600 times speed-up. The Non-recursive Plating Algorithm (NRPA), a non-recursive version of the plating algorithm for computing total exchange factors is then formulated. Due to low-complexity matrix multiplication algorithms, the NRPA has lower complexity than the PA does and it runs up to 750 times faster on the GPU by comparison to the CPU PA. On the other hand, an efficient GPU implementation for the Locally One Dimensional (LOD) finite difference split method for solving heat diffusion is presented, based on an optimiwed alternation between parallelization schemes and equation solvers, achieving accelerations from 75 to 250 times. Finally, all the methods are applied together for solving 3D heat transfer in a steel reheating furnace. A multi-grid approach is applied for MACZM and a zone-by zone computation for the NRPA. As a result, high precision and very fast computation time are achieved, making the methods of high interest for building precise and efficient control units.

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.