Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “boundary effects”

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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Investigation of the effect of structural defects from hydride precipitation on superconducting properties of high purity SRF cavity Nb using magneto-optical and electron imaging methods

Abstract Nb is an elemental superconductor with a critical temperature of 9.3 K and is widely used to fabricate superconducting radiofrequency (SRF) cavities for particle accelerators. However, microstructural defects in Nb, such as grain boundaries (GBs) and dislocations, can act as pinning centers for magnetic flux that can degrade SRF cavity performance. Hydrogen contamination is also detrimental to SRF cavity performance due to the formation of normal conducting hydrides during cool down. In this study, disc shaped Nb bi-crystals extracted from high-purity large-grain Nb slices were investigated to study the effects of GBs, hydrogen, and dislocations on superconducting properties. Grain orientation and GB misorientation were measured using Laue x-ray diffraction and electron backscattered diffraction (EBSD) analyses. Cryogenic magneto-optical imaging was used to directly observe magnetic flux penetration below T c = 9.3 K. Damage caused by low temperature precipitation of hydrides and their dissolution upon reheating after cryogenic cycles was examined using electron channeling contrast imaging, and EBSD. The relationships between hydride formation, dislocation content, GBs, cryo-cooling, heat treatment (HT), and flux penetration indicate that both GB character and hydrogen content affect magnetic flux penetration. Such flux penetration could be facilitated by dislocation structures and low angle GBs resulting from hydride precipitation and HT.

36 MATERIALS SCIENCE↗
DOE OSTI2026

A mechanistic, multiscale model for predicting Pd penetration in TRISO fuels using BISON

TRistructural ISOtropic (TRISO) particles use silicon carbide (SiC) as the primary structural member and barrier against metallic fission product (FP) release. palladiums (PDs), produced by fission in the fuel kernel, can diffuse to and chemically interact with the SiC layer, degrading its structural integrity and ability to contain radioactive FPs. Existing temperature-dependent correlations for predicting Pd penetration rely on experimental data with significant scatter due to varying conditions, potentially complicating ongoing fuel qualification and licensing efforts for advanced reactors that would subject TRISO fuels to operating conditions outside of those examined in the experiments. A mechanistic model of Pd production, transport, and reaction is developed in this work to better understand and predict PD attack of SiC in TRISO particles. molecular dynamicss (MDs) simulations are utilized to calculate the Pd diffusivity in SiC grain bulk and grain boundaries. A mesoscale phase-field diffusion model, informed by the MD diffusivities, is used to develop a reduced order model (ROM) for the effect of SiC microstructure on the Pd penetration rate. The engineering scale BISON model calculates Pd production and transport, utilizing the ROM to predict penetration rates consistent with experimental data. This novel mechanistic ROM captures the effect of temperature, microstructure, and irradiation history on Pd penetration. In conclusion, these new capabilities are expected to support ongoing qualification and licensing efforts associated with near-term TRISO-fueled reactor applications.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Local inversion-symmetry breaking in a bismuthate high-Tc superconductor

Abstract The doped perovskite BaBiO 3 exhibits a maximum superconducting transition temperature ( T c ) of 34 K and was the first high- T c oxide to be discovered, yet pivotal questions regarding the nature of both the metallic and superconducting states remain unresolved. Although it is generally thought that superconductivity in the bismuthates is of the conventional s -wave type, the pairing mechanism is still debated, with strong electron-phonon coupling and bismuth valence or bond disproportionation possibly playing a role. Here we use diffuse x-ray scattering and Monte Carlo modeling to study the local structure of Ba 1-x K x BiO 3 across its insulator-metal boundary. We find no evidence for either long- or short-range disproportionation, which resolves a major conundrum, as disproportionation and the related polaronic effects are likely not relevant for the metallic and superconducting states. Instead, we uncover nanoscale structural correlations that break inversion symmetry, with far-reaching implications for the electronic physics. This unexpected finding furthermore establishes that the bismuthates belong to the broader classes of materials with hidden spin-orbit coupling and a tendency towards inversion-breaking displacements.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Photometry, Centroid and Point-spread Function Measurements in the LSST Camera Focal Plane Using Artificial Stars

