Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Measurement of $d^2 \sigma/d|\vec{q}|d E_{avail}$ and 2p2h contribution using charged current $\nu_\mu$ interactions in the NOvA Near Detector

This Thesis presents the analysis methods and measurement of the $d^2 \sigma/d|\vec{q}|d E_{avail}$-CC inclusive cross section that describes neutrino scattering in a predominantly hydrocarbon medium. Studies that are prerequisite for the analysis including, variable determinations, event selections, efficiency, purity determinations, and systematic uncertainty estimations, are summarized. Tests of the analysis are also detailed, and shown to reproduce the input distributions. The double-differential cross-section measurement provides the foundation for determination of the inclusive rate for the excess induced by 2p2h processes together with nuclear medium effects that are not described by the Fermi gas model. The methodology for the extraction of this rate, defined as the excess observed relative to the rate estimated for known single-nucleon interactions, is described in detail. The following new measurements are reported in this Thesis: (1) The $\nu_\mu-$CC inclusive double-differential cross section as a function of three-momentum transfer and available hadronic energy, for an average $E_\nu$ of 1.8 GeV, is obtained. The double-differential cross section value that obtains when the bin widths are taken to be the dimensions of the entire analysis domain, e.g. $0.2 \le |\vec{q}| \le 2.0 GeV/c$ and $0.0 \le E_{avail} \le 2.0$ GeV, is inclusive cross section is $(6.05 \pm 0.75)\times10^{-39}$ cm$^2$/GeV/GeV/$c$/nucleon. An event excess attributed to 2p2h-MEC processes is observed; most of the rate occurs in a contiguous phase space region $0.3 \le |\vec{q}| \le 1.0$ GeV$/c$ $0.0 \le E_{avail} \le 0.35$ GeV. This rate represents $12.0\pm 6.5\%$ of the observed CC inclusive cross section; its cross section ratio relative to CCQE scattering is estimated to be $44.3 \pm 23.9\%$.

Olson, Travis Grant↗
DOE OSTI2021

An Assessment of Global Positioning System Velocity Uncertainty in California

We analyze data from 580 continuous global positioning system (GPS) stations in California to quantify differences in published velocity estimates from five analysis centers. Horizontal and vertical rates for individual stations can differ up to 5 mm/yr, with systematic differences in some areas comparable to deformation rates. Published velocity uncertainties vary between analysis centers and are systematically underreported in the horizontal relative to empirical uncertainties calculated from the scatter of analysis center velocities. In the vertical, published velocity uncertainties are both over and underreported and vary more widely between centers. An interpolated ensemble vertical velocity field shows high-subsidence regions in the Central Valley and Salton Trough have the largest empirical uncertainties, while station density has a modest impact on uncertainties. Applications that rely on subcentimeter GPS accuracy should consider the possibility that formal errors published with velocity rate estimates understate true velocity uncertainties in both the horizontal and vertical.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A mixed, unified forward/inverse framework for earthquake problems: fault implementation and coseismic slip estimate

SUMMARY We introduce a new finite-element (FE) based computational framework to solve forward and inverse elastic deformation problems for earthquake faulting via the adjoint method. Based on two advanced computational libraries, FEniCS and hIPPYlib for the forward and inverse problems, respectively, this framework is flexible, transparent and easily extensible. We represent a fault discontinuity through a mixed FE elasticity formulation, which approximates the stress with higher order accuracy and exposes the prescribed slip explicitly in the variational form without using conventional split node and decomposition discrete approaches. This also allows the first order optimality condition, that is the vanishing of the gradient, to be expressed in continuous form, which leads to consistent discretizations of all field variables, including the slip. We show comparisons with the standard, pure displacement formulation and a model containing an in-plane mode II crack, whose slip is prescribed via the split node technique. We demonstrate the potential of this new computational framework by performing a linear coseismic slip inversion through adjoint-based optimization methods, without requiring computation of elastic Green’s functions. Specifically, we consider a penalized least squares formulation, which in a Bayesian setting—under the assumption of Gaussian noise and prior—reflects the negative log of the posterior distribution. The comparison of the inversion results with a standard, linear inverse theory approach based on Okada’s solutions shows analogous results. Preliminary uncertainties are estimated via eigenvalue analysis of the Hessian of the penalized least squares objective function. Our implementation is fully open-source and Jupyter notebooks to reproduce our results are provided. The extension to a fully Bayesian framework for detailed uncertainty quantification and non-linear inversions, including for heterogeneous media earthquake problems, will be analysed in a forthcoming paper.

