Search NASA⌕ Search

SEARCH · Search NASA

Results for “Stochastic dynamical 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.

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

At least 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

A minimally invasive, efficient method for propagation of full-field uncertainty in solid dynamics

In this work, we present a minimally invasive method for forward propagation of material property uncertainty to full-field quantities of interest in solid dynamics. Full-field uncertainty quantification enables the design of complex systems where quantities of interest, such as failure points, are not known a priori. The method, motivated by the well-known probability density function (PDF) propagation method of turbulence modeling, uses an ensemble of solutions to provide the joint PDF of desired quantities at every point in the domain. A small subset of the ensemble is computed exactly, and the remainder of the samples are computed with approximation of the evolution equations based on those exact solutions. Although the proposed method has commonalities with traditional interpolatory stochastic collocation methods applied directly to quantities of interest, it is distinct and exploits the parameter dependence and smoothness of the driving term of the evolution equations. The implementation is model independent, storage and communication efficient, and straightforward. We demonstrate its efficiency, accuracy, scaling with dimension of the parameter space, and convergence in distribution with two problems: a quasi-one-dimensional bar impact, and a two material notched plate impact. For the bar impact problem, we provide an analytical solution to PDF of the solution fields for method validation. With the notched plate problem, we also demonstrate good parallel efficiency and scaling of the method.

42 ENGINEERING↗
DOE OSTI2023

Regional biogeography versus intra-annual dynamics of the root and soil microbiome

Abstract Background Root and soil microbial communities constitute the below-ground plant microbiome, are drivers of nutrient cycling, and affect plant productivity. However, our understanding of their spatiotemporal patterns is confounded by exogenous factors that covary spatially, such as changes in host plant species, climate, and edaphic factors. These spatiotemporal patterns likely differ across microbiome domains (bacteria and fungi) and niches (root vs. soil). Results To capture spatial patterns at a regional scale, we sampled the below-ground microbiome of switchgrass monocultures of five sites spanning > 3 degrees of latitude within the Great Lakes region. To capture temporal patterns, we sampled the below-ground microbiome across the growing season within a single site. We compared the strength of spatiotemporal factors to nitrogen addition determining the major drivers in our perennial cropping system. All microbial communities were most strongly structured by sampling site, though collection date also had strong effects; in contrast, nitrogen addition had little to no effect on communities. Though all microbial communities were found to have significant spatiotemporal patterns, sampling site and collection date better explained bacterial than fungal community structure, which appeared more defined by stochastic processes. Root communities, especially bacterial, were more temporally structured than soil communities which were more spatially structured, both across and within sampling sites. Finally, we characterized a core set of taxa in the switchgrass microbiome that persists across space and time. These core taxa represented < 6% of total species richness but > 27% of relative abundance, with potential nitrogen fixing bacteria and fungal mutualists dominating the root community and saprotrophs dominating the soil community. Conclusions Our results highlight the dynamic variability of plant microbiome composition and assembly across space and time, even within a single variety of a plant species. Root and soil fungal community compositions appeared spatiotemporally paired, while root and soil bacterial communities showed a temporal lag in compositional similarity suggesting active recruitment of soil bacteria into the root niche throughout the growing season. A better understanding of the drivers of these differential responses to space and time may improve our ability to predict microbial community structure and function under novel conditions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Stochastic and mixed density functional theory within the projector augmented wave formalism for simulation of warm dense matter

Stochastic density functional theory (DFT) and mixed stochastic-deterministic DFT are burgeoning approaches for the calculation of the equation of state and transport properties in materials under extreme conditions. In the intermediate warm dense matter regime, a state between correlated condensed matter and kinetic plasma, electrons can range from being highly localized around nuclei to delocalized over the whole simulation cell. The plane-wave basis pseudopotential approach is thus the typical tool of choice for modeling such systems at the DFT level. Unfortunately, stochastic DFT methods scale as the square of the maximum plane-wave energy in this basis. To reduce the effect of this scaling and improve the overall description of the electrons within the pseudopotential approximation, we present stochastic and mixed DFT approaches developed and implemented within the projector augmented wave formalism. In conclusion, we compare results between the different DFT approaches for both single-point and molecular dynamics trajectories and present calculations of self-diffusion coefficients of solid density carbon from 1 to 50 eV.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Ab initio molecular dynamics on quantum computers

