SEARCH · Search NASA

Results for “uncertainity”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 343 records · Page 19

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Fuel performance uncertainty quantification and sensitivity analysis in the presence of epistemic and aleatoric sources of uncertainties

Fuel performance modeling and simulation includes many uncertain parameters from models to boundary conditions, manufacturing parameters and material properties. These parameters exhibit large uncertainties and can have an epistemic or aleatoric nature, something that renders fuel performance code-to-code and code-to-measurements comparisons for complex phenomena such as the pellet cladding mechanical interaction (PCMI) very challenging. Additionally, PCMI and other complex phenomena found in fuel performance modeling and simulation induce strong discontinuities and non-linearities that can render difficult to extract meaningful conclusions form uncertainty quantification (UQ) and sensitivity analysis (SA) studies. In this work, we develop and apply a consistent treatment of epistemic and aleatoric uncertainties for both UQ and SA in fuel performance calculations and use historical benchmark-quality measurement data to demonstrate it. More specifically, the developed methodology is applied to the OECD/NEA Multi-physics Pellet Cladding Mechanical Interaction Validation benchmark. A cold ramp test leading to PCMI is modeled. Two measured quantities of interest are considered: the cladding axial elongation during the irradiations and the cladding outer diameter after the cold ramp. The fuel performance code used to perform the simulation is FAST. The developed methodology involves various steps including a Morris screening to decrease the number of uncertain inputs, a nested loop approach for propagating the epistemic and aleatoric sources of uncertainties, and a global SA using Sobol indices. The obtained results indicate that the fuel and cladding thermal conductivities as well as the cladding outer diameter uncertainties are the three inputs having the largest impact on the measured quantities. More importantly, it was found that the epistemic uncertainties can have a significant impact on the measured quantities and can affect the outcome of the global sensitivity analysis.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2021

Data-Driven Reduced-Order Modeling of Convective Heat Transfer in Porous Media

This work presents a data-driven Reduced-Order Model (ROM) for parametric convective heat transfer problems in porous media. The intrusive Proper Orthogonal Decomposition aided Reduced-Basis (POD-RB) technique is employed to reduce the porous medium formulation of the incompressible Reynolds-Averaged Navier–Stokes (RANS) equations coupled with heat transfer. Instead of resolving the exact flow configuration with high fidelity, the porous medium formulation solves a homogenized flow in which the fluid-structure interactions are captured via volumetric flow resistances with nonlinear, semi-empirical friction correlations. A supremizer approach is implemented for the stabilization of the reduced fluid dynamics equations. The reduced nonlinear flow resistances are treated using the Discrete Empirical Interpolation Method (DEIM), while the turbulent eddy viscosity and diffusivity are approximated by adopting a Radial Basis Function (RBF) interpolation-based approach. The proposed method is tested using a 2D numerical model of the Molten Salt Fast Reactor (MSFR), which involves the simulation of both clean and porous medium regions in the same domain. For the steady-state example, five model parameters are considered to be uncertain: the magnitude of the pumping force, the external coolant temperature, the heat transfer coefficient, the thermal expansion coefficient, and the Prandtl number. For transient scenarios, on the other hand, the coastdown-time of the pump is the only uncertain parameter. The results indicate that the POD-RB-ROMs are suitable for the reduction of similar problems. The relative L2 errors are below 3.34% for every field of interest for all cases analyzed, while the speedup factors vary between 54 (transient) and 40,000 (steady-state).

German, Péter (ORCID:0000000307285283)↗

DOE OSTI2025

Aircraft in-situ measurements from SOCRATES constrain the anthropogenic perturbations of cloud droplet number

Aerosol-cloud interactions (ACI) in warm clouds alter reflected shortwave radiation by influencing cloud microphysical and macrophysical properties. The variable of state controlling ACI is the cloud droplet number concentration (N d ). Here, we examine the perturbations in N d due to anthropogenic aerosols (ΔN d, PD-PI ) using a perturbed parameter ensemble (PPE) hosted in the sixth Community Atmosphere Model (CAM6). Surrogate models are created for the CAM6 PPE outputs and are used to generate 1 million model variants of CAM6 by sampling 45 sources of parameter uncertainty. The range of uncertain physical parameters related to ACI are constrained with observations of aerosol and cloud properties from SOCRATES. The likely range of uncertain parameters and the associated range of ΔN d, PD-PI are more strongly constrained with observations of N d relative to observations of cloud condensation nuclei. We conduct sensitivity tests of how constraints on ΔN d, PD-PI are affected by systematic uncertainties in observations and our limitations in our surrogate models created for CAM6 PPE outputs. Based on this, we provide guidance on the impact of reducing systematic uncertainty in airborne microphysical observations and in surrogate models.

