Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “information diffusion”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Pulsed Thermal Tomography Nondestructive Examination of Additively Manufactured Reactor Materials (Second Annual Progress Report)

Additive manufacturing (AM) is an emerging method for cost-efficient fabrication of nuclear reactor parts. AM of metallic structures for nuclear energy applications is currently based on laser powder bed fusion (LPBF) process, which can introduce internal material flaws, such as pores and anisotropy. Integrity of AM structures needs to be evaluated nondestructively because material flaws could lead to premature failures due to exposure to high temperature, radiation and corrosive environment in a nuclear reactor. Thermal tomography (TT) provides a capability for non-destructive evaluation of sub-surface defects in arbitrary size structures. Thermal tomography is a computational method for heat diffusion-based imaging of solids, which provides 3D visualization of data from flash thermography measurements. We investigate thermal tomography imaging and nondestructive evaluation of stainless steel and nickel super alloy metallic structures produced with laser powder bed fusion (LPBF) additive manufacturing (AM) process. Metallic structures produced with LPBF contain defects, and there are limited capabilities to evaluate these structures non-destructively. Thermal tomography reconstruction of 3D apparent spatial effusivity provides information about AM structure geometry and internal material flaws. We study performance of thermal tomography in imaging of metallic structures through COMSOL computer simulations of transient heat transfer, and through reconstruction of data obtained from experimental measurements. Reconstruction of internal defects is investigated using a stainless steel 316L specimen with flat bottom hole (FBH) indentations, and Inconel 718 plate produced with laser powder bed fusion (LPBF) method, which contains imprinted hemispherical shape low density regions containing non-sintered metallic powder. The FBH’s have the same sizes as the imprinted defects in the LPBF specimens, but offer better imaging contrast. Thermal tomography reconstructions provide visualizations of internal defects, and allow for estimation of their sizes and locations. Detection sensitivity of TT is limited by noises. We investigate separation of signal from noise in thermography images using several machine learning (ML) methods, including new spatio-temporal blind source separation (STBSS) and spatio-temporal sparse dictionary learning (STSDL) methods. Performance of the ML methods is benchmarked using thermography data obtained from imaging stainless steel 316L and Inconel 718 specimens produced LPBF method with imprinted calibrated porosity defects. The ML methods are ranked by F-score and execution runtime. Finally, we investigate TT of AM stainless steel 316L specimen with imprinted internal porosity defects using relatively low-cost, small form factor infrared (IR) camera based on uncooled micro bolometer detector. Sparse coding related K-means singular value decomposition (SVD) machine learning, image processing algorithms are developed to improve quality of TT images through removal of Additive white Gaussian noise without blurring the images. Following initial qualification of an AM component for deployment in a nuclear reactor, a compact TT can also be used for in-service nondestructive evaluation (NDE). With capability to perform in-service NDE of the AM component lifecycle, TT data can be used for development of a component digital twin.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Predictive Large-Eddy Simulation of Supercritical-Pressure Reactive Flows in the Cold Ignition Regime

This report describes a one-year study performed under DOE sponsorship, continuing the investigation of high-pressure turbulent reactive flows. The interest is in the effect of the chemical species distribution in high-pressure turbulent flows in the presence of strong temperature gradients as would occur during reactions in realistic flows. The prime example of such flows are boundary layers in which the wall is at a lower temperature than that of the fluid, as would be the case in Diesel engines. Previous DOE BES work in the program further highlighted the importance of the boundary layer as a configuration for fundamental studies: soot formation in boundary layers is still a problem poorly understood and depends on the availability of particular chemical species at that location, and when boundary layer Large Eddy Simulations results were compared with measurements, the agreement was unfavorable, showing that this important ‘unit’ problem is not well understood. To understand this unit problem, modeling and Direct Numerical Simulations of this unit problem were conducted for the simplest possible multispecies system, that is a binary-species system. The results discovered a new phenomenon, that is, Soret effect induced uphill diffusion. The far-reaching implication is that through an imposed wall/fluid temperature difference it is possible to control the distribution of the species in the boundary layer. Because the results have been documented in one paper published in the refereed literature, and also in conference papers, this final report is in the form of an Executive Summary succinctly describing the results and putting them in perspective with respect to existing information. The refereed and conference papers published are individually listed as Appendices and attached to this report. One manuscript is still in progress and is thus not listed.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Computational materials reliability assessment of hydrogen fueled gas turbine power generation engines

