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 811 records · Page 45

An efficient hybrid downscaling framework to estimate high-resolution river hydrodynamics

Flow depth and velocity are the most important hydrodynamic variables that govern various river functions, including water resources, navigation, sediment transport, and biogeochemical cycling. Existing high-resolution flow depth simulations rely on either computationally expensive river hydrodynamic models (RHMs) or data-driven models with formidable training costs, whereas data-driven modeling of flow velocity has rarely been explored. Here, using the hybrid Low-fidelity, Spatial analysis, and Gaussian process learning (LSG) model, we developed a downscaling approach to construct high-resolution flow depth and velocity from a two-dimensional (2-D) RHM simulation at coarse resolution. The LSG models were trained and tested in an urban watershed in Houston using two different hurricane-driven flood events. The high-resolution (as fine as 30 m resolution) and low-resolution (mostly 1000 m resolution) meshes include 664 724 and 14 536 grid cells, respectively. The results showed that through downscaling, the simulation errors were reduced to less than one-fourth and one-third of the errors of the low-resolution 2-D RHM for flow depth and velocity, respectively. Our analysis further revealed that the dominant uncertainty sources of the downscaled hydrodynamics are different, with flow velocity dominated by the dimensionality reduction error, which we reduced by using a regionalized training procedure. The downscaling approach achieves an 84-fold acceleration in computational time compared to the high-resolution 2-D RHM, making high-fidelity ensemble flood modeling feasible. More importantly, the developed method provides an opportunity to couple large-scale hydrodynamical processes with local physical, chemical, and biological processes in river models.

