Search NASA⌕ Search

SEARCH · Search NASA

Results for “Model 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 1,171 records · Page 65

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↗

Reduced‐Order Modeling of Energetic Materials Using Physics‐Aware Recurrent Convolutional Neural Networks in a Latent Space (LatentPARC)

Physics-aware deep learning (PADL) has gained popularity for use in spatiotemporal dynamics simulations, such as those in computational modeling of energetic materials (EM). We show that the challenge PADL methods face while learning complex field evolution problems can be simplified and accelerated by decoupling it into two tasks: learning complex geometric features in evolving fields and modeling dynamics over these features in a lower-dimensional feature space. We build upon our previous work on physics-aware recurrent convolutional neural networks (PARC). PARC embeds knowledge of underlying physics into its neural network architecture for more robust and accurate prediction of evolving physical fields. PARC was shown to effectively learn complex nonlinear features such as the formation of hotspots and coupled shock fronts in various initiation scenarios of EMs, as a function of microstructures, serving effectively as a microstructure-aware burn model. Here, we further accelerate PARC and reduce its computational cost by projecting the original dynamics onto a lower-dimensional invariant manifold, or “latent space.” The projected latent representation encodes the complex geometry of evolving fields (e.g., temperature and pressure) in a set of data-driven features. The reduced dimension of this latent space allows us to learn the dynamics during the initiation of EM with a lighter and more efficient model. We observe a significant decrease in training and inference time while maintaining results comparable to PARC at inference. This work takes steps towards enabling rapid prediction of EM thermomechanics at larger scales and characterization of EM structure–property–performance linkages at a full application scale.

Mathematics and Computing↗

Machine Learning Models of Plastic Flow Based on Representation Theory

We use machine learning (ML) to infer stress and plastic flow rules using data from representative polycrystalline simulations. In particular, we use so-called deep (multilayer) neural networks (NN) to represent the two response functions. The ML process does not choose appropriate inputs or outputs, rather it is trained on selected inputs and output. Likewise, its discrimination of features is crucially connected to the chosen input-output map. Furthermore, we draw upon classical constitutive modeling to select inputs and enforce well-accepted symmetries and other properties. With these developments, we enable rapid model building in real-time with experiments, and guide data collection and feature discovery.

42 ENGINEERING↗

Forecasting Dark Matter Subhalo Constraints from Stellar Streams using Implicit Likelihood Inference

The evidence for dark matter (DM) remains compelling, although attempts to understand its particle nature remain inconclusive. One promising method to study DM is detecting DM subhalos through their gravitational interactions with stellar streams. In this study, we apply Neural Posterior Estimation (NPE) to constrain subhalo interaction parameters, including mass, scale radius, velocity, and encounter geometry, from stellar stream kinematics. We generate particle spray simulations based on the Lagrange Cloud stripping technique, focusing on the ATLAS-Aliqa Uma stream as a test case. We train multiple NPE models across multiple observational scenarios, quantifying how kinematic completeness affects inference and forecasting constraints from upcoming surveys including LSST, 4MOST, and 10-year Gaia data. Our results demonstrate that NPE can produce accurate and well-calibrated posteriors. In the idealized case with full 6D coordinates, we achieve subhalo mass uncertainties of 15-20% for a $10^7 \, \mathrm{M_\odot}$ subhalo, with 5D coordinates (excluding radial velocities) achieving similar performance. Under realistic observational conditions, mass uncertainties range from 50% (present-day) to 20-40% (future scenarios), with comparable performance between the photometric-only LSST sample and a smaller sample that includes Gaia proper motions and 4MOST radial velocities. Most notably, we find that velocity bimodality emerges when phase space is poorly sampled, whether due to missing kinematic information or limited stellar tracers. Combining large photometric samples with targeted spectroscopic follow-up can effectively resolves this degeneracy. These results demonstrate the power of implicit likelihood inference for optimizing stellar stream observational strategies and forecasting DM subhalo constraints from upcoming surveys.

