Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Vacancy engineering”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Negative thermal expansion coefficient and amorphization in defective 4H-SiC

This paper presents thermal expansion coefficient (TEC) and amorphization in 4H-SiC containing point defects at different concentrations. We considered vacancy defects, interstitial defects, and Frenkel pair defects and investigated the thermomechanical response of the lattice over a wide range of temperatures using classical molecular dynamics simulations. The results show that 4H-SiC with vacancy defects exhibits a negative TEC above a critical defect density of around 9% (irrespective of the temperature). With interstitial defects, it exhibits a positive TEC (regardless of the defect density), and with Frenkel pair defects it shows a transition from positive TEC to negative TEC for a defect density greater than 8%. The coupling between temperature-induced expansion and defect-introduced stress in the lattice forms the mechanistic basis for the observed variation in TEC. Furthermore, the specific heat decreases rapidly with an increase in defect density at room temperature, with the highest sensitivity of the lattice observed for the Frenkel pair defects followed by interstitial defects and then by vacancy defects. Finally, these findings highlight the critical implications of defects on thermal expansion behavior of 4H-SiC with applications in radiation environments.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Defect-Driven Redox Interplay on Anatase TiO 2 : Surface-Structure Dependent Activation for CO 2 Hydrogenation Catalysis

Titanium dioxide (TiO 2 ) is one of the most extensively studied oxides as an active catalyst or catalyst support, particularly in energy and environmental applications, but the atomistic mechanisms governing its dynamic response to reactive environments and their correlation to reactivity remain largely elusive. Using in situ environmental transmission electron microscopy (ETEM), synchrotron X-ray diffraction (XRD), ambient-pressure X-ray photoelectron spectroscopy (AP-XPS), temperature-programmed reduction (TPR), reactivity measurements, and theoretical modeling, we reveal the dynamic interplay between oxygen loss and replenishment of anatase TiO 2 under varying reactive conditions. Under H 2 exposure, anatase TiO 2 undergoes surface reduction via lattice oxygen loss, forming Ti 3 O 5 . In contrast, CO 2 exposure induces oxygen replenishment, reversing stoichiometry. In mixed H2/CO 2 environments, the reverse water–gas shift (RWGS) reaction proceeds selectively on stepped and high-indexed TiO 2 surfaces, whereas the thermodynamically stable TiO 2 (101) surface remains inactive and intact. Critically, H 2 pretreatment generates oxygen vacancies on TiO 2 (101), transforming it into an active Ti 3 O 5 or defect-rich surface that catalyzes RWGS. By correlating surface structure, defect dynamics, and gas-phase interactions, this work deciphers the competition between H 2 -driven reduction and CO 2 -driven oxidation pathways at the atomic scale. Furthermore, these insights establish defect engineering as a strategic lever to activate inert TiO 2 facets, advancing the design of adaptive catalysts for sustainable fuel synthesis technologies.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Enhancing the Multifunctional Photocatalytic Activity of Sustainable Magnetic Nanoparticles

The development of efficient photocatalytic materials has intensified in response to increasing emphasis on sustainable energy conversion and environmental restoration. However, the excessive use and indiscriminate release of heavy metals from photocatalytic nanoparticles pose potential environmental risks. This study provides insights into optimizing visible-light-active photocatalysts to enhance their photocatalytic properties and stability. Specifically, cobalt doping and controlled pH modulation are employed on the modified Fe 3 O 4 , forming two distinct samples: Co@Fe 3 O 4 -B (base-treated) and Co@Fe 3 O 4 -A (acid-treated) nanoparticles. Microscopic characterization reveals that Co@Fe 3 O 4 -A undergoes a phase transition to hematite, whereas Co@Fe 3 O 4 -B retains its mixed-phase configuration. Spectroscopic analyses confirm that the cobalt dopants decorate the outskirts of each nanoparticle, forming a core–shell structure. However, Co@Fe 3 O 4 -B exhibit Co 2+ states with a high oxygen vacancy content, whereas Co@Fe 3 O 4 -A contain mixed Co 2+/3+ states. The high density of defect states in the Co@Fe 3 O 4 -B results in superior photocatalytic efficiency, achieving near-complete oxidation of furfuraldehyde (97% conversion) and 5-hydroxymethylfurfural (94% conversion), as well as effective degradation of toluene (78% conversion). This study demonstrates that the combined approach of doping and pH treatment is promising for the surface-defect engineering of photocatalysts, enhancing their multifunctional performance and reusability for sustainable energy conversion.

