Power System Protection Parameter Sensitivity Analysis with Integrated Inverter Based Resources.
Abstract not provided.
SEARCH · Engineering Papers
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.
Abstract not provided.
Tensor decomposition models play an increasingly important role in modern data science applications. One problem of particular interest is fitting a low-rank Canonical Polyadic (CP) tensor decomposition model when the tensor has sparse structure and the tensor elements are nonnegative count data. SparTen is a high-performance C++ library which computes a low-rank decomposition using different solvers: a first-order quasi-Newton or a second-order damped Newton method, along with the appropriate choice of runtime parameters. Since default parameters in SparTen are tuned to experimental results in prior published work on a single real-world dataset conducted using MATLAB implementations of these methods, it remains unclear if the parameter defaults in SparTen are appropriate for general tensor data. Furthermore, it is unknown how sensitive algorithm convergence is to changes in the input parameter values. This report addresses these unresolved issues with large-scale experimentation on three benchmark tensor data sets. Experiments were conducted on several different CPU architectures and replicated with many initial states to establish generalized profiles of algorithm convergence behavior.
Abstract not provided.
Abstract not provided.
Sensitivity analysis and reliability assessment are two important aspects of structural and system safety. Epistemic uncertainty with respect to probabilistic model of input parameters due to lack of knowledge is present in many scarce-data applications and complicates the characterization of uncertainty in model response. In this article, we present two importance measures to evaluate the impact of distribution parameters on the probability distribution function (PDF) of the output and the failure probability. The epistemic uncertainty associated with the distribution parameters is modeled as random variables. Additionally, a modified extended polynomial chaos expansion (MEPCE) approach is introduced in which aleatory and epistemic random variables are modeled and propagated simultaneously while allowing the separate assessment for any single epistemic variable. A MEPCE-based kernel density estimation (KDE) construction provides a composite map from each epistemic variable to the response PDF. The functional global sensitivity index of the PDF with respect to the distribution parameters is thus derived, as a function of output, which is both more informative and more efficient than standard scalar sensitivity measures. Reliability sensitivity indices can be readily evaluated by integrating the global sensitivity index function over the failure zone. Three illustrative examples are used to demonstrate the proposed methodology.
Many problems in engineering and sciences require the solution of large scale optimization constrained by partial differential equations (PDEs). Though PDE-constrained optimization is itself challenging, most applications pose additional complexity, namely, uncertain parameters in the PDEs. Uncertainty quantification (UQ) is necessary to characterize, prioritize, and study the influence of these uncertain parameters. Sensitivity analysis, a classical tool in UQ, is frequently used to study the sensitivity of a model to uncertain parameters. In this article, we introduce "hyper-differential sensitivity analysis" which considers the sensitivity of the solution of a PDE-constrained optimization problem to uncertain parameters. Our approach is a goal-oriented analysis which may be viewed as a tool to complement other UQ methods in the service of decision making and robust design. We formally define hyper-differential sensitivity indices and highlight their relationship to the existing optimization and sensitivity analysis literatures. Assuming the presence of low rank structure in the parameter space, computational efficiency is achieved by leveraging a generalized singular value decomposition in conjunction with a randomized solver which converts the computational bottleneck of the algorithm into an embarrassingly parallel loop. Two multi-physics examples, consisting of nonlinear steady state control and transient linear inversion, demonstrate efficient identification of the uncertain parameters which have the greatest influence on the optimal solution.
The growing adoption of renewable energy is driving the prevalence of inverter-based resources (IBRs) within power grids. Future power grids will integrate both grid-following IBRs (GFM-IBRs) and grid-forming IBRs (GFL-IBRs) alongside synchronous generators. Therefore, it is crucial to perform stability studies that account for all components and especially control interactions related to IBRs. Extensive research has performed to study the IBR-related stability, however, the sensitivity study of IBRs' control parameters on system stability has not been adequately studied yet, especially from a systematic way. Therefore, this paper conducts a small signal stability analysis for a generic grid with multiple types of resources and develops an analytical framework for assessing the sensitivity of control parameters affecting stability margins. To achieve that, the non-autonomous reduced-order non-linear dynamic model is developed for a generic power system with multiple synchronous generator-based resources (SGBRs), GFM-IBRs, and GFL-IBRs. Based on the analytic model, a systematic framework for parametric sensitivity on systems' asymptotic stability is developed. A parameter sensitivity analysis based on eigenvalue methods is proposed. The impact of the droop controllers of GFM-IBRs, PQ-dispatch and the PLL controller of GFL-IBR on the system asymptotic stability is discussed. This sensitivity study is aiming to provide deep insights on control parameters' impact on system stability, and gives direction for parameter tuning in case of instability.
To address the need for an approach to evaluate cumulative ecological impacts of wind energy development, particularly those on critical wildlife habitats and species, the authors developed a proof-of-concept, landscape-based, spatially explicit individual-based modeling (IBM) framework for populations of the greater sage-grouse (Centrocercus urophasianus) in Albany County, Wyoming, based on published and other available information on the life history of the species. This sage-grouse IBM estimates the spatial and temporal movements of sage-grouse based on modeled effects of habitat characteristics and developments on sage-grouse condition, survivorship, and reproduction. After completing the initial version of the model, we conducted a series of sensitivity analyses of key model parameters and validation of lek occurrences and spatial distributions to better understand performance of the model relative to observed values of important life history characteristics. The parameter sensitivity analysis was conducted as a series of one-factor-at-a-time tests in which the values of selected parameters were individually varied to evaluate the corresponding effects on model output. We examined the model’s response to changes in parameter values to gain insight into the overall robustness of the model, identify key driver parameters, and initiate qualitative validation. Five functions selected for sensitivity analyses were (1) the rate of change in body condition in relationship to habitat suitability index (HSI) of occupied habitat (condition change rate), (2) the effect of competition between individuals with overlapping home ranges expressed as a percentage reduction in HSI (competition factor), (3) value of condition at which survivorship probability is equal to zero (zero survivorship condition), (4) the relationship between female body condition and clutch size, and (5) the relationship between female body condition and nest success. We varied the parameter values to be tested for the sensitivity analysis around the original value used in the model (i.e., tested values were higher and lower than the model’s original value) to examine the effect of this variation on the predicted output of the model. Model outputs evaluated in the sensitivity analysis included (1) predicted age-class distribution (as measured by the percentage of immature individuals within the population); (2) predicted sex ratio; (3) predicted age-sex class distribution (percent of population in each of six age-sex classes, i.e., male and female juveniles, yearlings, and adults); (4) lifespan; (5) population size; and (6) spatial distribution. Of the five parameters analyzed, condition change rate, survivorship, clutch size, and nest success did not substantially affect any of the evaluated model outputs. For the predicted age-class distribution, the percentage of immature individuals predicted by the model for these four parameters ranged from 61.9% to 65.0%, which was reasonable with respect to the reference values of 51.4 to 57.8% in published studies. The ratio of females to males ranged from 1.09 to 1.14, which was comparable to the reference values of 1.2 to 3.0 presented in an existing study. The predicted mean lifespan ranged from 1.62 to 1.74 years, which appears to be reasonable with respect to the reference values of 0.9 to 1.1 years for sharp-tail grouse and greater prairie chickens presented by published studies. The predicted population size ranged from 1,877 to 3,140. This prediction would be reasonable with respect to our estimate of 5,000 yearlings and adults, which was based on the USFWS estimate for the state of Wyoming scaled to the number of leks in the county. The predicted spatial distribution was not substantially affected by the parameter values examined. In our sensitivity analyses, the only tested parameter that had a noticeable effect on model output was competition factor. Over the range of values tested for this parameter, predicted population size exhibited the largest range (from 1,833 to 4,737) of all parameters tested. As the competition factor increased, the predicted population size decreased. This reflects the effect competition had on the number of individuals that could co occur in high-quality habitat patches. Despite this apparently significant effect, the range of values appears to be reasonable with respect to our estimate of 5,000 yearlings and adults, which was based on the USFWS estimate for the state of Wyoming scaled to the number of leks in the county.
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.
Phenological transitions determine the timing of changes in land surface properties and the seasonality of exchanges of biosphere-atmosphere energy, water, and carbon. Accurate mechanistic modeling of phenological processes is therefore critical to understand and correctly predict terrestrial ecosystem feedbacks with changing atmospheric and climate conditions. However, the phenological components in the land model of the US Department of Energy's (DOE) Energy Exascale Earth System Model (ELM of E3SM) were previously unable to accurately capture the observed phenological responses to environmental conditions in a well-studied boreal peatland forest. In this research, we introduced new seasonal-deciduous phenology schemes into version 1.0 of ELM and evaluated their performance against the PhenoCam observations at the Spruce and Peatland Responses Under Changing Environments (SPRUCE) experiment in northern Minnesota from 2015 to 2018. We found that phenology simulated by the revised ELM (i.e., earlier spring onsets and stronger warming responses of spring onsets and autumn senescence) was closer to observations than simulations from the original algorithms for both the deciduous conifer (Larix laricina) and mixed shrub layers. Moreover, the revised ELM generally produced higher carbon and water fluxes (e.g., photosynthesis and evapotranspiration) during the growing season and stronger flux responses to warming than the default ELM. A parameter sensitivity analysis further indicated the significant contribution of phenology parameters to uncertainty in key carbon and water cycle variables, underscoring the importance of precise phenology parameterization. Furthermore, this phenological modeling effort demonstrates the potential to enhance the E3SM representation of land-climate interactions at broader spatiotemporal scales, especially under anticipated elevated CO 2 and warming conditions.
A multiple input multiple output (MIMO) power line communication (PLC) model for industrial facilities was developed that uses the physics of a bottom-up model but can be calibrated like top-down models. The PLC model considers 4-conductor cables (three-phase conductors and a ground conductor) and has several load types, including motor loads. The model is calibrated to data using mean field variational inference with a sensitivity analysis to reduce the parameter space. The results show that the inference method can accurately identify many of the model parameters, and the model is accurate even when the network is modified.
Raw and processed output from the Multi-Assumption Architecture & Testbed (MAAT, https://github.com/walkeranthonyp/MAAT, tag: v1.2.1_walker2020; Walker et al., 2018). The dataset contains leaf-scale carbon assimilation data from three MAAT simulations and their post-processing and plotting scripts. Simulation 1, labelled "FvC_4proc_proc," is a four-process, process sensitivity analysis (Dai et al., 2017) of carbon assimilation based on the Farquhar et al. (1980) photosynthesis model, and its variant proposed by Collatz et al. (1991). Simulation 2, labelled "FvC_4proc_salt," is a 14 parameter sensitivity analysis (Saltelli et al., 2010) of carbon assimilation based on the Farquhar and Collatz photosynthesis models. Simulation 3, labelled "FvC," is a factorial combination of the alternative process representations and parameter values presented in those original papers across a factorial combination of incident photosynthetically active radiation (50-2000 µmol m-2 s-1, in increments of 50) and atmospheric CO2 concentration (50-1500 µmol m-2 s-1, in increments of 50). Additional scripts/functions for post-processing and plotting are part of the MAAT repository. These data were generated to evaluate and compare the Farquhar and Collatz models of leaf-scale photosynthesis across a range of CO2 and incident light conditions. Analysis of results can be found in Walker et al. (2020). This dataset is comprised of R scripts to process data and produce figures; .RDS files, a binary format specific to R and read with the 'readRDS' function, to store raw MAAT output; R and bash scripts (.bs) that help to run a MAAT simulation, .XML files that are output from MAAT describing a run, and .csv files which contain the calculated sensitivity indexes published in Walker et al. (2020).
The parameterization of key photosynthesis parameters is one of the key uncertain sources in modeling ecosystem gross primary productivity (GPP). Solar-induced chlorophyll fluorescence (SIF) offers a good proxy for GPP since it marks the actual process of photosynthesis; while machine learning (ML) provides a robust approach to model the GPP-SIF relationship. Here, we trained the boosted regressing tree (BRT) and the Random Forest ML models with Greenhouse Gases Observing Satellite SIF data and in situ GPP observations from 49 eddy covariance towers. These trained ML GPP-SIF models were fed into the Energy Exascale Earth System Model (E3SM) Land Model (ELM) to generate ELM-simulated global SIF estimates, which were then benchmarked against satellite SIF observations with a surrogate modeling approach. Our results indicated good modeling performance of the ML-based GPP-SIF relationship. The ELM model when fed with the ML GPP-SIF models also can well predict the spatial-temporal variations in SIF. We also found high model accuracy for the surrogate modeling. Model parameter sensitivity analysis suggested that the fraction of leaf nitrogen in RuBisCO (flnr) is the most sensitive parameter to the SIF; other sensitive parameters include the Ball-Berry stomatal conductance slope (mbbopt) and the vcmax entropy (vcmaxse). The posterior uncertainty in simulated GPP was greatly reduced after benchmarking, and the model produced improved spatial patterns of mean GPP relative to FLUXCOM GPP. Our integrated approach provides a new avenue for improving land models and using remote-sensing SIF, which can be further improved in the future with more ground- and satellite-based observations.
Abstract. Major multi-reservoir cascades represent a primary mechanism for dealing with hydrologic variability and extremes within institutionally complex river basins worldwide. These coordinated management processes fundamentally reshape water balance dynamics. Yet, multi-reservoir coordination processes have been largely ignored in the increasingly sophisticated representations of reservoir operations within large-scale hydrological models. The aim of this paper is twofold, namely (i) to provide evidence that the common modeling practice of parameterizing each reservoir in a cascade independently from the others is a significant approximation and (ii) to demonstrate potential unintended consequences of this independence approximation when simulating the dynamics of hydrological extremes in complex reservoir cascades. We explore these questions using the Water Balance Model, which features detailed representations of the human infrastructure coupled to the natural processes that shape water balance dynamics. It is applied to the Upper Snake River basin in the western US and its heavily regulated multi-reservoir cascade. We employ a time-varying sensitivity analysis that utilizes the method of Morris factor screening to explicitly track how the dominant release rule parameters evolve both along the cascade and in time according to seasonal high- and low-flow events. This enables us to address aim (i) by demonstrating how the progressive and cumulative dominance of upstream releases significantly dampens the ability of downstream reservoir rules' parameters to influence flow conditions. We address aim (ii) by comparing simulation results with observed reservoir operations during critical low-flow and high-flow events in the basin. Our time-varying parameter sensitivity analysis with the method of Morris clarifies how independent single-reservoir parameterizations and their tacit assumption of independence leads to reservoir release behaviors that generate artificial water shortages and flooding, whereas the observed coordinated cascade operations avoided these outcomes for the same events. To further explore the role of (non-)coordination in the large deviations from the observed operations, we use an offline multi-reservoir water balance model in which adding basic coordination mechanisms drawn from the observed emergency operations is sufficient to correct the deficiencies of the independently parameterized reservoir rules from the hydrological model. These results demonstrate the importance of understanding the state–space context in which reservoir releases occur and where operational coordination plays a crucial role in avoiding or mitigating water-related extremes. Understanding how major infrastructure is coordinated and controlled in major river basins is essential for properly assessing future flood and drought hazards in a changing world.
Abstract. Vegetation plays a key role in the global carbon cycle and thus is an important component within Earth system models (ESMs) that project future climate. Many ESMs are adopting methods to resolve plant size and ecosystem disturbance history, using vegetation demographic models. These models make it feasible to conduct more realistic simulation of processes that control vegetation dynamics. Meanwhile, increasing understanding of the processes governing plant water use, and ecosystem responses to drought in particular, has led to the adoption of dynamic plant water transport (i.e., hydrodynamic) schemes within ESMs. However, the extent to which variations in plant hydraulic traits affect both plant water stress and the risk of mortality in trait-diverse tropical forests is understudied. In this study, we report on a sensitivity analysis of an existing hydrodynamic scheme (HYDRO) model that is updated and incorporated into the Functionally Assembled Terrestrial Ecosystem Simulator (FATES) (FATES–HYDRO V1.0). The size- and canopy-structured representation within FATES is able to simulate how plant size and hydraulic traits affect vegetation dynamics and carbon–water fluxes. To better understand this new model system, and its functionality in tropical forest systems in particular, we conducted a global parameter sensitivity analysis at Barro Colorado Island, Panama. We assembled 942 observations of plant hydraulic traits on 306 tropical plant species for stomata, leaves, stems, and roots and determined the best-fit statistical distribution for each trait, which was used in model parameter sampling to assess the parametric sensitivity. We showed that, for simulated leaf water potential and loss of hydraulic conductivity across different plant organs, the four most important traits were associated with xylem conduit taper (buffers increasing hydraulic resistance with tree height), stomatal sensitivity to leaf water potential, maximum stem hydraulic conductivity, and the partitioning of total hydraulic resistance above vs. belowground. Our analysis of individual ensemble members revealed that trees at a high risk of hydraulic failure and potential tree mortality generally have a lower conduit taper, lower maximum xylem conductivity, lower stomatal sensitivity to leaf water potential, and lower resistance to xylem embolism for stem and transporting roots. We expect that our results will provide guidance on future modeling studies using plant hydrodynamic models to predict the forest responses to droughts and future field campaigns that aim to better parameterize plant hydrodynamic models.
Abstract. Organic matter (OM) composition plays a central role in microbial respiration of dissolved organic matter and subsequent biogeochemical reactions. Here, a direct connection of organic matter chemistry and thermodynamics to reactive transport simulators has been achieved through the newly developed Lambda-PFLOTRAN workflow tool that succinctly incorporates carbon chemistry data generated from Fourier transform ion cyclotron resonance mass spectrometry (FTICR-MS) into reaction networks to simulate organic matter degradation and the resulting biogeochemistry. Lambda-PFLOTRAN is a Python-based workflow, executed through a Jupyter notebook interface, that digests raw FTICR-MS data, develops a representative reaction network based on substrate-explicit thermodynamic modeling (also termed lambda modeling due to its key thermodynamic parameter λ used therein), and completes a biogeochemical simulation with the open source, reactive flow and transport code PFLOTRAN. The workflow consists of the following five steps: configuration, thermodynamic (lambda) analysis, sensitivity analysis, parameter estimation, and simulation output and visualization. Two test cases are provided to demonstrate the functionality of the Lambda-PFLOTRAN workflow. The first test case uses laboratory incubation data of temporal oxygen depletion to fit lambda parameters (i.e., maximum utilization rate and microbial carrying capacity). A slightly more complex second test case fits multiple lambda formulation and soil organic matter release parameters to temporal greenhouse gas generation measured during a soil incubation. Overall, the Lambda-PFLOTRAN workflow facilitates upscaling by using molecular-scale characterization to inform biogeochemical processes occurring at larger scales.
Many advanced nuclear reactor concepts currently being developed are targeting higher operating temperatures relative to the current fleet of light water nuclear reactors, for efficiency gains and other operational considerations. The design of high temperature structural components with reliable long-term operational performance will depend on material models that accurately capture the inelastic deformation mechanisms active in these environments. In this work, we perform a detailed parameter sensitivity analysis of two unified elasto-viscoplastic Grade 91 material models capable of capturing long term high temperature creep deformation. The first model is a phenomelogical material model from the Nuclear Engineering Material Library (NEML) developed at Argonne National Lab. The NEML model parameters and their uncertainty were fit to a range of Grade 91 experimental data using Bayesian Markov Chain Monte Carlo analysis. The second model is a LAROMance data-driven surrogate material model developed at Los Alamos National Lab. The LAROMance model is fit to a large database of responses produced by a mechanistic crystal plasticity based polycrystal model. Parameters for the LAROMance surrogate material model reflect the pedigree of the Grade 91 microstructure. Both material models have been integrated into the Grizzly code, based on the open-source MOOSE multiphysics simulation framework, to simulate both the progression of aging mechanisms and the effects of that aging on nuclear power plant structures. Grizzly is used analyze a three-dimensional Grade 91 piping system to compare the long-term inelastic response predicted by these two fundamentally different models and assess the sensitivity of the material model input parameters on this quantity of interest.