Plusieurs observations récentes semblent montrer l'existence d'effets relativistes dans l'émission des trous noirs (raie de fluorescence du Fer large, variabilité). Je présente des outils numériques pour calculer les effets relativistes (le décalage gravitationnel, la courbure de l'espace-temps) de manière cohérente dans les processus de rayonnement X et gamma autour des trous noirs. Ces codes basés sur des méthodes Monte-Carlo permettent de prendre en compte la dépendance angulaire de l'émission de la source. Les méthodes d'optimisation des codes sont aussi présentées car les temps de calcul peuvent devenir considérablement élevés. Ces outils sont appliqués pour étudier le comportement des photons et des leptons en métrique de Kerr, l'influence de la position de la source par rapport au trou noir et de l'angle d'émission, ainsi que les effets observables sur l'émission d'un disque localisé dans le plan équatorial. Le modèle axisymétrique utilisé comprend une source X en forme d'anneau située au dessus d'un disque d'accrétion autour d'un trou noir en rotation. Les spectres observés sont la superposition de la composante primaire, observée directement de la source, et de la composante réfléchie par le disque. La forme, la normalisation et la dépendance angulaire des deux composantes sont étudiées en faisant varier le rayon de l'anneau et la distance entre la source et le disque. Les effets de la relativité générale sont mis en évidence par comparaison avec les résultats obtenus en métrique Newtonienne et en relativité restreinte. Les réflexions multiples sont prises en compte, pour les photons réfléchis qui retournent vers le disque d'accrétion en raison de la courbure de l'espace, et leur influence quantifiée. Finalement, afin de confronter nos calculs aux observations, nous avons regardé la relation prédite par notre modèle, entre la composante directe et la composante réfléchie, quand la position de la source varie. Les résultats obtenus reproduisent bien les données reportées pour la source NGC4051, et donnent alors des informations sur la variation de la position de la source. J'ai donc démontré que les effets relativistes sont très importants dans la région proche du trou noir et ne peuvent être négligés dans les spectres d'émission à haute énergie. Je montre également comment les outils numériques correspondants (incluant des méthodes d'optimisation) peuvent être incorporés dans les codes de transfert de rayonnement

Les processus markoviens déterministes par morceaux (PMDM) ont été introduits dans la littérature par M.H.A. Davis en tant que classe générale de modèles stochastiques non-diffusifs. Les PMDM sont des processus hybrides caractérisés par des trajectoires déterministes entrecoupées de sauts aléatoires. Dans cette thèse, nous développons des méthodes numériques adaptées aux PMDM en nous basant sur la quantification d'une chaîne de Markov sous-jacente au PMDM. Nous abordons successivement trois problèmes : l'approximation d'espérances de fonctionnelles d'un PMDM, l'approximation des moments et de la distribution d'un temps de sortie et le problème de l'arrêt optimal partiellement observé. Dans cette dernière partie, nous abordons également la question du filtrage d'un PMDM et établissons l'équation de programmation dynamique du problème d'arrêt optimal. Nous prouvons la convergence de toutes nos méthodes (avec le plus souvent des bornes de la vitesse de convergence) et les illustrons par des exemples numériques.

La simulation numérique est devenue un outil central dans la modélisation de nombreux systèmes physiques tels que la dynamique des fluides, les plasmas, l'électromagnétisme, etc. L'existence de phénomènes multi-échelles rend l'intégration numérique de ces modèles très difficile du point de vue de la précision et du temps de calcul. En effet, dans les méthodes classiques d'intégration temporelle, le pas de temps est limité par la taille des plus petites mailles au travers d'une relation de type CFL. De plus, la forte disparité entre le pas de temps effectif et la condition CFL favorise les phénomènes de diffusion numérique. Dans la littérature, des nombreux algorithmes à pas de temps locaux (LTS) ont été développés. Pour la plupart des algorithmes LTS, les pas de temps locaux doivent être choisis parmi les fractions du pas de temps global. Nous présentons dans cette thèse une méthode asynchrone pour l'intégration explicite des équations différentielles multi-échelles. Cette méthode repose sur l'utilisation de critères de stabilité locaux, critères déterminés non pas globalement mais à partir de conditions CFL locales. De plus, contrairement aux schémas LTS, l'algorithme asynchrone permet la sélection de pas de temps indépendants pour chaque cellule de maillage. Cette thèse comporte plusieurs volets. Le premier concerne l'étude mathématique des propriétés du schéma asynchrone. Le deuxième a pour objectif d'étudier la montée en ordre, à la fois temporelle et spatiale, des méthodes asynchrones. De nombreux développements dans le cadre des méthodes de haute précision en temps ou en espace, telles que les méthodes de type Galerkin Discontinu, peuvent offrir un cadre naturel pour l'amélioration de la précision des méthodes asynchrones. Toutefois, les estimations garantissant l'ordre de précision de ces méthodes peuvent ne pas être directement compatibles avec l'aspect asynchrone. L'objectif de cette thèse est donc de développer un schéma numérique asynchrone d'ordre élevé mais qui permet également de limiter la quantité de calculs à effectuer. Le troisième volet de cette thèse se focalise sur l'application numérique puisqu'il concerne la mise en oeuvre de la méthode asynchrone dans la simulation des cas-tests représentatifs de problèmes multi-échelles. / Numerical simulation has become a central tool for the modeling of many physical systems (Fluid dynamics, plasmas, electromagnetism, etc). Multi-scale phenomena make the integration of these physical systems difficult in terms of accuracy and computational time. Numerical time-stepping integration techniques used for modeling such problems generally fall into two categories : explicit and implicit schemes. In the explicit schemes, all unknown variables are computed at the current time level from quantities already available. The time step used is then limited by the most restrictive CFL condition over the whole computation domain. In the implicit method the time step is no longer limited by the CFL conditions. However the scheme is generally not suitable for strongly coupled problems. To solve such problems, a number of local time-stepping (LTS) approaches have been developed. These methods are restricted by a local CFL condition rather than the traditional global CFL condition. For most of these LTS algorithms, local time steps are usually selected to be fractions of the global time step so that regular meeting points in time exist, and only little work is available on LTS methods with independent time steps. We present in this thesis an asynchronous method for the explicit integration of multi-scale partial differential equations. This method is restricted by a local CFL condition rather than the traditional global CFL condition. Moreover, contrary to other LTS methods, the asynchronous algorithm permits the selection of independent time steps in each mesh element. Our work consists of several components. The first one concerns the mathematical study of the properties of the asynchronous method. the objective of the second part is to study the improvement of the convergence rate for asynchronous methods. Many approaches in the context of high precision methods in time or in space, such as the Discontinuous Galerkin methods, may offer a natural setting to improve the precision of the asynchronous methods. However, the estimates ensuring the order of the accuracy of the method may not be directly compatible with the asynchronous aspect. Then, the objective is to develop a high order asynchronous numerical scheme which also preserves the computational time reduction. Finally, the third part is focused on the implementation of the asynchronous method and illustrate the advantages of the method on test-cases representative of multiscale problems.

Dans un assemblage combustible de réacteur à neutrons rapides, les aiguilles combustibles sont insérées en faisceau dans un tube hexagonal avec des jeux de montage. Lors de situations transitoires accidentelles de type séisme ou dans des phases de transport et de manutention des assemblages, les aiguilles sont donc susceptibles de se déplacer légèrement dans le faisceau et de s’entrechoquer les unes avec les autres, plusieurs centaines de milliers de contacts auront alors lieu dans le faisceau. L’objectif de cette thèse est de développer une méthode numérique capable simuler le comportement dynamique du faisceau avec jeu, afin d’évaluer les efforts de contact pour pouvoir dimensionner mécaniquement les aiguilles combustibles pour ce type de sollicitation. Le sodium n’est pas représenté, on s’intéresse donc ici au comportement de l’assemblage en air. Le comportement dynamique de l’aiguille a été étudié expérimentalement sur un banc d’essai spécialement conçu dans le cadre de cette thèse et il sera présenté dans un premier temps. L’analyse des essais de vibration libre a permis de mettre en évidence un amortissement non-linéaire dont la phénoménologie sera détaillée. La méthode numérique choisie sera ensuite décrite et comparée à des références analytiques et expérimentales, en particulier afin de justifier des critères de choix des paramètres numériques qui ont été identifiés. Enfin, on s’intéressera à la modélisation d’un assemblage combustible complet, avec différents types de représentation et à la comparaison des résultats à des essais de chocs ce qui conclura la présentation / Understanding of phenomena taking place in a structure with multiple clearances is an industrial challenge. The fuel pellets in the SFR (Sodium-cooled Fast Reactor) prototype ASTRID are enclosed in small and long pins which form a bundle inside a hexagonal assembly. The assessment of stresses in the pins during dynamic loadings is essential for the safety studies of the project. Experimental tests on the test bed CARNAC have been conducted to understand the dynamical behavior of fuel pin with their pellets. Then a numerical model has been chosen to simulate the release of an assembly against a stop. The difficulty is to simulate the dynamical behavior of a structure with a huge number of internal contacts. Numerical method has been validated on a basic problem with a reference semi-analytical method. Simplified models of the pin bundle are created to understand dynamical phenomena of a multicontact system. Then, the whole assembly is modeled. Sub-structuring to accelerate computation and a precise contact law representative of the pin to pin contacts are used. Displacements, energy and contacts force are analyzed with or without clearance. Numerical results are confronted to a previous experiment made in the CEA and provide a very good fit. The average kinetic behavior of assembly is well approximated by a beam structure, if pins are linked. Contact forces are well assess with conservatism using simplified model of a pin row

CETTE THESE RASSEMBLE DIVERSES CONTRIBUTIONS MATHEMATIQUES ET NUMERIQUES A LA CHIMIE QUANTIQUE. LE CHAPITRE 1 EST CONSACRE A UNE PRESENTATION DE L'ESPRIT ET DES MODELES DE LA CHIMIE QUANTIQUE. LE CHAPITRE 2 TRAITE DE LA CONVERGENCE D'ALGORITHMES POUR LA RESOLUTION DES EQUATIONS DE HARTREE-FOCK. LES CHAPITRES SUIVANTS PORTENT SUR DES PROBLEMES SPECIFIQUES AUX SYSTEMES MOLECULAIRES IN SITU, C'EST-A-DIRE EN INTERACTION AVEC UN ENVIRONNEMENT EXTERIEUR. UNE PREMIERE APPROCHE POUR SIMULER LES EFFETS D'ENVIRONNEMENT CONSISTE A TRAITER L'INTERACTION ENTRE LE SYSTEME MOLECULAIRE ET LE MILIEU EXTERIEUR COMME UNE PERTURBATION. AU CHAPITRE 3, ON ETEND LA THEORIE DES PERTURBATIONS DES OPERATEURS LINEAIRES AU CADRE NON LINEAIRE DU MODELE DE HARTREE-FOCK. L'INTERACTION D'UN SYSTEME MOLECULAIRE AVEC UN ENVIRONNEMENT EST SOUVENT UN PROCESSUS DYNAMIQUE. C'EST LE CAS BIEN EVIDEMMENT DES QU'ON ETUDIE UNE REACTION CHIMIQUE. LE CHAPITRE 4 CONSISTE EN L'ANALYSE MATHEMATIQUE D'UNE DES APPROXIMATIONS DE L'EQUATION DE SCHODINGER DEPENDANT DU TEMPS QUI DECRIT LA DYNAMIQUE DU SYSTEME : LE MODELE DE HARTREE-FOCK NON ADIABATIQUE. LA QUASI-TOTALITE DES REACTIONS CHIMIQUES INTERESSANT L'INDUSTRIE OU LES SCIENCES DE LA VIE SE DEROULENT EN PHASE LIQUIDE, OU LES EFFETS DE SOLVANTS JOUENT UN ROLE DETERMINANT. LES CHAPITRES 5, 6 ET 7 CONCERNENT LA RESOLUTION NUMERIQUE DES MODELES DE CONTINUUM QUI SONT LES MODELES DE SOLVATATION OFFRANT A L'HEURE ACTUELLE LE MEILLEUR COMPROMIS ENTRE QUALITE DES RESULTATS ET TEMPS DE CALCUL.

Réalisé en majeure partie sous la tutelle de feu le Professeur Paul Arminjon. Après sa disparition, le Docteur Aziz Madrane a pris la relève de la direction de mes travaux. / Après un bref rappel sur les équations hyperboliques et les méthodes numériques classiques qui serviront pour les calculs de base dans nos simulations (MCS -- Main Computation Scheme), nous revenons sur le principe de compression artificielle de A. Harten, et la construction de la méthode de compression artificielle (ACM -- Artificial Compression Method). Nous insistons dans cette partie sur la nécessité d'associer l'ACM à un détecteur de discontinuité (DoD -- Detector of Discontinuities). Le triplet MCS/DoD/ACM ainsi formé peut être optimisé module par module. Nous proposons ensuite plusieurs extensions de l'ACM à différents degrés et différents ordres, jusqu'à l'ACM(2,2). Nous formulons aussi une méthodologie de construction de ces extensions, notées ACM(d,r), pour un degré d et un ordre r. Ces extensions sont moins sensibles aux discontinuités des dérivées d'une solution continue et elles perturbent moins les calculs menés par le MCS hors des discontinuités. Ceci les rend plus modulaires et utilisables sur une plus grande variété de problèmes, même avec un DoD moins fiable. Puis, nous passons à l'extension bidimensionnelle de l'ACM(1,1). D'abord pour des maillages cartésien, nous explorons une approche utilisant la séparation en espace, et une autre utilisant une méthode amont classique, appelée DCU (Donor-Cell Upwind). Nous proposons aussi une approche originale, spécialement construite pour l'ACM, que nous qualifions de directionnelle. Les deux dernières approches (DCU et directionnelle) sont ensuite reprises pour un maillage triangulaire non-structuré, ou plus exactement sur les deux maillages duaux de Voronoï (cellules barycentriques et cellules en diamant) issus d'un maillage en triangles, avec pour MCS le schéma d'Arminjon-Viallon-Madrane. Enfin, nous abordons le problème de la détection des discontinuités et proposons un DoD basé sur une propriété physique des chocs: la production d'entropie. Les raisonnements mis en place dans cette partie font intervenir deux maillages complémentaires, pouvant être traités par deux unités de calcul différentes, tirant ainsi pleinement parti des nouvelles technologies de processeurs multi-coeurs. Cet outil est construit pour les maillages unidimensionnels ainsi que bidimensionnels cartésiens et triangulaires non-structurés. / After a short review about hyperbolic equations and the classical numerical methods that will be used as a basis for computations in our simulations (we call them MCS -- Main Computation Scheme), we come back to the principle of artificial compression, as defined by A. Harten, and the construction of the Artificial Compression Method (ACM). In this part, we insist on the necessity to associate the ACM with a Detector of Discontinuities (DoD). The trio MCS/DoD/ACM thus formed may be optimised one module at a time. Then we propose some extensions of ACM up to different degrees and different orders, ending with ACM(2,2). We also formulate a way to build these extensions, noted ACM(d,r), up to a degree d and an order r. These extensions are less sensitive to discontinuities in the derivatives of a continuous solution and they less perturbate the computations performed by the MCS outside of discontinuities. This makes them more modular and usable with a wider variety of problems, even with a less reliable DoD. We follow up with the extension in two space dimensions of the ACM(1,1). First for cartesian grids, we explore an approach using space-splitting, and another using the classical Donor-Cell Upwind scheme (DCU). We also propose an original approach, specially designed for the ACM, that we call directional. These two last approaches (DCU and directional) are then applied to triangular unstructured grids, and more precisely to the two dual Voronoï meshes (barycentric cells and diamond cells) built from the triangular elements, together with the Arminjon-Viallon-Madrane extension of the Nessyahu-Tadmor scheme as MCS. Last but not least, we turn our attention to the problem of detecting discontinuities, and propose a DoD based on a physical property of shocks: the entropy production. The ideas developped in this part imply the use of two complementary meshes, that can be treated independently by two computing units, thus taking full advantage of the new technology of multi-core processors. This tool is made available for one-dimensional meshes, as for two-dimensional cartesian and unstructured triangular grids.

Les structures nanophotoniques sont généralement simulées par des méthodes de volumes, comme la méthode des différences finies dans le domaine temporel (FDTD), ou la méthode des éléments finis (FEM). Toutefois, pour les grandes structures, ou des structures plasmoniques métalliques qui nécessitent, la mémoire et le temps de calcul requis peuvent augmenter de façon spectaculaire.Les méthodes de surface, comme la méthode des éléments de frontière (BEM) ont été développées afin de réduire le nombre d'éléments de maillage. Ces méthodes consistent à exprimer le champ formé dans tout l'espace en fonction des courants électrique et magnétique à la surface de l’objet. Combinées avec la méthode multipôle rapide (FMM) qui permet une accélération du calcul de l'interaction entre les éléments lointain du maillage, de grands systèmes peuvent ainsi être manipulés.Nous avons développé, pour la première fois à notre connaissance, une FMM sur un nouveau formalisme BEM, basé sur les potentiels scalaire et vectoriel au lieu de courants électriques et magnétiques. Cette méthode a été montrée pour permettre une simulation précise des systèmes plasmoniques métalliques, tout en offrant une réduction significative des besoins de calcul. Des systèmes nanophotoniques complexes ont été simulés, comme une lentille plasmonique composé d'un ensemble de nanotubes d'or. / Nanophotonic structures are generally simulated by volume methods, as Finite-difference time-domain (FDTD) method, or Finite element method (FEM). However, for large structures, or metallic plasmonic structures, the memory and time computation required can increase dramatically, and make proper simulation infeasible.Surface methods, like the boundary element method (BEM) have been developed to reduce the number of mesh elements. These methods consist in expressing the electromagnetic filed in whole space as a function of electric and magnetic currents at the surface of scatterers. Combined with the fast multipole method (FMM) that enables a huge acceleration of the calculation of interaction between far mesh elements, very large systems can thus be handled.What we performed is the development of an FMM on a new BEM formalism, based on scalar and vector potentials instead of electric and magnetic currents, for the first time to our knowledge. This method was shown to enable accurate simulation of metallic plasmonic systems, while providing a significant reduction of computation requirements, compared to BEM-alone. Several thousands of unknowns could be handled on a standard computer. More complex nanophotonic systems have been simulated, such as a plasmonic lens consisting of a collection of gold nanorods.

Dans la première partie, on fait une approche du verrouillage en élasticité non linéaire. Nous étudions déplacements d'une coque cylindrique peu profonde de longueur infinie. Sous des hypothèses raisonnables on peut supposer que le problème est modélisé par une arche circulaire soumise à une pression F. On montre comment calculer analytiquement, toutes ses solutions, moyennant une paramétrisation des solutions. On donne deux critères pour mettre en évidence le phénomène de verrouillage en élasticité non linéaire ; le premier est basé sur l'utilisation d'une abscisse curviligne sur les branches de solutions, et le second sur la résolution d'un problème de minimisation en utilisant l'inverse de Moore-Penrose. Après avoir retenu le critère de Moore-Penrose, pour caractériser le verrouillage en non linéaire, nous avons approché le problème par une méthode d'éléments finis, combinée à la méthode de Moore-Penrose. Ceci nous a permis d'établir différents diagrammes de convergence montrant le phénomène de verrouillage et l'efficacité des méthodes à le traiter. Dans la deuxième partie, on étudie le flambage d'une coque cylindrique verticale, soumise à un champ de forces volumiques de densité f et en sa partie supérieure, soumise à une densité de forces surfaciques λ. Nous utilisons un modèle de coques minces élastiques, en théorie non linéaire, en se référant aux travaux de Koiter (1966). On s'intéresse uniquement aux solutions axisymétriques de ce modèle. On montre l'existence de solution pour toute valeur des charges et on donne une condition suffisante d'unicité. On fait une étude détaillée de la bifurcation de ce problème. Cette étude qualitative est complétée par une étude quantitative

Les solveurs de calcul en mécanique des fluides numérique (solveurs CFD) ont atteint leur maturité en termes de précision et d'efficacité de calcul. Toutefois, des progrès restent à faire pour les écoulements instationnaires surtout lorsqu'ils sont régis par de grandes structures cohérentes. Pour ces écoulements, les solveurs CFD actuels n'apportent pas de solutions assez précises à moins d'utiliser des maillages très fins. De plus, la haute précision est une caractéristique cruciale pour l'application des stratégies avancées de simulation de turbulence, comme la Simulation des Grandes Echelles (LES). Afin d'appliquer les méthodes d'ordre élevé pour les écoulements instationnaires complexes plusieurs points doivent être abordés dont la robustesse numérique et la capacité à gérer des géométries complexes.Dans cette thèse, nous étudions une famille d'approximations compactes qui offrent une grande précision non pour chaque dérivée spatiale traitée séparement mais pour le résidu r complet, c'est à dire la somme de tous les termes des équations considérées. Pour des problèmes stationnaires résolus par avancement temporelle, r est le résidu à l'état stationnaire ne comprenant que des dérivées spatiales; pour des problèmes instationnaires r comprend également la dérivée temporelle. Ce type de schémas sont appelés schémas Compacts Basés sur le Résidu (RBC). Plus précisément, nous développons des schémas RBC d'ordre élevé pour des écoulements instationnaires compressibles, et menons une étude approfondie de leurs propriétés de dissipation. Nous analysons ensuite les erreurs de dissipation et la dispersion introduites par les schémas RBC afin de quantifier leur capacité à résoudre une longueur d'onde donnée en utilisant un nombre minimal de points de maillage. Les capacités de la dissipation de RBC à drainer seulement l'énergie aux petites échelles sous-résolues sont également examinées en vue de l'application des schémas RBC pour des simulations LES implicites (ILES). Enfin, les schémas RBC sont étendus à la formulation de type volumes finis (FV) afin de gérer des géométries complexes. Une formulation FV des schémas RBC d'ordre trois préservant une précision d'ordre élevé sur des maillages irréguliers est présentée et analysée. Des applications numériques, dont la simulation d'écoulements instationnaires complexes de turbomachines régis par les équations de Navier-Stokes moyennées et des simulations ILES d'écoulements turbulents dominés par des structures cohérentes dynamiques ou en décroissance, confirment les résultats théoriques. / Computational Fluid Dynamics (CFD) solvers have reached maturity in terms of solution accuracy as well as computational efficiency. However, progress remains to be done for unsteady flows especially when governed by large, coherent structures. For these flows, current CFD solvers do not provide accurate solutions unless very fine mesh are used. Moreover, high-accuracy is a crucial feature for the application of advanced turbulence simulation strategies, like Large Eddy Simulation (LES). In order to apply high-order methods to complex unsteady flows several issues needs to be addressed among which numerical robustness and the capability of handling complex geometries.In the present work, we study a family of compact approximations that provide high accuracy not for each space derivative treated apart but for the complete residual r, i.e. the sum of all of the terms in the governing equations. For steady problems solved by time marching, r is the residual at steady state and it involves space derivatives only; for unsteady problems, r also includes the time derivative. Schemes of this type are referred-to as Residual-Based Compact (RBC). Precisely, we design high-order finite difference RBC schemes for unsteady compressible flows, and provide a comprehensive study of their dissipation properties. The dissipation and dispersion errors introduced by RBC schemes are investigated to quantify their capability of resolving a given wave length using a minimal number of grid-points. The capabilities of RBC dissipation to drain energy only at small, ill-resolved scales are also discussed in view of the application of RBC schemes to implicit LES (ILES) simulations. Finally, RBC schemes are extended to the Finite Volume (FV) framework in order to handle complex geometries. A high-order accuracy preserving FV formulation of the third-order RBC scheme for general irregular grids is presented and analysed. Numerical applications, including complex Reynolds-Averaged Navier-Stokes unsteady simulation of turbomachinery flows and ILES simulations of turbulent flows dominated by coherent structure dynamics or decay, support the theoretical results.

L'objectif de cette thèse est l'étude de la reconnexion magnétique dans les plasmas non-collisionels et relativistes. De tels plasmas sont présents dans divers objets astrophysiques (MQs, AGNs, GRBs...), où la reconnexion pourrait expliquer la production de particules et de radiation de haute énergie, un chauffage, ou des jets. Une compréhension fondamentale de la reconnexion n'est cependant toujours pas acquise, en particulier dans les plasmas relativistes ion-électron. Nous présentons d'abord les bases de la reconnexion magnétique. Nous démontrons des résultats particuliers à la physique des plasmas relativistes, concernant par exemple la distribution de Maxwell-Jüttner. Ensuite, nous réalisons une étude détaillée de l'outil numérique utilisé : les simulations particle-in-cell (PIC). Le fait que le plasma réel contienne beaucoup plus de particules que le plasma PIC a des conséquences importantes (collisionalité, relaxation, bruit) que nous décrivons. Enfin, nous étudions la reconnexion magnétique dans les plasmas ion-électron et relativistes à l'aide de simulations PIC. Nous soulignons des points spécifiques : loi d'Ohm (l'inertie de bulk dominante), zone de diffusion, taux de reconnexion (et sa normalisation relativiste). Les ions et les électrons produisent des lois de puissance, avec un index qui dépend de la vitesse d'Alfvén et de la magnétisation, et qui peut être plus dur que dans le cas des chocs non-collisionels. De plus, les ions peuvent avoir plus ou moins d'énergie que les électrons selon la valeur du champ guide. Ces résultats fournissent une base solide à des modèles d'objets astrophysiques qui, jusque là, supposaient a priori ces résultats. / The purpose of this thesis is to study magnetic reconnection in collisionless and relativistic plasmas. Such plasmas can be encountered in various astrophysical objects (microquasars, AGNs, GRBs...), where reconnection could explain high-energy particle and photon production, plasma heating, or transient large-scale outflows. However, a first principle understanding of reconnection is still lacking, especially in relativistic ion-electron plasmas. We first present the basis of reconnection physics. We derive results relevant to relativistic plasma physics, including properties of the Maxwell-Jüttner distribution. Then, we provide a detailed study of our numerical tool, particle-in-cell simulations (PIC). The fact that the real plasma contains far less particles than the PIC plasma has important consequences concerning relaxation times or noise, that we describe. Finally, we study relativistic reconnection in ion-electron plasmas with PIC simulations. We stress outstanding properties: Ohm's law (dominated by bulk inertia), structure of the diffusion zone, energy content of the outflows (thermally dominated), reconnection rate (and its relativistic normalization). Ions and electrons produce power law distributions, with indexes that depend on the inflow Alfvén speed and on the magnetization of the corresponding species. They can be harder than those produced by collisionless shocks. Also, ions can get more or less energy than the electrons, depending on the guide field strength. These results provide a solid ground for astrophysical models that, up to now, assumed with no prior justification the existence of such distributions or of such ion/electron energy repartition.