SEARCH · Engineering Papers

Results for “Stochastic systems”

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.

SourceResearch connectionFrom yearThrough year

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

At least 1,261 records · Page 70

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Stochastic density functional theory: Real- and energy-space fragmentation for noise reduction

Stochastic density functional theory (sDFT) is becoming a valuable tool for studying ground-state properties of extended materials. The computational complexity of describing the Kohn-Sham orbitals is replaced by introducing a set of random (stochastic) orbitals leading to linear and often sub-linear scaling of certain ground-state observables at the account of introducing a statistical error. Schemes to reduce the noise are essential, for example, for determining the structure using the forces obtained from sDFT. Recently, we have introduced two embedding schemes to mitigate the statistical fluctuations in the electron density and resultant forces on the nuclei. Both techniques were based on fragmenting the system either in real space or slicing the occupied space into energy windows, allowing for a significant reduction in the statistical fluctuations. For chemical accuracy, further reduction of the noise is required, which could be achieved by increasing the number of stochastic orbitals. However, the convergence is relatively slow as the statistical error scales as 1 / N χ according to the central limit theorem, where N χ is the number of random orbitals. In this paper, we combined the embedding schemes mentioned above and introduced a new approach that builds on overlapped fragments and energy windows. The new approach significantly lowers the noise for ground-state properties, such as the electron density, total energy, and forces on the nuclei, as demonstrated for a G-center in bulk silicon.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Thermal fluctuations in the Landau-Lifshitz-Bloch model

A formulation for thermal noise in the stochastic form of the Landau-Lifshitz-Bloch equation used for modeling the magnetization dynamics at elevated temperatures is presented here. The diffusion coefficients for thermal fluctuations are obtained via the Fokker-Plank equation using the mean-field approximation of the field in defining the free energy. The presented model leads to a mean magnetization consistent with the equilibrium magnetization for small and large particles. The distribution of the magnetization magnitude is of the Maxwell-Boltzmann type. The presented model was tested by studying the equilibrium magnetization in macrospin particles at high temperatures. The model is appealing for multiscale modeling, such as modeling heat assisted magnetic recording systems and all-optical magnetization reversal.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRSDate not supplied

Supersonic Free-Flight Dynamics Testing in the Stratosphere

A new technique for obtaining stratospheric free-flight dynamics data for atmospheric entry capsules is described. The Stratospheric Projectile Experiment of Entry Dynamics (SPEED) represents a new approach to characterizing the free-flight dynamics of vehicles in the supersonic and transonic regime of flight. The SPEED test architecture leverages a stratospheric balloon and 3D-printed flight system in a novel two-stage configuration to deliver test articles to supersonic conditions in the atmosphere which achieve dynamic similitude with a full-scale vehicle. Designed to address the limitations of existing test facilities, SPEED captures the complete time evolution of a vehicle’s dynamic state while producing a statistically significant number of observed flight trajectories to address the stochastic nature of the wake-driven dynamic stability phenomenon. SPEED was developed over two years at NASA Ames Research Center. The inaugural flight of the SPEED test platform was conducted in the summer of 2024 where scaled capsules of the Mars Sample Return Earth Entry System and Dragonfly entry vehicle were tested. This paper describes the motivation for a new test architecture, operational constraints and achievable flight test envelopes, the design and development of the SPEED concept, and results from the demonstration flight. Techniques for estimating the dynamic aerodynamic characteristics and atmospheric conditions for each test article individually and as a collective are also discussed.

Entry Vehicle↗

DOE OSTI2021

Z 2 topological order and first-order quantum phase transitions in systems with combinatorial gauge symmetry

We study a generalization of the two-dimensional transverse-field Ising model, combining both ferromagnetic and antiferromagnetic two-body interactions, that hosts exact global and local Z 2 gauge symmetries. Using exact diagonalization and stochastic series expansion quantum Monte Carlo methods, we confirm the existence of the topological phase in line with previous theoretical predictions. Our simulation results show that the transition between the confined topological phase and the deconfined paramagnetic phase is of first order, in contrast to the conventional Z 2 lattice gauge model in which the transition maps onto that of the standard Ising model and is continuous. We further generalize the model by replacing the transverse field on the gauge spins with a ferromagnetic X X interaction while keeping the local gauge symmetry intact. We find that the Z 2 topological phase remains stable, while the paramagnetic phase is replaced by a ferromagnetic phase. The topological-ferromagnetic quantum phase transition is also of first order. For both models, we discuss the low-energy spinon and vison excitations of the topological phase and their avoided level crossings associated with the first-order quantum phase transitions.

2-dimensional systems↗

DOE OSTI2025

Data-Driven Compositional Optimization in Misspecified Regimes

