Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Stochastic simulations”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

Electron acceleration in stochastic double layers

Transversely localized double layers evolve randomly in turbulent regions of strongly magnetized plasma carrying current along the magnetic field. Results from numerical simulations and spacecraft observations in the auroral plasma indicate that the parallel electric field in such regions is microscopically intermittent or stochastic. The implications of stochastic double layer fields on electron acceleration will be discussed in terms of a statistical process involving ensemble averages over test particle motion. A Fokker-Planck equation can be derived for the electron phase space density, which depends on the mean and rms amplitudes of the double layers, the mean double layer density, and the initial electron velocity distribution. It is shown that the resulting electron acceleration is very sensitive to the ratio of the initial electron energy to the rms double layer amplitude. When this ratio is large, the acceleration process differs little from that expected in a dc electric field. When it is small, stochastic heating competes with directed acceleration. Evidence for both cases can be found in the auroral ionosphere in association with so-called inverted-V precipitation and collimated edge precipitation.

Lotko, William↗
DOE OSTI2022

Effects of Improved Simulation of Precipitation on Evapotranspiration and Its Partitioning Over Land

Evapotranspiration (ET) is a key component of the global hydrological cycle, which is strongly modulated by the occurrence of different rainfall intensities. Global climate models (GCMs) commonly suffer from “too much light rain” and a negative bias in the ratio of transpiration (T) to ET (T/ET). It is unclear whether these biases are related. Here, in this work, we show that with the improved simulation of probability density functions of rainfall intensity by suppressing light-rain occurrence using a stochastic convection parameterization in the NCAR CESM1.2, the canopy T increases in tropical forests while evaporation from canopy interception and bare soil decreases. The simulated T/ET is increased by 2.5% globally and up to 8% regionally, primarily attributable to reduced fraction of wet leaves due to less frequent light rain despite its weak intensity. These results imply that excessive light rain is an important cause of the negative T/ET bias in GCMs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Stochastic Maxwell-Bloch equations for modeling amplified spontaneous emission

An approach for studying atom-radiation interaction has been developed, associating quantum operators with stochastic variables governed by discrete Heisenberg equations. This framework models general multilevel atomic systems using the paraxial approximation in both single- and multipass configurations. Simulation results for x-ray laser oscillators with parameters of interest are presented. Finally, the stochastic model is efficient to solve, as the required computational resources scale linearly with the number of emitters, and it may be applicable to various other quantum systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2005

A Comprehensive Study of Three Delay Compensation Algorithms for Flight Simulators

This paper summarizes a comprehensive study of three predictors used for compensating the transport delay in a flight simulator; The McFarland, Adaptive and State Space Predictors. The paper presents proof that the stochastic approximation algorithm can achieve the best compensation among all four adaptive predictors, and intensively investigates the relationship between the state space predictor s compensation quality and its reference model. Piloted simulation tests show that the adaptive predictor and state space predictor can achieve better compensation of transport delay than the McFarland predictor.

Guo, Liwen↗
DOE OSTI2024

Bounds on galaxy stochasticity from halo occupation distribution modeling

The joint probability distribution of matter overdensity and galaxy counts in cells is a powerful probe of cosmology, and the extent to which variance in galaxy counts at fixed matter density deviates from Poisson shot noise is not fully understood. The lack of informed bounds on this stochasticity is currently the limiting factor in constraining cosmology with the galaxy–matter probability distribution function (PDF). We investigate stochasticity in the conditional distribution of galaxy counts along lines of sight with fixed matter density, and we present a halo occupation distribution (HOD)-based approach for obtaining plausible ranges for stochasticity parameters. To probe the high-dimensional space of possible galaxy–matter connections, we derive a set of HODs that conserve the galaxies’ linear bias and number density to produce RED M A G I C-like galaxy catalogs within the A BACUS S UMMIT suite of N -body simulations. We study the impact of individual HOD parameters and cosmology on stochasticity and perform a Monte Carlo search in HOD parameter space subject to the constraints on bias and density. In mock catalogs generated by the selected HODs, shot noise in galaxy counts spans both sub-Poisson and super-Poisson values, ranging from 80% to 133% of Poisson variance for cells with mean matter density. Nearly all of the derived HODs show a positive relationship between local matter density and stochasticity. For galaxy catalogs with higher stochasticity, modeling galaxy bias to second order is required for an accurate description of the conditional PDF of galaxy counts at fixed matter density. The presence of galaxy assembly bias also substantially extends the range of stochasticity in the super-Poisson direction. This HOD-based approach leverages degrees of freedom in the galaxy–halo connection to obtain informed bounds on nuisance model parameters and can be adapted to study other parametrizations of shot noise in galaxy counts, in particular to motivate prior ranges on stochasticity for cosmological analyses.