58 GEOSCIENCES↗
DOE OSTI2021

Critical Parameter Identification of Fuel-Cell Models Using Sensitivity Analysis

Numerical modeling has been a vital tool in proton-exchange-membrane fuel-cell (PEMFC) analysis; however, the predictive capabilities of these models depend on the input physical parameters, several of which are either not experimentally measured or have large scatter in measured values. This article presents an uncertainty propagation-based sensitivity analysis to identify the model parameters that impact the model predictions most. A comprehensive 2-D membrane electrode assembly (MEA) model is used to perform local sensitivity analysis at multiple operating conditions, which encompass the range of environments and operating conditions a cell can encounter. While at lower humidities, cathode kinetics and membrane-ohmic-loss related parameters are crucial, gas transport and porous-media saturation behavior are more important at humidified conditions. Several of these findings are different from previous studies presented in literature. Identifying the crucial parameters helps focus future material and cell optimization studies as well as experimental studies to quantify these parameters with higher accuracy.

25 ENERGY STORAGE↗
DOE OSTI2026

A Joint Search for Muon Neutrino Disappearance with the Short-Baseline Neutrino Program Using the SPINE Deep Learning-Based Reconstruction Package

We present the status of a joint search for muon neutrino disappearance in the Booster Neutrino Beam at Fermilab using the Short-Baseline Neutrino (SBN) Program's two-detector configuration, SBND and ICARUS. Charged-current interactions consistent with muon neutrinos and containing only a muon and at least one proton in the final state are reconstructed and selected using the SPINE deep learning-based particle reconstruction package. To exploit proton multiplicity information and enhance sensitivity to modeling effects, the selected sample is partitioned into exclusive channels with exactly one reconstructed proton and with more than one reconstructed proton. Comprehensive systematic uncertainties from the neutrino flux, interaction, and detector response are incorporated into the analysis, and the coverage of these systematic uncertainty models is validated using data from both detectors, including checks of near-far consistency in key kinematic and topology-sensitive observables. This analysis is intended for inclusion in SBN's first oscillation result.

Mueller, Justin [Fermilab]↗
DOE OSTI2025

Improved Joint FHC+RHC Sterile Analysis with Nu-on-E Sample

The NOvA experiment measures neutrino oscillations in a nearly pure muon neutrino beam. The neutral current (NC) and the muon neutrino charged current (CC) interaction disappearance are studied to search for active-to-sterile neutrino oscillations. NOvA consists of two functionally identical detectors placed 14 mrad off-axis of Fermilab's NuMI beamline. The Near Detector (ND) is located 100 m underground on-site at Fermilab, while the Far Detector (FD) sits on the surface 810 km away at Ash River, Minnesota. The NuMI beam is capable of running in forward or reverse horn current modes, creating a predominantly neutrino or antineutrino beam, respectively. For this analysis, a simultaneous near and far detector joint fit approach was utilized. Systematic uncertainties constrain the ND while the statistical uncertainties constrain the FD. For the upcoming analysis, we plan on including a dedicated neutrino-on-electron elastic scattering sample. This interaction is advantageous as it is subject to minimal cross-section uncertainties, which will help reduce degeneracies between other uncertainty categories. This thesis presents the sterile neutrino oscillation sensitivity results using 27.2e+20 POT in FHC mode, 12.6e+20 in RHC mode.

Haejun, Stella Oh [Cincinnati U.]↗
DOE OSTI2026

Inclusive Search for Anomalous Single-Photon Production in MicroBooNE