Fe3 O4 nanoparticles↗
DOE OSTI2025

Origin of enhanced performance when Mn-rich rocksalt cathodes transform to δ -DRX

Most Mn-rich cathodes are known to undergo phase transformation into structures resembling spinel-like ordering upon electrochemical cycling. Recently, the irreversible transformation of Ti-containing Mn-rich disordered rock-salt cathodes into a phase — named δ — with nanoscale spinel-like domains has been shown to increase energy density, capacity retention, and rate capability. However, the nature of the boundaries between domains and their relationship with composition and electrochemistry are not well understood. In this work, we discuss how the transformation into the multi-domain structure results in eight variants of Spinel domains, which is crucial for explaining the nanoscale domain formation in the δ -phase. We study the energetics of crystallographically unique boundaries and the possibility of Li-percolation across them with a fine-tuned CHGNet machine learning interatomic potential. Energetics of 16 d vacancies reveal a strong affinity to segregate to the boundaries, thereby opening Li-pathways at the boundary to enhance long-range Li-percolation in the δ structure. Defect calculations of the relatively low-mobility Ti show how it can influence the extent of Spinel ordering, domain morphology and size significantly; leading to guidelines for engineering electrochemical performance through changes in composition.

Anand, Shashwat↗
DOE OSTI2026

Entropy-defect synergy for dual luminescence mechanism in spinel: Time-resolved anti-counterfeiting and fingerprint visualization

Multimodal luminescent materials, while promising for anti-counterfeiting, often lack dynamic time-dependent responses and controllable spatial distribution, limiting their encryption capabilities in the spatiotemporal dimension. Here, this work presents a coordinated control strategy based on entropy and defect engineering, and uses a backpropagation (BP) neural network for material screening to successfully prepare spinel Mg 0.8 (Fe 0.04 Co 0.04 Ni 0.04 Cu 0.04 Zn 0.04 )Cr 2 O 4 (MgA 5 CO) phosphors with time-dependent dynamic luminescence behavior. This phosphor simultaneously activated the d-d transition luminescence (∼618 nm) derived from Co 2+ /Cr 3+ and the defect luminescence (∼398 nm) related to zinc vacancies (V Zn ) in a single-phase solid solution. The phosphor exhibits a time-dependent color evolution from pink to purple under fixed-wavelength excitation, due to the different excited-state dynamics and decay lifetimes associated with the d-d transition and defect luminescence. Structural characterization and spectral analysis confirmed the existence of V Zn and its significant role in defect luminescence process. The fluorescent and dynamic luminescent properties of entropy-based spinel oxide enable its use in advanced anti-counterfeiting applications like fingerprint recognition and color-changing dedicated anti-counterfeiting mark, showing promise in high-end and time-dynamic anti-counterfeiting fields. This research not only developed a new type of fluorescent dynamic anti-counterfeiting material, but also provided a new idea for constructing advanced optical functional materials with multiple luminescence mechanisms.

Defect project↗
DOE OSTI2024

Effect of Carbon on Void Nucleation in Iron

The study reports the significance of carbon presence in affecting void nucleation in Fe. Without carbon, void nucleation rates decrease gradually at high temperatures but remain significantly high and almost saturated at low temperatures. With carbon present, even at 1 atomic parts per million, void nucleation rates show a low-temperature cutoff. With higher carbon levels, the nucleation temperature window becomes narrower, the maximum nucleation rate becomes lower, and the temperature of maximum void nucleation shifts to a higher temperature. Fundamentally, this is caused by the change in effective vacancy diffusivity due to the formation of carbon-vacancy complexes. The high sensitivity of void nucleation to carbon comes from the high sensitivity of void nucleation to the vacancy arrival rate in a void. The void nucleation is calculated by first obtaining the effective vacancy diffusivity considering the carbon effect, then calculating the defect concentration and defect flux change considering both carbon effects and pre-existing dislocations, and finally calculating the void nucleation rate based on the recently corrected homogeneous void nucleation theory. The study is important not only in the fundamental understanding of impurity effects in ion/neutron irradiation but also in alloy engineering for judiciously introducing impurities to increase swelling resistance, as well as in the development of simulation and modeling methodologies applicable to other metals.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Ligation of Single-Site Ruthenium within Perovskite Oxides for Efficient Conversion of Thermodynamically Stable Molecules

