Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “data dependencies”

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,369 records · Page 76

Interpreting and reporting 40 Ar/ 39 Ar geochronologic data

The 40Ar/39Ar dating method is among the most versatile of geochronometers, having the potential to date a broad variety of K-bearing materials spanning from the time of Earth’s formation into the historical realm. Measurements using modern noble-gas mass spectrometers are now producing 40Ar/39Ar dates with analytical uncertainties of ~1%, thereby providing precise time constraints for a wide range of geologic and extraterrestrial processes. Analyses of increasingly smaller subsamples have revealed age dispersion in many materials, including some minerals used as neutron fluence monitors. Accordingly, interpretive strategies are evolving to address observed dispersion in dates from a single sample. Moreover, inferring a geologically meaningful “age” from a measured “date” or set of dates is dependent on the geological problem being addressed and the salient assumptions associated with each set of data. We highlight requirements for collateral information that will better constrain the interpretation of 40Ar/39Ar data sets, including those associated with single-crystal fusion analyses, incremental heating experiments, and in situ analyses of microsampled domains. To ensure the utility and viability of published results, we emphasize previous recommendations for reporting 40Ar/39Ar data and the related essential metadata, with the amendment that data conform to evolving standards of being findable, accessible, interoperable, and reusable (FAIR) by both humans and computers. Our examples provide guidance for the presentation and interpretation of 40 Ar/ 39 Ar dates to maximize their interdisciplinary usage, reproducibility, and longevity.

58 GEOSCIENCES↗

Comparison of Conventional and Constrained Variational Methods for Computing Large-Scale Budgets and Forcing Fields

Analyses of atmospheric heat and moisture budgets serve as an effective tool to study convective characteristics over a region and to provide large-scale forcing fields for various modeling applications. This paper examines two popular methods for computing large-scale atmospheric budgets: the conventional budget method (CBM) using objectively gridded analyses based primarily on radiosonde data and the constrained variational analysis (CVA) approach which supplements vertical profiles of atmospheric fields with measurements at the top of the atmosphere and at the surface to conserve mass, water, energy, and momentum. Successful budget computations are dependent on accurate sampling and analyses of the thermodynamic state of the atmosphere and the divergence field associated with convection and the large-scale circulation that influences it. Utilizing analyses generated from data taken during Dynamics of the Madden-Julian Oscillation (DYNAMO) field campaign conducted over the central Indian Ocean from October to December 2011, we evaluate the merits of these budget approaches and examine their limitations. While many of the shortcomings of the CBM, in particular effects of sampling errors in sounding data, are effectively minimized with CVA, accurate large-scale diagnostics in CVA are dependent on reliable background fields and rainfall constraints. Finally, for the DYNAMO analyses examined, the operational model fields used as the CVA background state provided wind fields that accurately resolved the vertical structure of convection in the vicinity of Gan Island. However, biases in the model thermodynamic fields were somewhat amplified in CVA resulting in a convective environment much weaker than observed.

54 ENVIRONMENTAL SCIENCES↗

Matrix Multiply Performance of GPUs on Exascale-class HPE/Cray Systems

The computation of dense matrix-matrix products (GEMMs) is central to many modeling and simulation workloads as well as AI/ML deep learning campaigns. In fact, millions of dollars are spent annually on computing GEMMs, and large model training demands are increasing exponentially. Specialized processors such as GPUs are designed to perform well for these operations. However, the performance of GEMMs on GPUs can exhibit complex behaviors depending on many factors, making it challenging to optimize the performance of GEMMs on these processors. In this study we undertake an examination of GEMM performance on several leading GPU models taken from product lines of GPUs to be deployed in forthcoming exascale computing systems. We show results to illustrate the many factors that can affect performance of GEMMs on GPUs. We then present data collected from a large number of test runs for an example GEMM operation to show the dependence behaviors of GEMM rate on matrix dimensions. Finally, we show results from machine learning-based performance models using novel feature engineering methods to fit the measured performance, providing a potential basis for GEMM performance tuning and autotuning methods for GPUs. Recommendations are also given for how to achieve high GEMM performance on modern GPUs.

