Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Stochastic”

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 811 records · Page 45

Multi-Period Active Distribution Network Planning Using Multi-Stage Stochastic Programming and Nested Decomposition by SDDIP

This paper presents a multi-period active distribution network planning (ADNP) with distributed generation (DG). The objective of the proposed ADNP is to minimize the total planning cost, subject to both investment and operation constraints. The paper proposes a multi-stage stochastic optimization model to address DG uncertainties over several periods, in which the decisions are made sequentially by only using the present-stage information. A nested decomposition method is proposed which applies the stochastic dual dynamic integer programming (SDDIP) method to address computational intractabilities of the proposed ADNP approach. The presented numerical results and discussions on a 33-bus distribution system and a large-scale 906-bus system verify the effectiveness of the proposed ADNP method and its solution method.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Sparse Control Synthesis for Uncertain Responsive Loads With Stochastic Stability Guarantees

In this report, recent studies have demonstrated the potential of flexible loads in providing frequency response services, predominantly due to their availability and cost-effectiveness. However, uncertainty and variability in various weather-related and end-use behavioral factors often impact the reliability of demand-side control performance. This work addresses this problem with the design of a demand-side control to achieve frequency response under load uncertainties. Our approach involves modeling the load uncertainties via stochastic processes that appear as both multiplicative and additive in the power system dynamics. Recently developed mean square exponential stability (MSES) results for continuous-time linear stochastic systems are applied to pose the control synthesis problem which results in an LMI-based optimization problem. Additional costs and constraints are added to the LMI-based controller synthesis to ensure MSES, improve closed-loop transient performance, maximize tolerable uncertainties, and promote sparsity in the controller. Additionally, the fundamental limitations between the tolerable uncertainties and control efforts while ensuring MSES are discussed. Further, the control synthesis problem for the case of the full-state measurement is generalized to the case of partial-state measurements. The proposed control synthesis is illustrated on an IEEE 39 bus system with rigorous studies to demonstrate the role of sparsity, closed-loop transient performance, tolerable uncertainties, and control efforts while ensuring MSES and achieving frequency response.

42 ENGINEERING↗

A Stochastic Multi-Commodity Logistic Model for Disaster Preparation in Distribution Systems

In this paper, we propose a stochastic optimization approach for disaster preparation in distribution systems. For an upcoming storm, utilities should have a preparation plan that includes warehousing restoration supplies, securing staging sites (depots), and prepositioning crews and equipment. Pre-storm planning enables faster and more efficient post-disaster deploy ment of crews and equipment resources to damage locations. To assist utilities in making this important preparation, this paper develops a two-stage stochastic mixed integer linear program. The first stage determines the depots, number of crews in each site, and the amount of equipment. The second stage is the recourse action that deals with acquiring new equipment and assigning crews to repair damages in realized scenarios. The objective of the developed model is to minimize the costs of depots, crews, equipment, and penalty costs associated with delays in obtain ing equipment and restoration. We consider the uncertainties of damaged lines, number and type of equipment required, and expected repair times. The model is validated on modified IEEE 123-bus distribution test system.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Reserve Model of Energy Storage in Day-ahead Joint Energy and Reserve Markets: A Stochastic UC Solution

With many favorable advantages including fast response ability in particular, utility-level energy storage systems (ESS) are being integrated into energy and reserve markets to help mitigate uncertain renewable resources and fluctuant demands. This paper discusses a stochastic unit commitment (UC) model to explore capabilities of ESSs in providing valuable grid services by simultaneously joining energy and reserve markets. The proposed reserve model of ESSs presents the following features: (i) two constraints are proposed to formulate ESS’s reserve provision ability in each hour via six operation modes, namely increasing/reducing the level of charge, switching to discharge, increasing/reducing the level of discharge, switching to charge; (ii) as an energy-limited asset, constraints on ESS’s reserve deliverability across multiple hours is introduced to hold enough state of charge (SOC) headroom and floor room, guaranteeing that reserves cleared in individual hours are continuously deliverable over multiple successive hours without violating operation limits. Furthermore, a scenario-based UC model, which addresses uncertainties of renewables and demands, is adopted to evaluate reserves provided by ESSs and generators. To effectively solve the stochastic UC problem, the progressive hedging algorithm with heuristic approaches is discussed. A 6-bus system and a modified IEEE 118-bus system are used to illustrate effectiveness of the proposed approaches.

25 ENERGY STORAGE↗

A Stochastic Multi-Criteria Decision-Making Algorithm for Dynamic Load Prioritization in Grid-Interactive Efficient Buildings

