Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Variational inference”

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

SDSS-IV MaNGA: Bayesian analysis of the star formation history of low-mass galaxies in the local Universe

ABSTRACT We measure the star formation histories (SFHs) of a sample of low-mass galaxies with M* < 109 M⊙ from the SDSS-IV MaNGA survey. The large number of IFU spectra for each galaxy are either combined to reach a high signal to noise ratio or used to investigate spatial variations. We use Bayesian inferences based on full spectrum fitting. Our analysis based on Bayesian evidence ratio indicates a strong preference for a model that allows the presence of an old stellar population, and that an improper model for the SFH can significantly underestimate the old population in these galaxies. The addition of near-infrared photometry to the constraining data can further distinguish between different SFH model families and significantly tighten the constraints on the mass fraction in the old population. On average more than half of the stellar mass in present-day low-mass galaxies formed at least 8 Gyr ago, while about 30 per cent within the past 4 Gyr. Satellite galaxies on average have formed their stellar mass earlier than central galaxies. The radial dependence of the SFH is quite weak. Our results suggest that most of the low-mass galaxies have an early episode of active star formation that produces a large fraction of their present stellar mass.

Zhou, Shuang↗
DOE OSTI2022

Multifidelity Ensemble Kalman Filtering Using Surrogate Models Defined by Theory-Guided Autoencoders

Data assimilation is a Bayesian inference process that obtains an enhanced understanding of a physical system of interest by fusing information from an inexact physics-based model, and from noisy sparse observations of reality. The multifidelity ensemble Kalman filter (MFEnKF) recently developed by the authors combines a full-order physical model and a hierarchy of reduced order surrogate models in order to increase the computational efficiency of data assimilation. The standard MFEnKF uses linear couplings between models, and is statistically optimal in case of Gaussian probability densities. This work extends the MFEnKF into to make use of a broader class of surrogate model such as those based on machine learning methods such as autoencoders non-linear couplings in between the model hierarchies. We identify the right-invertibility property for autoencoders as being a key predictor of success in the forecasting power of autoencoder-based reduced order models. We propose a methodology that allows us to construct reduced order surrogate models that are more accurate than the ones obtained via conventional linear methods. Numerical experiments with the canonical Lorenz'96 model illustrate that nonlinear surrogates perform better than linear projection-based ones in the context of multifidelity ensemble Kalman filtering. We additionality show a large-scale proof-of-concept result with the quasi-geostrophic equations, showing the competitiveness of the method with a traditional reduced order model-based MFEnKF.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Dataset of Generative AI Workload Power Profiles

This dataset provides a collection of high-resolution (5/10 Hz or every 0.2/0.1 seconds) power consumption profiles for generative artificial intelligence (GenAI) workloads executed on NLR's High Performance Computing (HPC) platform Kestrel. The dataset also includes examples of representative whole-facility power profiles generated using a bottom-up, event-driven, data center energy model . This dataset is designed to support research in energy modeling, infrastructure planning, energy system integration, and sustainability analysis for AI-driven computing systems. The dataset captures time-resolved electrical power measurements across a diverse set of configurations, including variations in job type (inference vs. training), workload (LLM vs. image generation), datasets, and number of compute nodes. Power traces are provided in a standardized format and include both raw/instantaneous and aggregated files. Each profile is accompanied by metadata describing workload parameters, enabling reproducibility and cross-study comparison. The dataset is intended for use in applications such as data center infrastructure planning, energy modeling, demand response and grid impact studies, and development and validation of system-level simulation tools. By making these workload-specific power profiles publicly available, this dataset aims to address the current lack of open, empirical energy data for generative AI systems and to facilitate transparent, reproducible research on the energy and environmental impacts of large-scale AI deployment. If you use this dataset, please cite the associated publication: Vercellino et al., “Measurement of Generative AI Workload Power Profiles for Whole-Facility Data Center Infrastructure Planning,” arXiv:2604.07345 (2026).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

A Graphical Model for Fusing Diverse Microbiome Data