Tan, Zeli [Pacific Northwest National Laboratory (↗

Value of Geothermal Energy Storage for Supply-Side and Demand-Side Applications

This report presents the results of a study examining the value potential for geothermal energy storage (GES), a long-duration energy storage resource that stores thermal and/or geomechanical energy in the subsurface. GES could benefit the overall U.S. power system by temporally shifting electricity generation (supply-side) or meeting building heating and cooling load (demand-side). This report analyzes supply-side and demand-side opportunities independently because of differences in applications and models. Currently there is significant uncertainty about the development costs for GES, with only a limited number of demonstration plants for electric energy storage and building heating and cooling storage developments. In this report, we estimate the value of supply-side and demand-side GES to the bulk power system in the contiguous United States. Because of the significant uncertainty about GES development costs, this analysis does not consider GES deployment costs but instead focuses on the value of GES to the U.S. electricity system. The estimated values of GES provide reference points for economically competitive commercial cost targets. Supply-side GES is modeled as part of an enhanced geothermal system (EGS) generation plant in NREL's Regional Energy Deployment System (ReEDS) capacity expansion model (Ho et al. 2021). In contrast to conventional geothermal plants, which generate constant power, EGS plants have unique features that may allow for in-reservoir energy storage for flexible generation. Demand-side GES for heating and cooling, including seasonal hot and cold storage and short-duration heat pump storage, is incorporated into a price-taker model using Cambium electricity marginal cost projections. To establish an upper bound for the value of GES, analysis focused on favorable scenarios for storage with high generation from zero marginal cost, variable renewable energy resources. High penetrations of variable renewable energy generation can increase hourly electricity price variability, which increases the value of temporal energy arbitrage for storage technologies like GES.

15 GEOTHERMAL ENERGY↗

Addressing Uncertainty in MultiSector Dynamics Research

MultiSector Dynamics (MSD) research explores the dynamics and co-evolutionary pathways of human and Earth systems and the emerging interdependent sectors of energy, water, agriculture, and transportation among others. The interactions between these sectors are central to MSD science, as they capture how processes and feedbacks across Earth, environmental, infrastructure, and societal systems shape transitions, socioeconomic risks, and the provision of services. By definition, MSD research requires deep integration across diverse scientific disciplines, ranging from the natural to the social sciences and engineering. All these disciplines apply a variety of numerical simulation models to study and understand their underlying systems of focus. The utility of these models hinges on the fidelity with which they represent the real systems and their ability to produce novel insights about systems and their interactions. The coupled human-natural systems typically represented are shaped by a multitude of interdependent human and natural processes which, when modeled, translate to highly complex, non-linear, interacting behaviors. This ever-increasing complexity massively expands the uncertainty space of a model and can be found in model inputs, processes and parameters. This is further amplified when several models are combined to answer multisectoral questions, as additional uncertainty regarding coupling relationships and interactions is introduced.

IM3, Uncertainty and Sensitivity Analysis, multise↗

Bidding Curve Design for Hybrid Power Plants with Uncertain Solar Forecast: Preprint

This paper presents a novel bidding curve design algorithm tailored for hybrid power plants (HPPs) to participate in the wholesale electricity market. Utilizing forecasts for photovoltaic (PV) generation and available battery power, our algorithm strategically computes the bidding curve to maximize HPP profit while adeptly managing the inherent uncertainty associated with PV power generation. In addition, the introduction of the penalty cost in HPP bidding curves provides the system operator a tool to effectively manage the system-level uncertainty that caused by HPPs. Numerical analysis through Monte Carlo simulations confirms that our bidding curve methodology outperforms the benchmark across various scenarios.

bidding curve↗

Assessing the sensitivity of land-atmosphere coupling strength to boundary and surface layer parameters in the WRF model over Amazon

Modeling tools can be used to diagnose regional land-atmosphere (L-A) coupling strength in the absence of sufficient observations, but subject to uncertainties associated with parameters in model physical parameterizations. Different sensitivity analysis (SA) approaches may lead to different conclusions about the underlying sensitivities. In this study, we quantify simulation uncertainties related to parameter perturbations, and use different approaches to conduct parameter SA on the WRF model pertaining to L-A coupling strength for simulations over the Amazon region. A total of twenty parameters from the Yonsei University (YSU) planetary boundary layer (PBL) and the revised MM5 surface layer (SL) schemes were selected in this analysis. Three different SA methods, the Morris One-at-A-Time (MOAT) method, the Multivariate Adaptive Regression Splines (MARS) method, and the Sobol’ method, were employed to analyze seven WRF-simulated variables and five L-A coupling metrics. Results show that 1) parameter perturbations cause large simulation uncertainties which are comparable to those in the observations; 2) three different SA methods give consistent L-A coupling strength outcomes; 3) six out of the twenty parameters contribute 80%–95% of the total variance in the metrics analyzed, and first-order effects dominate over interaction effects; 4) the twelve variables/metrics of interest show similar sensitivity patterns to the selected parameters, which is consistent across all the methods used. Physical mechanisms for how the sensitive parameters act in determining the L-A coupling strength and associated variables also are illustrated. Our results will help quantifying L-A coupling strength and establishing a basis for parameter calibration over the Amazon region.

Wang, Chen↗

Complete quantification of parametric uncertainties in ( d , p ) transfer reactions

Deuteron-induced transfer reactions are a popular probe in nuclear structure and nuclear astro physics studies. The interpretation of these transfer measurements relies on reaction theory that takes as input effective interactions between the nucleons and the target nucleus. Previous work quantified the uncertainty associated with the optical potentials between the nucleons and the target. In this study, we extend that work by also including the parameters of the mean field associated with the overlap function of the final bound state, thus obtaining the full parametric uncertainty on transfer observables. We use Bayesian Markov Chain Monte Carlo simulations to obtain parameter posterior distributions. We use elastic-scattering cross sections to constrain the optical potential parameters and use the asymptotic normalization coefficient of the final state to constrain the bound state interaction. We then propagate these posteriors to the transfer angular distributions and obtain confidence intervals for this observable. We study (d,p) reactions on 14 C, 16 O, and 48 Ca at energies in the range $E$ $d$ = 7 - 24 MeV. Our results show a strong reduction in uncertainty by using the asymptotic normalization coefficient as a constraint, particularly for those reactions most sensitive to ambiguities in the mean field. For those reactions, the importance of constraining the bound state interaction is equal to that of constraining the optical potentials. Furthermore, the case of 14 C is an outlier because the cross section is less sensitive to the nuclear interior. When minimal constraints are used on the parameters of the nucleon-target interaction, the 1σ uncertainties on the differential cross sections are large (~ 500-100%). However, if elastic-scattering data and the asymptotic normalization coefficient are used in the analysis, with an error of 10% (5%), this uncertainty reduces to ~ 30% (~ 15%).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

