Search NASA⌕ Search

SEARCH · Search NASA

Results for “Input Sensitivity”

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 829 records · Page 46

Employing Sensitivity Derivatives for Robust Optimization under Uncertainty in CFD

A robust optimization is demonstrated on a two-dimensional inviscid airfoil problem in subsonic flow. Given uncertainties in statistically independent, random, normally distributed flow parameters (input variables), an approximate first-order statistical moment method is employed to represent the Computational Fluid Dynamics (CFD) code outputs as expected values with variances. These output quantities are used to form the objective function and constraints. The constraints are cast in probabilistic terms; that is, the probability that a constraint is satisfied is greater than or equal to some desired target probability. Gradient-based robust optimization of this stochastic problem is accomplished through use of both first and second-order sensitivity derivatives. For each robust optimization, the effect of increasing both input standard deviations and target probability of constraint satisfaction are demonstrated. This method provides a means for incorporating uncertainty when considering small deviations from input mean values.

Newman, Perry A.↗

Energy input and HI spin temperatures in low pressure regions

Two recent (unpublished) HI emission/absorption studies carried out with good sensitivity using the Arecibo 21 cm beam are discussed. One study (Colgan, Salpeter and Terzian) looked for high velocity clouds of our own Galaxy in absorption in the directions of 63 of the brightest continuum sources reachable with the Arecibo telescope. HI emission mapping in the neighborhood of these directions was also carried out. The other study (Corbelli and Schneider) looked for absorption along lines of sight to about 50 weaker sources which pass within a few diameters of nearby disk galaxies. Neither study detected any absorption.

Corbelli, E.↗

Local and Global Sensitivity Analysis of a Reactive Transport Model Simulating Floodplain Redox Cycling

Reactive transport models (RTMs) are essential tools that simulate the coupling of advective, diffusive, and reactive processes in the subsurface, but their complexity makes them difficult to understand, develop and improve without accompanying statistical analyses. Although global sensitivity analysis (SA) can address these issues, the computational cost associated with most global SA techniques limits their use with RTMs. In this study, we apply distance-based generalized sensitivity analysis (DGSA), a novel and computationally efficient method of global SA, to a floodplain-scale RTM and compare DGSA results to those from local SA. Our test case focuses on the impact of 17 uncertain environmental parameters on spatially and temporally variable redox conditions within a floodplain aquifer. The input parameters considered include flow and diffusion rates, geochemical reaction rates, and the spatial distribution of sediment facies. Sensitivity was evaluated for three distinct components of the model response, encompassing both multidimensional and categorical output. Parameter rankings differ between local SA and DGSA, due to nonlinear effects of individual parameters and interaction effects between parameters. DGSA results show that fluid residence time, which is controlled by aquifer permeability, generally exerts a stronger control on redox conditions than do geochemical reaction rates. Sensitivity indices also demonstrate that sulfate reduction is key for establishing and maintaining reducing conditions throughout the aquifer. Furthermore, these results provide insights into the key drivers of heterogeneous redox processes within floodplain aquifers, as well as the main sources of uncertainty when modeling complex subsurface systems.

54 ENVIRONMENTAL SCIENCES↗

Field Sensitivity Analysis of Turbulence Model Parameters for Flow Over a Wing

Reynolds-averaged-Navier-Stokes (RANS) turbulence models are a critical tool in computational-fluid-dynamics simulations of aerodynamic systems, but simulation results can be highly sensitive to RANS-model parameter choices. Sensitivity analysis can be used to quantify these impacts, and the objective of this study is to demonstrate field sensitivity analysis with respect to ten parameters in the 2003 Menter shear-stress-transport (SST) turbulence model. Here, the analysis is demonstrated for an application relevant to wind energy, namely, flow over a NACA 0015 wing at 12 degree angle of attack and a Reynolds number of 1.5 x 10 6 . We quantify sensitivity using Sobol indices and the mean-squared gradient, which are estimated using polynomial chaos and active subspace models, respectively. Our results indicate that there are substantial spatial variations in parameter sensitivities, with different sets of most-sensitive parameters near the wing, as well as in the downstream wake, consistent with the physical interpretations of the turbulence model inputs. We show that, for this particular turbulence model and flow, simultaneous dimension reduction is possible across all quantities of interest, enabling efficient exploration of model outcomes. Ultimately, this analysis provides new insights into turbulence model parameter sensitivities in incompressible flows, and also demonstrates the implementation of field sensitivity analysis for applications relevant to aerodynamics simulations.

