Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “stochastic analysis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9

The structure of dilute combusting sprays

An experimental and theoretical study of drop processes in a turbulent flame is described. The experiments involved a monodisperse (105 and 180 micro m initial diameter) stream of methanol drops injected at the base of a turbulent methane-fueled diffusion flame burning in still air. The following measurements were made: mean and fluctuating phase velocities, mean drop number flux, drop-size distributions and mean gas-phase temperatures. Measurements were compared with predictions of two separated flow models: (1) deterministic separated flow, where drop-turbulence interactions are ignored; and (2) stochastic separated flow, where drop-turbulence interactions are considered using random-walk computations. The stochastic separated flow analysis yielded best agreement with measurements, since it provides for turbulent dispersion of drops which was important for present test conditions (and probably for most combusting sprays as well). Distinguishing the presence or absence of envelope flames around the drops, however, was relatively unimportant for present test conditions, since the drops spent most of their lifetime in fuel-rich regions of the flow where this distinction is irrelevant.

Shuen, J. S.↗

Generating Sequential PV Deployment Scenarios for High Renewable Distribution Grid Planning: Preprint

This paper introduces a novel approach for generating solar photovoltaic (PV) plant deployment scenarios for grid integration planning. The approach guarantees consistency among scenarios of the same deployment by ensuring that higher penetration scenarios contain PV units deployed in lower penetration scenarios. It also constrains the size and spatial distribution of the PV plants and considers three placement types. A case study on a real-world distribution system proves that the precepts of scenario consistency, deployment diversity, and placement are met. The study further investigates the impact of the resulting scenarios via a stochastic hosting capacity analysis. Results indicate that the ratio between PV and load sizes, referred to as the nodal PV penetration factor (NPPF), is a key driver of the grid integration impact. By reducing the NPPF from 5 to 2, the maximum hosting capacity increased by at least 112%. The study also reveals that scenarios under random placement can lead to higher hosting capacity values.

ENERGY PLANNING, POLICY, AND ECONOMY,SOLAR ENERGY↗

A modular methodology for time-domain stochastic seismic wave propagation

Presented here is a modular methodology for time-domain stochastic seismic wave propagation analysis. Presented methodology is designed to analyse uncertain seismic motions as an input, propagating through uncertain material. Traditional approach for uncertain wave propagation relies on models that include deep bedrock, local soil site, and their random process and random field information. Such models can become quite large and computationally intractable. The modular approach proposed herein features two step approach that allows separate consideration of the deep bedrock and local site along with corresponding random field information. In this work, the first step considers an auxiliary stochastic motions problem in the bedrock. Stochastic local site response can then be simulated in a reduced domain within certain depth from the surface. Application of uncertain seismic motions at depth, for local uncertain site response is done using stochastic effective forces developed through the Domain Reduction Method. By using Hermite polynomial chaos expansion to represent the non-Gaussian random field of material parameters and non-stationary random process of seismic motion, the proposed modular methodology is formulated using intrusive stochastic Galerkin approach, as seen in the Stochastic Elastic–Plastic Finite Element Method (SEPFEM). Developed modular methodology is illustrated using a 1-D stochastic seismic wave propagation analysis with three cases, and simulation results are also verified with results from conventional approach.

58 GEOSCIENCES↗

Flexural rigidity of the Basin and Range-Colorado Plateau-Rocky Mountain transition from coherence analysis of gravity and topography

Stochastic inversion for flexural loads and flexural rigidity of the continental elastic layer can be accomplished most effectively by using the coherence of gravity and topography. However, the spatial resolution of coherence analysis has been limited by use of two-dimensional periodogram spectra from very large (greater than 10(exp 5)sq km) windows that generally include multiple tectonic features. Using a two-dimensional spectral estimator based on the maximum entropy method, the spatial resolution of flexural proerties can be enhanced by a factor of 4 or more, enabling more detailed analysis at the scale of individual tectonic features. This new approach is used to map the spatial variation of flexural rigidity along the Basin and Range transition to the Colorado Plateau and Middle Rocky Mountains physiographic provinces. Large variations in flexural isostatic responses are found, with rigidities ranging from as low as 8.7 x 10(exp 20) N m (elastic thickness (T(sub e) = 4.6 km) in the Basin and Range to as high as 4.1 x 10(exp 24) N m T(sub e) = 77 km) in the Middle Rocky Mountains. These results compare favorably woith independent determinations of flexural rigidity in the region. Areas of low flexural rigidity correlate strongly with areas of high surface heat flow, as is expected from the contingence of flexural rigidity on a temperature-dependent flow law. Also, late Cenozoic normal faults with large displacements are found primarily in area of low flexural rigidity region. The highest flexural rigidity is found within the Archean Wyoming craton, where evidence suggests that deeply rooted cratonic lithosphere may play a role in determining the distribution of tectonism at the surface.