Abstract The Vera C. Rubin Observatory’s LSST Camera (LSSTCam) pixel response has been characterized using laboratory measurements with a grid of artificial stars. We quantify the contributions to photometry, centroid, point-spread function size, and shape measurement errors due to small anomalies in the LSSTCam CCDs. The main sources of those anomalies are quantum efficiency variations and pixel area variations induced by the amplifier segmentation boundaries and “tree-rings”—circular variations in silicon doping concentration. This laboratory study using artificial stars projected on the sensors shows overall small effects. The residual effects on point-spread function (PSF) size and shape are below 0.1%, meeting the ten-year LSST survey science requirements. However, the CCD mid-line presents distortions that can have a moderate impact on PSF measurements. This feature can be avoided by masking the affected regions. Effects of tree-rings are observed on centroids and PSFs of the artificial stars and the nature of the effect is confirmed by a study of the flat-field response. Nevertheless, further studies of the full-focal plane with stellar data should more completely probe variations and might reveal new features, e.g., wavelength-dependent effects. The results of this study can be used as a guide for the on-sky operation of LSSTCam.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Mapping the unoccupied state dispersions in Ta 2 NiSe 5 with resonant inelastic x-ray scattering

The transition metal chalcogenide Ta 2 NiSe 5 undergoes a second-order phase transition at T c = 328 K involving a small lattice distortion. Below T c , a band gap at the center of its Brillouin zone increases up to about 0.35 eV. In this work, we study the electronic structure of Ta 2 NiSe 5 in its low-temperature semiconducting phase, using resonant inelastic x-ray scattering (RIXS) at the Ni L 3 -edge. In addition to a weak fluorescence response, we observe a collection of intense Raman-like peaks that we attribute to electron-hole excitations. Using density functional theory calculations of its electronic band structure, we identify the main Raman-like peaks as interband transitions between valence and conduction bands. Further, by performing angle-dependent RIXS measurements, we uncover the dispersion of these electron-hole excitations that allows us to extract the low-energy boundary of the electron-hole continuum. From the dispersion of the valence band measured by angle-resolved photoemission spectroscopy, we derive the effective mass of the lowest unoccupied conduction band.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Electronic and magnetic phase diagrams of the Kitaev quantum spin liquid candidate Na 2 Co 2 TeO 6

The 3⁢d 7 Co 2+ -based insulating magnet Na 2 ⁢Co 2 ⁢TeO 6 has recently been reported to have strong Kitaev interactions on a honeycomb lattice and is thus being considered as a Kitaev quandum spin liquid candidate. However, due to the existence of other types of interactions, a spontaneous long-range magnetic order occurs. This order is suppressed by applied magnetic fields leading to a succession of phases and ultimately saturation of the magnetic moments. The precise phase diagram, the nature of the phases, and the possibility that one of the field-induced phases is a Kitaev quantum spin liquid phase are still a matter of debate. Here, in this study, we measured an extensive set of physical properties to build the complete temperature-field phase diagrams to magnetic saturation at 10 T for magnetic fields along the a and a* axes, and a partial phase diagram up to 60 T along c. We probe the phases using magnetization, specific heat, magnetocaloric effect, magnetostriction, dielectric constant, and electric polarization, which is a symmetry-sensitive probe. With these measurements, we identify all the previously incomplete phase boundaries and find additional high-field phase boundaries. We find strong magnetoelectric coupling in the dielectric constant and moderate magnetostrictive coupling at several phase boundaries. Furthermore, we detect the symmetry of the magnetic order using electrical polarization measurements under magnetic fields. Based on our analysis, the absence of electric polarization under zero or finite magnetic field in any of the phases or after any combination of magnetic/electric field cooling suggests that a zigzag spin structure is more likely than a triple-Q spin structure at zero field. Finally, we investigate the hysteresis and first- or second-order nature of each phase transition and its entropy changes. With this information, we establish a map of the magnetic phases of this compound and its magnetic, thermodynamic, and magnetoelectric properties, and discuss where spin liquid or other phases may be sought in future studies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