Ab initio molecular dynamics (AIMD) is a valuable technique for studying molecules and materials at finite temperatures where the nuclei evolve on potential energy surfaces obtained from accurate electronic structure calculations. In this work, we present an approach to running AIMD simulations on noisy intermediate-scale quantum (NISQ)-era quantum computers. The electronic energies are calculated on a quantum computer using the variational quantum eigensolver (VQE) method. Algorithms for computation of analytical gradients entirely on a quantum computer require quantum fault-tolerant hardware, which is beyond NISQ-era. Therefore, we compute the energy gradients numerically using finite differences, the Hellmann-Feynman theorem, and a correlated sampling technique. This method only requires additional classical calculations of electron integrals for each degree of freedom without any additional computations on a quantum computer beyond the initial VQE run. As a proof of concept, AIMD simulations are demonstrated for the H-2 molecule on IBM quantum devices. In addition, we demonstrate the validity of the method for larger molecules using full configuration interaction wave functions. As quantum hardware and noise mitigation techniques continue to improve, the method can be utilized for studying larger molecular systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Time-resolved turbulent dynamo in a laser plasma

Understanding magnetic-field generation and amplification in turbulent plasma is essential to account for observations of magnetic fields in the universe. A theoretical framework attributing the origin and sustainment of these fields to the so-called fluctuation dynamo was recently validated by experiments on laser facilities in low-magnetic-Prandtl-number plasmas (Pm<1). However, the same framework proposes that the fluctuation dynamo should operate differently when Pm≳1, the regime relevant to many astrophysical environments such as the intracluster medium of galaxy clusters. This paper reports an experiment that creates a laboratory Pm≳1 plasma dynamo. We provide a time-resolved characterization of the plasma’s evolution, measuring temperatures, densities, flow velocities, and magnetic fields, which allows us to explore various stages of the fluctuation dynamo’s operation on seed magnetic fields generated by the action of the Biermann-battery mechanism during the initial drive-laser target interaction. The magnetic energy in structures with characteristic scales close to the driving scale of the stochastic motions is found to increase by almost three orders of magnitude and saturate dynamically. It is shown that the initial growth of these fields occurs at a much greater rate than the turnover rate of the driving-scale stochastic motions. Our results point to the possibility that plasma turbulence produced by strong shear can generate fields more efficiently at the driving scale than anticipated by idealized magnetohydrodynamics (MHD) simulations of the nonhelical fluctuation dynamo; this finding could help explain the large-scale fields inferred from observations of astrophysical systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A Continuum Dislocation Dynamics Crystal Plasticity Approach to Irradiated Body-Centered Cubic α-Iron

Radiation-induced embrittlement of reactor pressure vessel (RPV) steels can potentially limit the operating life of nuclear power plants. Over extended exposure to radiation doses, these body-centered cubic (BCC) irons demonstrate irradiation damage. Here, we present a continuum dislocation density (CDD) crystal plasticity model to capture the interaction among dislocations and self-interstitial atom (SIA) loops in α-iron. We demonstrate the importance of modeling cross slip using a combined stochastic Monte Carlo approach and the role of slip system strength anisotropy in capturing stochastic cross slip interactions. Through these captured interactions, the CDD crystal plasticity model can capture both the stress response and the physical evolution of dislocations on different slip system planes. Single-crystal verification experiments are used to calibrate the CDD crystal plasticity model, and a set of simplified polycrystalline simulations demonstrates the model’s ability to capture the stress response from tensile experiments on α-iron.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Time-resolved turbulent dynamo in a laser plasma

Understanding magnetic-field generation and amplification in turbulent plasma is essential to account for observations of magnetic fields in the universe. A theoretical framework attributing the origin and sustainment of these fields to the so-called fluctuation dynamo was recently validated by experiments on laser facilities in low-magnetic-Prandtl-number plasmas (Pm < 1). However, the same framework proposes that the fluctuation dynamo should operate differently when Pm ≳ 1, the regime relevant to many astrophysical environments such as the intracluster medium of galaxy clusters. This paper reports an experiment that creates a laboratory Pm ≳ 1 plasma dynamo. We provide a time-resolved characterization of the plasma’s evolution, measuring temperatures, densities, flow velocities, and magnetic fields, which allows us to explore various stages of the fluctuation dynamo’s operation on seed magnetic fields generated by the action of the Biermann-battery mechanism during the initial drive-laser target interaction. The magnetic energy in structures with characteristic scales close to the driving scale of the stochastic motions is found to increase by almost three orders of magnitude and saturate dynamically. It is shown that the initial growth of these fields occurs at a much greater rate than the turnover rate of the driving-scale stochastic motions. Here, our results point to the possibility that plasma turbulence produced by strong shear can generate fields more efficiently at the driving scale than anticipated by idealized magnetohydrodynamics (MHD) simulations of the nonhelical fluctuation dynamo; this finding could help explain the large-scale fields inferred from observations of astrophysical systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Global decline in ocean memory over the 21st century

