Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Bayesian uncertainty 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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Reduced scale stripline platform to extend accessible pressures on the Z machine: Shockless compression of platinum to 650 GPa

Reaching astrophysically relevant high energy density (HED) material states in the laboratory is an ongoing effort at multiple experimental facilities. We have developed a new dynamic compression platform for the Z Pulsed Power Facility that allows for sample sizes 100s of [Formula: see text]m in thickness that accommodate multiple grains in order to fully capture bulk properties, such as material strength. A pair of experiments compressed platinum (Pt) to HED conditions and conventional inverse Lagrangian analysis as well as a recent Bayesian calibration technique were used to determine the principal isentrope to 650 GPa with density uncertainties of <2%. These low uncertainties are calculated for single sample experiments, presenting the possibility of even smaller experimental uncertainties with multiple samples the platform allows. Our new platform extends the accessible Pt ramp pressures achievable on the Z machine to over 80% of the pressure recently achieved using the National Ignition Facility planar Hohlraum platform. This new capability, the next generation evolution of the stripline platform, was made possible by advancements in both our understanding of the Z pulsed power driver and our overall magnetohydrodynamic modeling capabilities.

Porwitzky, Andrew↗
DOE OSTI2022

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↗
DOE OSTI2025

Elastic Bayesian Model Calibration

Functional data are ubiquitous in scientific modeling. For instance, quantities of interest are modeled as functions of time, space, energy, density, etc. Uncertainty quantification methods for computer models with functional response have resulted in tools for emulation, sensitivity analysis, and calibration that are widely used. However, many of these tools do not perform well when the computer model’s parameters control both the amplitude variation of the functional output and its alignment (or phase variation). This paper introduces a framework for Bayesian model calibration when the model responses are misaligned functional data. The approach generates two types of data out of the misaligned functional responses: (1) aligned functions so that the amplitude variation is isolated and (2) warping functions that isolate the phase variation. These two types of data are created for the computer simulation data (both of which may be emulated) and the experimental data. The calibration approach uses both types so that it seeks to match both the amplitude and phase of the experimental data. The framework is careful to respect constraints that arise, especially when modeling phase variation, and is framed in a way that it can be done with readily available calibration software. In conclusion, we demonstrate the techniques on two simulated data examples and on two dynamic material science problems: a strength model calibration using flyer plate experiments and an equation of state model calibration using experiments performed on the Sandia National Laboratories’ Z-machine.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Bayesian inference for plasmonic nanometrology

Here, we introduce a Bayesian method for the characterization of plasmonic nanoparticles, which is applicable to both near- and far-field problems. Designed to combine data generated from any photon-plasmon interaction experiment with physically motivated theoretical models, our approach leverages state-of-the-art Markov chain Monte Carlo sampling techniques and returns parameter estimates on nanometric scales. Simulated spectral data sets, describing resonant scattering of photons from ellipsoidal and toroidal nanoparticles, are explored as concrete examples of our approach, with the resulting Bayesian estimates showing excellent agreement with the ground truth, even under conditions of high statistical noise. By incorporating Bayes factors into the method as well, we reveal how model selection can determine which one of competing geometric shapes better explains the observed data. Our comprehensive nanometrology procedure can be tailored to a variety of light-particle interaction models, and its reliance on Bayesian inference furnishes automatic uncertainty quantification. In addition to applicability to a host of plasmonic configurations such as nanoparticle dimers, trimers, and array studies, it is proposed that the presented analysis can be extended to the quantum regime, where nonclassical photon statistics may provide additional insight for inference of scatterer properties.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Quantifying uncertainty in Pareto estimates of global lake area

Abstract Size is a critical factor determining the rate and occurrence of specific lake processes such as carbon sequestration and greenhouse gas emissions and emerging evidence suggests that small lakes in particular have particularly large CO 2 flux rates. Because we do not have a complete census of all lakes, upscaling estimates of such processes to small lakes at broad spatial scales requires the use of lake size‐abundance distributions rather than empirical measurements of area. Existing lake census efforts are incomplete such that as lakes become smaller, they are more likely to be omitted either because they are too small to be resolved from remote sensing products or because of limited ground surveying effort (i.e., “censoring” of small lakes relative to large lakes). The present study explores one potential shortcoming of prior approaches estimating global lake area using lake size‐abundance distributions. Namely, that these prior approaches rely on frequentist curve fitting techniques combined with an ad‐hoc cutoff determination strategy (visual inspection to determine a likely censoring point). This yields an over‐exact lake area estimate that is typically reported with no uncertainty bounds. I show how these shortcomings can be addressed with a Bayesian model that produces larger estimates of lake area uncertainty relative to the typical approach. When used as part of a sensitivity analysis, such an approach has the potential to enable more robust intercomparisons among studies of aquatic processes upscaling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Boosting efficiency and reducing graph reliance: Basis adaptation integration in Bayesian multi-fidelity networks

