Search NASA⌕ Search

SEARCH · Search NASA

Results for “dynamic 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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Lift & Learn: Physics-informed machine learning for large-scale nonlinear dynamical systems

In this work, we present Lift & Learn, a physics-informed method for learning low-dimensional models for large-scale dynamical systems. The method exploits knowledge of a system’s governing equations to identify a coordinate transformation in which the system dynamics have quadratic structure. This transformation is called a lifting map because it often adds auxiliary variables to the system state. The lifting map is applied to data obtained by evaluating a model for the original nonlinear system. This lifted data is projected onto its leading principal components, and low-dimensional linear and quadratic matrix operators are fit to the lifted reduced data using a least-squares operator inference procedure. Analysis of our method shows that the Lift & Learn models are able to capture the system physics in the lifted coordinates at least as accurately as traditional intrusive model reduction approaches. This preservation of system physics makes the Lift & Learn models robust to changes in inputs. Numerical experiments on the FitzHugh–Nagumo neuron activation model and the compressible Euler equations demonstrate the generalizability of our model.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Unlocking hidden information in sparse small-angle neutron scattering measurements

Hypothesis Small-Angle Neutron Scattering (SANS) is a powerful technique for studying soft matter systems such as colloids, polymers, and lyotropic phases, providing nanoscale structural insights. However, its effectiveness is limited by low neutron flux, leading to long acquisition times and noisy data. Here, we hypothesize that Bayesian statistical inference using Gaussian Process Regression (GPR) can reconstruct high-fidelity scattering data from sparse measurements by leveraging intensity smoothness and continuity. Experiments and Simulations The method was benchmarked computationally and validated through SANS experiments on various soft matter systems, including wormlike micelles, colloidal suspensions, polymeric structures, and lyotropic phases. GPR-based inference was applied to both experimental and synthetic data to evaluate its effectiveness in noise reduction and intensity reconstruction. Findings GPR significantly enhances SANS data quality and therefore reducing measurement times by up to two orders of magnitude. This cost-effective approach maximizes experimental efficiency, enabling high-throughput studies and real-time monitoring of dynamic systems. It is particularly beneficial for weakly scattering and time-sensitive studies. Beyond SANS, this framework applies to other low-SNR techniques, including laboratory-based small-angle X-ray scattering and various dynamical scattering methods. Furthermore, it offers transformative potential for compact neutron sources, enhancing their viability for structural analysis in resource-limited settings.

Small angle neutron scattering↗
DOE OSTI2022

Sparsifying priors for Bayesian uncertainty quantification in model discovery

We propose a probabilistic model discovery method for identifying ordinary differential equations governing the dynamics of observed multivariate data. Our method is based on the sparse identification of nonlinear dynamics (SINDy) framework, where models are expressed as sparse linear combinations of pre-specified candidate functions. Promoting parsimony through sparsity leads to interpretable models that generalize to unknown data. Instead of targeting point estimates of the SINDy coefficients, we estimate these coefficients via sparse Bayesian inference. The resulting method, uncertainty quantification SINDy (UQ-SINDy), quantifies not only the uncertainty in the values of the SINDy coefficients due to observation errors and limited data, but also the probability of inclusion of each candidate function in the linear combination. UQ-SINDy promotes robustness against observation noise and limited data, interpretability (in terms of model selection and inclusion probabilities) and generalization capacity for out-of-sample forecast. Sparse inference for UQ-SINDy employs Markov chain Monte Carlo, and we explore two sparsifying priors: the spike and slab prior, and the regularized horseshoe prior. UQ-SINDy is shown to discover accurate models in the presence of noise and with orders-of-magnitude less data than current model discovery methods, thus providing a transformative method for real-world applications which have limited data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Estimating ionization states and continuum lowering from ab initio path integral Monte Carlo simulations for warm dense hydrogen