Nguyen, Tri [Northwestern U. (main); SkAI, Chicago↗

Baycal

Bayesian Model Calibration (BayCal) toolkit is a software plugin for Risk Analysis Virtual Environment (RAVEN) framework, arming at inversely quantifying the uncertainties associated with simulation model parameters based on available experiment data. BayCal seeks statistical inference of the uncertain input parameters that are consistent with the available measurement data or observed data. The unique feature of BayCal is the capability to be linked with RAVEN to build corresponding calibration workflows for complex multi-physics simulations. In addition to be able to use the machine learning capability of RAVEN to significantly reduce the computational cost of expensive simulation models, another distinctive feature of BayCal is the capability to deal with high-dimensional correlated model outputs, such as time series observations at multiple locations, via principal component analysis (PCA) technique.

Wang, Congjian↗

Simultaneous diagnosis of radial profiles and mix in NIF ignition-scale implosions via X-ray spectroscopy

In a National Ignition Facility implosion, hydrodynamic instabilities may cause the cold material from the imploding shell to be injected into the hot-spot (hot-spot mix), enhancing the radiative and conductive losses, which in turn may lead to a quenching of the ignition process. The bound-bound features of the spectrum emitted by high-Z ablator dopants that get mixed into the hot-spot have been previously used to infer the total amount of mixed mass; however, the typical errorbars are larger than the maximum tolerable mix. We present in this paper an improved 2D model for mix spectroscopy which can be used to retrieve information on both the amount of mixed mass and the full imploded plasma profile. By performing radiation transfer and simultaneously fitting all of the features exhibited by the spectra, we are able to constrain self-consistently the effect of the opacity of the external layers of the target on the emission, thus improving the accuracy of the inferred mixed mass. The model's predictive capabilities are first validated by fitting simulated spectra arising from fully characterized hydrodynamic simulations, and then, the model is applied to previously published experimental results, providing values of mix mass in agreement with previous estimates. Finally, we show that the new self consistent procedure leads to better constrained estimates of mix and also provides insight into the sensitivity of the hot-spot spectroscopy to the spatial properties of the imploded capsule, such as the in-flight aspect ratio of the cold fuel surrounding the hotspot.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Joint Bayesian analysis of large angular scale CMB temperature anomalies

Cosmic microwave background measurements show an agreement with the concordance cosmology model except for a few notable anomalies: Power Suppression , the lack of large scale power in the temperature data compared to what is expected in the concordance model, and Cosmic Hemispherical Asymmetry , a dipolar breakdown of statistical isotropy. An expansion of the CMB covariance in Bipolar Spherical Harmonics naturally parametrizes both these large-scale anomalies, allowing us to perform an exhaustive, fully Bayesian joint analysis of the power spectrum and violations of statistical isotropy up to the dipole level. Our analysis sheds light on the scale dependence of the Cosmic Hemispherical Asymmetry. Assuming a scale-dependent dipole modulation model with a two-parameter power law form, we explore the posterior pdf of amplitude A($l$ = 16) and the power law index α and find the maximum a posteriori values A * (l = 16) = 0.064 ± 0.022 and $α_*$ = -0.92 ± 0.22. The maximum a posteriori direction associated with the Cosmic Hemispherical Asymmetry is ($l,b$) = (247.8 o , -19.6 o ) in Galactic coordinates, consistent with previous analyses. We evaluate the Bayes factor BSI-DM to compare the Cosmic Hemispherical Asymmetry model with the isotropic model. The data prefer but do not substantially favor the anisotropic model (BSI-DM=0.4). We consider several priors and find that this evidence ratio is robust to prior choice. The large-scale power suppression does not soften when jointly inferring both the isotropic power spectrum and the parameters of the asymmetric model, indicating no evidence that these anomalies are coupled.

79 ASTRONOMY AND ASTROPHYSICS↗

Dynamic logistic regression and variable selection: Forecasting and contextualizing civil unrest

Civil unrest can range from peaceful protest to violent furor, and researchers are working to monitor, forecast, and assess such events to allocate resources better. Twitter has become a real-time data source for forecasting civil unrest because millions of people use the platform as a social outlet. Additionally, daily word counts are used as model features, and predictive terms contextualize the reasons for the protest. To forecast civil unrest and infer the reasons for the protest, we consider the problem of Bayesian variable selection for the dynamic logistic regression model and propose using penalized credible regions to select parameters of the updated state vector. This method avoids the need for shrinkage priors, is scalable to high-dimensional dynamic data, and allows the importance of variables to vary in time as new information becomes available. A substantial improvement in both precision and F1-score using this approach is demonstrated through simulation. Finally, we apply the proposed model fitting and variable selection methodology to the problem of forecasting civil unrest in Latin America. Our dynamic logistic regression approach shows improved accuracy compared to the static approach currently used in event prediction and feature selection.

97 MATHEMATICS AND COMPUTING↗

Milky Way Satellite Census -- II. Galaxy-Halo Connection Constraints Including the Impact of the Large Magellanic Cloud

The population of Milky Way (MW) satellites contains the faintest known galaxies, and thus provides essential insight into galaxy formation and dark matter microphysics. Here, we combine a model of the galaxy--halo connection with newly derived observational selection functions based on searches for satellites in photometric surveys over nearly the entire high-Galactic-latitude sky. In particular, we use cosmological zoom-in simulations of MW-like halos that include realistic Large Magellanic Cloud (LMC) analogs to fit the position-dependent MW satellite luminosity function. We report decisive evidence for the statistical impact of the LMC on the MW satellite population due to an estimated $6.5\pm 1.5$ observed LMC-associated satellites, consistent with the number of LMC satellites inferred from $\textit{Gaia}$ proper motion measurements, confirming the predictions of cold dark matter models for the existence of satellites within satellite halos. Moreover, we infer that the LMC fell into the MW within the last $2\ \rm{Gyr}$ at high confidence. Based on our detailed full-sky modeling, we find that the faintest observed satellites inhabit halos with peak virial masses below $2.2\times 10^{8}\ M_{\rm{\odot}}$ at $95\%$ confidence, and we place the first robust constraints on the fraction of halos that host galaxies in this regime. We predict that the faintest detectable satellites occupy halos with peak virial masses above $10^{6}\ M_{\rm{\odot}}$, highlighting the potential for powerful galaxy formation and dark matter constraints from future dwarf galaxy searches.

79 ASTRONOMY AND ASTROPHYSICS↗

Impact of petrophysical uncertainty on Bayesian hydrogeophysical inversion and model selection

Quantitative hydrogeophysical studies rely heavily on petrophysical relationships that link geophysical properties to hydrogeological properties and state variables. Coupled inversion studies are frequently based on the questionable assumption that these relationships are perfect (i.e., no scatter). Using synthetic examples and crosshole ground-penetrating radar (GPR) data from the South Oyster Bacterial Transport Site in Virginia, USA, we investigate the impact of spatially-correlated petrophysical uncertainty on inferred posterior porosity and hydraulic conductivity distributions and on Bayes factors used in Bayesian model selection. Our study shows that accounting for petrophysical uncertainty in the inversion (I) decreases bias of the inferred variance of hydrogeological subsurface properties, (II) provides more realistic uncertainty assessment and (III) reduces the overconfidence in the ability of geophysical data to falsify conceptual hydrogeological models.

58 GEOSCIENCES↗

Multi‐hypothesis comparison of Farquhar and Collatz photosynthesis models reveals the unexpected influence of empirical assumptions at leaf and global scales

Abstract Mechanistic photosynthesis models are at the heart of terrestrial biosphere models (TBMs) simulating the daily, monthly, annual and decadal rhythms of carbon assimilation ( A ). These models are founded on robust mathematical hypotheses that describe how A responds to changes in light and atmospheric CO 2 concentration. Two predominant photosynthesis models are in common usage: Farquhar (FvCB) and Collatz (CBGB). However, a detailed quantitative comparison of these two models has never been undertaken. In this study, we unify the FvCB and CBGB models to a common parameter set and use novel multi‐hypothesis methods (that account for both hypothesis and parameter variability) for process‐level sensitivity analysis. These models represent three key biological processes: carboxylation, electron transport, triose phosphate use (TPU) and an additional model process: limiting‐rate selection. Each of the four processes comprises 1–3 alternative hypotheses giving 12 possible individual models with a total of 14 parameters. To broaden inference, TBM simulations were run and novel, high‐resolution photosynthesis measurements were made. We show that parameters associated with carboxylation are the most influential parameters but also reveal the surprising and marked dominance of the limiting‐rate selection process (accounting for 57% of the variation in A vs. 22% for carboxylation). The limiting‐rate selection assumption proposed by CBGB smooths the transition between limiting rates and always reduces A below the minimum of all potentially limiting rates, by up to 25%, effectively imposing a fourth limitation on A . Evaluation of the CBGB smoothing function in three TBMs demonstrated a reduction in global A by 4%–10%, equivalent to 50%–160% of current annual fossil fuel emissions. This analysis reveals a surprising and previously unquantified influence of a process that has been integral to many TBMs for decades, highlighting the value of multi‐hypothesis methods.

54 ENVIRONMENTAL SCIENCES↗

Ice Nucleating Particles, Aerosols and Clouds over the Higher Latitude Southern Ocean

Improved understanding of aerosol-cloud-precipitation interactions is a key goal of the DOE-ASR program. Ice Nucleating Particles (INPs) are rare aerosol particles that trigger ice formation in clouds. By doing so they may initiate precipitation in clouds, and this impacts cloud lifetime. Nowhere may the role of INPs be manifested so dramatically as in the Southern Ocean (SO). This vast region was the focus of the DOE-ARM MARCUS (Measurement of Aerosols, Radiation and CloUds over the Southern Oceans) and MICRE (Macquarie Island Cloud and Radiation Experiment) campaigns. In the SO, a deficit of INPs has been a leading hypothesis for why liquid clouds persist, which underlays the bias of global climate models in predicting excess shortwave radiation reaching the ocean surface between 55°S and Antarctica. In this study, we used the 6-month acquisition of INP data on four MARCUS cruises from Hobart, Tasmania, to Antarctica, and the single-location full annual INP cycle observed in MICRE, to gain a holistic, quantitative picture of the INP number concentrations as a function of temperature and their variability over SO latitudes from 43°S to Antarctica. Augmenting initial processing of filter collections of aerosol particles during MARCUS and MICRE, additional freezing studies following thermal treatment (to removes INPs associated with microbes) and peroxide treatment (to remove all organics), as well as ionic chemistry and total organic carbon analysis, were applied to samples to improve time and space resolution, and differentiate marine from terrestrial contributions. We also performed Next Generation Sequencing to determine that bacterial composition can tag periods of enhanced and degraded marine organic INPs and the relative lack of occurrence of land-sourced bio-particles. Through these and other analyses, INP data were categorized into representative source types (e.g., marine versus terrestrially-influenced). The products of our analyses are being used to parameterize ice formation for use in numerical modeling studies to determine if, and when, ocean-derived INPs control the microphysical composition and radiative balance of SO clouds in a manner not presently captured by global models. The specific objectives and approaches proposed were largely accomplished, including: 1) Additional and value-added chemical and biological analyses of archived samples of aerosols were completed to categorize and quantify INP types and concentrations over the SO. 2) Analyses were completed to place INP data in meteorological and aerosol context. 3) Sea spray aerosol INPs were demonstrated to dominate the SO region in all but episodic events. 4) New parameterizations for INP sources relevant to the SO were constructed. 5) Collaborative modeling studies demonstrating the crucial role of INPs in determining cloud radiative and precipitation properties were conducted with MARCUS and MICRE partners. Collaborative publication submissions in this regard ensued within the two year study. This work will advance the science of interactions of aerosols, clouds and precipitation, to improve their representation in regional and global climate models. Results have informed representation of primary ice crystal nucleation, give inference to the role of secondary processes, and improve investigations of aerosol influences on clouds via ice nucleation over SO high latitudes and similar vast ocean regions. Application will improve representation of such interactions for clouds in both regional and global climate models. The data base and analyses methods applied will continue to serve research studies in the future.

