Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Variability 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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Fire–precipitation interactions amplify the quasi-biennial variability in fires over southern Mexico and Central America

Fires have great ecological, social, and economic impact. However, fire prediction and management remain challenges due to a limited understanding of their roles in the Earth system. Fires over southern Mexico and Central America (SMCA) are a good example of this, greatly impacting local air quality and regional climate. Here we report that the spring peak (April–May) of fire activities in this region has a distinct quasi-biennial signal based on multiple satellite datasets measuring different fire characteristics. The variability is initially driven by quasi-biennial variations in precipitation. Composite analysis indicates that strong fire years correspond to suppressed ascending motion and weakened precipitation over the SMCA. The anomalous precipitation over the SMCA is further found to be mostly related to the East Pacific–North Pacific (EP-NP) pattern 2 months prior to the fire season. The positive phase of the EP-NP leads to enhanced precipitation over the eastern US but suppressed precipitation over the SMCA, similar to the spatial pattern of precipitation differences between strong and weak fire years. Meanwhile, the quasi-biennial signals in precipitation and fires appear to be amplified by their interactions through a positive feedback loop at short timescales. Model simulations show that in strong fire years, more aerosol particles are released and transported downstream over the Gulf of Mexico and the eastern US, where suspended light-absorbing aerosols warm the atmosphere and cause the ascending motion of the air aloft. Subsequently, a compensating downward motion is formed over the region of the fire source and ultimately suppresses precipitation and intensifies fires. Statistical analysis shows the different durations of the two-way interaction, where the fire suppression effect of precipitation lasts for more than 20 d, while fire leads to a decrease in precipitation at shorter timescales (3–5 d). This study demonstrates the importance of fire–climate interactions in shaping the fire activities on an interannual scale and highlights how precipitation–fire interactions at short timescales contribute to the interannual variability in both fire and precipitation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Projecting climate change in South America using variable‐resolution Community Earth System Model: An application to Chile

Abstract We introduce variable‐resolution enabled Community Earth System Model (VR‐CESM) results simulating historical and future climate conditions at 28 km over South America and 14 km over the Andes. Three 30‐year simulations are performed: a historic (1985–2014), a near future (2030–2059), and an end‐century (2070–2099) simulation under the RCP8.5 scenario. Historic results compare favourably to several temperature and precipitation reanalysis products, though local biases are present, particularly during austral summer. Future simulations highlight broad warming patterns (+3–6°C by end‐century) and heterogeneous precipitation responses across South America that qualitatively agree with prior modelling efforts. Our results reveal that the interaction between temperature and precipitation changes produce shifts in several Köppen–Geiger climates. Notable changes include the near‐elimination of the Andean Tundra or Alpine climates, a 15% decrease in Tropical Rainforests and a Tropical Savannah expansion of 20%. To provide a regionally focused analysis of projected climate change and to illustrate the benefits of variable resolution modelling, we analyse changes in the magnitude and trend in seasonal and daily temperature and precipitation in Chile. We also examined several metrics [e.g., snow water equivalent (SWE), temperatures on wet days, and days below 0°C] to evaluate potential impacts of climate change on the Chilean cryosphere between the end‐of‐century and historic periods, finding wide‐ranging indications of cryospheric decline. These changes are interpreted through reductions in the timing (1–2.5 months earlier peak SWE) and magnitude (200–1,000 mm SWE decreases) of water stored as snow in the Andes, a 10–30% decrease in number of cool season wet days with temperatures below 1°C, and 50–200 fewer days (annually) with minimum temperatures below 0°C. Our aim in producing a high‐resolution dataset of climate projections from VR‐CESM is to support analyses of climate change throughout South America but especially in vulnerable montane regions and to provide additional results for comparison with previous, ongoing, and upcoming modelling efforts.

variable-resolution global climate models↗
DOE OSTI2021

FLUXNET-CH 4 : a global, multi-ecosystem dataset and analysis of methane seasonality from freshwater wetlands