Warm dense matter (WDM) is an active field of research, with applications ranging from astrophysics to inertial confinement fusion. Ionization degree and continuum lowering are important quantities to understand how materials behave under these conditions, but can be difficult to diagnose since experimental campaigns are limited and often require model-dependent analysis. This is especially true for hydrogen, which has a comparably low scattering cross section, making high-quality data particularly difficult to obtain. Consequently, building equation of state tables often relies on simulations in combination with untested approximations to extract properties from experiments. Here, we investigate an approach for extracting the ionization potential depression and ionization degree—quantities which are otherwise not directly accessible from the physical model—from first-principles path integral Monte Carlo (PIMC) simulations utilizing a chemical model. In contrast to experimental measurements, where noise and nonequilibrium effects add to the uncertainty of the inferred parameters, PIMC simulations provide a clean signal with well-defined thermodynamic conditions. Comparisons against commonly used models show a qualitative agreement, but we find deviations primarily for the high-density and high-temperature cases. We also demonstrate the decreasing sensitivity of the dynamic structure factor with respect to both ionization and continuum lowering for increasing scattering angles in x-ray Thomson scattering experiments. Our work has important implications for the design of future experiments, but also offers qualitative understanding of structure factors and the imaginary-time correlation function obtained from first-principles quantum Monte Carlo simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Indicator-directed Dynamic Power Management for Iterative Workloads on GPU-Accelerated Systems

Modern high-performance and warehouse computing centers show strong interest in minimizing system power consumption while satisfying customers’ quality of service (QoS). Dynamic voltage and frequency scaling (DVFS) is effective for achieving this goal. Nevertheless, automating the process online and making it transparent to users must address three major challenges: (1) Complexity — today’s hardware components (e.g., CPUs, GPUs, memory, network, etc.) can be configured in several or dozens of frequency/voltage states for satisfying divergent system demands. Given their combination and the emergence of heterogeneity, searching the optimal configuration in the design space online can be timing consuming. (2) QoS guarantee — user-defined objectives such as power constraint and performance target must be monitored, predicted and ensured at the best effort. (3) Adaptability — various known and unknown workloads run on systems. Workloads characteristics should be quickly determined and configurations dynamically adjusted in accord with workloads and QoS. In this work, we focus on applications exhibiting an interesting feature – iterative or periodic, which is common among conventional HPC and emerging machine learning workloads. We propose an online dynamic power-performance (ODPP) management framework to dynamically adjust GPU DVFS configurations to meet performance and power objectives and constraints, without any code annotation or intrusion. Particularly, ODPP extracts the performance and power indicators for applications from their resources utilization profiles in a short episode. It further automatically constructs an accurate model that infers from the indicators how the application's performance and power vary with GPU core and memory frequencies. Aided with the model, for both seen and unseen applications, ODPP can quickly determine the most appropriate DVFS configuration for their execution. We evaluate ODPP on an NVIDIA GPU using multiple exascale computing (ECP) and deep learning applications.

Zou, Pengfei↗
DOE OSTI2022

Initial soil organic carbon stocks govern changes in soil carbon: Reality or artifact?

Abstract Changes in soil organic carbon (SOC) storage have the potential to affect global climate; hence identifying environments with a high capacity to gain or lose SOC is of broad interest. Many cross‐site studies have found that SOC‐poor soils tend to gain or retain carbon more readily than SOC‐rich soils. While this pattern may partly reflect reality, here we argue that it can also be created by a pair of statistical artifacts. First, soils that appear SOC‐poor purely due to random variation will tend to yield more moderate SOC estimates upon resampling and hence will appear to accrue or retain more SOC than SOC‐rich soils. This phenomenon is an example of regression to the mean. Second, normalized metrics of SOC change—such as relative rates and response ratios—will by definition show larger changes in SOC at lower initial SOC levels, even when the absolute change in SOC does not depend on initial SOC. These two artifacts create an exaggerated impression that initial SOC stocks are a major control on SOC dynamics. To address this problem, we recommend applying statistical corrections to eliminate the effect of regression to the mean, and avoiding normalized metrics when testing relationships between SOC change and initial SOC. Careful consideration of these issues in future cross‐site studies will support clearer scientific inference that can better inform environmental management.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A Step in Understanding the Hubble Tension

