Search NASA⌕ Search

SEARCH · Search NASA

Results for “global 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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Design and analysis of numerical experiments

Methods for the design and analysis of numerical experiments that are especially useful and efficient in multidimensional parameter spaces are presented. The analysis method, which is similar to kriging in the spatial analysis literature, fits a statistical model to the output of the numerical model. The method is applied to a fully nonlinear, global, equivalent-barotropic dynamical model. The statistical model also provides estimates for the uncertainty of predicted numerical model output, which can provide guidance on where in the parameter space to conduct further experiments, if necessary. The method can provide significant improvements in the efficiency with which numerical sensitivity experiments are conducted.

Bowman, Kenneth P.↗
NASA NTRS1997

Technical Report Series on Global Modeling and Data Assimilation: Interannual Variability and Potential Predictability in Reanalysis Products - Volume 13

The Data Assimilation Office (DAO) at Goddard Space Flight Center and the National Center for Environmental Prediction and National Center for Atmospheric Research (NCEP/NCAR) have produced multi-year global assimilations of historical data employing fixed analysis systems. These "reanalysis" products are ideally suited for studying short-term climatic variations. The availability of multiple reanalysis products also provides the opportunity to examine the uncertainty in the reanalysis data. The purpose of this document is to provide an updated estimate of seasonal and interannual variability based on the DAO and NCEP/NCAR reanalyses for the 15-year period 1980-1995. Intercomparisons of the seasonal means and their interannual variations are presented for a variety of prognostic and diagnostic fields. In addition, atmospheric potential predictability is re-examined employing selected DAO reanalysis variables.

Min, Wei↗
DOE OSTI2026

Uncertainty quantification for competing failure mechanisms in unidirectionally reinforced carbon–carbon composites

Microstructure-informed finite element models play a key role in the carbon–carbon composite design process. Variability in manufacturing process parameters and experimental limitations introduce model parameter uncertainty. This study quantifies the effect of model parameter uncertainty on transverse tensile fracture behavior and proposes a methodology to predict the failure mode based on competing microscale damage mechanisms. Finite element simulations incorporate fiber–matrix interface debonding with cohesive zones and matrix damage with a smeared crack band approach in a unidirectional carbon–carbon composite. Results from a variance-based global sensitivity analysis identifies interfacial and matrix damage parameters as the primary source of variability in fracture behavior. Sobol’ indices indicate that matrix and cohesive zone strengths contribute 94% of the variance in the effective ultimate stress. A local analysis elucidates the relationship between these constituent strength parameters and failure mode by estimating the probability of cohesive, matrix, and mixed-mode dominated failure. Based on the results for 4000 simulations, 93% exhibit mixed-mode or interfacial dominated failure, which underscores the crucial role of fiber–matrix interface debonding in the transverse tensile failure of carbon–carbon composites. These uncertainty quantification results facilitate more efficient model calibration and provide a framework for microstructure-informed failure predictions in the face of manufacturing-induced uncertainty.

36 MATERIALS SCIENCE↗
DOE OSTI2022

LES of HCCI combustion of iso-octane/air in a flat-piston rapid compression machine

Homogeneous Charge Compression Ignition (HCCI) engines promise better efficiency and cleaner emissions than conventional piston engines, but can be challenging to control. Rapid compression machines (RCM) provide a simplified configuration for investigating HCCI combustion behavior, which is necessary for effective control of engine ignition timing and peak pressures. Here, in this study, we assess the utility of large eddy simulations (LES) for predicting HCCI combustion in a 3-D configuration. To this end, LES with finite-rate chemistry employing a 99-species iso-octane/air mechanism of two RCM operating conditions are performed. The RCM configuration under consideration was designed by Strozzi et al. with a flat piston to introduce large amounts of thermal stratification representative of realistic HCCI engine conditions, through the generation of corner vortices. It is shown that the simulation provides reasonable agreement with temperature fluctuations (7% difference), as well as ignition delay in the short ignition case (1 ms difference), while the long ignition case (35 ms difference) highlights more substantial deficiencies that are still within the expected uncertainty from the employed chemical mechanism. Flame propagation modes predicted by LES agree with experimental observations: spontaneous ignition is seen in the short ignition case, while deflagration is more prominent in the long ignition case. Analysis of global and local quantities classify the short ignition case in a mixed ignition regime, and the long ignition case in the mild ignition regime. These results demonstrate the utility of FRC-LES for investigations of multimode combustion regimes of HCCI combustion in a 3-D configuration.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