Lowery, Anthony R.↗

Particle-laden swirling free jets: Measurements and predictions

A theoretical and experimental investigation of single-phase and particle-laden weakly swirling jets was conducted. The jets were injected vertically downward from a 19 mm diameter tube with swirl numbers ranging from 0 to 0.33. The particle-laden jets had a single loading ratio (0.2) with particles having a SMD of 39 microns. Mean and fluctuating properties of both phases were measured using nonintrusive laser based methods while particle mass flux was measured using an isokinetic sampling probe. The continuous phase was analyzed using both a baseline kappa-epsilon turbulence model and an extended version with modifications based on the flux Richardson number to account for effects of streamline curvature. To highlight effects of interphase transport rates and particle/turbulence interactions, effects of the particles were analyzed as follows: (1) locally homogeneous flow (LHF) analysis, where interphase transport rates are assumed to be infinitely fast; (2) deterministic separated flow (DSF) analysis, where finite interphase transport rates are considered but particle/turbulence interactions are ignored; and (3) stochastic separated flow (SSF) analysis, where both effects are considered using random-walk computations.

Bulzan, D. L.↗

Particle-laden swirling free jets - Measurements and predictions

A theoretical and experimental investigation of single-phase and particle-laden weakly swirling jets was conducted. The jets were injected vertically downward from a 19 mm diameter tube with swirl numbers ranging from 0 to 0.33. The particle-laden jets had a single loading ratio (0.2) with particles having an SMD of 39 micrometers. Mean and fluctuating properties of both phases were measured using nonintrusive laser based methods while particle mass flux was measured using an isokinetic sampling probe. The continuous phase was analyzed using both a baseline k-epsilon turbulence model and an extended version with modifications based on the flux Richardson number to account for effects of streamline curvature. To highlight effects of interphase transport rates and particle/turbulence interactions, effects of the particles were analyzed as follows: (1) locally homogeneous flow (LHF) analysis, where interphase transport rates are assumed to be infinitely fast; (2) deterministic separated flow (DSF) analysis, where finite interphase transport rates are considered but particle/turbulence interactions are ignored; and (3) stochastic separated flow (SSF) analysis, where both effects are considered using random-walk computations.

Bulzan, D. L.↗

Structure of nonevaporating sprays. II - Drop and turbulence properties

This is the second part of a study reporting structure measurements in the dilute portion of axisymmetric nonevaporating sprays. Measurements are compared with predictions of three typical methods for analyzing sprays: (1) locally homogeneous flow (LHF) analysis, where slip between the phases is neglected; (2) deterministic separated flow (DSF) analysis, where slip is considered but effects of drop interactions with turbulence are ignored; and (3) stochastic separated flow (SSF) analysis, where both slip and effects of drop interactions with turbulence are considered. This part of the study reports measurements of mean and fluctuating drop velocities, the variation of Sauter mean diameter, and gas-phase turbulence properties in the dilute portion of the sprays. Best agreement between predictions and measurements was obtained with the SSF analysis. For present measurements in the dilute region (void fraction greater than 99.1 percent), effects of drops on gas-phase turbulence properties (turbulence modulation) were small. However, as the dense spray regions near the injector were approached, the measurements indicated modification of turbulence properties by drop motion.

Solomon, A. S. P.↗

A Novel Simulation-Based Framework for Sensor Error Impact Analysis in Smart Building Systems: A Case Study for a Demand-Controlled Ventilation System