The Pantheon+ Analysis: Cosmological Constraints

We present constraints on cosmological parameters from the Pantheon+ analysis of 1701 light curves of 1550 distinct Type Ia supernovae (SNe Ia) ranging in redshift from z = 0.001 to 2.26. This work features an increased sample size from the addition of multiple cross-calibrated photometric systems of SNe covering an increased redshift span, and improved treatments of systematic uncertainties in comparison to the original Pantheon analysis, which together result in a factor of 2 improvement in cosmological constraining power. For a flat ΛCDM model, we find Ω M = 0.334 ± 0.018 from SNe Ia alone. For a flat w 0 CDM model, we measure w 0 = –0.90 ± 0.14 from SNe Ia alone, H 0 = 73.5 ± 1.1 km s –1 Mpc –1 when including the Cepheid host distances and covariance (SH0ES), and w 0 = $-{0.978}_{-0.031}^{+0.024}$ when combining the SN likelihood with Planck constraints from the cosmic microwave background (CMB) and baryon acoustic oscillations (BAO); both w 0 values are consistent with a cosmological constant. We also present the most precise measurements to date on the evolution of dark energy in a flat w 0 w a CDM universe, and measure w a = $-{0.1}_{-2.0}^{+0.9}$ from Pantheon+ SNe Ia alone, H 0 = 73.3 ± 1.1 km s –1 Mpc –1 when including SH0ES Cepheid distances, and w a = $-{0.65}_{-0.32}^{+0.28}$ when combining Pantheon+ SNe Ia with CMB and BAO data. Finally, we find that systematic uncertainties in the use of SNe Ia along the distance ladder comprise less than one-third of the total uncertainty in the measurement of H 0 and cannot explain the present "Hubble tension" between local measurements and early universe predictions from the cosmological model.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Discovering the Multisectoral Impacts of Global Energy Sector Outcomes Through Multiple Ensemble Aggregation Measures

Understanding complex human-Earth system interactions often involves analyzing large scenario ensembles that encompass a wide range of plausible futures. These ensembles often require aggregation to summarize information based on specific criteria or conditions. However, previous research using global change scenario ensembles has largely overlooked how the choice of aggregation method influences the interpretation of results. To address this gap, we leverage a large ensemble data set designed to capture broad energy system dynamics generated using the Global Change Analysis Model. We first explore how energy-related uncertainties are propagated to both global and regional water-energy-food sectors. We then conduct a rank correlation analysis across seven ensemble aggregation measures and demonstrate the need to consider multiple measures in global change scenarios. Our results suggest that global water and food sector outcomes in the 21st century vary widely depending on different scenario assumptions. The global energy productivity is projected to improve by the end of the century across all scenarios. Moreover, regions facing water scarcity challenges in 2100 do not always overlap with those facing extreme energy and food sector outcomes. Although rank correlations across seven aggregation measures are relatively stable across sectors, we identify cases where relying on a single measure leads to losing critical information in the full ensemble. Reliance on a single aggregation measure can distort the interpretation of global change scenario outcomes. Instead, adopting multiple ensemble aggregation measures provides a more holistic understanding of global change scenario ensembles.

Kim, Gijoo↗

Sensitivity of an integrated experiment to uncertainty in the high explosive equations of state

Traditionally, hydrodynamics simulations are performed with a single equation-of-state (EOS) to describe each material. These EOSs typically have a physics-informed functional form with adjustable parameters that are calibrated in order to replicate small-scale data. However, because the calibration data have uncertainty and there are typically inherent degeneracies in fitting the EOS, there are actually multiple EOSs that might be consistent with calibration data. In this work, we perform uncertainty quantification (UQ) for the reactant and product equations of state for the high explosive PBX 9501 to yield an ensemble of EOSs that match the uncertain small-scale calibration data. We then simulate an experiment of an explosively formed penetrator repeatedly with different EOSs to both validate the UQ analysis and determine the effects of EOS uncertainty on the prediction of quantities of interest in the experiment. In general, we find good agreement between the simulation predictions and the experimental measurements, and we identify an EOS variable that contributes most directly to the spread in the predictions as the EOSs are varied.

