Cette thèse porte sur l'étude des points entiers et rationnels de certaines courbes et variétés modulaires. Après une brève introduction décrivant les motivations et le cadre de ce genre d'études ainsi que les résultats principaux de la thèse, le manuscrit se divise en trois parties. Le premier chapitre s'intéresse aux Q-courbes, et aux morphismes Gal(Q/Q) -> PGL2(Fp) qu'on peut leur associer pour tout p premier. Nous montrons que sous de bonnes hypothèses, pour p assez grand par rapport au discriminant du corps de définition de la Q-courbe, ce morphisme est surjectif, ce qui résout un cas particulier du problème d'uniformité de Serre (toujours ouvert en général). Les outils principaux du chapitre sont la méthode de Mazur (basée ici sur des résultats d'Ellenberg), la méthode de Runge et des théorèmes d'isogénie, suivant la structure de preuve de Bilu et Parent. Le second chapitre consiste en des estimations analytiques de sommes pondérées de valeurs de fonctions L de formes modulaires, dans l'esprit de techniques développées par Duke et Ellenberg. La motivation de départ d'un tel résultat est l'application de la méthode de Mazur dans le premier chapitre. Le troisième chapitre est consacré à la recherche de généralisations de la méthode de Runge pour des variétés de dimension supérieure. Nous y redémontrons un résultat de Levin inspiré de cette méthode, avant d'en prouver une forme assouplie dite "de Runge tubulaire", plus largement applicable. Dans l'optique de recherche de points entiers de variétés modulaires, nous en donnons enfin un exemple d'utilisation à la réduction d'une surface abélienne en produit de courbes elliptiques. / This thesis concerns the study of integral and rational points on some modular curves and varieties. After a brief introduction which describes the motivation and the setting of this topic as well as the main results of this thesis, the manuscript follows a threefold development. The first chapter focuses on Q-curves, and on the morphisms Gal(Q/Q) -> PGL2(Fp) that we can build with a Q-curve for every prime p. We prove that, under good hypotheses, for p large enough with respect to the discriminant of the definition field of the Q-curve, such a morphism is surjective, which solves a particular case of Serre's uniformity problem (still open in general). The main tools of the chapter are Mazur's method (based here on results of Ellenberg), Runge's method, and isogeny theorems, following the strategy of Bilu and Parent. The second chapter covers analytic estimates of weighted sums of L-function values of modular forms, in the fashion of techniques designed by Duke and Ellenberg. The initial goal of such a result is the application of Mazur's method in the first chapter. The third chapter is devoted to the search for generalisations of Runge's method for higherdimensional varieties. Here we prove anew a result of Levin inspired by this method, before proving an enhanced version called "tubular Runge", more generally applicable. In the perspective of studying integral points of modular varieties, we finally give an example of application of this theorem to the reduction of an abelian surface in a product of elliptic curves.

Cette thèse étudie le problème de minoration de la hauteur canonique sur les courbeselliptiques. Son résultat diophantien principal utilise des méthodes d’intersectionarithmétique pour retrouver un résultat de Laurent, qui démontrait la conjecturede Lehmer pour les courbes elliptiques à multiplications complexes à un exposant" près, tout en explicitant complètement sa dépendance en divers paramètres liésà la courbe elliptique ; une telle démarche peut être motivée par la conjecture deLang, qui présage une minoration possible de la hauteur canonique proportionnelle,essentiellement, à la hauteur de Faltings de la courbe.Notre dissertation commence toutefois par une partie dédiée à l’explicitation duthéorème de densité de Chebotarev, qui reprend les grandes lignes d’un travail deLagarias et Odlyzko, et s’avère être cruciale dans notre approche du problème deLehmer elliptique. On obtient également des majorations des zéros de Siegel et de lanorme du plus petit idéal premier entrant en jeu dans le théorème de Chebotarev. / In this thesis we consider the problem of lower bounds for the canonical height onelliptic curves, aiming for the conjecture of Lehmer. Our main diophantine result isan explicit version of a theorem of Laurent (who proved this conjecture for ellipticcurves with CM up to a " exponent) using arithmetic intersection, enlightening thedependence with parameters linked to the elliptic curve ; such a result can be motivatedby the conjecture of Lang, hoping for a lower bound proportional to, roughly,the Faltings height of the curve.Nevertheless, our dissertation begins with a part dedicated to a completely explicitversion of the density theorem of Chebotarev, along the lines of a previous workdue to Lagarias and Odlyzko, which will be crucial to investigate the elliptic Lehmerproblem. We also obtain upper bounds for Siegel zeros, and for the smallest primeideal whose Frobenius is in a fixed conjugacy class.

