Search NASA⌕ Search

SEARCH · Search NASA

Results for “conditional probability distribution”

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 109 records · Page 6

Sensitivity of Probable Maximum Flood in a Changing Environment

With likely increases in probable maximum precipitation (PMP) in a changing environment, critical infrastructures such as major reservoirs and nuclear power plants are subject to elevated risk. To understand how factors such as PMP variability, climate change, land use land cover (LULC) change, antecedent soil moisture conditions, and reservoir storage may individually or jointly affect the magnitude of probable maximum flood (PMF), we conducted integrated hydrometeorological simulations involving both the Weather Research Forecasting model and the distributed hydrologic model (DHSVM) over the Alabama–Coosa–Tallapoosa (ACT) River Basin in the southeastern United States. A total of 120 relative humidity–maximized PMP storms under historic and projected future climate conditions were used to drive DHSVM in current and projected future LULC conditions. Overall, PMP and PMF are projected to increase significantly over the ACT River Basin. Sources of meteorological forcing data sets and climate change were found to be the most sensitive factors affecting PMF, followed by antecedent soil moisture, reservoir storage, and then LULC change. In conclusion, the ensemble of PMP and PMF simulations, along with their sensitivity, allows us to better quantify the potential risks associated with hydroclimatic extreme events to critical infrastructures for energy–water security.

54 ENVIRONMENTAL SCIENCES↗

Minimizing Uncertainties Impact in Decision Making with an Applicability Study for Economic Power Dispatch

This report summaries the work performed under the LDRD project on the preliminary study on knowledge automation, where specific focus has been made on the investigation of the impact of uncertainties of human decision making onto the optimization of the process operation. At first the statistics on signals from the Brain-Computing Interface (BCI) is analyzed so as to obtain the uncertainties characterization of human operators during the decision making phase using the electroencephalogram (EEG) signals. This is then followed by the discussions of an architecture that reveals the equivalence between optimization and closed loop feedback control design, where it has been shown that all the optimization problems can be transferred into the control design problem for closed loop systems. This has led to a “closed loop” framework, where the structure of the decision making is shown to be subjected to both process disturbances and controller’s uncertainties. The latter can well represent the uncertainties or randomness occurred during human decision making phase. As a result, a stochastic optimization problem has been formulated and a novel solution has been proposed using probability density function (PDF) shaping for both the cost function and the constraints using stochastic distribution control concept. A sufficient condition has been derived that guarantees the convergence of the optimal solution and discussions have been made for both the total probabilistic solution and chanced constrained optimization which have been well-studied in optimal power flows (OPF) area. A simple case study has been carried out for the economic dispatch of powers for a grid system when there are distributed energy resources (DERs) in the system, and encouraging results have been obtained showing that a significant savings on the generation cost can be expected.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗

A probabilistic creep model incorporating test condition, initial damage, and material property uncertainty

Uncertainty is prevalent in the creep resistance of alloys, where at elevated temperature and low pressure, rupture can range across logarithmic decades. In this study, a probabilistic continuum-damage-mechanics (CDM)-based model is derived to capture the uncertainty of creep resistance. To meet this objective, creep data for alloy 304 Stainless Steel is gathered. A constitutive model, “Sinh”, is calibrated deterministically to determine the statistical variability of the material properties. Three sources of uncertainty are injected into the model: test condition (stress and temperature), initial damage, and material properties. Probabilistic simulations are carried out by (a) calibrating probability distribution functions (pdfs) for each source of uncertainty (b) randomly sampling the pdfs using Monte Carlo methods and (c) executing simulations to replicate the uncertain creep behavior. A sensitivity analysis is performed to evaluate the relative effect of each source of uncertainty. In full probabilistic simulations, the cumulative uncertainty of creep behavior is evaluated. The probabilistic model accurately predicts the creep deformation and rupture of the available experiments. The probabilistic model is validated for interpolation but lacks extrapolation ability. Several future works are proposed to further improve the model.

36 MATERIALS SCIENCE↗

Environmental Selection, Dispersal, and Organism Interactions Shape Community Assembly in High-Throughput Enrichment Culturing