17 WIND ENERGY↗

Coarse woody debris decomposition assessment tool: Model development and sensitivity analysis

Coarse woody debris (CWD) is an important component in forests, hosting a variety of organisms that have critical roles in nutrient cycling and carbon (C) storage. We developed a process-based model using literature, field observations, and expert knowledge to assess woody debris decomposition in forests and the movement of wood C into the soil and atmosphere. The sensitivity analysis was conducted against the primary ecological drivers (wood properties and ambient conditions) used as model inputs. The analysis used eighty-nine climate datasets from North America, from tropical (14.2° N) to boreal (65.0° N) zones, with large ranges in annual mean temperature (26.5°C in tropical to -11.8°C in boreal), annual precipitation (6,143 to 181 mm), annual snowfall (0 to 612 kg m -2 ), and altitude (3 to 2,824 m above mean see level). The sensitivity analysis showed that CWD decomposition was strongly affected by climate, geographical location and altitude, which together regulate the activity of both microbial and invertebrate wood-decomposers. CWD decomposition rate increased with increments in temperature and precipitation, but decreased with increases in latitude and altitude. CWD decomposition was also sensitive to wood size, density, position (standing vs downed), and tree species. The sensitivity analysis showed that fungi are the most important decomposers of woody debris, accounting for over 50% mass loss in nearly all climatic zones in North America. The model includes invertebrate decomposers, focusing mostly on termites, which can have an important role in CWD decomposition in tropical and some subtropical regions. The role of termites in woody debris decomposition varied widely, between 0 and 40%, from temperate areas to tropical regions. Woody debris decomposition rates simulated for eighty-nine locations in North America were within the published range of woody debris decomposition rates for regions in northern hemisphere from 1.6° N to 68.3° N and in Australia.

54 ENVIRONMENTAL SCIENCES↗

Estimation of 3D Woven Design Sensitivities Using a Rapid Multiscale Analysis Technique

Highly-refined finite element models of three-dimension (3D) woven composite systems currently require excessive computational demands that limit their use in sensitivity analysis, uncertainty quantification, and optimization. An alternative analysis methodology was developed using the NASA Multiscale Analysis Tool (NASMAT) where multiscale models of a 3D woven composite (including inter-tow matrix voids and constituent failure) can be completed on a single central processing unit (CPU) on the order of ~30s. To develop inputs and validation data for the NASMAT model, coupon and acid-digesting testing and x-ray computed tomography were performed. The NASMAT inputs were parameterized using a set of 25 input variables and distributions. These inputs were randomly sampled to generate a total of 100,000 NASMAT analyses that could be used to understand the influence of different material and geometric properties on the warp and weft-direction stiffness and strength. These analyses (including pre/post-processing) were performed in less than eight hours on a 120 CPU cluster. The computational efficiency of the NASMAT model enabled a sensitivity analysis to be performed, and dominant input variables were able to be identified. Key results were consistent with theoretical and experimental observations for the specific 3D woven system studied in this work.

NASMAT↗

Dynamic, Intermediate Soil Carbon Pools May Drive Future Responsiveness to Environmental Change

