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,585 records · Page 88

Estimates of Southern Hemispheric Gravity Wave Momentum Fluxes across Observations, Reanalyses, and Kilometer-Scale Numerical Weather Prediction Model

Abstract Gravity waves (GWs) are among the key drivers of the meridional overturning circulation in the mesosphere and upper stratosphere. Their representation in climate models suffers from insufficient resolution and limited observational constraints on their parameterizations. This obscures assessments of middle atmospheric circulation changes in a changing climate. This study presents a comprehensive analysis of stratospheric GW activity above and downstream of the Andes from 1 to 15 August 2019, with special focus on GW representation ranging from an unprecedented kilometer-scale global forecast model (1.4 km ECMWF IFS), ground-based Rayleigh lidar (CORAL) observations, modern reanalysis (ERA5), to a coarse-resolution climate model (EMAC). Resolved vertical flux of zonal GW momentum (GWMF) is found to be stronger by a factor of at least 2–2.5 in IFS compared to ERA5. Compared to resolved GWMF in IFS, parameterizations in ERA5 and EMAC continue to inaccurately generate excessive GWMF poleward of 60°S, yielding prominent differences between resolved and parameterized GWMFs. A like-to-like validation of GW profiles in IFS and ERA5 reveals similar wave structures. Still, even at ∼1 km resolution, the resolved waves in IFS are weaker than those observed by lidar. Further, GWMF estimates across datasets reveal that temperature-based proxies, based on midfrequency approximations for linear GWs, overestimate GWMF due to simplifications and uncertainties in GW wavelength estimation from data. Overall, the analysis provides GWMF benchmarks for parameterization validation and calls for three-dimensional GW parameterizations, better upper-boundary treatment, and vertical resolution increases commensurate with increases in horizontal resolution in models, for a more realistic GW analysis. Significance Statement Gravity wave–induced momentum forcing forms a key component of the middle atmospheric circulation. However, complete knowledge of gravity waves, their atmospheric effects, and their long-term trends are obscured due to limited global observations, and the inability of current climate models to fully resolve them. This study combines a kilometer-scale forecast model, modern reanalysis, and a coarse-resolution climate model to first compare the resolved and parameterized momentum fluxes by gravity waves generated over the Andes, and then evaluate the fluxes using a state-of-the-art ground-based Rayleigh lidar. Our analysis reveals shortcomings in current model parameterizations of gravity waves in the middle atmosphere and highlights the sensitivity of the estimated flux to the formulation used.

Meteorology & Atmospheric Sciences↗

Uncertainty evaluation in reference material characterization by double isotope dilution inductively coupled plasma mass spectrometry (ID-ICP-MS)

Double isotope dilution inductively coupled plasma mass spectrometry (ID-ICP-MS) is one of the most important measurement methods in the establishment of measurement standards for inorganic analysis. However, there is still no good consensus on approaches for uncertainty evaluation in double ID-ICP-MS. Particularly, approaches to treating the instrumental isotopic fractionation (IIF) phenomenon and correlation between uncertainty sources are diverse and often unclear in the literature. Here, an approach for uncertainty evaluation in double ID-ICP-MS measurements is described. Using three certified reference materials (CRMs), it is shown that the IIF effects cancel out and should not be included in the uncertainty budget. Further, a consistent way of considering correlated sources of uncertainty to prevent double counting of uncertainties is described in detail along with an example of measuring the mass fraction of magnesium in an oyster powder CRM.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Solar energy system economic evaluation: IBM System 2, Togus, Maine

The economic analysis of the solar energy system, is developed for Torgus and four other sites typical of a wide range of environmental and economic conditions in the continental United States. This analysis is accomplished based on the technical and economic models in the f-chart design procedure with inputs taken on the characteristics of the installed system and local conditions. The results are expressed in terms of the economic parameters of present worth of system cost over a projected twenty year life, life cycle savings, year of positive savings and year of payback for the optimized solar energy system at each of the analysis sites. The sensitivity of the economic evaluation to uncertainties in constituent system and economic variables is also investigated. Results demonstrate that the solar energy system is economically viable at all of the five sites for which the analysis was conducted.