Ocean memory, the persistence of ocean conditions, is a major source of predictability in the climate system beyond weather time scales. We show that ocean memory, as measured by the year-to-year persistence of sea surface temperature anomalies, is projected to steadily decline in the coming decades over much of the globe. This global decline in ocean memory is predominantly driven by shoaling of the upper-ocean mixed layer depth in response to global surface warming, while thermodynamic and dynamic feedbacks can contribute substantially regionally. As the mixed layer depth shoals, stochastic forcing becomes more effective in driving sea surface temperature anomalies, increasing high-frequency noise at the expense of persistent signals. Reduced ocean memory results in shorter lead times of skillful persistence-based predictions of sea surface thermal conditions, which may present previously unknown challenges for predicting climate extremes and managing marine biological resources under climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

The Dynamics of Interstellar Asteroids and Comets within the Galaxy: An Assessment of Local Candidate Source Regions for 1I/’Oumuamua and 2I/Borisov

The low velocity of interstellar asteroid 1I/’Oumuamua with respect to our galaxy’s local standard of rest implies that it is young. Adopting the young age hypothesis, we assess possible origin systems for this interstellar asteroid and for 2I/Borisov, though the latter’s higher speed means it is unlikely to be young. First, their past trajectories are modeled under gravitational scattering by galactic components (“disk heating”) to assess how far back one can trace them. The stochastic nature of disk heating means that a back integration can only expect to be accurate to within 15 pc and 2 km s{sup −1} at −10 Myr, dropping steeply to 400 pc and 10 km s{sup −1} at −100 Myr, sharply limiting our ability to determine a precise origin. Nevertheless, we show that ’Oumuamua’s origin system likely is currently within 1 kpc of Earth, in the local Orion Arm. Second, we back integrate ’Oumuamua’s trajectory in order to assess source regions, emphasizing young systems and moving groups. Though disk heating allows for only a statistical link to source regions, ’Oumuamua passed through a considerable subset of the Carina and Columba moving groups when those groups were forming. This makes them perhaps the most plausible source region, if ’Oumuamua was ejected during planet formation or via intra-cluster interactions. We find three stars in the Ursa Major group, one brown dwarf, and seven other stars to have plausible encounters with 2I/Borisov, within 2 pc and 30 km s{sup −1}. These encounters’ high relative speeds mean none are likely to be the home of 2I/Borisov.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Improving Estimation of the Koopman Operator with Kolmogorov–Smirnov Indicator Functions

It has become common to perform kinetic analysis using approximate Koopman operators that transform high-dimensional timeseries of observables into ranked dynamical modes. The key to the practical success of the approach is the identification of a set of observables that form a good basis on which to expand the slow relaxation modes. Good observables are, however, difficult to identify a priori and suboptimal choices can lead to significant underestimations of characteristic time scales. Leveraging the representation of slow dynamics in terms of Hidden Markov Models (HMM), we propose a simple and computationally efficient clustering procedure to infer surrogate observables that form a good basis for slow modes. Here, we apply the approach to an analytically solvable model system as well as on three protein systems of different complexities. We consistently demonstrate that the inferred indicator functions can significantly improve the estimation of the leading eigenvalues of Koopman operators and correctly identify key states and transition time scales of stochastic systems, even when good observables are not known a priori.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Aging phenomena in the two-dimensional complex Ginzburg-Landau equation

The complex Ginzburg-Landau equation with additive noise is a stochastic partial differential equation that describes a remarkably wide range of physical systems which include coupled non-linear oscillators subject to external noise near a Hopf bifurcation instability and spontaneous structure formation in non-equilibrium systems, e.g., in cyclically competing populations or oscillatory chemical reactions. Here, we employ a finite-difference method to numerically solve the noisy complex Ginzburg-Landau equation on a two-dimensional domain with the goal to investigate its non-equilibrium dynamics when the system is quenched into the “defocusing spiral quadrant”. We observe slow coarsening dynamics as oppositely charged topological defects annihilate each other, and characterize the ensuing aging scaling behavior. We conclude that the physical aging features in this system are governed by non-universal aging scaling exponents. We also investigate systems with control parameters residing in the “focusing quadrant”, and identify slow aging kinetics in that regime as well. Finally, we provide heuristic criteria for the existence of slow coarsening dynamics and physical aging behavior in the complex Ginzburg-Landau equation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Adaptive tau-leaping methods for microscopic-lattice kinetic Monte Carlo simulations