Britt, Dylan (ORCID:000000019905601X)↗
DOE OSTI2022

ASPIRES: Airport Shuttle Planning and Improved Routing Event-driven Simulation

Most of the existing traffic simulation packages require significant calibration work to be able to reflect reality. To evaluate special operations including emerging technologies, a microscopic simulation that tracks detailed interactions of all the elements of the traffic systems is usually needed. This type of simulation is usually computationally demanding. This work developed an Airport Shuttle Planning and Improved Routing Event-driven Simulation (ASPIRES) package to simulate and evaluate current, potential, and future airport shuttle operations. Here, the simulation was driven by data and thus did not require much calibration effort. The discrete-event simulation nature of ASPIRES makes the simulation computationally efficient. Simulating 1 day of shuttle operations takes less than 2s. The study site of this work is the Dallas/Fort Worth International Airport in the U.S. The shuttle service that connects the five terminals of the airport and the rental car center was studied. Travel times, dwell times, and passenger arrivals were simulated using empirical distributions derived mainly from real data to capture the stochastic nature of the rental car center shuttle bus operations. Data on bus miles traveled, bus energy consumption, passenger wait times, and passengers left behind at stops were collected to study the trade-off between energy use and passenger experience. Electric bus and on-demand bus operations were also included. The simulation outputs can show passengers statistics at terminals, shuttles statistics, and charging station statistics. ASPIRES cannot be used to model a generic traffic system but is well-suited for fleet systems.

33 ADVANCED PROPULSION SYSTEMS↗
NASA NTRS1994

Random vibration of nonlinear beams by the new stochastic linearization technique

In this paper, the beam under general time dependent stationary random excitation is investigated, when exact solution is unavailable. Numerical simulations are carried out to compare its results with those yielded by the conventional linearization techniques. It is found that the modified version of the stochastic linearization technique yields considerably more accurate results for the mean square displacement of the beam than the conventional equivalent linearization technique, especially in the case of large nonlinearity.

Fang, J.↗
DOE OSTI2021

Tutorial: Langevin Dynamics methods for aerosol particle trajectory simulations and collision rate constant modeling

The Langevin Dynamics (LD) method (also known in the literature as Brownian Dynamics) is routinely used to simulate aerosol particle trajectories for transport rate constant calculations as well as to understand aerosol particle transport in internal and external fluid flows. This tutorial intends to explain the methodological details of setting up a LD simulation of a population of aerosol particles and to deduce rate constants from an ensemble of classical trajectories. We discuss the applicability and limitations of the translational Langevin equation to model the combined stochastic and deterministic motion of particles in fields of force or fluid flow. The drag force and stochastic “diffusion” force terms that appear in the Langevin equation are discussed elaborately, along with a summary of common forces relevant to aerosol systems (electrostatic, gravity, van der Waals, …); a commonly used first order and a fourth order Runge-Kutta time stepping schemes for linear stochastic ordinary differential equations are presented. A MATLAB® implementation of a LD code for simulating particle settling under gravity using the first order scheme is included for illustration. Scaling analysis of aerosol transport processes and the selection of timestep and domain size for trajectory simulations are demonstrated through two specific aerosol processes: particle diffusion charging and coagulation. Fortran® implementations of the first order and fourth order time-stepping schemes are included for simulating the 3D motion of a particle in a periodic domain. Lastly, potential applications and caveats to the usage of LD are included as a summary.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Stochastic finite volume method for uncertainty quantification of transient flow in gas pipeline networks

We develop a weakly intrusive framework to simulate the propagation of uncertainty in solutions of generic hyperbolic partial differential equation systems on graph-connected domains with nodal coupling and boundary conditions. The method is based on the Stochastic Finite Volume (SFV) approach and can be applied for uncertainty quantification (UQ) of the dynamical state of fluid flow over actuated transport networks. The numerical scheme has specific advantages for modeling intertemporal uncertainty in time-varying boundary parameters, which cannot be characterized by strict upper and lower (interval) bounds. We describe the scheme for a single pipe, and then formulate the controlled junction Riemann problem (JRP) that enables the extension to general network structures. In conclusion, we demonstrate the method's capabilities and performance characteristics using a standard benchmark test network.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