Source record↗

Physics-based Entry, Descent and Landing Risk Model

A physics-based risk model was developed to assess the risk associated with thermal protection system failures during the entry, descent and landing phase of a manned spacecraft mission. In the model, entry trajectories were computed using a three-degree-of-freedom trajectory tool, the aerothermodynamic heating environment was computed using an engineering-level computational tool and the thermal response of the TPS material was modeled using a one-dimensional thermal response tool. The model was capable of modeling the effect of micrometeoroid and orbital debris impact damage on the TPS thermal response. A Monte Carlo analysis was used to determine the effects of uncertainties in the vehicle state at Entry Interface, aerothermodynamic heating and material properties on the performance of the TPS design. The failure criterion was set as a temperature limit at the bondline between the TPS and the underlying structure. Both direct computation and response surface approaches were used to compute the risk. The model was applied to a generic manned space capsule design. The effect of material property uncertainty and MMOD damage on risk of failure were analyzed. A comparison of the direct computation and response surface approach was undertaken.

descent and landing↗

Solar energy system economic evaluation final report for SEMCO-Loxahatchee, Loxahatchee National Wildlife refuge, Palm Beach County, Florida

Economic analysis of the solar energy system installed at Loxahatchee, was developed for Loxahatchee and four other sites typical of a wide range of environmental and economic conditions in the continental United States. This analysis was accomplished based on the technical and economic models in the f Chart design procedure with inputs based on the characteristics of the installed system and local conditions. The results are expressed in terms of the economic parameters of present worth of system costs over a projected twenty year life, life cycle savings, year of positive savings and year of payback for the optimized solar energy system at each of the analysis sites. The sensitivity of the economic evaluation to uncertainties in constituent system and economic variables was also investigated. The results demonstrate that the solar energy system is economically viable at all of the five sites for which the analysis was conducted.

Source record↗

New Perspectives on Spectroscopic Factor Quenching from Reactions

The evolution of single-particle strengths as the neutron-to-proton asymmetry changes informs us of the importance of short- and long-range correlations in nuclei and has therefore been extensively studied for the last two decades. Surprisingly, the strong asymmetry dependence of these strengths and their extreme values for highly asymmetric nuclei inferred from knockout reaction measurements on a target nucleus are not consistent with what is extracted from electron-induced, transfer, and quasi-free reaction data, constituting a two-decade old puzzle. This work presents the first consistent analysis of one-nucleon transfer and one-nucleon knockout data, in which theoretical uncertainties associated with the nucleon-nucleus effective interactions considered in the reaction models are quantified using a Bayesian analysis. Here, our results demonstrate that, taking into account these uncertainties, the spectroscopic strengths of loosely bound nucleons extracted from both probes agree with each other and, although there are still discrepancies for deeply bound nucleons, the slope of the asymmetry dependence of the single-particle strengths inferred from transfer and knockout reactions are consistent within 1 σ . Both probes are consistent with a small asymmetry dependence of these strengths. The uncertainties obtained in this work represent a lower bound and are already significantly larger than the original estimates.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Uncertainty reduction in residual stress measurements by an optimised inverse solution using nonconsecutive polynomials

Many destructive methods for measuring residual stresses such as the slitting method require an inverse analysis to solve the problem. The accuracy of the result as well as an uncertainty component (the model uncertainty) depends on the basis functions used in the inverse solution. The use of a series expansion as the basis functions for the inverse solution was analysed in a previous work for the particular case where functions orders grew consecutively. The present work presents a new estimation of the model uncertainty and a new improved methodology to select the final basis functions for the case where the basis is composed of polynomials. Including nonconsecutive polynomial orders in the basis generates a larger space of possible solutions to be evaluated and allows the possibility to include higher-order polynomials. The paper includes a comparison with two other inverse analyses methodologies applied to synthetically generated data. With the new methodology, the final error is reduced and the uncertainty estimation improved.

36 MATERIALS SCIENCE↗