The use of blended fuel sources in land based gas turbine engines drives variations in the resulting operational profile (temperatures and pressures) which can impact engine reliability. Furthermore, variability in the manufacture of components affects the resulting microstructure which directly impacts material performance and reliability. Currently, data-driven models are typically used for maintaining and inspecting fleets of engines. Without explicitly capturing material and operational sources of variability conservatism must be used in developing component-level reliability models. Therefore, there exists an opportunity to use information from materials-scale physics models to better inform reliability modeling and reduce conservatism; the impact is more cost-efficient operation and maintenance of current and future fleets. Specifically, this work establishes a computational framework for evaluating the probabilistic high temperature creep performance of hot-section Ni-based superalloys where uncertainty comes from both microstructural and operational variability. A novel high-fidelity physics model which phenomenologically captures grain-boundary sensitive phenomena has been established. A probabilistic calibration procedure was used to calibrate the model and capture uncertainty in the parameterized model coefficients. A design of experiments methodology was established for identifying informative microstructural digital representations for suitable for forward model evaluation. Results show that training a machine-learning surrogate using this design criteria outperforms random selection of microstructural representations. Finally, two surrogate models were developed: (1) a deterministic surrogate model which predicts the local field response given microstructure, constitutive model parameters, and operating conditions (stress, temperature) and (2) a probabilistic model, where uncertainty comes from constitutive law uncertainty, built using denoising diffusion probabilistic models which samples responses given (1) microstructure and (2) operating conditions. These surrogate models enable partner Siemens Energy to rapidly perform UQ analysis specific to creep deformation across a range of microstructures and operating conditions. The impact is that these ML and physics codes can be used to establish more advanced reliability models for the inspection, servicing, and maintenance of land based gas turbine engines.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Modeling and simulation to investigate the electrification potential of medium- and heavy-duty vehicle fleets

This project involves developing and integrating new modeling tools to simulate the dynamics of electric medium- and heavy-duty fleet vehicle adoption. A technical and economic modeling tool, combining a data-driven hardware cost model with a cost-optimal charging strategy microsimulation, enables tailored analysis of the costs and benefits of electrifying individual fleets. Next, a novel text synthesis process, applied to a curated corpus of literature, quantifies trade-offs between technical, economic, and other factors in the fleet vehicle procurement decision. The outcomes of these tasks combine with knowledge from recent literature on fleet decision processes to specify the vehicle procurement model used by fleets in an agent-based model of the medium- and heavy-duty electric vehicle market. This model embodies an especially disaggregated approach to adoption modeling, internalizing factors and dynamics that conventional adoption models externalize. In particular, explicitly modeling the formation and diffusion of opinions among agents enables experiments that conventional models cannot support. Demonstrations show, for example, that increasing the extent of interactions between populations with different proclivities to electric vehicles has an asymmetrical outcome. High-proclivity electric vehicle adoption is generally unaffected as interactions increase, but low-proclivity adoption is accelerated. By representing individual fleets' requirements and costs at a high level of detail, incorporating an adoption decision model informed by a wide body of empirical research, and broadening the array of variables and dynamics available for experimentation, this integrated model offers a new way to understand the urgent challenge of eliminating emissions from the most emissions-intensive transportation sectors.

Trinko, David A.↗
DOE OSTI2025

Dynamical Signatures of Thermotoga maritima Maltose-Binding Proteins Affected by Ligand Binding