This paper develops a Bayesian graphical model for fusing disparate types of count data. The motivating application is the study of bacterial communities from diverse high-dimensional features, in this case, transcripts, collected from different treatments. In such datasets, there are no explicit correspondences between the communities and each corresponds to different factors, making data fusion challenging. We introduce a flexible multinomial-Gaussian generative model for jointly modeling such count data. This latent variable model jointly characterizes the observed data through a common multivariate Gaussian latent space that parameterizes the set of multinomial probabilities of the transcriptome counts. The covariance matrix of the latent variables induces a covariance matrix of co-dependencies between all the transcripts, effectively fusing multiple data sources. We present a computationally scalable variational Expectation-Maximization (EM) algorithm for inferring the latent variables and the parameters of the model. Here, the inferred latent variables provide a common dimensionality reduction for visualizing the data and the inferred parameters provide a predictive posterior distribution. In addition to simulation studies that demonstrate the variational EM procedure, we apply our model to a bacterial microbiome dataset.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Microphysical Modeling of Fault Slip and Stability Transition in Hydrothermal Conditions

Abstract Field and laboratory observations indicate that the frictional behaviors of faults depend on hydrothermal conditions. We extend the microphysical Chen‐Niemeijer‐Spiers (CNS) model to hydrothermal conditions by using the observed temperature variation of indentation hardness to infer the temperature dependence of a microphysical parameter . This parameter is assumed constant in previous versions of the CNS model. A simple spring‐slider system is used to simulate the fault system and investigate the steady‐state frictional behaviors of wet granite gouges. Our numerical results quantitatively reproduce experimental data showing the frictional‐plastic transition. The results also describe the transition from velocity‐strengthening at low temperatures (<160°C), to velocity‐weakening at intermediate temperatures (160°C–370°C), then back to velocity‐strengthening at high temperatures (>370°C). In our extended CNS model, these results suggest that the dominant shear deformation mechanism does transition from frictional granular flow to fully plastic creep with increasing temperature.

58 GEOSCIENCES↗
DOE OSTI2024

Stomatal Parameters in a Changing Environment

Here, we recommend that stomatal slope parameters (g 1 ) be inferred by inversion so that variations in g 1 may be attributed to variations physiological and environmental conditions. Understanding g1 will advance predictions of plant gas exchange and performance under global climate.

Stomatal conductance modeling↗
DOE OSTI2025

The DESI Y1 RR Lyrae catalog I: Empirical modeling of the cyclic variation of spectroscopic properties and a chemodynamical analysis of the outer halo

We present the catalog of RR Lyrae stars observed in the first year of operations of the Dark Energy Spectroscopic Instrument (DESI) survey. This catalog contains 6,240 RR Lyrae stars out to $\sim100$ kpc from the Galactic center and over 12,000 individual epochs with homogeneously-derived stellar atmospheric parameters. We introduce a novel methodology to model the cyclical variation of the spectroscopic properties of RR Lyrae from single-epoch measurements. We employ this method to infer the radial velocity and effective temperature variation of fundamental mode and first-overtone RR Lyrae stars and to determine their systemic velocities and mean temperatures. For fundamental mode pulsators, we obtain radial velocity curves with amplitudes of $\sim$30--50 km s$^{-1}$ and effective temperature curves with 300--1,000 K variations, whereas for first-overtone pulsators these amplitudes are $\sim20$ km s$^{-1}$ and $\sim 600$ K, respectively. We use our sample to study the metallicity distribution of the halo and its dependence on Galactocentric distance ($R_{\rm GC}$). Using a radius-dependent mixture model, we split the data into chemodynamically distinct components and find that our inner halo sample ($R_{\rm GC}\lesssim50$ kpc) is predominantly composed of stars with [Fe/H] $\sim-1.5$ dex and largely radial orbits (with an anisotropy parameter $\beta\sim0.94$), that we associate with the Gaia-Sausage-Enceladus merger event. Stars in the outer halo field exhibit a broader and more metal-poor [Fe/H] distribution with more circular orbits ($\beta\sim0.39$). The metallicity gradient of the metal-rich and the metal-poor components is found to be $0.005$ and $0.010$ dex kpc$^{-1}$, respectively. Our catalog highlights DESI's tantalizing potential for studying the Milky Way and the pulsation properties of RR Lyrae stars in the era of large spectroscopic surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

When action is not least for systems with action-dependent Lagrangians