Probabilistic Predictions for Fastener Failure in the Sandia Mechanics Challenge Using the Discrete-Direct Uncertainty Quantification Approach

This paper documents the blind and post-blind analysis predictions for the 2023 Sandia Mechanics Challenge (SMC), which involved predicting the behavior of a threaded fastener joint structure subjected to shock loading. Utilizing repeat sets of fastener calibration data from various experimental configurations including tension, double shear, and joint tension, we developed a library of calibrated models which were propagated through the application model using the Discrete-Direct (DD) uncertainty quantification (UQ) approach. Although the initial blind predictions did not incorporate spare-sample processing to quantify fastener failure probabilities, the analyses yielded reasonable conclusions aligned with experimental results. In the post-blind analysis phase, we focused on enhancing the fidelity of the aluminum constitutive model and innovating the DD approach to obtain probabilistic predictions for fastener failure, particularly when quantities of interest (QoIs) approach their bounds. The improved aluminum model captures the behavior of the cantilever under shock loading more accurately, predicting both partial and complete cracks, although it tends to underpredict failure propagation. The enhanced DD approach facilitates probabilistic predictions that reflect the interdependent failure mechanisms of the fasteners and the cantilever, revealing that while certain fasteners are more likely to fail, the failure does not necessarily follow a progressive pattern. Overall, the post-blind analyses significantly improved the predictive capabilities of the model, providing valuable insights into the SMC application and establishing a robust foundation for informed engineering decisions. The methodology demonstrates a cost-effective and extensible approach suitable for a wide range of applications, highlighting the importance of uncertainty quantification to provide context for engineering decision making.

42 ENGINEERING↗

Machine Learning Driven Sensitivity Analysis of E3SM Land Model Parameters for Wetland Methane Emissions

Methane (CH 4 ) is globally the second most critical greenhouse gas after carbon dioxide, contributing to 16%–25% of the observed atmospheric warming. Wetlands are the primary natural source of methane emissions globally. However, wetland methane emission estimates from biogeochemistry models contain considerable uncertainty. One of the main sources of this uncertainty arises from the numerous uncertain model parameters within various physical, biological, and chemical processes that influence methane production, oxidation, and transport. Sensitivity Analysis (SA) can help identify critical parameters for methane emission and achieve reduced biases and uncertainties in future projections. This study performs SA for 19 selected parameters responsible for critical biogeochemical processes in the methane module of the Energy Exascale Earth System Model (E3SM) land model (ELM). The impact of these parameters on various CH 4 fluxes is examined at 14 FLUXNET- CH 4 sites with diverse vegetation types. Given the extensive number of model simulations needed for global variance-based SA, we employ a machine learning (ML) algorithm to emulate the complex behavior of ELM methane biogeochemistry. We found that parameters linked to CH 4 production and diffusion generally present the highest sensitivities despite apparent seasonal variation. Comparing simulated emissions from perturbed parameter sets against FLUXNET-CH 4 observations revealed that better performances can be achieved at each site compared to the default parameter values. This presents a scope for further improving simulated emissions using parameter calibration with advanced optimization techniques.

54 ENVIRONMENTAL SCIENCES↗

Efficient Probabilistic Visualization of Local Divergence of 2D Vector Fields with Independent Gaussian Uncertainty

This work focuses on visualizing uncertainty of local divergence of two-dimensional vector fields. Divergence is one of the fundamental attributes of fluid flows, as it can help domain scientists analyze potential positions of sources (positive divergence) and sinks (negative divergence) in the flow. However, uncertainty inherent in vector field data can lead to erroneous divergence computations, adversely impacting downstream analysis. While Monte Carlo (MC) sampling is a classical approach for estimating divergence uncertainty, it suffers from slow convergence and poor scalability with increasing data size and sample counts. Thus, we present a two-fold contribution that tackles the challenges of slow convergence and limited scalability of the MC approach. (1) We derive a closed-form approach for highly efficient and accurate uncertainty visualization of local divergence, assuming independently Gaussian-distributed vector uncertainties. (2) We further integrate our approach into Viskores, a platform-portable parallel library, to accelerate uncertainty visualization. In our results, we demonstrate significantly enhanced efficiency and accuracy of our serial analytical (speed-up up to 1946×) and parallel Viskores (speed-up up to 19698×) algorithms over the classical serial MC approach. We also demonstrate qualitative improvements of our probabilistic divergence visualizations over traditional mean-field visualization, which disregards uncertainty. We validate the accuracy and efficiency of our methods on wind forecast and ocean simulation datasets.

