Search NASA⌕ Search

SEARCH · Search NASA

Results for “Stochastic Modeling and Simulation”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

A Mechanistic Model for Cooperative Behavior of Co-transcribing RNA Polymerases

In fast-transcribing prokaryotic genes, such as an rrn gene in Escherichia coli, many RNA polymerases (RNAPs) transcribe the DNA simultaneously. Active elongation of RNAPs is often interrupted by pauses, which has been observed to cause RNAP traffic jams; yet some studies indicate that elongation seems to be faster in the presence of multiple RNAPs than elongation by a single RNAP. We propose that an interaction between RNAPs via the torque produced by RNAP motion on helically twisted DNA can explain this apparent paradox. We have incorporated the torque mechanism into a stochastic model and simulated transcription both with and without torque. Simulation results illustrate that the torque causes shorter pause durations and fewer collisions between polymerases. Our results suggest that the torsional interaction of RNAPs is an important mechanism in maintaining fast transcription times, and that transcription should be viewed as a cooperative group effort by multiple polymerases.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Bill Savings vs. Backup Power: Evaluating operational tradeoffs for home solar+storage systems [Slides]

This study explores tradeoffs between the use of home solar+storage systems for backup power versus day-to-day utility bill savings. The study focuses specifically on the “reserve setting” available with most home battery storage systems, which allow the customer to maintain some minimum level of storage in reserve in case of an unforeseen power interruption. The more capacity that is held in reserve, the greater the customer’s ability to ride-through possible power interruptions, but less capacity is then available to manage utility bills on a day-to-day basis. This study evaluates this operational tradeoff across a diverse set of locations and residential electricity tariff structures, relying on Berkeley Lab’s PRESTO model to stochastically simulate power interruption events, and exploring a range of sensitivities, including variations in customer value of lost load (VoLL), interruption frequency, and other key drivers. The results show that, in most circumstances, the opportunity cost of holding storage capacity in reserve, in terms of foregone bill saving, tends to outweigh any gains in reliability value associated with mitigated power interruptions. This finding is robust across tariff structures and across most of the sensitivities considered, including those related to rate level, customer load level, and storage sizing. There are a limited set of circumstances where raising the reserve setting improves the overall customer value (comprised of bill savings plus reliability value). Specifically, that exception occurs when all of the following conditions apply: (a) the customer resides in a location with exceptionally poor reliability, (b) the customer has exceptionally high VoLL; (c) the customer is on a net billing rate or on a TOU rate that allows grid discharging but not grid charging; and (d), depending on the location, the price arbitrage differential on that rate is relatively small. In all other circumstances analyzed, total customer value declines with reserve level.

14 SOLAR ENERGY↗
DOE OSTI2014

Synthetic Biology Outside the Cell: Linking Computational Tools to Cell-Free Systems

As mathematical models become more commonly integrated into the study of biology, a common language for describing biological processes is manifesting. Many tools have emerged for the simulation of in vivo synthetic biological systems, with only a few examples of prominent work done on predicting the dynamics of cell-free synthetic systems. At the same time, experimental biologists have begun to study dynamics of in vitro systems encapsulated by amphiphilic molecules, opening the door for the development of a new generation of biomimetic systems. In this review, we explore both in vivo and in vitro models of biochemical networks with a special focus on tools that could be applied to the construction of cell-free expression systems. We believe that quantitative studies of complex cellular mechanisms and pathways in synthetic systems can yield important insights into what makes cells different from conventional chemical systems.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1986

Communications satellite business ventures - Measuring the impact of technology programmes and related policies

An economic evaluation and planning procedure which assesses the effects of various policies on fixed satellite business ventures is described. The procedure is based on a stochastic financial simulation model, the Domsat II, which evaluates spacecraft reliability, market performance, and cost uncertainties. The application of the Domsat II model to the assessment of NASA's ion thrusters for on-orbit propulsion and GaAs solar cell technology is discussed. The effects of insurance rates and the self-insurance option on the financial performance of communication satellite business ventures are investigated. The selection of a transportation system for placing the satellites into GEO is analyzed.

Greenberg, J. S.↗
DOE OSTI2017

An agent-based stochastic Occupancy Simulator

Occupancy has significant impacts on building performance. However, in current building performance simulation programs, occupancy inputs are static and lack diversity, contributing to discrepancies between the simulated and actual building performance. This work presents an Occupancy Simulator that simulates the stochastic behavior of occupant presence and movement in buildings, capturing the spatial and temporal occupancy diversity. Each occupant and each space in the building are explicitly simulated as an agent with their profiles of stochastic behaviors. The occupancy behaviors are represented with three types of models: (1) the status transition events (e.g., first arrival in office) simulated with probability distribution model, (2) the random moving events (e.g., from one office to another) simulated with a homogeneous Markov chain model, and (3) the meeting events simulated with a new stochastic model. A hierarchical data model was developed for the Occupancy Simulator, which reduces the amount of data input by using the concepts of occupant types and space types. Finally, a case study of a small office building is presented to demonstrate the use of the Simulator to generate detailed annual sub-hourly occupant schedules for individual spaces and the whole building. The Simulator is a web application freely available to the public and capable of performing a detailed stochastic simulation of occupant presence and movement in buildings. Future work includes enhancements in the meeting event model, consideration of personal absent days, verification and validation of the simulated occupancy results, and expansion for use with residential buildings.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

