Search NASA⌕ Search

SEARCH · Search NASA

Results for “Bayesian Analysis”

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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
DOE OSTI2010

LOCAL NON-GAUSSIANITY IN THE COSMIC MICROWAVE BACKGROUND THE BAYESIAN WAY

We introduce an exact Bayesian approach to search for non-Gaussianity of local type in cosmic microwave background (CMB) radiation data. Using simulated CMB temperature maps, the newly developed technique is compared against the conventional frequentist bispectrum estimator. Starting from the joint probability distribution, we obtain analytic expressions for the conditional probabilities of the primordial perturbations given the data, and for the level of non-Gaussianity, f{sub NL}, given the data and the perturbations. We propose Hamiltonian Monte Carlo sampling as a means to derive realizations of the primordial fluctuations from which we in turn sample f{sub NL}. Although computationally expensive, this approach allows us to construct exactly the full target posterior probability distribution. When compared to the frequentist estimator, applying the Bayesian method to Gaussian CMB maps provides consistent results. For the analysis of non-Gaussian maps, however, the error bars on f{sub NL} do not show excess variance within the Bayesian framework. This finding is of particular relevance in the light of upcoming high-precision CMB measurements obtained by the Planck satellite mission.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Planck intermediate results: XVI. Profile likelihoods for cosmological parameters

We explore the 2013 Planck likelihood function with a high-precision multi-dimensional minimizer (Minuit). This allows a refinement of the CDM best-fit solution with respect to previously-released results, and the construction of frequentist confidence intervals using profile likelihoods. The agreement with the cosmological results from the Bayesian framework is excellent, demonstrating the robustness of the Planck results to the statistical methodology. We investigate the inclusion of neutrino masses, where more significant differences may appear due to the non-Gaussian nature of the posterior mass distribution. By applying the Feldman-Cousins prescription, we again obtain results very similar to those of the Bayesian methodology. However, the profile-likelihood analysis of the cosmic microwave background (CMB) combination ( Planck +WP+highL) reveals a minimum well within the unphysical negative-mass region. We show that inclusion of the Planck CMB-lensing information regularizes this issue, and provide a robust frequentist upper limit Σ$m_v$ ≤0:26 eV (95% confidence) from the CMB+lensing+BAO data combination.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

First observation of increased DT yield over prediction due to addition of hydrogen

Here, in a number of reported instances, implosions utilizing fuel mixtures have resulted in anomalously low fusion yields below those predicted by radiation-hydrodynamics simulations. Inter-species ion diffusion has been suggested as a possible cause of the observed yield degradation in fuel mixture implosions. An experimental platform utilizing hydro-equivalent deuterium–tritium (DT), deuterium–tritium–hydrogen (DTH), and deuterium-tritium-helium3 (DT 3 He) capsule implosions was developed to determine whether the inter-species ion diffusion theory may describe the resulting fuel mixture implosion behavior. The implosion experiments were performed at the Omega laser facility. X-ray images and shell areal density diagnostics results show that the hydro-equivalent three capsules (DT, DTH, and DT 3 He) have similar compression behavior. However, nuclear yield deviation was observed from the scaling determined using a fusion yield formula. In the DT 3 He mixture, a reduced yield of a factor of 0.65 ± 0.13 was observed, which is similar to a yield reduction observed in D 3 He mixture by Rygg et al. (i.e., Rygg effect). In contrast, in the DTH mixture, a factor of 1.17 ±0.15 yield increase was observed, which we named the inverse Rygg effect. The yield increase observed in the DTH mixture is consistent with the inter-species ion diffusion theory where lighter H diffuses away from the core and concentrated DT in the core produces higher yield. An inter-species ion diffusion model, the Zimmerman–Paquette–Kagan–Zhdanov model, implemented in a Lagrangian radiation-hydrodynamics fluid code, was also used to analyze the present data, without the need to assume hydrodynamic equivalence of the capsules, but it does not completely explain the DTH or DT 3 He capsules although its effects are in the correct direction. Simulation-based Bayesian inference was used in the latter analysis to quantify the uncertainty in the numerical simulations. The simulation-based analysis resulted in an inferred Rygg-effect yield decrease factor of 0.91 ± 0.02 for the DT 3 He mixture, and an inferred inverse-Rygg yield increase factor of 1.21 ± 0.04 for the DTH mixture, based on simulations ignoring ion diffusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Analysis of Repeatability and Reproducibility Studies With Ordinal Measurements