Accurately capturing dynamic soil response to disturbance effects in agroecosystem models remains elusive, thereby limiting projections of climate change mitigation potential. Perennial grasses cultivated in zero‐tillage management systems hold promise as sustainable agroecosystems. High‐yielding tropical C 4 grasses often have extensive rooting systems, and the belowground processes of root turnover, aggregate formation, and mineral stabilization drove rapid C accumulation after cultivation in a recent study. We sought (i) to understand and constrain the size and responsiveness of dynamic, intermediate‐cycling C pools contributing to the observed C accrual rates, and (ii) to simulate C stocks over time under the disturbance of elevated temperature using soil incubation at multiple temperatures and physical fractionation via density and sonication. Three‐pool transfer modeling of soil incubations revealed small pools of readily available (i.e., days to months) microbial substrate that were responsive to temperature, time since cultivation, and inputs. Larger, kinetically slow‐cycling pools were more indicative of long‐term (i.e., years to decades) changes in C stock and strongly connected to measured changes in physical fractions. Combining the sensitivity of readily available microbial substrate with three‐pool transfer modeling of the physical fractions over time since cultivation revealed that dynamic transfers of inputs occurred between the free organic and aggregate‐protected fractions, and from these fractions to the mineral‐associated dense fraction. Under 5°C temperature elevation, increased transfer rates outweighed elevated decomposition losses to sustain soil C accrual into the future. To effectively plan managed landscapes and monitor sustainable agroecosystems for climate change mitigation, tools must incorporate the complexity of soil response to change. Core Ideas Dynamic aspects of soil C are inadequately represented in most models. Intermediate‐cycling pools are more responsive to disturbances than thought. Root inputs are rapidly moved from particulate debris to mineral‐bound stable pools. Temperature increased the transfer rate to stable pools and sustained C accumulation Climate change mitigation projections must capture dynamic responses in agroecosystems.

Crow, Susan E.↗

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

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

36 MATERIALS SCIENCE↗

Utilizing In Situ Directional Hyperspectral Measurements to Validate Bio-Indicator Simulations for a Corn Crop Canopy

Two radiative transfer canopy models, SAIL and the two-layer Markov-Chain Canopy Reflectance Model (MCRM), were coupled with in situ leaf optical properties to simulate canopy-level spectral band ratio vegetation indices with the focus on the photochemical reflectance index in a cornfield. In situ hyperspectral measurements were made at both leaf and canopy levels. Leaf optical properties were obtained from both sunlit and shaded leaves. Canopy reflectance was acquired for eight different relative azimuth angles (psi) at three different view zenith angles (Theta (sub v)), and later used to validate model outputs. Field observations of photochemical reflectance index (PRI) for sunlit leaves exhibited lower values than shaded leaves, indicating higher light stress. Canopy PRI expressed obvious sensitivity to viewing geometry, as a function of both Theta (sub v) and psi . Overall, simulations from MCRM exhibited better agreements with in situ values than SAIL. When using only sunlit leaves as input, the MCRM-simulated PRI values showed satisfactory correlation and RMSE, as compared to in situ values. However, the performance of the MCRM model was significantly improved after defining a lower canopy layer comprised of shaded leaves beneath the upper sunlit leaf layer. Four other widely used band ratio vegetation indices were also studied and compared with the PRI results. MCRM simulations were able to generate satisfactory simulations for these other four indices when using only sunlit leaves as input; but unlike PRI, adding shaded leaves did not improve the performance of MCRM. These results support the hypothesis that the PRI is sensitive to physiological dynamics while the others detect static factors related to canopy structure. Sensitivity analysis was performed on MCRM in order to better understand the effects of structure related parameters on the PRI simulations. Leaf area index (LAI) showed the most significant impact on MCRM-simulated PRI among the parameters studied. This research shows the importance of hyperspectral and narrow band sensor studies, and especially the necessity of including the green wavelengths (e.g., 531 nm) on satellites proposing to monitor carbon dynamics of terrestrial ecosystems.

Cheng, Yen-Ben↗

Estimation of 3D Woven Design Sensitivities Using a Rapid Multiscale Analysis Technique