A central goal of microbial ecology is to identify and quantify the forces that lead to observed population distributions and dynamics. However, these forces, which include environmental selection, dispersal, and organism interactions, are often difficult to assess in natural environments. Here in this paper, we present a method that links microbial community structures with selective and stochastic forces through highly replicated subsampling and enrichment of a single environmental inoculum. Specifically, groundwater from a well-studied natural aquifer was serially diluted and inoculated into nearly 1,000 aerobic and anaerobic nitrate-reducing cultures, and the final community structures were evaluated with 16S rRNA gene amplicon sequencing. We analyzed the frequency and abundance of individual operational taxonomic units (OTUs) to understand how probabilistic immigration, relative fitness differences, environmental factors, and organismal interactions contributed to divergent distributions of community structures. We further used a most probable number (MPN) method to estimate the natural condition-dependent cultivable abundance of each of the nearly 400 OTU cultivated in our study and infer the relative fitness of each. Additionally, we infer condition-specific organism interactions and discuss how this high-replicate culturing approach is essential in dissecting the interplay between overlapping ecological forces and taxon-specific attributes that underpin microbial community assembly. IMPORTANCE Through highly replicated culturing, in which inocula are subsampled from a single environmental sample, we empirically determine how selective forces, interspecific interactions, relative fitness, and probabilistic dispersal shape bacterial communities. These methods offer a novel approach to untangle not only interspecific interactions but also taxon-specific fitness differences that manifest across different cultivation conditions and lead to the selection and enrichment of specific organisms. Additionally, we provide a method for estimating the number of cultivable units of each OTU in the original sample through the MPN approach.

16S RNA↗

Review of the WECC EDT phase 2 EIM benefits analysis and results report.

A region-wide Energy Imbalance Market (EIM) was recently proposed by the Western Electricity Coordinating Council (WECC). In order for the Western Area Power Administration (Western) to make more informed decisions regarding its involvement in the EIM, Western asked Argonne National Laboratory (Argonne) to review the EIM benefits study (the October 2011 revision) performed by Energy and Environmental Economics, Inc. (E3). Key components of the E3 analysis made use of results from a study conducted by the National Renewable Energy Laboratory (NREL); therefore, we also reviewed the NREL work. This report examines E3 and NREL methods and models used in the EIM study. Estimating EIM benefits is very challenging because of the complex nature of the Western Interconnection (WI), the variability and uncertainty of renewable energy resources, and the complex decisions and potentially strategic bidding of market participants. Furthermore, methodologies used for some of the more challenging aspects of the EIM have not yet matured. This review is complimentary of several components of the EIM study. Analysts and modelers clearly took great care when conducting detailed simulations of the WI using well-established industry tools under stringent time and budget constraints. However, it is our opinion that the following aspects of the study and the interpretation of model results could be improved upon in future analyses. The hurdle rate methodology used to estimate current market inefficiencies does not directly model the underlying causes of sub-optimal dispatch and power flows. It assumes that differences between historical flows and modeled flows can be attributed solely to market inefficiencies. However, flow differences between model results and historical data can be attributed to numerous simplifying assumptions used in the model and in the input data. We suggest that alternative approaches be explored in order to better estimate the benefits of introducing market structures like the EIM. In addition to more efficient energy transactions in the WI, the EIM would reduce the amount of flexibility reserves needed to accommodate forecast errors associated with variable production from wind and solar energy resources. The modeling approach takes full advantage of variable resource diversity over the entire market footprint, but the projected reduction in flexibility reserves may be overly optimistic. While some reduction would undoubtedly occur, the EIM is only an energy market and would therefore not realize the same reduction in reserves as an ancillary services market. In our opinion the methodology does not adequately capture the impact of transmission constraints on the deployment of flexibility reserves. Estimates of flexibility reserves assume that forecast errors follow a normal distribution. Improved estimates could be obtained by using other probability distributions to estimate up and down reserves to capture the underlying uncertainty of these resources under specific operating conditions. Also, the use of a persistence forecast method for solar is questionable, because solar insolation follows a deterministic pattern dictated by the sun's path through the sky. We suggest a more rigorous method for forecasting solar insolation using the sun's relatively predictable daily pattern at specific locations. The EIM study considered only one scenario for hydropower resources. While this scenario is within the normal range over the WI footprint, it represents a severe drought condition in the Colorado River Basin from which Western schedules power. Given hydropower's prominent role in the WI, we recommend simulating a range of hydropower conditions since the relationship between water availability and WI dispatch costs is nonlinear. Also, the representation of specific operational constraints faced by hydropower operators in the WI needs improvements. The model used in the study cannot fully capture all of the EIM impacts and complexities of power system operations. In particular, a primary benefit of the EIM is a shorter dispatch interval; namely, 5 minutes. However, the model simulates the dispatch hourly. Therefore it cannot adequately measure the benefits of a more frequent dispatch. A tool with a finer time resolution would significantly improve simulation accuracy. When the study was conducted, the rules for the EIM were not clearly defined and it was appropriate to estimate societal benefits of the EIM assuming a perfect market without a detailed specification of the market design. However, incorporating a more complete description of market rules will allow for better estimates of EIM benefits. Furthermore, performing analyses using specific market rules can identify potential design flaws that may be difficult and expensive to correct after the market is established. Estimated cost savings from a more efficient dispatch are less than one percent of the total cost of electricity production.