A New Shape Optimization Approach for Lightweighting Electric Machines Inspired by Additive Manufacturing

Minimizing the mass in electric machines while maintaining superior performance is a new requirement for the advancement of drivetrains used in wind energy and electric mobility. Topology optimization (TO) for lightweighting electric machines using traditional approaches typically explores a restricted design space allowed by standard parametrizable geometry and manufacturing, while advanced methods, such as cell-based density approaches, suffer from a lack of robust manufacturability constraints during the optimization process. To overcome these drawbacks, we explore a grid-independent, boundary optimization where the outer shape of the magnet is parameterized using Bezier curves. We conduct a design of experiments (DOE) to study the effect of different magnet shapes on machine performance by varying the control points on the Bezier curves. A machine-learning-based surrogate model is constructed using the data from the DOE to quantify the relationship between the control points, air-gap torque, and mass. The control points are then optimized to maximize the torque density. The approach is used for minimizing electrical steel mass in the International Energy Agency (IEA) 15-MW radial flux direct-drive wind turbine generator. The new approach to shape optimization resulted in smooth and concise shapes that can be easily additively manufactured with up to a 20-ton reduction in electrical steel mass.

additive manufacturing↗
DOE OSTI2023

Cross-Scale Land-Atmosphere Experiment (CSLAEX) (Final Report)

It is now well recognized that land-atmosphere interactions play an important role in both weather and climate. These interactions take place over a wide range of spatial and temporal scales. Most of the current formulations of the momentum, heat and moisture transport laws are informed by observations with a limited temporal or spatial scale, and hence may have limited applicability. The objective of this work is to improve the multi-scale representation of the near-surface heat transfer in the ground and in the atmosphere using a combination of modeling and field observations. Our knowledge and predictive capacity in land-atmosphere coupling is limited by our capacity to observe the full-scale spectrum at play. In order to bridge this fundamental gap in our observation capacity I will work to: create the first long-term Integrated Cross-Scale land-Atmosphere Experiment (CSLAEX) since I will work to: measure the components of the surface energy balance at a spatial scale relevant to remote sensing and meteorological applications (5km), along with observations of the nested subscale processes with Distributed Temperature Sensors (DTS). This proposal aims to fill the gaps in modeling and observing land-atmosphere interactions across scales via a novel Cross-Scale Land-Atmosphere Experiment (CSLAEX) to be implemented as part of the Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) site near Lamont, Oklahoma. The proposed research effort aims at exploring land-atmosphere processes across multiple spatial scales and to investigate the effect of landscape heterogeneity using Distributed Temperature Sensor measurements. The novel instrumentation of the Surface Energy Balance and simultaneous observations of the atmospheric state (boundary layer, cloud, convection) available at the DOE ARM SGP site will be used to create CSLAEX. This dataset will be used to frame and evaluate new transport laws and parameterizations for the boundary layer and the land surface, which can be implemented into the Weather and Research Forecasting (WRF) model and Community Atmospheric Model (CAM). The result of this study should impact multiple fields: meteorology, climate, and hydrology.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Cybersecurity Lessons Learned from Vehicle to Grid Engagement