Ouermi, Timbwaoga [University of Utah]↗

Reaction History Uncertainty Propagation from Flux to α and from Time to α [Slides]

We discuss uncertainty propagation in reaction history data from underground nuclear tests. γ reaction history detectors measured flux as a function of time. After data analysis we report α as a function of time. Flux uncertainty and time uncertainty propagation into α uncertainty are discussed in detail. Discretization of the formulas is discussed.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

The right conditions for high-precision dynamic temperature and heat capacity measurement via pyrometry and conductivity

The pursuit of accurate bulk temperature T under extreme conditions has been a long-standing goal of the high pressure science community, complicated by a lack of data to inform models. To reach these extremely high-pressure, high-temperature (high P − T) conditions, a combination of dynamic and heated static experiments (e.g., diamond or gem anvil cel experiments) are used. For example, in a diamond anvil cell (DAC) experiment, a sample placed in the DAC is first pressurized. Following pressurization, the sample T is increased either by heating the entire DAC (usually using resistive heating, and limited to ∼1000K) or by applying intense laser power to the sample surfaces. In a dynamic experiment, the process of pressurizing the sample also heats it. In the case of shock physics experiments, such heating is substantial, easily reaching thousands of Kelvin; in our work we have seen T ∼17000K. Most methods of measuring temperature at ambient are not compatible with experiments under these high-pressure, high-temperature conditions: thermocouples break, melt, or have conductivity properties that differ from ambient where they are calibrated; thermometers would melt; both are too slow. As a result most methods are based on non-contact techniques such as x-ray diffraction broadening, neutron scattering, or optical methods. Of these, optical methods using the visible and near-infrared region of the spectrum are the most commonly used as the sources and detectors are readily available. In the case of optical methods the optical depth, and therefore the measurement location, is limited to the surface. When a window or anvil material is used, heat flows from the sample into the window/anvil. Likewise, if the sample undergoes a change in thermodynamic state, such as expansion upon release, different T may be expected. As a result, the surface or apparent temperature T app measurement will differ from the bulk or interior temperature that is desired. This surface measurement must be related to the bulk measurement using thermal transport models and material models. While it is tempting to conclude that one should just use x-ray methods that directly probe the interior, even these methods have been shown to depend on thermal transport and material models. Regardless of the method used to create the high P − T condition, therefore, we must understand the role of thermal transport and material models upon our interpretation of the T measurement, as well as the errors and uncertainties associated with the choice of models used in the analysis. This is a substantial area of research and this paper is by no means a complete survey of the relevant sources of uncertainty. For example, we have yet to begin to address alternate transport models in a detailed manner (e.g., Tan-Ahrens), or the many models that use additional layers to approximate melting, turbulence, or epitaxial phenomena). Likewise, we have not explored the impact upon uncertainty of thermal models that use temperature-dependent thermal transport coefficients, or the wide range of material models that can be applied. Instead, this paper focuses on using one simple model, the Urtiew-Grover model, to understand the sources of error in T measurement so that we may identify how best to focus future research efforts to return the best improvements and avoid working on over-optimizing a single type of measurement. To this end, we work through some of the best and worst case scenarios for T measurement.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Benchmarking FFTF LOFWOS Test# 13 using SAM code: Baseline model development and uncertainty quantification