14 SOLAR ENERGY↗

The use of copulas to practical estimation of multivariate stochastic differential equation mixed effects models

A system of stochastic differential equations (SDE) with mixed-effects parameters and multivariate normal copula density function were used to develop tree height model for Scots pine trees in Lithuania. A two-step maximum likelihood parameter estimation method is used and computational guidelines are given. After fitting the conditional probability density functions to outside bark diameter at breast height, and total tree height, a bivariate normal copula distribution model was constructed. Predictions from the mixed-effects parameters SDE tree height model calculated during this research were compared to the regression tree height equations. The results are implemented in the symbolic computational language MAPLE.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Estimating the influence of field inventory sampling intensity on forest landscape model performance for determining high-severity wildfire risk

Abstract Historically, fire has been essential in Southwestern US forests. However, a century of fire-exclusion and changing climate created forests which are more susceptible to uncharacteristically severe wildfires. Forest managers use a combination of thinning and prescribed burning to reduce forest density to help mitigate the risk of high-severity fires. These treatments are laborious and expensive, therefore optimizing their impact is crucial. Landscape simulation models can be useful in identifying high risk areas and assessing treatment effects, but uncertainties in these models can limit their utility in decision making. In this study we examined underlying uncertainties in the initial vegetation layer by leveraging a previous study from the Santa Fe fireshed and using new inventory plots from 111 stands to interpolate the initial forest conditions. We found that more inventory plots resulted in a different geographic distribution and wider range of the modelled biomass. This changed the location of areas with high probability of high-severity fires, shifting the optimal location for management. The increased range of biomass variability from using a larger number of plots to interpolate the initial vegetation layer also influenced ecosystem carbon dynamics, resulting in simulated forest conditions that had higher rates of carbon uptake. We conclude that the initial forest layer significantly affects fire and carbon dynamics and is dependent on both number of plots, and sufficient representation of the range of forest types and biomass density.

Science & Technology - Other Topics↗

Delineating Facies Spatial Distribution by Integrating Ensemble Data Assimilation and Indicator Geostatistics with Level Set Transformation.

A new approach is developed to delineate the spatial distribution of discrete facies (geological units that have unique distributions of hydraulic, physical, and/or chemical properties) conditioned not only on direct data (measurements directly related to facies properties, e.g., grain size distribution obtained from borehole samples) but also on indirect data (observations indirectly related to facies distribution, e.g., hydraulic head and tracer concentration). Our method integrates for the first time ensemble data assimilation with traditional transition probability-based geostatistics. The concept of level set is introduced to build shape parameterization that allows transformation between discrete facies indicators and continuous random variables. The spatial structure of different facies is simulated by indicator models using conditioning points selected adaptively during the iterative process of data assimilation. To evaluate the new method, a two-dimensional semi-synthetic example is designed to estimate the spatial distribution and permeability of two distinct facies from transient head data induced by pumping tests. The example demonstrates that our new method adequately captures the spatial pattern of facies distribution by imposing spatial continuity through conditioning points. The new method also reproduces the overall response in hydraulic head field with better accuracy compared to data assimilation with no constraints on spatial continuity on facies.