Functional segregation among protein isoforms depends on the interplay of their overall structures and the molecular dynamics of these structures. Thermotoga maritima maltose-binding protein (tmMBP) isoforms show size-dependent differential binding of maltose and malto-oligomers while maintaining remarkable fold conservation. This differential behavior needs detailed characterization in native-like aqueous conditions to understand the effects of protein dynamics on ligand binding and recognition. Small-angle neutron scattering (SANS), neutron spin echo (NSE) spectroscopy, and dynamic light scattering (DLS) were used in conjunction with previously published computational molecular dynamics (MD) simulations to understand the dynamic behavior of tmMBPs experimentally. SANS provided information on the overall structure of the molecules, while NSE was used to determine the dynamics in the nanosecond time scale. Both tmMBP2 and tmMBP3 have a bidomain architecture linked with a flexible hinge, with the binding pocket sitting in the cleft between the two domains. tmMBP2 and tmMBP3 showed different solution dynamics, with the translational and rotational components dominating the dynamics of both systems, resulting in a clear differentiation of their diffusion pattern. A faster dynamics component was also observed and was attributed to segmental dynamics. Differences observed between the ligand-free (apo) and ligand-bound (holo) states of the two proteins are attributed to conformational entropy. Our results highlight the intricacies of how structure and dynamics can together shape binding to a repertoire of substrates in structurally similar proteins.

Diffusion↗
DOE OSTI2020

Thermal tomography 3D imaging of additively manufactured metallic structures

Thermal tomography is a computational method for heat diffusion-based imaging of solids, which provides 3D visualization of data from flash thermography measurements. We investigate thermal tomography imaging and nondestructive evaluation of stainless steel and nickel super alloy metallic structures produced with the laser powder bed fusion (LPBF) additive manufacturing (AM) process. Metallic structures produced with LPBF contain defects, and there are limited capabilities to evaluate these structures non-destructively. Thermal tomography reconstruction of 3D apparent spatial effusivity provides information about AM structure geometry and internal material flaws. We study performance of thermal tomography in imaging of metallic structures through COMSOL computer simulations of transient heat transfer and through reconstruction of data obtained from experimental measurements. Thermal tomography reconstructions of structure shape and dimensions are shown for the Inconel 718 AM structure which has variations in the horizontal plane but is uniform along the depth dimension. Reconstruction of internal defects is investigated using a stainless steel 316L specimen with flat bottom hole (FBH) indentations, and the Inconel 718 plate is produced with the LPBF method, which contains imprinted hemispherical shape low density regions containing non-sintered metallic powder. The FBHs have the same sizes as the imprinted defects in the LPBF specimens but offer better imaging contrast. Thermal tomography reconstructions provide visualizations of internal defects and allow for estimation of their sizes and locations. Results of this study demonstrate that thermal tomography can be used for visualization and quality control in AM.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Water, Ion, and Proton Dynamics and Interactions in Nanoconfined Systems

The research was directed at understanding the role that nanoconfinement has on the dynamics and structure of liquids confined in nanometer size structures. Experiments were conducted using ultrafast infrared methods. The methods were two dimensional infrared (2D IR) spectroscopy and polarization pump-probe (PSPP) experiments. In addition FT-IR spectroscopy was used. To augment the experimental work, high level molecular dynamics simulations of the systems were conducted. The system used for confinement were nanoporous silica with essentially monodispersed pores. Pores of different diameters were studies, i.e., 2.5 nm to 8 nm. Water and organic liquids were investigated. The experiments were performed on a vibrational probe molecule, i.e., the CN stretch of the anion SeCN-. The 2D IR experiments measured spectral diffusion, which directly reported on the structural dynamics of the confined liquids. 2D IR was also used to perform chemical exchange measurements, which observed the formation and dissociation of molecular complexes with surface silanols. The PSPP experiments measured orientational relaxation of the probe molecule as well as it vibrational lifetime. All of the observables depended on the diameter of the pores and the type of liquid. The simulations of the observables provided detailed structural information on the structure of the confined liquids and there interactions with the surfaces. In addition to the study of the confinement in nanopores, water in hydrogels was also studied. The same experimental methods were used along with MD simulations. It was found the influence of hydrogels on water did not come from confinement of the water in the hydrogel polymer network, but rather from the interactions of the water with the individual monomers that compose the polymer. In addition, the temperature dependence of a polymer/water system was tracked from below to above the lower critical solution temperature to observe the nature of the phase separated polymer and how water responded to the phase separation.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Conversion between triplet pair states is controlled by molecular coupling in pentadithiophene thin films