With a manifold growth in the scale and intricacy of systems, the challenges of parametric misspecification become pronounced. These concerns are further exacerbated in compositional settings, which emerge in problems complicated by modeling risk and robustness. In “Data-Driven Compositional Optimization in Misspecified Regimes,” the authors consider the resolution of compositional stochastic optimization problems, plagued by parametric misspecification. In considering settings where such misspecification may be resolved via a parallel learning process, the authors develop schemes that can contend with diverse forms of risk, dynamics, and nonconvexity. They provide asymptotic and rate guarantees for unaccelerated and accelerated schemes for convex, strongly convex, and nonconvex problems in a two-level regime with extensions to the multilevel setting. Surprisingly, the nonasymptotic rate guarantees show no degradation from the rate statements obtained in a correctly specified regime and the schemes achieve optimal (or near-optimal) sample complexities for general T-level strongly convex and nonconvex compositional problems.

Business & Economics↗

DOE OSTI2024

Collisional simulations of the modulator section in coherent electron cooling

The first section of any coherent electron cooling (CeC) system is the modulator, where the density of the electron beam is modulated by the copropagating ion beam. This density modulation is a result of Coulomb collisions between the individual particles of the two beams. The pairwise, stochastic part of the interactions impacts the overall performance of the CeC process. We present the first simulations of the density modulations of the electron beams from a collisional picture of the dynamics, considering the proof-of-principle CeC experiments at Brookhaven National Laboratory. These simulations were performed using PHAD, which is the first efficient, large-scale collisional numerical method in beam physics that we have previously developed and benchmarked. Realistic beam distributions and external fields have been optimized to provide strong modulation signals necessary for variations of coherent electron cooling systems. Cooling performance limits and potential collisionless simulation pitfalls are pointed out. Published by the American Physical Society 2024

43 PARTICLE ACCELERATORS↗

DOE OSTI2020

Multigrid deflation for Lattice QCD

Computing the trace of the inverse of large matrices is typically addressed through statistical methods. Deflating out the lowest eigenvectors or singular vectors of the matrix reduces the variance of the trace estimator. This work summarizes our efforts to reduce the computational cost of computing the deflation space while achieving the desired variance reduction for Lattice QCD applications. Previous efforts computed the lower part of the singular spectrum of the Dirac operator by using an eigensolver preconditioned with a multigrid linear system solver. Despite the improvement in performance in those applications, as the problem size grows the runtime and storage demands of this approach will eventually dominate the stochastic estimation part of the computation. In this work, we propose to compute the deflation space in one of the following two ways. First, by using an inexact eigensolver on the Hermitian, but maximally indefinite, operator. Second, by exploiting the fact that the multigrid prolongator for this operator is rich in components toward the lower part of the singular spectrum. We show experimentally that the inexact eigensolver can approximate the lower part of the spectrum even for ill-conditioned operators. Also, the deflation based on the multigrid prolongator is more efficient to compute and apply, and, despite its limited ability to approximate the fine level spectrum, it obtains similar variance reduction on the trace estimator as deflating with approximate eigenvectors from the fine level operator.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

A Kaczmarz-inspired approach to accelerate the optimization of neural network wavefunctions

Neural network wavefunctions optimized using the variational Monte Carlo method have been shown to produce highly accurate results for the electronic structure of atoms and small molecules, but the high cost of optimizing such wavefunctions prevents their application to larger systems. We propose the Subsampled Projected-Increment Natural Gradient Descent (SPRING) optimizer to reduce this bottleneck. SPRING combines ideas from the recently introduced minimum-step stochastic reconfiguration optimizer (MinSR) and the classical randomized Kaczmarz method for solving linear least-squares problems. We demonstrate that SPRING outperforms both MinSR and the popular Kronecker-Factored Approximate Curvature method (KFAC) across a number of small atoms and molecules, given that the learning rates of all methods are optimally tuned. For example, on the oxygen atom, SPRING attains chemical accuracy after forty thousand training iterations, whereas both MinSR and KFAC fail to do so even after one hundred thousand iterations.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

High-throughput bend-strengths of ultra-small polysilicon MEMS components

The strength distribution of polysilicon bend specimens, approximately 10 μm in size, is measured using a high-throughput microelectromechanical system fabrication and testing method. The distribution is predicted from reference tests on tensile specimens and finite element analysis of the bend specimen geometry incorporated into a stochastic extreme-value strength framework. Agreement between experiment and prediction suggests that the ultra-small specimens may be at the limit of extreme-value scaling and contain only one strength-controlling flaw/specimen.

Cook, Robert F. (ORCID:0000000304228881)↗

DOE OSTI2021

Effect of diffusion constant on the morphology of dendrite growth in lithium metal batteries