The computational cost of high-fidelity numerical models makes outer-loop analysis, which requires repeated interrogation of the model such as uncertainty quantification, computationally demanding. Multi-fidelity methods, which construct a surrogate model using data from an ensemble of models of varying cost and accuracy, can substantially reduce the cost of outer-loop analysis. However, these methods can be difficult to apply when the model ensemble does not admit a clear hierarchy a priori and the correlations between models are low. Consequently, in this paper, we present a multi-fidelity method that leverages dimension reduction to enhance the correlation between models, thereby reducing the amount of data needed to train a surrogate from an unordered ensemble of models. Our method utilizes basis adaptation to build low-dimensional polynomial chaos expansions of each model and employs Multi-fidelity Networks to encode the relationships among models. We show that the resulting method exhibit two notable advantages over its counterpart: (1) enhanced accuracy (both reduced bias and variance); and (2) reduced dependency on the graph structure encoding relationships among models. We demonstrate the approach on an analytical test problem and a challenging finite element model for a spent nuclear fuel. Our method produces a surrogate model that is significantly more accurate than either a single-fidelity surrogate or a multi-fidelity surrogate constructed without basis adaptation.

42 ENGINEERING↗
DOE OSTI2021

Analysis and Optimization of Seismo-Acoustic Monitoring Networks with Bayesian Optimal Experimental Design

The Bayesian optimal experimental design (OED) problem seeks to identify data, sensor configurations, or experiments which can optimally reduce uncertainty. The goal of OED is to find an experiment that maximizes the expected information gain (EIG) about quantities of interest given prior knowledge about expected data. Therefore, within the context of seismic monitoring, we can use Bayesian OED to configure sensor networks by choosing sensor locations, types, and fidelity in order to improve our ability to identify and locate seismic sources. In this work, we develop the framework necessary to use Bayesian OED to optimize the ability to locate seismic events from arrival time data of detected seismic phases. In order to do utilize Bayesian OED we must develop four elements:1. A likelihood function that describes the uncertainty of detection and travel times; 2. A Bayesian solver that takes a prior and likelihood to identify the posterior; 3. An algorithm to compute EIG; and, 4. An optimizer that finds a sensor network which maximizes EIG. Once we have developed this framework, we can explore many relevant questions to monitoring such as: how and what multiphenomenology data can be used to optimally reduce uncertainty, how to trade off sensor fidelity and earth model uncertainty, and how sensor types, number, and locations influence uncertainty

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Impact of multiplicity fluctuations on entropy scaling across system size

The initial state is one of the greatest uncertainties in heavy-ion collisions. A model-agnostic approach is taken in the phenomenological Trento framework which constrains parameters using Bayesian analysis. However, the color-glass condensate (CGC) effective theory predicts initial energy densities that lie outside the recent Bayesian analyses due, in part, to the assumption in Trento of event-by-event multiplicity fluctuations following a Γ distribution. We compare the Trento-preferred √T A T B scaling to CGC-like T A T B scaling coupled with log-normal fluctuations in AuAu and dAu collisions and find there is a significant impact on the multiplicity distributions and on the eccentricities, which may affect the extraction of viscosity in small systems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Atikokan Digital Twin: Machine learning in a biomass energy system

The Atikokan Generating Station, operated by Ontario Power Generation, has a 200 MW, biomass-fired tower boiler that operates on a dispatch schedule with a five-minute cycle. The boiler is generally operated in the range of 40–100 MW using two of five burner levels. In order to optimize boiler performance, we propose the implementation of a unique digital twin. Our digital twin abstraction couples Bayesian inference from science-based models and from observations (machine learning) with decision theory to predict operating-variable set points that optimize the physical asset (the boiler) in the presence of uncertainty (artificial intelligence). We focus this paper on the continuous Bayesian machine learning part of the Atikokan Digital Twin; we discuss decision theory in a companion paper. We identify and learn about 12 operational, model, and measured-output parameters and their uncertainties from high-fidelity, science-based simulations of the Atikokan boiler and from the observed measurements at the power plant. Since the goal of the Atikokan Digital Twin is to implement it online in real time, we require fast function evaluations for the quantities of interest extracted from the simulations in the Bayesian analysis. We use Gaussian process regression/interpolation to create accurate, robust surrogate models. We define the Bayesian priors and likelihood function and solve for the posterior distributions of the 12 parameters. Here we then propagate these distributions (i.e., parameters with uncertainty) into the predicted distributions of 790 quantities of interest to learn about the relative importance of various sources of error including experimental, model, and operating-parameter errors.

