Search NASA⌕ Search

SEARCH · Search NASA

Results for “inference methods”

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

Operator inference for non-intrusive model reduction of systems with non-polynomial nonlinear terms

Here in this work we present a non-intrusive model reduction method to learn low-dimensional models of dynamical systems with non-polynomial nonlinear terms that are spatially local and that are given in analytic form. In contrast to state-of-the-art model reduction methods that are intrusive and thus require full knowledge of the governing equations and the operators of a full model of the discretized dynamical system, the proposed approach requires only the non-polynomial terms in analytic form and learns the rest of the dynamics from snapshots computed with a potentially black-box full-model solver. The proposed method learns operators for the linear and polynomially nonlinear dynamics via a least-squares problem, where the given non-polynomial terms are incorporated on the right-hand side. The least-squares problem is linear and thus can be solved efficiently in practice. The proposed method is demonstrated on three problems governed by partial differential equations, namely the diffusion–reaction Chafee–Infante model, a tubular reactor model for reactive flows, and a batch-chromatography model that describes a chemical separation process. The numerical results provide evidence that the proposed approach learns reduced models that achieve comparable accuracy as models constructed with state-of-the-art intrusive model reduction methods that require full knowledge of the governing equations.

42 ENGINEERING↗
DOE OSTI2025

Application of Indirect Quantification of 133mXe to Calibration of HPGe Detector

The INL Noble Gas Laboratory provides intercomparison samples for the noble gas analysis laboratories as part of the CTBTO PrepCom IMS. Xe-133m is one of the four relevant radionuclides in nuclear explosion monitoring. Without commercially available Xe-133m calibration standards laboratories must create and improve calibration methods. Improvements in calibration methods at the INL NGL benefit the CTBTO PrepCom through better certified values for Xe-133m intercomparison samples. Calibration of High Purity Germanium detectors for Xe-133m quantification is complicated by the coexistence of Xe-133 in samples under analysis. Xe-133 is typically produced in larger quantities, has higher gamma emission probabilities, and its gammas are detected more efficiently than Xe-133m. Xe-133m activity of samples can be indirectly inferred through the 133:133m activity ratio of a batch of material, and the Xe-133 counts in the assay of a small aliquot of the same material. This indirect quantification method can be leveraged to perform detector calibrations for quantification of Xe-133m. Calibrations can be performed by inferring the Xe-133m to certify the sample, and direct counting to determine detector efficiency. A comparison of method results will be shown.

133mXe↗
DOE OSTI2019

Estimating redshift distributions using hierarchical logistic Gaussian processes

ABSTRACT This work uses hierarchical logistic Gaussian processes to infer true redshift distributions of samples of galaxies, through their cross-correlations with spatially overlapping spectroscopic samples. We demonstrate that this method can accurately estimate these redshift distributions in a fully Bayesian manner jointly with galaxy-dark matter bias models. We forecast how systematic biases in the redshift-dependent galaxy-dark matter bias model affect redshift inference. Using published galaxy-dark matter bias measurements from the Illustris simulation, we compare these systematic biases with the statistical error budget from a forecasted weak gravitational lensing measurement. If the redshift-dependent galaxy-dark matter bias model is mis-specified, redshift inference can be biased. This can propagate into relative biases in the weak lensing convergence power spectrum on the 10–30 per cent level. We, therefore, showcase a methodology to detect these sources of error using Bayesian model selection techniques. Furthermore, we discuss the improvements that can be gained from incorporating prior information from Bayesian template fitting into the model, both in redshift prediction accuracy and in the detection of systematic modelling biases.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1994

On the rms anisotropy at 7 deg and 10 deg observed in the COBE-DMR two year sky maps

The frequency-independent rms temperature fluctuations determined from the Cosmic Background Explorer-Differential Microwave Radiometer (COBE-DMR) two-year sky maps are used to infer the parameter Q(sub rms-PS), which characterizes the normalization of power-law models of primordial cosmological temperature anisotropy, for a forced fit to a scale-invariant Harrison-Zel'dovich (n = 1) spectral model. Using a joint analysis of the 7 deg and 10 deg 'cross'-rms derived from both the 53 and 90 GHz sky maps, we find Q(sub rms-PS) = 17.0(sub -2.1 sup +2.5) micro Kelvin when the low quadrupole is included, and Q(sub rms-PS) = 19.4(sub -2.1 sup +2.3) micro Kelvin excluding the quadrupole. These results are consistent with the n = 1 fits from more sensitive methods. The effect of the low quadrupole derived from the COBE-DMR data on the inferred Q(sub rms-PS) normalization is investigated. A bias to lower Q(sub rms-PS) is found when the quadrupole is included. The higher normalization for a forced n = 1 fit is then favored by the cross-rms technique.