Melesse Vergara, Veronica↗

Structural Evidence of Interanionic Hydrogen Bonding in Phosphoric Acid Solutions

Interanionic hydrogen bonding (IAHB) is a noncovalent interaction between like-charged ions that challenges conventional electrostatic understanding. This study provides direct structural evidence of IAHB in concentrated aqueous phosphoric acid (PA) solutions, which exhibit >60% dissociation under these conditions. Oxygen K-edge X-ray absorption fine structure spectroscopy, combined with electron affinity time-dependent density functional theory calculations, reveals the formation of stable, cyclic phosphate-phosphate IAHB dimers at PA concentrations ≥7 M. Extended X-ray absorption fine structure data show distinct long-range ordering consistent with these dimers, and near-edge X-ray absorption fine structure spectra confirm a concentration-dependent transition from monomeric to dimeric species. Energy decomposition analysis through density functional theory shows that the formation of solution-phase IAHB is energetically favored and is attributed to polarization of and the charge transfer between the two fragments driven by the surrounding solvent molecules, in addition to the permanent electrostatics. These findings offers crucial structural insights into the H-bonded networks in concentrated PA, highlighting the critical role of solvent in facilitating anion–anion association.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Accelerating multiscale electronic stopping power predictions with time-dependent density functional theory and machine learning

Knowing the rate at which particle radiation releases energy in a material, the “stopping power,” is key to designing nuclear reactors, medical treatments, semiconductor and quantum materials, and many other technologies. While the nuclear contribution to stopping power, i.e., elastic scattering between atoms, is well understood in the literature, the route for gathering data on the electronic contribution has for decades remained costly and reliant on many simplifying assumptions, including that materials are isotropic. We establish a method that combines time-dependent density functional theory (TDDFT) and machine learning to reduce the time to assess new materials to hours on a supercomputer and provide valuable data on how atomic details influence electronic stopping. Our approach uses TDDFT to compute the electronic stopping from first principles in several directions and then machine learning to interpolate to other directions at a cost of 10 million times fewer core-hours. We demonstrate the combined approach in a study of proton irradiation in aluminum and employ it to predict how the depth of maximum energy deposition, the “Bragg Peak,” varies depending on the incident angle—a quantity otherwise inaccessible to modelers and far outside the scales of quantum mechanical simulations. The lack of any experimental information requirement makes our method applicable to most materials, and its speed makes it a prime candidate for enabling quantum-to-continuum models of radiation damage. The prospect of reusing valuable TDDFT data for training the model makes our approach appealing for applications in the age of materials data science.

36 MATERIALS SCIENCE↗

Exploring the Three-Dimensional Momentum Distribution of Longitudinally Polarized Quarks in the Proton

By analyzing experimental data on semi-inclusive deep inelastic scattering off longitudinally polarized targets, we extract the transverse momentum dependence of the quark helicity distribution, i.e., the difference between the three-dimensional motion of quarks with polarization parallel or antiparallel to the longitudinal polarization of the parent hadron. We perform the analysis at next-to-leading and next-to-next-to-leading perturbative accuracy. The quality of the fit is very good for both cases, reaching a χ 2 per number of data points equal to 1.11 and 1.09, respectively. Although the limited number of data points leads to significant uncertainties, the data are consistent with an interpretation in which the helicity distribution is narrower in transverse momentum than the unpolarized distribution. Published by the American Physical Society 2025

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Constraints on the IBEX Ribbon’s Origin from Its Evolution over a Solar Cycle

