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 1,207 records · Page 67
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

BeyondPlanck: XII. Cosmological parameter constraints with end-to-end error propagation

We present cosmological parameter constraints estimated using the Bayesian BEYONDPLANCK analysis framework. This method supports seamless end-to-end error propagation from raw time-ordered data onto final cosmological parameters. As a first demonstration of the method, we analyzed time-ordered Planck LFI observations, combined with selected external data (WMAP 33–61 GHz, Planck HFI DR4 353 and 857 GHz, and Haslam 408 MHz) in the form of pixelized maps that are used to break critical astrophysical degeneracies. Overall, all the results are generally in good agreement with previously reported values from Planck 2018 and WMAP, with the largest relative difference for any parameter amounting about 1σ when considering only temperature multipoles between 30 ≤ ℓ ≤ 600. In cases where there are differences, we note that the BEYONDPLANCK results are generally slightly closer to the high-ℓ HFI-dominated Planck 2018 results than previous analyses, suggesting slightly less tension between low and high multipoles. Using low-ℓ polarization information from LFI and WMAP, we find a best-fit value of τ = 0.066 ± 0.013, which is higher than the low value of τ = 0.052 ± 0.008 derived from Planck 2018 and slightly lower than the value of 0.069 ± 0.011 derived from the joint analysis of official LFI and WMAP products. Most importantly, however, we find that the uncertainty derived in the BEYONDPLANCK processing is about 30 % greater than when analyzing the official products, after taking into account the different sky coverage. We argue that this uncertainty is due to a marginalization over a more complete model of instrumental and astrophysical parameters, which results in more reliable and more rigorously defined uncertainties. We find that about 2000 Monte Carlo samples are required to achieve a robust convergence for a low-resolution cosmic microwave background (CMB) covariance matrix with 225 independent modes, and producing these samples takes about eight weeks on a modest computing cluster with 256 cores.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Numerical Simulations of Flow-Induced Deflections in MITR LEU Fuel Plate Due to Channel Size Disparity

The hydromechanical stability of the fuel plates in parallel coolant channels of a Materials Testing Reactor (MTR) fuel element design is of great importance to the safety of research and test reactors. Previous analytical, experimental, and numerical efforts focused on parallel channels with the same or similar size; also, in the prior numerical simulations, the fuel plate was often assumed to be perfectly flat. This work presents the results of a fluid-structure interaction simulation performed to evaluate the flow-induced deflections of the fuel plates in the low-enriched uranium (LEU, <20 wt% 235 U) fuel element design for the conversion (from highly enriched uranium) of the Massachusetts Institute of Technology Reactor (MITR-II, also referred to as MITR). Various manufacturing and assembly tolerances of the MITR LEU elements are considered in the analysis, and the effects of channel size disparity, nonideal plate shape, and flow rate uncertainty are investigated. Results show that, for all cases analyzed, the deflection occurs toward the larger channel, and the change in any channel stripe remains small (less than 0.021 mm) compared to fabrication tolerances. In addition to simulation work, a hydraulic performance test of the MITR LEU fuel element is currently planned to support conversion to the use of LEU fuel.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Angular systematics-free cosmological analysis of galaxy clustering in configuration space

ABSTRACT Galaxy redshift surveys are subject to incompleteness and inhomogeneous sampling due to the various constraints inherent to spectroscopic observations. This can introduce systematic errors on the summary statistics of interest, which need to be mitigated in cosmological analysis to achieve high accuracy. Standard practices involve applying weighting schemes based on completeness estimates across the survey footprint, possibly supplemented with additional weighting schemes accounting for density-dependent effects. In this work, we concentrate on pure angular systematics and describe an alternative approach consisting in analysing the galaxy two-point correlation function where angular modes are nulled. By construction, this procedure removes all possible known and unknown sources of angular observational systematics, but also part of the cosmological signal. We use a modified Landy–Szalay estimator for the two-point correlation function that relies on an additional random catalogue where angular positions are randomly drawn from the galaxy catalogue, and provide an analytical model to describe this modified statistic. We test the model by performing an analysis of the full anisotropic clustering in mock catalogues of luminous red and emission-line galaxies at 0.43 < z < 1.1. We find that the model fully accounts for the modified correlation function in redshift space, without introducing new nuisance parameters. The derived cosmological parameters from the analysis of baryon acoustic oscillations and redshift-space distortions display slightly larger statistical uncertainties, mostly for the growth rate of structure parameter fσ8 that exhibits a $50{{\ \rm per\ cent}}$ statistical error increase, but free from angular systematic error.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Joint Determination of Reactor Antineutrino Spectra from U 235 and Pu 239 Fission by Daya Bay and PROSPECT