Banday, A. J.↗
NASA NTRS1992

On the Retrieval and Analysis of Multilevel Clouds

An accurate satellite retrieval of cloud properties depends upon the detection and analysis of multilayered, overlapping cloud systems that surface observations show to be common. Multiple cloud layers are often found, for instance, in frontal situations, where cirrus overlays boundary layer convective cloud or low-to mid-level stratus cloud. Surface observers (Hahan et al., 1982) indicate that over ocean in the Northern Hemisphere between 30 deg. N and 60 deg. N, 51 percent of observations are of multilevel clouds. A satellite analysis by Coakley (1983) over the Pacific Ocean finds that more than 50 percent of 500 (250 sq km) frames exhibit evidence of multilayered cloud systems. The questions addressed in this study are the following: What error is introduced when inferring the cloud pressure from a Field-Of-View (FOV) that contains some arbitrary amount of transparent cloud overlaying a lower-level black cloud, such as stratus, by making the assumption that there is only a single cloud layer in the FOV, and what may be done to improve the cloud retrieval? The CO2 slicing methods (e.g. McCleese and Wilson, 1976; Smith and Platt, 1978; Chahine, 1974) have been shown to provide accurate means of inferring cirrus cloud altitude from passive infrared radiance measurements. The CO2 techniques have been applied to radiometric data from several instruments, notably the High Resolution Infrared Radiometric Sounder (HIRS/2, hereafter referred to as HIRS), the VISSR Atmospheric Sounder (VAS) (e.g., Menzel et al., 1983; Wylie and Menzel, 1989), and most recently to the High Resolution Interferometer Sounder (HIS) (Smith and Frey, 1990). The methods take advantage of the fact that infrared CO2 sounding channels spaced closely in wavenumber each have varying opacity to CO2, thereby causing each channel to be sensitive to a different level in the atmosphere. The techniques have been shown to be effective for single-layered, nonblack, mid- to high-level clouds such as cirrus, but are generally applied operationally to any given cloud occurrence. The CO2 slicing algorithms are most accurate for clouds than occur in a single, well-defined layer, or for multi-layered cloud cases in which the uppermost cloud layer is nearly black. Significant cloud height retrieval errors may ensue if the HIRS Field-Of-View (FOV) is cotaminated with low cloud. McCleese and Wilson (1976) have shown that the retrieved cloud height for the case of multiple cloud layers is a weighted average of the cloud heights actually present. The weight is approximately proportional to the product of the cloud heigt and the effective cloud amount. The effect of their result is that the uppermost cloud layer dominates the cloud pressure retrieval. Beyond stating that the higher cloud dominates the cloud pressure retrieval, there is no quantitative information to provide a way of estimating the errors in cloud pressure retrieval one should expect for certain common multilevel cloud situations or any suggestions on how to reduce the errors. In this paper we estimate the magnitude of the errors and use a simple algorithm to reduce the errors in optically thin cloud height retrival.

Baum, Bryan A.↗
DOE OSTI2022

Dissociating the phononic, magnetic and electronic contributions to thermal conductivity: a computational study in alpha-iron

Computational tools to study thermodynamic properties of magnetic materials have, until recently, been limited to phenomenological modeling or to small domain sizes limiting our mechanistic understanding of thermal transport in ferromagnets. Herein, we study the interplay of phonon and magnetic spin contributions to the thermal conductivity in a-iron utilizing non-equilibrium molecular dynamics simulations. It was observed that the magnetic spin contribution to the total thermal conductivity exceeds lattice transport for temperatures up to two-thirds of the Curie temperature after which only strongly coupled magnon-phonon modes become active heat carriers. Characterizations of the phonon and magnon spectra give a detailed insight into the coupling between these heat carriers, and the temperature sensitivity of these coupled systems. Comparisons to both experiments and ab initio data support our inferred electronic thermal conductivity, supporting the coupled molecular dynamics/spin dynamics framework as a viable method to extend the predictive capability for magnetic material properties.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Towards stacking fault energy engineering in FCC high entropy alloys

Stacking Fault Energy (SFE) is an intrinsic alloy property that governs much of the plastic deformation mechanisms observed in fcc alloys. While SFE has been recognized for many years as a key intrinsic mechanical property, its inference via experimental observations or prediction using, for example, computationally intensive first-principles methods is challenging. This difficulty precludes the explicit use of SFE as an alloy design parameter. In this work, we combine DFT calculations (with necessary configurational averaging), machine-learning (ML) and physics-based models to predict the SFE in the fcc CoCrFeMnNiV-Al high-entropy alloy space. The best-performing ML model is capable of accurately predicting the SFE of arbitrary compositions within this 7-element system. Finally, this efficient model along with a recently developed model to estimate intrinsic strength of fcc HEAs is used to explore the strength–SFE Pareto front, predicting new-candidate alloys with particularly interesting mechanical behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A Robust Estimate of Continental-Scale Terrestrial Carbon Sinks Using GOSAT XCO 2 Retrievals