Ru cation ligation in SrTiO 3 perovskite and their migration to the surface through exsolution are investigated for the dry reforming of methane (DRM), a chemistry that requires activation of two thermodynamically stable molecules, CH 4 and CO 2 . Compared to a supported 1 wt % Ru/SrTiO 3 benchmark, doped and exsolved Ru-SrTiO 3 demonstrate ≥3× higher CH 4 turnover rates (873 K), with isolated, ligated Ru exhibiting higher reactivity. Reactor studies assert that CH 4 and CO 2 activation are both kinetically relevant for CH 4 turnover rates on exsolved Ru-SrTiO 3 , unlike for supported Ru systems, where H-abstraction from CH 4 is the sole kinetically relevant step. As such, exsolved and doped Ru architectures are responsive toward co-reactant activation strategies, with their consequent reaction networks showing marked departures from those established for supported Ru catalysts. In situ spectroscopy and kinetic analyses propose distinct sensitivity toward CO 2 on Ru-SrTiO 3 systems, where higher CH4 turnover rates result from O-assisted C─H bond activation pathways. These pathways occur on paired Ru-oxygen vacancy sites that are inherent to the perovskite structure and are not readily accessible on supported Ru catalysts. Here, CO 2 is activated on oxygen vacancies directly adjacent to Ru active sites, which facilitates CH 4 C─H bond activation through a surface methoxy intermediate. The high reactivity of Ru-SrTiO 3 systems enables stable CH 4 turnover rates at milder reaction temperatures (673 K), conditions under which the supported counterpart, Ru/SrTiO 3 , is inactive. Overall, this work demonstrates that ligation of catalytically active cations in perovskite oxides facilitates site engineering toward atom-efficient activation of thermodynamically stable feedstocks.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Revealing Robust Room Temperature Ferromagnetism in Gd‐Doped Few‐Layered MoS 2 Thin Films

2D MoS 2 holds great promise for spintronics, yet is limited by intrinsic diamagnetism. This study demonstrates inducing ferromagnetic behavior in MoS 2 films doped with 0.47% Gd, achieving an ultrahigh saturation magnetization of 454 emu/cm 3 in a few‐layered film over 11‐times higher than bulk films (40 nm). Raman spectroscopy, X‐ray photoelectron spectroscopy, X‐ray magnetic circular dichroism, and density functional theory (DFT) calculations reveal an interplay between Gd dopants and Mo, S vacancies (V 1Mo+2S ), leading to the formation of bound magnetic polarons (BMPs) that drive ferromagnetic ordering. H 2 S annealing and DFT calculations reveal that defect healing reduces the saturation magnetization by 83%. High sulfur migration barrier in few‐layered films helps preserve BMPs, thereby sustaining ferromagnetism, whereas lower migration barriers in bulk films lead to suppression. These findings highlight the synergy between Gd doping and defect engineering in achieving ultrahigh room‐temperature ferromagnetism, offering a scalable strategy for developing high‐performance 2D magnetic materials for spintronic applications.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Synthesis of nanodiamonds encapsulated by zeolitic imidazole framework-8 for quantum sensing applications

Nitrogen vacancy (NV)-containing nanodiamonds are widely used in quantum sensing applications due to their high sensitivity to magnetic fields, relatively low cost, and ability to be initialized, manipulated, and read out at room temperature. Quantum sensing techniques such as optically detected magnetic resonance (ODMR) and spin relaxometry have exploited the sensitivity of the NV nanodiamonds to magnetic fields to detect a range of analytes, such as pH, metal ions, and biomolecules. However, diversifying the sensing targets accessible by NV diamond quantum sensors typically requires careful engineering of the diamond surface chemistry with stimuli-responsive functional groups. Here, a simple protocol for coating NV nanodiamonds with the zeolitic imidazole framework 8 (ZIF-8), a widely used metal-organic framework, is presented. ZIF-8 is a highly porous material that has been used as a selective sensor for gasses, metal ions, and other analytes. The material is well-characterized by x-ray diffraction, transmission electron microscopy, scanning electron microscopy, x-ray photoelectron spectroscopy, and luminescence spectroscopy. Encapsulation of NV nanodiamonds with a porous scaffold such as ZIF-8 provides a promising method for improving the selectivity for the quantum sensing of various analytes. Importantly, the ZIF-8 coating does not impact the luminescence properties of the NV diamond, which is a key readout in ODMR and spin relaxometry sensing approaches. Indeed, the ODMR spectra with and without the ZIF-8 shell is nearly identical. Moreover, the ZIF-8 coating increases the longitudinal spin relaxation time of the NV nanodiamond by a factor of 4 relative to aggregated diamond, a desirable outcome for spin relaxation-based quantum sensing. Metal-organic framework composites with nanodiamonds thus are an exciting strategy for enhancing NV nanodiamond performance in applications such as quantum sensing and quantum-enhanced nuclear magnetic resonance spectroscopy.