Increasing deployment of advanced sensing, controls, and communication infrastructure enables buildings to provide services to the power grid, leading to the concept of grid-interactive efficient buildings. Since occupant activities and preferences primarily drive the availability and operational flexibility of building devices, there is a critical need to develop occupant-centric approaches that prioritize devices for providing grid services, while maintaining the desired end-use quality of service. In this paper, we present a decision-making framework that facilitates a building owner/operator to effectively prioritize loads for curtailment service under uncertainties, while minimizing any adverse impact on the occupants. The proposed framework uses a stochastic (Markov) model to represent the probabilistic behavior of device usage from power consumption data, and a load prioritization algorithm that dynamically ranks building loads using a stochastic multi-criteria decision-making algorithm. The proposed load prioritization framework is illustrated via numerical simulations in a residential building use-case, including plug-loads, air-conditioners, and plug-in electric vehicle chargers, in the context of load curtailment as a grid service. Suitable metrics are proposed to evaluate the closed-loop performance of the proposed prioritization algorithm under various scenarios and design choices. Scalability of the proposed algorithm is established via computational analysis, while time-series plots are used for intuitive explanation of the ranking choices.

24 POWER TRANSMISSION AND DISTRIBUTION↗

A stochastic model of cell survival for high-Z nanoparticle radiotherapy

Purpose: The authors present a stochastic framework for the assessment of cell survival in gold nanoparticle radiotherapy. Methods: The authors derive the equations for the effective macroscopic dose enhancement for a population of cells with nonideal distribution of gold nanoparticles (GNP), allowing different number of GNP per cell and different distances with respect to the cellular target. They use the mixed Poisson distribution formalism to model the impact of the aforementioned physical factors on the effective dose enhancement. Results: The authors show relatively large differences in the estimation of cell survival arising from using approximated formulae. They predict degeneration of the cell killing capacity due to different number of GNP per cell and different distances with respect to the cellular target. Conclusions: The presented stochastic framework can be used in interpretation of experimental cell survival or tumor control probability studies.

60 APPLIED LIFE SCIENCES↗

MO-DE-BRA-03: TOPAS-edu: A Window Into the Stochastic World Through the TOPAS Tool for Particle Simulation

Purpose: The stochastic nature of the subatomic world presents a challenge for physics education. Even experienced physicists can be amazed at the varied behavior of electrons, x-rays, protons, neutrons, ions and the any short-lived particles that make up the overall behavior of our accelerators, brachytherapy sources and medical imaging systems. The all-particle Monte Carlo particle transport tool, TOPAS Tool for Particle Simulation, originally developed for proton therapy research, has been repurposed into a physics teaching tool, TOPAS-edu. Methods: TOPAS-edu students set up simulated particle sources, collimators, scatterers, imagers and scoring setups by writing simple ASCII files (in the TOPAS Parameter Control System format). Students visualize geometry setups and particle trajectories in a variety of modes from OpenGL graphics to VRML 3D viewers to gif and PostScript image files. Results written to simple comma separated values files are imported by the student into their preferred data analysis tool. Students can vary random seeds or adjust parameters of physics processes to better understand the stochastic nature of subatomic physics. Results: TOPAS-edu has been successfully deployed as the centerpiece of a physics course for master’s students at Queen’s University Belfast. Tutorials developed there takes students through a step by step course on the basics of particle transport and interaction, scattering, Bremsstrahlung, etc. At each step in the course, students build simulated experimental setups and then analyze the simulated results. Lessons build one upon another so that a student might end up with a full simulation of a medical accelerator, a water-phantom or an imager. Conclusion: TOPAS-edu was well received by students. A second application of TOPAS-edu is currently in development at Zurich University of Applied Sciences, Switzerland. It is our eventual goal to make TOPAS-edu available free of charge to any non-profit organization, along with associated tutorial materials developed by the TOPAS-edu community. Work supported in part by the U.S. Department of Energy under contract number DE-AC02-76SF00515. B. Villagomez-Bernabe is supported by CONACyT (Mexican Council for Science and Technology) project 231844.

60 APPLIED LIFE SCIENCES↗

Propagation of excitation in long 1D chains: Transition from regular quantum dynamics to stochastic dynamics

The quantum dynamics problem for a 1D chain consisting of 2N + 1 sites (N Much-Greater-Than 1) with the interaction of nearest neighbors and an impurity site at the middle differing in energy and in coupling constant from the sites of the remaining chain is solved analytically. The initial excitation of the impurity is accompanied by the propagation of excitation over the chain sites and with the emergence of Loschmidt echo (partial restoration of the impurity site population) in the recurrence cycles with a period proportional to N. The echo consists of the main (most intense) component modulated by damped oscillations. The intensity of oscillations increases with increasing cycle number and matrix element C of the interaction of the impurity site n = 0 with sites n = {+-}1 (0 < C {<=} 1; for the remaining neighboring sites, the matrix element is equal to unity). Mixing of the components of echo from neighboring cycles induces a transition from the regular to stochastic evolution. In the regular evolution region, the wave packet propagates over the chain at a nearly constant group velocity, embracing a number of sites varying periodically with time. In the stochastic regime, the excitation is distributed over a number of sites close to 2N, with the populations varying irregularly with time. The model explains qualitatively the experimental data on ballistic propagation of the vibrational energy in linear chains of CH{sub 2} fragments and predicts the possibility of a nondissipative energy transfer between reaction centers associated with such chains.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Crossing resonance of stochastically interacting wave fields