Highly-refined finite element models of three-dimension (3D) woven composite systems currently require excessive computational demands that limit their use in sensitivity analysis, uncertainty quantification, and optimization. An alternative analysis methodology was developed using the NASA Multiscale Analysis Tool (NASMAT) where multiscale models of a 3D woven composite (including inter-tow matrix voids and constituent failure) can be completed on a single central processing unit(CPU)on the order of ~30 s. To develop inputs and validation data for the NASMAT model, coupon and acid-digesting testing and x-ray computed tomography were performed. The NASMAT inputs were parameterized using a set of 25 input variables and distributions. These inputs were randomly sampled to generate a total of 100,000 NASMAT analyses that could be used to understand the influence of different material and geometric properties on the warp and weft-direction stiffness and strength. These analyses (including pre/post-processing) were performed in less than eight hours on a 120 CPU cluster. The computational efficiency of the NASMAT model enabled a sensitivity analysis to be performed, and dominant input variables were able to be identified. Key results were consistent with theoretical and experimental observations for the specific 3D woven system studied in this work.

NASMAT↗

A Fresh Take: Seasonal Changes in Terrestrial Freshwater Inputs Impact Salt Marsh Hydrology and Vegetation Dynamics

Abstract Salt marshes exist at the terrestrial-marine interface, providing important ecosystem services such as nutrient cycling and carbon sequestration. Tidal inputs play a dominant role in salt marsh porewater mixing, and terrestrially derived freshwater inputs are increasingly recognized as important sources of water and solutes to intertidal wetlands. However, there remains a critical gap in understanding the role of freshwater inputs on salt marsh hydrology, and how this may impact marsh subsurface salinity and plant productivity. Here, we address this knowledge gap by examining the hydrologic behavior, porewater salinity, and pickleweed ( Sarcocornia pacifica also known as Salicornia pacifica) plant productivity along a salt marsh transect in an estuary along the central coast of California. Through the installation of a suite of hydrometric sensors and routine porewater sampling and vegetation surveys, we sought to understand how seasonal changes in terrestrial freshwater inputs impact salt marsh ecohydrologic processes. We found that salt marsh porewater salinity, shallow subsurface saturation, and pickleweed productivity are closely coupled with elevated upland water level during the winter and spring, and more influenced by tidal inputs during the summer and fall. This seasonal response indicates a switch in salt marsh hydrologic connectivity with the terrestrial upland that impacts ecosystem functioning. Through elucidating the interannual impacts of drought on salt marsh hydrology, we found that the severity of drought and historical precipitation can impact contemporary hydrologic behavior and the duration and timing of the upland-marsh hydrologic connectivity. This implies that the sensitivity of salt marshes to climate change involves a complex interaction between sea level rise and freshwater inputs that vary at seasonal to interannual timescales.

54 ENVIRONMENTAL SCIENCES↗

A Sensitivity-informed, Process Model Abstraction Implementation in GoldSim - 19276

System-level Performance Assessment (PA) models rely on specific process-level models for some inputs, such as infiltration. For probabilistic PA, the value of the process model output depends on the values of input parameters that also vary in the PA model. The simplest conceptual solution is to draw values of the input parameters from their distributions, run the process level model at these values, and then use the output where it is needed in the PA model. However, this solution is unsatisfactory for multiple reasons: (1) the PA Model should be set up such that new runs can be made and distributions can change without having to go back and run the process-level model; (2) calling the process model from the system model is possible, but technical difficulties remain in coupling the two; and (3) there is no uncertainty in the process-level model results incorporated into the PA model. Therefore, it is desirable to develop a method that maps the values of PA Model inputs to process level model outputs, in real time within the PA Model set up and using as simple of a map as possible without losing key features of the relationship. This has been termed 'model abstraction' because information about process model runs is 'abstracted' to be used with the PA model. A simple 'nearest neighbor' abstraction approach is presented to overcome these problems, in which a sensitivity-informed distance metric is used to map a system-level model realization to the most appropriate process-level model realization. (authors)