A joint determination of the reactor antineutrino spectra resulting from the fission of 235U and 239Pu has been carried out by the Daya Bay and PROSPECT collaborations. This Letter reports the level of consistency of 235U spectrum measurements from the two experiments and presents new results from a joint analysis of both data sets. The measurements are found to be consistent. The combined analysis reduces the degeneracy between the dominant 235U and 239Pu isotopes and improves the uncertainty of the 235U spectral shape to about 3%. The 235U and 239Pu antineutrino energy spectra are unfolded from the jointly deconvolved reactor spectra using the Wiener-SVD unfolding method, providing a data-based reference for other reactor antineutrino experiments and other applications. This is the first measurement of the 235U and 239Pu spectra based on the combination of experiments at low- and highly enriched uranium reactors.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Multiprobe constraints on early and late time dark energy

We perform a multiprobe analysis combining cosmic microwave background (CMB) data from Planck and the Atacama Cosmology Telescope (ACT), ACT CMB lensing, and large-scale structure (LSS) measurements from the Dark Energy Spectroscopic Instrument (DESI), including DESI Legacy Imaging Survey (LS) galaxies and baryon acoustic oscillations (BAOs). We present the first 5×2 pt analysis of ACT DR6 lensing, DESI LS, and Planck Integrated Sachs-Wolfe (ISW). Within Λ⁢CDM, this yields 𝑆 8 =𝜎 8 ⁢(Ω 𝑚 /0.3) 0.5 =0.819 ±0.016, in good agreement with primary CMB inferences and provides a sound-horizon-free Hubble constant constraint of 𝐻 0 =70.0±4.4 km s −1 Mpc −1 . Then, combining with CMB primary and BAO, we reconfirm a CMB–BAO discrepancy in the Ω 𝑚 –$\frac{𝐷𝑣}{𝑟𝑑}$ plane, which is heightened when combining BAO with the 5 ×2 pt data vector. We explore two dark-energy extensions that may reconcile this: an early-time modification, early dark energy (EDE), and late-time dynamical dark energy (DDE) parametrized by 𝑤 0 ⁢𝑤 𝑎 . For CMB primary +BAO +5 ×2 pt, we find a 3.3⁢𝜎 preference for DDE over Λ⁢CDM, while EDE is modestly favored at 2.3⁢𝜎. The models address different shortcomings of Λ⁢CDM: DDE relaxes the neutrino mass bound (𝑀 𝜈 <0.17 eV vs <0.050 eV under Λ⁢CDM), making it compatible with neutrino oscillation measurements, while EDE raises the Hubble constant to 𝐻 0 =70.5±1.2 km s −1 Mpc −1 , easing the discrepancy with SH0ES. However, neither model resolves both issues simultaneously. Our analysis indicates that both DDE and EDE remain viable extensions of Λ⁢CDM within current uncertainties and demonstrates the capacity of combined probes to place increasingly stringent constraints on cosmological parameters.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Low Energy Analysis of NuclEar Reactions v0.0.1

This is a lightweight, easy to use Python analysis package for analyzing public low-energy nuclear reaction data relevant to understanding basic nuclear processes, often of relevance to Solar Fusion, to improve our theoretical understanding of them, and enable connections with input from lattice QCD. It performs the analysis in a Bayesian Framework and supports Bayesian Model Averaging to provide a robust uncertainty quantification.