Abstract In 2009, the Interstellar Boundary Explorer (IBEX) discovered a narrow “ribbon” of energetic neutral atom emissions across the sky with properties correlated with the solar wind latitudinal structure and the interstellar magnetic field draped around the heliosphere. It is widely believed that the ribbon is formed from the escape of heliospheric ENAs into the local interstellar medium and their eventual return as secondary ENAs. However, there is no consensus on the rate of pitch angle scattering of these PUIs before they become secondary ENAs. We test two opposing limits of scattering rates (“weak” versus “strong”) by solving a time-dependent model of the ribbon that evolves with the solar cycle, and we compare them to IBEX observations over 2009–2019. First, we find that both models qualitatively reproduce the evolution of IBEX fluxes for most of the data set, with a few exceptions, although the strong (or “spatial retention”) scattering model greatly underestimates the observed fluxes. Regardless, time dependence of fluxes cannot distinguish these models. Second, the ribbon’s geometric properties, i.e., its center and radius, are significantly different between the models. The spatial retention model reproduces the observed ribbon centers as a function of energy and time slightly better than the weak scattering model, and the spatial retention model reproduces the observed ribbon radius over energy and time almost perfectly, whereas the weak scattering model compares poorly. Our analysis favors the spatial retention mechanism as the source of the IBEX ribbon, but it requires modification to increase the flux of ENAs observed at 1 au.

Zirnstein, E. J. (ORCID:0000000172400618)↗

Determining reference standard strength for neutron-irradiated reduced activation ferritic/martensitic steel F82H by Bayesian method

The deterministic approach widely adopted in the design of structural components relies on systematically defined design limits using empirically determined safety factors. However, this approach is not always appropriate because structures are subjected to a variety of loads in the practical environment, which may result in excessively conservative design limits. In recent years, a more rigorous probabilistic approach that incorporates material strength distributions has become an important solution. In the probabilistic approach, the probability density functions of material strength properties underpin the design criteria. Here, the objective of this study is to identify the density distribution functions that best describe tensile properties of irradiated F82H to define a reference strength for DEMO design. Due to the limited number of existing data, this study specifically employs a Bayesian prediction method based on Monte Carlo simulations to determine a material reference value with statistical reliability and to investigate its effectiveness. For example, the dependence of tensile properties of 300 °C irradiated materials on irradiation damage and the range predicted by 95% Bayesian estimation was evaluated. As a statistical model for the dose dependence of statistical parameters, the normal distribution exhibited a better fit for 0.2% proof strength and tensile strength, whereas the distribution of total elongation data gave comparable reference values for both the normal and Weibull distribution models. Both models gave comparable criteria for the distribution of total elongation data. The Weibull model also gave better results for uniform elongation. The function best describing the model was a logarithmic law for both 0.2% proof strength and tensile strength, while a power law for both total and uniform elongation, which allowed for more comprehensive data prediction of irradiation data with statistical accuracy for DEMO reactor design.

36 MATERIALS SCIENCE↗

Energy dependent calculations of fission product, prompt, and delayed neutron yields for neutron induced fission on 235 U, 238 U, and 239 Pu

In this article, we perform energy-dependent calculations of independent and cumulative fission product yields for 235 U, 238 U, and 239 Pu in the first chance fission region. Starting with the primary fission fragment distributions taken from available experimental data and analytical functions based on assumptions for the excitation energy and spin-parity distributions, the Hauser-Feshbach statistical decay treatment for fission fragment de-excitation is applied to more than 1,000 fission fragments for the incident neutron energies up to 5 MeV. The calculated independent yields are then used as an input of β-decay calculations to produce the cumulative yield, and summation calculations are performed. Model parameters in these procedures are adjusted by applying the Bayesian technique at the thermal energy for 235 U and 239 Pu and in the fast energy range for 238 U. The calculated fission observable quantities, such as the energy-dependent cumulative yields, and prompt and delayed neutron yields, are compared with available experimental data. We also study the possible impact of the second chance fission opening on the energy dependence of the delayed neutron yield by on the opening of second chance fission on the energy dependence calculation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Deep glassy state dynamic data challenge glass models: Configurational entropy models

