Engineering topics

Beeler, Benjamin W.

Publications and source records attributed to Beeler, Benjamin W..

10 records

DOE OSTI2024

A Finite Difference informed Random Walk solver for simulating radiation defect evolution in polycrystalline structures with strongly inhomogeneous diffusivity

Diffusivity of species and defects on grain boundaries is usually several orders of magnitude larger than that inside grains. Such strongly inhomogeneous diffusivity requires prohibitively high computational demands for modeling microstructural evolution. Here, this paper presents a highly-efficient numerical solver, combining the Finite Difference method and Random Walk model, designed for accurately modeling strongly inhomogeneous diffusion within polycrystalline structures. The proposed solver, termed Finite Difference informed Random Walk (FDiRW), integrates a customized Finite Difference (cFD) scheme tailored for fast diffusion along thin grain boundaries represented by a single-layer of nodes. Numerical experiments demonstrate that the FDiRW solver achieves an impressive efficiency gain of 1560x compared to traditional Finite Difference methods while maintaining accuracy, making it feasible for personal computer machines to handle diffusional systems with strongly inhomogeneous diffusivity across static polycrystalline microstructures. The model has been successfully applied to simulate radiation defect evolution, showcasing its scalability to engineering scales in both length and time dimensions.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Assessment of effective elastic constants of U-10Mo fuel: A multiscale modeling and homogenization study

The significant microstructural changes that U-Mo fuel undergoes during operation degrades its mechanical properties and structural integrity. Microstructural evolution entails the formation, evolution, and redistribution of porosity in conjunction with grain refinement. In the present paper, we employ numerical approaches to assess the impact of the various microstructural features—grains, nanoscale intragranular fission gas bubbles, and mesoscale intergranular voids—on the degradation of elastic constants. Phase-field microstructure models are combined with the asymptotic expansion homogenization technique in order to derive the effective elastic constants as a function of porosity and fission density. Here the results are verified and compared against theoretical bounds. Using this approach, elastic degradation in operating nuclear fuels can be quantified when the distributions of microstructural features are known.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2023

Microstructural and phase changes in alpha uranium investigated via in-situ studies and molecular dynamics

A deeper knowledge of thermally induced microstructural and phase evolution in nuclear metallic fuel can be obtained using novel in-situ microscopic analyses. Such studies can provide information on the dynamics of phase transitions which is not possible with conventional postmortem characterization (post-irradiation examination). In this work, the behavior of alpha uranium (α-U) was investigated via in-situ heating tests in a transmission electron microscope. Here, the main objective is to understand the microstructural and phase changes, such as defect annihilation and ß phase formation and retention, observed in reactor in-pile transient studies at the Transient Reactor Test facility. Indeed, defect migration and rearrangement were observed within the α phase starting at 673 K; α→ß phase transition was observed between 773 K and 1,073 K during the heating ramp (which is in the temperature window reported for α→ß transition temperatures). Recrystallization and formation of nano grains was observed at high temperatures (over 1,073 K). Such recrystallization was possibly related to the formation of the γ phase. Finally, it was indeed observed that the ß phase (but not γ phase) was retained at room temperature upon rapid cooling. Molecular dynamics studies support these experimental results and shows that the γ phase of pure uranium cannot be retained at room temperature if not stabilized with the addition of an alloying element.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2022

An ab initio molecular dynamics investigation of the thermophysical properties of molten NaCl-MgCl 2

Molten salts have many applications in the nuclear and solar energy industries for thermal storage and heat transfer applications. However, there is a knowledge gap in molten salt thermophysical properties which hinders the technical readiness level of molten salt applications, especially in the nuclear industry. A common method of investigating new materials is through ab initio Molecular Dynamics (AIMD) simulations which is an effective tool to investigate structural and thermophysical properties at realistic temperatures. NaCl-MgCl 2 is an inexpensive salt that is a good candidate for use as a heat transfer medium in solar power applications or in the secondary loop of a nuclear reactor. In this article, the thermophysical properties of NaCl-MgCl 2 are computed via AIMD calculations to supplement the limited experimental studies in the literature. Here a wide range of compositions and temperatures for the pseudo-binary NaCl-MgCl 2 were used to calculate the density, heat capacity, compressibility, enthalpy of mixing, and volumetric thermal expansion coefficient. AIMD is shown to accurately model the densities of molten NaCl-MgCl 2 as there is good agreement with the available literature. This work observed a transition to a monotonic increase of the density with respect to MgCl 2 composition occurring above 1100 K. The heat capacity values increase uniformly with respect to concentration of MgCl 2 at a rate of 2.85 J/mol-K per 10 mol% of MgCl 2 . Select thermophysical properties are fit to a Redlich-Kister expansion for utilization in multiphysics simulations.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Phase-field simulations of fission gas bubble growth and interconnection in U-(Pu)-Zr nuclear fuel

