1 |
Simulation of thermomechanical properties of U-PuO2 nuclear fuel under irradiation / Simulations des propriétés thermomécaniques du combustible nucléaire (U,PuO2) sous irradiationBalboa lópez, Hector 10 December 2018 (has links)
L’objectif de cette thèse a été d’utiliser une approche numérique pour étudier l’impact des dommages d'irradiation sur la microstructure du combustible nucléaire composé d'un mélange d’oxyde de plutonium et d’uranium (MOX). Cette approche, réalisée à l'échelle atomique, repose sur l'utilisation de potentiels empiriques développés dans la littérature pour l'oxyde mixte $(U,Pu)O_2$ dans l’approximation des interactions d’ions rigides.Une première étape a été une analyse critique des propriétés structurales, thermodynamiques et mécaniques prédites par 5 potentiels de la littérature. Les calculs de dynamique moléculaire ont été réalisés à l'aide du code LAMMPS, sur l'ensemble de la gamme de composition de $UO_2$ à $PuO_2$ et à des températures comprises entre 300 K et le point de fusion. Nous montrons que les potentiels les plus satisfaisants sont ceux développés par Cooper et Potashnikov. Ces deux potentiels reproduisent correctement la stabilité thermodynamique des phases ainsi que l'évolution en température des paramètres de maille et de la chaleur spécifique. Cependant le comportement mécanique des oxydes est différent selon le potentiel choisi. Tout d'abord, on note que les constantes élastiques et le facteur d'anisotropie obtenus par le potentiel Cooper sont plus fidèles aux recommandations de la littérature. Ensuite, nous montons que la propagation d'une fracture induit une transformation de phase en pointe de fissure avec le potentiel de Cooper alors qu'une fissuration fragile est observée avec le potentiel Potashnikov.Une seconde étape a été l'étude des dommages d'irradiation dans les oxydes mixtes en se limitant à l'utilisation des deux potentiels de Cooper et de Potashnikov. Des calculs de dynamique moléculaire de cascades de déplacement ont été réalisés à différentes énergies, température et compositions en plutonium. Ensuite, afin d'évaluer les dommages d'irradiation sur des temps plus longs que ceux accessibles en dynamique moléculaire, la méthode d’accumulation de défauts a été utilisée. Pour les deux potentiels, l’évolution des dommages primaires avec l’augmentation de la dose suit globalement les mêmes étapes que celles trouvées précédemment dans $UO_2$. Tout d’abord, les défauts ponctuels sont créés. Par la suite, ils se regroupent et forment de petites boucles de Frank, qui après se transforment en dislocations parfaites. Cependant, la cinétique de la recombinaison de défauts ponctuels est significativement plus lente avec le potentiel de Cooper ce qui conduit à la création de petites régions désordonnées pour les cascades d'énergie élevée. L’effet de la teneur en plutonium est analysé. Nous montrons en particulier que la densité de dislocations crée diminue lorsque la teneur en plutonium augmente.Bien que la dynamique moléculaire a été décrite comme un microscope moléculaire en raison de sa capacité à décrire avec précision des systèmes atomiques, elle présente un inconvénient majeur, celui lié aux temps de l’ordre de la femto-seconde nécessaires pour résoudre les vibrations atomiques. Cela limite le temps total de simulation approximativement quelques nanosecondes. Afin de stimuler les processus, tels que la diffusion de cations, un autre outil de calcul est nécessaire. Les techniques de Monte Carlo (KMC) atomiques peuvent simuler de plus longtemps que la dynamique moléculaire . Cependant, pour que KMC fonctionne avec précision, il est nécessaire de connaître a priori tous les mécanismes de transition entre les états possibles. Pour cette raison, la method de Monte Carlo cinétique adaptative (aKMC) est choisie pour surmonter ces limitations. Cette méthode détermine les états de transition disponibles pendant la simulation. De cette manière, elle entraîne le système dans des états imprévus via des mécanismes complexes. La puissance de cet outil s'est révélée efficace pour découvrir la recombinaison de cations sur de plus longues périodes de temps que la DM. / The objective of this doctoral research is to use a numerical approach to study the impact of irradiation damage on the microstructure of the mixed uranium-plutonium oxide fuel (MOX). This numerical approach comprises mainly the use of Molecular Dynamics (MD) using empirical potential. Several empirical potentials for $(U,Pu)O_2$ can be found in the literature. The results of these potentials can exhibit significant differences. For this reason an extensive assessment of the main empirical potential found in the literature had to be performed.Five empirical interatomic potentials were assessed in the approximation of rigid ions and pair interactions for the $(U_{1-y},Pu_y)O_2$ solid solution. Simulations were carried out on the structural, thermodynamics and mechanical properties over the full range of plutonium composition, meaning from pure $UO_2$ to pure $PuO_2$ and for temperatures ranging from 300 K up to the melting point. The best results are obtained by potentials referred as Cooper and Potashnikov. The first one reproduces more accurately recommendations for the thermodynamics and mechanical properties exhibiting ductile-like behaviour during crack propagation, whereas the second one gives brittle behaviour at low temperature.From our results from the empirical potentials assessment, we can move to the radiation damage using only two potentials (Cooper and Potashnikov). In order to know the main source of defect during irradiation, MD displacement cascades were simulated. This revealed the damage created due to varying projectile energies. In addition, the Frenkel pair accumulation method was chosen to investigate the dose effect. This method circumvents the highly computing time demanding accumulation of displacement cascade by directly creating their final states, i.e. mainly point defects. Overall, results obtained with both potentials show the same trend. However, kinetics of point defect recombination are significantly slower with Cooper potential implying creation of small disordered region with high energy displacement cascades. The evolution of the primary damage with increasing dose follows the same steps as those found previously in pure $UO_2$. First, point defects are created. Subsequently, they cluster and form small Frank loops, which in turn transform and grow into unfaulted loops. We demonstrate also that increasing temperatures accelerate the production of dislocations shifting their creation to lower doses. The effect of the plutonium content is also evidenced, especially with Cooper potential. It shows that the dislocation density decreases when the plutonium content increases.Although, MD has been described as a molecular microscope due to its ability to discribe accuratily systems of atoms, it has a large drawback that is the short time steps of the order of femto-seconds needed to resolve the atomic vibrations. This limits the time typically few microsecond. In order to invetigate processess, such as, cation diffusion and rare-event annihilation of defects after cascaces, another computational tool is required. Atomistic or object kinetic Monte Carlo (KMC) techniques can run for longer timescales than MD. However, for KMC to work accurately, all of the possible inter-state transition mechanisms and their associated rates need to be known a priori. For this reason, the adaptive kinetic Monte Carlo (AKMC) is chosen to overcome these limitations. This method determines the available transition states during simulation. In this way, it takes the system into unforeseen states via complex mechanisms. The power and range of this tool proved to be efficient to discover cation Frenkel pair recombination over a longer periods of time than MD.
|
2 |
Gradient-damage modeling of dynamic brittle fracture : variational principles and numerical simulations / Analyse de la rupture dynamique fragile via les modèles d'endommagement à gradient : principes variationnels et simulations numériquesLi, Tianyi 06 October 2016 (has links)
Une bonne tenue mécanique des structures du génie civil en béton armé sous chargements dynamiques sévères est primordiale pour la sécurité et nécessite une évaluation précise de leur comportement en présence de propagation dynamique de fissures. Dans ce travail, on se focalise sur la modélisation constitutive du béton assimilé à un matériau élastique-fragile endommageable. La localisation des déformations sera régie par un modèle d'endommagement à gradient où un champ scalaire réalise une description régularisée des phénomènes de rupture dynamique. La contribution de cette étude est à la fois théorique et numérique. On propose une formulation variationnelle des modèles d'endommagement à gradient en dynamique. Une définition rigoureuse de plusieurs taux de restitution d'énergie dans le modèle d'endommagement est donnée et on démontre que la propagation dynamique de fissures est régie par un critère de Griffith généralisé. On décrit ensuite une implémentation numérique efficace basée sur une discrétisation par éléments finis standards en espace et la méthode de Newmark en temps dans un cadre de calcul parallèle. Les résultats de simulation de plusieurs problèmes modèles sont discutés d'un point de vue numérique et physique. Les lois constitutives d'endommagement et les formulations d'asymétrie en traction et compression sont comparées par rapport à leur aptitude à modéliser la rupture fragile. Les propriétés spécifiques du modèle d'endommagement à gradient en dynamique sont analysées pour différentes phases de l'évolution de fissures : nucléation, initiation, propagation, arrêt, branchement et bifurcation. Des comparaisons avec les résultats expérimentaux sont aussi réalisées afin de valider le modèle et proposer des axes d'amélioration. / In civil engineering, mechanical integrity of the reinforced concrete structures under severe transient dynamic loading conditions is of paramount importance for safety and calls for an accurate assessment of structural behaviors in presence of dynamic crack propagation. In this work, we focus on the constitutive modeling of concrete regarded as an elastic-damage brittle material. The strain localization evolution is governed by a gradient-damage approach where a scalar field achieves a smeared description of dynamic fracture phenomena. The contribution of the present work is both theoretical and numerical. We propose a variationally consistent formulation of dynamic gradient damage models. A formal definition of several energy release rate concepts in the gradient damage model is given and we show that the dynamic crack tip equation of motion is governed by a generalized Griffith criterion. We then give an efficient numerical implementation of the model based on a standard finite-element spatial discretization and the Newmark time-stepping methods in a parallel computing framework. Simulation results of several problems are discussed both from a computational and physical point of view. Different damage constitutive laws and tension-compression asymmetry formulations are compared with respect to their aptitude to approximate brittle fracture. Specific properties of the dynamic gradient damage model are investigated for different phases of the crack evolution: nucleation, initiation, propagation, arrest, kinking and branching. Comparisons with experimental results are also performed in order to validate the model and indicate its further improvement.
|
Page generated in 0.0606 seconds