97 MATHEMATICS AND COMPUTING↗

Agent-based simulation and child protection systems: Rationale, implementation, and verification

Simulation models are an important tool used in health care and other disciplines to support operational research and decision-making. In the child protection literature, simulation models are an under-utilized source of research evidence. Here, in this paper, we describe the rationale for and the development of an agent-based simulation of a child protection system in the US. Using the investigation, prevention service, and placement histories of 600,000 children served in an urban child welfare system, we walk the reader through the development of a prototype known as OSPEDALE. The governing equations built into OSPEDALE probabilistically simulate the onset of investigations. Then, drawing from empirical survival distributions, the governing equations trace the probability of subsequent interactions with the system (recurrence of maltreatment, service referrals, and placement) conditional on the characteristics of children, their assessed risk level, and prior child protection system involvement. As an initial test of OSPEDALE's utility, we compare empirical admission counts with counts generated from OSPEDALE. Though the verification step is admittedly simple, the comparison shows that OSPEDALE replicates the empirical count of new admissions closely enough to justify further investment in OSPEDALE. Management of public child protection systems is increasingly research evidence-dependent. The emphasis on research evidence as a decision-support tool has elevated evidence acquired through randomized clinical trials. Though important, the evidence from clinical trials represents only one type of research evidence. Properly specified, simulation models are another source of evidence with real-world relevance.

60 APPLIED LIFE SCIENCES↗

Cellulose nanofibrils and nanocrystals in confined flow: Single-particle dynamics to collective alignment revealed through scanning small-angle x-ray scattering and numerical simulations

Nanostructured materials made through flow-assisted assembly of proteinaceous or polymeric nanosized fibrillar building blocks are promising contenders for a family of high-performance biocompatible materials in a wide variety of applications. Optimization of these processes relies on improving our knowledge of the physical mechanisms from nano- to macroscale and especially understanding the alignment of elongated nanoparticles in flows. Here, we study the full projected orientation distributions of cellulose nanocrystals (CNCs) and nanofibrils (CNFs) in confined flow using scanning microbeam SAXS. For CNCs, we further compare with a simulated system of dilute Brownian ellipsoids, which agrees well at dilute concentrations. However, increasing CNC concentration to a semidilute regime results in locally arranged domains called tactoids, which aid in aligning the CNC at low shear rates, but limit alignment at higher rates. Similarly, shear alignment of CNF at semidilute conditions is also limited owing to probable bundle or flock formation of the highly entangled nanofibrils. Finally, this work provides a quantitative comparison of full projected orientation distributions of elongated nanoparticles in confined flow and provides an important stepping stone towards predicting and controlling processes to create nanostructured materials on an industrial scale.

36 MATERIALS SCIENCE↗

A Resilient Condition Assessment Monitoring System

An architecture and supporting methods are presented for the implementation of a resilient condition assessment monitoring system that can adaptively accommodate both cyber and physical anomalies to a monitored system under observation. In particular, the architecture includes three layers: information, assessment, and sensor selection. The information layer estimates probability distributions of process variables based on sensor measurements and assessments of the quality of sensor data. Based on these estimates, the assessment layer then employs probabilistic reasoning methods to assess the plant health. The sensor selection layer selects sensors so that assessments of the plant condition can be made within desired time periods. Resilient features of the developed system are then illustrated by simulations of a simplified power plant model, where a large portion of the sensors are under attack.

97 MATHEMATICS AND COMPUTING↗

A computationally efficient ductile damage model accounting for nucleation and micro-inertia at high triaxialities