nitrogen vacancy nanodiamond↗
DOE OSTI2024

3D‐Mapping and Manipulation of Photocurrent in an Optoelectronic Diamond Device

Abstract Establishing connections between material impurities and charge transport properties in emerging electronic and quantum materials, such as wide‐bandgap semiconductors, demands new diagnostic methods tailored to these unique systems. Many such materials host optically‐active defect centers which offer a powerful in situ characterization system, but one that typically relies on the weak spin‐electric field coupling to measure electronic phenomena. In this work, charge‐state sensitive optical microscopy is combined with photoelectric detection of an array of nitrogen‐vacancy (NV) centers to directly image the flow of charge carriers inside a diamond optoelectronic device, in 3D and with temporal resolution. Optical control is used to change the charge state of background impurities inside the diamond on‐demand, resulting in drastically different current flow such as filamentary channels nucleating from specific, defective regions of the device. Conducting channels that control carrier flow, key steps toward optically reconfigurable, wide‐bandgap optoelectronics are then engineered using light. This work might be extended to probe other wide‐bandgap semiconductors (SiC, GaN) relevant to present and emerging electronic and quantum technologies.

Wood, Alexander A.↗
DOE OSTI2025

Atomic-Scale Insights into Carbon Dissolution in α-, γ-, and θ-Al 2 O 3 : Phase-dependent Transport Dynamics from First-Principles Calculations

α-Al 2 O 3 exhibits superior carburizing corrosion resistance compared to metastable γ-Al 2 O 3 and θ-Al 2 O 3 phases in high-temperature CO 2 environments, yet its atomic-scale origins remain unclear. Using first-principles density functional theory, we systematically investigate carbon dissolution and diffusion in α-Al 2 O 3 , γ-Al 2 O 3 , and θ-Al 2 O 3 , including the effects of oxygen (O) and aluminum (Al) vacancies. Our results show that α-Al 2 O 3 consistently exhibits higher carbon solution enthalpies than γ-Al 2 O 3 and θ-Al 2 O 3 in both pristine and defective structures, indicating lower intrinsic carbon solubility in α-Al 2 O 3 . Vacancies significantly enhance carbon incorporation: O vacancies reduce solution enthalpy, while Al vacancies further amplify this effect, with a strong preference for carbon at Al vacancy sites. Carbon diffusion barriers are also highest in α-Al 2 O 3 , reflecting slower carbon mobility. Al vacancies increase diffusion barriers across all phases, while O vacancies raise barriers in α- and γ-Al 2 O 3 but slightly lower them in θ-Al 2 O 3 . These results reveal a dual mechanism behind the carburizing resistance of α-Al 2 O 3 : reduced carbon solubility and elevated diffusion barriers. Furthermore, this work provides atomic-scale insights to guide the design of alumina-based materials with improved carburizing resistance through phase selection and defect engineering.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Defect Engineering of Bi 2 SeO 2 Thermoelectrics