Extreme-scale stochastic optimization and simulation via learning-enhanced decomposition and parallelization (Final Technical Report)

Stochastic optimization and simulation models ubiquitously arise in designing and operating complex service/engineering systems. They can be extreme in scale due to high-dimensional data and decisions, and can also involve decisions made sequentially in response to newly revealed data, both causing significant computational challenge. The objective of this research is to explore a unified framework that integrates machine learning with discrete optimization and risk-averse modeling, to improve the efficiency of decomposition paradigms for stochastic optimization and simulations at extreme scale. The models we consider represent a broad class of complex decision-making problems, where 0-1 or continuous decisions are made before and/or after knowing multiple sources of uncertainties that could be correlated. We will employ machine learning methods to dynamically decide and prioritize computational procedures, including cut generation, branching, and bounding of the optimal objective. Furthermore, the research will shed new lights on the traditional decomposition algorithms for extreme-scale computing. Deliverables of the research include new modeling and computational methods for advancing the state-of-the-art research in optimization and simulation, bringing many relevant risk-averse, data-driven optimization problems in practice within the range of tractability. Examples include distributed computing server scheduling and sensor deployment for monitoring critical infrastructures. Success in this effort will enable progress in solving multiple extreme-scale problems in the complex system design and operations arising from DoE missions in energy, environment, and national security.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2014

Estimating the fraction of progeny virions that must incorporate APOBEC3G for suppression of productive HIV-1 infection

The contest between the host factor APOBEC3G (A3G) and the HIV-1 protein Vif presents an attractive target of intervention. The extent to which the A3G–Vif interaction must be suppressed to tilt the balance in favor of A3G remains unknown. We employed stochastic simulations and mathematical modeling of the within-host dynamics and evolution of HIV-1 to estimate the fraction of progeny virions that must incorporate A3G to render productive infection unsustainable. Using three different approaches, we found consistently that a transition from sustained infection to suppression of productive infection occurred when the latter fraction exceeded ∼0.8. The transition was triggered by A3G-induced hypermutations that led to premature stop codons compromising viral production and was consistent with driving the basic reproductive number, R{sub 0}, below unity. The fraction identified may serve as a quantitative guideline for strategies targeting the A3G–Vif axis. - Highlights: • We perform simulations and mathematical modeling of the role of APOBEC3G in suppressing HIV-1 infection. • In three distinct ways, we estimate that when over 80% of progeny virions carry APOBEC3G, productive HIV-1 infection would be suppressed. • Our estimate of this critical fraction presents quantitative guidelines for strategies targeting the APOBEC3G–Vif axis.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Accelerating population balance-Monte Carlo simulation for coagulation dynamics from the Markov jump model, stochastic algorithm and GPU parallel computing

This paper proposes a comprehensive framework for accelerating population balance-Monte Carlo (PBMC) simulation of particle coagulation dynamics. By combining Markov jump model, weighted majorant kernel and GPU (graphics processing unit) parallel computing, a significant gain in computational efficiency is achieved. The Markov jump model constructs a coagulation-rule matrix of differentially-weighted simulation particles, so as to capture the time evolution of particle size distribution with low statistical noise over the full size range and as far as possible to reduce the number of time loopings. Here three coagulation rules are highlighted and it is found that constructing appropriate coagulation rule provides a route to attain the compromise between accuracy and cost of PBMC methods. Further, in order to avoid double looping over all simulation particles when considering the two-particle events (typically, particle coagulation), the weighted majorant kernel is introduced to estimate the maximum coagulation rates being used for acceptance–rejection processes by single-looping over all particles, and meanwhile the mean time-step of coagulation event is estimated by summing the coagulation kernels of rejected and accepted particle pairs. The computational load of these fast differentially-weighted PBMC simulations (based on the Markov jump model) is reduced greatly to be proportional to the number of simulation particles in a zero-dimensional system (single cell). Finally, for a spatially inhomogeneous multi-dimensional (multi-cell) simulation, the proposed fast PBMC is performed in each cell, and multiple cells are parallel processed by multi-cores on a GPU that can implement the massively threaded data-parallel tasks to obtain remarkable speedup ratio (comparing with CPU computation, the speedup ratio of GPU parallel computing is as high as 200 in a case of 100 cells with 10 000 simulation particles per cell). These accelerating approaches of PBMC are demonstrated in a physically realistic Brownian coagulation case. The computational accuracy is validated with benchmark solution of discrete-sectional method. The simulation results show that the comprehensive approach can attain very favorable improvement in cost without sacrificing computational accuracy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Comparison of the kinetics of different Markov models for ligand binding under varying conditions