97 MATHEMATICS AND COMPUTING↗

Machine learning-assisted early ignition prediction in a complex flow

Machine learning methods are used to improve the efficiency by which turbulence-resolved simulations predict whether a hydrogen jet in air crossflow will successfully ignite. The flush-mounted jet issues perpendicularly from the wall of a low-speed wind tunnel into a turbulent boundary layer wherein a laser-induced optical breakdown (LIB) hotspot is deposited. A detailed hydrogen chemical mechanism is used to model the radicals and any subsequent chemical reactions. A dielectric-barrier discharge actuator generates body forces and hydrogen radicals near the jet orifice. Here we focus on the success or not of the ignition based on LIB location. A challenge is that definitive determination of this requires long simulations, up to 440 µs after the LIB deposition. This is particularly expensive since multiple simulations are required to find the threshold. To reduce the computational effort, three short-time (91 µs) criteria are proposed, evaluated, and compared: a constructed criterion based on detailed observations of radicals near the stoichiometric surface and two machine learning approaches, each trained on 38 realizations. The constructed criterion provides a low-cost estimate of the ignition boundary that is unambiguous in 45 of the 50 training and test trials, so only 10% of them would need to be simulated longer. The trained neural networks correctly predict outcomes in all cases evaluated, with the more automated procedure—a convolutional neural network (CNN) trained on two-dimensional images—providing the most definitive outcome prediction. From the CNN, a sensitivity analysis is used to determine which kernel features from the two-dimensional input data, as defined by intermediate-layer network weights, are important for identifying ignition.

42 ENGINEERING↗

Sediment Characteristics and Methane Ebullition in Three Subarctic Lakes

Ebullition (bubbling) from climate-sensitive northern lakes remains an unconstrained source of atmospheric methane (CH 4 ). Although the focus of many recent studies, ebullition is rarely linked to the physical characteristics of lakes. In this study we analyze the sediments of subarctic postglacial lakes and investigate how sediment properties relate to the large spatial variation in CH 4 bubble flux, quantified over multiple years using bubble traps. The results show that the sediments from our lakes are rich in total organic carbon, containing 37 kg/m 3 on average. This number is roughly 40% higher than the average for yedoma deposits, which have been identified as high CH 4 emitters. However, the quantity of total organic carbon is not a useful indicator of high emissions from the study lakes. Neither is the amount of CH 4 in the sediment a reliable measure of ebullition potential. Instead, our data point to coarse detritus, partly from buried submerged aquatic vegetation and redeposited peat as spatial controls on fluxes, often in combination with previously established effects of incoming solar radiation and water depth. The results once again highlight the climate sensitivity of northern lakes, indicating that biological responses to warmer waters and increased energy input and heating of organic sediments during longer ice-free seasons can substantially alter future CH 4 emissions.

58 GEOSCIENCES↗

Impacts of anthropogenic aerosols on fog in North China Plain

Fog poses a sever environmental problem in the North China Plain (NCP), China, which has been witnessing increases anthropogenic emission since the early 1980s. Here, this work first uses the WRF/Chem model coupled with the local anthropogenic emissions to simulate and evaluate a severe fog event occurring in NCP. Comparison of the simulations against observations shows that WRF/Chem well reproduces the general features of temporal evolution of PM 2.5 mass concentration, fog spatial distribution, visibility, and vertical profiles of temperature, water vapor content and relative humidity in the planetary boundary layer throughout the whole period of the fog event. Sensitivity studies are then performed with five different levels of anthropogenic emission as model inputs to systematically examine the comprehensive impacts of aerosols on fog microphysical, macrophysical, radiative, and dynamical properties. The results show that as aerosol concentration increases, fog droplet number concentration and liquid water content all increase nonlinearly; but effective radius decreases. Macrophysical properties (fog fraction, fog duration, fog height, and liquid water path) also increase nonlinearly with increasing aerosol concentration, with rates of changes smaller than microphysical properties. Further analysis reveals distinct aerosol effects on thermodynamic and dynamical conditions during different stages of fog evolution: increasing aerosols invigorate fog formation and development by enhancing longwave-induced instability, fog droplet condensation accompanying latent heat release, and thus turbulence, but delay fog dissipation by reducing surface solar radiation, surface sensible and latent heat fluxes, and thus suppressing turbulence during the dissipation stage.