Walker-Loud, André↗
DOE OSTI2024

Analysis of the Challenges in Developing Sample-Based Multi-fidelity Estimators for Non-deterministic Models

Multifidelity (MF) uncertainty quantification (UQ) seeks to leverage and fuse information from a collection of models to achieve greater statistical accuracy with respect to a single-fidelity counterpart, while maintaining an efficient use of computational resources. Despite many recent advancements in MF UQ, several challenges remain and these often limit its practical impact in certain application areas. In this manuscript, we focus on the challenges introduced by nondeterministic models to sampling MF UQ estimators. Nondeterministic models produce different responses for the same inputs, which means their outputs are effectively noisy. MF UQ is complicated by this noise since many state-of-the-art approaches rely on statistics, e.g., the correlation among models, to optimally fuse information and allocate computational resources. Here, we demonstrate how the statistics of the quantities of interest, which impact the design, effectiveness, and use of existing MF UQ techniques, change as functions of the noise. With this in hand, we extend the unifying approximate control variate framework to account for nondeterminism, providing for the first time a rigorous means of comparing the effect of nondeterminism on different multifidelity estimators and analyzing their performance with respect to one another. Numerical examples are presented throughout the manuscript to illustrate and discuss the consequences of the presented theoretical results.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Comparative study of decay heat calculations with FLUKA and MCNP/CINDER2008

In designing a high-power tungsten target, decay heat driven temperature rise in the spallation volume is a safety concern during maintenance and in loss of coolant accidents. As tungsten hydrates and becomes volatile in steam at above 800 ◦C, it is important to keep the target temperature below this threshold when active cooling is unavailable. Decay heat in a target is calculated with particle transport simulation codes combined with transmutation codes. The calculated decay heat usually differs depending on the nuclear cross sections and the decay particle transport models built in the code architecture. In this paper, we calculated decay heat of a water-cooled tungsten target using popular particle transport codes, FLUKA and MCNP6® paired with CINDER2008 and compared the results. The target-moderator-reflector (TMR) system is modeled with a water-cooled solid tungsten target, water premoderators, liquid hydrogen cold moderators and beryllium reflectors. Water-cooled stainless-steel shielding is modeled around the TMR system. The tungsten volume is clad with a thin layer of erosion/corrosion resistant material. This study provides information about the uncertainty range in decay heat prediction of high-power spallation targets for hazard analysis.

Lee, Yong Joong↗
DOE OSTI2025

SSAPy - Space Situational Awareness for Python

SSAPy is a fast and flexible orbit modeling and analysis tool for orbits spanning from low-Earth into the cislunar regime. Orbits can be flexibly specified from common input formats such as Keplerian elements or two-line element (TLE) data files. SSAPy allows users to model satellites and specify parameters such as satellite area, mass, and drag coefficients. SSAPy includes a customizable force-propagation with a range of Earth, Lunar, radiation, atmospheric, and maneuvering models. SSAPy makes use of various community integration methods and can calculate time-evolved orbital quantities, including satellite magnitudes and state vectors. Users can specify various space- and ground-based observation models with support for multiple coordinate and reference frames. SSAPy also supports orbit analysis and propagation methods such as multiple hypothesis tracking and has built-in uncertainty quantification. The majority of SSAPy’s methods are vectorized and parallelizable, allowing for effective use of high-performance computer (HPC) systems. Finally, SSAPy has plotting functionality, allowing users to visualize orbits and trajectories. Examples are shown in Figure 1 and Figure 2.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Accounting for uncertainty sources in pyrometric measurements: Variables affecting Surface Temperature measurements using the LLNL Streak Spectrometers

Here we present a detailed analysis of the error sensitivities of the LLNL temperature analysis package used in analyzing streak spectrometer data. These sensitivities include noise, temperature, wavelength, emissivity, missing bands, and missing data points. The analysis package was tested for its ability to accurately reconstruct randomly generated noisy data and found to perform within specifications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Progress on the R-matrix Analysis for the n+ 181 Ta Evaluation [Slides]