Uncertainties in greenhouse gas emission factors: A comprehensive analysis of switchgrass‐based biofuel production

Abstract This study investigates uncertainties in greenhouse gas (GHG) emission factors related to switchgrass‐based biofuel production in Michigan. Using three life cycle assessment (LCA) databases—US lifecycle inventory (USLCI) database, GREET, and Ecoinvent—each with multiple versions, we recalculated the global warming intensity (GWI) and GHG mitigation potential in a static calculation. Employing Monte Carlo simulations along with local and global sensitivity analyses, we assess uncertainties and pinpoint key parameters influencing GWI. The convergence of results across our previous study, static calculations, and Monte Carlo simulations enhances the credibility of estimated GWI values. Static calculations, validated by Monte Carlo simulations, offer reasonable central tendencies, providing a robust foundation for policy considerations. However, the wider range observed in Monte Carlo simulations underscores the importance of potential variations and uncertainties in real‐world applications. Sensitivity analyses identify biofuel yield, GHG emissions of electricity, and soil organic carbon (SOC) change as pivotal parameters influencing GWI. Decreasing uncertainties in GWI may be achieved by making greater efforts to acquire more precise data on these parameters. Our study emphasizes the significance of considering diverse GHG factors and databases in GWI assessments and stresses the need for accurate electricity fuel mixes, crucial information for refining GWI assessments and informing strategies for sustainable biofuel production.

Kim, Seungdo↗
NASA NTRS1989

Trends in total column ozone measurements

It is important to ensure the best available data are used in any determination of possible trends in total ozone in order to have the most accurate estimates of any trends and the associated uncertainties. Accordingly, the existing total ozone records were examined in considerable detail. Once the best data set has been produced, the statistical analysis must examine the data for any effects that might indicate changes in the behavior of global total ozone. The changes at any individual measuring station could be local in nature, and herein, particular attention was paid to the seasonal and latitudinal variations of total ozone, because two dimensional photochemical models indicate that any changes in total ozone would be most pronounced at high latitudes during the winter months. The conclusions derived from this detailed examination of available total ozone can be split into two categories, one concerning the quality and the other the statistical analysis of the total ozone record.

Rowland, F. S.↗
DOE OSTI2021

An efficient method to identify uncertainties of WRF-Solar variables in forecasting solar irradiance using a tangent linear sensitivity analysis

Uncertainty in predicting solar energy resources introduces major challenges in power system management and necessitates the development of reliable probabilistic solar forecasts. As the first part of the development of probabilistic forecasts based on the Weather Research and Forecasting model with solar extensions (WRF-Solar), this study presents a tangent linear approach to identify input variables responsible for the largest uncertainties in predicting surface solar irradiance and clouds. A tangent linear analysis is capable of efficiently investigating sensitivities of output variables with respect to various input variables of WRF-Solar because this approach avoids the computational burden of perturbing the initial conditions of individual input variables. We develop tangent linear models (TLMs) for six WRF-Solar physics packages that control the formation and dissipation of clouds and solar radiation, and we evaluate the validity of TLMs using a linearity test. The tangent linear sensitivity analysis is conducted under various scenarios based on satellite observations and model simulations to consider realistic input conditions. A simple method is used to quantify the impact of the uncertainty of input variables on the output variables from the TLMs. The results demonstrate that uncertainties in the output variables that are the focus of this study—including global horizontal irradiance, direct normal irradiance, cloud mixing ratio, cloud tendency, cloud fraction, and sensible and latent heat fluxes—are highly sensitive to uncertainties in 14 input variables. This study indicates that the tangent linear method can identify key variables of physics modules in WRF-Solar that can be stochastically perturbed to generate ensemble-based probabilistic forecasts.

