Engineering Papers⌕ Search

Engineering topics

Beeler, Benjamin

Publications and source records attributed to Beeler, Benjamin.

At least 37 records · Page 2

Density functional theory calculations of the thermodynamic and kinetic properties of point defects in β -U

Density functional theory (DFT) calculations of the thermodynamic and kinetic properties of point defects in the β phase of uranium are reported. Defect energies and entropies were calculated using 2 x 2 x 2 supercells and the Generalized Gradient Approximation (GGA) of Perdew-Burke-Ernzerhof (PBE) for the exchange-correlation potential. Due to computational cost, calculations of the vibrational properties governing entropies were performed by only displacing atoms within (roughly) the 3rd or 4th nearest neighbor shell of the defect, which implicitly assumes that atoms beyond this distance are unaffected by the defect. Migration barriers were estimated by nudged elastic band (NEB) calculations. The low symmetry of the β-U phase (the unit cell is tetragonal and contains 30 atoms) results in many point defect configurations and even more migration pathways. A connectivity map, starting from the most stable point defects, was developed in order to identify the rate-limiting step controlling the net diffusion rate in each crystallographic direction. The uranium self-diffusivity tensor was calculated by combining the defect formation energies, entropies, migration barriers and attempt frequencies. The fastest diffusion rate was determined to be a vacancy mechanism in the z crystallographic direction. The predicted uranium self-diffusivity for this mechanism agrees well with available experimental data. Finally, the diffusion mechanisms and rates identified in this study will inform fuel performance models of swelling and gas evolution.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Irradiation-enhanced diffusion and diffusion-limited creep in U 3 Si 2

U 3 Si 2 is an advanced fuel candidate due to its relatively high fissile density and attractive thermal properties. Compared to standard UO 2 fuel, there are significant data gaps for the thermophysical and thermomechanical properties of U 3 Si 2 . Point defect concentrations and mobilities under irradiation govern a number of important fuel performance properties, such as creep and fission gas release. In this work, we utilized density functional theory (DFT) data to inform a cluster dynamics framework to predict point defect concentrations in U 3 Si 2 under irradiation. Molecular dynamics (MD) simulations were used to examine the contribution of atomic mixing during ballistic cascades to diffusion, as well as the diffusivity of U and Si at grain boundaries. These atomic scale models for diffusivity were then used to inform a creep model based on bulk (Nabarro-Herring) and grain boundary (Coble) diffusional creep, and climb-limited dislocation creep. The model compares well against available experimental data and has been implemented in the BISON fuel performance code. Finally, a demonstration case using simple power profiles has been carried out, showing that negligible creep occurs due to the low temperatures experienced by U 3 Si 2 in-reactor, a consequence of its high thermal conductivity.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Determination of Thermal Expansion, Defect Formation Energy, and Defect-Induced Strain of α-U Via ab Initio Molecular Dynamics

Uranium (U) is often alloyed with molybdenum (Mo) or zirconium (Zr) in order to stabilize its high-temperature body-centered cubic phase for use in nuclear reactors. However, in all metallic fuel forms, the α phase of U remains in some fraction. This phase decomposition due to temperature or compositional variance can play an outsized role on fuel performance and microstructural evolution. Relatively little is known about fundamental point defect properties in α -U at non-zero temperatures, from either computational or experimental studies. This work performs the first thorough evaluation of the α phase of U via ab initio molecular dynamics (AIMD). A number of thermophysical properties are calculated as a function of temperature, including equilibrium lattice parameters, thermal expansion, and heat capacity. These results indicate a two-region behavior, with the transition at 400 K. The thermal expansion/contraction in the a / b direction occurs rapidly from 100 up to 400 K, after which a more linear and gradual change in the lattice constant takes place. The volumetric expansion matches experiments quantitatively, but the individual lattice constant expansion only matches experiments qualitatively. Point defect formation energies and induced lattice strains are also determined as a function of temperature, providing insight on defect populations and the fundamentals of irradiation growth in α -U. Interstitials induce significantly more strain than vacancies, and the nature of that strain is highly dependent on the individual lattice directions. The direction of point defect-induced lattice strain is contrary to the irradiation growth behavior of α -U. This work shows that isolated point defects cannot be the primary driving force responsible for the significant irradiation-induced growth of α -U observed experimentally.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Gas Bubble Evolution in Polycrystalline UMo Fuels Under Elastic-Plastic Deformation: A Phase-Field Model With Crystal-Plasticity