This work was the result of the student's (Dongjie Chen) qualifying exam response. In it the experimental observation from our labs that the relaxation times of ultrastable glasses (a 20 million year old amber and a vapor deposited amorphous teflon) do not follow the expected Vogel-Fulcher or WLF types of divergence as the temperatures fall below the glass transition temperature. In the work we evaluated several theories by considering them in different ways. We used them as generally used and show that they do not capture the non-diverging behavior. Because the models were all entropy-based models, we also considered the possibility that the temperature dependence of the entropy used in the models is incorrect and chose two ways of estimating the entropy that do not follow the classical theoretical models. First we used the actually measured calorimetric response for the two ultrastable glasses, and we used the Milchev model, which does not go to zero above 0 Kelvin. Both of these approaches led to improved agreement between the theories and the experimental data. Finally, we also asked what the temperature dependence of the entropy would have to be so that the models agreed with the observed data. In some cases the entropy so-calculated was reasonable, while in others it was unphysical. In sum, the important point of the work is that the observation of non-diverging time-scales that suggest that there is no ideal glass transition can be reasonably captured by current theories if they take more realistic models of the entropy than is usually the case.

36 MATERIALS SCIENCE↗

Constraining the Charge-Sign and Rigidity-Dependence of Solar Modulation

Our ability to identify the sources of cosmic rays and understand how these particles propagate through the interstellar medium is hindered by the combined effects of the solar wind and its embedded magnetic field, collectively known as solar modulation. In this paper, we build upon our previous work to model and constrain the effects of solar modulation on the cosmic-ray spectrum, using data from AMS-02 and BESS Polar II collected between 2007 and 2012, during which the heliospheric magnetic field was in a state of negative polarity. Our model uses measurements of the heliospheric magnetic field and the tilt angle of the heliospheric current sheet to accurately predict the effects of solar modulation as a function of time, charge, and rigidity. By incorporating data from a period of negative polarity, we have been able to robustly observe and constrain the charge-dependent effects of solar modulation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Stable Methane Isotopologues From Northern Lakes Suggest That Ebullition Is Dominated by Sub‐Lake Scale Processes

Abstract Stable isotopes have emerged as popular study targets when investigating emission of methane (CH 4 ) from lakes. Yet little is known on how isotopic patterns conform to variations in emission magnitudes—a highly relevant question. Here, we present a large multiyear data set on stable isotopes of CH 4 ebullition (bubbling) from three small adjacent subarctic lakes. The δ 13 C‐CH 4 and δD‐CH 4 range from −78.4‰ to −53.1‰ and from −369.8‰ to −218.8‰, respectively, and vary greatly among the lakes. The signatures suggest dominant hydrogenotrophic methanogenesis, particularly in the deep zones, but there are also signals of seemingly acetoclastic production in some high fluxing shallow areas, possibly fueled by in situ vegetation, but in‐sediment anaerobic CH 4 oxidation cannot be ruled out as an alternative cause. The observed patterns, however, are not consistent across the lakes. Neither do they correspond to the spatiotemporal variations in the measured bubble CH 4 fluxes. Patterns of acetoclastic and hydrogenotrophic production plus oxidation demonstrate that gains and losses of sediment CH 4 are dominated by sub‐lake scale processes. The δD‐CH 4 in the bubbles was significantly different depending on measurement month, likely due to evaporation effects. On a larger scale, our isotopic data, combined with those from other lakes, show a significant difference in bubble δD‐CH 4 between postglacial and thermokarst lakes, an important result for emission inventories. Although this characteristic theoretically assists in source partitioning studies, most hypothetical future shifts in δD‐CH 4 due to high‐latitude lake area or production pathway are too small to lead to atmospheric changes detectable with current technology.

54 ENVIRONMENTAL SCIENCES↗

A data-driven multiscale model for reactive wetting simulations