58 GEOSCIENCES↗

The Dark World: A Tale of WASP-43b in Reflected Light with HST WFC3/UVIS

Optical reflected light eclipse observations provide a direct probe of exoplanet scattering properties, such as from aerosols. We present here the photometric reflected light observations of WASP-43b using the Hubble Space Telescope (HST) WFC3/UVIS instrument with the F350LP filter (346–822 nm) encompassing the entire optical band. This is the first reflected light photometric eclipse using UVIS in scanning mode; as such, we further detail our scanning extraction and analysis pipeline Arctor. Our HST WFC3/UVIS eclipse light curve for WASP-43b derived a 3σ upper limit of 67 ppm on the eclipse depth, which implies that WASP-43b has a very dark dayside atmosphere. With our atmospheric modeling campaign, we compared our reflected light constraints with predictions from global circulation and cloud models benchmarked with HST and Spitzer observations of WASP-43b. We infer that we do not detect clouds on the dayside within the pressure levels probed by HST WFC3/UVIS with the F350LP filter (P > 1 bar). This is consistent with the general circulation model predictions based on previous WASP-43b observations. Dayside emission spectroscopy results from WASP-43b with HST and Spitzer observations are likely to not be significantly affected by contributions from cloud particles.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Machine learning coarse grained models for water