La reconstruction en tomographie discrète de la microstructure de l’os joue un role très important pour le diagnostic de l’ostéoporse, une maladie des os très fréquente. Le diagnostic clinique est basé sur l’absortiométrie duale de rayons X. Avec la tomographie de rayons X, une résolution spatiale élevée avec des images reconstruites in vivo requiert une dose d’irradiation élevée et un temps de balayage long, ce qui est dangereux pour le patient. Une des méthodes pour résoudre ce problème est de limiter le nombre de projections. Cependant, avec cette méthode le problème de reconstruction devient mal posé. Deux types de régularisation par Variation Totale minimisées avec la méthode Alternate Direction of Minimization Method (ADMM) et deux schémas basés sur les méthodes de régularisation Level-set sont appliquées à deux images d’os expérimentales acquises avec un synchrotron (pixel size: 15 μm). Des images de tailles variées et avec différents niveaux de bruit Gaussien additifs ajoutés aux projections sont utlisées pour étudier l’efficacité des méthodes de régularisation. Des minima locaux sont obtenus avec ces méthodes déterministes. Une approche globale d’optimisation est nécessaire pour améliorer les résultats. Des perturbations stochastiques peuvent être un moyen très utile pour échapper aux minima locaux. Dans une première approche, une équation différentielle stochastique basée sur la régularisation level-set est étudiée. Cette méthode améliore les résultats de reconstruction mais ne modifie que les frontières entre les régions 0 et 1. Ensuite une équation aux dérivées partielles stochastique est obtenue avec la régularisation TV pour améliorer la méthode stochastique level-set. A la fin de notre travail, nous avons étendu la méthode de régularisation à des images 3D avec des données réelles. Cette algorithme a été implémenté avec RTK. Nous avons aussi étendu l’approche level-set utilisée pour la tomographie binaire au cas multi-level. / Discrete tomography reconstruction of bone microstructure is important in diagnosis of osteoporosis. One way to reduce the radiation dose and scanning time in CT imaging is to limit the number of projections. This method makes the reconstruction problem highly ill-posed. A common solution is to reconstruct only a finite number of intensity levels. In this work, we investigate only binary tomography reconstruction problem. First, we consider variational regularization methods. Two types of Total Variation (TV) regularization approaches minimized with the Alternate Direction of Minimization Method (ADMM) and two schemes based on Level-set (LS) regularization methods are applied to two experimental bone cross-section images acquired with synchrotron micro-CT. The numerical experiments have shown that good reconstruction results were obtained with TV regularization methods and that level-set regularization outperforms the TV regularization for large bone image with complex structures. Yet, for both methods, some reconstruction errors are still located on the boundaries and some regions are lost when the projection number is low. Local minima were obtained with these deterministic methods. Stochastic perturbations is a useful way to escape the local minima. As a first approach, a stochastic differential equation based on level-set regularization was studied. This method improves the reconstruction results but only modifies the boundaries between the 0 and 1 regions. Then partial stochastic differential equation obtained with the TV regularization semi-norm were studied to improve the stochastic level-set method. The random change of the boundary are performed in a new way with the gradient or wavelet decomposition of the reconstructed image. Random topological changes are included to find the lost regions in the reconstructed images. At the end of our work, we extended the TV regularization method to 3D images with real data on RTK (Reconstruction Toolkit). And we also extended the level-set to the multi-level cases.