As cosmological data have improved, tensions have arisen. One such tension is the difference between the locally measured Hubble constant H 0 and the value inferred from the cosmic microwave background (CMB). Interacting radiation has been suggested as a solution, but studies show that conventional models are precluded by high-ℓ CMB polarization data. It seems at least plausible that a solution may be provided by related models that distinguish between high- and low-ℓ multipoles. When interactions of strongly-coupled radiation are mediated by a force-carrier that becomes non-relativistic, the dark radiation undergoes a "step" in which its relative energy density increases as the mediator deposits its entropy into the lighter species. If this transition occurs while CMB-observable modes are inside the horizon, high- and low-ℓ peaks are impacted differently, corresponding to modes that enter the horizon before or after the step. Here, these dynamics are naturally packaged into the simplest supersymmetric theory, the Wess-Zumino model, with the mass of the scalar mediator near the eV-scale. We investigate the cosmological signatures of such "Wess-Zumino Dark Radiation" (WZDR) and find that it provides an improved fit to the CMB alone, favoring larger values of H 0 . If supernovae measurements from the SH0ES collaboration are also included in the analysis, the inferred value of H 0 is yet larger, but the preference for dark radiation and the location of the transition is left nearly unchanged. Utilizing a standardized set of measures, we compare to other models and find that WZDR is among the most successful at addressing the H 0 tension and the best of those with a Lagrangian formulation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Impact of mid- Z gas fill on dynamics and performance of shock-driven implosions at the OMEGA laser

Shock-driven implosions with 100% deuterium (D 2 ) gas fill compared to implosions with 50:50 nitrogen-deuterium (N 2 ⁢D 2 ) gas fill have been performed at the OMEGA laser facility to test the impact of the added mid-Ζ fill gas on implosion performance. Ion temperature (Τ ion ) as inferred from the width of measured DD-neutron spectra is seen to be 34%±6% higher for the N 2⁢ D 2 implosions than for the D 2 -only case, while the DD-neutron yield from the D 2 -only implosion is 7.2±0.5 times higher than from the N 2⁢ D 2 gas fill. The T ion enhancement for N 2 ⁢D 2 is observed in spite of the higher Z, which might be expected to lead to higher radiative loss, and higher shock strength for the D 2 -only versus N 2 ⁢D 2 implosions due to lower mass, and is understood in terms of increased shock heating of N compared to D, heat transfer from N to D prior to burn, and limited amount of ion-electron-equilibration-mediated additional radiative loss due to the added higher-Z material. Further, this picture is supported by interspecies equilibration timescales for these implosions, constrained by experimental observables. The one-dimensional (1D) kinetic Vlasov-Fokker-Planck code ifp and the radiation hydrodynamic simulation codes hyades (1D) and xrage [1D, two-dimensional (2D)] are brought to bear to understand the observed yield ratio. Comparing measurements and simulations, the yield loss in the N 2 ⁢D 2 implosions relative to the pure D 2 -fill implosion is determined to result from the reduced amount of D 2 in the fill (fourfold effect on yield) combined with a lower fraction of the D 2 fuel being hot enough to burn in the N 2 ⁢D 2 case. The experimental yield and T ion ratio observations are relatively well matched by the kinetic simulations, which suggest interspecies diffusion is responsible for the lower fraction of hot D 2 in the N 2 ⁢D 2 relative to the D 2 -only case. The simulated absolute yields are higher than measured; a comparison of 1D versus 2D XRAGE simulations suggest that this can be explained by dimensional effects. The hydrodynamic simulations suggest that radiative losses primarily impact the implosion edges, with ion-electron equilibration times being too long in the implosion cores. The observations of increased T ion and limited additional yield loss (on top of the fourfold expected from the difference in D content) for the N 2 ⁢D 2 versus D 2 -only fill suggest it is feasible to develop the platform for studying CNO-cycle-relevant nuclear reactions in a plasma environment.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