Abstract Bi 2 SeO 2 is a promising n ‐type semiconductor to pair with p ‐type BiCuSeO in a thermoelectric (TE) device. The TE figure of merit zT and, therefore, the device efficiency must be optimized by tuning the carrier concentration. However, electron concentrations in self‐doped n ‐type Bi 2 SeO 2 span several orders of magnitude, even in samples with the same nominal compositions. Such unsystematic variations in the electron concentration have a thermodynamic origin related to the variations in native defect concentrations. In this study, first‐principles calculations are used to show that the selenium vacancy, which is the source of n ‐type conductivity in Bi 2 SeO 2 , varies by 1–2 orders of magnitude depending on the thermodynamic conditions. It is predicted that the electron concentration can be enhanced by synthesizing under more Se‐poor conditions and/or at higher solid‐state reaction temperatures ( T SSR ), which promote the formation of selenium vacancies without introducing extrinsic dopants. The computational predictions are validated through solid‐state synthesis of Bi 2 SeO 2 . More than two orders of magnitude increase are observed in the electron concentration simply by adjusting the synthesis conditions. Additionally, a significant effect of grain boundary scattering on the electron mobility in Bi 2 SeO 2 is revealed, which can also be controlled by adjusting T SSR . By simultaneously optimizing the electron concentration and mobility, a zT of ≈0.2 is achieved at 773 K for self‐doped n ‐type Bi 2 SeO 2 . The study highlights the need for careful control of thermodynamic growth conditions and demonstrates TE performance improvement by varying synthesis parameters according to thermodynamic guidelines.