An accurate and computationally efficient molecular level description of mesoscopic behavior of ice-water systems remains a major challenge. Here, we introduce a set of machine-learned coarse-grained (CG) models (ML-BOP, ML-BOP dih , and ML-mW) that accurately describe the structure and thermodynamic anomalies of both water and ice at mesoscopic scales, all at two orders of magnitude cheaper computational cost than existing atomistic models. In a significant departure from conventional force-field fitting, we use a multilevel evolutionary strategy that trains CG models against not just energetics from first-principles and experiments but also temperature-dependent properties inferred from on-the-fly molecular dynamics (~ 10’s of milliseconds of overall trajectories). Our ML BOP models predict both the correct experimental melting point of ice and the temperature of maximum density of liquid water that remained elusive to-date. Our ML workflow navigates efficiently through the high-dimensional parameter space to even improve upon existing high-quality CG models (e.g. mW model).

42 ENGINEERING↗

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↗

A Correct Flat Field Model For DARHT

The forward modeling approach has been utilized to obtain the density distribution during dynamic radiographic experiments using the Bayesian Inference Engine (BIE) at the Dual Axis Radiographic Hydrodynamic Test Facility (DARHT) for over two decades. The forward modeling approach attempts to incorporate all of the known physics to obtain the transmission, from an unknown density field. The density field is then optimized in an attempt to match the flat-fielded radiographic image using a Chi-squared comparison. A depiction of the forward modeling approach implemented in the BIE is given.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Dark Matter Constraints from a Unified Analysis of Strong Gravitational Lenses and Milky Way Satellite Galaxies