Song, Ci [Univ. of Wyoming, Laramie, WY (United St↗

DOE OSTI2025

Constraining aerosol–cloud adjustments by uniting surface observations with a perturbed parameter ensemble

Abstract. Aerosol–cloud interactions (ACIs) are the largest source of uncertainty in inferring the magnitude of future warming consistent with the observational record. The effective radiative forcing due to ACI (ERFaci) is dominated by liquid clouds and is composed of two terms: the change in cloud albedo due to redistributing liquid over a larger number of cloud droplets (Nd) and the change in cloud macrophysical properties due to changes in cloud microphysics. These terms are, respectively, referred to as the radiative forcing due to ACI (RFaci) and aerosol–cloud adjustments. While the magnitude of RFaci is uncertain, its sign is confidently negative and results in a cooling in the historical record. In contrast, the adjustment of cloud liquid water path (LWP) to enhanced Nd and associated radiative forcing is uncertain in sign. Increased LWP in response to increased Nd is consistent with precipitation suppression, while decreased LWP in response to increased Nd is consistent with enhanced evaporation from cloud top. Observational constraints of these processes are poor in part because of causal ambiguity in the relationship between Nd and LWP. To better understand this relationship, precipitation (P), Nd, and LWP surface observations from the Eastern North Atlantic (ENA) atmospheric observatory are combined with the output from a perturbed parameter ensemble (PPE) hosted in the Community Atmosphere Model version 6 (CAM6). This allows for causal interpretation of observed covariability. Observations of precipitation and cloud from ENA constrain the range of possible LWP aerosol–cloud adjustments relative to the prior from the PPE by 15 %, resulting in a global value that is confidently positive (a historical cooling) ranging from 2.1 to 6.9 g m−2. It is found that observed covariability between Nd and LWP is dominated by coalescence scavenging and that this observed covariability is not strongly related to aerosol–cloud adjustments.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Evaluating alternative ebullition models for predicting peatland methane emission and its pathways via data–model fusion

Abstract. Understanding the dynamics of peatland methane (CH4) emissions and quantifying sources of uncertainty in estimating peatland CH4 emissions are critical for mitigating climate change. The relative contributions of CH4 emission pathways through ebullition, plant-mediated transport, and diffusion, together with their different transport rates and vulnerability to oxidation, determine the quantity of CH4 to be oxidized before leaving the soil. Notwithstanding their importance, the relative contributions of the emission pathways are highly uncertain. In particular, the ebullition process is more uncertain and can lead to large uncertainties in modeled CH4 emissions. To improve model simulations of CH4 emission and its pathways, we evaluated two model structures: (1) the ebullition bubble growth volume threshold approach (EBG) and (2) the modified ebullition concentration threshold approach (ECT) using CH4 flux and concentration data collected in a peatland in northern Minnesota, USA. When model parameters were constrained using observed CH4 fluxes, the CH4 emissions simulated by the EBG approach (RMSE = 0.53) had a better agreement with observations than the ECT approach (RMSE = 0.61). Further, the EBG approach simulated a smaller contribution from ebullition but more frequent ebullition events than the ECT approach. The EBG approach yielded greatly improved simulations of pore water CH4 concentrations, especially in the deep soil layers, compared to the ECT approach. When constraining the EBG model with both CH4 flux and concentration data in model–data fusion, uncertainty of the modeled CH4 concentration profiles was reduced by 78 % to 86 % in comparison to constraints based on CH4 flux data alone. The improved model capability was attributed to the well-constrained parameters regulating the CH4 production and emission pathways. Our results suggest that the EBG modeling approach better characterizes CH4 emission and underlying mechanisms. Moreover, to achieve the best model results both CH4 flux and concentration data are required to constrain model parameterization.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

A multi-model assessment of the early last deglaciation (PMIP4 LDv1): a meltwater perspective

Abstract. The last deglaciation (∼20–11 ka BP) is a period of a major, long-term climate transition from a glacial to interglacial state that features multiple centennial- to decadal-scale abrupt climate variations whose root cause is still not fully understood. To better understand this time period, the Paleoclimate Modelling Intercomparison Project (PMIP) has provided a framework for an internationally coordinated endeavour in simulating the last deglaciation whilst encompassing a broad range of models. Here, we present a multi-model intercomparison of 17 transient simulations of the early part of the last deglaciation (∼20–15 ka BP) from nine different climate models spanning a range of model complexities and uncertain boundary conditions and forcings. The numerous simulations available provide the opportunity to better understand the chain of events and mechanisms of climate changes between 20 and 15 ka BP and our collective ability to simulate them. We conclude that the amount of freshwater forcing and whether it follows the ice sheet reconstruction or induces an inferred Atlantic meridional overturning circulation (AMOC) history, heavily impacts the deglacial climate evolution for each simulation rather than differences in the model physics. The course of the deglaciation is consistent between simulations except when the freshwater forcing is above 0.1 Sv – at least 70 % of the simulations agree that there is warming by 15 ka BP in most places excluding the location of meltwater input. For simulations with freshwater forcings that exceed 0.1 Sv from 18 ka BP, warming is delayed in the North Atlantic and surface air temperature correlations with AMOC strength are much higher. However, we find that the state of the AMOC coming out of the Last Glacial Maximum (LGM) also plays a key role in the AMOC sensitivity to model forcings. In addition, we show that the response of each model to the chosen meltwater scenario depends largely on the sensitivity of the model to the freshwater forcing and other aspects of the experimental design (e.g. CO2 forcing or ice sheet reconstruction). The results provide insight into the ability of our models to simulate the first part of the deglaciation and how choices between uncertain boundary conditions and forcings, with a focus on freshwater fluxes, can impact model outputs. We can use these findings as helpful insight in the design of future simulations of this time period.

Snoll, Brooke (ORCID:0000000230914394)↗

DOE OSTI2023

Rainbows and Climate Change: A tutorial on climate model diagnostics and parameterization

Earth System Models (ESMs) must represent processes below the grid scale of a model using representations (parameterizations) of physical and chemical processes. As a tutorial exercise to understand diagnostics and parameterization, this work presents a representation of rainbows for an ESM: the Community Earth System Model version 2 (CESM2). Using the 'state' of the model, basic physical laws, and some assumptions, we generate a representation of this unique optical phenomena as a diagnostic output. Rainbow occurrence and it's possible changes are related to cloud occurrence and rain formation which are critical uncertainties for climate change prediction. The work highlights issues which are typical of many diagnostics parameterizations such as assumptions, uncertain parameters and the difficulty of evaluation against uncertain observations. Results agree qualitatively with limited available global 'observations' of Rainbows. Rainbows are seen in expected locations in the sub-tropics over the ocean where broken clouds and frequent precipitation occurs. The diurnal peak is in the morning over ocean and in the evening over land. The representation of rainbows is found to be quantitatively sensitive to the assumed amount of cloudiness and the amount of stratiform rain. Rainbows are projected to have decreased, mostly in the Northern Hemisphere, due to aerosol pollution effects increasing cloud coverage since 1850. In the future, continued climate change is projected to decrease cloud cover, associated with a positive cloud feedback. As a result the rainbow diagnostic projects that rainbows will increase in the future, with the largest changes at mid-latitudes. The diagnostic may be useful for assessing cloud parameterizations, and is an exercise in how to build and test parameterizations of atmospheric phenomena.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

ISMIP6-based Antarctic projections to 2100: simulations with the BISICLES ice sheet model

The contribution of the Antarctic Ice Sheet is one of the most uncertain components of sea level rise to 2100. Ice sheet models are the primary tool for projecting future sea level contribution from continental ice sheets. The Ice Sheet Model Intercomparison for the Coupled Model Intercomparison Phase 6 (ISMIP6) provided projections of the ice sheet contribution to sea level over the 21st century, quantifying uncertainty due to ice sheet model, climate model, emission scenario, and uncertain parameters. We present simulations following the ISMIP6 framework with the BISICLES ice sheet model and new experiments extending the ISMIP6 protocol to more comprehensively sample uncertainties in future climate, ice shelf sensitivity to ocean melting, and their interactions. These results contributed to the land ice projections of Edwards et al. (2021), which formed the basis of sea level projections for the Sixth Assessment Report of the Intergovernmental Panel on Climate Change (AR6). Our experiments show the important interplay between surface mass balance processes and ocean-driven melt in determining Antarctic sea level contribution. Under higher-warming scenarios, high accumulation offsets more ocean-driven mass loss when sensitivity to ocean-driven melt is low. Conversely, we show that when sensitivity to ocean warming is high, ocean melting drives increased mass loss despite high accumulation. Overall, we simulate a sea level contribution range across our experiments from 2 to 178 mm. Finally, we show that collapse of ice shelves due to surface warming increases sea level contribution by 25 mm relative to the no-collapse experiments, for both moderate and high sensitivity of ice shelf melting to ocean forcing.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

A comparison of eight optimization methods applied to a wind farm layout optimization problem

Abstract. Selecting a wind farm layout optimization method is difficult. Comparisons between optimization methods in different papers can be uncertain due to the difficulty of exactly reproducing the objective function. Comparisons by just a few authors in one paper can be uncertain if the authors do not have experience using each algorithm. In this work we provide an algorithm comparison for a wind farm layout optimization case study between eight optimization methods applied, or directed, by researchers who developed those algorithms or who had other experience using them. We provided the objective function to each researcher to avoid ambiguity about relative performance due to a difference in objective function. While these comparisons are not perfect, we try to treat each algorithm more fairly by having researchers with experience using each algorithm apply each algorithm and by having a common objective function provided for analysis. The case study is from the International Energy Association (IEA) Wind Task 37, based on the Borssele III and IV wind farms with 81 turbines. Of particular interest in this case study is the presence of disconnected boundary regions and concave boundary features. The optimization methods studied represent a wide range of approaches, including gradient-free, gradient-based, and hybrid methods; discrete and continuous problem formulations; single-run and multi-start approaches; and mathematical and heuristic algorithms. We provide descriptions and references (where applicable) for each optimization method, as well as lists of pros and cons, to help readers determine an appropriate method for their use case. All the optimization methods perform similarly, with optimized wake loss values between 15.48 % and 15.70 % as compared to 17.28 % for the unoptimized provided layout. Each of the layouts found were different, but all layouts exhibited similar characteristics. Strong similarities across all the layouts include tightly packing wind turbines along the outer borders, loosely spacing turbines in the internal regions, and allocating similar numbers of turbines to each discrete boundary region. The best layout by annual energy production (AEP) was found using a new sequential allocation method, discrete exploration-based optimization (DEBO). Based on the results in this study, it appears that using an optimization algorithm can significantly improve wind farm performance, but there are many optimization methods that can perform well on the wind farm layout optimization problem, given that they are applied correctly.

17 WIND ENERGY↗

DOE OSTI2025

Data for Unraveling the Drivers of Water Shortage across Spatial Scales and Sectors in Colorado's West Slope River Basins

Results from Sunkara et al. submitted to Earth's Future. All code to reproduce the experiment and make the figures can be found here: https://github.com/ssaiveena/Sunkara-etal_2025_WestSlope For a detailed guide to data in this repository, see the README.txt file. Shortage output from .xdd files generated by StateMod was compressed into .parquet files in the Uncertain_scenarios folder. Reservoir output from .xre files can be found in the "Reservoir" directory within the Uncertain_scenarios folder. We adopted the CDSS naming convention: cm = Upper Colorado River Basin gm = Gunnison River Basin ym = Yampa River Basin wm = White River Basin sj = Southwest Basin

Climate Change↗

DOE OSTI2024

Recent Advances of PyROS: A Pyomo Solver for Nonconvex Two-Stage Robust Optimization in Process Systems Engineering

The document presents recent algorithmic and implementation advances of the two-stage robust optimization (RO) solver PyROS, and a benchmarking study which demonstrates the utility of PyROS for two-stage RO problems. The advances include extensions of the scope of PyROS to models with uncertain variable bounds, improvements to the initializations of the subproblems used by the underlying cutting set algorithm, and extensions of the uncertainty set interfaces. The benchmarking study is performed on a library of over 8,500 instances, with variations in the nonlinearities, degree-of-freedom partitioning, uncertainty sets, and polynomial decision rule approximations. An amine-based CO2 capture case study is presented to demonstrate the utility of PyROS for large-scale process models. Overall, the results highlight the effectiveness of PyROS for obtaining robust solutions to optimization problems with uncertain equality constraints.

Sherman, Jason↗

NASA NTRS1972

Minimax design of Kalman-like filters in the presence of large parameter uncertainties.

An attempt is made to guide the designer by means of a systematic procedure which will yield a filter that performs in an acceptable fashion when there is some uncertainty in the various system parameters. Very general restrictions are placed upon these uncertainties, and it is shown that they are applicable in most practical situations. The minimax approach results in a unique fixed filter design which places a least upper bound on a given sensitivity measure over the assumed range of uncertain parameters. The resulting filters are identical in form to the Kalman filter, provide nearly optimal performance over the entire range of uncertain statistics, and are independent of the actual noise statistics. Examples are given.

Hutchinson, C. E.↗

NASA NTRS1976

Computer program for single input-output, single-loop feedback systems

Additional work is reported on a completely automatic computer program for the design of single input/output, single loop feedback systems with parameter uncertainly, to satisfy time domain bounds on the system response to step commands and disturbances. The inputs to the program are basically the specified time-domain response bounds, the form of the constrained plant transfer function and the ranges of the uncertain parameters of the plant. The program output consists of the transfer functions of the two free compensation networks, in the form of the coefficients of the numerator and denominator polynomials, and the data on the prescribed bounds and the extremes actually obtained for the system response to commands and disturbances.

Source record↗

NASA NTRS1986

Integrated control/structure design and robustness

When a flexible structure is to be controlled actively, optimum performance is obtained by integrated, or simultaneous, design of the structure and the controller, as opposed to the common practice of designing the structure independently of control consideration and then designing a controller for a fixed structure. The primary design objective from the structural point of view usually is to minimize weight, while the control design objectives depend on the application. An important requirement for a practical control system is robustness with respect to uncertain plant parameters. This paper discusses simultaneous control/structure design when the overall design objective combines the weight of the structure and the robustness of the closed-loop control system. For numerical optimization, robustness is represented by the sensitivity of the closed-loop eigenvalues with respect to uncertain parameters. An example illustrates the optimal design of a flexible structure along with a robust compensator.

Adamian, A.↗

NASA NTRS1988

The mean density and two-point correlation function for the CfA redshift survey slices

The effect of large-scale inhomogeneities on the determination of the mean number density and the two-point spatial correlation function were investigated for two complete slices of the extension of the Center for Astrophysics (CfA) redshift survey (de Lapparent et al., 1986). It was found that the mean galaxy number density for the two strips is uncertain by 25 percent, more so than previously estimated. The large uncertainty in the mean density introduces substantial uncertainty in the determination of the two-point correlation function, particularly at large scale; thus, for the 12-deg slice of the CfA redshift survey, the amplitude of the correlation function at intermediate scales is uncertain by a factor of 2. The large uncertainties in the correlation functions might reflect the lack of a fair sample.

De Lapparent, Valerie↗

NASA NTRS1988

Analysis and use of VAS satellite data

A series of interrelated investigations has examined the analysis and use of VAS (VISSR Atmospheric Sounder) satellite data. A case study of VAS-derived mesoscale stability parameters suggested that they would have been a useful supplement to conventional data in the forecasting of thunderstorms on the day of interest. However, the meteorological significance of small or short lived stability features was uncertain. A second investigation examined the roles of first guess and VAS radiometric data in producing sounding retrievals. The radiance data often did not have a decisive influence on the final satellite soundings. Broad-scale patterns of the first guess, radiances, and retrievals frequently were similar, whereas small scale retrieval features, especially in the dew points, were often of uncertain origin.

Fuelberg, Henry E.↗

NASA NTRS1989

AGARD standard aeroelastic configurations for dynamic response

Since emphasis is on the transonic speed range, special importance is placed on configurations for which available data are sufficient to define accurately a transonic flutter boundary. Only configurations with clean, smooth surfaces are considered suitable. Segmented models or models with surface-slope discontinuities are inappropriate. Excluded also, in general, are configurations and data sets that involve behavior that is uncertain or not well understood, uncertain model properties, or know sensitivities to small variations in model properties. In order to assess the suitability of configurations already tested and the associated data for designation as standard, a survey of AGARD member countries was conducted to seek candidates for the prospective set. The results of that survey are given and summarized along with the initial selection of a standard configuration.

Yates, E. Carson, Jr.↗

NASA NTRS1989

Computational analysis of a stability robustness margin for structured real-parameter perturbations

An efficient computational method is presented for stability robustness analysis with structured real-parameter perturbations. A generic model of a class of uncertain dynamical systems is used as an example. The parameter uncertainty is characterized by a real scalar, epsilon. Multilinearity of the closed-loop characteristic polynomial is exploited to permit application of the mapping theorem to calculate the stability robustness margin. It is found that sensitive geometry of the stability boundary in the epsilon, omega-plane renders problematic the calculation of the minimum epsilon as a function of omega. This difficulty is avoided by calculating the minimum distance to the image of the uncertainty domain over omega as a function of epsilon. It is also shown that a certain class of uncertain dynamical systems has the required multilinearity property and are thus amenable to the proposed technique.

Wedell, Evan↗