We recently derived a Markov model for macromolecular ligand binding dynamics from few physical assumptions and showed that its stationary distribution is the grand canonical ensemble [J. W. R. Martini, M. Habeck, and M. Schlather, J. Math. Chem. 52, 665 (2014)]. The transition probabilities of the proposed Markov process define a particular Glauber dynamics and have some similarity to the Metropolis-Hastings algorithm. Here, we illustrate that this model is the stochastic analog of (pseudo) rate equations and the corresponding system of differential equations. Moreover, it can be viewed as a limiting case of general stochastic simulations of chemical kinetics. Thus, the model links stochastic and deterministic approaches as well as kinetics and equilibrium described by the grand canonical ensemble. We demonstrate that the family of transition matrices of our model, parameterized by temperature and ligand activity, generates ligand binding kinetics that respond to changes in these parameters in a qualitatively similar way as experimentally observed kinetics. In contrast, neither the Metropolis-Hastings algorithm nor the Glauber heat bath reflects changes in the external conditions correctly. Both converge rapidly to the stationary distribution, which is advantageous when the major interest is in the equilibrium state, but fail to describe the kinetics of ligand binding realistically. To simulate cellular processes that involve the reversible stochastic binding of multiple factors, our pseudo rate equation model should therefore be preferred to the Metropolis-Hastings algorithm and the Glauber heat bath, if the stationary distribution is not of only interest.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Simulating the optical properties of soot using a stochastic soot model

KL extinction measurements are common in literature concerning soot producing sprays because of their nonintrusiveness. Unfortunately, these measurements often rely on bold assumptions of uniform and monodisperse spherical particles. In this work, a Spray A case is simulated using a highly detailed 3D stochastic soot model and the measured and simulated soot volume fraction are compared. Although the simulated volume fraction magnitude is approximately correct, features are much sharper than in the measured data. Optical models based on particle shape are applied to the simulated soot particles to nd the local optical thickness. The modeled optical thickness aligns with the measured volume fraction much more closely than the simulated volume fraction indicating KL extinction measurements are not simply related to soot volume fraction.

Strickland, Tyler↗
DOE OSTI2024

Computational Modeling, Simulation, and Potential Applications of Optical Stochastic Cooling