The MicroBooNE experiment is an 85-ton active volume liquid argon time projection chamber (LArTPC) neutrino detector situated in the Fermilab Booster Neutrino Beam. One of the main goals for MicroBooNE is to investigate the “low energy excess” (LEE) of electromagnetic events observed by the MiniBooNE experiment by leveraging the unique capabilities of the LArTPC technology to distinguish photons from electron electromagnetic showers. In this poster, we present an inclusive search for anomalous production of single-photon events from neutrino interactions in the MicroBooNE experiment, motivated by the MiniBooNE LEE. The analysis uses the Wire-Cell reconstruction framework to select a sample of inclusive single-photon final-state interactions and leverages simultaneous measurements of sidebands to constrain signal and background predictions and reduce uncertainties. We perform a blind analysis using a dataset collected from February 2016 to July 2018 from the Booster Neutrino Beam (BNB) at Fermilab. In the full signal region, we observe agreement between the data and the prediction, with a goodness-of-fit p-value of 0.11. We then isolate a sub-sample of these events containing no visible protons, and observe 93 $/pm$ 22(stat.) $\pm$ 35(syst.) data events above prediction, corresponding to just above 2$\sigma$ local significance, concentrated at shower energies below 600 MeV.

Yandel, Erin [Los Alamos] (ORCID:0000000277123709)↗
DOE OSTI2026

Determination of proton PDF uncertainties with Markov chain Monte Carlo

We present an analysis of parton distribution functions (PDFs) of the proton using Markov chain Monte Carlo (MCMC) methods. The MCMC approach naturally implements Bayes’ theorem and, thus, provides a means to directly sample the underlying probability distribution—in this case, the probability distribution of the PDF parameters. This allows for a straightforward propagation of the resulting uncertainties into any PDF-dependent observable, preserving their simple probabilistic interpretation. In our analysis we include a broad set of deep inelastic scattering data from HERA, BCDMS and NMC experiments along with the Drell-Yan, 𝑊 and 𝑍 boson data from LHC and Tevatron experiments, which combined with theoretical calculations at next-to-next-to-leading order in QCD allow for realistic determination of PDFs. The main focus of this analysis is to explore alternative methods for PDF uncertainty estimation that are more firmly grounded in statistical principles. We show that the flexibility of the Bayes framework, allowing one, e.g., to account for non-Gaussianity or inconsistencies of datasets, is crucial to extract realistic uncertainties when such assumptions are not fulfilled. We also demonstrate that MCMC allows one to determine the Δ⁢𝜒 2 value corresponding to a given confidence level in the sample, which can, in turn, be used as a statistically well-founded tolerance criterion used in the Hessian method, thus addressing one of its main long-standing drawbacks.