The dynamics of some non-conservative and dissipative systems can be derived by calculating the first variation of an action-dependent action according to the variational principle of Herglotz. This is directly analogous to the variational principle of Hamilton commonly used to derive the dynamics of conservative systems. In a similar fashion, just as the second variation of a conservative system’s action can be used to infer whether that system’s possible trajectories are dynamically stable, so too can the second variation of the action-dependent action be used to infer whether the possible trajectories of non-conservative and dissipative systems are dynamically stable. In this paper, I show, generalizing earlier analyses of the second variation of the action for conservative systems, how to calculate the second variation of the action-dependent action and how to apply it to two physically important systems: a time-independent harmonic oscillator and a time-dependent harmonic oscillator.

Physics↗
DOE OSTI2021

Surveying Grassland Islands: the genetics and performance of Appalachian switchgrass (Panicum virgatum L.) collections

The interior Southeastern United States could contain novel germplasm for the bioenergy crop switchgrass due to its diverse habitats and geographic location between genetic subpopulations (Atlantic, Midwest, and Gulf). Collections from this region could accelerate breeding progress, contribute to conservation efforts, and improve understanding of isolated grasslands in the region. This study located 22 sites in the Midsouth region and obtained 1,521,210 single nucleotide polymorphism markers of 202 individuals through genotype-by-sequencing. Here, individuals were evaluated for flowering time, winter survival and tiller number. Comparison to a national diversity panel revealed that branches of two major subpopulations occur in the region with two levels of polyploidy: Atlantic tetraploids and Midwest octoploids. Two locations contained admixed octoploid individuals with Midwest and Gulf genetics. Field performance of the Midwest octoploids conformed with prior reported performance of the Midwest subpopulation, although three sites contained promising late flowering traits. The Atlantic tetraploids had moderate winter survival, short stature, and anomalously early flowering. Atlantic populations mostly occurred in marginal sites and their morphological and flowering time adaptations may be a resource conservation strategy. Demographic inference of historical effective population size variation in a subset of tetraploid locations indicated a widespread recent decline in effective population size. This pattern is consistent with isolation of these switchgrass communities from larger populations and is further supported by evidence of inbreeding within the populations (F I = 0.18). The populations documented in this study contain novel genetic diversity and adaptations to a range of marginal habitats. Therefore, this study provides a new source of germplasm for future breeding and conservation programs.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Toward In Situ Synchrotron Mapping of Crystal Selection Processes during Crystal Growth

In this work, we present an automated and rapid method for non-destructive mapping of crystal grains in a rod-shaped sample. The approach was designed for application to in situ float-zone crystal growth experiments at an x-ray synchrotron source, but could be useful in other applications. The methods have been tested on a TiO 2 boule grown in an optical float zone furnace. The approach applies a statistical filter to polycrystalline diffraction patterns on 2D detectors to rapidly determine the degree of powder-quality 1of the signal. When larger crystals emerge in the growth, their position, size and shape can be tracked using an automated blob-tracking algorithm that follows individual Bragg peaks as a function of position in a grid-scan, even when multiple crystals are contributing spots to diffraction images. This method is found to be robust as the same crystal shape can be independently reconstructed using different sets of Bragg reflections. Image segmentation methods are then used to map out the polycrystalline grains. We also note that other information about crystal quality, such as mosaicity or strain state, may be inferred and mapped from the intensity variation of the Bragg peaks at different locations within the sample.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Spin skyrmion gaps as signatures of strong-coupling insulators in magic-angle twisted bilayer graphene

Abstract The flat electronic bands in magic-angle twisted bilayer graphene (MATBG) host a variety of correlated insulating ground states, many of which are predicted to support charged excitations with topologically non-trivial spin and/or valley skyrmion textures. However, it has remained challenging to experimentally address their ground state order and excitations, both because some of the proposed states do not couple directly to experimental probes, and because they are highly sensitive to spatial inhomogeneities in real samples. Here, using a scanning single-electron transistor, we observe thermodynamic gaps at even integer moiré filling factors at low magnetic fields. We find evidence of a field-tuned crossover from charged spin skyrmions to bare particle-like excitations, suggesting that the underlying ground state belongs to the manifold of strong-coupling insulators. From the spatial dependence of these states and the chemical potential variation within the flat bands, we infer a link between the stability of the correlated ground states and local twist angle and strain. Our work advances the microscopic understanding of the correlated insulators in MATBG and their unconventional excitations.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Planck intermediate results. LVI. Detection of the CMB dipole through modulation of the thermal Sunyaev-Zeldovich effect: Eppur si muove II