This presentation discusses how nuclear data evaluations can influence the accuracy of computational nuclear design calculations. It also discusses the motivation for a new evaluation of 181 Ta. The goals of the evaluations are to generate a set of resonance parameters in the Resolved Resonance Region (RRR) that extends up to 2.6 keV, to generate a set of average resonance parameters in the Unresolved Resonance Region (URR) up to 100 keV (derived from the RRR parameters), and to generate covariance information for both RRR and URR. The presentation discusses how several different experimental data sets will be fit simultaneously in the evaluation of 181 Ta. It states that the level density of 181 Ta can make R-matrix fitting a challenge as the width of the levels increase with energy (as well as decreasing experimental resolution) and the measured resonances can overlap. The presentation concludes by discussing R-matrix analysis and covariance information as the latter is important when propagated to quantify the uncertainty on output responses of nuclear design calculations such as the multiplication factor k .

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2024

Integrating Atmospheric Specifications into Seismoacoustic Event Localization

This report investigates the integration of infrasound and seismic data to improve event localization accuracy, specifically focusing on a surface explosion at the Utah Training and Testing Range (UTTR). Utilizing the Seismoacoustic Bayesian Event Locator (SABEL) framework, we incorporated atmospheric specifications derived from Ground to Space (G2S) profiles to enhance celerity-range priors. Our analysis revealed that while the combination of infrasound and seismic observations significantly reduced localization uncertainty, challenges remained, particularly with returns at distances less than 200 km from the source and the influence of specific observations on location estimates. The results indicate that broader celerity distributions, such as those from Blom et al. (2020), facilitate better alignment with ground truth locations compared to narrower models. Overall, this work demonstrates the promise of seismoacoustic approaches in refining event localization and highlights the need for further exploration of celerity-range models to ensure reliable outcomes.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2024

2024 Milestone Report: Site 300 mPDV Optical Fiducials

In the early 2000s Photon Doppler Velocimetry (PDV) replaced the Fabry-Perot many beam system1 and has since become ubiquitous across dynamic experimental platforms to measure velocity and event times such as shock breakout, key variables for high pressure physics research. The advent of optically multiplexed oscilloscope channels to create multiplexed PDV (mPDV) increased portability and reduced price per data point. However, these advantages came with the cost of additional fielding complexity and thermally induced timing drift, which directly affects the ability to use PDV for high precision time measurement. Controlling the temperature or incorporating optical fiducials allows researchers to characterize, reduce, and correct this thermal drift in analysis. Our implemented optical fiducial or “timing marker” allows us to reduce thermal drift uncertainty from the 10s of nanoseconds down to the 100s of picoseconds with minimal added complexity to existing systems. The implementation of this timing marker normalizes uncertainty across optical delays, brings our facility cross timing into the sub nanosecond regime and lets us identify/study anomalies in our data. This improvement increases experiment reliability and quality enabling a new class of high precision experiments at S300.

42 ENGINEERING↗
DOE OSTI2022

SNIa Cosmology Analysis Results from Simulated LSST Images: From Difference Imaging to Constraints on Dark Energy