Lithium dendrites can lead to a short circuit and battery failure, and developing strategies for their suppression is of considerable importance. In this work, we study the growth of dendrites in a simple model system where the solvent is a continuum and the lithium ions are hard spheres that can deposit by sticking to existing spheres or the electrode surface. Using stochastic dynamics simulations, we investigate the effect of applied voltage and diffusion constant on the growth of dendrites. We find that the diffusion constant is the most significant factor, and the inhomogeneity of the electric field does not play a significant role. Additionally, the growth is most pronounced when the applied voltage and diffusion constant are both low. We observe a structural change from broccoli to cauliflower shape as the diffusion constant is increased. The simulations suggest that a control of electrolyte parameters that impact lithium diffusion might be an attractive route to controlling dendrite growth.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

On-lattice kinetic Monte Carlo approaches for modeling molecular anisotropy in resveratrol crystallization

Stilbenes are a class of organic compounds with broad-ranging pharmaceutical and agricultural applications, which are typically isolated and purified through recrystallization. We are motivated by reducing experimental waste and optimizing yield via developing predictive simulations for processing-dependent crystal morphologies. Using resveratrol as a model stilbene system, we have developed an approach for simulating crystallization with molecular resolution using on-lattice kinetic Monte Carlo. In this work, we highlight modifications to the Stochastic Parallel PARticle Kinetic Simulator (SPPARKS) software package, which were essential to this application. Key enhancements include the incorporation of non-orthogonal cell shapes and monomer anisotropy approximations using bound hard spheres. This new SPPARKS application has been applied to resveratrol with attachment energy libraries obtained from density functional theory, resulting in excellent agreement with experimental morphology prediction.

crystallization↗

DOE OSTI2023

Joint Resource Modeling and Assessment for Hybrid Distributed Solar and Wind Systems

The inherent variability and uncertainty in distributed energy resources can presents myriad challenges to the planning and operations of power systems. These risks are poised to become larger as the penetration of renewable energy sources rises in the power generation mix. Hybrid solar-wind energy systems are able to mitigate some of these risks by their complementary resource availability. Surface solar and wind fields are coupled and correlated in both space and time. Appropriately estimating the hybrid solar wind energy system requires simulating the spatio-temporal structure of these fields that can be produced for each time horizon. We introduce a novel joint spatio-temporal stochastic differential equation (SPDE) approach that captures the spatio-temporal dynamics of solar and wind fields and their joint dependency over a domain for each time step. In the case study on Colorado, we consider nonstationary three-level hierarchical spatio temporal models for both hourly solar irradiance data and wind speed data in Colorado. Dependence between the solar irradiance data and wind speed data is captured by a shared spatio-temporal random effect. Our approach performs well in terms of the prediction score criterion.

joint modeling↗

NASA NTRS1976

Trends in modern system theory

The topics considered are related to linear control system design, adaptive control, failure detection, control under failure, system reliability, and large-scale systems and decentralized control. It is pointed out that the design of a linear feedback control system which regulates a process about a desirable set point or steady-state condition in the presence of disturbances is a very important problem. The linearized dynamics of the process are used for design purposes. The typical linear-quadratic design involving the solution of the optimal control problem of a linear time-invariant system with respect to a quadratic performance criterion is considered along with gain reduction theorems and the multivariable phase margin theorem. The stumbling block in many adaptive design methodologies is associated with the amount of real time computation which is necessary. Attention is also given to the desperate need to develop good theories for large-scale systems, the beginning of a microprocessor revolution, the translation of the Wiener-Hopf theory into the time domain, and advances made in dynamic team theory, dynamic stochastic games, and finite memory stochastic control.

Athans, M.↗

NASA NTRS1983

Approximation of the optimal compensator for a large space structure

This paper considers the approximation of the optimal compensator for a Large Space Structure. The compensator is based upon a solution to the Linear Stochastic Quadratic Regulator problem. Colocation of sensors and actuators is assumed. A small gain analytical solution for the optimal compensator is obtained for a single input/single output system, i.e., certain terms in the compensator can be neglected for sufficiently small gain. The compensator is calculated in terms of the kernel to a Volterra integral operator using a Neumann series. The calculation of the compensator is based upon the C sub 0 semigroup for the infinite dimensional system. A finite dimensional approximation of the compensator is, therefore, obtained through analysis of the infinite dimensional compensator which is a compact operator.

Mackay, M. K.↗

NASA NTRS2006

Chance of Necessity: Modeling Origins of Life