Traditional Kinetic Monte Carlo (KMC) approaches, rooted in Gillespie’s stochastic simulation algorithm, become computationally demanding in systems with a large range of timescales. The goal of this work is to propose and study new adaptive lattice-KMC time integration strategies for spatially non-uniform systems. To that end, two novel adaptive tau-leaping methods and their corresponding time integration strategies are developed based on the idea of the “n-fold” direct KMC method. These strategies allow for the simultaneous execution of multiple reactions, advancing time by adaptively selected coarse increments. We present numerical experiments comparing the proposed methods with existing approaches in a catalytic surface kinetics application involving ammonia decomposition.

Bimolecular reactions↗
DOE OSTI2026

Understanding and Controlling Photoexcited Molecules in Complex Environments (Final Technical Report)

Predicting the outcome of a photochemical reaction is complicated by factors outside the scope of tradition theoretical chemistry. On such small scales, fluctuations of the environment render outcomes stochastic, nonadiabatic dynamics blur the separation between nuclear and electronic motion, and ultrafast relaxation renders typical assumptions of equilibrium inappropriate. This project aimed to elucidate and control nonadiabatic molecular dynamics in condensed phases. We developed accurate, efficient simulation tools to study irreversible molecular processes computational across a range of scales. We applied these methodological advances to study electron transfer under exotic conditions where spin and topology opened new pathways for reactivity. We also brought these ideas to bare on light harvesting systems, where energy and charge transport are coupled, and where long lifetimes can be deleterious to function. These advances helped to establish general principles for chemical efficiency and guide the design of nanoscale energy materials, molecular machines, and related technologies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Structural Aspects of Neutron Survival Probabilities

The neutron survival probability (and related quantities including probabilities of extinction and initiation) is a central element of the broader stochastic theory of neutron populations and finds application in fields including reactor start-up, analysis of reactor power bursts and criticality accidents, and safeguards. In a full neutron transport formulation, the equation governing the single-neutron survival probability is a backward or adjoint-like integro-partial differential equation with the added complexity of being highly nonlinear. Analogous formulations of this equation exist in the context of many approximate theories of neutron transport, with the point kinetics formulation having received significant theoretical attention since the 1940s. This work continues this tradition by providing a novel analysis of the single-neutron survival probability equation using the tools of boundary layer theory. The analysis reveals that the “fully dynamic” solution of the single-neutron survival probability equation—and some key probability distributions derived from it—may be cast as a singular perturbation around the underlying quasi-static single-neutron probability of initiation. In this perturbation solution, the expansion parameter is the ratio of the neutron generation time to a macroscopic time scale characterizing the overall system evolution; this interpretation illuminates some of the fundamental structural aspects of neutron survival phenomena.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

A multi-stage stochastic programming model for adaptive biomass processing operation under uncertainty

Variations of physical and chemical characteristics of biomass reduce equipment utilization and increase operational costs of biomass processing. Biomass processing facilities use sensors to monitor the changes in biomass characteristics. Integrating sensory data into the operational decisions in biomass processing will increase its flexibility to the changing biomass conditions. In this paper, we propose a multi-stage stochastic programming model that minimizes the expected operational costs by identifying the initial inventory level and creating an operational decision policy for equipment speed settings. These policies take the sensory information data and the current biomass inventory level as inputs to dynamically adjust inventory levels and equipment settings according to the changes in the biomass' characteristics. We ensure that a prescribed target reactor utilization is consistently achieved by penalizing the violation of the target reactor feeding rate. A case study is developed using real-world data collected at Idaho National Laboratory's biomass processing facility. We show the value of multi-stage stochastic programming from an extensive computational experiment. Our sensitivity analysis indicates that updating the infeed rate of the system, the processing speed of equipment, and bale sequencing based on the moisture level of biomass improves the processing rate of the reactor and reduces operating costs.

09 BIOMASS FUELS↗
DOE OSTI2021

Uncertainty quantification for Multiphase-CFD simulations of bubbly flows: a machine learning-based Bayesian approach supported by high-resolution experiments