54 ENVIRONMENTAL SCIENCES↗

Distribution–Based Global Sensitivity Analysis in Hydrology

Global sensitivity analysis (GSA) is routinely used in academic setting to quantify the influence of input variability and uncertainty on predictions of a quantity of interest. Practical applications of GSA are hampered by its high computational cost, which arises from the need to run large (e.g., groundwater) models multiple times, and by its reliance on the analysis of variance, which formally requires input parameters to be uncorrelated. The former difficulty can be alleviated by replacing expensive models with inexpensive (e.g., polynomial) surrogates, while adoption of distribution-based (rather than variance-based) metrics can, in principle, overcome the latter but at significantly increased computational cost. To make use of distribution-based GSA feasible for regional-scale models with a large number of degrees of freedom, we supplement it with a surrogate model built with polynomial chaos expansions with analytically updated coefficients. Here, we demonstrate the computational efficiency of our algorithm on a case study dealing with evaluation of the effects of temperature variability on annual evapotranspiration at the regional scale.

42 ENGINEERING↗

Uncertainty Quantification Study of CTF for the OECD/NEA LWR Uncertainty Analysis in Modeling Benchmark

This work describes the results of a quantitative uncertainty analysis of the thermal-hydraulic subchannel code for nuclear engineering applications, Coolant Boiling in Rod Arrays-Three Field (COBRA-TF). CTF is used, which is a version of COBRA-TF developed in cooperation between the Consortium for Advanced Simulation of Light Water Reactors and North Carolina State University. Four steady-state cases from Phase II Exercise 3 of the Organisation for Economic Co-operation and Development/Nuclear Energy Agency Light Water Reactor Uncertainty Analysis in Modeling (UAM) Benchmark are analyzed using the statistical analysis tool, Design Analysis Kit for Optimization and Terascale Applications (Dakota). The input parameters include boundary condition, geometry, and modeling uncertainties, which are selected using a sensitivity study and then defined based on expert judgment. Here, a forward uncertainty quantification method with Latin hypercube sampling (LHS) is used, where the sample size is based on available computational resources.

42 ENGINEERING↗

Effects of Photovoltaic Module Materials and Design on Module Deformation Under Load

Quasi-static structural finite-element models of an aluminum-framed crystalline silicon photovoltaic module and a glass-glass thin-film module were constructed and validated against experimental measurements of deflection under uniform pressure loading. Specific practices in the computational representation of module assembly were identified as influential to matching experimental deflection observations. Additionally, parametric analyses using Latin hypercube sampling were performed to propagate input uncertainties related to module materials, dimensions, and tolerances into uncertainties in simulated deflection. Sensitivity analyses were performed on the uncertainty quantification datasets using linear correlation coefficients and variance-based sensitivity indices to elucidate key parameters influencing module deformation. Results identified edge tape and adhesive material properties as being strongly correlated to module deflection, suggesting that optimization of these materials could yield module stiffness gains at par with the conventionally structural parameters, such as glass thickness. This exercise verifies the applicability of finite-element models for accurately predicting mechanical behavior of solar modules and demonstrates a workflow for model-based parametric uncertainty quantification and sensitivity analysis. Finally, applications of this capability include the assessment of field environment loads, derivation of representative loading conditions for reduced-scale testing, and module design optimization, among others.

42 ENGINEERING↗