The largest temperature anisotropy in the cosmic microwave background (CMB) is the dipole, which has been measured with increasing accuracy for more than three decades, particularly with the Planck satellite. The simplest interpretation of the dipole is that it is due to our motion with respect to the rest frame of the CMB. Since current CMB experiments infer temperature anisotropies from angular intensity variations, the dipole modulates the temperature anisotropies with the same frequency dependence as the thermal Sunyaev-Zeldovich (tSZ) effect. Here, we present the first, and significant, detection of this signal in the tSZ maps and find that it is consistent with direct measurements of the CMB dipole, as expected. The signal contributes power in the tSZ maps, which is modulated in a quadrupolar pattern, and we estimate its contribution to the tSZ bispectrum, noting that it contributes negligible noise to the bispectrum at relevant scales.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Variational autoencoders for at-source data reduction and anomaly detection in high energy particle detectors

Detectors in next-generation high-energy physics experiments face several daunting requirements, such as high data rates, damaging radiation exposure, and stringent constraints on power, space, and latency. To address these challenges, machine learning in readout electronics can be leveraged for smart detector designs, enabling intelligent inference and data reduction at-source. Variational autoencoders (VAEs) offer a variety of benefits for front-end readout; an on-sensor encoder can perform efficient lossy data compression while simultaneously providing a latent space representation that can be used for anomaly detection. Results are presented from low-latency and resource-efficient VAEs for front-end data processing in a futuristic silicon pixel detector. Encoder-based data compression is found to preserve good performance of off-detector analysis while significantly reducing the off-detector data rate as compared to a similarly sized data filtering approach. Furthermore, the latent space information is found to be a useful discriminator in the context of real-time sensor defect monitoring. Together, these results highlight the multifaceted utility of autoencoder-based front-end readout schemes and motivate their consideration in future detector designs.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

A species’ response to spatial climatic variation does not predict its response to climate change

The dominant paradigm for assessing ecological responses to climate change assumes that future states of individuals and populations can be predicted by current, species-wide performance variation across spatial climatic gradients. However, if the fates of ecological systems are better predicted by past responses to in situ climatic variation through time, this current analytical paradigm may be severely misleading. Empirically testing whether spatial or temporal climate responses better predict how species respond to climate change has been elusive, largely due to restrictive data requirements. Here, we leverage a newly collected network of ponderosa pine tree-ring time series to test whether statistically inferred responses to spatial versus temporal climatic variation better predict how trees have responded to recent climate change. When compared to observed tree growth responses to climate change since 1980, predictions derived from spatial climatic variation were wrong in both magnitude and direction. This was not the case for predictions derived from climatic variation through time, which were able to replicate observed responses well. Future climate scenarios through the end of the 21st century exacerbated these disparities. These results suggest that the currently dominant paradigm of forecasting the ecological impacts of climate change based on spatial climatic variation may be severely misleading over decadal to centennial timescales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Benchmarking of Different Inverse Point Kinetics Implementations for an Autocorrected Reactimeter Algorithm

In November 2017, the Transient Reactor Test Facility returned to operation. Since that time, many transient test series have been completed; each has provided valuable data for materials performance, reactor safety that can be applied in future designs. During each experimental series, detector count rates provided important information on the core behavior during transients. However, a limitation of these data is that variations in the neutron distribution during experiments cause errors when attempting to infer reactivity evolution from detector signals. Neutron physics codes can be used to compute the flux shape variations. However, this is a poor solution when the experimental data is used to do verification, validation and uncertainty quantification (VVUQ) on codes. Indeed, if the output of the code is used both as a reference and to correct what the reference is compared to, the circular dependency limits the quality of the VVUQ approach. To overcome this problem, an Autocorrected Reactimeter Algorithm (ACRA) has been developed. This approach infers time-dependent reactivity evolution by testing different spatial corrections and selecting the one that minimizes reactivity variations when the core is in a frozen configuration (i.e. when there is no variation in parameters affecting reactivity). However, the scope of this method was limited to transients where there were negligible thermal feedback. Indeed, the core is never in a frozen configuration when the fuel temperature varies during the whole transient. This is our motivation for the development of an improved version of the ACRA which does not require frozen configurations

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Monitoring of ground water table depth and soil moisture at the Point Reyes field site