GradeIT (Road Grade Inference Tool)

The National Renewable Energy Laboratory (NREL) has developed the Road Grade Inference Tool, "GradeIT". This tool is a modular Python package the leverages the United States Geological Survey’s (USGS) 1/3 arc-second Digital Elevation Model (DEM), an open and free public dataset. Roadway gradient (slope) has a strong influence on vehicle energy consumption behavior and longitudinal dynamics. Therefore, any tools, models, or analyses that consider vehicle energy consumption or driving profiles is incomplete without accurate elevation and grade information. Despite the acute need for this data, it is rarely collected alongside GPS (latitude and longitude) data, and tools to post-process GPS data to append elevation and grade are difficult to come by. GradeIT satisfies this significant need in the research community. Utilization of the USGS DEM makes the underlying data available to any potential user and ensures continuous data anywhere in the continental United States. The USGS DEM is a great resource; however, in order to be useful for automotive and traffic engineering applications, significant interpolation and filtering is required to achieve an accurate representation of gradient on the roadway. GradeIT employs a two-dimensional interpolation to approximate elevation at a point on the road from the nearest nodes in the DEM. Then, GradeIT applies multiple filtering techniques to smooth the elevation profile to a realistic road topology. Finally, road grade is derived from the elevation profile and returned as the primary output.

Holden, Jacob↗
DOE OSTI2020

FY20 Proxy App Suite Release: Report for ECP Proxy App Project Milestone ADCD-504-10

Version 4.0 of the ECP Proxy App Suite is practically unchanged from the previous release. The current set of proxies has proven useful for many aspects of benchmarking and co-design and we see little reason to alter the suite. Although there have been few changes to the ECP suite, the team has been hard at work in other areas. In the area of Machine Learning (ML) we have now created a separate proxy suite dedicated to this scientific applications of ML. The suite includes: miniGAN (Generative Adversarial Networks), miniRL (Reinforcement Learning), CRADL (inline inference), Cosmoflow-Benchmark (Convolutional Neural Network), and MLPerf-DeepCam (Climate Segmentation Benchmark). Section 3 contains more information about these proxies as well as the principles that are guiding the development of the suite. We have surveyed available proxies for several application domains including Computational Fluid Dynamics, Quantum Chemistry, Quantum Computing Simulation, Molecular Dynamics, Monte Carlo Transport, and Density Functional Theory to identify gaps in proxy coverage. Several new proxy apps are either already available or will be released soon to fill these gaps. Section 4 provides full details. Finally, section 5 reports on our continued collaboration with the ECP Continuous Integration (CI) effort to use proxy apps to help identify problems and roadblocks to cross-lab CI. We have also assisted the El Capitan Center of Excellence (COE) to stand up a CI system that will be used to test releases of HPE and AMD software stacks. We hope that the COE effort can serve as a model for ECP by showing how the Proxy App Team can work with AD and ST teams identify critical features, kernels, patterns, etc. and incorporate them into a CI system that will help ensure that Frontier and Aurora will provide those needed capabilities.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Imaging electron angular distributions to assess a full-power petawatt-class laser focus