As the transportation industry continues to become electrified, introduction of additional digital devices within associated actions such as recharging bring additional potential for cybersecurity attacks. Devices that are designed, implemented, and operated with cybersecurity as a crucial consideration exacerbate these concerns by failing to provide strict boundaries on access to and use of the equipment. Emerging use cases such as Vehicle to Grid (V2G) charging may expand the potential physical effects of a cybersecurity attack by providing indirect access to electrical components of a building microgrid or portions of the larger power grid. This paper serves as an overview of findings and recommendations based on cybersecurity testing performed at a V2G implementation site operated by a member of the Memorandum of Understanding (MOU) to Establish the Vehicle-to-Everything (V2X) Collaboration [1]. The Department of Energy Office of Cybersecurity, Energy Security, and Emergency Response is a signatory of the MOU, and has funded this research paper and associated body of work regarding V2X cybersecurity. Sandia has a large background of previous research focused on Electric Vehicle (EV) cybersecurity, such as reference [2], which includes an overall survey of EV infrastructure cybersecurity and recommendations based on those findings. This report seeks to expand knowledge of EV cybersecurity status and needs by focusing on a specific implementation of V2G charging, and providing recommendations based on the relevant findings. This report serves as a publicly available, sanitized description of applied vulnerability testing on an operational V2G implementation. A more in-depth technical version of the report is provided to the MOU partner, but not available at the time of writing due to inclusion of proprietary information. V2G charging comes with many research problems that must be solved before the technology can securely implemented in sites with unrestricted public access or where cybersecurity attacks could have increased consequences, such as government offices. V2G charging requires many stakeholders such as end users, host sites, equipment vendors, and integrators, which all rely on operational safety and security as well as security and trustworthiness of any associated financial transactions.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2021

Exploring the origins of the indentation size effect at submicron scales

Significance Downscaling material structure into ever-decreasing levels for modern technological advancements demands small-scale mechanical characterizations of materials, such as nanoindentation. Thus, a comprehensive understanding of indentation size effect (ISE) has a far-reaching impact, but such understanding remains incomplete at the submicron scale despite a decades-long quest since the 1980s. Here, we explored the origins of this effect experimentally by linking dislocation behavior evolution with nanoindentation mechanical data at various depths. Importantly, the dislocation behavior is found to be strongly depth dependent and gives rise to subgrain formation progressively. Aided by subgrain boundaries, the ISE mechanism transitions from dislocation source starvation to dislocation interaction as indentation depth increases. The critical mechanism transition elucidates the distinctive ISE behavior at the submicron scale.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Strong and almost strong modes of Floquet spin chains in Krylov subspaces

Integrable Floquet spin chains are known to host strong zero and π modes which are boundary operators that respectively commute and anticommute with the Floquet unitary generating stroboscopic time evolution, in addition to anticommuting with a discrete symmetry of the Floquet unitary. Thus the existence of strong modes implies a characteristic pairing structure of the full spectrum. Weak interactions modify the strong modes to almost strong modes that almost commute or anticommute with the Floquet unitary. Manifestations of strong and almost strong modes are presented in two different Krylov subspaces. One is a Krylov subspace obtained from a Lanczos iteration that maps the time evolution generated by the Floquet Hamiltonian onto dynamics of a single particle on a fictitious chain with nearest-neighbor hopping. The second is a Krylov subspace obtained from the Arnoldi iteration that maps the time evolution generated directly by the Floquet unitary onto dynamics of a single particle on a fictitious chain with longer-range hopping. While the former Krylov subspace is sensitive to the branch of the logarithm of the Floquet unitary, the latter obtained from the Arnoldi scheme is not. In this work, the effective single-particle models in the Krylov subspace are discussed, and the topological properties of the Krylov chain that ensure stable zero and π modes at the boundaries are highlighted. The role of interactions is discussed. Expressions for the lifetime of the almost strong modes are derived in terms of the parameters of the Krylov subspace, and are compared with exact diagonalization.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

De Sitter space is sometimes not empty