A Bayesian inferential approach with a noninformative prior is introduced to analyze ordinal repeatability and reproducibility (R&R) data using the De Mast–Van Wieringen model. This approach is extended with a weakly informative prior and random effects to allow for the consideration of a population of raters and prediction of a new rater. This random-effects approach is also shown to result in partial pooling of estimates across raters. In addition, match-probability-based measures to decompose ordinal R&R study data into contributions due to repeatability and due to reproducibility are defined. All extensions involving Bayesian inference (for fixed or random effects) and measures are illustrated on real and simulated ordinal R&R study data and are applicable in business and industry settings. This methodology can be implemented using the supplemental R package ordinalRR available from CRAN. Additional supplementary material for this article is available online.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Probabilistic modeling and prediction of out-of-plane unidirectional composite lamina properties

Computational simulation provides an efficient means to predict the behavior of customized hybrid material configurations using validated, physics-based models. One limitation to this approach is the quality and quantity of available data to characterize the many constituent input properties. Therefore, a systematic approach to identify the most influential parameters on the hybrid behavior and quantify the corresponding uncertainty in predictive capabilities is required. Here, an approach using Bayesian multimodel inference and imprecise global sensitivity analysis is presented to investigate the effects of sparse constituent data on the prediction of composite material properties. The methodology allows the identification, using quantified uncertainties, of the most influential constituent material parameters for specified homogenized properties. This sensitivity analysis further enables a dimension reduction when assessing the influence of uncertainties on material properties and can be used to inform testing programs of the constituent properties that require additional testing/data collection in order to minimize uncertainty in macro-scale composite properties. The methodology is specifically demonstrated on the prediction and sensitivity analysis of out-of-plane mechanical properties of a unidirectional lamina.

36 MATERIALS SCIENCE↗
DOE OSTI2013

Constraints on the tensor-to-scalar ratio for non-power-law models

Recent cosmological observations hint at a deviation from the simple power-law form of the primordial spectrum of curvature perturbations. In this paper we show that in the presence of a tensor component, a turn-over in the initial spectrum is preferred by current observations, and hence non-power-law models ought to be considered. For instance, for a power-law parameterisation with both a tensor component and running parameter, current data show a preference for a negative running at more than 2.5σ C.L. As a consequence of this deviation from a power-law, constraints on the tensor-to-scalar ratio r are slightly broader. We also present constraints on the inflationary parameters for a model-independent reconstruction and the Lasenby and Doran (LD) model. In particular, the constraints on the tensor-to-scalar ratio from the LD model are: r{sub LD} = 0.11±0.024. In addition to current data, we show expected constraints from Planck-like and CMB-Pol sensitivity experiments by using Markov-Chain-Monte-Carlo sampling chains. For all the models, we have included the Bayesian Evidence to perform a model selection analysis. The Bayes factor, using current observations, shows a strong preference for the LD model over the standard power-law parameterisation, and provides an insight into the accuracy of differentiating models through future surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Statistical tools for a better optical model

Background: Modern statistical tools provide the ability to compare the information content of observables and provide a path to explore which experiments would be most useful to give insight into and constrain theoretical models. Purpose: Here we study three such tools in the context of nuclear reactions with the goal of constraining the optical potential. Method: The three statistical tools examined are (i) the principal component analysis, (ii) the sensitivity analysis based on derivatives, and (iii) the Bayesian evidence. We first apply these tools to a toy-model case, comparing the form of the imaginary part of the optical potential. Then we consider two different reaction observables, elastic angular distributions and polarization data for reactions on 48 Ca and 208 Pb at two different beam energies. Results: For the toy-model case, we find significant discrimination power in the sensitivities and the Bayesian evidence, showing clearly that the volume imaginary term is more useful to describe scattering at higher energies. When comparing between elastic cross sections and polarization data using realistic optical models, sensitivity studies indicate that both observables are roughly equally sensitive but the variability of the optical model parameters is strongly angle dependent. The Bayesian evidence shows some variability between the two observables, but the Bayes factor obtained is not sufficient to discriminate between angular distributions and polarization. Conclusions: From the cases considered, we conclude that, in general, elastic scattering angular distributions have similar impact in constraining the optical potential parameters compared with the polarization data. The angular ranges for the optimum experimental constraints can vary significantly with the observable considered.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Landmark-Warped Emulators for Models with Misaligned Functional Response