Here we present a technique to assess the focal volume of petawatt-class lasers at full power. Our approach exploits quantitative measurement of the angular distribution of electrons born in the focus via ionization of rarefied gas, which are accelerated forward and ejected ponderomotively by the field. We show that a bivariate (θ,Φ) angular distribution, which was obtained with image plates, not only enables the peak intensity to be extracted, but also reflects nonideality of the focal-spot intensity distribution. In our prototype demonstration at intensities of a few ×10 19 to a few ×10 20 W/cm 2 , an f/10 optic produced a focal spot in the paraxial regime. This allows a plane-wave parametrization of the peak intensity given by tanθ c =2/a 0 (a 0 being the normalized vector potential and θ c the minimum ejection angle) to be compared with our measurements. Qualitative agreement was found using an a 0 inferred from the pulse energy, pulse duration, and focal spot distribution with a modified parametrization, tanθ c =2η/a 0 ($η=2.02^{+0.26}_{–0.22}$). This highlights the need for (i) better understanding of intensity degradation due to focal-spot distortions and (ii) more robust modeling of the ejection dynamics. Using single-shot detection of electrons, we showed that while there is significant shot-to-shot variation in the number of electrons ejected at a given angular position, the average distribution scales with the pulse energy in a way that is consistent with that seen with the image plates. Finally, we note that the asymptotic behavior as θ→0° limits the usability of angular measurement. For 800 nm, this limit is at an intensity ~10 21 W/cm 2 .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Physics-based modeling and information-theoretic sensor and settings selection for tool wear detection in precision machining

Precision machining of metals is an energy intensive process with applications and impacts across the manufacturing industry. The energy efficiency, product yield, and maintenance of the precision machine require a digital twin that can assist with prognostics and health management. In this report a physics-based model is developed and validated against face milling data, and then used for the timely and precise inference of machining faults that cannot be measured directly. Computer numerical control (CNC) measurements of power and force are used through this physics-based machining model to predict deviations of the outputs of power consumption and cutting forces during normal operation. A model-based fault detection and isolation methodology is applied to determine the optimal (traditional and available) sensor suite and the test settings (admissible input values) that improve the inference of tool wear in face milling. The optimal sensor suite and input test settings are obtained by solving a mixed integer non-linear program that optimizes information-theoretic metrics relevant to the detection and isolation of tool wear from steady-state or transient machining measurements. Dynamic time warping and k—NN classification are then used to validate the robustness of the optimal design for fault detection test design, including the optimal sensor suite.

42 ENGINEERING↗
DOE OSTI2025

Disorder-induced spin excitation continuum and spin-glass ground state in the inverse spinel CuGa 2 O 4

Spinel-structured compounds serve as prototypical examples of highly frustrated systems and are promising candidates for realizing the long-sought quantum spin liquid (QSL) state. However, structural disorder is inevitable in many real QSL candidates and its impact remains a topic of intense debate. In this work, we conduct comprehensive investigations on CuGa 2 ⁢ O 4 , a spinel compound with significant structural disorder, focusing on its thermodynamic properties and spectroscopic behaviors. No long-range magnetic order is observed down to ∼80 mK, as evidenced by magnetic susceptibility, specific-heat, and elastic neutron scattering measurements. More intriguingly, inelastic neutron scattering experiments reveal a broad gapless continuum of magnetic excitations around the Brillouin zone boundary, resembling the magnetic excitation spectra expected for a QSL. Nevertheless, a spin-freezing transition at 𝑇 f ≈ 0.88 K is identified from the cusp in the dc susceptibility curves, where a bifurcation between zero-field-cooling and field-cooling curves occurs. Furthermore, ac susceptibility measurements show a peak close to 𝑇 f at low frequency, which shifts to higher temperature with increasing frequency. These results show that CuGa 2 ⁢ O 4 has a spin-glass ground state, consistent with the establishment of short-range order inferred from the specific-heat measurements. Collectively, these results illustrate the crucial role of disorder in defining the excitation spectrum out of the disordered ground state. Furthermore, our findings shed light onto the broader class of 𝐴⁢𝐵 2 ⁢O 4 spinels and advance our understanding of the spin dynamics in magnetically disordered systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Parameter estimation for X-ray scattering analysis with Hamiltonian Markov Chain Monte Carlo