La caméra Compton est un dispositif permettant d’imager les sources de rayonnement gamma. Ses avantages sont sa sensibilité (absence de collimateur mécanique) et la possibilité de reconstruire des images 3D avec un dispositif immobile. Elle également adaptée pour des sources à large spectre énergétique. Ce dispositif est un candidat prometteur en médecine nucléaire et en hadronthérapie. Ces travaux, financés par le projet européen ENVISION (European NoVel Imaging Systems for ION therapy) Coopération-FP7, portent sur le développement de méthodes de reconstruction d’images pour la caméra Compton pour la surveillance de la thérapie par ions. Celle-ci nécessite idéalement une reconstruction temps réel avec une précision millimétrique, même si le nombre de données acquises est relativement faible. Nous avons développé des méthodes analytiques et itératives. Leurs performances sont analysées dans le contexte d’acquisitions réalistes (géométrie de la caméra, nombre d’événements). Nous avons développé une méthode analytique de rétroprojection filtrée. Cette méthode est rapide mais nécessite beaucoup de données. Nous avons également développé des méthodes itératives utilisant un algorithme de maximisation de la fonction de vraisemblance. Nous avons proposé un modèle probabiliste pour l’estimation des éléments de la matrice système nécessaire à la reconstruction et nous avons développé différentes approches pour le calcul de ses éléments : l’une néglige les incertitudes de mesure sur l’énergie, l’autre les prend en compte en utilisant une distribution gaussienne. Nous avons étudié une méthode simplifiée utilisant notre modèle probabiliste. Plusieurs reconstructions sont menées à partir de données simulées, obtenues avec Geant4, mais provenant aussi de plusieurs prototypes simulés de caméra Compton proposés par l’Institut de Physique Nucléaire de Lyon (IPNL) et par le Centre de recherche de Dresde-Rossendorf en Allemagne. Les résultats sont prometteurs et des études plus poussées, à partir de données encore plus réalistes, viseront à les confirmer. / The Compton camera is a device for imaging gamma radiation sources. The advantages of the system lie in its sensitivity, due to the absence of mechanical collimator, and the possibility of imaging wide energy spectrum sources. These advantages make it a promising candidate for application in hadrontherapy. Funded by the european project ENVISION, FP7-Cooperation Work Program, this work deals with the development of image reconstruction methods for the Compton camera. We developed both analytical and iterative methods in order to reconstruct the source from cone-surface projections. Their performances are analyzed with regards to the context (geometry of the camera, number of events). We developped an analytical method using a Filtered BackProjection (FBP) formulation. This method is fast but really sensitive to the noise. We have also developped iterative methods using a List Mode-Maximum Likelihood Expectation Maximization (LM-MLEM) algorithm. We proposed a new probabilistic model for the computation of the elements of the system matrix and different approaches for the calculation of these elements neglecting or not the measurement uncertainties. We also implemented a simplified method using the probabilistic model we proposed. The novelty of the method also lies on the specific discretization of the cone-surface projections. Several studies are carried out upon the reconstructions of simulated data worked out with Geant4, but also simulated data obtained from several prototypes of Compton cameras under study at the Institut de Physique Nucléaire de Lyon (IPNL) and at the Research Center of Dresden-Rossendorf. Results are promising, and further investigations on more realistic data are to be done.