In this paper, we developed a machine learning-based Bayesian approach to inversely quantify and reduce the uncertainties of multiphase computational fluid dynamics (MCFD) simulations for bubbly flows. The proposed approach is supported by high-resolution two-phase flow measurements, including those by double-sensor conductivity probes, high-speed imaging, and particle image velocimetry. Local distributions of key physical quantities of interest (QoIs), including the void fraction and phasic velocities, are obtained to support the Bayesian inference. In the process, the epistemic uncertainties of the closure relations are inversely quantified while the aleatory uncertainties from stochastic fluctuations of the system are evaluated based on experimental uncertainty analysis. The combined uncertainties are then propagated through the MCFD solver to obtain uncertainties of the QoIs, based on which probability-boxes are constructed for validation. The proposed approach relies on three machine learning methods: feedforward neural networks and principal component analysis for surrogate modeling, and Gaussian processes for model form uncertainty modeling. The whole process is implemented within the framework of an open-source deep learning library PyTorch with graphics processing unit (GPU) acceleration, thus ensuring the efficiency of the computation. The results demonstrate that with the support of high-resolution data, the uncertainties of MCFD simulations can be significantly reduced. The proposed approach has the potential for other applications that involve numerical models with empirical parameters.

42 ENGINEERING↗
DOE OSTI2021

Transportation Hub Infrastructure Expansion: Decision Support Under Uncertainty

The Athena project (www.athena-mobility.org) has worked to investigate the relationship between the Dallas-Fort Worth Airport (DFW) and the greater Dallas area in order to better understand and therefore better inform future decision-making regarding the critical infrastructure that influence mobility between the airport and the city. Through this work, infrastructure related to curbside pickup and drop-off, parking, public transit, and the road network congestion were identified as critical to the operation of the DFW transportation hub. The infrastructure analysis and expansion aspect of the Athena project is focused on the restructuring of the CTA curb as a hierarchical curb and the building or repurposing of parking infrastructure as the interplay between these two areas. Many sources of uncertainty exist that may impact future airport and transportation hub operations, such as passenger volume growth, population demographic changes over time, electric vehicle (EV) adoption rates, and autonomous vehicle (AV) adoption rates. Due to these sources of uncertainty, we have selected for our research a modeling framework that can capture various types of uncertainty and hedge against those uncertainties in the optimization process. We analyze road network and curb congestion, the rise of transportation networking companies, trends in parking usage, existing policies around this infrastructure, airport revenue streams, and other contributing factors to enable infrastructure decision making with less uncertainty. To accomplish this wholistic analysis, we have developed a novel multi-stage, multi-period stochastic optimization model which considers the airport's decisions from 2025-2045 under different possible future macro trajectories and day-to-day variations in operational conditions captured as "annual representation of operations" scenarios with respective probabilities. This model has also been designed to leverage the outputs of various efforts under the Athena project to create a combined decision framework for infrastructure decisions. These various efforts include the route optimization model, the ASPIRES simulation, the mode choice model, and the SUMO traffic simulation. Our computational experiments of this system at scale have resulted in a working version of our infrastructure model which enables the explicit representation and consideration of various sources of uncertainty in the decision process to enable robust, flexible decision-making. This model has been effectively run on NREL's HPC system, Eagle, with large numbers of stochastic scenarios and shows promise as a scalable tool for robust consideration of uncertainties in airport planning. We have tested our model using 30,240 operational circumstances in total, resulting in a problem with more 200 million variables. This model was solved in several different configurations, and a workflow to simulate the performance of the infrastructure model results was developed and deployed. In general, our results indicate that a combination of remote parking, remote curb infrastructure, and dynamic pricing can generate revenue, reduce emissions, accommodate emerging technologies such as AVs and EVs, and manage airport passenger growth over time. We note the success of the proposed strategy depends on the data collection and forecasting abilities of DFW. We have also seen that the AV adoption by TNCs might necessitate larger amounts of remote curb. The results of this work inform strategies for airport infrastructure decision making, as well as demonstrate the value of an adaptable model, but also indicate that there are avenues remaining where further research would be of value.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Binding kinetics of harmonically confined random walkers

Diffusion-mediated binding of molecules under the influence of discrete spatially confining potentials is a commonly encountered scenario in systems subjected to explicit fields or implicit fields arising from tethering restraints. Here, in this work, we derive analytical expressions for the mean binding time of two random walkers geometrically confined by means of two harmonic potentials in one- and two-dimensional systems, which show excellent agreement with Brownian dynamics simulations. As a demonstration of its utility, we use this theory to maximize the communication speed in existing DNA walkers, obtaining quantitative agreement with previously reported experimental findings. The analytical expressions derived in this paper are broadly applicable to diverse systems, providing ways to characterize communication processes and optimize the rate of signal propagation for sensing and computing applications at the nanoscale.

36 MATERIALS SCIENCE↗