09 BIOMASS FUELS↗
DOE OSTI2022

BayesMT: A Probabilistic Bayesian Framework for the Seismic Moment Tensor

Moment tensors (MTs) have long been used in earthquake and explosion source analysis, and there has been a renewed interest in how they can inform us about the seismic source, particularly in the geophysical monitoring community due to its application in event identification and yield analysis. However, parameter uncertainties in seismic MT inversion are rarely available. The inverse procedure often does not quantify MT model errors such as event location, data noise and Earth model that are essential for estimating solution robustness. To address this need, we propose to adopt the Bayesian probabilistic framework to incorporate uncertainties in MT inversions. In this study, we present the theoretical background of a probabilistic Bayesian framework for MT inversion accounting for model and measurements errors and illustrate the implementation of the method using a synthetic example.

58 GEOSCIENCES↗
DOE OSTI2022

A Bayesian Deep Learning Approach to Near-Term Climate Prediction

Since model bias and associated initialization shock are serious shortcomings that reduce prediction skills in state-of-the-art decadal climate prediction efforts, we pursue a complementary machine-learning-based approach to climate prediction. The example problem setting we consider consists of predicting natural variability of the North Atlantic sea surface temperature on the interannual timescale in the pre-industrial control simulation of the Community Earth System Model. While previous works have considered the use of recurrent networks such as convolutional LSTMs and reservoir computing networks in this and other similar problem settings, we currently focus on the use of feedforward convolutional networks. In particular, we find that a feedforward convolutional network with a Densenet architecture is able to outperform a convolutional LSTM in terms of predictive skill. Next, we go on to consider a probabilistic formulation of the same network based on Stein variational gradient descent and find that in addition to providing useful measures of predictive uncertainty, the probabilistic (Bayesian) version improves on its deterministic counterpart in terms of predictive skill. Finally, we characterize the reliability of the ensemble of machine learning models obtained in the probabilistic setting by using analysis tools developed in the context of ensemble numerical weather prediction.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

A Bayesian inferencing framework for ultrasound wave speed measurements in metal additive manufacturing

Process-related changes during metal additive manufacturing introduce microstructural variability in the material properties of printed parts, directly affecting component reliability. Accurate estimation of these property variations with part performance are essential for quality assurance. Ultrasound testing offers a non-destructive means to estimate mechanical properties and detect defects; however, conventional analysis methods often neglect the influence of microstructural variability, limiting their effectiveness. Here, this research presents a Bayesian inference technique for quantifying wave speed uncertainty from ultrasound measurements of metal additive manufactured parts. By integrating prior ultrasound data with a Bayesian model, the proposed approach generates posterior density estimates of wave speed that systematically account for manufacturing-induced variability and uncertainty. The novelty of this research lies in applying a Bayesian framework to analyze experimental ultrasound measurements within the context of metal additive manufacturing variability. The method enhances the accuracy of wave speed estimation by 64%, defect position by 50% and increases confidence associated with wave speed variance by 30% across different porosity levels, thereby providing a robust foundation for improved decision-making and increased reliability in additively manufactured components.

Additive manufacturing↗
DOE OSTI2022

Compressed Natural Gas Component Leak Frequency Estimation

The frequency of unintended releases in a compressed natural gas system is an important aspect of the system quantitative risk assessment. The frequencies for possible release scenarios, along with engineering models, are utilized to quantify the risks for compressed natural gas facilities. This report documents component leakage frequencies representative of compressed natural gas components that were estimated as a function of the normalized leak size. A Bayesian statistical method was used which results in leak frequency distributions for each component which represent variation and uncertainty in the leak frequency. The analysis shows that there is high uncertainty in the estimated leak frequencies due to sparsity in compressed natural gas data. These leak frequencies may still be useful in compressed natural gas system risk assessments, as long as this high uncertainty is acknowledged and considered appropriately.

03 NATURAL GAS↗
DOE OSTI2020

Bayesian analysis of the Zn 70 ( d , He 3 ) Cu 69 transfer reaction

Transfer reactions provide information about the single-particle nature of nuclear levels. In particular, the differential cross sections from these measurements are sensitive to the angular momentum of the transferred particle and the spectroscopic factor of the populated level. However, the process of extracting these properties is subject to uncertainties, both from experimental and theoretical sources. By integrating the distorted wave Born approximation into a Bayesian model, we propagate these uncertainties through to the spectroscopic factors and orbital angular momentum values. We use previously reported data of the proton pickup reaction Zn 70 ( d , He 3 ) Cu 69 as an example. By accounting for uncertainties in the experimental data, optical model parameters, and reaction mechanism, we find that the extracted spectroscopic factors for low-lying states of 69 Cu are subject to large, asymmetric uncertainties ranging from 35 to 108%. Additionally, Bayesian model comparison is employed to assign probabilities to each of the allowed angular momentum transfers. This method confirms the assignments for many states, but suggests that the data for a state lying at 3.70 MeV are better characterized by an ℓ =3 transfer, rather than the previously reported ℓ =2 .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Uncertainty Quantification and Sensitivity Analysis for Simulation of Hostile Blast Events [Slides]