Abstract The Vera Rubin Observatory Legacy Survey of Space and Time (LSST) is expected to process ∼10 6 transient detections per night. For precision measurements of cosmological parameters and rates, it is critical to understand the detection efficiency, magnitude limits, artifact contamination levels, and biases in the selection and photometry. Here we rigorously test the LSST Difference Image Analysis (DIA) pipeline using simulated images from the Rubin Observatory LSST Dark Energy Science Collaboration Data Challenge (DC2) simulation for the Wide-Fast-Deep survey area. DC2 is the first large-scale (300 deg 2 ) image simulation of a transient survey that includes realistic cadence, variable observing conditions, and CCD image artifacts. We analyze ∼15 deg 2 of DC2 over a 5 yr time span in which artificial point sources from Type Ia supernova (SNIa) light curves have been overlaid onto the images. The magnitude limits per filter are u = 23.66 mag, g = 24.69 mag, r = 24.06 mag, i = 23.45 mag, z = 22.54 mag, and y = 21.62 mag. The artifact contamination levels are ∼90% of all detections, corresponding to ∼1000 artifacts deg –2 in g band, and falling to 300 deg –2 in y band. The photometry has biases <1% for magnitudes 19.5 < m < 23. Our DIA performance on simulated images is similar to that of the Dark Energy Survey difference-imaging pipeline on real images. We also characterize DC2 image properties to produce catalog-level simulations needed for distance bias corrections. We find good agreement between DC2 data and simulations for distributions of signal-to-noise ratio, redshift, and fitted light-curve properties. Applying a realistic SNIa cosmology analysis for redshifts z < 1, we recover the input cosmology parameters to within statistical uncertainties.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Detection of atmospheric rivers with inline uncertainty quantification: TECA-BARD v1.0.1

It has become increasingly common for researchers to utilize methods that identify weather features in climate models. There is an increasing recognition that the uncertainty associated with choice of detection method may affect our scientific understanding. For example, results from the Atmospheric River Tracking Method Intercomparison Project (ARTMIP) indicate that there are a broad range of plausible atmospheric river (AR) detectors and that scientific results can depend on the algorithm used. There are similar examples from the literature on extratropical cyclones and tropical cyclones. It is therefore imperative to develop detection techniques that explicitly quantify the uncertainty associated with the detection of events. We seek to answer the following question: given a “plausible” AR detector, how does uncertainty in the detector quantitatively impact scientific results? We develop a large dataset of global AR counts, manually identified by a set of eight researchers with expertise in atmospheric science, which we use to constrain parameters in a novel AR detection method. We use a Bayesian framework to sample from the set of AR detector parameters that yield AR counts similar to the expert database of AR counts; this yields a set of “plausible” AR detectors from which we can assess quantitative uncertainty. This probabilistic AR detector has been implemented in the Toolkit for Extreme Climate Analysis (TECA), which allows for efficient processing of petabyte-scale datasets. We apply the TECA Bayesian AR Detector, TECA-BARD v1.0.1, to the MERRA-2 reanalysis and show that the sign of the correlation between global AR count and El Niño–Southern Oscillation depends on the set of parameters used.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Evaluation of a high-resolution regional climate simulation for surface and hub-height wind climatology over North America

Assessing the availability of key wind resources requires augmenting observations to support the implementation of wind energy infrastructure. However, observations are limited, necessitating the development of high-resolution, long-term gridded datasets. This study presents a robust, dynamically downscaled climatological dataset, offering 20 years of hourly wind data at a 4 km spatial resolution across North America, and evaluates its performance against observations, including meteorological towers and automated surface-observing system (ASOS) stations, as well as coarse-resolution reanalysis data (the European Centre for Medium-Range Weather Forecasts (ECMWF) reanalysis version 5 (ERA5)). Results demonstrate that the downscaled high-resolution wind data outperform ERA5 in regions of complex terrain and coastal areas, with improved overlap coefficients for wind data distributions and reduced root mean square errors (RMSEs) for hub-height and near-surface diurnal wind patterns. The downscaled simulation also captures the synoptic drivers of seasonal wind direction patterns reasonably well, indicated by high wind rose similarity indices. This study also provides an analysis of interannual variability, utilizing the dataset's full 20-year period, and model uncertainty, generated by varying model initial conditions and physics parameterizations across 1-year ensemble members, which are key considerations for wind resource assessment in wind farm development.

17 WIND ENERGY↗
DOE OSTI2020

Status Update on the High Precision Isotopic Measurements on High Burnup LWR Fuel in 2020