Sensors are one of the fundamental components for sensor-rich controls in buildings but are prone to different errors. Existing studies show that sensor errors hold a place among top-priority faults in building systems. Before we take countermeasures to mitigate the sensor errors, it is vital to prioritize key sensors and quantify the collective impacts of concurrent sensor errors. In response to this, a simulation-based methodology is introduced to conduct a comprehensive sensor error impact analysis in building systems, which adds a stochastic sensor prioritization through a sensitivity analysis on top of a commonly used deterministic sensor error quantification. The synergies of these two parts help better interpret the sensor error impacts on building energy consumption, ventilation performance, thermal comfort, etc. A sensor-rich CO2-based Demand-Controlled Ventilation system is used as a case study to demonstrate the viability of the methodology as a proof-of-the-concept. The results show that the energy savings potential and ventilation performance are mostly influenced by the accuracy of the AHU outdoor airflow sensors. The accuracy of zone level airflow sensors has a negligible impact on both energy savings and ventilation performance. The accuracy of zone CO2 sensors has more influence on the ventilation performance compared with the accuracy of zone airflow sensors. Compared with the baseline case with zero errors, the largest deviation percentages could reach 16.90% and 94.32%, respectively, in terms of the Heating, Ventilation, and Air-Conditioning (HVAC) annual energy consumption and the Outdoor Air Ratio (OAR) when multiple key sensors suffer from normal error intensities simultaneously.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Risk-averse optimal control of semilinear elliptic PDEs

In this paper, we consider the optimal control of semilinear elliptic PDEs with random inputs. These problems are often nonconvex, infinite-dimensional stochastic optimization problems for which we employ risk measures to quantify the implicit uncertainty in the objective function. In contrast to previous works in uncertainty quantification and stochastic optimization, we provide a rigorous mathematical analysis demonstrating higher solution regularity (in stochastic state space), continuity and differentiability of the control-to-state map, and existence, regularity and continuity properties of the control-to-adjoint map. Our proofs make use of existing techniques from PDE-constrained optimization as well as concepts from the theory of measurable multifunctions. We illustrate our theoretical results with two numerical examples motivated by the optimal doping of semiconductor devices.

97 MATHEMATICS AND COMPUTING↗

Further studies using matched filter theory and stochastic simulation for gust loads prediction

This paper describes two analysis methods - one deterministic, the other stochastic - for computing maximized and time-correlated gust loads for aircraft with nonlinear control systems. The first method is based on matched filter theory; the second is based on stochastic simulation. The paper summarizes the methods, discusses the selection of gust intensity for each method and presents numerical results. A strong similarity between the results from the two methods is seen to exist for both linear and nonlinear configurations.

Scott, Robert C.↗

Further studies using matched filter theory and stochastic simulation for gust loads prediction

This paper describes two analysis methods -- one deterministic, the other stochastic -- for computing maximized and time-correlated gust loads for aircraft with nonlinear control systems. The first method is based on matched filter theory; the second is based on stochastic simulation. The paper summarizes the methods, discusses the selection of gust intensity for each method and presents numerical results. A strong similarity between the results from the two methods is seen to exist for both linear and nonlinear configurations.

Scott, Robert C.↗

Structure of nonevaporating sprays. I - Initial conditions and mean properties

Structure measurements were completed within the dilute portion of axisymmetric nonevaporating sprays. Measurements included: mean velocities, velocity fluctuations, and Reynolds stress of the gas phase; and mean velocities, fluctuating velocities, mass flux, and diameter distributions of the drop phase. The measurements were used to evaluate three typical methods of analyzing sprays: (1) a locally homogeneous flow (LHF) analysis, where slip between the phases is neglected; (2) a deterministic separated flow (DSF) analysis, where slip is considered but effects of drop interactions with turbulence are ignored; and (3) a stochastic separated flow (SSF) analysis, where effects of both slip and turbulence are considered using random-walk computations for drop motion. Measurements of initial conditions of both phases near the injector, mean gas-phase properties, and liquid flux distributions are described herein. Best agreement between predictions and measurements was obtained with the SSF model, which provided a reasonable representation of turbulent dispersion of drops. A companion paper presents additional measurements of drop continuous phase properties.

Solomon, A. S. P.↗

Methods for R&D Portfolio Analysis and Evaluation (Workshop Report)

The Workshop on Methods for R&D Portfolio Analysis and Evaluation convened on 17–18 July 2019 at the National Renewable Energy Laboratory in Golden, Colorado, and examined strengths and weaknesses of the various methodologies applicable to R&D portfolio modeling, analysis, and decision support, given pragmatic constraints such as data availability, uncertainties in estimating the impact of R&D spending, and practical operational overheads. Participants employed their deep expertise in approaches such as stochastic optimization, real options, Monte-Carlo analysis, Bayesian networks, decision theory, complex systems analysis, deep uncertainty, and technology-evolution modeling to critique the initial example models developed by the project’s core team and to conduct thought experiments grounded in real-life technology models, progress data, expert elicitation, and portfolio information. This engagement of participants’ methodological expertise with the practical requirements of real-life portfolio decision support yielded ideas for improved approaches, alternative methodological hypotheses, and hybridization of methodologies that are well-grounded theoretically, computationally sound, and realistically executable given data availability and other practical constraints.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Methods for R&D Portfolio Analysis and Evaluation (Workshop Report)