The development and deployment of advanced reactors, such as the sodium-cooled fast reactor (SFR), relies on sophisticated modeling tools to ensure the safety of the design under various transients. The predictive capability of these advanced modeling tools requires validation to garner trust in supporting the licensing of the advanced reactors. For this reason, the International Atomic Energy Agency (IAEA) initiated a coordinated research project (CRP) in 2018 for the analysis of the Fast Flux Test Facility (FFTF) Loss of Flow Without Scram (LOFWOS) Test #13.In this study, we present and discuss the benchmarking efforts of the modern system code SAM on the FFTF LOFWOS Test #13. Further, the SAM baseline model was developed according to the benchmark specification, which included a detailed core model with reactivity feedback. Generally, good agreement was observed between the baseline results and benchmark measurements; however, discrepancies persisted, particularly in predicted fuel assembly coolant outlet temperatures. Utilizing the baseline model, uncertainty quantification (UQ) and sensitivity analysis (SA) were conducted with the assistance of various statistical learning and machine learning methods, including kernel density estimation, Gaussian processes, and Sobol indices. Following the baseline model prediction and UQ and SA results, we discuss the reasons for the simulation discrepancies and propose further improvements to the model. This benchmarking effort adheres to the best-estimate plus uncertainty approach and can serve as a valuable example for supporting risk-informed licensing of advanced reactors.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

INL Poster - Juan Barrera Salazar

Generation IV nuclear reactors introduce several advantages and benefits in terms of safety and efficiency when compared with their predecessors from previous generation. This is due, among many things, to the use of innovative forms of fuel and coolant, different from the conventional ones used in the last decades. Given that these upcoming designs utilize emerging technologies, the related instrumentation is also in the process of being developed; therefore, it is necessary to establish the sensitivity requirements and the effects of uncertainty on different properties of the components and elements of the reactor designs. This report presents the results of simulations that quantify the impacts of the uncertainties of four thermophysical properties of the refrigerant salt (LiF-BeF2) for the Kairos Power benchmark model (g-FHR) for steady state making use of the Sobol’ method through polynomial chaos surrogate modeling. The properties of the salt to which uncertainty was evaluated were density, dynamic viscosity, thermal conductivity and heat capacity. This study was carried out using the Griffin/Pronghorn multiphysics model under the computational resources of the Idaho National Laboratory (INL) High Performance Computing (HPC). The results indicate a weak dependence of the uncertainty of thermal conductivity on the quantities of core pressure drop and core outlet temperature.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗

High pressure/high temperature multiphase simulations of dodecane injection to nitrogen: Application on ECN Spray-A

The present work investigates the complex phenomena associated with pressure/high temperature dodecane injection for the Engine Combustion Network (ECN) Spray-A case, employing more elaborate thermodynamic closures, to avoid well known deficiencies concerning density and speed of sound prediction using traditional cubic models. A tabulated thermodynamic approach is proposed here, based on log 10 (p)-T tables, providing very high accuracy across a large range of pressures, spanning from 0 to 2500 bar, with only a small number of interpolation points. The tabulation approach is directly extensible to any thermodynamic model, existing or to be developed in the future. Here NIST REFPROP properties are used, combined with PC-SAFT Vapor-Liquid-Equilibrium to identify the liquid in mixtures penetration, hence avoiding the use of an arbitrary threshold for mass fraction. Identified liquid and vapour penetration are compared against experimental data from the ECN database showing a good agreement, within approximately 3–8% for axial penetration of liquid, 2% for vapor axial penetration and within experimental uncertainty for radial distribution of mass fraction. Analysis of the vortex evolution indicates that driving mechanisms behind the jet break-up are vortex tilting/stretching, then baroclinic torque, leading to Rayleigh-Taylor instabilities, closely followed by vortex dilation and finally viscous effects.

42 ENGINEERING↗

A novel framework for hosting capacity analysis with spatio-temporal probabilistic voltage sensitivity analysis