Many computer models output functional data, and in some cases, these functional data have similar, but misaligned, shape characteristics. In this paper, we introduce a general approach for building emulators for computer models that output misaligned functional data when key values in the functional response (landmarks) can be easily identified. This approach has two main parts: modeling the aligned (using the landmarks) functional data, and modeling the functions that map the misaligned data to the aligned space (warping functions). As the warping functions are required to be monotonic, we give special attention to modeling monotonic functional response data. We discuss how our approach can be easily applied for a variety of typical emulators, such as Gaussian processes, Bayesian multivariate adaptive regression splines, and Bayesian additive regression trees, and how sensitivity analysis can be performed. We demonstrate our approach by building emulators for two applications: (1) a high-energy-density physics computer model used to simulate inertial confinement fusion ignition experiments, where model outputs are highly misaligned, and (2) a multiphysics continuum hydrocode used to simulate high-velocity impact experiments, where model outputs are only slightly misaligned. In case (1) traditional methods cannot be applied, while in (2) they can be applied, but the proposed method performs significantly better.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2013

Use of limited data to construct Bayesian networks for probabilistic risk assessment.

Probabilistic Risk Assessment (PRA) is a fundamental part of safety/quality assurance for nuclear power and nuclear weapons. Traditional PRA very effectively models complex hardware system risks using binary probabilistic models. However, traditional PRA models are not flexible enough to accommodate non-binary soft-causal factors, such as digital instrumentation&control, passive components, aging, common cause failure, and human errors. Bayesian Networks offer the opportunity to incorporate these risks into the PRA framework. This report describes the results of an early career LDRD project titled %E2%80%9CUse of Limited Data to Construct Bayesian Networks for Probabilistic Risk Assessment%E2%80%9D. The goal of the work was to establish the capability to develop Bayesian Networks from sparse data, and to demonstrate this capability by producing a data-informed Bayesian Network for use in Human Reliability Analysis (HRA) as part of nuclear power plant Probabilistic Risk Assessment (PRA). This report summarizes the research goal and major products of the research.

Groth, Katrina M.↗
DOE OSTI2013

UQTk version 2.0 user manual

The UQ Toolkit (UQTk) is a collection of libraries and tools for the quantification of uncertainty in numerical model predictions. Version 2.0 ffers intrusive and non-intrusive methods for propagating input uncertainties through computational models, tools for sensitivity analysis, methods for sparse surrogate construction, and Bayesian inference tools for inferring parameters from experimental data. This manual discusses the download and installation process for UQTk, provides pointers to the UQ methods used in the toolkit, and describes some of the examples provided with the toolkit.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

UQTk Version 3.0.3 User Manual

The UQ Toolkit (UQTk) is a collection of libraries and tools for the quantification of uncertainty in numerical model predictions. Version 3.0.3 offers intrusive and non-intrusive methods for propagating input uncertainties through computational models, tools for sen- sitivity analysis, methods for sparse surrogate construction, and Bayesian inference tools for inferring parameters from experimental data. This manual discusses the download and installation process for UQTk, provides pointers to the UQ methods used in the toolkit, and describes some of the examples provided with the toolkit.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗
DOE OSTI2016

UQTk (V. 3.0) User Manual

The UQ Toolkit (UQTk) is a collection of libraries and tools for the quantification of uncertainty in numerical model predictions. Version 3.0 offers intrusive and non-intrusive methods for propagating input uncertainties through computational models, tools for sensitivity analysis, methods for sparse surrogate construction, and Bayesian inference tools for inferring parameters from experimental data. This manual discusses the download and installation process for UQTk, provides pointers to the UQ methods used in the toolkit, and describes some of the examples provided with the toolkit.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

UQTk User Manual (V.3.1.0)

The UQ Toolkit (UQTk) is a collection of libraries and tools for the quantification of uncertainty in numerical model predictions. Version 3.1.0 offers intrusive and non-intrusive methods for propagating input uncertainties through computational models, tools for sensitivity analysis, methods for sparse surrogate construction, and Bayesian inference tools for inferring parameters from experimental data. This manual discusses the download and installation process for UQTk, provides pointers to the UQ methods used in the toolkit, and describes some of the examples provided with the toolkit.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

UQTk Version 3.1.1 User Manual

The UQ Toolkit (UQTk) is a collection of libraries and tools for the quantification of uncertainty in numerical model predictions. Version 3.1.1 offers intrusive and non-intrusive methods for propagating input uncertainties through computational models, tools for sensitivity analysis, methods for sparse surrogate construction, and Bayesian inference tools for inferring parameters from experimental data. This manual discusses the download and installation process for UQTk, provides pointers to the UQ methods used in the toolkit, and describes some of the examples provided with the toolkit.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

UQTk Version 3.0.4 User Manual

The UQ Toolkit (UQTk) is a collection of libraries and tools for the quantification of uncertainty in numerical model predictions. Version 3.0.4 offers intrusive and non-intrusive methods for propagating input uncertainties through computational models, tools for sensitivity analysis, methods for sparse surrogate construction, and Bayesian inference tools for inferring parameters from experimental data. This manual discusses the download and installation process for UQTk, provides pointers to the UQ methods used in the toolkit, and describes some of the examples provided with the toolkit.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