A mesh convergence study shows that shock arrival times converge for basic flow shots. Sensitivity study reveals a surprisingly high sensitivity to the placement of the HE package. Bayesian inference tools used on the experimental datasets suggest this is a real effect. Sensitivity analysis of ambient conditions such as temperature and pressure reveal minimal effect. Porting of pressure data from QUINOA simulations to perform structural analysis on the aeroshell geometry was successful, and shows a strong dependence of payload response to the angle of attack.

42 ENGINEERING↗
DOE OSTI2020

Integrated Data Analysis to Improve Te and Zeff Measurement on NSTX-U

The multi-energy soft x-ray (ME-SXR) camera is a new concept for a fusion plasma diagnostic based on a pixelated detector with adjustable photon cutoff thresholds. It is a single, compact, solid-state detector which combines spatial and spectral sensitivity and can be configured remotely at run time, making it ideal for a fusion environment. The detector can be configured for application to a wide range of physics problems including impurity transport, electron temperature profiles, and non-thermal effects. The detector has been installed on the Madison Symmetric Torus (MST) as part of a collaboration between Princeton Plasma Physics Laboratory and the University of Wisconsin–Madison to develop the technology and analysis techniques. A full diagnostic forward model of the ME-SXR diagnostic has been developed to enable analysis of the raw photon count data. This model incorporates the atomic physics underlying the plasma emission, diagnostic geometry, and detector response characteristics. Employing this model, a Bayesian analysis approach is used to extract electron temperature from the raw data. This analysis approach uses all available information, clearly quantifies prior information, and handles uncertainties transparently. The ME-SXR camera is being used to study the 3D (helical, stellarator-like) state that forms in the MST device. This helical core structure has improved energy confinement. The transition to the helical state is accompanied by a rise in electron temperature in the central part of the helical structure, but then this enhancement fluctuates. The characteristics of this fluctuation are being compared to a magneto-hydrodynamics model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Harnessing Uncertainty through Functional Data Analysis in Gas Breakthrough Data

Detecting subsurface explosions from radionuclide gas migration through rock fractures is an effective way to identify nuclear activity. Los Alamos National Laboratory (LANL) has developed simulation methods, based on data from the 1962 Hardhat underground nuclear test, to predict gas breakthrough times at the surface. However, these methods rely on an imperfect understanding of the relationship between rock damage and fracture permeability. Our clinic project studies methods for predicting breakthrough curves that characterize total mass produced as a function of time, as well as quantifying the uncertainty associated with these predictions. The model that is currently employed to relate damage to permeability uses an empirically motivated power-law expression, with a range of parameter values that are compatible with the experimental Hardhat data. We develop emulators, built from functional data analysis techniques and trained on simulation data, that rapidly predict the gas breakthrough curve given a damage field and given the parameter values of the power-law equation. Using Bayesian regression, we address the problem of uncertainty quantification in our emulators. Finally, in order to test the robustness of the model, we further validate it on a damage field representing different physical conditions.

58 GEOSCIENCES↗
DOE OSTI2023

Absolute cross section of the C 12 ( p , γ ) N 13 reaction

Solar neutrino measurements have recently reached a level of sensitivity such that CNO fluxes can now be experimentally determined. While these first measurements are still only sensitive to the higher energy neutrinos resulting from the β + decays of 15 Ο produced by the 14 Ν ⁡( p ,$\gamma$)⁢ 15 O reaction, future measurements will work towards detection of neutrinos from the β + decay of 13 N from the 12 C ⁡( p ,$\gamma$)⁢ 13 N reaction. Here, this paper reports on a recent measurement of the 12 C ⁡( p ,$\gamma$)⁢ 13 N reaction covering a broad laboratory energy range between 1.0 and 2.5 MeV. The measurement was made to better determine the overall normalization of the absolute cross section and to explore the interference effects between the two broad, overlapping resonances at proton energies of 0.460 and 1.689 MeV and the direct capture to the ground state of 13 N in the framework of a multichannel R -matrix analysis. This work takes into account previous radiative capture as well as elastic 12 C ( p , p ) 12 C scattering data, making uncertainty estimations using a Bayesian framework, to determine a reliable extrapolation of the low energy S factor towards the stellar energy range of CNO hydrogen burning. These new experimental results, and a detailed investigation of the past literature data, suggest that the resonant component of the cross section should be 30% lower than previously accepted.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