14 SOLAR ENERGY↗
NASA NTRS2016

Influence of Measurement Uncertainties on Fractional Solubility of Iron in Mineral Aerosols over the Oceans

The atmospheric supply of mineral dust iron (Fe) plays a crucial role in the Earths biogeochemical cycle and is of specific importance as a micronutrient in the marine environment. Observations show several orders of magnitude variability in the fractional solubility of Fe in mineral dust aerosols, making it hard to assess the role of mineral dust for the global ocean biogeochemical Fe cycle. In this study we compare the operational solubility of mineral dust aerosol Fe associated with the flow-through leaching protocol to the results of the global 3-D chemical transport model GEOS-Chem. According to the protocol, aerosol Fe is defined as soluble by first deionized water leaching of mineral dust through a 0.45 m pore size membrane followed by acidification and storage of the leachate over a long period of time prior to analysis. To estimate the uncertainty in soluble Fe results introduced by the flow-through leaching protocol, we prescribe an average 50 (range of 30 to 70) fractional solubility to sub-0.45 m sized mineral dust particles that may inadvertently pass the filter and end up in the acidified (at pH1.7) leachate for a couple of month period. In the model, the fractional solubility of Fe is either explicitly calculated using a complex mineral aerosol Fe dissolution module, or prescribed to be 1 and 4. Calculations show that the fractional solubility of Fe derived through the flow-through leaching is higher compared to the model results. The largest differences (40) are predicted to occur farther away from the dust source regions, over the areas where sub-0.45 m sized mineral dust particles contribute a larger fraction of the total mineral dust mass. This study suggests that different methods used in soluble Fe measurements and inconsistences in the operational definition of filterable Fe in marine environment and soluble Fe in atmospheric aerosols are likely to contribute to the wide range of fractional solubility of aerosol Fe reported in the literature.

Soluble Iron↗
DOE OSTI2021

A probabilistic economic and environmental impact assessment of a cyanobacteria-based biorefinery

Microbial based biofuels represent a potential promising solution as an environmentally favorable transportation fuel. Cyanobacteria have many of the same advantages as microalgae: ability for rapid growth in otherwise non-arable regions, suitability for genetic engineering, and simple nutritional needs. Additionally, cyanobacteria can be engineered to secrete valuable co-products that can be harvested independent from the produced biomass. However, little work has been done to identify the processes and the economic and environmental impacts associated with a large-scale cyanobacteria-to-fuels facility. The present study is a concurrent techno-economic and life cycle assessment of a facility that generates fuels and methyl laurate, an oleochemical, from the cyanobacterial species Synechocystis sp. PCC 6803. Here, the biorefinery model includes all aspects of cultivation, separation of the secreted methyl laurate, biomass harvesting and fuel processing via hydrothermal liquefaction (HTL) of the dewatered biomass. The assessments leverage Monte Carlo analysis (MCA) to address uncertainty and variability inherent in the most significant input parameters, replacing them with probabilistic functions. For the facility configuration producing both fuels and the oleochemical co-product, the MCA average minimum fuel selling price (MFSP) is $\$2.47$ per decimeter (dm 3 ) or $\$9.34$ per gallon of gasoline equivalent (gge) with the corresponding average global warming potential determined to be 118 g CO 2 -eq-MJ -1 . The case producing only fuels results in an MCA average MFSP of $\$2.01$-(dm 3 ) -1 ($\$7.60$-gge -1 ) and an average environmental impact of 100 g CO 2 -eq-MJ -1 . These results are compared to static optimistic and conservative scenario analysis estimates, illustrating the over- and under-estimation of outcomes associated with non-stochastic methods. Suggested facility improvements include increases in pond productivity of both the biomass and methyl laurate oil production, as well as improvements to carbon utilization and bio-crude yield from HTL processing.

09 BIOMASS FUELS↗
DOE OSTI2021

The updated ITPA global H-mode confinement database: description and analysis