Multiple lines of evidence suggest that the Hilbert space of an isolated de Sitter universe is one dimensional but can appear larger when probed by a gravitating observer. To test this idea, we compute the von Neumann entropy of a field theory in a two-dimensional de Sitter universe which is entangled in a thermal-like state with the same field theory on a disjoint, asymptotically anti-de Sitter (AdS) black hole. Previously, it was shown that the replica trick for computing the entropy of such entangled gravitating systems requires the inclusion of a non-perturbative effect in quantum gravity — novel wormholes connecting the two spaces. Here we show that: (a) the expected wormholes connecting de Sitter and AdS universes exist, avoiding a no-go theorem via the presence of sources on the AdS boundary; (b) the entanglement entropy vanishes if the nominal entropy of the de Sitter cosmological horizon (S ds = $A^{dS}_{horizon}$/4G N ) is less than the entropy of the AdS black hole horizon (S BH = $A^{AdS}_{horizon}$/4G N ), i.e., S dS < S BH ; (c) the entanglement entropy is finite when S dS > S BH . Thus, the de Sitter Hilbert space is effectively nontrivial only when S dS > S BH . The AdS black hole we introduce can be regarded as an “observer” for de Sitter space. In this sense, our result is a non-perturbative generalization of the recent perturbative argument that the algebra of observables on the de Sitter static patch becomes nontrivial in the presence of an observer.

2D gravity↗
DOE OSTI2022

MECHANISTIC MULTIPHYSICS MODELING OF CLADDING RUPTURE IN NUCLEAR FUEL RODS DURING LOSS-OF-COOLANT ACCIDENT CONDITIONS

The Loss of Coolant Accident (LOCA) is a design basis accident that is included as part of the safety analysis of nuclear power plants. As the nuclear industry desires to increase the cycle length and discharge burnup of existing nuclear power plants they must demonstrate safe operation during a LOCA on high burnup fuel. During a LOCA transient on high burnup fuel rods, the rods may undergo a process known as fuel fragmentation, relocation, and dispersal (FFRD). To permit dispersal, the cladding encapsulating the fuel must rupture with an opening size large enough to allow the fragmented fuel particles to release. Current licensing tools used by industry and the United States Nuclear Regulatory Commission are limited in geometric fidelity and materials that can be analyzed. These simulation tools generally employ a quasi-two-dimensional (1.5D or Layered1D) or 2D-RZ axisymmetric geometric representations exclusively. While a valid approach under some instances, there are times when important physics have an asymmetric behavior in the fuel rod. Examples include fuel fragmentation, thermal-hydraulic boundary conditions, and cladding rupture, all of which are important for LOCA analysis. As industry pursues burnup extensions it must be demonstrated that fuel dispersal can be mitigated or eliminated. To do this, an understanding of the rupture opening after cladding failure is required. This work presents the development of a model for predicting the size and location of the rupture opening in failed fuel rods during LOCA conditions using advanced modeling and simulation tools. In order to supply the rupture model with appropriate boundary conditions, improvements to fuel fragmentation, axial relocation and oxidation modeling were required. First, the eXtended Finite Element Method (XFEM) is used to mechanistically predict the number of fuel fragments that form due to material strength randomization, criteria for strength randomization, mesh density, power ramping rates and irradiation effects. These predictions with associated uncertainty were compared to empirical correlations developed for UO2 verifying that they can be used with increased confidence in subsequent axial relocation analyses. Secondly, a new first-of-its-kind Layered2D computational framework was developed that provides the ability to apply azimuthally varying boundary conditions while providing discrete layers to track fuel movement during the LOCA. An existing fuel axial relocation model developed for Layered1D was extended to work within the Layered2D framework. A large sensitivity study was performed on the initial version of the model to identify modeling parameters of particular importance, with the emissivity used for radiation after blowdown being the primary source of uncertainty. Then, a simplistic approach to incorporate mechanical degradation of the cladding due to oxidation was also developed to investigate the impact of reduced cladding thickness on predictions of the time to failure of cladding tubes. It was found that the cladding will typically rupture prior to a reduction in thickness that is sufficient to impact the rupture behavior. A model was then developed for predicting cladding rupture that transfers the cladding surface temperatures, rod internal and external pressures, fast neutron fluence, and fast neutron flux from a more detailed Layered1D, Layered2D, or 2D-RZ analysis to a 3D cladding only analysis. Comparisons of the rupture model to a few experiments indicate reasonable predictions. The rupture model was then applied to two accident tolerant fuel concepts (FeCrAl and Cr-coated Zircaloy) where it predicted that both ATF concepts would have smaller rupture openings and delayed rupture times than the standard Zircaloy-4 cladding material under identical loading conditions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