In singlet fission (SF) the initially formed correlated triplet pair state, 1 (TT), may evolve toward independent triplet excitons or higher spin states of the (TT) species. The latter result is often considered undesirable from a light harvesting perspective but may be attractive for quantum information sciences (QIS) applications, as the final exciton pair can be spin-entangled and magnetically active with relatively long room temperature decoherence times. In this study we use ultrafast transient absorption (TA) and time-resolved electron paramagnetic resonance (TR-EPR) spectroscopy to monitor SF and triplet pair evolution in a series of alkyl silyl-functionalized pentadithiophene (PDT) thin films designed with systematically varying pairwise and long-range molecular interactions between PDT chromophores. The lifetime of the (TT) species varies from 40 ns to 1.5 μs, the latter of which is associated with extremely weak intermolecular coupling, sharp optical spectroscopic features, and complex TR-EPR spectra that are composed of a mixture of triplet and quintet-like features. On the other hand, more tightly coupled films produce broader transient optical spectra but simpler TR-EPR spectra consistent with significant population in 5 (TT) 0 . These distinctions are rationalized through the role of exciton diffusion and predictions of TT state mixing with low exchange coupling J versus pure spin substate population with larger J. The connection between population evolution using electronic and spin spectroscopies enables assignments that provide a more detailed picture of triplet pair evolution than previously presented and provides critical guidance for designing molecular QIS systems based on light-induced spin coherence.

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

Formation and Microwave Losses of Hydrides in Superconducting Niobium Thin Films Resulting from Fluoride Chemical Processing

Abstract Superconducting niobium (Nb) thin films have recently attracted significant attention due to their utility for quantum information technologies. In the processing of Nb thin films, fluoride‐based chemical etchants are commonly used to remove surface oxides that are known to affect superconducting quantum devices adversely. However, these same etchants can also introduce hydrogen to form Nb hydrides, potentially negatively impacting microwave loss performance. Here, comprehensive materials characterization of Nb hydrides formed in Nb thin films as a function of fluoride chemical treatments is presented. In particular, secondary‐ion mass spectrometry, X‐ray scattering, and transmission electron microscopy reveal the spatial distribution and phase transformation of Nb hydrides. The rate of hydride formation is determined by the fluoride solution acidity and the etch rate of Nb 2 O 5, which acts as a diffusion barrier for hydrogen into Nb. The resulting Nb hydrides are detrimental to Nb superconducting properties and lead to increased power‐independent microwave loss in coplanar waveguide resonators. However, Nb hydrides do not correlate with two‐level system loss or device aging mechanisms. Overall, this work provides insight into the formation of Nb hydrides and their role in microwave loss, thus guiding ongoing efforts to maximize coherence time in superconducting quantum devices.

36 MATERIALS SCIENCE↗
DOE OSTI2024

From Bulk to Interface: Solvent Exchange Dynamics and Their Role in Ion Transport and the Interfacial Model of Rechargeable Magnesium Batteries

Multivalent battery chemistries have been explored in response to the increasing demand for high-energy rechargeable batteries utilizing sustainable resources. Solvation structures of working cations have been recognized as a key component in the design of electrolytes; however, most structure–property correlations of metal ions in organic electrolytes usually build upon favorable static solvation structures, often overlooking solvent exchange dynamics. We here report the ion solvation structures and solvent exchange rates of magnesium electrolytes in various solvents by using multimodal nuclear magnetic resonance (NMR) analysis and molecular dynamics/density functional theory (MD/DFT) calculations. These magnesium solvation structures and solvent exchange dynamics are correlated to the combined effects of several physicochemical properties of the solvents. Moreover, Mg 2+ transport and interfacial charge transfer efficiency are found to be closely correlated to the solvent exchange rate in the binary electrolytes where the solvent exchange is tunable by the fraction of diluent solvents. Our primary findings are (1) most battery-related solvents undergo ultraslow solvent exchange coordinating to Mg 2+ (with time scales ranging from 0.5 μs to 5 ms), (2) the cation transport mechanism is a mixture of vehicular and structural diffusion even at the ultraslow exchange limit (with faster solvent exchange leading to faster cation transport), and (3) an interfacial model wherein organic-rich regions facilitate desolvation and inorganic regions promote Mg 2+ transport is consistent with our NMR, electrochemistry, and cryogenic X-ray photoelectron spectroscopy (cryo-XPS) results. In conclusion, this observed ultraslow solvent exchange and its importance for ion transport and interfacial properties necessitate the judicious selection of solvents and informed design of electrolyte blends for multivalent electrolytes.

