Engineering topics

Beeler, Benjamin

Publications and source records attributed to Beeler, Benjamin.

6 records

DOE OSTI2024

Grain boundary self-diffusion and point defect interactions in α -U via molecular dynamics

Though metallic U-Zr fuel has been used in nuclear reactors since the 1960s, many of its fundamental and thermodynamic properties are still unknown. The a-U phase, which has a highly anisotropic crystal structure and physical properties, is present in U-Zr fuel. The character and behavior of a-U grain boundaries will strongly impact fuel thermophysical performance under irradiation. Here, we study the interaction of point defects with grain boundaries, diffusion along grain boundaries, and the predicted diffusional creep behavior of a-U via molecular dynamics. We calculate the segregation energy of vacancies and interstitials to grain boundaries and quantify the biased sink strength of the grain boundaries, and observe that this sink strength is not strongly dependent on the grain boundary orientation. We also find that grain boundary diffusivity is strongly dependent on the grain boundary energy and grain boundary orientation. The presence of point defects within the grain boundary can induce diffusion in grain boundaries with low formation energies and can enhance diffusion in high-energy grain boundaries. We also find that diffusional creep of a-U at prototypical metallic fuel operation conditions is extremely high and could help explain observed metallic fuel swelling behaviors.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2024

Vacancy-mediated transport and segregation tendencies of solutes in fcc nickel under diffusional creep: A density functional theory study

The Nabarro-Herring (NH) diffusional creep theory postulates the vacancy-mediated transport of atoms under a stress gradient as the creep mechanism under low-stress and high-temperature conditions. In multicomponent alloys, we premise that this stress-assisted flow of vacancies to and from grain boundaries will produce elemental segregation. An observation of such segregation, validated with theoretical predictions, can provide the necessary experimental evidence for the occurrence of NH creep. Theoretical calculations of the segregation tendencies via analyzing the dominant solute diffusion mechanisms and the difference in diffusivities of the elements are therefore essential. To this end, this study applies density functional theory calculations of migration barriers and solute-vacancy binding energies as input to the self-consistent mean-field theory to assess the vacancy-mediated diffusion mechanisms, transport coefficients, and segregation tendencies of Co, Cr, Mo, Re, Ta, and W solutes in face-centered-cubic Ni. We find Co, Re, and W to be slow diffusers at high temperatures and Cr, Mo, and Ta to be fast diffusers. Further analysis shows that the slow diffusers tend to always enrich at vacancy sinks over a wide range of temperatures. In contrast, the fast diffusers show a transition from depletion to enrichment as the temperature lowers. Furthermore, our analysis of the segregation tendencies under tensile hydrostatic strains shows that slow diffusers are largely unaffected by the strain and favor enrichment. On the other hand, the fast diffusers exhibit high sensitivity to strain and their segregation tendency can transition from depletion to enrichment at a given temperature. Finally, the transport coefficients calculated in this work are expected to serve as input to mesoscale microstructure models to provide a more rigorous assessment of solute segregation under NH creep conditions.

36 MATERIALS SCIENCE↗

DOE OSTI2024

First-principles investigation of lanthanides diffusion in HCP zirconium via vacancy-mediated transport

The diffusion of lanthanide fission products plays an important role in the growth of the fuel-cladding chemical interaction (FCCI) region in metallic fuels. The use of a Zr interdiffusion barrier may mitigate the transport of lanthanides from the fuel to the cladding, but the efficacy of such a liner is not yet known. In this paper, the stability and vacancy-mediated diffusion of La, Ce, Pr, and Nd in hexagonal close-packed (HCP) Zr is investigated via density functional theory (DFT) calculations and self-consistent mean field (SCMF) analysis. DFT is used to calculate the formation, binding, and migration energies of vacancies and vacancy-solute pairs. The DFT energetics are used in the KineCluE code to calculate the Onsager transport coefficients. La is found to be the fastest diffusing species in HCP Zr and experiences an almost isotropic diffusion behavior. The other three species (Ce, Pr, and Nd) demonstrate anisotropic diffusion where the diffusion in the basal planes is significantly faster than that along the c-axis. Further, the calculated lanthanide diffusivities in HCP Zr are fitted to an Arrhenius relation and the activation energies and prefactors are reported for the first time. Furthermore, the vacancy drag and the segregation tendencies were analyzed using the calculated off-diagonal transport coefficients. According to our vacancy-mediated diffusion model, lanthanides are expected to be enriched at vacancy sinks at low temperatures, while at high temperatures, lanthanides are depleted at sinks and will preferably diffuse into the bulk. The enrichment/depletion transition temperature depends on the diffusion direction (basal or axial) and hence will be controlled by the grain texture and orientation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2024