36 MATERIALS SCIENCE↗

Processing Aleatory and Epistemic Uncertainties in Experimental Data From Sparse Replicate Tests of Stochastic Systems for Real-Space Model Validation

This paper presents a practical methodology for propagating and processing uncertainties associated with random measurement and estimation errors (that vary from test-to-test) and systematic measurement and estimation errors (uncertain but similar from test-to-test) in inputs and outputs of replicate tests to characterize response variability of stochastically varying test units. Also treated are test condition control variability from test-to-test and sampling uncertainty due to limited numbers of replicate tests. These aleatory variabilities and epistemic uncertainties result in uncertainty on computed statistics of output response quantities. The methodology was developed in the context of processing experimental data for “real-space” (RS) model validation comparisons against model-predicted statistics and uncertainty thereof. The methodology is flexible and sufficient for many types of experimental and data uncertainty, offering the most extensive data uncertainty quantification (UQ) treatment of any model validation method the authors are aware of. It handles both interval and probabilistic uncertainty descriptions and can be performed with relatively little computational cost through use of simple and effective dimension- and order-adaptive polynomial response surfaces in a Monte Carlo (MC) uncertainty propagation approach. A key feature of the progressively upgraded response surfaces is that they enable estimation of propagation error contributed by the surrogate model. Sensitivity analysis of the relative contributions of the various uncertainty sources to the total uncertainty of statistical estimates is also presented. Finally, the methodologies are demonstrated on real experimental validation data involving all the mentioned sources and types of error and uncertainty in five replicate tests of pressure vessels heated and pressurized to failure. Simple spreadsheet procedures are used for all processing operations.

97 MATHEMATICS AND COMPUTING↗

Stability of Cool Cores during Galaxy Cluster Growth: A Joint Chandra/SPT Analysis of 67 Galaxy Clusters along a Common Evolutionary Track Spanning 9 Gyr

We present the results of a joint analysis of Chandra X-ray and South Pole Telescope (SPT) Sunyaev–Zel’dovich (SZ) observations targeting the first sample of galaxy clusters at 0.3 < z < 1.3, selected to be the progenitors of well-studied nearby clusters based on their expected accretion rate. We develop a new procedure in order to tackle the analysis challenge that is estimating the intracluster medium (ICM) properties of low-mass and high-redshift clusters with ~150 X-ray counts. One of the dominant sources of uncertainty on the ICM density profile estimated with a standard X-ray analysis with such shallow X-ray data is due to the systematic uncertainty associated with the ICM temperature obtained through the analysis of the background-dominated X-ray spectrum. We show that we can decrease the uncertainty on the density profile by a factor varying between 2 and 8 with a joint deprojection of the X-ray surface brightness profile measured by Chandra and the SZ-integrated Compton parameter available in the SPT cluster catalog. We apply this technique to the whole sample of 67 clusters in order to track the evolution of the ICM core density during cluster growth. We confirm that the evolution of the gas density profile is well modeled by the combination of a fixed core and a self-similarly evolving non-cool-core profile. We show that the fraction of cool cores in this sample is remarkably stable with redshift although clusters have gained a factor of ~4 in total mass over the past ~9 Gyr.

79 ASTRONOMY AND ASTROPHYSICS↗

Field Measurement and Analysis of Frequency and RoCoF for Low-Inertia Power Systems