Methane (CH 4 ) emissions from natural landscapes constitute roughly half of global CH4 contributions to the atmosphere, yet large uncertainties remain in the absolute magnitude and the seasonality of emission quantities and drivers. Eddy covariance (EC) measurements of CH4 flux are ideal for constraining ecosystem-scale CH 4 emissions due to quasi-continuous and high-temporal-resolution CH 4 flux measurements, coincident carbon dioxide, water, and energy flux measurements, lack of ecosystem disturbance, and increased availability of datasets over the last decade. Here, we (1) describe the newly published dataset, FLUXNET-CH 4 Version 1.0, the first open-source global dataset of CH 4 EC measurements (available at https://fluxnet.org/data/fluxnet-ch4-community-product/, last access: 7 April 2021). FLUXNET-CH 4 includes half-hourly and daily gap-filled and non-gap-filled aggregated CH4 fluxes and meteorological data from 79 sites globally: 42 freshwater wetlands, 6 brackish and saline wetlands, 7 formerly drained ecosystems, 7 rice paddy sites, 2 lakes, and 15 uplands. Then, we (2) evaluate FLUXNET-CH 4 representativeness for freshwater wetland coverage globally because the majority of sites in FLUXNET-CH4 Version 1.0 are freshwater wetlands which are a substantial source of total atmospheric CH 4 emissions; and (3) we provide the first global estimates of the seasonal variability and seasonality predictors of freshwater wetland CH4 fluxes. Our representativeness analysis suggests that the freshwater wetland sites in the dataset cover global wetland bioclimatic attributes (encompassing energy, moisture, and vegetation-related parameters) in arctic, boreal, and temperate regions but only sparsely cover humid tropical regions. Seasonality metrics of wetland CH 4 emissions vary considerably across latitudinal bands. In freshwater wetlands (except those between 20°S to 20°N) the spring onset of elevated CH 4 emissions starts 3 d earlier, and the CH 4 emission season lasts 4 d longer, for each degree Celsius increase in mean annual air temperature. On average, the spring onset of increasing CH 4 emissions lags behind soil warming by 1 month, with very few sites experiencing increased CH 4 emissions prior to the onset of soil warming. In contrast, roughly half of these sites experience the spring onset of rising CH 4 emissions prior to the spring increase in gross primary productivity (GPP). The timing of peak summer CH 4 emissions does not correlate with the timing for either peak summer temperature or peak GPP. Our results provide seasonality parameters for CH 4 modeling and highlight seasonality metrics that cannot be predicted by temperature or GPP (i.e., seasonality of CH 4 peak). FLUXNET-CH 4 is a powerful new resource for diagnosing and understanding the role of terrestrial ecosystems and climate drivers in the global CH 4 cycle, and future additions of sites in tropical ecosystems and site years of data collection will provide added value to this database. All seasonality parameters are available at https://doi.org/10.5281/zenodo.4672601 (Delwiche et al., 2021). Additionally, raw FLUXNET-CH 4 data used to extract seasonality parameters can be downloaded from https://fluxnet.org/data/fluxnet-ch4-community-product/ (last access: 7 April 2021), and a complete list of the 79 individual site data DOIs is provided in Table 2 of this paper.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Characterization of Electronic Stress-Induced Changes in Multilayer MoS 2

Transition metal dichalcogenides like molybdenum disulfide (MoS 2 ) are compelling for next-generation electronic devices. In this work, we investigate the impact of electronic stress on MoS 2 to illustrate that observational and phenomenological information on multiple devices can be useful to describe changes in the device, and caution against the rationalization of paltry results as representative or correlative to device behavior. Here, we stress MoS 2 by applying a sustained 20 V DC bias to study the material’s response. Post-stress electronic characterization revealed nonuniform shifts in current–voltage (I–V) behavior alongside microscale changes. Complementary mechanical, spectroscopic, and scanning microwave impedance measurements showed that stress-induced features locally modulate stiffness, surface potential, Raman intensity, and charge carrier density. We correlated I–V behavior with morphological features (wrinkles, tears, folds, height) and device-level geometry (MoS 2 overlap with electrodes, channel area, contact length) on 50 test structures across five chips to move beyond anecdotal conclusions. We found no universal correlations before DC stress. However, device-level geometry was correlated with I–V behavior after DC stress, suggesting that electrode contacts play a more dominant role than morphology in determining performance. Delamination and thinning induced by DC stress led to localized reductions in charge carrier density within the affected regions. Further, delamination and thinning appear to map to I–V device performance in a few samples, but the correlation is lost when a larger sample size is considered. This suggests significant sample-to-sample variability in surface electronic states of the test structures. We also discuss how environmental factors introduced during fabrication may contribute to the observed heterogeneous device response. Progress will require high-resolution, multimodal analysis across many samples constructed under controlled, clean conditions. By building data sets that capture variability, we can better identify the true drivers of performance.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A novel random forest approach to revealing interactions and controls on chlorophyll concentration and bacterial communities during coastal phytoplankton blooms

Abstract Increasing occurrence of harmful algal blooms across the land–water interface poses significant risks to coastal ecosystem structure and human health. Defining significant drivers and their interactive impacts on blooms allows for more effective analysis and identification of specific conditions supporting phytoplankton growth. A novel iterative Random Forests (iRF) machine-learning model was developed and applied to two example cases along the California coast to identify key stable interactions: (1) phytoplankton abundance in response to various drivers due to coastal conditions and land-sea nutrient fluxes, (2) microbial community structure during algal blooms. In Example 1, watershed derived nutrients were identified as the least significant interacting variable associated with Monterey Bay phytoplankton abundance. In Example 2, through iRF analysis of field-based 16S OTU bacterial community and algae datasets, we independently found stable interactions of prokaryote abundance patterns associated with phytoplankton abundance that have been previously identified in laboratory-based studies. Our study represents the first iRF application to marine algal blooms that helps to identify ocean, microbial, and terrestrial conditions that are considered dominant causal factors on bloom dynamics.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

What's in a Name? Developing a Standardized Taxonomy for HVAC System Faults

Faults occurring in heating, ventilation and air-conditioning (HVAC) systems have significantly negative impacts on building energy consumption, occupant comfort, and indoor air quality. In the past thirty years, extensive research has been conducted on fault detection and diagnostics (FDD) methods, and there are now dozens of commercially available FDD software tools. Growing adoption of FDD tools has the potential to generate a massive and useful data set on fault characteristics. However, the lack of a unifying taxonomy is a significant barrier to efficient analysis and evaluation of FDD outputs. Therefore, there is a strong need to develop a robust taxonomy which can better represent and interpret FDD output data. This paper documents the development of a unifying taxonomy for HVAC system faults in commercial buildings, with initial focus on air handling units, variable air volume terminal units, and roof top unit systems. The developed fault taxonomy employs both a physical hierarchy of HVAC equipment and a cause-effect relationship model as tools to better understand and support root cause analysis for HVAC faults. A variable air volume terminal unit is used as an example to demonstrate the application of the developed fault taxonomy. The taxonomy has short-term application in a major U.S. study on fault prevalence, and promises longer term benefits to FDD software developers and building operators by creating a foundation for improved approaches to identifying and resolving HVAC faults.

Chen, Yimin↗
DOE OSTI2021

3-D Geologic Controls of Hydrothermal Fluid Flow at Brady geothermal field, Nevada, USA

In many hydrothermal systems, fracture permeability along faults provides pathways for groundwater to transport heat from depth. Faulting generates a range of deformation styles that cross-cut heterogeneous geology, resulting in complex patterns of permeability, porosity, and hydraulic conductivity. Vertical connectivity (a throughgoing network of permeable areas that allows advection of heat from depth to the shallow subsurface) is rare and is confined to relatively small volumes that have highly variable spatial distribution. This local compartmentalization of connectivity represents a significant challenge to understanding hydrothermal circulation and for exploring, developing, and managing hydrothermal resources. Here, we present an evaluation of the geologic characteristics that control this compartmentalization in hydrothermal systems through 3-D analysis of the Brady geothermal field in western Nevada. A published 3-D geologic map of the Brady area is used as a basis to develop structural and geological variables that are hypothesized to control or effect permeability or connectivity. The 3-D distribution of these variables is compared to the distribution of productive and non-productive fluid flow intervals along production wells and non-productive wells via principal component analysis (PCA). This comparison elucidates which geologic and structural variables are most closely associated with productive fluid flow intervals. Results indicate that production intervals at Brady are located: (1) within or near to known and stress-loaded macro-scale faults, and (2) in areas of high fault and fracture density.

58 GEOSCIENCES↗
DOE OSTI2021

3-D Geologic Controls of Hydrothermal Fluid Flow at Brady Geothermal Field, Nevada using PCA

In many hydrothermal systems, fracture permeability along faults provides pathways for groundwater to transport heat from depth. Faulting generates a range of deformation styles that cross-cut heterogeneous geology, resulting in complex patterns of permeability, porosity, and hydraulic conductivity. Vertical connectivity (a through going network of permeable areas that allows advection of heat from depth to the shallow subsurface) is rare and is confined to relatively small volumes that have highly variable spatial distribution. This local compartmentalization of connectivity represents a significant challenge to understanding hydrothermal circulation and for exploring, developing, and managing hydrothermal resources. Here, we present an evaluation of the geologic characteristics that control this compartmentalization in hydrothermal systems through 3-D analysis of the Brady geothermal field in western Nevada. A published 3-D geologic map of the Brady area is used as a basis to develop structural and geological variables that are hypothesized to control or effect permeability or connectivity. The 3-D distribution of these variables is compared to the distribution of productive and non-productive fluid flow intervals along production wells and non-productive wells via principal component analysis (PCA). This comparison elucidates which geologic and structural variables are most closely associated with productive fluid flow intervals. Results indicate that production intervals at Brady are located: (1) within or near to known and stress-loaded macro-scale faults, and (2) in areas of high fault and fracture density. This submission includes the published journal article detailing this work, the published 3-D geologic map of the Brady Geothermal Area used as a basis to develop structural and geological variables that are hypothesized to control or effect permeability or connectivity, 3-D well data, along which geologic data were sampled for PCA analyses, and associated metadata file. This work was done using existing R programs.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Unsupervised learning of ferroic variants from atomically resolved STEM images

An approach for the analysis of atomically resolved scanning transmission electron microscopy data with multiple ferroic variants in the presence of imaging non-idealities and chemical variabilities based on a rotationally invariant variational autoencoder (rVAE) is presented. We show that an optimal local descriptor for the analysis is a sub-image centered at specific atomic units, since materials and microscope distortions preclude the use of an ideal lattice as a reference point. The applicability of unsupervised clustering and dimensionality reduction methods is explored and is shown to produce clusters dominated by chemical and microscope effects, with a large number of classes required to establish the presence of rotational variants. Comparatively, the rVAE allows extraction of the angle corresponding to the orientation of ferroic variants explicitly, enabling straightforward identification of the ferroic variants as regions with constant or smoothly changing latent variables and sharp orientational changes. This approach allows further exploration of the chemical variability by separating the rotational degrees of freedom via rVAE and searching for remaining variability in the system. The code used in this article is available at https://github.com/saimani5/ferroelectric_domains_rVAE.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2021

Identifying Periodic Variable Stars and Eclipsing Binary Systems with Long-term Las Cumbres Observatory Photometric Monitoring of ZTF J0139+5245

We present the results of our search for variable stars using the long-term Las Cumbres Observatory (LCO) monitoring of white dwarf ZTF J0139+5245 with the two 1.0 m telescope nodes located at McDonald Observatory using the Sinistro imaging instrument. In this search, we find 38 variable sources, of which 27 are newly discovered or newly classified (71%) based on comparisons with previously published catalogs, thereby increasing the number of detections in the field of view under consideration by a factor of ≈2.5. We find that the improved photometric precision per exposure due to longer exposure time for LCO images combined with the greater time sampling of LCO photometry enables us to increase the total number of detections in this field of view. Each LCO image covers a field of view of 26' × 26' and observes a region close to the Galactic plane (b = -9.°4) abundant in stars with an average stellar density of ≈8 arcmin{sup -2}. We perform aperture photometry and Fourier analysis on over 2000 stars across 1560 LCO images spanning 537 days to find 28 candidate BY Draconis variables, three candidate eclipsing binaries of type EA, and seven candidate eclipsing binaries of type EW. In assigning preliminary classifications to our detections, we demonstrate the applicability of the Gaia color–magnitude diagram as a powerful classification tool for variable-star studies.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Towards data-driven constitutive modelling for granular materials via micromechanics-informed deep learning

The analytical description of path-dependent elastic-plastic responses of a granular system is highly complicated because of continuously evolving microstructures and strain localisation within the system undergoing deformation. This study offers an alternative to the current analytical paradigm by developing micromechanics-informed machine-learning based constitutive modelling approaches for granular materials. A set of critical variables associated with the constitutive behaviour of granular materials are identified through an incremental stress-strain relationship analysis. Depending on the strategy to exploit the priori micromechanical knowledge, three different training strategies are explored. The first model uses only the measurable external variables to make stress predictions; the second model utilises a directed graph to link all the external strain sequences and internal microstructural evolution variables into a single prediction model comprised of a series of sub-mappings, and the third model explicitly integrates the physically important non-temporal properties with external strain paths into training through an enhanced Gated Recurrent Unit (GRU). These three models show satisfactory agreement with unseen test specimens based on multi-directional loading cases. The basic features and potential applications of each model are explained. Lastly, the key factors for constitutive training and limitations of the current work are also discussed in detail.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Modeling Variable Renewable Energy and Storage in the Power Sector

The emergence of variable renewable energy and battery storage technologies have fundamentally transformed the electric power sector and generated demand for analysis to understand their roles in future energy systems. Although unique characteristics of these resources are well-recognized and require more sophisticated methodologies to capture effectively, guidance is limited on best practices and research gaps. This paper selectively reviews recent literature and draws upon our collective modeling experience to offer recommendations to analysts and consumers of model outputs on approaches for modeling variable renewable energy and storage in long-term electric sector models. We focus on regional- and national-scale models with technological, temporal, and spatial detail given their prevalence in planning and policy analysis, though insights are applicable in other settings. The review highlights how the research frontier has advanced in representing renewables and energy storage over the past decade; however, given the many considerations involved with appropriately capturing salient economic and operational characteristics of renewables, there is a gap between commonly used models and state-of-the-art methods. Model simplifications can materially impact policy analysis associated with power sector decarbonization and high renewables deployment, and improving model representations of variable renewables can enhance insights for policymakers and other stakeholders. This review can point the way for improved methods for established models and designs for emerging ones.

energy and environmental policy↗
DOE OSTI2022

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

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

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

System Analysis of an Internal Combustion Engine (ICE) - Solid Oxide Fuel Cell (SOFC) Hybrid Cycle

The variability of renewable energy sources poses challenges for reliable grid operation. Conventional thermal power sources, though reliable, often lack operational flexibility. Hybrid energy systems that integrate Solid Oxide Fuel cells (SOFC) with Internal Combustion Engines (ICE) offer a promising solution by achieving relatively higher efficiency and grid-following capability. This study investigates performance of a 100-kW pressurized SOFC-ICE hybrid cycle. In this configuration, unutilized SOFC fuel is used to drive the engine, with a turbocharger providing air supply and an external reformer generating syngas. System components were numerically modeled using MATLAB/SIMULINK for the SOFC and reformer, and EBSILON® for the ICE and balance of plant. Parametric studies varied fuel utilization (70-90%), reformer temperature (600 – 1000K), anode off-gas recirculation (0 – 70%), and current density (0.2 – 0.55 A/cm2). Results show that the SOFC and ICE operate as thermally independent topping and bottoming cycles, achieving peak efficiency of 62% under optimized conditions.

hybrid↗
DOE OSTI2020

Probabilistic Look-ahead Contingency Analysis Integration with Commercial Tool and Practical Data

This paper presents an initial effort of integrating a smart sampling-based probabilistic look-ahead contingency analysis algorithm with General Electric (GE) Grid Solutions’ commercial energy management system (EMS) tool as a proof-of-concept for a seamless research tool integration using real world large-scale grid data. With the increasing impact of random forces such as variable generation and load, their stochastic behaviors cannot be ignored. However, the current practices are still dominated by deterministic tools. They are becoming increasingly inadequate for the future grid. The developed look-ahead contingency analysis algorithm incorporates forecast errors of variable energy and load to address the challenges brought by the increasing uncertainty of power system. The algorithm can reveal the potential violations caused by the variance of variable energy and load that are not normally detected by traditional deterministic approaches. To test its performance under practical environments ( real data with real commercial tool), significant efforts have been made to prepare test cases, modify GE EMS tool, and adapt an extreme value distribution algorithm to analyze the GE EMS’s violation-only outputs. The test results clearly demonstrate the effectiveness of the developed algorithm as new transformer violations that were not previously detected have been identified. This performance provides better situational awareness to engineers for their decision-making process under uncertainty. Moreover, with the discussion of computational performance and future work, this paper has shown a clear path for integrating the probabilistic algorithm with commercial tools to make us better equipped for the changing power system.

Modeling and simulation of power systems, constrai↗
DOE OSTI2024

Comparative Analyses of Bioequivalence Assessment Methods for In Vitro Permeation Test Data

ABSTRACT For topical, dermatological drug products, an in vitro option to determine bioequivalence (BE) between test and reference products is recommended. In particular, in vitro permeation test (IVPT) data analysis uses a reference‐scaled approach for two primary endpoints, cumulative penetration amount (AMT) and maximum flux ( J max ), which takes the within donor variability into consideration. In 2022, the Food and Drug Administration (FDA) published a draft IVPT guidance that includes statistical analysis methods for both balanced and unbalanced cases of IVPT study data. This work presents a comprehensive evaluation of various methodologies used to estimate critical parameters essential in assessing BE. Specifically, we investigate the performance of the FDA draft IVPT guidance approach alongside alternative empirical and model‐based methods utilizing mixed‐effects models. Our analyses include both simulated scenarios and real‐world studies. In simulated scenarios, empirical formulas consistently demonstrate robustness in approximating the true model, particularly in effectively addressing treatment–donor interactions. Conversely, the effectiveness of model‐based approaches heavily relies on precise model selection, which significantly influences their results. The research emphasizes the importance of accurate model selection in model‐based BE assessment methodologies. It sheds light on the advantages of empirical formulas, highlighting their reliability compared to model‐based approaches and offers valuable implications for BE assessments. Our findings underscore the significance of robust methodologies and provide essential insights to advance their understanding and application in the assessment of BE, employed in IVPT data analysis.

Leon, Sami↗
DOE OSTI2023

Causal interaction in high frequency turbulence at the biosphere–atmosphere interface: Structure–function coupling

At the biosphere–atmosphere interface, nonlinear interdependencies among components of an ecohydrological complex system can be inferred using multivariate high frequency time series observations. Information flow among these interacting variables allows us to represent the causal dependencies in the form of a directed acyclic graph (DAG). Here, we use high frequency multivariate data at 10 Hz from an eddy covariance instrument located at 25 m above agricultural land in the Midwestern US to quantify the evolutionary dynamics of this complex system using a sequence of DAGs by examining the structural dependency of information flow and the associated functional response. We investigate whether functional differences correspond to structural differences or if there are no functional variations despite the structural differences. We base our analysis on the hypothesis that causal dependencies are instigated through information flow, and the resulting interactions sustain the dynamics and its functionality. To test our hypothesis, we build upon causal structure analysis in the companion paper to characterize the information flow in similarly clustered DAGs from 3-min non-overlapping contiguous windows in the observational data. We characterize functionality as the nature of interactions as discerned through redundant, unique, and synergistic components of information flow. Through this analysis, we find that in turbulence at the biosphere–atmosphere interface, the variables that control the dynamic character of the atmosphere as well as the thermodynamics are driven by non-local conditions, while the scalar transport associated with CO and H 2 O is mainly driven by short-term local conditions.

58 GEOSCIENCES↗
DOE OSTI2020

Imprecise global sensitivity analysis using bayesian multimodel inference and importance sampling

Global Sensitivity Analysis (GSA) aims to understand the relative importance of uncertain input variables to model response. Conventional GSA involves calculating sensitivity (Sobol’) indices for a model with known model parameter distributions. However, model parameters are affected by aleatory and epistemic uncertainty, with the latter often caused by lack of data. In this paper, we propose a new framework to quantify uncertainty in probability model-form and model parameters resulting from small datasets and integrate these uncertainties into Sobol’ index estimates. First, the process establishes, through Bayesian multimodel inference, a set of candidate probability models and their associated probabilities. Imprecise Sobol’ indices are calculated from these probability models using an importance sampling reweighting approach. This results in probabilistic Sobol’ indices, whose distribution characterizes uncertainty in the sensitivity resulting from small dataset size. The imprecise Sobol’ indices thus provide a measure of confidence in the sensitivity estimate and, moreover, can be used to inform data collection efforts targeted to minimize the impact of uncertainties. Through an example studying the parameters of a Timoshenko beam, we show that these probabilistic Sobol’ indices converge to the true/deterministic Sobol’ indices as the dataset size increases and hence, distribution-form uncertainty reduces. The approach is then applied to assess the sensitivity of the out-of-plane properties of an E-glass fiber composite material to its constituent properties. This second example illustrates the approach for an important class of materials with wide-ranging applications when data may be lacking for some input parameters.

42 ENGINEERING↗