The multi-machine ITPA Global H-mode Confinement Database has been upgraded with new data from JET with the ITER-like wall and ASDEX Upgrade with the full tungsten wall. This paper describes the new database and presents results of regression analysis to estimate the global energy confinement scaling in H-mode plasmas using a standard power law. Various subsets of the database are considered, focusing on type of wall and divertor materials, confinement regime (all H-modes, ELMy H or ELM-free) and ITER-like constraints. Apart from ordinary least squares, two other, robust regression techniques are applied, which take into account uncertainty on all variables. Regression on data from individual devices shows that, generally, the confinement dependence on density and the power degradation are weakest in the fully metallic devices. Using the multi-machine scalings, predictions are made of the confinement time in a standard ELMy H-mode scenario in ITER. The uncertainty on the scaling parameters is discussed with a view to practically useful error bars on the parameters and predictions. One of the derived scalings for ELMy H-modes on an ITER-like subset is studied in particular and compared to the IPB98(y,2) confinement scaling in engineering and dimensionless form. Transformation of this new scaling from engineering variables to dimensionless quantities is shown to result in large error bars on the dimensionless scaling. Regression analysis in the space of dimensionless variables is therefore proposed as an alternative, yielding acceptable estimates for the dimensionless scaling. The new scaling, which is dimensionally correct within the uncertainties, suggests that some dependencies of confinement in the multi- machine database can be reconciled with parameter scans in individual devices. This includes vanishingly small dependence of confinement on line-averaged density and normalized plasma pressure (β), as well as a noticeable, positive dependence on effective atomic mass and plasma triangularity. Extrapolation of this scaling to ITER yields a somewhat lower confinement time compared to the IPB98(y, 2) prediction, possibly related to the considerably weaker dependence on major radius in the new scaling (slightly above linear). Further studies are needed to compare more flexible regression models with the power law used here. In addition, data from more devices concerning possible ‘hidden variables’ could help to determine their influence on confinement, while adding data in sparsely populated areas of the parameter space may contribute to further disentangling some of the global confinement dependencies in tokamak plasmas.

Database↗
DOE OSTI2024

A global urban heat island intensity dataset: Generation, comparison, and analysis

The urban heat island (UHI) effect, a phenomenon of local warming over urban areas, is the most well-known impact of urbanization on climate. Globally consistent estimates of the UHI intensity (UHII) are crucial for examining this phenomenon across time and space. However, publicly available UHII datasets are limited and have several constraints: (1) they are for clear-sky surface UHII, not all-sky surface UHII and canopy (air temperature) UHII; (2) the estimation methods often neglect anthropogenic disturbance, introducing uncertainties in the estimated UHII. To address these issues, this study proposes a new dynamic equal-area (DEA) method that can minimize the influence of various confounding factors on UHII estimates through a dynamic cyclic process. Utilizing the DEA method and leveraging various gridded temperature data, we develop a global-scale (>10,000 cities), long-term (over 20 years by month), and multi-faceted (clear-sky surface, all-sky surface, and canopy) UHII dataset. Further, based on these estimates, we provide a comprehensive analysis of the UHII and its trends in global cities. The UHII is found to be greater than zero in >80% of cities, with global annual average magnitudes around 1.0 °C (day) and 0.8 °C (night) for surface UHII, and close to 0.5 °C for canopy UHII. Furthermore, an interannual upward trend in UHII is observed in >60% of cities, with global annual average trends exceeding 0.1 °C/decade (day) and over 0.06 °C/decade (night) for surface UHII, and slightly surpassing 0.03 °C/decade for canopy UHII. Notably, there exists a positive correlation between the magnitude and trend of UHII, suggesting that cities with stronger UHII tend to experience faster growth in UHII. Additionally, discrepancies in UHII are found between different temperature data, stemming not only from distinctions in data types (surface or air temperature) but also from differences in data acquisition times (Terra or Aqua), weather conditions (clear-sky or all-sky), and processing methodologies (with or without gap filling). Overall, our proposed method, dataset, and analysis results have the potential to provide valuable insights for future urban climate studies. The UHII dataset is publicly available at https://doi.org/10.6084/m9.figshare.24821538.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Quantitative assessment of eddy viscosity rans models for turbulent mixed convection in a differentially heated plane channel