Here, we describe a data-driven, multiscale technique to model reactive wetting of a silver–aluminum alloy on a Kovar™ (Fe-Ni-Co alloy) surface. We employ molecular dynamics simulations to elucidate the dependence of surface tension and wetting angle on the drop’s composition and temperature. A design of computational experiments is used to efficiently generate training data of surface tension and wetting angle from a limited number of molecular dynamics simulations. The simulation results are used to parameterize models of the material’s wetting properties and compute the uncertainty in the models due to limited data. The data-driven models are incorporated into an engineering-scale (continuum) model of a silver–aluminum sessile drop on a Kovar™ substrate. Model predictions of the wetting angle are compared with experiments of pure silver spreading on Kovar™ to quantify the model-form errors introduced by the limited training data versus the simplifications inherent in the molecular dynamics simulations. The paper presents innovations in the determination of “convergence” of noisy MD simulations before they are used to extract the wetting angle and surface tension, and the construction of their models which approximate physio-chemical processes that are left unresolved by the engineering-scale model. Together, these constitute a multiscale approach that integrates molecular-scale information into continuum scale models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Tools to perform data reduction of time-of-flight data taken on BL4A (MagRef)—a Jupyter notebook approach

A collection of modules or tools to produce reflectivity profiles from data acquired with unpolarized and polarized neutron beams using the BL4A (MagRef) instrument at the Spallation Neutron Source. The Jupyter notebook approach offers transparency to the user on what is done with the data, and once the event mode data are stored as NumPy binary files empowers the user to perform the steps from their own laptop without the need for Mantid. The notebook offers two options for background subtraction, applies corrections for the wavelength dependence of the polarizer and flipper efficiencies, accounts for Zeeman energy splitting of the spin-flip channels, Qz-bins the data along contours of constant Qz, propagates errors and writes results in text format suitable for analysis with tools such as GenX.