The computational formulation of a micro-mechanical material model for the dynamic failure of ductile metals is presented in this paper. The statistical nature of porosity initiation is accounted for by introducing an arbitrary probability density function which describes the pores nucleation pressures. Each micropore within the representative volume element is modeled as a thick spherical shell made of plastically incompressible material. The treatment of porosity by a distribution of thick-walled spheres also allows for the inclusion of micro-inertia effects under conditions of shock and dynamic loading. The second order ordinary differential equation governing the microscopic porosity evolution is solved with a robust implicit procedure. A new Chebyshev collocation method is employed to approximate the porosity distribution and remapping is used to optimize memory usage. The adaptive approximation of the porosity distribution leads to a reduction of computational time and memory usage of up to two orders of magnitude. Moreover, the proposed model affords consistent performance: changing the nucleation pressure probability density function and/or the applied strain rate does not reduce accuracy or computational efficiency of the material model. The numerical performance of the model and algorithms presented is tested against three problems for high density tantalum: single void, one-dimensional uniaxial strain, and two-dimensional plate impact. Here, the results using the integration and algorithmic advances suggest a significant improvement in computational efficiency and accuracy over previous treatments for dynamic loading conditions.

42 ENGINEERING↗

Distributionally Robust Building Load Control to Compensate Fluctuations in Solar Power Generation

This paper investigates the use of a collection of dispatchable heating, ventilation and air conditioning (HVAC) systems to absorb low-frequency fluctuations in renewable energy sources, especially in solar photo-voltaic (PV) generation. Given the uncertain and time-varying nature of solar PV generation, its probability distribution is difficult to be estimated perfectly, which poses a challenging problem of how to optimally schedule a fleet of HVAC loads to consume as much as local PV generation. We formulate a distributionally robust chance-constrained (DRCC) model to ensure that PV generation is consumed with a desired probability for a family of probability distributions, termed as an ambiguity set, built upon mean and covariance information. We benchmark the DRCC model with a deterministic optimization model and a stochastic programming model in a one-day simulation. We show that the DRCC model achieves constantly good performance to consume most PV generation even in the case with the presence of probability distribution ambiguity.

Zhang, Yiling↗

The creation of a comprehensive metocean data set for offshore wind turbine simulations: Comprehensive metocean data set

A database of meteorological and ocean conditions is presented for use in offshore wind energy research and design. The original data are from 23 ocean sites around the USA and were obtained from the National Data Buoy Center run by the National Oceanic and Atmospheric Administration. The data are presented in a processed form that includes the variables of interest for offshore wind energy design: wind speed, significant wave height, wave peak-spectral period, wind direction and wave direction. For each site, a binning process is conducted to create conditional probability functions for each of these variables. The sites are then grouped according to geographic location and combined to create three representative sites, including a West Coast site, an East Coast site and a Gulf of Mexico site. Both the processed data and the probability distribution parameters for the individual and representative sites are being hosted on a publicly available domain by the National Renewable Energy Laboratory, with the intent of providing a standard basis of comparison for meteorological and ocean conditions for offshore wind energy research worldwide.

17 WIND ENERGY↗

Uncertainty quantification tools for multiphase gas-solid flow simulations using MFIX

Computational fluid dynamics (CFD) has been widely studied and used in the scientific community and in the industry. Various models were proposed to solve problems in different areas. However, all models deviate from reality. Uncertainty quantification (UQ) process evaluates the overall uncertainties associated with the prediction of quantities of interest. In particular it studies the propagation of input uncertainties to the outputs of the models so that confidence intervals can be provided for the simulation results. In the present work, a non-intrusive quadrature-based uncertainty quantification (QBUQ) approach is proposed. The probability distribution function (PDF) of the system response can be then reconstructed using extended quadrature method of moments (EQMOM) and extended conditional quadrature method of moments (ECQMOM). The report first explains the theory of QBUQ approach, including methods to generate samples for problems with single or multiple uncertain input parameters, low order statistics, and required number of samples. Then methods for univariate PDF reconstruction (EQMOM) and multivariate PDF reconstruction (ECQMOM) are explained. The implementation of QBUQ approach into the open-source CFD code MFIX is discussed next. At last, QBUQ approach is demonstrated in several applications. The method is first applied to two examples: a developing flow in a channel with uncertain viscosity, and an oblique shock problem with uncertain upstream Mach number. The error in the prediction of the moment response is studied as a function of the number of samples, and the accuracy of the moments required to reconstruct the PDF of the system response is discussed. The QBUQ approach is then demonstrated by considering a bubbling fluidized bed as example application. The mean particle size is assumed to be the uncertain input parameter. The system is simulated with a standard two-fluid model with kinetic theory closures for the particulate phase implemented into MFIX. The effect of uncertainty on the disperse-phase volume fraction, on the phase velocities and on the pressure drop inside the fluidized bed are examined, and the reconstructed PDFs are provided for the three quantities studied. Then the approach is applied to a bubbling fluidized bed with two uncertain parameters, particle-particle and particle-wall restitution coefficients. Contour plots of the mean and standard deviation of solid volume fraction, solid phase velocities and gas pressure are provided. The PDFs of the response are reconstructed using EQMOM with appropriate kernel density functions. The simulation results are compared to experimental data provided by the 2013 NETL small-scale challenge problem. Lastly, the proposed procedure is demonstrated by considering a riser of a circulating fluidized bed as an example application. The mean particle size is considered to be the uncertain input parameter. Contour plots of the mean and standard deviation of solid volume fraction, solid phase velocities, and granular temperature are provided. Mean values and confidence intervals of the quantities of interest are compared to the experiment results. The univariate and bivariate PDF reconstructions of the system response are performed using EQMOM and ECQMOM.