Turbulent mixed convection between two vertical, infinite parallel plates at different temperatures is studied using various two-equation turbulence models. The numerical simulations are performed at a turbulent Reynolds number of Re τ = 150 and a Grashof number of Gr = 9.6 × 10 5 . Comparisons are made against the highly trusted DNS results. Consistent with the DNS approach, the current simulations are performed using constant properties and the Boussinesq approximation to predict the influence of buoyancy. Previous studies have provided assessments of two-equation turbulence models for various scenarios, but often rely on a qualitative “eye” test in order to determine the most appropriate model to predict a given flow. This study aims to provide a new form of quantitative assessment that accounts for both the physics captured by the turbulence model as well as the magnitude of the system response quantities (SRQ) using a modified symmetric mean absolute percent error (SMAPE) method. This method is designed to be approachable to researchers at any level and can be applied to system response quantities from multiple research fields. Uncertainty quantification is also performed to determine the discretization error for each turbulence model. Recommendations are made as to which turbulence models best capture the physics – hydrodynamically and thermally – using both local and global validation metrics. Lastly, a sensitivity analysis is performed on the damping functions used in the most accurate models. This underpins the potential of model developments and adjustments most worth pursuing for buoyant flows. Finally, this framework provides a more physics-based comparative analysis of the selected turbulence models.

42 ENGINEERING↗
DOE OSTI2024

Analysis of parton distributions in a pion with Bézier parametrizations

We explore the role of parametrizations for nonperturbative QCD functions in global analyses, with a specific application to extending a phenomenological analysis of the parton distribution functions (PDFs) in the charged pion realized in the xFitter fitting framework. The parametrization dependence of PDFs in our pion fits substantially enlarges the uncertainties from the experimental sources estimated in the previous analyses. We systematically explore the parametrization dependence by employing a novel technique to automate generation of polynomial parametrizations for PDFs that makes use of Bézier curves. This technique is implemented in a ++ module that is included in the xFitter program. Our analysis reveals that the sea and gluon distributions in the pion are not well disentangled, even when considering measurements in leading-neutron deep inelastic scattering. For example, the pion PDF solutions with a vanishing gluon and large quark sea are still experimentally allowed, which elevates the importance of ongoing lattice and nonperturbative QCD calculations, together with the planned pion scattering experiments, for conclusive studies of the pion structure. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

New physics in rare B decays after Moriond 2021

Abstract The anomalies in rare B decays endure. We present results of an updated global analysis that takes into account the latest experimental input – in particular the recent results on $$R_K$$ R K and BR $$(B_s \rightarrow \mu ^+\mu ^-)$$ ( B s → μ + μ - ) – and that qualitatively improves the treatment of theory uncertainties. Fit results are presented for the Wilson coefficients of four-fermion contact interactions. We find that muon specific Wilson coefficients $$C_9 \simeq -0.73$$ C 9 ≃ - 0.73 or $$C_9 = -C_{10} \simeq -0.39$$ C 9 = - C 10 ≃ - 0.39 continue to give an excellent description of the data. If only theoretically clean observables are considered, muon specific $$C_{10} \simeq 0.60$$ C 10 ≃ 0.60 or $$C_9=-C_{10} \simeq -0.35$$ C 9 = - C 10 ≃ - 0.35 improve over the Standard Model by $$\sqrt{\Delta \chi ^2} \simeq 4.7\sigma $$ Δ χ 2 ≃ 4.7 σ and $$\sqrt{\Delta \chi ^2} \simeq 4.6\sigma $$ Δ χ 2 ≃ 4.6 σ , respectively. In various new physics scenarios we provide predictions for lepton flavor universality observables and CP asymmetries that can be tested with more data. We update our previous combination of ATLAS, CMS, and LHCb data on BR $$(B_s \rightarrow \mu ^+\mu ^-)$$ ( B s → μ + μ - ) and BR $$(B^0\rightarrow \mu ^+\mu ^-)$$ ( B 0 → μ + μ - ) taking into account the full two-dimensional non-Gaussian experimental likelihoods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

An uncertainty-aware strategy for plasma mechanism reduction with directed weighted graphs