Fitzsimmons, Michael [Oak Ridge Associated Univ., ↗

Experimental Validation of Thermal Scattering Evaluations

In order to test the performance of new neutron thermal scattering law (TSL) evaluations it is desirable to have experimental data that is highly sensitive to the TSL and provides high fidelity information on the energy dependent performance of TSL evaluations. Three relevant experiments are discussed including: accurate thermal total cross section measurements, thermal neutron die-away experiments, and neutron leakage experiments. The experimental setups and results are reviewed, and examples provided for some moderators including polyethylene, Plexiglas, and YHx. For the experiments preformed thus far, there is generally good agreement between the measured total cross section and simulations using current TSL evaluations, however in certain energy ranges differences were observed. Similarly, neutron die-away and leakage measurements for samples at room temperature are in good agreement with data computed from TSLs, however leakage measurements for polyethylene at 29K show discrepancies with TSL evaluations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Sap Velocity Data for Urban Trees in Chicago, Illinois (2024-2025)

This dataset contains uncorrected sap velocity measurements using the heat ratio method (HRM) collected using ICT International SFM1x sensors at five urban sites in Chicago, Illinois, as part of the DOE CROCUS project. The data includes continuous monitoring of sap velocity from various tree species, including Maples (Acer spp.): Sugar Maple (Acer saccharum), Silver Maple (Acer saccharinum), and Red Maple (Acer rubrum); Oaks (Quercus spp.): Swamp White Oak (Quercus bicolor); American Elm (Ulmus americana); Honey Locust (Gleditsia triacanthos); Cottonwood (Populus deltoides); and Tree of Heaven (Ailanthus altissima) across Chicago State University (CSU), Northeastern Illinois University (NEIU), Northwestern University (NU), University of Illinois Chicago (UIC), and West Woodlawn "Blacks in Green" (BIG). These include both street trees and those in urban park locations. Measurements were collected at 15-20 minute intervals, depending on the sensor, and transmitted via Long Range Wide Area Network (LoRaWAN) protocols. The wireless data was collected by Sage Network (https://sagecontinuum.org/) nodes. The dataset includes sensor ID, Global Positioning System (GPS) coordinates, tree species (common and scientific names), tree identification number, diameter at breast height (DBH in cm), uncorrected sap velocity measurements (cm/hr) from both inner and outer probes, and Sage Node identifiers so the data can be mapped to related variables such as air quality and wind speed that were collected on the Sage nodes. All timestamps are in local Chicago time (CDT/CST). Quality control flags are provided using a 3-bit binary system indicating physical range violations (< -10 or > 60 cm/hr), step spikes (absolute difference > 36 cm/hr), and stuck sensor conditions (> 10 consecutive identical values). These are raw data, not corrected for wood anatomy or species-specific characteristics. Data is provided in comma separated (CSV) format. This dataset is part of a larger collection of CROCUS environmental monitoring data, including linked datasets from Air Quality Transmitter (AQT) sensors, Weather Transmitter (WXT) sensors, and Multi-Function Research LoRaWAN (MFR) Nodes. DOIs for the supporting data are provided as part of this data package.

Chicago↗

High-pressure ignition delay time measurements of a four-component gasoline surrogate and its high-level blends with ethanol and methyl acetate

Ignition delay times of a gasoline surrogate (iso-octane/toluene/n-heptane/1-hexene at 55/25/15/5% by liquid volume, developed in McCormick et al. 2017) with high-level bioblendstock-gasoline surrogate blends (50% and 85% biofuel by liquid volume of each ethanol and methyl acetate) were collected behind reflected shock waves. Post-reflected-shock temperatures ranged from 968 to 1361 K at pressures of about 4, 10, 25, and 50 atm. Data were collected for real fuel–air mixtures at fuel lean and stoichiometric conditions (φ = 0.5, 1.0) with focus on the more-practical, fuel-lean conditions. Ignition delay times were measured from OH* chemiluminescence around 307 nm using an endwall diagnostic. The RON and MON of the surrogate were 90.3 and 84.7, respectively. Ethanol was chosen due to its wide use in flex-fuel vehicles and availability, while also providing an increased octane rating. Methyl acetate was chosen for its especially high octane rating to investigate the effect of an extreme case. To validate the gasoline surrogate, the data are compared to real gasoline (RD387) and several gasoline surrogate experiments from the literature. Using the wide range of pressures studied, the gasoline surrogate’s pressure dependence was quantified to account for test-to-test variations using regression analysis. Similarly, a global correlation for gasoline and its surrogates was developed using all available data from the literature. Two modern chemical kinetics models targeting gasoline and its surrogates are compared to the ignition delay time measurements. These new high-pressure, high-bioblendstock concentration tests provide required chemical kinetic data for optimizing fuel and engine design.

42 ENGINEERING↗

Examining the EMC effect using the ${F}_{2}^{n}$ neutron structure function

The persistently mysterious deviations from unity of the ratio of nuclear target structure functions to those of deuterium as measured in deep inelastic scattering (often termed the “EMC effect”) have become the canonical observable for studies of nuclear medium modifications to free nucleon structure in the valence regime. The structure function of the free proton is well known from numerous experiments spanning decades. The free neutron structure function, however, has remained difficult to access. Recently it has been extracted in a systematic study of the global data within a parton distribution function extraction framework and is available from the CTEQ–Jefferson Lab (CJ) Collaboration. Here, we leverage the latter to introduce a new method to study the EMC effect in nuclei by reexamining existing data in light of the the magnitude of the medium modifications to the free neutron and proton structure functions independently. From the extraction of the free neutron from world data, it is possible to examine the nuclear effects in deuterium and their contribution to our interpretation of the EMC effect. In this study, we observe that the ratio of the deuteron to the sum of the free neutron and proton structure functions has some x B and Q 2 dependencies that impact the magnitude of the EMC effect as typically observed. Specifically, different EMC slopes are obtained when data from different x B and Q 2 values are utilized. Furthermore, while a linear correlation persists between the EMC and short range correlation effects, the slope is modified when deuteron nuclear effects are removed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