The Workshop on Methods for R&D Portfolio Analysis and Evaluation convened on 17-18 July 2019 at the National Renewable Energy Laboratory in Golden, Colorado, and examined strengths and weaknesses of the various methodologies applicable to R&D portfolio modeling, analysis, and decision support, given pragmatic constraints such as data availability, uncertainties in estimating the impact of R&D spending, and practical operational overheads. Participants employed their deep expertise in approaches such as stochastic optimization, real options, Monte-Carlo analysis, Bayesian networks, decision theory, complex systems analysis, deep uncertainty, and technology-evolution modeling to critique the initial example models developed by the project’s core team and to conduct thought experiments grounded in real-life technology models, progress data, expert elicitation, and portfolio information. This engagement of participants’ methodological expertise with the practical requirements of real-life portfolio decision support yielded ideas for improved approaches, alternative methodological hypotheses, and hybridization of methodologies that are well-grounded theoretically, computationally sound, and realistically executable given data availability and other practical constraints.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Dynamics retrieval from stochastically weighted incomplete data by low-pass spectral analysis

Time-resolved serial femtosecond crystallography (TR-SFX) provides access to protein dynamics on sub-picosecond timescales, and with atomic resolution. Due to the nature of the experiment, these datasets are often highly incomplete and the measured diffracted intensities are affected by partiality. To tackle these issues, one established procedure is that of splitting the data into time bins, and averaging the multiple measurements of equivalent reflections within each bin. This binning and averaging often involve a loss of information. Here, we propose an alternative approach, which we call low-pass spectral analysis (LPSA). In this method, the data are projected onto the subspace defined by a set of trigonometric functions, with frequencies up to a certain cutoff. This approach attenuates undesirable high-frequency features and facilitates retrieving the underlying dynamics. A time-lagged embedding step can be included prior to subspace projection to improve the stability of the results with respect to the parameters involved. Subsequent modal decomposition allows to produce a low-rank description of the system's evolution. Using a synthetic time-evolving model with incomplete and partial observations, we analyze the LPSA results in terms of quality of the retrieved signal, as a function of the parameters involved. We compare the performance of LPSA to that of a range of other sophisticated data analysis techniques. We show that LPSA allows to achieve excellent dynamics reconstruction at modest computational cost. Finally, we demonstrate the superiority of dynamics retrieval by LPSA compared to time binning and merging, which is, to date, the most commonly used method to extract dynamical information from TR-SFX data.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Analysis and improvement of the vector quantization in SELP (Stochastically Excited Linear Prediction)

The Stochastically Excited Linear Prediction (SELP) algorithm is described as a speech coding method employing a two-stage vector quantization. The first stage uses an adaptive codebook which efficiently encodes the periodicity of voiced speech, and the second stage uses a stochastic codebook to encode the remainder of the excitation signal. The adaptive codebook performs well when the pitch period of the speech signal is larger than the frame size. An extension is introduced, which increases its performance for the case that the frame size is longer than the pitch period. The performance of the stochastic stage, which improves with frame length, is shown to be best in those sections of the speech signal where a high level of short-term correlations is present. It can be concluded that the SELP algorithm performs best during voiced speech where the pitch period is longer than the frame length.

Kleijn, W. B.↗

Stochastic Gradients for Large-Scale Tensor Decomposition

Tensor decomposition is a well-known tool for multiway data analysis. This work proposes using stochastic gradients for efficient generalized canonical polyadic (GCP) tensor decomposition of large-scale tensors. GCP tensor decomposition is a recently proposed version of tensor decomposition that allows for a variety of loss functions such as Bernoulli loss for binary data or Huber loss for robust estimation. Here, the stochastic gradient is formed from randomly sampled elements of the tensor and is efficient because it can be computed using the sparse matricized-tensor times Khatri--Rao product tensor kernel. For dense tensors, we simply use uniform sampling. For sparse tensors, we propose two types of stratified sampling that give precedence to sampling nonzeros. Numerical results demonstrate the advantages of the proposed approach and its scalability to large-scale problems.

97 MATHEMATICS AND COMPUTING↗