The frequency measurement and rate of change of frequency (RoCoF) calculations are important for active power, frequency control, and frequency-related protection in power system operations. With the rapid decrease in inertia under modern power systems, however, multiple power variations introduced by variable renewable energy bring large RoCoF and introduce oscillations in frequency measurements. To shed light on this issue, this article proposes a frequency measurement and RoCoF analysis technology, as well as the sensitivity analysis for the real low-inertia power grid. First, a real-world event in an island power grid is analyzed as an example, and an event numerical model is estimated considering four uncertainties, including event RoCoF, oscillation magnitude, oscillation frequency, and underfrequency load-shedding threshold. Then, a sensitivity analysis is performed using the Sobol indices method and quasi-Monte Carlo simulations to identify the dominant uncertainties based on the event model. Here, the dominant uncertainties are compared to provide a guideline for the frequency estimation algorithm selection under different event cases. This algorithm is compared with an enhanced zero-crossing algorithm and gives a 3.41% better RoCoF error rate. Finally, an experiment is conducted using universal grid analyzers to validate the simulation results and the theoretical analysis.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Measurement of Helium and Xenon Knudsen permeability in graphite IG-11

Here, understanding the flow and diffusion (transport) of gases through nuclear graphite is of interest in both the prismatic and pebble bed high-temperature gas reactors (HTGRs), for normal operation and under accident conditions, for example, as related to graphite oxidation. Both the laminar and turbulent gaseous transport are of interest. It has also been noted that nuclear graphite can have pore sizes that can possibly lead to non-continuum transport, as the pore size can be comparable to the gaseous mean free path. We report measurements and analyzes of experiments performed with IG-11 graphite in the pressure range of 10 Pa to 14,700 Pa, 293 K, using Helium, and research grade xenon (99.999%), respectively. These values are in fair agreement with those reported for different graphites and gases by other investigators, but are also higher by six orders of magnitude than those reported for Kr permeability in a developmental very low permeability graphite-HXT-90. Our measured ratio of the Helium and Xenon permeabilities is about 4, different from $\sqrt{\frac{m_{Xe}}{m_{He}}}$ ≈ 5.72 that corresponds to diffuse reflection at the graphite surface for both gases. Simple analysis indicates that the two gases have slightly different accommodation coefficients (about 25%) with the graphite surface. However, there could also be other reasons, such as experimental uncertainties or the need for a more detailed analysis.

36 MATERIALS SCIENCE↗

Round Robin Measurements of Molten Salt Properties for LiF-NaF-KF (FLiNaK) and NaCl-KCl Mixtures

The development, operation, and regulation of nuclear reactors that utilize molten salts as fuel or as heat transfer media require knowledge of the thermal properties of the salt systems and quantification of the corresponding uncertainties. Knowledge of molten salt properties is also necessary for applications in material synthesis, processing, separations, solar thermal power generation, and energy storage. A round robin was conducted with national laboratory and university participants from twenty-one laboratories in five countries to compare property measurements, to better understand uncertainties, and to identify possible best practices. Two salt mixtures, each from a common batch, were distributed to participants for evaluation: equimolar NaCl-KCl and 45.0LiF-13.7NaF-41.3KF mol % (FLiNaK). Measurements were performed to determine the major constituent composition, oxygen content, density, thermal expansivity, melting point, and thermal conductivity. Error analysis was performed on each measurement for uncertainty quantification for each type of property that was explored. Finally, the resulting discussion of the methodologies used in this work is meant to lay the groundwork for the development of standard methods and reference materials for future high-temperature property measurements on halide melts.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Uncertainty Quantification Framework for Predicting Material Response with Large Number of Parameters: Application to Creep Prediction in Ferritic-Martensitic Steels Using Combined Crystal Plasticity and Grain Boundary Models