UQTk Version 3.1.2 User Manual

The UQ Toolkit (UQTk) is a collection of libraries and tools for the quantification of uncertainty in numerical model predictions. Version 3.1.2 offers intrusive and non-intrusive methods for propagating input uncertainties through computational models, tools for sensitivity analysis, methods for sparse surrogate construction, and Bayesian inference tools for inferring parameters from experimental data. This manual discusses the download and installation process for UQTk, provides pointers to the UQ methods used in the toolkit, and describes some of the examples provided with the toolkit.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Luminous Infrared Galaxies with the Submillimeter Array. V. Molecular Gas in Intermediate to Late-stage Mergers

We present new high-resolution ALMA ({sup 13}CO J = 1−0 and J = 2−1) and CARMA ({sup 12}CO and {sup 13}CO J = 1−0) observations of two luminous infrared galaxies (LIRGs), Arp 55 and NGC 2623. The new data are complementary to published and archival submillimeter array observations of {sup 12}CO J = 2−1 and J = 3−2. We perform a Bayesian likelihood non-local thermodynamic equilibrium analysis to constrain the molecular gas physical conditions such as temperature, column, and volume densities and the [{sup 12}CO]/[{sup 13}CO] abundance ratio. For Arp 55, an early/intermediate-staged merger, the line measurements are consistent with cold (∼10–20 K), dense (>10{sup 3.5} cm{sup −3}) molecular gas. For NGC 2623, the molecular gas is warmer (∼110 K) and less dense (∼10{sup 2.7} cm{sup −3}). Because Arp 55 is an early/intermediate stage merger, while NGC 2623 is a merger remnant, the difference in physical conditions may be an indicator of merger stage. Comparing the temperature and volume density of several LIRGs shows that the molecular gas, averaged over ∼kiloparsec scales, of advanced mergers is in general warmer and less dense than early/intermediate stage mergers. We also find that the [{sup 12}CO]/[{sup 13}CO] abundance ratio of NGC 2623 is unusually high (>250) when compared with the Milky Way; however, it follows a trend seen with other LIRGs in the literature. This high [{sup 12}CO]/[{sup 13}CO] value is very likely due to stellar nucleosynthesis enrichment of the interstellar medium. On the other hand, Arp 55 has a more Galactic [{sup 12}CO]/[{sup 13}CO] value with the most probable [{sup 12}CO]/[{sup 13}CO] value being 20–30. We measure the CO-to-H{sub 2} conversion factor, α{sub CO}, to be ∼0.1 and ∼0.7 (3 × 10{sup −4}/x{sub CO}) M {sub ⊙} (K km s{sup −1} pc{sup 2}){sup −1} for Arp 55 and NGC 2623, respectively. Because Arp 55 is an early/intermediate-stage merger, this suggests that the transition from a Galactic conversion factor to a LIRG value happens at an even earlier merger stage.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

A novel application of cosmic rays for spent nuclear fuel monitoring - 15738

Monitoring nuclear waste and controlling spent nuclear fuel at its source is one of the main strategies to minimize the risks of nuclear proliferation and reduce potential homeland threats. Storage of high level waste from a fission reactor typically involves placing spent fuel in sealed dry casks. After the spent nuclear fuel has been placed inside the dry cask, the cask is welded, not allowing for visual inspection. Conventional methods for examining the interior of materials e.g., X-rays, are limited by the fact that they cannot penetrate very dense well-shielded objects. Instead, cosmic ray muons have the potential to allow for assessment of spent nuclear fuel stored within sealed dense dry casks. Objectives: 1.Exploit the ability of cosmic ray muons to penetrate high density materials allowing the distribution of materials within the object to be inferred from the muon tracks. 2. Develop a monitoring system that tracks the path of cosmic ray muons, including their deflection and attenuation, to identify the contents of a closed dry cask. Methods: - Monte-Carlo simulations using Geant4; - Bayesian decision theory for signal classification; - Wavelet analysis for noise cancellation; - Chaotic Neural Networks for signal encryption. Results: 1. Preliminary statistical analysis proved that the method is feasible and has the ability to identify the presence of spent nuclear fuel in a dry cask. 2. Using monoenergetic muons, a dry cask with one fuel assembly missing can be correctly identified with error rate <1%. Broader merit: Once completed, the proposed project for cosmic ray muon transport will enable a new and potentially transformative approach towards efficient, inexpensive and remote cosmic ray muon applications in archaeology, medicine and safeguarding of nuclear materials.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