A low-cost and realistic noisy light system for studying photosynthesis

Unlike the light conditions commonly used to grow photosynthetic organisms in the research laboratory, the light intensity in real environments is dynamic. A simple and low-cost system is described in which a commercial dimmable LED panel is controlled to simulate a sinusoidal function representing daylight hours and overlaid with stochastic shading events. The output closely resembles light intensity measurements on Earth’s surface on partly cloudy days or in lower levels of plant canopies. Here, this tool may be useful to researchers studying photosynthetic acclimation responses.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

A Machine Learning Assisted Development of a Model for the Populations of Convective and Stratiform Clouds

Abstract Traditional parameterizations of the interaction between convection and the environment have relied on an assumption that the slowly varying large‐scale environment is in statistical equilibrium with a large number of small and short‐lived convective clouds. They fail to capture nonequilibrium transitions such as the diurnal cycle and the formation of mesoscale convective systems as well as observed precipitation statistics and extremes. Informed by analysis of radar observations, cloud‐permitting model simulation, theory, and machine learning, this work presents a new stochastic cloud population dynamics model for characterizing the interactions between convective and stratiform clouds, with the goal of informing the representation of these interactions in global climate models. Fifteen wet seasons of precipitating cloud observations by a C‐band radar at Darwin, Australia are fed into a machine learning algorithm to obtain transition functions that close a set of coupled equations relating large‐scale forcing, mass flux, the convective cell size distribution, and the stratiform area. Under realistic large‐scale forcing, the derived transition functions show that, on the one hand, interactions with stratiform clouds act to dampen the variability in the size and number of convective cells and therefore in the convective mass flux. On the other, for a given convective area fraction, a larger number of smaller cells is more favorable for the growth of stratiform area than a smaller number of larger cells. The combination of these two factors gives rise to solutions with a few convective cells embedded in a large stratiform area, reminiscent of mesoscale convective systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Concurrent multi-peak Bragg coherent x-ray diffraction imaging of 3D nanocrystal lattice displacement via global optimization

Abstract In this paper we demonstrated a method to reconstruct vector-valued lattice distortion fields within nanoscale crystals by optimization of a forward model of multi-reflection Bragg coherent diffraction imaging (MR-BCDI) data. The method flexibly accounts for geometric factors that arise when making BCDI measurements, is amenable to efficient inversion with modern optimization toolkits, and allows for globally constraining a single image reconstruction to multiple Bragg peak measurements. This is enabled by a forward model that emulates the multiple Bragg peaks of a MR-BCDI experiment from a single estimate of the 3D crystal sample. We present this forward model, we implement it within the stochastic gradient descent optimization framework, and we demonstrate it with simulated and experimental data of nanocrystals with inhomogeneous internal lattice displacement. We find that utilizing a global optimization approach to MR-BCDI affords a reliable path to convergence of data which is otherwise challenging to reconstruct.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Charging-driven coarsening and melting of a colloidal nanoparticle monolayer at an ionic liquid–vacuum interface

In this work, we induce and investigate the coarsening and melting dynamics of an initially static nanoparticle colloidal monolayer at an ionic liquid–vacuum interface, driven by a focused, scanning electron beam. Coarsening occurs through grain interface migration and larger-scale motions such as grain rotations, often facilitated by sliding dislocations. The progressive decrease in area fraction that drives melting of the monolayer is explained using an electrowetting model whereby particles at the interface are solvated once their accumulating charge recruits sufficient counterions to subsume the particle. Subject to stochastic particle removal from the monolayer, melting is recapitulated in simulations with a Lennard-Jones potential. This new driving mechanism for colloidal systems, whose dynamical timescales we show can be controlled with the accelerating voltage, opens the possibility to manipulate particle interactions dynamically without need to vary particle intrinsic properties or surface treatments. Furthermore, the decrease in particle size availed by electron imaging presents opportunities to observe force and time scales in a lesser-explored regime intermediate between typical colloidal and molecular systems.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Simulating pitch angle scattering using an explicitly solvable energy-conserving algorithm