The US Department of Energy (DOE) Office of Nuclear Energy (NE) is currently investigating the feasibility of directly disposing dual-purpose (storage and transportation) canisters (DPCs) in a spent nuclear fuel (SNF) repository. Criticality during the repository performance period (10,000 years or more) is one of the major concerns related to direct disposal of DPCs, specifically as the system undergoes degradation in the repository environment and timeframe. Oak Ridge National Laboratory (ORNL) is developing an as-loaded criticality analysis methodology using full (actinides + fission products) burnup credit that exploits the inherent criticality margin associated with actual canister-specific loading configuration. Burnup credit criticality analysis requires validation of the depletion/decay codes used to generate the burned isotopic inventory of an assembly by comparing the code-predicted inventory with experimentally determined isotopic data. Currently, isotopic measurement data for boiling water reactors (BWR) SNF are limited, and additional measurements will be highly beneficial for BWR burnup credit analysis, which is essential to demonstrate disposability of BWR DPCs. Moreover, new pressurized water reactor (PWR) samples of isotopic measurements will expand the PWR sample population and consequently will greatly improve sample statistics. In turn, this will reduce uncertainty in the computational determination of the isotopic composition of commercial SNF and eliminate additional penalties currently used for lack of data. As part of the isotopic analysis task, eight diverse, high burnup (HBU) samples from several PWR rods are being dissolved and isotopically analyzed to provide high-quality measurement data for PWR SNF to reducing the uncertainties associated with PWR isotopic depletion validation. The isotopics of interest focus on nuclides important to burnup credit, shielding, and decay heat in PWR SNF. Additionally, eight high burnup BWR fuel samples from the Limerick nuclear power plant have been identified for measurement to support reducing uncertainties associated with BWR isotopic depletion validation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Bioproducts from high-protein algal biomass: An economic and environmental sustainability review and risk analysis

This report has High-protein algal biomass as an important bio-commodity that has the potential to provide a new source of sustainable protein products. Herein is a critical review that identifies 1) the most relevant sustainability findings related to the processing of proteinaceous algal biomass to higher value protein products and 2) the potential pathways to improve life cycle assessment (LCA) and techno-economic analysis (TEA) metrics, including life-cycle carbon dioxide equivalent (CO 2 eq), life cycle energy, and minimum selling price (MSP) of these products. The critical review of the literature revealed a large variation in model input parameters relating to these metrics. Therefore, a Monte Carlo analysis was conducted to assess the risk associated with these input variations. To understand the uncertainties that propagate into high-protein algae to products’ systems, we reviewed more than 20 state-of-the-art unit operations for algal biomass processing., including cell disruption, protein solubilization, protein precipitation and purification, and protein concentration. We evaluated displacement of proteinaceous products by algal-bioproducts, including ruminant feed, aquaculture feed, protein tablets, and biopolymers and biopolyesters, with prices in the market ranging from 1.9 to 120 $\$$ kg -1 protein. This review realized that the MSP of ruminant and non-ruminant feed ranges from 0.65±0.56 to 2.9±1.1 $\$$ kg -1 protein, and bioplastics’ MSP ranges from 0.97 to 7.0 $\$$ kg -1 protein. Regarding LCA metrics, there is limited research on life cycle energy in proteinaceous biomass concentration and bioproduct systems, reported at 32.7 MJ kg protein -1 , for animal feed displacement. Animal feed emissions in the literature report negative fluxes, representing environmental benefits, as low as -3.7 kgCO 2 eq kg -1 protein and positive fluxes, i.e., global warming potential, as high as 12.8 kgCO 2 eq kg -1 protein. There is limited research on bioplastics life cycle emissions reported at 0.6 kgCO 2 eq kg -1 protein. In general, the studies to date of algae-derived protein bioproducts showed similar life cycle emissions to soybean meals, nylon, polymers, and polystyrenes. Our risk analysis realized that more than 50% of scenarios can result in negative-net life cycle CO 2 eq emissions. This review and risk analysis assess and demonstrate the scenarios that improve economic and environmental sustainability metrics in high-protein algal bioproducts systems.

09 BIOMASS FUELS↗