The dynamic susceptibilities (Green's functions) of the system of two interacting wave fields of different physical natures with a stochastically inhomogeneous coupling parameter between them with zero mean value have been examined. The well-known self-consistent approximation taking into account all diagrams with noncrossing correlation/interaction lines has been generalized to the case of stochastically interacting wave fields. The analysis has been performed for spin and elastic waves. The results obtained taking into account the processes of multiple scattering of waves from inhomogeneities are significantly different from those obtained for this situation earlier in the Bourret approximation [R.C. Bourret, Nuovo Cimento 26, 1 (1962)]. Instead of frequencies degeneracy removal in the wave spectrum and the splitting of resonance peaks of dynamic susceptibilities, a wide single-mode resonance peak should be observed at the crossing point of the unperturbed dispersion curves. The fine structure appears at vertices of these wide peaks in the form of a narrow resonance on the Green's-function curve of one field and a narrow antiresonance on the vertex of the Green's-function curve of the other field.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Stochastic variational quantization and maximum entropy principle

In this note, we discuss the stochastic variational method (SVM) to derive the action principle of Schrödinger equation. In this formalism, the wave function is introduced as one of convenient ways to express the action. This work is a reformulation of Nelson—Yasue’s stochastic quantization scheme and its variational formulation by Yasue, and will furnish possible way to deal with more complicated system such as quantum field theory, suggesting an origin of quantum mechanics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

An Asynchronous Bundle-Trust-Region Method for Dual Decomposition of Stochastic Mixed-Integer Programming

Here, we present an asynchronous bundle-trust-region algorithm within the context of Lagrangian dual decomposition for stochastic mixed-integer programs. The approach solves the Lagrangian master problem by using a bundle method with a trust-region constraint. This scheme enables asynchronous computations and can thus help mitigate severe load imbalance issues (associated with the solution of scenario subproblems) and improve parallel efficiency. We provide a convergence analysis and an implementation of the proposed scheme. We also present extensive numerical results on eighty instances of a large-scale stochastic unit commitment problem, and demonstrate that the proposed approach provides significant reductions in solution time and achieves strong scaling.

97 MATHEMATICS AND COMPUTING↗

Adaptive Uncertainty Quantification for Stochastic Hyperbolic Conservation Laws

Here, we propose a predictor-corrector adaptive method for the study of hyperbolic partial differential equations (PDEs) under uncertainty. Constructed around the framework of stochastic finite volume (SFV) methods, our approach circumvents sampling schemes or simulation ensembles while also preserving fundamental properties, in particular hyperbolicity of the resulting systems and conservation of the discrete solutions. Furthermore, we augment the existing SFV theory with a priori convergence results for statistical quantities, in particular push-forward densities, which we demonstrate through numerical experiments. By linking refinement indicators to regions of the physical and stochastic spaces, we drive anisotropic refinements of the discretizations, introducing new degrees of freedom where deemed profitable. To illustrate our proposed method, we consider a series of numerical examples for nonlinear hyperbolic PDEs based on Burgers’ and Euler’s equations.

97 MATHEMATICS AND COMPUTING↗

The Optical Signatures of Stochastic Processes in Many-Body Exciton Scattering

We review our recent quantum stochastic model for spectroscopic lineshapes in the presence of a coevolving and nonstationary background population of excitations. Starting from a field theory description for interacting bosonic excitons, we derive a reduced model whereby optical excitons are coupled to an incoherent background via scattering as mediated by their screened Coulomb coupling. The Heisenberg equations of motion for the optical excitons are then driven by an auxiliary stochastic population variable, which we take to be the solution of an Ornstein–Uhlenbeck process. Here, we present an overview of the theoretical techniques we have developed as applied to predicting coherent nonlinear spectroscopic signals. We show how direct (Coulomb) and exchange coupling to the bath give rise to distinct spectral signatures and discuss mathematical limits on inverting spectral signatures to extract the background density of states.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A stochastic electricity market clearing formulation with consistent pricing properties

We argue that deterministic market clearing formulations introduce arbitrary distortions between day-ahead and expected real-time prices that bias economic incentives. We extend and analyze a previously proposed stochastic clearing formulation in which the social surplus function induces penalties between day-ahead and real-time quantities. We prove that the formulation yields price bounded price distortions, and we show that adding a similar penalty term to transmission flows and phase angles ensures boundedness throughout the network. We prove that when the price distortions are zero, day-ahead quantities equal a quantile of their real-time counterparts. The undesired effects of price distortions suggest that stochastic settings provide significant benefits over deterministic ones that go beyond social surplus improvements. Finally, we propose additional metrics to evaluate these benefits.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Stochastic Optimization of Complex Systems

This project focused on methodologies for the solution of stochastic optimization problems based on relaxation and penalty methods, Monte Carlo simulation, parallel processing, and inverse optimization. The main results of the project were the development of a convergent method for the solution of models that include expectation constraints as in equilibrium models, improvement of Monte Carlo convergence through the use of a new method of sample batch optimization, the development of new parallel processing methods for stochastic unit commitment models, and the development of improved methods in combination with parallel processing for incorporating automatic differentiation methods into optimization.

97 MATHEMATICS AND COMPUTING↗

NSR&D Program Fiscal Year 2015 Funded Research Stochastic Modeling of Radioactive Material Releases Final Report

Nonreactor nuclear facilities operating under the approval authority of the U.S. Department of Energy use unmitigated hazard evaluations to determine if potential radiological doses associated with design basis events challenge or exceed dose evaluation guidelines. Unmitigated design basis events that sufficiently challenge dose evaluation guidelines or exceed the guidelines for members of the public or workers, merit selection of safety structures, systems, or components or other controls to prevent or mitigate the hazard. Idaho State University, in collaboration with Idaho National Laboratory, has developed a portable and simple to use software application called SODA (Stochastic Objective Decision-Aide) that stochastically calculates the radiation dose distribution associated with hypothetical radiological material release scenarios. Rather than producing a point estimate of the dose, SODA produces a dose distribution result to allow a deeper understanding of the dose potential. SODA allows users to select the distribution type and parameter values for all of the input variables used to perform the dose calculation. Users can also specify custom distributions through a user defined distribution option. SODA then randomly samples each distribution input variable and calculates the overall resulting dose distribution. In cases where an input variable distribution is unknown, a traditional single point value can be used. SODA, developed using the MATLAB coding framework, has a graphical user interface and can be installed on both Windows and Mac computers. SODA is a standalone software application and does not require MATLAB to function. SODA provides improved risk understanding leading to better informed decision making associated with establishing nuclear facility material-at-risk limits and safety structure, system, or component selection. It is important to note that SODA does not replace or compete with codes such as MACCS or RSAC; rather it is viewed as an easy to use supplemental tool to help improve risk understanding and support better informed decisions. The SODA development project was funded through a grant from the DOE Nuclear Safety Research and Development Program.

42 ENGINEERING↗

R&D toward a proof-of-principle experiment on optical stochastic cooling at the Fermilab’s IOTA ring

The stochastic-cooling method was invented in the late 70s and played a key role in the accumulation of intense beams of antiprotons for a head-on collision with counter rotating proton beams at 500 GeV in the Super Proton Synchrotron (SPS) at CERN. Such an improvement led to the discovery of the W and Z bosons in 1983. In brief, stochastic cooling relies on recording information (position) associated with short longitudinal slices within the bunch. The radiation is encoded on a radiofrequency (RF) signal, manipulated and used to correct the bunch slices at a later time.The repetitive correction is applied as the beam circulates in the ring ultimately cool the beam.Stochastic cooling has been successfully implemented in several machines for particle cooling and accumulation. Cooling of bunched beams in proton-(anti)proton colliders was challenging due to the high phase-space density of the bunches but was recently demonstrated at the Relativistic Heavy Ion Collider (RHIC).

43 PARTICLE ACCELERATORS↗

Multi-Fidelity Stochastic Economic Dispatch for Operating Low Carbon Power Grids

Grid operators can address the inherently stochastic nature of renewables by solving a two-stage stochastic programming model that minimizes the cost of dispatch decisions while accounting for the complex grid dynamics. It is common to use a sample average approximation to estimate the expectation of the second stage costs in this model. However, the large sample count needed for numerical accuracy makes effective modeling large-scale electric grids computationally intractable. We introduce a control variate multi-fidelity estimator for the second-stage recourse that enables high quality dispatch decisions in real-time with a reduced computational burden. We obtain a hierarchy of model fidelities by linearizing the AC power flow system representation to DC power flow, and by relaxing transmission and voltage network constraints. We evaluate the performance of our proposed method on a synthetic grid with 73 buses against a deterministic baseline with persistence forecast and a high-fidelity reference. Our analysis shows a computational speed-up of 7.62x with a minimal loss in accuracy. The multi-fidelity method is well suited to fidelity combinations that use a simplified network topology in their lower fidelity model and is an attractive option for applications where accurate grid modeling needed on a limited computational budget.

24 POWER TRANSMISSION AND DISTRIBUTION↗