Bayesian-inference-based approaches, in particular the random-walk Markov Chain Monte Carlo (MCMC) method, have received much attention recently for X-ray scattering analysis. Hamiltonian MCMC, a state-of-the-art development in the field of MCMC, has become popular in recent years. It utilizes Hamiltonian dynamics for indirect but much more efficient drawings of the model parameters. We described the principle of the Hamiltonian MCMC for inversion problems in X-ray scattering analysis by estimating high-dimensional models for several motivating scenarios in small-angle X-ray scattering, reflectivity, and X-ray fluorescence holography. Hamiltonian MCMC with appropriate preconditioning can deliver superior performance over the random-walk MCMC, and thus can be used as an efficient tool for the statistical analysis of the parameter distributions, as well as model predictions and confidence analysis.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

libEnsemble: A complete Python toolkit for dynamic ensembles of calculations

Almost all science and engineering applications eventually stop scaling: their runtime no longer decreases as available computational resources increase. Therefore, many applications will struggle to efficiently use emerging extreme-scale high-performance, parallel, and distributed systems. libEnsemble is a complete Python toolkit and workflow system for intelligently driving ensembles of experiments or simulations at massive scales. It enables and encourages multidisciplinary design, decision, and inference studies portably running on laptops, clusters, and supercomputers.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

ARM-IRL: Adaptive Resilience Metric Quantification Using Inverse Reinforcement Learning

The resilience of safety-critical systems is gaining importance due to the rise in cyber and physical threats, especially within critical infrastructure. Traditional static resilience metrics may not capture dynamic system states, leading to inaccurate assessments and ineffective responses to cyber threats. This work aims to develop a data-driven, adaptive method for resilience metric learning. We propose a data-driven approach using inverse reinforcement learning (IRL) to learn a single, adaptive resilience metric. The method infers a reward function from expert control actions. Unlike previous approaches using static weights or fuzzy logic, this work applies adversarial inverse reinforcement learning (AIRL), training a generator and discriminator in parallel to learn the reward structure and derive an optimal policy. The proposed approach is evaluated on multiple scenarios: optimal communication network rerouting, power distribution network reconfiguration, and cyber–physical restoration of critical loads using the IEEE 123-bus system. The adaptive, learned resilience metric enables faster critical load restoration in comparison to conventional RL approaches.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Bayesian calibration of bubble size dynamics applied to CO 2 gas fermenters

To accelerate the scale-up of gaseous CO 2 fermentation reactors, computational models need to predict gas-to-liquid mass transfer which requires capturing the bubble size dynamics, i.e. bubble breakup and coalescence. However, the applicability of existing models beyond air–water mixtures remains to be established. Here, an inverse modeling approach, accelerated with a neural network surrogate, calibrates the breakup and coalescence closure models, that are used in class methods for population balance modeling (PBM). The calibration is performed based on experimental results obtained in a CO 2 -air–water-coflowing bubble column reactor. Bayesian inference is used to account for noise in the experimental dataset and bias in the simulation results. To accurately capture gas holdup and interphase mass transfer, the results show that the breakage rate needs to be increased by one order of magnitude. In conclusion, the inferred model parameters are then used on a separate configuration and shown to also improve bubble size distribution predictions.

09 BIOMASS FUELS↗
DOE OSTI2025

Array-Based Seismic Measurements of OSIRIS-REx’s Re-Entry

The return home of the OSIRIS-REx spacecraft in September 2023 marked only the fifth time that an artificial object entered the Earth’s atmosphere at interplanetary velocities. Although rare, such events serve as valuable analogs for natural meteoroid re-entries; enabling study of hypersonic dynamics, shock wave generation, and acoustic-to-seismic coupling. Here, in this study, we report on the signatures recorded by a dense (100 m scale) 11-station array located almost directly underneath the capsule’s point of peak atmospheric heating in northern Nevada. Seismic data are presented, which allow inferences to be made about the shape of the shock wave’s footprint on the surface, the capsule’s trajectory, and its flight parameters.

58 GEOSCIENCES↗