Les terres agricoles ont été durablement contaminées à la fois par les pesticides mis à la disposition des agriculteurs pour lutter contre les charançons et autres insectes nuisibles, et par les engrais azotées pour augmenter la productivité chez les plantes.Des recherches récentes concernent des cultures alternatives écologiques utilisant les plantes de service qui fournissent les nutriments aux plantes principales. Ce travail de thèse s'inscrit dans cette perspective, d'un point de vue modélisation.L'accent est mis sur la résolution de problèmes de contrôle du phénomène d'absorption de nutriments, par les racines dans la rhizosphère (partie proche des racines), en considérant les deux cas de sols : sol sain et sol pollué.Ces phénomènes d'absorption sont modélisés par des systèmes d'advection-diffusion de type Nye-Tinker-Barber (NTB). La concentration de nutriments absorbée, solution du problème, est une fonction du temps et de l'espace.On étudie l'existence de solution du système NTB dans les deux cas où la fonction d'absorption de nutriments à la frontière (surface de la racine) appelée fonction de Michealis-Menten, est linéaire et/ou non linéaire, à l’aide des outils d’analyse fonctionnelle. On étudie ensuite les problèmes de contrôle optimal associés au système NTB, en considérant les deux cas linéaire et non linéaire, en application pour les deux cas d’absorption de nutriments en sol non pollué puis en sol pollué. Pour le premier cas, on utilise les techniques classiques de recherche d'un contrôle pour les systèmes distribués, tandis que, pour le second cas, on fait appel aux notions de contrôle sans regret et contrôle à moindres regrets de J.-L. Lions. Les contrôles obtenus pour les différents problèmes sont caractérisés chacun par un système d'optimalité (SO) cas sans pollution, et système d’optimalité singulier (SOS) dans le cas avec pollution.= / Agriculture soils were highly contaminated for a long time by pesticides which were widely used by producers to fight against weevils. Soils where also contaminated by the use of fertilizers to increase the plant development. An ecological alternative using service plants is encouraged following recent research. The aim of this work is to give a mathematical and a modelling point of view as we study the mecha- nisms of nutrient transfer to plants using the mathematical analysis and optimal control theories. The two cases of polluted and non-polluted soils are considered. The nutrient transfer and uptake processes are modeled by an advection-diffusion system derived from the Nye-Tinker-Barber (NTB) model. The absorbed nutrient concentration represented by the Michaelis-Mention function at the root surface of the principal plant, depends on time and space. We study the existence of a solution for the linear and nonlinear NTB systems, then we characterize the opti- mal control which corresponds to the added nutrients from the service plant. For the pollution case, we use the concept of low-regret and no-regret control of J.-L. Lions.

Nous étudions l’analyse asymptotique en relativité générale sous deux aspects: le peeling et le scattering (diffusion) conforme. Le peeling est construit pour les champs scalaires linéaire et non-linéaires et pour les champs de Dirac en espace-temps de Kerr (qui est non-stationnaire et à symétrie simplement axiale), généralisant les travaux de L. Mason et J-P. Nicolas (2009, 2012). La méthode des champs de vecteurs (estimations d’énergie géométriques) et la technique de compactification conforme sont développées. Elles nous permettent de formuler les définitions du peeling à tous ordres et d’obtenir les données initiales optimales qui assurent ces comportements. Une théorie de la diffusion conforme pour les équations de champs sans masse de spîn n/2 dans l’espace-temps de Minkowski est construite.En effectuant les compactifications conformes (complète et partielle), l’espace-temps est complété en ajoutant une frontière constituée de deux hypersurfaces isotropes représentant respectivement les points limites passés et futurs des géodésiques de type lumière. Le comportement asymptotique des champs s’obtient en résolvant le problème de Cauchy pour l’équation rééchelonnée et en considérant les traces des solutions sur ces bords. L’inversibilité des opérateurs de trace, qui associent le comportement asymptotique passé ou futur aux données initiales, s’obtient en résolvant le problème de Goursat sur le bord conforme. L’opérateur de diffusion conforme est alors obtenu par composition de l’opérateur de trace futur avec l’inverse de l’opérateur de trace passé. / This work explores two aspects of asymptotic analysis in general relativity: peeling and conformal scattering.On the one hand, the peeling is constructed for linear and nonlinear scalar fields as well as Dirac fields on Kerr spacetime, which is non-stationary and merely axially symmetric. This generalizes the work of L. Mason and J-P. Nicolas (2009, 2012). The vector field method (geometric energy estimates) and the conformal technique are developed. They allow us to formulate the definition of the peeling at all orders and to obtain the optimal space of initial data which guarantees these behaviours. On the other hand, a conformal scattering theory for the spin-n/2 zero rest-mass equations on Minkowski spacetime is constructed. Using the conformal compactifications (full and partial), the spacetime is completed with two null hypersurfaces representing respectively the past and future end points of null geodesics. The asymptotic behaviour of fields is then obtained by solving the Cauchy problem for the rescaled equation and considering the traces of the solutions on these hypersurfaces. The invertibility of the trace operators, that to the initial data associate the future or past asymptotic behaviours, is obtained by solving the Goursat problem on the conformal boundary. The conformal scattering operator is then obtained by composing the future trace operator with the inverse of the past trace operator.