97 MATHEMATICS AND COMPUTING↗

Index for Wind Power Variability

Variability of large scale wind power generation is dependent on several factors: characteristics of installed wind power plants, size of the area where the plants are installed, geographic dispersion within that area and its weather regime(s). Variability can be described by ramps in power generation, i.e. changes from time period to time period. Given enough data points, it can be described with a probability density function. This approach focuses on two dimensions of variability: duration of the ramp and probability distribution. This paper proposes an index based on these two dimensions to enable comparisons and characterizations of variability under different conditions. The index is tested with real, large scale wind power generation data from several countries. Considerations while forming an index are discussed, as well as the main results regarding what the drivers of variability experienced for different data.

17 WIND ENERGY↗

Learning the temporal evolution of multivariate densities via normalizing flows

In this work, we propose a method to learn multivariate probability distributions using sample path data from stochastic differential equations. Specifically, we consider temporally evolving probability distributions (e.g., those produced by integrating local or nonlocal Fokker–Planck equations). Here, we analyze this evolution through machine learning assisted construction of a time-dependent mapping that takes a reference distribution (say, a Gaussian) to each and every instance of our evolving distribution. If the reference distribution is the initial condition of a Fokker–Planck equation, what we learn is the time-T map of the corresponding solution. Specifically, the learned map is a multivariate normalizing flow that deforms the support of the reference density to the support of each and every density snapshot in time. We demonstrate that this approach can approximate probability density function evolutions in time from observed sampled data for systems driven by both Brownian and Lévy noise. We present examples with two- and three-dimensional, uni- and multimodal distributions to validate the method.

97 MATHEMATICS AND COMPUTING↗

Microfield dynamics in dense hydrogen plasmas with high- Z impurities

We use large-scale classical molecular dynamics to determine microfield properties for several dense plasma mixtures. By employing quantum statistical potentials (QSPs) to regularize the Coulomb interaction, our simulations follow motions of electrons as well as ions for times long enough to track relaxation phenomena involving both types of particles. Coulomb coupling, relative to temperature, of different pairs of species in the hot, dense matter being simulated ranges from weak to strong. We first study the effect of such coupling differences, along with composition and QSP differences, on the roles of electrons and various mixture components in determining probability distributions of instantaneous, total microfields experienced by the ions. Then, we address two important dynamical questions: (1) How is the quasistatic part of the total field to be extracted from the time-dependent simulation data? (2) Under what conditions does the commonly used approximation of ions with fixed Yukawa-like screening by free electrons accurately describe quasistatic fields? We identify a running, short-time average of the total field at each ion as its slowly evolving, quasistatic part. We consider several ways to specify the averaging interval, and note the influence of ion dynamics in this issue. When all species are weakly coupled, the quasistatic fields have probability distributions agreeing well with those we obtain from simulations of Yukawa-screened ions. However, agreement deteriorates as the coupling between high- Z ions increases well beyond unity, principally because the Yukawa model tends to underestimate the true screening of close high- Z pairs. Lastly, examples of this fact are given, and some consequences for the high-field portions of probability distributions are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