Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interpretable models”

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 1,351 records · Page 75

Distinguishing Electron and Hole Dynamics in Functionalized CdSe/CdS Core/Shell Quantum Dots Using Complementary Ultrafast Spectroscopies and Kinetic Modeling

The evolution of excitation energy and photogenerated charges in semiconductor quantum dots (QDs) functionalized with molecular acceptors can be probed on ultrafast time scales using techniques such as transient absorption (TA) spectroscopy. However, historical interpretations that the 1S(e)-1S 3/2 (h) transition in Cd-chalcogenide QDs is fully attributable to electrons may be misleading, and multiexponential models used to fit TA kinetics do not correspond directly to specific photophysical processes. Here, we present visible-wavelength and mid-IR TA and time-resolved photoluminescence measurements to inform a comprehensive kinetic model of the photoexcited CdSe/CdS core/shell QDs functionalized with passivating oleic acid (OA), hole-accepting ferrocene, or electron-accepting naphthalene bisimide (NBI). We show that ~30% of the 1S signal and 72% of the IR signal can originate from holes in well-passivated core/shell QDs. We also demonstrate evidence of electron trapping in OA-capped core/shell QDs, with additional electron transfer and hole trapping in the QDs functionalized with NBI. Electron (hole) trapping and detrapping occur in 450 ± 100 ps (430 ± 70 ps) and 340 ± 100 ps (1.1 ± 0.4 ns) respectively, while the time constant for electron transfer to NBI is ~1.8 ns. The comprehensive picture of photophysical processes provided by the complementary ultrafast techniques and kinetic modeling can accelerate both the fundamental science and application development of nanostructured and molecular systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Initial Fracture Propagation Modeling of Graphite Components with Grizzly

Graphite has historically been extensively used in power reactor cores and will be used in multiple types of advanced reactors currently under development. These graphite structural components can experience significant stresses due to nonuniform volumetric strains induced by irradiation and thermal expansion, which can lead to fracture. Robust tools for predicting fracture initiation and propagation in graphite structural components in nuclear reactors are important for evaluating component integrity, developing design standards, and interpreting experimental results to characterize graphite performance. The U.S. Department of Energy’s Nuclear Energy Advanced Modeling and Simulation program has been developing degradation models for other structural components in nuclear reactors within the Grizzly and BlackBear codes. This report documents an effort to develop initial capabilities for modeling graphite fracture within these codes, building on prior efforts to model fracture in other materials. Major elements of this effort include developing a new system for modeling fracture nucleation and growth in two dimensions using the extended finite element method and incorporating a damage and plasticity model. These capabilities are applied here to model a representative graphite component and a splitting disc experiment used to obtain tensile strength.

36 MATERIALS SCIENCE↗

Fiber Bundle Morphisms as a Framework for Modeling Many-to-Many Maps