25 ENERGY STORAGE↗
DOE OSTI2022

Self-mixing in microtubule-kinesin active fluid from nonuniform to uniform distribution of activity

Active fluids have applications in micromixing, but little is known about the mixing kinematics of systems with spatiotemporally-varying activity. To investigate, UV-activated caged ATP is used to activate controlled regions of microtubule-kinesin active fluid and the mixing process is observed with fluorescent tracers and molecular dyes. At low Péclet numbers (diffusive transport), the active-inactive interface progresses toward the inactive area in a diffusion-like manner that is described by a simple model combining diffusion with Michaelis-Menten kinetics. At high Péclet numbers (convective transport), the active-inactive interface progresses in a superdiffusion-like manner that is qualitatively captured by an active-fluid hydrodynamic model coupled to ATP transport. Results show that active fluid mixing involves complex coupling between distribution of active stress and active transport of ATP and reduces mixing time for suspended components with decreased impact of initial component distribution. This work will inform application of active fluids to promote micromixing in microfluidic devices.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Reduced models for ETG transport in the tokamak pedestal

This paper reports on the development of reduced models for electron temperature gradient (ETG) driven transport in the pedestal. Model development is enabled by a set of 61 nonlinear gyrokinetic simulations with input parameters taken from pedestals in a broad range of experimental scenarios. The simulation data have been consolidated in a new database for gyrokinetic simulation data, the multiscale gyrokinetic database (MGKDB), facilitating the analysis. The modeling approach may be considered a generalization of the standard quasilinear mixing length procedure. The parameter η, the ratio of the density to temperature gradient scale length, emerges as the key parameter for formulating an effective saturation rule. With a single order-unity fitting coefficient, the model achieves an error of 15%. A similar model for ETG particle flux is also described. We also present simple algebraic expressions for the transport informed by an algorithm for symbolic regression.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Enhanced accuracy through ensembling of randomly initialized auto-regressive models for dynamical systems