Abstract The growth and interconnection of fission gas bubbles in the hotter central regions of U-(Pu)-Zr nuclear fuel has been simulated with a phase-field model. The Cahn-Hilliard equation was used to represent the two-phase microstructure, with a single defect species. The volume fraction of the bubble phase and surface area of the bubble-matrix interface were determined during growth and interconnection. Surface area increased rapidly during the initial stages of growth, then slowed and finally decreased as bubble interconnection began and coarsening acted to reduce surface area. The fraction of the bubbles vented to a simulation domain boundary, f V , was quantified as a measure of the microstructure’s interconnectivity and plotted as a function of porosity p . The defect species diffusivity was varied; although changes in diffusivity significantly affected the microstructure, the plots of f V vs. p did not change significantly. The percolation threshold p c was calculated to be approximately 0.26, depending on the assumed diffusivity and using an initial bubble number density based on experimental observations. This is slightly smaller than the percolation threshold for continuum percolation of overlapping 3D spheres. The simulation results were used to parameterize two different engineering-scale swelling models for U-(Pu)-Zr in the nuclear fuel performance code BISON.

Aagesen, Larry K. (ORCID:000000034936676X)↗

DOE OSTI2021

Radiation driven diffusion in γU-Mo

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). A critical phenomenon of interest with U-Mo monolithic fuel is the large amount of swelling that takes place during operation, particularly at high fission densities. The accurate prediction of fuel evolution under irradiation requires implementation of correct thermodynamic and kinetic properties into mesoscale and continuum level fuel performance modeling codes. One such property where there exists incomplete data is the diffusion of relevant species under irradiation. Fuel performance swelling predictions rely on an accurate representation of diffusion in order to determine the rate of fission gas swelling and the local microstructural evolution. In this work, we present molecular dynamics simulations of the radiation driven diffusion of U, Mo and Xe in U-Mo nuclear fuels. Diffusion coefficients for each species are determined over a range of temperatures and compositions. In this work, updated diffusion coefficients are presented that are applicable under irradiation that incorporate both intrinsic and radiation driven diffusion.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2021

Microstructural-Level Fuel Performance Modeling of U-Mo Monolithic Fuel

As the physics that governs the microstructural evolution of nuclear fuel span various time and spatial scales, to fully understand the fuel behavior inevitably involves atomic to mesoscale resolution that can be difficult to determine experimentally. Microstructural-level modeling and simulations can be used to develop physics-based materials models that can provide physical understanding to inform fabrication process control, as well as a valuable feedback mechanism between post-irradiation examination (PIE) results and fabrication parameters. In accordance with the program schedule, the primary goals of the microstructure modeling effort are to: 1. Address critical microstructural questions and provide practical guidance to the fabricator via the fuel product specification; 2. Provide mechanistic inputs for the existing fuel performance code to improve its descriptive and predictive capability at the macroscopic scale. In fiscal year (FY)-21, the work scope consisted of six main facets: (1) the effect of carbides on fuel performance; (2) gas diffusivity in different phases; (3) integration of microstructural fuel performance modeling; (4) property degradation; (5) irradiation creep; and (6) historical analysis of microstructure data. Brief summaries of each are included below.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2020

Phase-field simulations of intergranular fission gas bubble behavior in U 3 Si 2 nuclear fuel

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. Finally, 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↗

DOE OSTI2020

AFIP6-MkII and RERTR-12 Porosity Data Collection and Analysis for Modeling and Simulation

Gathering data for the improvement of nuclear fuel modeling and simulation efforts is the primary driver for this work. Mechanistic models allow for a better understanding of the material on a micro- and macrostructural level while saving time and money over traditional experiment efforts. Historically, summarized data and correlations are the inputs for empirical material models and model validation. When improving these models for nuclear fuels with experimental results, there is a lack of reliable data readily available. Experiments - RERTR-12 and AFIP6-MkII - were conducted to understand the irradiation behavior of metallic U-10Mo monolithic fuels for use in extreme reactor environments such as research reactors like the Advanced Test Reactor (ATR) or the High Flux Isotope Reactor (HFIR). Microstructural characteristics of fission gas pores (FGP) in each experiment are collected using an automated image analysis technique developed at the University of Florida and presented here. A series of statistical tests are performed to explore the reliability of the results, as well as understand where the data is lacking and what future data collection is necessary to provide sufficient information to assist modeling efforts. The focus is on the porosity, pore size, and eccentricity of FGPs formed during irradiation in three AFIP6-MkII samples and one RERTR-12 sample. From the analysis, it is clear there are substantial impacts of fission density on the pore structure, but there also exist also underlying connections between each sample and the behavior observed in the pores. Further analyses of the pre- and post-irradiation microstructure are needed to improve the understanding of these connections. An early method for microstructural data analysis is presented within and is currently being expanded to include other microstructure data.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2019

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.↗