La modélisation du transport réactif du contaminant en milieu poreux est un problème complexe cumulant les difficultés de la modélisation du transport avec celles de la modélisation de la chimie et surtout du couplage entre les deux. Cette modélisation conduit à un système d'équations aux dérivées partielles et algébriques dont les inconnues sont les quantités d'espèces chimiques. Une approche possible, déjà utilisée par ailleurs, est de choisir la méthode globale DAE : l'utilisation d'une méthode de lignes, correspondant à la discrétisation en espace seulement, conduit à un système différentiel algébrique (DAE) qui doit être résolu par un solveur adapté. Dans notre cas, on utilise le solveur IDA de Sundials qui s'appuie sur une méthode implicite, à ordre et pas variables, et qui requiert à chaque pas de temps la résolution d'un grand système non linéaire associé à une matrice jacobienne. Cette méthode est implémentée dans un logiciel qui s'appelle GRT3D (Transport Réactif Global en 3D). Le présent travail présente une amélioration de la méthode GDAE, du point de vue de la performance, de la stabilité et de la robustesse. Nous avons ainsi enrichi les possibilités de GRT3D, par la prise en compte complète des équations de précipitation-dissolution permettant l'apparition ou la disparition d'une espèce précipitée. En complément de l'étude de la méthode GDAE, nous présentons aussi une méthode séquentielle non itérative (SNIA), qui est une méthode basée sur le schéma d'Euler explicite : à chaque pas de temps, on résout explicitement l'équation de transport et on utilise ces calculs comme données pour le système chimique, résolu dans chaque maille de façon indépendante. Nous présentons aussi une comparaison entre cette méthode et l'approche GDAE. Des résultats numériques pour deux cas tests, celui proposé par l'ANDRA (cas-test 2D) d'une part, celui proposé par le groupe MoMas (Benchmark "easy case") d'autre part, sont enfin présentés, commentés et analysés. / Modeling reactive transport of contaminants in porous media is a complex problem combining the difficulties of modeling the trasport with those of modeling the chemistry and especially the coupling between the two .This model leads to a system of partial differential equations and algebraic equations whose unknowns are the quantities of chemical species. One approach , already used elsewhere , is choosing the global DAE method : using the method of lines, discretization in space only, leads to a differential algebraic system (DAE ) to be solved by a suitable solver . In our case , the solver IDA Sundials relies on an implicit method, order is used but not variables, and requires at each time solving a large nonlinear system associated with a Jacobian matrix . This method is implemented in a software called GRT3D (Global Reactive Transport in 3D). This paper presents an improved GDAE method , from the standpoint of performance, the stability and robustness. We have enriched the possibilities of GRT3D , by taking full account of the equations of dissolution – precipitation for the appearance or disappearance of precipitated species. In addition to the study of the GDAE method, we also present a non-iterative sequential method ( SNIA ) which is a method based on the explicit Euler scheme : at each time step, we explicitly solve the transport equation and we use these calculations as data for the chemical system which is resolved in each cell independently. We also present a comparison between this method and GDAE approach . Numerical results for two test cases , one proposed by ANDRA ( 2D test case ) on one hand and one proposed by the group MOMAS ( Benchmark "easy case" ) on the other hand, are finally presented , discussed and analyzed.

