Anders David Ragnar Andersson

Effect of point defects on the thermal conductivity of UO2: molecular dynamics simulations

The thermal conductivity of uranium dioxide (UO2) fuel is an important materials property that affects fuel performance since it is a key parameter determining the temperature distribution in the fuel, thus governing, e.g., dimensional changes due to thermal expansion, fission gas release rates, etc. [1] The thermal conductivity of UO2 nuclear fuel is also affected by fission gas, fission products, defects, and microstructural features such as grain boundaries. Here, molecular dynamics (MD) simulations are carried out to determine quantitatively, the effect of irradiation induced point defects on the thermal conductivity of UO2, as a function of defect concentrations, for a range of temperatures, 300 – 1500 K. The results will be used to develop enhanced continuum thermal conductivity models for MARMOT and BISON by INL. These models express the thermal conductivity as a function of microstructure state-variables, thus enabling thermal conductivity models with closer connection to the physical state of the fuel [2].

Thermal transport in UO2 with defects and fission products by molecular dynamics simulations

The importance of the thermal transport in nuclear fuel has motivated a wide range of experimental and modelling studies. In this report, the reduction of thermal transport in UO2 due to defects and fission products has been investigated using non-equilibrium MD simulations, with two sets of empirical potentials for studying the degregation of UO2 thermal conductivity including a Buckingham type interatomic potential and a recently developed EAM type interatomic potential. Additional parameters for U5+ and Zr4+ in UO2 have been developed for the EAM potential. The thermal conductivity results from MD simulations are then corrected for the spin-phonon scattering through Callaway model formulations. To validate the modelling results, comparison was made with experimental measurements on single crystal hyper-stoichiometric UO2+x samples.

Kinetic Interaction of Uranium Vacancies and Dislocations in UO2

Understanding how point defects and defect clusters interact with dislocations in urania is important for a number of phenomena. For example, dislocations are one (of many) trap sites in the fuel where fission gases may accumulate and ultimately nucleate fission gas bubbles. Further, some creep mechanisms are governed by the flow of point defects to dislocations. Thus, for a variety of reasons, it is important to examine how dislocations attract and accelerate the kinetics of point defects.

The thermal conductivity of mixed fuel UxPu1-xO2: molecular dynamics simulations

Mixed oxides (MOX), in the context of nuclear fuels, are a mixture of the oxides of heavy actinide elements such as uranium, plutonium and thorium. The interest in the UO2-PuO2 system arises from the fact that these oxides are used both in fast breeder reactors (FBRs) as well as in pressurized water reactors (PWRs). The thermal conductivity of UO2 fuel is an important material property that affects fuel performance since it is the key parameter determining the temperature distribution in the fuel, thus governing, e.g., dimensional changes due to thermal expansion, fission gas release rates, etc. For this reason it is important to understand the thermal conductivity of MOX fuel and how it differs from UO2. Here, molecular dynamics (MD) simulations are carried out to determine quantitatively, the effect of mixing on the thermal conductivity of UxPu1-xO2, as a function of PuO2 concentrations, for a range of temperatures, 300 – 1500 K. The results will be used to develop enhanced continuum thermal conductivity models for MARMOT and BISON by INL. These models express the thermal conductivity as a function of microstructure state-variables, thus enabling thermal conductivity models with closer connection to the physical state of the fuel.

Derivation of effective fission gas diffusivities in UO2 from lower length scale simulations and implementation of fission gas diffusion models in BISON

This report summarizes the development of new fission gas diffusion models from lower length scale simulations and assessment of these models in terms of annealing experiments and fission gas release simulations using the BISON fuel performance code. Based on the mechanisms established from density functional theory (DFT) and empirical potential calculations, continuum models for diffusion of xenon (Xe) in UO2 were derived for both intrinsic conditions and under irradiation. The importance of the large XeU3O cluster (a Xe atom in a uranium + oxygen vacancy trap site with two bound uranium vacancies) is emphasized, which is a consequence of its high mobility and stability. These models were implemented in the MARMOT phase field code, which is used to calculate effective Xe diffusivities for various irradiation conditions. The effective diffusivities were used in BISON to calculate fission gas release for a number of test cases. The results are assessed against experimental data and future directions for research are outlined based on the conclusions.

Simulation of xenon, uranium vacancy and interstitial diffusion and grain boundary segregation in UO2

In light water reactor fuel, gaseous fission products segregate to grain boundaries, resulting in the nucleation and growth of large intergranular fission gas bubbles. Based on the mechanisms established from density functional theory (DFT) and empirical potential calculations 1, continuum models for diffusion of xenon (Xe), uranium (U) vacancies and U interstitials in UO2 have been derived for both intrinsic conditions and under irradiation. Segregation of Xe to grain boundaries is described by combining the bulk diffusion model with a model for the interaction between Xe atoms and three different grain boundaries in UO2 ( Σ5 tilt, Σ5 twist and a high angle random boundary),as derived from atomistic calculations. All models are implemented in the MARMOT phase field code, which is used to calculate effective Xe and U diffusivities as well as redistribution for a few simple microstructures.

Simulations of Xe and U diffusion in UO2

Diffusion of xenon (Xe) and uranium (U) in UO{sub 2} is controlled by vacancy mechanisms and under irradiation the formation of mobile vacancy clusters is important. Based on the vacancy and cluster diffusion mechanisms established from density functional theory (DFT) calculations, we derive continuum thermodynamic and diffusion models for Xe and U in UO{sub 2}. In order to capture the effects of irradiation, vacancies (Va) are explicitly coupled to the Xe and U dynamics. Segregation of defects to grain boundaries in UO{sub 2} is described by combining the bulk diffusion model with models of the interaction between Xe atoms and vacancies with grain boundaries, which were derived from atomistic calculations. The diffusion and segregation models were implemented in the MOOSE-Bison-Marmot (MBM) finite element (FEM) framework and the Xe/U redistribution was simulated for a few simple microstructures.

MARMOT simulations of Xe segregation to grain boundaries in UO2

Diffusion of Xe and U in UO{sub 2} is controlled by vacancy mechanisms and under irradiation the formation of mobile vacancy clusters is important. We derive continuum thermodynamic and diffusion models for Xe and U in UO{sub 2} based on the vacancy and cluster diffusion mechanisms established from recent density functional theory (DFT) calculations. Segregation of defects to grain boundaries in UO{sub 2} is described by combining the diffusion model with models of the interaction between Xe atoms and vacancies with grain boundaries derived from separate atomistic calculations. The diffusion and segregation models are implemented in the MOOSE/MARMOT (MBM) finite element (FEM) framework and we simulate Xe redistribution for a few simple microstructures. In this report we focus on segregation to grain boundaries. The U or vacancy diffusion model as well as the coupled diffusion of vacancies and Xe have also been implemented, but results are not included in this report.