SEARCH · Engineering Papers
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.
Presentation for my Nuclear Engineering II course at Idaho State University about critical data extrapolation of ISU's AGN-201m reactor.
Nuclear nonproliferation serves as a key goal, being undertaken by the International Atomic Energy Agency (IAEA). To recognize proliferation there are two pathways that states, who intend to use nuclear material for malicious purposes can take, diversion can misuse. Diversion is when fissile nuclear material is declared to the IAEA for non-weapon purposes, but then covertly removed. If the source of nuclear material, that is not declared and not fissionable, is placed inside the reactor core to create fissile material used to create weapons then the state is using the second pathway of proliferation, misuse. With the emerging development in areas of simulation and machine learning the creation of virtual models of reactor systems, digital twins, serve as a potential method to identify proliferation through detecting anomalous behavior in the reactor. A digital twin for a physical nuclear reactor has never been developed, as digital twins serve as an emerging technology. To investigate the process for development and use of a digital twin for a nuclear reactor Idaho State University’s AGN-201m serves as the nuclear reactor used for development of this digital twin. A data acquisition system has been installed to the reactor system allowing for the transfer of collected data from a reactor operation to Idaho National Laboratory’s Deeplynx data warehouse. When utilizing data to train reactor physics and machine learning models, a significant challenge encountered is the initial state of the data. Nuclear proliferation will have the capacity to be detected when the reactor immediately starts up, nor will it occur after the reactor shuts down. Generally, it will be detected when the reactor is operating at some desired power over a sufficient period for that specific reactor design. For the AGN-201m this will be when the reactor is critical (generally 1 mW or above) for a timespan that is within or less than the range of a regular business day. Datasets sent to Deeplynx have had to be manually cut to when the reactor is critical based on plots of power levels. This method is inefficient and laborious, especially when using multiple datasets at once to train a model. To provide a more streamlined approach an automated critical data extrapolator is developed, with capabilities of recognizing when the reactor operation first reaches criticality, and when the reactor undergoes a SCRAM and is shutdown.
This report provides a method for evaluating coordinate measuring machine (CMM) data and extrapolating the data to estimate the dimensions of a cylindrical component beyond the sampled measurements. The theoretical basis is provided, and a calculational method is presented. The report also provides detailed, formatted, printed, and plotted representations.
Measuring the spin structure of protons and neutrons tests our understanding of how they arise from quarks and gluons, the fundamental building blocks of nuclear matter. At long distances, the coupling constant of the strong interaction becomes large, requiring non-perturbative methods to calculate quantum chromodynamics processes, such as lattice gauge theory or effective field theories. In this work, we report proton spin structure measurements from scattering a polarized electron beam off polarized protons. The spin-dependent cross-sections were measured at large distances, corresponding to the region of low momentum transfer squared between 0.012 and 1.0GeV 2 . This kinematic range provides unique tests of chiral effective field theory predictions. Our results show that a complete description of the nucleon spin remains elusive, and call for further theoretical works, for example, in lattice quantum chromodynamics. Finally, our data extrapolated to the photon point agree with the Gerasimov–Drell–Hearn sum rule, a fundamental prediction of quantum field theory that relates the anomalous magnetic moment of the proton to its integrated spin-dependent cross-sections.
The data used to produce all of the figures and tables in the manuscript titled "Highly Accurate Many-Body Theory Reaches 2D Materials" can be found here. This data set includes: -SHDMC results for graphene -selected CI with and without re-normalized second-order perturbation (rPT2) theory corrections for graphene -data demonstrating that SHDMC displays an exponential rate of convergence -data used for sCI + rPT2 complete basis set extrapolation -data used to extrapolate SHDMC energies to the infinite basis limit -data used to demonstrate compactness of SHDMC wavefunction
Data to reproduce the paper "Data efficiency and extrapolation trends in neural network interatomic potentials" (LLNL-JRNL-845001). Contains: plots, results, and visualization scripts. Preprint available at https://doi.org/10.48550/arXiv.2302.05823.
The purpose of this study was to develop a process to convert input signals from one facility into another by reflecting geometric and environmental settings. The Dynamic Energy Transport and Integration Laboratory (DETAIL) is a research facility in development. Its aim is to emulate the daily interactions among power production industry systems and receive real-time data from those systems as inputs. To convert signals and ensure that the temporal sequences and magnitudes reflect laboratory settings, the ability to scale and project data is essential. To demonstrate this ability, Dynamical System Scaling (DSS) and Hierarchical Two-Tiered Scaling (H2TS) (methodologies that enable systems to scale and project or extrapolate data sets to desired environments while conserving the observed behavior based on first principles) were applied to DETAIL’s thermocline thermal storage system in the Thermal Energy Distribution System (TEDS) facility and solid-oxide electrolysis cell in the High Temperature Hydrogen Electrolysis (HTHE) facility. Both thermocline and electrolysis cell systems were successfully scaled, and test cases were conducted to generate a doubly accelerated energy charge and discharge in reference to past experimental data from the facilities. The research results represented a case for the thermocline system that required signals to be accelerated without altering the stored energy. To enhance the quality of the accelerated data, error propagation analyses were conducted on DSS post-processing terms to determine the consequences of raw-data-associated errors.
In this paper, we analyze the numerical aspects of the inherent multi reference density matrix renormalization group (DMRG) calculations on top of the periodic Kohn-Sham density functional theory using the complete active space approach. The potential of the framework is illustrated by studying hexagonal boron nitride nanoflakes embedding a charged single boron vacancy point defect by revealing a vertical energy spectrum with a prominent multireference character. We investigate the consistency of the DMRG energy spectrum from the perspective of sample size, basis size, and active space selection protocol. Results obtained from standard quantum chemical atom-centered basis calculations and plane-wave based counterparts show excellent agreement. Furthermore, we also discuss the spectrum of the periodic sheet which is in good agreement with extrapolated data of finite clusters. These results pave the way toward applying the DMRG method in extended correlated solid-state systems, such as point defect qubit in wide band gap semiconductors.
Abstract The immiscible alloy Cu–Ta has the potential for enhanced mechanical performance in applications as a functional coating. To establish baseline mechanical properties, four Cu–Ta films were co-sputtered at the temperatures 23, 400, 600, and 800 °C and tested with nanoindentation at strain rates 5 $$\times $$ × 10 −3 s −1 to 10 s −1 . Each film had a unique microstructure morphology. The hardness and elastic modulus of the four films were insensitive to strain rate changes. Instead, the measured properties were spatially dependent, particularly in the 600 and 800 °C films. In those two films, there is a bimodal deformation behavior due to Cu-agglomeration under protruding grains and planar Ta-rich regions. Increasing the indentation depth revealed shear band suppression which is related to a homogenous distribution of flow stresses for all four microstructure morphologies. Finally, the Cu–Ta hardness appeared to follow a rule-of-mixtures when compared to extrapolated data of Cu and Ta monolithic films.
Abstract Human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) have been widely used in the Comprehensive in vitro Proarrhythmia Assay (CiPA). The notable difference of the electrophysiological (EP) responses of hiPSC-CMs in serum and serum-free media (SFM) is puzzling and may impact regulatory decision-making on the cardiac safety of candidate drugs in inducing QT prolongation and torsade de pointes (TdP). In this study, we compared the EP responses of hiPSC-CMs to 10 CiPA compounds and moxifloxacin in serum and SFM; explained the potential reason behind the different EP responses—abiotic compound loss to plastic tubes/plates of hydrophobic compounds prepared in SFM; and investigated the impact of compound preparation methods on drug bioavailability in exposure media, which affects the TdP risk prediction of drugs tested in serum-containing and SFM. For assays to be conducted in SFM, awareness of abiotic compound loss of hydrophobic compounds in serum-free preparations is critical for delay repolarization evaluation and data extrapolation from in vitro to in vivo.
The $^{12}$C($\alpha $,$\gamma $)$^{16}$O reaction rate is crucial in determining the carbon-to-oxygen abundance ratio in stellar nucleosynthesis. Measuring this reaction’s cross section at stellar energies is challenging due to its extremely small value, approximately 10 -17 barn at E c.m. = 300 keV. To address this, R-matrix calculations are employed to extrapolate data to lower energies, requiring a comprehensive understanding of each contribution to the cross section. The dominant contributions to the cross section at stellar energies arise from electric dipole (E1) and electric quadrupole (E2) transitions to the ground state of 16 O, along with a significant cascade contribution. Traditionally, these contributions have been separated using the γ-ray angular distribution. In this work, we propose a novel technique using the energy distribution of the 16 O recoils at the focal plane. This method involves a neural network trained on detailed Monte Carlo simulations of the energy distribution of recoils transported through the recoil mass separator ERNA. This approach enables the simultaneous determination of all three contributions with errors around 10% in the energy range E c.m. = 1.0–2.2 MeV. In conclusion, by employing this new technique, we aim to significantly improve the accuracy of determining the cross section of the $^{12}$C($\alpha $,$\gamma $)$^{16}$O reaction at astrophysical energies.
Battery lifetime models are used to extrapolate data from accelerated aging tests to simulate degradation in real-world applications such as electric vehicles and battery energy storage systems. Methods developed at NREL utilize both expert domain-knowledge and machine-learning to identify models, using statistical methods such as cross-validation and bootstrap resampling to interrogate model performance and quantify uncertainty. These models can be utilized in systems level simulations to predict battery performance or technoeconomic models to estimate the lifetime cost of battery systems.
The internal structure of protons and neutrons provides insight into both the dynamical behavior of the constitute quarks and gluons, and emergent properties of the nucleons (such as mass, spin, and electromagnetic distributions). Elastic electron-nucleon scattering can probe the elastic electromagnetic form factors of the nucleon. The electric and magnetic form factors, respectively, encode information about the internal charge and magnetization distributions within the nucleon. Precision data for these form factors, over a broad range of the four-momentum transfer squared, Q^2, can benchmark theoretical models describing the strong interaction of nuclear physics. The Super BigBite Spectrometer (SBS) program in Hall A at Jefferson Lab, is a series of high-precision experiments which seek to significantly extend the Q^2 reach of previous data for the nucleon electromagnetic form factors. The first two experiments of this program are known as G_M^n and the neutron Two Photon Exchange (nTPE) and the data were collected from October 2021 to February 2022. Both experiments were conducted with the simultaneous measurement of D(e,e'n) and D(e,e'p) reactions for quasi-elastic electron-deuteron scattering. The scattered electrons were detected in the BigBite Spectrometer, which features multiple large-acceptance Gas Electron Multiplier (GEM) detectors. The Super BigBite Spectrometer provided simultaneous detection of scattered nucleons, and utilized a large acceptance dipole magnet and Hadron Calorimeter (HCal). The G_M^n experiment provides precision measurements of the neutron magnetic form factor, via the ratio method, over a Q^2 range of 3.0 to 13.5 (GeV/c)2. From this data analysis, preliminary values for G_M^n/µ_n G_D are extracted. For Q^2=4.48 (GeV/c)2 we find G_M^n/µ_n G_D=0.9546±0.0132 and for Q^2=4.476 (GeV/c)2 we find G_M^n/µ_n G_D=0.9563±0.0110. These preliminary G_M^n/µ_n G_D values are more precise than existing world data in this Q^2 regime and are consistent with the most recent parameterization of the G_M^n/µ_n G_D world data. The nTPE experiment provides a first measurement of the neutron Rosenbluth Slope and seeks to quantify the two-photon exchange(TPE) contribution to elastic electron-neutron scattering at a fixed Q^2=4.5 (GeV/c)2 with two different beam energies and scattering angle values. For data of the proton form factor ratio, µ_p G_E^p/G_M^p, significant discrepancies exist between values obtained from Rosenbluth Separation and polarization transfer measurement, particularly at large Q^2, and TPE contributions are thought to resolve this discrepancy. The impacts of TPE contributions have not yet been experimentally established for the neutron. From the data analysis presented in this dissertation, a preliminary result for the neutron Rosenbluth Slope is found as S^n=(G_E^n )^2/t_n (G_M^n )^2=0.0916±0.0476 for Q^2=4.48 (GeV/c)2. This value of the neutron Rosenbluth Slope is consistent with the world data extrapolation and the absence of large TPE corrections.
Class B Materials Code Case to include: FY25 - High temperature design methodology, overview of ASME Section III (Division 5), Design for elevated temperature service, Goal to address current gaps, evaluation of proposed rules, allowable stress development, Class A materials (existing ASME approach), variable confidence index, data extrapolation, bounding value for SEE, and next steps.
Abstract Recently, key architectural advances have been proposed for neural network interatomic potentials (NNIPs), such as incorporating message-passing networks, equivariance, or many-body expansion terms. Although modern NNIP models exhibit small differences in test accuracy, this metric is still considered the main target when developing new NNIP architectures. In this work, we show how architectural and optimization choices influence the generalization of NNIPs, revealing trends in molecular dynamics (MD) stability, data efficiency, and loss landscapes. Using the 3BPA dataset, we uncover trends in NNIP errors and robustness to noise, showing these metrics are insufficient to predict MD stability in the high-accuracy regime. With a large-scale study on NequIP, MACE, and their optimizers, we show that our metric of loss entropy predicts out-of-distribution error and data efficiency despite being computed only on the training set. This work provides a deep learning justification for probing extrapolation and can inform the development of next-generation NNIPs.
The connector/adapter assembly (CAA) is one component of a higher-level assembly. A simplified representation of the CAA is illustrated in Figure 1. The location of each CAA is required to meet a specific positional requirement. Several CAAs have been measured and found to be out of tolerance. The connector and the adapter are welded to each other such that the cylindrical axes of the two components are coaxial. The positional tolerance condition requires that all points on the outer surface of the CAA must be within a defined maximum diameter, as illustrated in Figure 1. Two different methods have been developed to measure the CAA positions from coordinate measuring machine (CMM) data. One method (RevD) extrapolates the data from two circles approximately midway from the weld region to each end. The second method (RevF) measures points that are much closer to the actual ends.
The mean excitation energy (I-value) of liquid argon is a critical input for energy estimation in neutrino oscillation experiments. It is measured to be 199 +4 -5 eV using the range of 402.2 MeV protons from the Fermilab Linac. This compares to the author's recent evaluation of (197 ± 7) eV based on a combination of an oscillator strength distribution analysis, gaseous argon range measurements, sparse stopping power data on solid argon, and an extrapolation of data on the effect of phase from other substances. Using all sources of information, we recommend a value of (198 ± 4) eV for liquid argon, which is significantly higher than 188 eV, from ICRU-37's gaseous argon evaluation, commonly used in Monte Carlo codes such as Geant4.
Representing real-time data as a sum of complex exponentials provides a compact form that enables both denoising and extrapolation. As a fully data-driven method, the Estimation of Signal Parameters via Rotational Invariance Techniques (ESPRIT) algorithm is agnostic to the underlying physical equations, making it broadly applicable to various observables and experimental or numerical setups. In this work, we consider applications of the ESPRIT algorithm primarily to extend real-time dynamical data from simulations of quantum systems. We evaluate ESPRIT's performance in the presence of noise and compare it to other extrapolation methods. We demonstrate its ability to extract information from short-time dynamics to reliably predict long-time behavior and determine the minimum time interval required for accurate results. We discuss how this insight can be leveraged in numerical methods that propagate quantum systems in time, and we show how ESPRIT can predict infinite-time values of dynamical observables, offering a purely data-driven approach to characterizing quantum phases.