The fundamental nature of processes that led to the emergence of life has been a subject of long-standing debate. One view holds that the origin of life is an event governed by chance, and the result of so many random events is unpredictable. This view was eloquently expressed by Jacques Monod in his book Chance or Necessity. In an alternative view, the origin of life is considered a deterministic event. Its details need not be deterministic in every respect, but the overall behavior is predictable. A corollary to the deterministic view is that the emergence of life must have been determined primarily by universal chemistry and biochemistry rather than by subtle details of environmental conditions. In my lecture I will explore two different paradigms for the emergence of life and discuss their implications for predictability and universality of life-forming processes. The dominant approach is that the origin of life was guided by information stored in nucleic acids (the RNA World hypothesis). In this view, selection of improved combinations of nucleic acids obtained through random mutations drove evolution of biological systems from their conception. An alternative hypothesis states that the formation of protocellular metabolism was driven by non-genomic processes. Even though these processes were highly stochastic the outcome was largely deterministic, strongly constrained by laws of chemistry. I will argue that self-replication of macromolecules was not required at the early stages of evolution; the reproduction of cellular functions alone was sufficient for self-maintenance of protocells. In fact, the precise transfer of information between successive generations of the earliest protocells was unnecessary and could have impeded the discovery of cellular metabolism. I will also show that such concepts as speciation and fitness to the environment, developed in the context of genomic evolution also hold in the absence of a genome.

Pohorille, Andrew↗

DOE OSTI2025

Realizing Steady-State Microbunching with Optical Stochastic Crystallization

Optical Stochastic Cooling (OSC) is a state-of-the-art beam cooling technology first demonstrated in 2021 at the IOTA storage ring at Fermilab's FAST facility. A second phase of the research program is planned to run in 2026 and will incorporate an optical amplifier to enable significantly increased cooling rates and greater operational flexibility. In addition to beam cooling, an OSC system can be configured to enable advanced control over the phase space of the beam. An example operational mode could enable crystallization, where the particles in a bunch are locked into a self-reinforcing, regular microstructure at the OSC fundamental wavelength; we refer to this as Optical Stochastic Crystallization (OSX). OSX represents a new path toward Steady-State Microbunching (SSMB), which may enable light sources combining the high brightness of a free-electron laser with the high repetition rate of a storage ring. Such a source has applications from the terahertz to the extreme ultraviolet (EUV), including high-power EUV generation for semiconductor lithography. This contribution will discuss the integration of OSX development as part of the OSC program at IOTA. The design of an accelerator lattice to enable the mechanism and associated high fidelity simulations will be shown, and a path to realizing an experimental demonstration will be discussed.

Wallbank, Michael James [Fermilab]↗

DOE OSTI2025

Steady-State Microbunching using Optical Stochastic Cooling

Optical Stochastic Cooling (OSC) is a state-of-the-art beam cooling technology first demonstrated in 2021 at the IOTA storage ring at Fermilab's FAST facility. A second phase of the research program is planned to run in 2026 and will incorporate an optical amplifier to enable significantly increased cooling rates and greater operational flexibility. In addition to beam cooling, an OSC system can be configured to enable advanced control over the phase space of the beam. An example operational mode could enable crystallization, where the particles in a bunch are locked into a self-reinforcing, regular microstructure at the OSC fundamental wavelength; we refer to this as Optical Stochastic Crystallization (OSX). OSX represents a new path toward Steady-State Microbunching (SSMB), which may enable light sources combining the high brightness of a free-electron laser with the high repetition rate of a storage ring. Such a source has applications from the terahertz to the extreme ultraviolet (EUV), including high-power EUV generation for semiconductor lithography. This contribution will discuss the integration of OSX development as part of the OSC program at IOTA. The design of an accelerator lattice to enable the mechanism and associated high fidelity simulations demonstrating the beam dynamics will be shown, and a path to realizing an experimental demonstration will be discussed.

Wallbank, Michael James [Fermilab]↗

DOE OSTI2024

Realizing Steady-State Microbunching with Optical Stochastic Crystallization

Optical Stochastic Cooling (OSC) is a state-of-the-art beam cooling technology first demonstrated in 2021 at the IOTA storage ring at Fermilab's FAST facility. A second phase of the research program is planned to run in early 2025 and will incorporate an optical amplifier to enable significantly increased cooling rates and greater operational flexibility. In addition to beam cooling, an OSC system can be configured to enable advanced control over the phase space of the beam. An example operational mode could enable crystallization, where the particles in a bunch are locked into a self-reinforcing, regular microstructure at the OSC fundamental wavelength; we refer to this as Optical Stochastic Crystallization (OSX). OSX represents a new path toward Steady-State Microbunching (SSMB), which may enable light sources combining the high brightness of a free-electron laser with the high repetition rate of a storage ring. Such a source has applications from the terahertz to the extreme ultraviolet (EUV), including high-power EUV generation for semiconductor lithography. This contribution will discuss the status of the OSC experimental program and its potential to achieve the first demonstration of SSMB during the upcoming experimental run.

Wallbank, M.↗