In monolithic UMo fuels, the interaction between the Al cladding and large gas bubble volumetric swelling causes both elastic-plastic and creep deformation. In this work, a phase-field model of gas bubble evolution in polycrystalline UMo under elastic-plastic deformation was developed for studying the dynamic interaction between evolving gas bubble/voids and deformation. A crystal plasticity model, which assumes that the plastic strain rate is proportional to resolved shear stresses of dislocation slip systems on their slip planes, was used to describe plastic deformation in polycrystalline UMo. Xe diffusion and gas bubble evolution are driven by the minimization of chemical and deformation energies in the phase-field model, while evolving gas bubble structure was used to update the mechanical properties in the crystal plasticity model. With the developed model, we simulated the effect of gas bubble structures (different volume fractions and internal gas pressures) on stress-strain curves and the effect of local stresses on gas bubble evolution. The results show that 1) the effective Young’s modulus and yield stress decrease with the increase of gas bubble volume fraction; 2) the hardening coefficient increases with the increase of gas bubble volume fraction, especially for gas bubbles with higher internal pressure; and 3) the pressure dependence of Xe thermodynamic and kinetic properties in addition to the local stress state determine gas bubble growth or shrinkage. The simulated results can serve as a guide to improve material property models for macroscale fuel performance modeling.

36 MATERIALS SCIENCE↗

Evaluation of the anisotropic grain boundaries and surfaces of α -U via molecular dynamics

Alloys based on uranium-zirconium are gaining renewed interest as fuels for the Versatile Test Reactor and a number of microreactor designs. Implementing metallic fuel in reactors creates the need for robust descriptive and predictive fuel performance modeling. Here, the current state of metallic fuel performance modeling relies on empirical equations derived from historical experiments, which may be unreliable when applied outside of their temperature, power, and composition phase space. One area where such data is lacking is the irradiation behavior of alpha -U, specifically tearing and porosity formation at the early stages of irradiation. While grain boundaries likely play a key role in this fuel behavior, relatively little is known about grain boundaries in alpha-U. Thus, we evaluate the grain boundary, surface energy, and work of adhesion of alpha -U utilizing molecular dynamics. Symmetric tilt grain boundaries (STGBs) are analyzed with the tilt plane oriented along each major crystallographic axis, for a total of eighty unique grain boundaries. The effect of temperature, tilt plane, and misorientation angle on interfacial energies are analyzed. The interfacial energies typically increase with temperature and there is significant variance as a function of misorientation angle, irrespective of the tilt plane. At 500 K, the average surface energy (1.23 J/m 2 ) is approximately 1.5 times the grain boundary energy (0.79 J/m 2 ), and the work of adhesion is approximately twice the grain boundary energy (1.68 J/m 2 ). Orientations for the likely formation of twins and likely failure planes are identified.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Defect cluster and nonequilibrium gas bubble associated growth in irradiated UMo fuels – a cluster dynamics and phase field model

Irradiation examination shows that gas bubble swelling kinetics is much faster after irradiation-induced recrystallization than that prior recrystallization in UMo fuels. It implies that gas bubbles in coarse grains and small recrystallized grains have different growth behavior. In this work, a phase-field model of gas bubble evolution integrating microstructure dependent cluster dynamics has been developed, for the first time, to study the gas bubble swelling behavior in the recrystallization zone of UMo fuels. Generation, diffusion, reaction, sink, emission and clustering of vacancies and interstitials are described by the cluster dynamics model while a phase-field model is used to describe the evolution of non-equilibrium gas bubbles including nucleation and growth. With the coupled model, the effect of defect generation rate, clustering rate, interstitial emission and sink rates on grain boundaries on the gas bubble evolution are systematically simulated. A set of model parameters (defect generation rate, clustering rate, interstitial emission and sink rates) is determined by comparing measured and simulated gas bubble swelling kinetics. The results demonstrate that interstitial clustering is one of the important physical mechanisms which results in a fast gas bubble swelling kinetics in the recrystallization zone. The developed model can also be extended to study the associated growth of defect and second phase precipitates often observed in irradiated materials.

Hu, Shenyang↗

First-principles study of the surface properties of uranium carbides