Particle distribution functions evolving under the Lorentz operator can be simulated with the Langevin equation for pitch angle scattering. This approach is frequently used in particle based Monte-Carlo simulations of plasma collisions, among others. However, most numerical treatments do not guarantee energy conservation, which may lead to unphysical artifacts such as numerical heating and spectra distortions. In this paper, we present a novel structure-preserving numerical algorithm for the Langevin equation for pitch angle scattering. Similar to the well-known Boris algorithm, the proposed numerical scheme takes advantage of the structure-preserving properties of the Cayley transform when calculating the velocity-space rotations. The resulting algorithm is explicitly solvable, while preserving the norm of velocities down to machine precision. We demonstrate that the method has the same order of numerical convergence as the traditional stochastic Euler-Maruyama method. The numerical scheme is benchmarked by simulating the pitch angle scattering of a particle beam, and comparing with the analytical solution. Benchmark results show excellent agreement with theoretical predictions, showcasing the remarkable long time accuracy of the proposed algorithm.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Reliability Analysis of Power Grids Considering Component Failures of Variable Energy Resources

This paper proposes an improved model for the reliability assessment of power systems considering component failures of variable energy resources (VER). The inherent intermittency of VER such as solar photovoltaic (PV) and wind farms, along with their susceptibility to component failures, present significant challenges to reliable system operation. These issues, combined with power grid operation and network constraints, complicate the reliable operation of VER-integrated power systems. Here, to address these concerns, this paper introduces a reliability assessment framework that considers VER input variability, its impact on component availability, and their resulting impact on overall system reliability. Stochastic models based on discrete Markov processes are developed to incorporate variable irradiance, wind speeds, and their effects on PV and wind component failure rates. A next-event and state transition-based approach is then developed to integrate the stochastic models into a mixed-timing sequential Monte Carlo simulation framework for composite reliability assessment. Case studies on the RTS-GMLC system demonstrate the effectiveness of the proposed model in evaluating the reliability of VER-integrated systems.

Pandit, Dilip [Sandia National Laboratories (SNL-N↗
DOE OSTI2023

Reducing Urban Traffic Congestion Using Deep Learning and Model Predictive Control

This article proposes a deep learning (DL)-based control algorithm—DL velocity-based model predictive control (VMPC)—for reducing traffic congestion with slowly time-varying traffic signal controls. This control algorithm consists of system identification using DL and traffic signal control using VMPC. For the training process of DL, we established a modeling error entropy loss as the criteria inspired by the theory of stochastic distribution control (SDC) originated by the fourth author. Herein, simulation results show that the proposed algorithm can reduce traffic congestion with a slowly varying traffic signal control input. Results of an ablation study demonstrate that this algorithm compares favorably to other model-based controllers in terms of prediction error, signal varying speed, and control effectiveness.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Utah FORGE 5-2428: Fracture Permeability Impact on Reservoir Stress and Seismic Slip Behavior - Workshop Presentation

This is a presentation on the Fracture Permeability Impact on Seismic Slip Behavior project by Lawrence Livermore National Laboratory, presented by Dr. Kayla A. Kroll. The project's objective is to develop, apply and validate a holistic thermal, hydrologic, mechanical, and chemical (THMC) workflow that includes evaluation of induced seismic slip in EGS reservoirs. This presentation was featured in the Utah FORGE R&D Annual Workshop on September 8, 2023. The workshop provided a valuable opportunity to explore the progress made in each of the 17 Research and Development projects funded under Solicitation 2020-1 which aim to enhance our understanding of the crucial factors influencing the development of Enhanced Geothermal Systems (EGS) reservoirs and resources.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

Understanding and Strategies for Controlled Interfacial Phenomena in Lithium-Ion Batteries and Beyond

Electrolyte chemistry and properties and electrode structure and chemical properties particularly at the interfaces are crucial for the development of advanced battery components. This project focused on elucidating the roles of the electrolyte and electrode on the formation and evolution of the SEI layer and cell electrochemical performance on silicon and on Li metal anodes. For Si nano/microstructures we evaluated lithiation, volume changes, reactivity, and chemo-mechanical transformations as functions of nanoparticle size, shape, presence of coatings, electrolyte composition, electron leakage to the electrolyte, and cycling. For Li metal anodes we provided detailed analyses of the electrode and electrolyte effects on interfacial reactivity and as a function of microstructural evolution and underlying stochasticity as dendrites nucleate and grow. We thoroughly investigated and simulated how the chemistry of the various components of the electrolyte and the electrode architecture may affect the electrochemical reactions as well as cell degradation. A multiscale modeling approach was utilized, where atomistic simulations informed about the microscopic behavior of the system, and intermediate time and length scales were investigated with mesoscopic models. Results were also tested against selected experiments from our collaborators. The combined theoretical-experimental strategy saves costs by utilizing predictions from computational analyses to guide experimentation.

25 ENERGY STORAGE↗