An ion specific continuum model on ionic surfactant's binary phase diagram in aqueous solution

Here, a continuum aggregation model is proposed to account for specific ion effects, enabling accurate phase diagrams calculations for the micellar, cylindrical and lamellar aggregates of sodium/potassium carboxylate surfactants in aqueous solution across a range of temperatures. Three groups of concentrations at distinctive temperatures are fitted to build empirical temperature dependence given the limited availability of relevant experimental measurements. The specific ion effects are manifested in the aggregates' surface tension as well as in the distributions of counter-ions' concentrations in the vicinity of the aggregates. The aggregates' geometric sizes are well-reproduced. The differential evolution algorithm is applied to address boundary conditions of the electrostatic potential, aggregate size optimization as well as the equilibrium of monomers transferring between the aggregate and aqueous region.

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

Exascale granular microstructure reconstruction in 3D volumes of arbitrary geometries with generative learning

Reconstructing 3D granular microstructures within volumes of arbitrary geometries from limited 2D image data is crucial for predicting the material properties, as well as performances of structural components accounting for material microstructural effects. We present a novel generative learning framework that enables exascale reconstruction of granular microstructures within complex 3D geometric volumes. Building upon existing transfer learning techniques using pre-trained convolutional neural networks (CNN), we introduce several key innovations to overcome the difficulties inherent in arbitrary geometries. Our framework incorporates periodic boundary conditions using circular padding techniques, ensuring continuity and representativeness of the reconstructed microstructures. We also introduce a novel seamless transition reconstruction (STR) method that creates statistically equivalent transition zones to integrate multiple pre-existing 3D microstructure volumes. Based on STR, we propose a cost-effective strategy for reconstructing microstructures within complex geometric volumes, minimizing computational waste. Validation through numerical experiments using kinetic Monte Carlo simulations demonstrates accurate reproduction of grain statistics, including grain size distributions and morphology. A case study involving the reconstruction of a 4-blade propeller microstructure illustrates the method’s capability to efficiently handle complex geometries. In conclusion, the proposed framework significantly reduces computational demands while maintaining high reconstruction quality, paving the way for scalable microstructure reconstruction in materials design and analysis.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Modeling the effect of wind speed and direction shear on utility‐scale wind turbine power production

Abstract Wind speed and direction variations across the rotor affect power production. As utility‐scale turbines extend higher into the atmospheric boundary layer (ABL) with larger rotor diameters and hub heights, they increasingly encounter more complex wind speed and direction variations. We assess three models for power production that account for wind speed and direction shear. Two are based on actuator disc representations, and the third is a blade element representation. We also evaluate the predictions from a standard power curve model that has no knowledge of wind shear. The predictions from each model, driven by wind profile measurements from a profiling LiDAR, are compared to concurrent power measurements from an adjacent utility‐scale wind turbine. In the field measurements of the utility‐scale turbine, discrete combinations of speed and direction shear induce changes in power production of −19% to +34% relative to the turbine power curve for a given hub height wind speed. Positive speed shear generally corresponds to over‐performance and increasing magnitudes of direction shear to greater under‐performance, relative to the power curve. Overall, the blade element model produces both higher correlation and lower error relative to the other models, but its quantitative accuracy depends on induction and controller sub‐models. To further assess the influence of complex, non‐monotonic wind profiles, we also drive the models with best‐fit power law wind speed profiles and linear wind direction profiles. These idealized inputs produce qualitative and quantitative differences in power predictions from each model, demonstrating that time‐varying, non‐monotonic wind shear affects wind power production.

Energy & Fuels↗
DOE OSTI2026

Embedded Ferroelectric Nanoclusters can drive Polarization Reversal in a Non-Ferroelectric Polar Film via the Proximity Effect

