Simulation de la cinétique d’absorption des défauts ponctuels par les dislocations et amas de défauts

Cluster Dynamics (CD) is a mean field simulation method which makes it possible to predict the materials evolution under irradiation. In this work, we focus on sink strengths, which are the parameters used in CD to represent the capacity of sinks (dislocations, cavities…) to absorb point defects (vacancies, self-interstitials). To calculate the sink strengths, object kinetic Monte Carlo (OKMC) simulations are performed. The energy of point defects at stable and saddle points is described through their elastic dipoles. These elements are computed in an aluminum crystal using density functional theory. In a first part, the sink strengths of the main objects found in irradiated microstructures (dislocations, cavities and dislocation loops) are calculated on simple configurations. This study reveals the importance of elastic interactions, and enables us to identify the saddle point anisotropy of point defects as an important parameter, as it modifies both the point defects trajectories and the sink strength values. Then, we focus on the role of the relative position of sinks in their capacity to absorb point defects. Microstructures containing a large number of dislocation loops are generated by OKMC simulations and the absorption of point defects by those microstructures is measured. It is shown that the neighborhood of a sink modifies noticeably its capacity to absorb the point defects and this behavior is rationalized through the Voronoi volume associated with each sink. This study leads to the proposal of a new sink strength expression, and of a new formalism for CD in which clusters are described by their size and their Voronoi volume. The results show that this formalism makes it possible to strongly improve the prediction of the evolution of cluster size distributions during irradiation.La dynamique d'amas (DA) est une méthode de simulation en champ moyen permettant la prédiction de l'évolution des matériaux sous irradiation. Dans cette thèse, on s'intéresse aux forces de puits, qui sont les paramètres utilisés en DA pour représenter la capacité des puits (dislocations, cavités...) à absorber les défauts ponctuels (lacunes, interstitiels). Pour calculer les forces de puits, on réalise des simulations Monte Carlo cinétique sur objets (OKMC) dans lesquelles l'énergie des défauts ponctuels au point stable et au point col est décrite à l'aide de dipôles élastiques. Ces quantités sont préalablement calculées pour un cristal d'aluminium en utilisant la théorie de la fonctionnelle de la densité. Une première partie de ce travail vise à évaluer les forces des principaux puits présents dans les microstructures sous irradiation (cavité, dislocation droite et boucle de dislocation) dans des configurations simples. L'étude met en évidence l'importance des interactions élastiques, et permet d'identifier l'anisotropie des dipôles élastiques au point col comme un paramètre très important qui modifie à la fois les trajectoires des défauts ponctuels et les valeurs de forces de puits. On s'intéresse dans une deuxième partie à l'importance de la position relative des puits sur leur capacité à absorber les défauts ponctuels. Des microstructures contenant un grand nombre de boucles de dislocation sont générées par des simulations OKMC et l'absorption des défauts ponctuels par ces microstructures est mesurée. Il est montré que le voisinage d'un puits modifie sensiblement sa capacité à absorber les défauts ponctuels et ce comportement est rationalisé en utilisant le volume de Voronoï associé à chaque puits. Cette étude conduit à proposer une nouvelle expression de force de puits, ainsi qu'un nouveau formalisme de DA dans lequel les amas sont caractérisés par leur taille et par le volume de leur cellule de Voronoï. Il est montré que ce formalisme permet d'améliorer fortement la prédiction de l'évolution des distributions de tailles d’amas lors d'une irradiation