In this work, we present a framework for the analysis and reduction of plasma mechanisms by means of weighted directed graphs, in which reactions and species are both treated as nodes. The methodology consists of two distinct analyses. The first, which is qualitative, relies on graph spatializations via force-directed algorithms to discover the predominant global patterns in the chemical model. The second ranks the reactions based on their shortest paths' lengths from/to the species of interest and their relative contributions to the power balance. Further, this quantitative investigation enables a strategy for mechanism reduction that is fully automatized, as it does not require any expert knowledge, highly effective, as it generates reduced mechanisms that are highly accurate while relying on a small number of processes, and easily interpretable, as the algorithm justifies the importance of the retained reactions by outputting their related chemical pathways. Additionally, the work proposes a methodology extension that employs ensembles of graphs to improve the robustness of the reduced mechanism to reaction parameter uncertainties. The approach, here tested for steady-state predictions of a plasma system characterizing negative hydrogen ion sources, is general and can be used in a wide variety of applications outside the particular nuclear fusion context demonstrated in this work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRSDate not supplied

NASA’s Prototype Spectral Water Inversion Processor and Emulator (SWIPE): Towards Global Coastal and Inland Water Quality and Algal Biodiversity Monitoring

Degradation of Earth’s inland water resources due to anthropogenic perturbations and climate anomalies at both local and global scales continues to place human health at substantial risk. There is now a growing necessity to develop pragmatic approaches that allow timely and effective extrapolation of local processes, to spatially resolved global products, and to promote operational and sustainable resource policy management. This presentation will provide updates on NASA’s prototype open-source aquatic modeling platform, Spectral Water Inversion Processor and Emulator (SWIPE), which is a comprehensive, multi-faceted modeling platform for both forward and inverse modeling of diverse aquatic ecosystems from the benthos to top-of-atmosphere (TOA). SWIPE provides a cohesive application which leverages recent advancements in particle modeling, Big Data analytics, and machine learning to develop a high-fidelity synthetic training ground for sensitivity studies and algorithm development for multispectral or upcoming hyperspectral missions. Some of the prominent features of SWIPE to be discussed include: 1. Advanced hyperspectral modeling of globally diverse algal and non-algal particles using a novel two-layer coated sphere scattering model and radiative transfer modeling, 2. Massive, highly detailed synthetic spectral libraries of Analysis-Ready-Data (ARD) which include spectral libraries of particle microphysics, water biogeophysical and optical properties, as well as surface and TOA reflectances at 1 nm resolution, 3. An ensemble of pre-built analytic, machine learning, and deep learning inversion algorithms for various water quality and biodiversity related retrieval parameters and uncertainty quantification, 4. Sensor-agnostic water quality inversion at wide ranging spatial and spectral resolutions including a codebase for seamless application in the Google Earth Engine and NASA Earth Exchange (NEX) for planetary scale analysis. SWIPE will be a fully open-source platform based in python with comprehensive documentation, tutorials, and options for distributed computing on high performance computing clusters or on single, local machines. Further, we will discuss how we envision SWIPE contributing towards a global analysis of coastal and inland water quality dynamics.

top-of-atmosphere (TOA)↗
DOE OSTI2021

Assessing the predictive impact of factor fixing with an adaptive uncertainty-based approach .

Despite widespread use of factor fixing in environmental modeling, its effect on model predictions has received little attention and is instead commonly presumed to be negligible. We propose a proof-of-concept adaptive method for systematically investigating the impact of factor fixing. The method uses Global Sensitivity Analysis methods to identify groups of sensitive parameters, then quantifies which groups can be safely fixed at nominal values without exceeding a maximum acceptable error, demonstrated using the 21-dimensional Sobol’ G-function. Furthermore, three error measures are considered for quantities of interest, namely Relative Mean Absolute Error, Pearson Product-Moment Correlation and Relative Variance. Results demonstrate that factor fixing may cause large errors in the model results unexpectedly, when preliminary analysis suggests otherwise, and that the default value selected affects the number of factors to fix. To improve the applicability and methodological development of factor fixing, a new research agenda encompassing five opportunities is discussed for further attention.

54 ENVIRONMENTAL SCIENCES↗