Ground water table (GWT) depth and soil moisture (SM) have been monitored at several locations at the Point Reyes field site (Californian coastal grassland) from 2021 to 2024. Monitoring is still on-going and data may be added to this archive at later time. The SM data have been acquired using Teros 12 Meter soil moisture sensors placed at 10, 30, 60 and 90 cm depth at 5 locations along a small hillslope. These sensors also collect soil temperature and bulk conductance. In addition, some collocated sensors provide pore pressure and Photochemical Reflectance Index (PRI). The GWT depth has been inferred from various type of Onset pressure transducers. The pressure measurements have been corrected for atmospheric pressure variations and sensor position relative to the ground surface to infer GWT depth, as well as with RTK GPS data to infer GWT elevation. The GWT data have been acquired at 5 distinct locations from 2020 to 2024 with the sensors placed at about 4 m depth. In addition, GWT data has been acquired for the 2023-2024 period with sensors located in 1 m deep shallow wells installed near each deeper well. This data is intended to evaluate possibly different dynamic in shallow (perched) and deep aquifer. The datasets are all provided in csv format. Please note that the interpretation of the GWT data needs to be done with consideration of environmental and well characteristics at the site and uncertainty in various variables. For more information on GWT and SM data, please contact the author.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

(Doublon) Benchmarking of Different Inverse Point Kinetics Implementations for an Autocorrected Reactimeter Algorithm

In November 2017, the Transient Reactor Test Facility returned to operation. Since that time, many transient test series have been completed, such as the Transient Heatsink Overpower Response capsule (THOR), the Transient Water Irradiation System for TREAT (TWIST), and Sirius. Each has provided valuable data for materials performance and reactor safety that can be applied in future designs. During each experimental series, detector count rates provided important information on the core behavior during transients. However, a limitation of these data is that variations in the neutron distribution during experiments can cause errors when attempting to infer reactivity evolution from detector signals. Neutron physics codes can be used to compute the flux shape variations. However, this is a poor solution when the experimental data is used for code verification, validation and uncertainty quantification. Indeed, if the output of the code is used both as a reference and to correct what the reference is compared to, the circular dependency limits the quality of the verification, validation and uncertainty quantification approach. To overcome this problem, the autocorrected reactimeter algorithm (ACRA) has been developed. This approach infers a time-dependent reactivity evolution by testing different spatial corrections and selecting the one that minimizes reactivity variations when the core is in a frozen configuration (i.e., when there is no variation in parameters affecting reactivity). However, the scope of this method was limited to transients where there were negligible thermal feedback. Indeed, the core is never in a frozen configuration when the fuel temperature varies during the whole transient. This is our motivation for developing an improved version of the ACRA that does not require frozen configurations. To develop this new algorithm, we need a precise and unbiased implementation of the inverse point kinetic equations (IPKEs) as any error in the reactivity evaluation will be propagated into the choice of the optimal spatial correction. Indeed, the previous reactimeter algorithm would use approximations, such as a negligible flux amplitude derivative, to focus on rapidity. For the numerical validation of ACRA, we aim at absolute error under for reactivity derived from signals similar to the one of this study. In this summary, we test eight different IPKE implementations. Each will process a mockup signal built for this study, similar to those that the future ACRA will process. Each reactivity output will be compared to the reference reactivity that has been used to generate the mockup signal. The implementation minimizing the difference with the reference reactivity will be used in the development of a new ACRA formulation.

73 - NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Your Clean Graphene is Still Not Clean

Researchers working with thin samples, such as monolayer graphene, are consistently struggling against contamination. Indeed, the problem of hydrocarbon contamination is known from the earliest days of electron microscopy and efforts to reduce this problem are ubiquitous to almost all high-vacuum experiments. Accurate knowledge of the behavior of such contamination is essential for electron beam (e-beam) based atomic fabrication, where it is aspired to select and control matter on an atom-by-atom basis. Here, the vexing question of hydrocarbon contamination on graphene is taken up. Image intensity is used to directly reveal the presence of diffusing hydrocarbons on ostensibly clean graphene. These diffusing hydrocarbons are previously inferred but not directly observed. Surprising dynamic variations of the concentration of these hydrocarbons impels questions about their origin. Here, some possible explanations are presented and some tentative conclusions are drawn. This work updates the conceptual model of “clean graphene” and offers refinements to the description of e-beam induced hydrocarbon deposition.

atomic fabrication↗