Dans cette thèse, nous présentons des résultats d’homogénéisation périodique d’un problème spectral et de l’équation d’ondes de Bloch. Il permet de modéliser les ondes à basse et haute fréquences. La partie modèle à basse fréquence est bien connu et n’est pas donc abordée dans ce travail. A contrario ; la partie à haute fréquence du modèle, qui représente des oscillations aux échelles microscopiques et macroscopiques, est un problème laissé ouvert. En particulier, les conditions aux limites de l’équation macroscopique à hautes fréquences établies dans [36] n’étaient pas connues avant le début de la thèse. Ce dernier travail apporte trois contributions principales. Les deux premières contributions, portent sur le comportement asymptotique du problème d’homogénéisation périodique du problème spectral et de l’équation des ondes en une dimension. La troisième contribution consiste en une extension du modèle du problème spectral posé dans une bande bi dimensionnelle et bornée. Le résultat d’homogénéisation comprend des effets de couche limite qui se produisent dans les conditions aux limites de l’équation macroscopique à haute fréquence. / In this dissertation, we present the periodic homogenization of a spectral problem and the waveequation with periodic rapidly varying coefficients in a bounded domain. The asymptotic behavioris addressed based on a method of Bloch wave homogenization. It allows modeling both the lowand high frequency waves. The low frequency part is well-known and it is not a new point here.In the opposite, the high frequency part of the model, which represents oscillations occurringat the microscopic and macroscopic scales, was not well understood. Especially, the boundaryconditions of the high-frequency macroscopic equation established in [36] were not known prior to thecommencement of thesis. The latter brings three main contributions. The first two contributions, areabout the asymptotic behavior of the periodic homogenization of the spectral problem and waveequation in one-dimension. The third contribution consists in an extension of the model for thespectral problem to a thin two-dimensional bounded strip Ω = (0; _) _ (0; ") _ R2. The homogenizationresult includes boundary layer effects occurring in the boundary conditions of the high-frequencymacroscopic equation.

L’objectif de ce mémoire est d’introduire une nouvelle méthode smoothed particle hydrodynamics (SPH) implicite purement lagrangienne, pour la résolution des équations de Navier- Stokes incompressibles bidimensionnelles en présence d’une surface libre. Notre schéma de discrétisation est basé sur celui de Kéou Noutcheuwa et Owens [19]. Nous avons traité la surface libre en combinant la méthode multiple boundary tangent (MBT) de Yildiz et al. [43] et les conditions aux limites sur les champs auxiliaires de Yang et Prosperetti [42]. Ce faisant, nous obtenons un schéma de discrétisation d’ordre $\mathcal{O}(\Delta t ^2)$ et $\mathcal{O}(\Delta x ^2)$, selon certaines contraintes sur la longueur de lissage $h$. Dans un premier temps, nous avons testé notre schéma avec un écoulement de Poiseuille bidimensionnel à l’aide duquel nous analysons l’erreur de discrétisation de la méthode SPH. Ensuite, nous avons tenté de simuler un problème d’extrusion newtonien bidimensionnel. Malheureusement, bien que le comportement de la surface libre soit satisfaisant, nous avons rencontré des problèmes numériques sur la singularité à la sortie du moule. / The objective of this thesis is to introduce a new implicit purely lagrangian smoothed particle hydrodynamics (SPH) method, for the resolution of the two-dimensional incompressible Navier-Stokes equations in the presence of a free surface. Our discretization scheme is based on that of Kéou Noutcheuwa et Owens [19]. We have treated the free surface by combining Yildiz et al. [43] multiple boundary tangent (MBT) method and boundary conditions on the auxiliary fields of Yang et Prosperetti [42]. In this way, we obtain a discretization scheme of order $\mathcal{O}(\Delta t ^2)$ and $\mathcal{O}(\Delta x ^2)$, according to certain constraints on the smoothing length $h$. First, we tested our scheme with a two-dimensional Poiseuille flow by means of which we analyze the discretization error of the SPH method. Then, we tried to simulate a two-dimensional Newtonian extrusion problem. Unfortunately, although the behavior of the free surface is satisfactory, we have encountered numerical problems on the singularity at the output of the die.

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.