Uranium carbides have attracted renewed interest as advanced nuclear fuels for Generation IV reactors. As an important property required for gas bubble modeling in nuclear fuels, the surface energy of uranium carbides is scarce in literature. In this work, we study the surface properties of uranium carbides by first-principles density functional theory calculations. Surface orientations with maximum Miller index up to 3, 2 and 2 are investigated for UC, U 2 C 3 and α-UC 2 , respectively. By studying the effects of surface termination and chemical potential on surface energy, we identify the factors that determines the surface stability. From the calculated surface energies, the surface properties of uranium carbide single crystals are obtained from Wulff construction, including equilibrium morphology, dominant surface orientation and area weighted surface energy.

36 MATERIALS SCIENCE↗

Mesoscale modeling to inform Bison models of accident tolerant fuel concepts

U 3 Si 2 and doped UO 2 have been under investigation in recent years as potential accident-tolerant fuel concepts. In this report, lower-length scale studies of these fuel concepts carried out during Fiscal Year 2020 are detailed. U 3 Si 2 is a potential accident-tolerant fuel that shows promise due to its high thermal conductivity and higher uranium density relative to UO 2 . However, its swelling and fission gas release behavior in light water reactor (LWR) conditions is relatively unknown. To provide mechanistic insight and determine parameters for engineering-scale fuel performance modeling of pellet-form U 3 Si 2 , phase-field simulations of the growth, interconnection, and venting of intergranular fission gas bubbles were performed. The fractional coverage of the grain boundary and the fraction of bubble area that is vented were calculated as a function of time. From the simulation data, the fractional grain boundary coverage at saturation, an important parameter needed in engineering-scale modeling of swelling and fission gas release, was determined. Multiple simulations were run to determine the uncertainty in the calculated value. The effect of model assumptions and input parameters that are not well known was evaluated. Simulation results are compared to related theoretical and computational work. Based on the simulation results, a value of 0.60 for the fractional grain boundary coverage at saturation is recommended for U 3 Si 2 fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Microstructure, Thermal, and Mechanical Properties Relationships in U and UZr Alloys (Final Report)

Uranium-zirconium (U-Zr) alloys are candidate fuel systems for transmutation based reactors that can be used to burn long-lived minor actinides and fission products in fast spectrum reactors. Metallic fuels have also been gaining more recent attention for applications as accident tolerant fuels, as well as for use in small modular reactors. This research focused on a “science-based” approach to capture the connections between U and U-Zr alloys’ three-dimensional (3-D) microstructure, thermal properties, and mechanical properties through closely coordinated experiments and modeling efforts from the unirradiated to the irradiated fuels. Advanced characterization and modeling techniques were used to understand irradiation-induced microstructural evolution and its direct impact on the thermal and mechanical properties of U and U-Zr fuel. Closely coordinated experiments and modeling were performed to provide crucial data that does not currently exist. Overall, this research spanned multiple length and time scales within the models and experiments. The scope of the research encompassed the understanding of the irradiation effects in U and various U-Zr alloys with particular attention paid to three task areas: (1) microstructural evolution, (2) in-situ/ex-situ thermal and mechanical properties, and (3) multiscale modeling connections to microstructure, thermal, and mechanical properties. This research resulted in (1) the 3-D characterization of neutron irradiated U-Zr fuel in multiple phase regions to better understand fission gas swelling and constituent redistribution, (2) development of a microstructural model linking thermal and mechanical properties via in situ Raman and nanoindentation, (3) and mesoscale phase field modeling was coupled with the AEH method in the MOOSE framework was used to calculate the effective thermal conductivities of U-Zr fuels consisting of α-U and δ-UZr 2 heterogeneous microstructures.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Research Needs for Uranium-Zirconium- Based Metallic Fuels

The purpose of this report is to summarize the needs and present recommendations related to the future direction for U-Zr-based metallic fuel research (including binary U-Zr and ternary U-Pu-Zr alloys). These needs and recommendations were determined by subject matter experts from various institutions during a two-day workshop held at the University of Florida in November, 2019. During open-floor discussions, the highest priority gaps in our understanding of U-Zr-based fuels were down-selected, and near- and long-term needs that directly impact the implementation of these metallic fuels were identified. The identified near-term needs include investigation of the following phenomena: i) swelling and fission gas release, ii) fuel-cladding chemical interaction, iii) phase evolution/constituent redistribution, and iv) thermal properties of the fuels. The long-term needs are: i) investigation of fuel creep and plasticity and ii) fission product (lanthanide) transport. In addition, there was general agreement that all institutions and subject matter experts would benefit from an open-source metallic fuels database with thermophysical property and microstructural data, along with fuel operation/irradiation history, which should be regularly updated with vetted information from new experimental and computational investigations and used to advance metallic fuels research and development. Finally, we recommend that metallic fuel research should be ongoing and that fuel qualification and fuel optimization should be equally prioritized; research combining experiments with modeling and simulation has the largest potential impact.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A improved equation of state for Xe gas bubbles in gamma U-Mo fuels