Risse, Peter Clemens [Universität Münster (Germany↗
DOE OSTI2024

An assessment of China’s methane mitigation potential and costs and uncertainties through 2060

Abstract China, the world’s largest methane emitter, is increasingly focused on methane mitigation in support of its climate goals, but gaps exist in the understanding of key methane sources, as well as mitigation opportunities and their associated uncertainties. We use a bottom-up modeling approach with updated methane emission projections and abatement cost analysis to account for additional sources, uncertainties, and mitigation measures in China’s energy and agricultural sectors. Here we show the significant cost-effective potential for reducing methane emissions in China by 2030, with 660 million tonnes of carbon dioxide equivalent possible with average negative abatement costs of US$6.40 per tonne CO 2 e. Most of this potential exists in the energy sector, particularly coal mining, but the greater potential will shift towards agriculture by 2060. Aquaculture and biochar applications in rice cultivation have net economic benefits but need greater support for deployment, while new mitigation measures will be needed for remaining emissions from enteric fermentation, rice cultivation, and wastewater.

Science & Technology - Other Topics↗
DOE OSTI2021

Gluon distributions and their applications to Ioffe-time distributions

We investigate unpolarized and polarized gluon distributions and their applications to the Ioffe-time distributions, which are related to lattice QCD calculations of parton distribution functions. Guided by the counting rules based on the perturbative QCD at large momentum fraction x and the color coherence of gluon couplings at small x, we parametrize gluon distributions in the helicity basis. By fitting the unpolarized gluon distribution, the inferred polarized gluon distribution from our parametrization agrees with the one from global analysis. A simultaneous fit to both unpolarized and polarized gluon distributions is also performed to explore the model uncertainty. The agreement with the global analysis supports the (1 - x) power suppression of the helicity-antialigned distribution relative to the helicity-aligned distribution. The corresponding Ioffe-time distributions and their asymptotic expansions are calculated from the gluon distributions. Our results of the Ioffe-time distributions can provide guidance to the extrapolation of lattice QCD data to the region lacking precise gluonic matrix elements. Therefore, they can help regulate the ill-posed inverse problem associated with extracting the gluon distributions from discrete data from first-principle calculations, which are available in a limited range of the nucleon momentum and the spatial separation between the gluonic currents. Given various limitations in obtaining lattice QCD data at large Ioffe time, phenomenological approaches can provide important complementary information to extract the gluon distributions in the entire momentum fraction region, especially at small x. The possibility of investigating higher-twist effects and other systematic uncertainties in the contemporary first-principle calculations of parton distributions from phenomenologically well-determined Ioffe-time distributions in the large Ioffe-time region is also discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Gamma Spectrometry Code Rodeo for Uranium Enrichment—FY 2022 Report

In the first two quarters of FY22, data acquisition continued at ORNL and LLNL using uranium sources of known enrichments. This was an FY21 task which could not be completed in FY21 because of problems encountered with the ORNL M400 CZT in Q4 of FY21, and the subsequent repairs. The detector was received back from H3D in the first of September 2021 , and the measurements resumed . Measurements using the repaired detector were completed in Q1 of FY22. The spectra were distributed by ORNL to the analyzing labs. Analysis results were received in Q2 of FY2022. The results from the various codes were intercompared and an ANOVA analysis was performed. Random and systematic uncertainties were established for each code. The ANOVA results and discussions were included in a revised version of FY21Annual Report issued in March 2022. A paper was presented at the INMM 2022Annual Conference, with the analysis results from the various isotopic codes, and the ANOVA table with random and systematic uncertainties. The Project Work Plan (PWP) for FY22 included a task to perform field testing of the M400 CZT and the analysis codes using UF 6 cylinder measurements at the Framatome Fuel Fabrication Facility in Richland, WA. PNNL was the lead for the field testing task. PNNL drafted a Field Test Plan, and refined it based on comments received from the team. PNNL coordinated with Framatome facility, the logistics of carrying out the field testing . A collimator and shield made out of TFlex (tungsten impregnated polymer) was designed and professionally manufactured. The collimators were used in the field test measurements. The measurements at Framatome were completed on April 21, 2022. A total of 34 Type 30B cylinders were measured using three M400 detectors (PNNL, LLNL, and ORNL detectors). Measurements using M400 were taken at two locations on the side of each cylinder and from the end-on bottom location. Additionally, HPGe measurements were taken at the end-on location to establish ground truth. Cylinder wall thickness measurements were also made at all three locations. To gain a better understanding of the effect of background from surrounding cylinders, the same five cylinders measured individually in low background locations were measured again in the cylinder storage yards. Due to inclement weather, manufacturer delays, shipping delays, and equipment failure, the measurement campaign spanned twice as long compared to the original timeline. Gamma-ray spectra from M400 and HPGe detectors, along with the cylinder data and photographs were organized and shared with the collaborators for further analysis. Spectra were analyzed by the participating laboratories. FY22 PWP also consists of tasks related to plutonium source measurements, adapting the codes to analyze plutonium spectra, and inter-comparison of the results from various codes. Plutonium spectra are being acquired at LANL, ORNL, and LLNL. LANL, SNL and LLNL are in the process of modifying FRAM, GADRAS, and CZTU, respectively. The plutonium related tasks will be completed in Q2 of FY23.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Estimation of process steady state with autoregressive models and Bayesian inference

To improve efficiency, separations engineers will typically design process circuits containing recirculating streams, which mix one or more of the process outputs with the feed material. Doing so can improve efficiency, but will cause a delay in the system reaching steady state conditions until the recirculating load mass flows stabilize. In testing separation circuits, engineers will often test a variety of factors and complete an analysis from sample results. Knowledge of if a process is at steady state, as well as the steady state conditions of a process, is essential for a valid techno-economic analysis. However, the definition of process steady state is often poorly defined, or does not include uncertainty quantification. If the performance of a process operating under two different sets of conditions are compared, an engineer who does not test for steady state or quantify steady state conditions risks producing a faulty analysis. In this work, a Bayesian statistical method for testing if all streams are at steady state is further motivated and then derived. Then after testing for steady state, the same model is used with a prior distribution that enforces a steady state assumption to estimate steady state conditions. Further, these methods were validated in a solvent extraction pilot plant where steady state conditions for all outflows were inferred with uncertainty quantification. Analysis is completed with functions available to the reader as part of the BayesMassBal (V 1.1.0) software package written in R.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2021

An Error Analysis Toolkit for Binned Counting Experiments

We introduce the MINERvA Analysis Toolkit (MAT), a utility for centralizing the handling of systematic uncertainties in HEP analyses. The fundamental utilities of the toolkit are the MnvHnD, a powerful histogram container class, and the systematic Universe classes, which provide a modular implementation of the many universe error analysis approach. These products can be used stand-alone or as part of a complete error analysis prescription. They support the propagation of systematic uncertainty through all stages of analysis, and provide flexibility for an arbitrary level of user customization. This extensible solution to error analysis enables the standardization of systematic uncertainty definitions across an experiment and a transparent user interface to lower the barrier to entry for new analyzers.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Policy support and technology development trajectory for renewable natural gas in the U.S.

Renewable natural gas (RNG) is a clean alternative to fossil natural gas, which can be used as transportation fuel, among other applications. This study projects the development trajectory of RNG and evaluates its impacts on the future U.S. transportation market using a hybrid computable general equilibrium model. This analysis considers various factors and uncertainties affecting RNG production, such as technology development, market conditions, competition with other advanced biofuels, and national and state policies. In 2050, RNG production will grow to 2.7 billion gallons (10 billion liters), mostly from swine manure, under current policy provisions. This will lead to a reduction in greenhouse gas (GHG) emissions by 58.56 million metric tonne of CO 2e in 2050. Analysis of different technology cases finds RNG from animal manure to be predominant, while RNG from corn stover and cellulosic ethanol are less competitive. Furthermore, a high mandatory target of 1 billion gallons will drive RNG production higher by 8–18 %, while an extended 2 nd -generation biofuel production tax credit will mostly increase cellulosic ethanol production. The model also finds RNG production being affected by uncertainties in market conditions, such as GDP growth, fossil fuel prices, and oil and gas supply.

Biomethane↗
DOE OSTI2025

Using Simulation to Interpret LArS's Results

Scintillation light analysis in liquid argon based neutrino detectors is restrained in capability due to uncertainty in fundamental constants critical to the analysis process. One such property is the Rayleigh scattering length of liquid argon. In the fall of 2023, the TallBo cryostat, located in the Noble Liquid Testing Facility (NLTF) at Fermilab, was used to study the scattering length of liquid argon in the Liquid Argon Scattering (LArS) experiment. Due to systematic errors unknown during measurement analysis, LArS s measurements were quite uncertain. By simulating the LArS experiment, we found that the downturn in detector count rate as a function of liquid argon height at low heights was caused by a misplaced silicon photo multiplier (SiPM). With concentrated effort, we may be able to successfully correct this effect by understanding the relationship between the specified Rayleigh scattering length and the measured attenuation length. With this information and further progression in analysis, we may be able to extract corrected measurements from the LArS data and attain a tangible experimental measurement of the scattering length of liquid argon.

Breaux, Auto D. [Tulane U.]↗