Heterogeneous nucleation from defects dominates the electric field required for polarization switching of ferroelectrics. Here, we consider the switching of a nominally non-switchable polar thin film of AlN due to the proximity effect arising from embedded ferroelectric nanoclusters of Al 1-x Sc x N. Using a Landau-Ginzburg-Devonshire thermodynamic approach and finite element modeling, we study the influence of nanocluster shape on polarization switching and domain nucleation emerging in AlN. The ferroelectric nanocluster boundary is modeled as a thin layer transitioning from Al 1-x Sc x N to AlN. We analyze the conditions under which polarization switching in the AlN film occurs at coercive fields significantly lower than its dielectric breakdown field. In the presence of spike-like Al 1-x Sc x N nanoclusters, the proximity effect enables switching of the spontaneous polarization in AlN and significantly reduces the corresponding coercive field. The internal field, which is depolarizing inside the AlN (due to its larger spontaneous polarization) and polarizing within the ferroelectric Al 1-x Sc x N nanoclusters (due to its smaller spontaneous polarization), lowers the potential barrier in the clusters and nucleates nanodomains at the Al 1-x Sc x N-AlN interface, forming localized regions of reversed polarization. Proximity effect can thus provide a pathway towards "thawing" previously "frozen" ferroelectrics through engineered nucleation for memory, actuation and optical technologies.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Cluster Dynamics Simulations of Intra-Granular Fission Gas Bubble Size and Pressure Evolution in UO 2

Fission gases such as xenon (Xe) play a critical role in determining the behavior and response of nuclear fuel. Given that Xe has little solubility in UO 2 , it accumulates and forms bubbles, which significantly impact fuel performance. Intra- and inter-granular bubble nucleation and growth can lead to fuel swelling, and once bubbles interconnect at grain boundaries, fission gas can be released into the plenum. At low temperatures, limited uranium vacancy mobility can restrict swelling, therefore causing the bubbles to become highly pressurized. Consequently, this can induce micro-cracking, promote fission gas release (increasing the likelihood of cladding failure), and even lead to fuel pulverization under accident conditions such as a loss of coolant accident. As bubble evolution is strongly influenced by local temperature and fission rate, markedly different behavior occurs across the radial profile of the fuel pellet. Capturing the mechanisms that underpin bubble evolution is therefore important to predict these behaviors in the fuel. Previous models describing important mechanisms informed by lower length scale simulations have been developed under the NEAMS program. These can describe the evolution of a single bubble type (i.e., single value for radius and pressure) at each position in the pellet, for instance using the Centipede cluster dynamic code. However, in reality, a full distribution in bubble sizes and pressures exists within the microstructure at a given position in the pellet. To address this the cluster dynamics code Xolotl, which can predict Xe and vacancy phase space (i.e., bubble distributions) for intra-granular bubbles, has been used before. Prior work benchmarked the Xolotl code against the Centipede cluster dynamics code to ensure compatibility and to verify that mobile defect properties are adequately transferred between the two codes, along with some physics improvements. In this work, we go further by introducing a physics-based set of improvements that will allow us to accurately predict bubble size distributions and internal bubble pressures under representative UO 2 irradiation conditions. The improvements include (i) coupling bubble-defect reaction energies to a virial equation of state (EOS), (ii) including a bubble surface tension contribution, (iii) incorporating radiation-induced re-solution of Xe and vacancies, (iv) enabling pressure-driven dislocation loop punching through an effective emission of interstitial clusters informed by interstitial loop energetics, (v) accounting for radiation induced athermal diffusion of Xe, and (vi) implementing a Booth-type grain boundary sink representation for all mobile defects and defect clusters. After these modifications, we observe good agreement of Xolotl fission gas bubble size and concentration predictions with legacy experimental measurements. Additionally, it allows the distribution of Xe bubble pressures and radius to also be predicted and compared to data produced through the Advanced Fuels Campaign (AFC) program. Here, we have done this by running simulations under conditions similar to the AFC post-irradiation examination (PIE) samples irradiated at North Anna 2 light water reactor (LWR). Our results shows excellent agreement with these experimental measurements.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