A monolithic fuel design based on a U-Mo alloy has been selected as the fuel type for conversion of the United States High-Performance Research Reactors (HPRRs). An issue with U-Mo monolithic fuel is the large amount of swelling that takes place during operation. The accurate prediction of fuel evolution under irradiation requires implementation of correct thermodynamic properties into mesoscale and continuum level fuel performance modeling codes. However, the thermodynamic properties of the fission gas bubbles (such as the relationship among bubble size, equilibrium Xe concentration, and bubble pressure) are not well known. This work studies Xe bubbles in ?U-Mo from a diameter of 3 nm up to 8.5 nm and from 400 K up to 700 K. The energetic relationship of Xe bubbles with regard to voids and Xe substitutional atoms is described. The transition is also determined for when a bubble becomes over-pressurized. Finally, an equation of state is fit to the pressure as a function of molar volume and temperature.

Beeler, Benjamin↗

Development of a grain growth model for U 3 Si 2 using experimental data, phase field simulation and molecular dynamics

The purpose of this work is to develop a model for normal grain growth in U 3 Si 2 . The average grain boundary energy was determined from previously published molecular dynamics simulations. The grain growth kinetics were quantified at various temperatures by annealing nanocrystalline samples. Here, the mobility was determined by comparing phase field grain growth simulations to the experimental data.

36 MATERIALS SCIENCE↗

Status of LEU U-10Mo Monolithic Fuel Testing and Qualification Efforts to Support US High Performance Research Reactor Conversions

The US High Performance Research Reactor Project has begun a new phase of fuel testing and post irradiation examination (PIE) of plate-type low-enriched uranium (LEU) U-10Mo Monolithic fuel. The testing campaigns are intended to first confirm that the fuel performance behavior of commercially produced fuel matches that of the lab-scale fabricated fuel tested during development, followed by qualification of the U-10Mo monolithic fuel system through the US Nuclear Regulatory Commission. Planned tests include miniature-plate, full-size plate, and fuel element irradiations, during which fuel test specimens fabricated with qualified processes will be irradiated under conditions that bound the NRC-regulated reactors in order to demonstrate that the fuel meets the overarching fuel performance requirements of 1) maintaining mechanical integrity 2) maintaining geometric stability and 3) exhibiting stable and predictable behavior. The MP-1 test, with the objective of confirming fuel performance behavior, is currently undergoing irradiation in the Advanced Test Reactor and PIE in INL hot cells. Status of upcoming qualification and licensing test design activities and highlights of initial as-fabricated characterization and testing results from the MP-1 test will be provided.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A Monte Carlo Model of Irradiation-Induced Recrystallization in Polycrystalline UMo Fuels

Experiments show that irradiation-induced recrystallization speeds up the swelling kinetics in U10wt%Mo fuels. However, recrystallization mechanisms and the effect of initial grain microstructures on recrystallization kinetics are still unclear. In this work a Monte Carlo model coupling the rate theory of defect evolution has been developed to study the irradiation-induced recrystallization. The rate theory is used to describe the spatial evolution of gas bubbles, interstitials and interstitial loops; First-Passage Kinetic Monte Carlo (FPKMC) approach is used to describe the fast and strongly anisotropic migration of interstitials, and a Cellular Automata method is used to model the formation of recrystallized grains. With the assumption that recrystallization may occur when the local interstitial loop density is larger than a given critical value, simulation results reveal that 1) recrystallized grains first nucleate on grain boundaries and the recrystallization zone front moves to the center of original coarse grains in the UMo matrix, and 2) recrystallization starts earlier in coarse polycrystalline structures, while the overall recrystallization kinetics decreases with increasing grain size. These results agree with experimental observations. The comparison of recrystallization kinetics obtained from experiments and modeling suggests that the interstitial loop accumulation leads to the recrystallization and the interstitial loop growth is suppressed inside coarse grains due to the over-pressured intra-granular gas bubbles.

Frazier, William E.↗