While it is not generally reflected in the `nice' datasets used for benchmarking machine learning algorithms, the real-world is full of processes that would be best described as many-to-many. That is, a single input can potentially yield many different outputs (whether due to noise, imperfect measurement, or intrinsic stochasticity in the process) and many different inputs can yield the same output (that is, the map is not injective). For example, imagine a sentiment analysis task where, due to linguistic ambiguity, a single statement can have a range of different sentiment interpretations while at the same time many distinct statements can represent the same sentiment. When modeling such a multivalued function $f: X \rightarrow Y$, it is frequently useful to be able to model the distribution on $f(x)$ for specific input $x$ as well as the distribution on fiber $f^{-1}(y)$ for specific output $y$. Such an analysis helps the user (i) better understand the variance intrinsic to the process they are studying and (ii) understand the range of specific input $x$ that can be used to achieve output $y$. Following existing work which used a fiber bundle framework to better model many-to-one processes we describe how morphisms of fiber bundles provide a template for building models which naturally capture the structure of many-to-many processes.

Coda, Elizabeth D.↗

Beyond Energy Efficiency: A clustering approach to embed demand flexibility into building energy benchmarking

The intermittency of carbon-free renewables and the demand changes associated with the widespread push for electrifying the transportation and building sectors provides an opportunity for buildings to go beyond energy efficiency and push towards providing demand flexibility to the electricity grid. The duality of energy efficiency and demand flexibility is necessary for success in a sustainable and reliable energy transition. Current building energy benchmarking models are limited in their ability to integrate concepts of demand flexibility and/or utilize granular smart meter data. Thus, current benchmarking methods are focused annual energy usage and fail to incorporate how the time of use of energy consumption impacts emissions in a quickly changing energy grid. Without a more comprehensive view of energy usage and associated real-time emissions, current benchmarking methods are unlikely to realize the full decarbonization potential of buildings. New emerging data streams provide an opportunity to develop a new generation of benchmarking energy models that embed dimensions of energy efficiency, grid interactivity, and demand flexibility into their analysis. In this paper, we propose a four-step method for embedding grid interactivity and demand flexibility into building benchmarking models that utilizes emerging building and time-series electricity data streams. We first engineer features to produce a mix-type dataset that encompasses many attributes of grid-interactive and efficient buildings, and then we apply K-medoids using Gower's Distance to produce peer-group clusters. We apply the method to a case study of 306 primary and secondary schools in southern California, USA. The results show that the method effectively clusters buildings by attributes of demand flexibility and energy efficiency. The clustering results reveal patterns in inefficient building operations and demand inflexibility at the building peer group level. In conclusion, the interpretation of clusters can serve as an integrated energy efficiency and demand flexibility benchmarking model and inform performance-specific policy targeting for buildings that go beyond traditional efficiency measures.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Graph network heterogeneity predicts interplant wake losses

Wind plants generate large-scale wakes, which can affect the performance of neighboring installations. Such wakes are challenging to estimate due to the inherent complexity in modeling wake interactions between large quantities of turbines at various distances. Weighted directed graph networks can inform complex models by linking turbine pairs into chains of upstream and downstream neighbors for a given wind direction. A novel interpretation of the graph network adjacency matrix is proposed where each element of the matrix represents the cumulative impact of upstream turbines on an individual. In this study, wake losses were estimated with an engineering wake model across a range of inflow conditions for nine parametric variations of a system containing two neighboring wind plants. The parametric nature of the study isolates turbine spacing within the plant, separation distance between plants, and wind direction as the main drivers of wake losses. Spatial heterogeneity is computed from the weighted average adjacency matrix of each plant arrangement. The proposed method is orders of magnitude faster than wake modeling and does not require detailed turbine information or atmospheric conditions. Furthermore, the weighted average adjacency matrix provides insight on the spatial organization of wake losses at various scales. Plant heterogeneity is correlated with wake losses within and among plants. Framing wind plant wake interaction in terms of graph network spatial heterogeneity provides an efficient approach for predicting wake losses within and among neighboring wind plants with applications to other complex systems where wake interactions are key factor.

17 WIND ENERGY↗

Anomalously abrupt switching of wurtzite-structured ferroelectrics: simultaneous non-linear nucleation and growth model

Ferroelectric polarization switching is one common example of a process that occurs via nucleation and growth, and understanding switching kinetics is crucial for applications such as ferroelectric memory. Here we describe and interpret anomalous switching dynamics in the wurtzite-structured nitride thin film ferroelectrics Al 0.7 Sc 0.3 N and Al 0.94 B 0.06 N using a general model that can be directly applied to other abrupt transitions that proceed via nucleation and growth. When substantial growth and impingement occur while nucleation rate is increasing, such as in these wurtzite-structured ferroelectrics under high electric fields, abrupt polarization reversal leads to very large Avrami coefficients ( e.g., n = 11), inspiring an extension of the KAI (Kolmogorov–Avrami–Ishibashi) model. We apply this extended model to two related but distinct scenarios that crossover between (typical) behavior described by sequential nucleation and growth and a more abrupt transition arising from significant growth prior to peak nucleation rate. This work therefore provides a more complete description of general nucleation and growth kinetics applicable to any system while specifically addressing the anomalously abrupt polarization reversal behavior in new wurtzite-structured ferroelectrics.

36 MATERIALS SCIENCE↗

Machine Learning an Ab-Initio Based Bond-Order Potential for Bismuthene

Bismuthene is a heavy 2D material whose strong spin–orbit coupling and recently observed single-element ferroelectricity have intensified interest in its structural, vibrational, and transport properties. Accurate modeling of these behaviors requires a short-range interatomic potential that can reproduce the underlying bonding physics at a fraction of the computational cost of first-principles methods. However, such a potential is currently unavailable. Here, in this work, we construct a Tersoff bond-order potential for β-bismuthene using a reinforcement-learning framework that integrates a continuous Monte Carlo Tree Search with a simplex-based local optimizer. The optimized parameter sets reproduce first-principles lattice constants, cohesive energy, the equation of state, elastic constants, and phonon dispersion. We validate the models by performing thermal-conductivity calculations and uniaxial fracture simulations our findings confirm the reliability of the resulting models across multiple thermomechanical regimes. Comparison of the three best solutions reveals how differences in pairwise interactions, angular terms, and bond-order behavior govern phonon features and mechanical responses. We demonstrate an interpretable and computationally efficient potential for bismuthene and demonstrate a general reinforcement-learning strategy for developing bond-order models in emerging 2D materials.

deformation↗

Ionization yield measurement in a germanium CDMSlite detector using photo-neutron sources

Two photo-neutron sources, Y 88 Be 9 and Sb 124 Be 9 , have been used to investigate the ionization yield of nuclear recoils in the CDMSlite germanium detectors by the SuperCDMS collaboration. In this study, we evaluate the yield for nuclear recoil energies between 1 and 7 keV at a temperature of ~ 50 mK . We use a geant4 simulation to model the neutron spectrum assuming a charge yield model that is a generalization of the standard Lindhard model and consists of two energy dependent parameters. We perform a likelihood analysis using the simulated neutron spectrum, modeled background, and experimental data to obtain the best fit values of the yield model. The ionization yield between recoil energies of 1 and 7 keV is shown to be significantly lower than predicted by the standard Lindhard model for germanium. There is a general lack of agreement among different experiments using a variety of techniques studying the low energy range of the nuclear recoil yield, which is most critical for interpretation of direct dark matter searches. This suggests complexity in the physical process that many direct detection experiments use to model their primary signal detection mechanism and highlights the need for further studies to clarify underlying systematic effects that have not been well understood up to this point.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Neutrino Flavor Model Building and the Origins of Flavor and CP Violation

The neutrino sector offers one of the most sensitive probes of new physics beyond the Standard Model of Particle Physics (SM). The mechanism of neutrino mass generation is still unknown. The observed suppression of neutrino masses hints at a large scale, conceivably of the order of the scale of a rand unified theory (GUT), which is a unique feature of neutrinos that is not shared by the charged fermions. The origin of neutrino masses and mixing is part of the outstanding puzzle of fermion masses and mixings, which is not explained ab initio in the SM. Flavor model building for both quark and lepton sectors is important in order to gain a better understanding of the origin of the structure of mass hierarchy and flavor mixing, which constitute the dominant fraction of the SM parameters. Recent activities in neutrino flavor model building based on non-Abelian discrete flavor symmetries and modular flavor symmetries have been shown to be a promising direction to explore. The emerging models provide a framework that has a significantly reduced number of undetermined parameters in the flavor sector. In addition, such a framework affords a novel origin of CP violation from group theory due to the intimate connection between physical CP transformation and group theoretical properties of non-Abelian discrete groups. Model building based on non-Abelian discrete flavor symmetries and their modular variants enables the particle physics community to interpret the current and anticipated upcoming data from neutrino experiments. Non-Abelian discrete flavor symmetries and their modular variants can result from compactification of a higher-dimensional theory. Pursuit of flavor model building based on such frameworks thus also provides the connection to possible UV completions: in particular, to string theory. We emphasize the importance of constructing models in which the uncertainties of theoretical predictions are smaller than, or at most compatible with, the error bars of measurements in neutrino experiments. While there exist proof-of-principle versions of bottom-up models in which the theoretical uncertainties are under control, it is remarkable that the key ingredients of such constructions were discovered first in top-down model building. We outline how a successful unification of bottom-up and top-down ideas and techniques may guide us towards a new era of precision flavor model building in which future experimental results can give us crucial insights into the UV completion of the SM.

Almumin, Yahya (ORCID:0000000196396772)↗

Effect of likelihood misspecification in Gaussian process-driven autonomous experimentation

In recent years, several groups have designed Autonomous Experiment (AE) models with the aim of using them as an alternative method for neutron scattering scanning. In an AE, Gaussian processes (GPs) are most frequently used due to their interpretability, their non-parametric nature, their universal approximation, and their closed-form predictive distribution. GPs have two key components, namely, the model for the likelihood of a neutron count knowing the underlying dynamic structure factor and the acquisition function. In this paper, we investigate the impact, on the quality of an AE, of the likelihood and acquisition function choices, in energy scans and (Q, ω) ones, with respect to the signal-over-noise ratio. While we hypothesized that the quality of GP predictions would decrease when the normal to Poisson likelihood approximation breaks down at low count rates, we found that the use of the correct Poisson likelihood does not improve the quality of the data collected, as well as yields very poor results in (Q, ω) scans at low count rates. In fact, the best results are obtained with a combination of normal likelihood, including the observation noise, and the change in variance acquisition function. In addition, we find that the performance, or quality of the predictive distribution, is a misleading measure of efficiency, that is, of the quality of the data collected.

Perryman, David Elliott [Inst. Laue-Langevin (ILL)↗

Measuring differential particle correlations in relativistic nuclear collisions

This study explores the transverse momentum ( p T ) dependencies of symmetric and asymmetric correlations (SC and ASC) with one and two particles of interest in Au + Au collisions at 200 GeV. Leveraging a multiphase transport (AMPT) model, the investigation delves into the sensitivity of these correlations to the final state effects, providing valuable insights into their potential for constraining the final state effects' p T dependencies. The HIJING model is employed as a benchmark for nonflow correlations, shedding light on their impact on interpreting SC and ASC data. Moreover, the study points out that differential SC and ASC with one and two particles of interest (POIs) typically incorporate contributions from event-plane angle fluctuations. Consequently, this work highlights the significance of SC and ASC with one and two POIs as valuable tools for investigating the p T nature of the final state effects and advocates for comprehensive experimental measurements across various beam energies and system sizes to enhance understanding and provide additional constraints for theoretical models. Published by the American Physical Society 2024

Physics↗

Simulation of Deuterium Oxidation in a Room with a Fire

A series of experiments where deuterium was released in trace amounts into a room with a fire were performed at Sandia National Laboratories. This report describes the corresponding effort to model the test series using SIERRA/FUEGO. The objective of this modeling effort was to produce a simulation test matrix that can be utilized to help interpret the corresponding experiments and be used to assess the credibility of using the SIERRA/FUEGO simulations as a surrogate for real tritium reaction data for mock fire and release studies. The simulations assessed multiple configurations of fire size and location and overall found that there was low conversion of deuterium to oxidized deuterium. There are some notable differences between the simulation predictions and the experimental results, such as the simulations generally predicting higher amounts of deuterium conversion compared to the experiment. A strategy to address modeling uncertainties and unknowns is also given in this report.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Stress intensity factor models using mechanics-guided decomposition and symbolic regression

The finite element method can be used to compute accurate stress intensity factors (SIFs) for cracks with complex geometries and boundary conditions. In contrast, handbook solutions act as surrogate SIF models that provide significantly faster evaluation times. However, the development of conventional surrogate SIF models relies on manual development based on low-order parameterizations. This limits surrogate model accuracy and generalizability. Here, in this paper, we develop a framework for the automated development of mechanics-guided handbook SIF solutions by using interpretable machine learning via genetic programming for symbolic regression (GPSR). Formalizing the mechanics-based approach of Raju and Newman, SIF training data is decomposed into multiple subsets. This decomposition enables parallel GPSR model development of subfunctions, each of which accounts for specific geometrical corrections with respect to a known analytical model. Using this mechanics-based approach with GPSR allows for equations to be learned with improved accuracy and reduced complexity relative to the Raju Newman equations while maintaining the inherent interpretability of mathematical expressions. In this paper, we present equations that match the complexity of the Raju Newman equations while having reduced error, as well as equations with similar errors and reduced complexity.

42 ENGINEERING↗

Hierarchical semi-Markov models with duration-aware dynamics for activity sequences

Residential electricity demand at granular scales is driven by what people do and for how long. Accurately forecasting this demand for applications like microgrid management and demand response therefore requires generative models for activities that can produce realistic daily activity sequences, capturing both the timing and duration of human behavior. This paper develops a generative model of human activity sequences using nationally representative time-use diaries at a 10-min resolution. We use this model to quantify which demographic factors are most critical for improving predictive performance. We propose a hierarchical semi-Markov framework that addresses two key modeling challenges. First, a time-inhomogeneous Markov router learns the patterns of “which activity comes next.” Second, a semi-Markov hazard component explicitly models activity durations, capturing “how long” activities realistically last. To ensure statistical stability when data are sparse, the model pools information across related demographic groups and time blocks. The entire framework is trained and evaluated using survey design weights to ensure our findings are representative of the U.S. population. On a held-out test set, we demonstrate that explicitly modeling durations with the hazard component provides a substantial and statistically significant improvement over purely Markovian models. Furthermore, our analysis reveals a clear hierarchy of demographic factors: Sex, Day-Type, and Household Size provide the largest predictive gains, while Region and Season, though important for energy calculations, contribute little to predicting the activity sequence itself. The result is an interpretable and robust generator of synthetic activity traces, providing a high-fidelity foundation for downstream energy systems modeling.

24 POWER TRANSMISSION AND DISTRIBUTION↗

EWPD in the SMEFT to dimension eight

We calculate the O H † H 2 / Λ 4 corrections to LEP electroweak precision data using the geometric formulation of the Standard Model Effective Field Theory (SMEFT). We report our results in simple-to-use interpolation tables that allow the interpretation of this data set to dimension eight for the first time. We demonstrate the impact of these previously unknown terms in the case of a general analysis in the SMEFT, and also in the cases of two distinct models matched to dimension eight. Neglecting such dimension-eight corrections to LEP observables introduces a theoretical error in SMEFT studies. We report some preliminary studies defining such a theory error, explicitly demonstrating the effect of previously unknown dimension-eight SMEFT corrections on LEP observables.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

How macroscopic limits on neutron-star baryon loss yield microscopic limits on non-standard-model baryon decay

We investigate how our baryon-loss limits from anomalous binary-pulsar period lengthening can be interpreted microscopically to yield specific constraints on the particle physics of baryon number violation within a neutron star. We focus on the possibility of anomalous baryon disappearance via dark baryon processes and on scenarios in which the produced dark-sector particles do not survive to influence the response of the star to baryon-number-violating effects. We flesh out the conditions for which this may occur, as well as other key assumptions. We then turn to the analysis of particle processes in the dense nuclear medium found at the core of a neutron star, employing the techniques of relativistic mean-field theory. Using our study of in-medium effects and limits on macroscopic baryon number violation, we extract limits on in-vacuum baryon-number-violating processes, and we determine them for various equations of state. We conclude by noting the implications of our results for models of dark-sector-enabled baryogenesis. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

SigTime: Learning and Visually Explaining Time Series Signatures

Understanding and distinguishing temporal patterns in time series data is essential for scientific discovery and decision-making. For example, in biomedical research, uncovering meaningful patterns in physiological signals can improve diagnosis, risk assessment, and patient outcomes. However, existing methods for time series pattern discovery face major challenges, including high computational complexity, limited interpretability, and difficulty in capturing meaningful temporal structures. Here, to address these gaps, we introduce a novel learning framework that jointly trains two Transformer models using complementary time series representations: shapelet-based representations to capture localized temporal structures and traditional feature engineering to encode statistical properties. The learned shapelets serve as interpretable signatures that differentiate time series across classification labels. Additionally, we develop a visual analytics system—SigTime—with coordinated views to facilitate exploration of time series signatures from multiple perspectives, aiding in useful insights generation. We quantitatively evaluate our learning framework on eight publicly available datasets and one proprietary clinical dataset. Additionally, we demonstrate the effectiveness of our system through two usage scenarios along with the domain experts: one involving public ECG data and the other focused on preterm labor analysis.

97 MATHEMATICS AND COMPUTING↗

A physics-informed and hierarchically regularized data-driven model for predicting fluid flow through porous media

This paper presents a new deep learning data-driven model for predicting structure dependent pore-fluid velocity fields in rock. The model is based on a Convolutional Auto-Encoder (CAE) artificial neural network capable of learning from image data generated by direct numerical simulations of fluid flow through pore-structures, such as by Lattice Boltzmann or molecular dynamics methods. The main novelty of the model in comparison to previous CAE-based data-driven approaches consists of three parts. The first is a methodology for decomposing the full-domain of the porous media into sub-regions, or “sub-domains”, in order to reduce the overall size of the CAE, batch process the sub-domains in parallel, and enable the CAE to learn local and generalizable nonlinear mappings of pore-fluid velocities. The second consists of embedding the finite difference solutions of the incompressible Navier-Stokes and continuity equations into convolutional layers prior to the CAE in order to provide the CAE with knowledge of fluid dynamics physics (PhyFlow). The third main novelty is that the training of the CAE is regularized with a hierarchical loss function that encourages the learning of fluid flow patterns (in a way similar to ranked modes in principal component analysis), ranking from most to least important. This is shown to increase the stability in learning, reduce over-fitting, and promote interpretability of the CAE neural network layers (HierCAE). The comprehensive new data-driven model, which we call the PhyFlow-HierCAE model, is shown to exhibit improved accuracy and generalizability of flow field predictions over conventional CAE models, attributable to the embedded physical knowledge and the hierarchical regularization, as well as realize orders of magnitude speed-ups in computation times as a surrogate for the direct numerical simulations. Examples of training and forward predictions on unseen pore-structures are provided and evaluated for data from Lattice Boltzmann and molecular dynamics simulations of pore-fluid flow. The model is shown to be a fast and accurate emulator (or “surrogate”) for predicting effective permeability of unseen pore-structures based on learning from relatively small direct numerical simulation datasets.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