This paper presents an uncertainty quantification (UQ) framework for the physics-based model prediction of material response with a large number of parameters. The application problem presented in this work is that of predicting creep in Grade 91 steel at 600°C. The material response is defined with a physically based microstructural model with constitutive equations emulating several observed phenomena in Grade 91 and embodied into an explicit geometry mesoscale finite element model for prior austenite grains and grain boundaries. Creep within the grains and in grain boundaries are represented by crystal plasticity for dislocation motion and a physics-based model for cavity growth and nucleation, respectively. The creep behavior of this material is influenced by several parameters, some of which have a wide range of variation based on experimental data. UQ combined with microstructural modeling can discover the core microstructural causes of experimental variability, leading to improved materials with lower variability in critical long-term material properties. In this study, we investigate the model's uncertainty to identify material properties that may be modified during production to increase creep life and analyze different components of the crystal plasticity model for improvements. For this purpose, a quantity of interest is defined as time to minimum creep rate, which correlates well to the creep failure of the material. A deep neural network model was trained and validated to be used as a surrogate for the finite element model. Then, a variance-based sensitivity analysis is performed on the surrogate model to find the Sobol indices of the input parameters in respect to the output quantity of interest. The Sobol indices are used to reduce the dimensionality of the model. Generalized polynomial chaos expansion is used on the reduced basis models to propagate the uncertainty from the input parameters to the quantity of interest using the deep neural network surrogate model. These results are benchmarked against uncertainty propagation using Monte Carlo simulations. In conclusion, the UQ performed through the reduced basis model captures almost all the uncertainty in the model with significantly fewer simulations, making it possible to perform the UQ directly via simulations with the finite element model rather than surrogate machine-learned models.

36 MATERIALS SCIENCE↗

Improving helium abundance determinations with Leo P as a case study

Currently, the primordial helium abundance is best estimated through spectroscopic observations of H II regions in metal-poor galaxies. However these determinations are limited by several systematic uncertainties which ultimately limit our ability to accurately ascertain the primordial abundance. In this study, we improve the methodologies for solving for the reddening, the emission contributions from collisional excitation of the H I atoms, the effects underlying absorption in the H I and He I emission lines, and the treatment of the blended H I and He I emission at λ3889 with the aim of lowering the systematic uncertainties in helium abundance determinations. To apply these methods, we have obtained observations of the He I λ10830 emission line in the brightest H II region in the extremely metal-poor (3% Z ⊙ ) galaxy Leo P with the LUCI1 instrument on the LBT. We combine this measurement with previous MODS/LBT observations to derive an improved helium abundance. In doing so, our present analysis results in a decrease in the uncertainty in the helium abundance of Leo P by approximately 70%. Furthermore, this result is combined with data from other observations to estimate the primordial helium mass fraction, Yp= 0.2453 ± 0.0034.

79 ASTRONOMY AND ASTROPHYSICS↗

Fast Computational Algorithms for Partial Differential Equations and Uncertainty Quantifications

This project concerned the construction, testing and analysis of computational algorithms for solving parameterized and stochastic partial differential equations. The study and understanding of equations of this type is of fundamental importance in numerous engineering and scientific applications. Examples include simulation of plasma dynamics in models of electric propulsion and nuclear fusion, simulation of multiphase flows, such as the flow of water, gas and oil in reservoirs, and structural analysis of the dependence of structures on materials. Parametrization is used in such settings when properties of the models such as viscosity of fluids or electric resistivity of materials are not precisely understood and instead are treated as random variables. The resulting solutions are themselves random, and having such solutions will enable engineers to use probabilistic methods to assess the likelihood of events, for example, whether a pollutant in a liquid will exceed a limit, and to use such analyses to develop ways to ensure positive outcomes. Construction of accurate (high resolution) computational solutions is expensive, requiring significant computer time and computational resources, and there is need to reduce computational cost to make simulation useful and effective. The aim of the project was to construct algorithms to efficiently compute surrogate solutions to parameterized problems to allow for efficient and accurate simulation. The technical approach used focused on two related strategies, based on rank-reduction methods and reduced-order models. These methods construct surrogate solutions of parameter-dependent models by projection or interpolation into low-dimensional approximation spaces. Cost savings are achieved if the low-dimensional spaces can be identified and constructed efficiently and if the resulting low-dimensional algebraic systems can be solved cheaply. Accomplishments include: Theoretical and empirical demonstration of the effectiveness of fast multigrid solution strategies for computing low-rank representations of parameter-dependent solutions to discrete partial differential equations, including the first proof establishing so-called textbook convergence properties for low-rank methods. Development of efficient solution algorithms for solving nonlinear parameter-dependent partial differential equations used in models of fluid dynamics. Developent of efficient algorithms for low-rank representation of solutions of time-dependent simulations of fluid dynamics using multi-dimensional tensor representations of solutions.

97 MATHEMATICS AND COMPUTING↗