Assessment of uranium nitride interatomic potentials

Uranium mononitride (UN) is a promising nuclear fuel due to its high fissile density, high thermal conductivity, and suitability for reprocessing. In this study, two uranium nitride interatomic potentials are assessed: Tseplyaev and Starikov's angular-dependent potential and Kocevski et al.'s embedded atom model potential. Predictions of the thermophysical and elastic properties of UN, UN 2 , and α- and β-U 2 N 3 computed using both potentials are assessed and compared to available experimental data. Notably, the Tseplyaev potential performs better with the energetic aspects of UN, e.g., specific heat capacity and point defect formation energies, whereas the Kocevski potential performs better with the structural aspects of UN, e.g., thermal expansion as well as with the elastic properties. The reasons why the Kocevski potential underestimates the UN specific heat are explained by examining the UN phonon properties modeled using both potentials. The Kocevski potential shows better identification of the mechanical stability ranges of UN, UN 2 , and α- and β-U 2 N 3 , reasonably predicting the melting point of UN and predicting stable structures for UN 2 and α- and β-U 2 N 3 . On the other hand, the Tseplyaev potential predicts a premature phase change of both UN and UN 2 and cannot stabilize α- nor β-U 2 N 3 . However, the Kocevski potential cannot predict a stable α-U phase and is thus not suitable for the calculation of formation energies for non-stoichiometric point defects.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2024

KCl-UCl 3 molten salts investigated by Ab Initio Molecular Dynamics (AIMD) simulations

Ab Initio Molecular Dynamics (AIMD) simulations are performed on molten KCl-UCl3 salt mixtures to determine energies, heat capacities, and densities. The density-dependent energy correction (DFT-dDsC), Grimme et al.'s DFT-D3, and Langreth & Lundqvist (vdW-cx) models are used for dispersion forces and combined with the Perdew-Burke-Ernzerhof (PBE) exchange-correlation potential with a Hubbard U parameter for the 5f electrons of uranium. After validating predictions for the end-member systems to literature data, KCl-UCl 3 mixtures are studied at select temperatures. Densities and energies both deviate from ideal solution behavior, with the maximum deviation occurring around 36% UCl 3 for mixing energies and slightly lower (29% UCl3) for densities. Compared to the NaCl-UCl 3 system, which was previously investigated using the same simulation methodologies, the KCl-UCl 3 density and mixing energy deviations from ideal solution behavior are larger by almost a factor of two. No deviation from ideal solution behavior for heat capacity was observed. The AIMD predictions for mixing energies and densities agree qualitatively with experimental data, though the spread in data obtained from the various dispersion force models utilized, measurements, and empirical estimates makes strong conclusions difficult. The dependence of thermodynamic and thermophysical properties on composition is correlated with the local chemistry of the solution phase, in particular, the tendency of UCl 3 to form network structures.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Understanding the effect of crystal anisotropy on grain growth, texturing and transport via the orthorhombic ?-U system [Slides]

The alpha phase of uranium exhibits the orthorhombic crystal structure, resulting in significant physical property anisotropy. These complex properties make the material challenging to work with in engineering settings but provide a rich arena to investigate the fundamental, multi-scale effects of anisotropy on physical behaviors such as grain growth and microstructure evolution under irradiation. We apply molecular dynamics and phase field modeling to study mass transport and grain growth in alpha-uranium. We characterize the mobilities of crystalline defects and the interfacial energy, and We find that the grain boundary energy is highly variable depending on the interface structure and that the relative rates of defect transport change with temperature. We also find that anisotropic thermal expansion has a significant impact on grain growth kinetics and texture development. Our results are used to explain experimental observations and to provide a basis for further hypotheses into the complex physical behavior of alpha-uranium.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