With the rising demand for intense particles beams, much research is being conducted in the area of particle beam cooling. One of these methods, called Stochastic Cooling (SC) (developed at CERN in the 1970's), delivered a feedback method to improve the quality and lifetime of circulating proton beams by reducing their 6D phase-space and has been widely implemented in a number of hadron machines. However, traditional stochastic cooling schemes are limited by the bandwidth of microwave frequency systems. Optical Stochastic Cooling (OSC) is a promising extension of the stochastic cooling beam cooling technique. OSC instead uses optical wavelengths which allows for improved control and increased cooling but creates its own technical challenges. This lays out work conducted toward the experimental demonstration of OSC at Fermilab's Integrable Optics Test Accelerator (IOTA) storage ring. This includes the design and characterization of parts of the optical delay system, the development and validation of a high-fidelity computational model of the OSC process, and the investigations into possible applications of the OSC mechanism to advanced beam manipulations.

43 PARTICLE ACCELERATORS↗
NASA NTRS1990

Comparison of effects of copropagated and precomputed atmosphere profiles on Monte Carlo trajectory simulation

A realization of a stochastic atmosphere model for use in simulations is presented. The model provides pressure, density, temperature, and wind velocity as a function of latitude, longitude, and altitude, and is implemented in a three degree of freedom simulation package. This implementation is used in the Monte Carlo simulation of an aeroassisted orbital transfer maneuver and results are compared to those of a more traditional approach.

Queen, Eric M.↗
NASA NTRS2003

Subgrid Combustion Modeling for the Next Generation National Combustion Code

In the first year of this research, a subgrid turbulent mixing and combustion methodology developed earlier at Georgia Tech has been provided to researchers at NASA/GRC for incorporation into the next generation National Combustion Code (called NCCLES hereafter). A key feature of this approach is that scalar mixing and combustion processes are simulated within the LES grid using a stochastic 1D model. The subgrid simulation approach recovers locally molecular diffusion and reaction kinetics exactly without requiring closure and thus, provides an attractive feature to simulate complex, highly turbulent reacting flows of interest. Data acquisition algorithms and statistical analysis strategies and routines to analyze NCCLES results have also been provided to NASA/GRC. The overall goal of this research is to systematically develop and implement LES capability into the current NCC. For this purpose, issues regarding initialization and running LES are also addressed in the collaborative effort. In parallel to this technology transfer effort (that is continuously on going), research has also been underway at Georgia Tech to enhance the LES capability to tackle more complex flows. In particular, subgrid scalar mixing and combustion method has been evaluated in three distinctly different flow field in order to demonstrate its generality: (a) Flame-Turbulence Interactions using premixed combustion, (b) Spatially evolving supersonic mixing layers, and (c) Temporal single and two-phase mixing layers. The configurations chosen are such that they can be implemented in NCCLES and used to evaluate the ability of the new code. Future development and validation will be in spray combustion in gas turbine engine and supersonic scalar mixing.

Menon, Suresh↗
DOE OSTI2017

Toward Development of a Stochastic Wake Model: Validation Using LES and Turbine Loads

Wind turbines within an array do not experience free-stream undisturbed flow fields. Rather, the flow fields on internal turbines are influenced by wakes generated by upwind unit and exhibit different dynamic characteristics relative to the free stream. The International Electrotechnical Commission (IEC) standard 61400-1 for the design of wind turbines only considers a deterministic wake model for the design of a wind plant. This study is focused on the development of a stochastic model for waked wind fields. First, high-fidelity physics-based waked wind velocity fields are generated using Large-Eddy Simulation (LES). Stochastic characteristics of these LES waked wind velocity field, including mean and turbulence components, are analyzed. Wake-related mean and turbulence field-related parameters are then estimated for use with a stochastic model, using Multivariate Multiple Linear Regression (MMLR) with the LES data. To validate the simulated wind fields based on the stochastic model, wind turbine tower and blade loads are generated using aeroelastic simulation for utility-scale wind turbine models and compared with those based directly on the LES inflow. The study’s overall objective is to offer efficient and validated stochastic approaches that are computationally tractable for assessing the performance and loads of turbines operating in wakes.

17 WIND ENERGY↗
DOE OSTI2023

Stochastic Models, Indices & Optimization Algorithms for Pricing & Hedging Reliability Risks in Modern Power Grids: Data Plan - Princeton

We collected and cleaned the synthetic grid data produced by NREL for the Texas and New York synthetic grids. We developed a high dimensional joint stochastic model for load at the zone level, and solar and wind power productions at the asset level, capturing the spatial and temporal dependencies between all the variables, and demonstrated how such a model could be fitted to historical data. We designed and implemented a simulation engine which can produce Monte Carlo scenarios for the hourly day-ahead values of load, and solar and wind power productions at the spatial and temporal resolutions of the historical data used to fit the model. Finally we developed an open-source Python package which can, from an input grid model, efficiently use forecasts and large numbers of Monte Carlo scenarios to provide unit commitment and economic dispatch for each of these scenarios. The high dimensional stochastic model and the subsequent Monte Carlo simulation engine were implemented in the package PGscen and the corresponding UC and ED optimization programs in the package Vatic.

14 SOLAR ENERGY↗
NASA NTRS1993

Simulated orbiter lift-off acceleration - A stochastic model

A simple procedure for generating synthetic orbiter liftoff acceleration is described. This stochastic method is based on a filtered white noise model which accounts for the nonstationary time evolution of amplitude and frequency content of orbiter lift-off accelerations. The example of Space Transportation System STS-41 Z-acceleration is carefully examined, and the corresponding synthetic orbiter lift-off accelerations model was developed. It is shown that the resulting simulated accelerations preserve significant properties of the original record. Statistically based acceleration response spectra curves for the simulated orbiter lift-off excitations are also evaluated and the results are compared with those for the actual flight data.

Lee-Glauser, Gina↗
DOE OSTI2016

Statistical characteristic in time-domain of direct current corona-generated audible noise from conductor in corona cage

The corona-generated audible noise (AN) has become one of decisive factors in the design of high voltage direct current (HVDC) transmission lines. The AN from transmission lines can be attributed to sound pressure pulses which are generated by the multiple corona sources formed on the conductor, i.e., transmission lines. In this paper, a detailed time-domain characteristics of the sound pressure pulses, which are generated by the DC corona discharges formed over the surfaces of a stranded conductors, are investigated systematically in a laboratory settings using a corona cage structure. The amplitude of sound pressure pulse and its time intervals are extracted by observing a direct correlation between corona current pulses and corona-generated sound pressure pulses. Based on the statistical characteristics, a stochastic model is presented for simulating the sound pressure pulses due to DC corona discharges occurring on conductors. The proposed stochastic model is validated by comparing the calculated and measured A-weighted sound pressure level (SPL). The proposed model is then used to analyze the influence of the pulse amplitudes and pulse rate on the SPL. Furthermore, a mathematical relationship is found between the SPL and conductor diameter, electric field, and radial distance.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