Satellite XCO 2 retrievals could improve the estimates of surface carbon fluxes, but it remains unknown on what scales these estimates are robust. Here, we use the time-dependent Bayesian synthesis top-down method and prior net ecosystem exchanges (NEEs) from 12 terrestrial biosphere models (TBMs) to infer the monthly carbon fluxes of 51 land regions with constraints by GOSAT XCO 2 retrievals. We find that the uncertainty (standard deviation of 12 TBMs) reduction rates (uncertainty reduction rate (URR)) decrease significantly at decreasing spatial scales. On the continental-scale, the mean URR is about 57%, and the annual and seasonal cycle estimates of NEE are rather robust. The evaluation shows that the posterior CO 2 concentrations are significantly improved at the continental scale. Our study suggests that the GOSAT XCO 2 can only promise a robust continental-scale NEE estimate and improving the XCO 2 accuracy is an effective way to achieve robust estimates on smaller scales under current spatial coverage.

58 GEOSCIENCES↗
DOE OSTI2018

Bayesian uncertainty calculation in neural network inference of ion and electron temperature profiles at W7-X

We make use of a Bayesian description of the neural network (NN) training for the calculation of the uncertainties in the NN prediction. Having uncertainties on the NN prediction allows having a quantitative measure for trusting the NN outcome and comparing it with other methods. Within the Bayesian framework, the uncertainties can be calculated under different approximations. The NN has been trained with the purpose of inferring ion and electron temperature profile from measurements of a X-ray imaging diagnostic at W7-X. The NN has been trained in such a way that it constitutes an approximation of a full Bayesian model of the diagnostic, implemented within the Minerva framework. Finally, the network has been evaluated using measured data and the uncertainties calculated under different approximations have been compared with each other, finding that neglecting the noise on the NN input can lead to an underestimation of the error bar magnitude in the range of 10%–30%.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Determining the electrical conductivity of metals using the 2 MA Thor pulsed power driver

We present the development of a pulsed power experimental technique to infer the electrical conductivity of metals from ambient to high energy density conditions. The method is implemented on Thor, a moderate scale (1–2 MA) pulsed power driver. The electrical conductivity of copper at elevated temperature (>4000 K) and pressure (>10 GPa) is determined, and a new tabular material model is developed, guided by density functional theory, which preserves agreement with existing experimental data. Minor modifications (<10%) are found to be necessary to the previous Lee–More–Desjarlais model isotherms in the vicinity of the melt transition in order to account for observed discrepancies with the new experimental data. An analytical model for magnetic direct drive flyer acceleration and Joule heating induced vaporization based on the Tsiolkovsky “rocket equation” is presented to assess sensitivity of the method to minor changes in electrical conductivity.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Structural inference of networked dynamical systems with universal differential equations

Networked dynamical systems are common throughout science in engineering; e.g., biological networks, reaction networks, power systems, and the like. For many such systems, nonlinearity drives populations of identical (or near-identical) units to exhibit a wide range of nontrivial behaviors, such as the emergence of coherent structures (e.g., waves and patterns) or otherwise notable dynamics (e.g., synchrony and chaos). Here, we seek to infer (i) the intrinsic physics of a base unit of a population, (ii) the underlying graphical structure shared between units, and (iii) the coupling physics of a given networked dynamical system given observations of nodal states. These tasks are formulated around the notion of the Universal Differential Equation, whereby unknown dynamical systems can be approximated with neural networks, mathematical terms known a priori (albeit with unknown parameterizations), or combinations of the two. We demonstrate the value of these inference tasks by investigating not only future state predictions but also the inference of system behavior on varied network topologies. The effectiveness and utility of these methods are shown with their application to canonical networked nonlinear coupled oscillators.

97 MATHEMATICS AND COMPUTING↗
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 OSTI2021

Anomaly Attribution with Likelihood Compensation