Joint analyses of small-scale cosmological structure probes have been relatively unexplored yet promise to advance measurements of microphysical dark matter properties using heterogeneous data. Here, we present a multidimensional analysis of dark matter substructure using strong gravitational lenses and the Milky Way (MW) satellite galaxy population, accounting for degeneracies in model predictions and using covariances in the constraining power of these individual probes for the first time. We simultaneously infer the projected subhalo number density and the half-mode mass describing the suppression of the subhalo mass function in thermal relic warm dark matter (WDM), $M_{\mathrm{hm}}$, using the semianalytic model \texttt{Galacticus} to connect the subhalo population inferred from MW satellite observations to the strong lensing host halo mass and redshift regime. Combining MW satellite and strong lensing posteriors in this parameter space yields $M_{\mathrm{hm}}<10^{7.0}\ M_{\mathrm{\odot}}$ (WDM particle mass $m_{\mathrm{WDM}}>9.7\ \mathrm{keV}$) at $95\%$ confidence and disfavors $M_{\mathrm{hm}}=10^{7.4}\ M_{\mathrm{\odot}}$ ($m_{\mathrm{WDM}}=7.4\ \mathrm{keV}$) with a 20:1 marginal likelihood ratio, improving limits on $m_{\mathrm{WDM}}$ set by the two methods independently by $\sim 30\%$. These results are marginalized over the line-of-sight contribution to the strong lensing signal, the mass of the MW host halo, and the efficiency of subhalo disruption due to baryons and are robust to differences in the disruption efficiency between the MW and strong lensing regimes at the $\sim 10\%$ level. This work paves the way for unified analyses of next-generation small-scale structure measurements covering a wide range of scales and redshifts.

79 ASTRONOMY AND ASTROPHYSICS↗

The Intrinsic Shapes of Low Surface Brightness Galaxies (LSBGs): A Discriminant of LSBG Galaxy Formation Mechanisms

We use the low surface brightness galaxy (LSBG) samples created from the Hyper Suprime-Cam Subaru Strategic Program (781 galaxies), the Dark Energy Survey (20977 galaxies), and the Legacy Survey (selected via H ι detection in the Arecibo Legacy Fast ALFA Survey, 188 galaxies) to infer the intrinsic shape distribution of the LSBG population. To take into account the effect of the surface brightness cuts employed when constructing LSBG samples, we simultaneously model both the projected ellipticity and the apparent surface brightness in our shape inference. We find that the LSBG samples are well characterized by oblate spheroids, with no significant difference between red and blue LSBGs. This inferred shape distribution is in good agreement with similar inferences made for ultra-diffuse cluster galaxy samples, indicating that environment does not play a key role in determining the intrinsic shape of LSBGs. In this study, we also find some evidence that LSBGs are more thickened than similarly massive high surface brightness dwarfs. We compare our results to intrinsic shape measures from contemporary cosmological simulations, and find that the observed LSBG intrinsic shapes place considerable constraints on the formation path of such galaxies. In particular, LSBG production via the migration of star formation to large radii produces intrinsic shapes in good agreement with our observational findings.

79 ASTRONOMY AND ASTROPHYSICS↗