Novitskii, Andrei [MANA, National Institute for Ma↗
DOE OSTI2024

Evaluation of δ-Phase ZrH1.4 to ZrH1.7 Thermal Neutron Scattering Laws Using Ab Initio Molecular Dynamics Simulations

Zirconium hydride is commonly used for next-generation reactor designs due to its excellent hydrogen retention capacity at temperatures below 1000 K. These types of reactors operate at thermal neutron energies and require accurate representation of thermal scattering laws (TSLs) to optimize moderator performance and evaluate the safety indicators for reactor design. In this work, we present an atomic-scale representation of sub-stoichiometric ZrH2−x(0.3≤x≤0.6), which relies on ab initio molecular dynamics (AIMD) in tandem with velocity auto-correlation (VAC) analysis to generate phonon density of states (DOS) for TSL development. The novel NJOY+NCrystal tool, developed by the European Spallation Source community, was utilized to generate the TSL formulations in the A Compact ENDF (ACE) format for its utility in neutron transport software. First, stoichiometric zirconium hydride cross sections were benchmarked with experiments. Then sub-stoichiometric zirconium hydride TSLs were developed. Significant deviations were observed between the new δ-phase ZrH2−x TSLs and the TSLs in the current ENDF release. It was also observed that varying the hydrogen vacancy defect concentration and sites did not cause as significant a change in the TSLs (e.g., ZrH1.4 vs. ZrH1.7) as was caused by the lattice transformation from ϵ- to δ-phase.

42 ENGINEERING↗
DOE OSTI2025

Pt Nanoparticle Disintegration at Oxide Interfaces Enhances CO Oxidation Catalysis

Understanding how supported metal nanoparticles dynamically evolve under reaction conditions is critical for controlling their catalytic function. Here, the mechanism behind the dynamic disintegration of Pt nanoparticles (NPs) supported on CeOx-TiO2 (CT) during CO oxidation is elucidated, leading to the formation of single atoms (SAs) and/or sub-nanometer clusters. Density functional theory (DFT) calculations reveal that strong Pt-CO interactions weaken Pt─Pt cohesion, while electronic coupling between Pt and Ce ions stabilizes Pt-CO* intermediates at the oxide interface. Surface oxygen vacancies kinetically trap Pt-CO*, but the vacancies are replenished under oxygen-rich conditions, enabling Pt-CO* surface diffusion and subsequent structural reorganization. In situ spectroscopic analyses confirm the oxygen-driven transformation of Pt NPs, correlating with a threefold increase in mass-specific activity at 150 °C. These findings demonstrate that interfacial oxygen dynamics and metal-support interactions can be leveraged to induce nanoparticle disintegration and optimize catalytic performance, highlighting the catalytic potential of interface-engineered Pt nanostructures.

CO oxidation↗
DOE OSTI2026

Unraveling Defect-Dependent Conductivity-Type Switching in CuFe 2 O 4 for Enhanced Photoelectrocatalytic Reduction of Benzaldehyde

Photoelectrocatalytic (PEC) reduction provides a sustainable route for upgrading biomass-derived feedstocks with reduced energy requirements, yet remains largely unexplored beyond hydrogen evolution and CO 2 reduction due to the scarcity of stable photocathodes. Here, we report a defect-engineered CuFe 2 O 4 photocathode that enables directly quantified PEC reduction of benzaldehyde to benzyl alcohol using a singlecomponent, earth-abundant oxide. By controlling annealing temperature and oxygen partial pressure, CuFe 2 O 4 is systematically tuned from n-type to p-type conductivity. Electrochemical measurements, X-ray and ultraviolet photoelectron spectroscopy, and firstprinciples defect calculations collectively show that oxygen-rich annealing conditions suppress deep donor-type oxygen vacancies while stabilizing shallow acceptor-type copper vacancies, resulting in enhanced hole concentration and improved charge transport. In a mixed acetonitrile/water electrolyte employing 1,4-benzoquinone as a redox mediator, the optimized CuFe 2 O 4 photocathode achieves stable photoelectrochemical operation over 18 h under continuous illumination with a benzyl alcohol production rate of 2.57 μmol/h at −0.50 V vs Ag/AgNO 3 , corresponding to a Faradaic efficiency of 51.3%. This PEC approach lowers the required applied potential by ∼1 V compared to traditional electrocatalytic methods, offering a more energy-efficient route for carbonyl reduction. These findings establish CuFe 2 O 4 as a viable photocathode platform for sustainable photoelectrocatalytic organic transformations under mild reaction conditions.

benzaldehyde reduction↗
DOE OSTI2024

Chemical Environment and Structural Variations in High Entropy Oxide Thin Film Probed with Electron Microscopy

For this work, we employ analytical transmission electron microscopy (TEM) to correlate the structural and chemical environment variations within a stacked epitaxial thin film of the high entropy oxide (HEO) Mg 0.2 Co 0.2 Ni 0.2 Cu 0.2 Zn 0.2 O (J14), with two layers grown at different substrate temperatures (500 and 200 °C) using pulsed laser deposition (PLD). Electron diffraction and atomically resolved STEM imaging reveal the difference in out-of-plane lattice parameters in the stacked thin film, which is further quantified on a larger scale using four-dimensional STEM (4D-STEM). In the layer deposited at a lower temperature, electron energy loss spectroscopy (EELS) mapping indicates drastic changes in the oxidation states and bonding environment for Co ions, and energy-dispersive X-ray spectroscopy (EDX) mapping detects more significant cation deficiency. Ab initio density functional theory (DFT) calculations validate that vacancies on the cation sublattice of J14 result in significant electronic and structural changes. The experimental and computational analyses indicate that low temperatures during film growth result in cation deficiency, an altered chemical environment, and reduced lattice parameters while maintaining a single phase. Our results demonstrate that the complex correlation of configurational entropy, kinetics, and thermodynamics can be utilized for accessing a range of metastable configurations in HEO materials without altering cation proportions, enabling further engineering of functional properties of HEO materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Meso-scale modeling of UO2 nuclear fuel to high burnup

To improve the economics of light water reactors for commercial nuclear energy generation, utility operators are seeking to obtain regulatory approval to run UO2 fuel to higher levels of burnup. One potential impediment to obtaining this approval is the phenomenon of fuel fragmentation, relocation, and dispersal (FFRD). FFRD can result when fuel experiences a rapid temperature transient, such as that occurring during a Loss Of Coolant Accident (LOCA). FFRD has historically been most associated with the rim region in UO2 fuel pellets, where the phenomenon of fragmentation is also referred to as pulverization due to the small size of the fragments. More recent evidence suggests that the so-called “dark zone” (due to its appearance in micrographs) that can be observed in the mid-radial regions of high burnup fuel is also susceptible to FFRD. Although empirical fuel performance models have been developed that can adequately predict pulverization in the rim region under typical LWR conditions, a scientific understanding of what underlies fuel restructuring and subsequent FFRD is lacking even in the rim region, and no models are currently available for the behavior the dark zone. To address these challenges, the U.S. Department of Energy’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program has employed a multi-scale modeling approach to improve scientific understanding and develop new fuel performance models. In this talk, I will focus on meso-scale efforts, which form a crucial link between atomic-scale and engineering-scale models. Phase-field modeling combined with cluster dynamics is used to predict the restructuring process in the rim region. Phase-field fracture modeling, informed by atomistic simulations, is used to predict the onset of pulverization in the rim region. Combining these techniques together allows the extent of rim pulverization to be predicted. The formation and evolution of the dark zone has also been simulated with the phase-field method, using an improved approach to vacancy source term parameterization. The work shows the important impact of microstructure on fuel performance.

fracture↗