This paper addresses the task of explaining anomalous predictions of a black-box regression model. When using a black-box model, such as one to predict building energy consumption from many sensor measurements, we often have a situation where some observed samples may significantly deviate from their prediction. It may be due to a sub-optimal black-box model, or simply because those samples are outliers. In either case, one would ideally want to compute a responsibility score indicative of the extent to which an input variable is responsible for the anomalous output. In this work, we formalize this task as a statistical inverse problem: Given model deviation from the expected value, infer the responsibility score of each of the input variables. We propose a new method called likelihood compensation (LC), which is founded on the likelihood principle and computes a correction to each input variable. To the best of our knowledge, this is the first principled framework that computes a responsibility score for real valued anomalous model deviations. We apply our approach to a real-world building energy prediction task and confirm its utility based on expert feedback.

Idé, Tsuyoshi↗
DOE OSTI2023

What is ab initio in nuclear theory?

Ab initio has been used as a label in nuclear theory for over two decades. Its meaning has evolved and broadened over the years. We present our interpretation, briefly review its historical use, and discuss its present-day relation to theoretical uncertainty quantification.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Evidence for Cosmic Acceleration Is Robust to Observed Correlations between Type Ia Supernova Luminosity and Stellar Age

Type Ia supernovae (SNe Ia) are powerful standardizable candles for constraining cosmological models and provided the first evidence of the accelerated expansion of the universe. Their precision derives from empirical correlations, now measured from >1000 SNe Ia, between their luminosities, light-curve shapes, colors, and most recently with the stellar mass of their host galaxy. As mass correlates with other galaxy properties, alternative parameters have been investigated to improve SN Ia standardization though none have been shown to significantly alter the determination of cosmological parameters. We re-examine a recent claim, based on 34 SN Ia in nearby passive host galaxies, of a 0.05 mag Gyr -1 dependence of standardized SN Ia luminosity on host age, which, if extrapolated to higher redshifts, would be a bias up to 0.25 mag, challenging the inference of dark energy. We reanalyze this sample of hosts using both the original method and a Bayesian hierarchical model and find after a fuller accounting of the uncertainties the significance of a dependence on age to be ≤2σ and ~1σ after the removal of a single poorly sampled SN Ia. Here, to test the claim that a trend seen in old stellar populations can be applied to younger ages, we extend our analysis to a larger sample that includes young hosts. We find the residual dependence of host age (after all standardization typically employed for cosmological measurements) to be consistent with zero for 254 SNe Ia from the Pantheon sample, ruling out the large but low significance trend seen in passive hosts.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1972

Ogo 5 magnetic-field data near the earth's bow shock - A correlation with theory.

Magnetic-field data obtained in the earth's bow-shock region with a high-resolution triaxial fluxgate magnetometer aboard the Ogo 5 satellite have been correlated with a theory of Tidman and Northrop (1968). These authors have shown that either of two hypotheses about the nature of low-frequency magnetic waves could be invoked to explain previous observations. We have observed exponentially decaying upstream waves that are consistent with only one of these hypotheses. This observation allows use of the theory to infer the local shock velocity and frequency of driving currents within the shock. This method of finding the shock velocity is less sensitive to errors in the plasma parameters than is the method based on the Rankine-Hugoniot relations.

Guha, J. K.↗
NASA NTRS1983

Endurance test of a 30-CM-diameter engineering model ion thruster. Task 12: Investigation of thin-film erosion monitors for ion thrusters

An investigation of short term measurement techniques for predicting the wearout of ion thrusters resulting from sputter erosion damage is described. The previously established laminar thin film techniques to provide high precision erosion rate data. However, the erosion rates obtained using this technique are generally substantially higher than those obtained during long term endurance tests (by virtue of the as deposited nature of the thin films), so that the results must be interpreted in a relative sense. Absolute measurements can be performed using a new masked substrate arrangement which was developed during this study. This new technique provides a means for estimating the lifetimes of critical discharge chamber components based on direct measurements of sputter erosion depths obtained during short duration (10 hour) tests. The method enables the effects on lifetime of thruster design and operating parameters to be inferred without the investment of the time and capital required to conduct long term (1000 hour) endurance tests. Results obtained using the direct measurement technique are shown to agree with sputter erosion depths calculated for the plasma conditions of the test and also with lifetest results. The direct measurement approach is shown to be applicable to both mercury and argon discharge plasma environments and should be useful in estimating the lifetimes of inert gas and extended performance mercury ion thrusters presently under development.

Beattie, J. R.↗
NASA NTRS1985

Stress rate and proof-testing of silicon wafers

Fracture mechanics test methods were applied to evaluate the proof-test characteristics of single-crystal silicon wafers. The results indicate that the strength distribution of silicon wafers is truncated by proof-testing. No subcritical crack growth occurred during proof-loading, as inferred from the lack of a stress-rate effect on strength. Mechanical proof-testing appears to be an effective method for eliminating weak samples before cell processing.

Chen, C. P.↗