Smart grids are envisioned to accommodate high penetration of distributed photovoltaic (PV) generation, which may cause adverse grid impacts in terms of voltage violations. Therefore, PV Hosting capacity is being used as a planning tool to determine the maximum PV installation capacity that causes the first voltage violation and above which would require infrastructure upgrades. Additionally, traditional methods of Hosting capacity analysis are scenario based and computationally complex as they rely on iterative load flow algorithms that require investigating a large number of scenarios for accurate assessment of PV impacts. Therefore, this paper presents a computationally efficient analytical approach to compute the probability distribution of voltage change due to random behavior of randomly located multiple distributed PVs. The proposed approach is based on Spatio-temporal probabilistic voltage sensitivity analysis that exploits both spatial and temporal uncertainties associated with PV injections. Thereafter, the derived distribution is used to quantify voltage violations for various PV penetration levels and subsequently determine the hosting capacity of the system without the need to examine large number of scenarios. Results of the proposed framework are validated via conventional load flow based simulation approach on the IEEE 37 and IEEE 123 node test systems.

42 ENGINEERING↗

The gallium anomaly

In order to test the end-to-end operations of gallium solar neutrino experiments, intense electron-capture sources were fabricated to measure the responses of the radiochemical SAGE and GALLEX/GNO detectors to known fluxes of low-energy neutrinos. Such tests were viewed at the time as a cross-check, given the many tests of 71 Ge recovery and counting that had been routinely performed, with excellent results. However, the four 51 Cr and 37 Ar source experiments yielded rates below expectations, a result commonly known as the Ga anomaly. As the intensity of the electron-capture sources can be measured to high precision, the neutrino lines they produce are fixed by known atomic and nuclear rates, and the neutrino absorption cross section on 71 Ga is tightly constrained by the lifetime of 71 Ge, no simple explanation for the anomaly has been found. To check these calibration experiments, a dedicated experiment BEST was performed, utilizing a neutrino source of unprecedented intensity and a detector optimized to increase statistics while providing some information on counting rate as a function of distance from the source. The results BEST obtained are consistent with the earlier solar neutrino calibration experiments, and when combined with those measurements, yield a Ga anomaly with a significance of approximately 4𝜎, under conservative assumptions. But BEST found no evidence of distance dependence and thus no explicit indication of new physics. In this review we describe the extensive campaigns carried out by SAGE, GALLEX/GNO, and BEST to demonstrate the reliability and precision of their experimental procedures, including 71 Ge recovery, counting, and analysis. We also describe efforts to define uncertainties in the neutrino capture cross section, which now include estimates of effects at the ≲ 0.5% level such as radiative corrections and weak magnetism. With the results from BEST, an anomaly remains even if one retains only the transition to the 71 Ge ground state, whose strength is fixed by the known lifetime of 71 Ge. We then consider the new-physics solution most commonly suggested to resolve the Ga anomaly, oscillations into a sterile fourth neutrino, 𝜈 𝑒 → 𝜈 𝑠 . We find such a solution generates substantial tension with several null experiments, owing to the large mixing angle required. While this does not exclude such solutions – the sterile sector might include multiple neutrinos as well as new interactions – it shows the need for more experimental constraints, if we are to make progress in resolving the Ga and other low-energy neutrino anomalies. We conclude by consider the role future low-energy electron-capture sources could play in this effort.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Gridded daily weather data for North America with comprehensive uncertainty quantification

Abstract Access to daily high-resolution gridded surface weather data based on direct observations and over long time periods is essential for many studies and applications including vegetation, wildlife, soil health, hydrological modelling, and as driver data in Earth system models. We present Daymet V4, a 40-year daily meteorological dataset on a 1 km grid for North America, Hawaii, and Puerto Rico, providing temperature, precipitation, shortwave radiation, vapor pressure, snow water equivalent, and day length. The dataset includes an objective quantification of uncertainty based on strict cross-validation analysis for temperature and precipitation results. The dataset represents several improvements from a previous version, and this data descriptor provides complete documentation for updated methods. Improvements include: reductions in the timing bias of input reporting weather station measurements; improvement to the three-dimensional regression model techniques in the core algorithm; and a novel approach to handling high elevation temperature measurement biases. We show cross-validation analyses with the underlying weather station data to demonstrate the technical validity of new dataset generation methods, and to quantify improved accuracy.

54 ENVIRONMENTAL SCIENCES↗