Computational mechanics simulations using traditional finite element methods (FEM) require prohibitively expensive computational resources for real-time engineering applications, design optimization, and digital twin implementations. While machine learning (ML) surrogate models offer significant computational speedups, autoregressive ML models for time-dependent mechanical systems suffer from error accumulation that compromises long-term prediction reliability - a critical concern for engineering applications where accuracy over extended time horizons is essential for safety and performance assessments. Here, we propose a deep ensemble framework specifically designed to address this challenge in computational mechanics applications, where multiple ML surrogate models with random weight initializations are trained in parallel and their predictions aggregated during inference. This approach leverages statistical diversity to maximize information gain from a fixed set of training data and to mitigate error propagation, while maintaining the computational efficiency that makes ML surrogates attractive for engineering practice. We validate the framework on three representative problems spanning critical areas of computational mechanics: stress field evolution in heterogeneous microstructures under complex loading (relevant to advanced materials design and composite analysis), planetary-scale shallow water dynamics (applicable to environmental and geotechnical engineering), and Gray-Scott reaction-diffusion systems (relevant to mass transport and chemical process engineering). Across all test cases, the ensemble approach demonstrates consistent error reduction of 15-33% compared to individual models. The codes for this work are available on GitHub (https://github.com/Graham-Brady-Research-Group/AutoregressiveEnsemble_SpatioTemporal_Evolution).

autoregressive prediction↗
DOE OSTI2022

Understanding the fundamental driver of semiconductor radiation tolerance with experiment and theory

Space is the operating environment of a multitude of systems that our society is heavily reliant on today, however, maintaining operability there necessitates special consideration of the electronic systems' tolerance of space radiation. Electronic systems are critically dependent on the electronic properties of their semiconductor components, which are modified by space radiation with an adverse impact on the space system performance. What innate property allows some semiconductors to sustain little damage while others accumulate defects rapidly with dose is poorly understood, which limits the extent to which radiation tolerance can be implemented as a design criterion. Here, to gain insight into what properties are drivers of semiconductor radiation tolerance, the first step is to generate a dataset of the relative radiation tolerance of a broad sampling of semiconductors. To accomplish this, Rutherford backscatter channeling experiments are used to compare the displaced lattice atom buildup in InAs, InP, GaP, GaN, ZnO, MgO, and Si as a function of stepwise alpha particle dose. With this experimental information on radiation-induced incorporation of interstitial defects in hand, hybrid density functional theory electron densities (and their derived quantities) are calculated and their gradient and Laplacian are evaluated to obtain key fundamental information about the interactions in each material. It is shown that simple, undifferentiated values (which are typically used to describe bond strength) are insufficient to predict radiation tolerance. Instead, the curvature of the electron density at bond critical points provides a measure of radiation tolerance consistent with the experimental results obtained. This curvature and associated forces surrounding bond critical points have the potential to disfavor the localization of displaced lattice atoms at these points, favoring their diffusion toward perfect lattice positions. With this criterion to predict radiation tolerance, simple density functional theory simulations can be conducted on potential new materials to gain insight into how they may operate in demanding high radiation environments.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Bayesian Analysis of TRISO Fuel: Quantifying Model Inadequacy, Incorporating Lower-Length-Scale Effects, and Developing Parallel Active Learning Capabilities

The U.S. Department of Energy (DOE)’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program aims to develop predictive capabilities by applying computational methods to the analysis and design of advanced reactor and fuel-cycle systems. This program has been providing engineering-scale support for the continued development of BISON, a high-fidelity, high-resolution fuel performance tool. Fuel behavior in nuclear reactors is governed by a complex network of mechanisms that interact with various other physics aspects in the reactor system. Any model developed to represent fuel behavior will likely be idealized, resulting in uncertainties when comparing their predictions against the observed data. In Fiscal Year (FY)-23, we initiated the Uncertainty Quantification (UQ) work by using Bayesian methods to establish a level of model trustworthiness and further improve it, with a particular emphasis on TRI-Structural isOtropic (TRISO) nuclear fuel. This year, we further expanded on that UQ work by investigating an approach to quantifying model inadequacy and accounting for lower-length scale (LLS) effects in TRISO silver (Ag) release modeling. Furthermore, we are implementing parallel active learning capabilities to reduce the computational cost (i.e., required computational resources and elapsed time) of performing UQ. Specifically, we utilized The Kennedy O’Hagan framework for Bayesian uncertainty quantification (KOH) to account for model inadequacy in TRISO Ag release predictions made by BISON. The KOH framework represents an improvement over the standard Bayesian framework used in FY-23. Explicitly accounting for model inadequacy in the Bayesian framework helps establish the level of experimental noise uncertainty in the Advanced Gas Reactor (AGR) data. We compared the inverse UQ results obtained from both the standard Bayesian and KOH frameworks in light of the AGR-2/3/4 data, and also compared the predictive UQ results obtained from these two frameworks in light of the AGR-1 data. Next, we investigated the impact of considering LLS effects in the Ag release simulations. We developed an expanded database of LLS simulated effective diffusivities for Ag, covering a wide range of microstructures and temperatures. Using this database, we developed a framework for incorporating LLS effects into the engineering-scale Ag release UQ. We developed both parametric and non-parametric approaches for bridging the length scales. We then investigated the inverse UQ results in light of the AGR-2/3/4 data and the predictive UQ results in light of the AGR-1 data, and compared the LLS-informed approach and the Arrhenius equation, which does not include microstructure information. Finally, we discussed implementing parallel active learning capabilities in the Multiphysics Object Oriented Simulation Environment (MOOSE)/BISON to reduce the computational cost (i.e., computational resources and elapsed time) of Bayesian UQ. For verification purposes, we first tested these new capabil ities on a species interaction problem. We then demonstrated them on the TRISO Ag release application, showing that parallel active learning capabilities can enhance the accuracy of UQ while also substantially reducing the computational cost in comparison to the reference methods developed in FY-23.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Templates for Risk Informed Assurance with Curvature Embeddings (TRACE)

We investigate recovery of geometric structure from networks embedded in manifolds with spatially varying curvature, extending the constant-curvature framework of Lubold et al. (2023). Our work supports cascade risk assessment in critical infrastructure through the Templates for Risk-informed Assurance with Curvature Embeddings (TRACE) framework. Simulations on a bi-modal Gaussian surface show that constant-curvature methods yield weighted averages shaped by clique patterns, while hierarchical clustering identifies distinct regimes. Localized estimation, however, reveals boundary contamination in transitional regions. To address heterogeneity, we develop distance metrics for graphs with edge and node features, proving their metric validity, and validate them via deterministic graph generation from canonical tilings. We further propose a diffusion-based anomaly detection approach that treats networks as glued manifolds, using curvature discontinuities to detect structural anomalies. Employing the carré-du-champ operator and scalar curvature, we achieve robust anomaly discrimination, demonstrated on the Singapore Water Treatment (SWaT) dataset with joint network-traffic and sensor features. Integration with TRACE reveals how curvature shapes cascade dynamics: positive curvature impedes, while negative curvature accelerates propagation. This geometric perspective provides interpretable risk metrics and visualization tools for critical infrastructure managers. While full validation remains ongoing, our contributions establish a rigorous foundation for geometric analysis of network resilience and cascade vulnerability.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Atomistic simulations and machine learning of solute grain boundary segregation in Mg alloys at finite temperatures

Understanding solute segregation thermodynamics is the first step in investigating grain boundary (GB) properties, such as strong yttrium (Y) effects on grain growth and texture evolution in micro-scale polycrystalline magnesium (Mg) alloys. To estimate the average GB segregation behavior in low-solute-concentration Mg alloys (e.g., 2 at.% Y), a state-of-the-art spectral approach is applied based on a per-site segregation energy spectrum for Y solute atoms at zero K obtained from molecular statistics (MS) simulations of ~10 4 GB sites in Mg symmetric tilt GBs (STGBs). Although selected MS simulation results are consistent with verification by density functional theory (DFT) calculations, estimates of average segregation tendency based on the zero-K energy spectrum deviate from experimental observations. To resolve this problem, thermodynamic integration (TI) methods based on molecular dynamics (MD) simulations are used to determine the per-site segregation free energies of Y at representative GB sites, which show contributions beyond harmonic approximations can be important for certain GB sites at high temperatures. A surrogate model of per-site segregation free energy is constructed from a small set of TI data points using stacking cross-validation regressors and physics-informed descriptors. This model is applied to predict the Y segregation free energy spectra for thousands of GB sites in Mg STGBs with uncertainty quantification. Finally, the average segregation tendency predicted by the spectral approach based on the free energy spectra agrees well (within the uncertainty range) with experimental observations for micro-scale polycrystalline Mg alloys at typical thermomechanical processing temperatures (500 ~ 800 K), where thermodynamic equilibrium states are likely to be achieved due to fast diffusion.

Atomistic simulations↗
DOE OSTI2026

The impact of chemistry on anion migration in bixbyite-structured lanthanide oxides

This manuscript describes atomistic calculations of oxygen vacancy and interstitial migration in bixbyite structured lanthanide oxides. We examine two types of compounds, one in which only one type of lanthanide cation is present and a second class in which two lanthanides are present in a 3:1 ratio as dictated by the symmetry of the bixbyite crystal structure. Using temperature accelerated dynamics and the nudged elastic band method, we quantify the role of chemistry on the energy barriers for the most important pathways for both vacancy and interstitial migration. We then analyze the impact of these variations on the overall diffusivity of each defect, quantifying the contribution of each pathway using the theory of kinosons. We find that vacancy mobility can vary by as much as three orders of magnitude through changes in chemistry at 500 K. Changes in interstitial mobility are more modest but can still vary by an order of magnitude. This points to the ability to tune the mass transport characteristics of these compounds through appropriate choices in chemistry. We have also included supplementary information containing the atomic structures of the relevant pathways.

36 MATERIALS SCIENCE↗