Search NASA⌕ Search

SEARCH · Search NASA

Results for “Random walk”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Devices and methods for increasing the speed and efficiency at which a computer is capable of modeling a plurality of random walkers using a particle method

A method for increasing a speed or energy efficiency at which a computer is capable of modeling a plurality of random walkers. The method includes defining a virtual space in which a plurality of virtual random walkers will move among different locations in the virtual space. The method also includes either assigning a corresponding set of ringed neurons in a spiking neural network to a corresponding virtual random walker, or assigning a corresponding set of ringed neurons to a point in the virtual space. Movement of a given virtual random walker is tracked by decoding differences between states of individual neurons in a corresponding given set of ringed neurons. A virtual random walk of the plurality of virtual random walkers is executed using the spiking neural network.

Aimone, James Bradley↗
DOE OSTI2022

Devices and methods for increasing the speed and efficiency at which a computer is capable of modeling a plurality of random walkers using a density method

A method for increasing a speed or energy efficiency at which a computer is capable of modeling a plurality of random walkers. The method includes defining a virtual space in which a plurality of virtual random walkers will move among different locations in the virtual space, wherein the virtual space comprises a plurality of vertices and wherein the different locations are ones of the plurality of vertices. A corresponding set of neurons in a spiking neural network is assigned to a corresponding vertex such that there is a correspondence between sets of neurons and the plurality of vertices, wherein a spiking neural network comprising a plurality of sets of spiking neurons is established. A virtual random walk of the plurality of virtual random walkers is executed using the spiking neural network, wherein executing includes tracking how many virtual random walkers are at each vertex at a given time increment.

Aimone, James Bradley↗
NASA NTRS1969

Solar protons E greater than 100 Mev incident over Antarctica during January- February 1967.

Commencing at 0825 +a -1 UT on January 28, 1967, a large and prolonged increase in the intensity of penetrating charged particles was observed by balloon-borne instruments floating over Byrd Station, Antarctica. (80°S, 120°W). A peak intensity of approximately 50 protons per cm 2 -sec-steradian with E> 100 MeV occurred at about 1230 UT on the 28th. The event was under observation almost continuously over a period of about 100 hours until the intensity decayed below cosmic-ray background on February 1. The initial decay was rapid but, some 40 hours after onset, went over into a slow exponential decay characterized by a 20 hour time-constant. The decay phase of an additional, though considerably less intense, event was observed on February 3 and 4. Presumably both events had their origins in major disturbances on the far side of the sun since nether event has been definitely linked to any feature which existed on the visible disk within an appropriate time interval. Results pertaining to the time-intensity profile and to the energy spectrum for protons E> 100 MeV are presented for the January 29 event. Comparison of the balloon results with neutron-monitor and satellite measurements and with models of interplanetary diffusion has led to some conclusions regarding the role of small-angle scattering by irregularities and by the random walk of magnetic lines of force relative to the mean interplanetary field within the orbit of earth.

Energy spectrum↗
DOE OSTI2019

Development of Upscaling Techniques and Construction of Calibrated Models for Fractured Rocks Using Discrete Fracture Network Approaches

The Pahute Mesa Corrective Action Unit on the Nevada National Security Site (NNSS) contains several fractured aquifers that can potentially provide high-permeability pathways for radionuclides to migrate away from underground nuclear testing locations. Although the geometric properties of fractures, such as length and orientation, can be found using geophysical methods and borehole image analyses, the hydraulic property of fractures (which are primarily influenced by apertures) are often unknown or have a high degree of uncertainty attached to their values. The aperture values have a profound impact on groundwater velocity and are critical to understanding transport scales at Pahute Mesa. Fracture aperture values and their distributions also have dominant control over radionuclide transport with potential for long-range correlations or scale-invariant dynamics. Understanding of the small-scale transport variability in fractures could lead to identification of distinct characteristics on larger scales. Furthermore, upscaling techniques need to be developed to adequately predict transport behavior at field scales by capturing the transport trends of smaller scale networks. This report explores these vital modeling needs to understand the role of fractures in Pahute Mesa. Discrete fracture network (DFN) models are the basis for constructing a well-calibrated model to ascertain aperture values of fractures in lava-flow aquifers. These models are also used to develop and compare the usefulness of several random-walk-based models to perform transport upscaling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Binary operations on neuromorphic hardware with application to linear algebraic operations and stochastic equations

Abstract Non-von Neumann computational hardware, based on neuron-inspired, non-linear elements connected via linear, weighted synapses—so-called neuromorphic systems—is a viable computational substrate. Since neuromorphic systems have been shown to use less power than CPUs for many applications, they are of potential use in autonomous systems such as robots, drones, and satellites, for which power resources are at a premium. The power used by neuromorphic systems is approximately proportional to the number of spiking events produced by neurons on-chip. However, typical information encoding on these chips is in the form of firing rates that unarily encode information. That is, the number of spikes generated by a neuron is meant to be proportional to an encoded value used in a computation or algorithm. Unary encoding is less efficient (produces more spikes) than binary encoding. For this reason, here we present neuromorphic computational mechanisms for implementing binary two’s complement operations. We use the mechanisms to construct a neuromorphic, binary matrix multiplication algorithm that may be used as a primitive for linear differential equation integration, deep networks, and other standard calculations. We also construct a random walk circuit and apply it in Brownian motion simulations. We study how both algorithms scale in circuit size and iteration time.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Odd Diffusivity of Chiral Random Motion

Diffusive transport is characterized by a diffusivity tensor which may, in general, contain both a symmetric and an antisymmetric component. Although the latter is often neglected, we derive Green-Kubo relations showing it to be a general characteristic of random motion breaking time-reversal and parity symmetries, as encountered in chiral active matter. In analogy with the odd viscosity appearing in chiral active fluids, we term this component the odd diffusivity. Furthermore, we show how odd diffusivity emerges in a chiral random walk model, and demonstrate the applicability of the Green-Kubo relations through molecular dynamics simulations of a passive tracer particle diffusing in a chiral active bath.

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

Semicoherent symmetric quantum processes: Theory and applications

Discovering pragmatic and efficient approaches to construct ε-approximations of quantum operators such as real (imaginary) time-evolution propagators in terms of the basic quantum operations (gates) is challenging. Prior ε-approximations are invaluable, in that they enable the compilation of classical and quantum algorithm modeling of, e.g., dynamical and thermodynamic quantum properties. In parallel, symmetries are powerful tools concisely describing the fundamental laws of nature; the symmetric underpinnings of physical laws have consistently provided profound insights and substantially increased predictive power. In this work, we consider the interplay between the ε-approximate processes and the exact symmetries in a semicoherent context—where measurements occur at each logical clock cycle. Here we draw inspiration from Pascual Jordan's groundbreaking formulation of nonassociative, but commutative, symmetric algebraic form. Our symmetrized formalism is then applied in various domains such as quantum random walks, real-time evolutions, variational algorithm ansatzes, and efficient entanglement verification. Our work paves the way for a deeper understanding and greater appreciation of how symmetries can be used to control quantum dynamics in settings where coherence is a limited resource.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Multiscale modeling of solute diffusion in triblock copolymer membranes

We develop a multiscale simulation model for diffusion of solutes through porous triblock copolymer membranes. The approach combines two techniques: self-consistent field theory (SCFT) to predict the structure of the self-assembled, solvated membrane and on-lattice kinetic Monte Carlo (kMC) simulations to model diffusion of solutes. Solvation is simulated in SCFT by constraining the glassy membrane matrix while relaxing the brush-like membrane pore coating against the solvent. The kMC simulations capture the resulting solute spatial distribution and concentration-dependent local diffusivity in the polymer-coated pores; we parameterize the latter using particle-based simulations. We apply our approach to simulate solute diffusion through nonequilibrium morphologies of a model triblock copolymer, and we correlate diffusivity with structural descriptors of the morphologies. We also compare the model’s predictions to alternative approaches based on simple lattice random walks and find our multiscale model to be more robust and systematic to parameterize. Furthermore, our multiscale modeling approach is general and can be readily extended in the future to other chemistries, morphologies, and models for the local solute diffusivity and interactions with the membrane.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Dynamical properties of a room temperature ionic liquid: Using molecular dynamics simulations to implement a dynamic ion cage model [plus supplemental information]

The transport behavior of ionic liquids (ILs) is pivotal for a variety of applications, especially when ILs are used as electrolytes. Many aspects of the transport dynamics of ILs remain to be understood. Here, a common ionic liquid, 1-butyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide (BmimNTf 2 ), was studied with molecular dynamics simulations. The results show that BmimNTf 2 displays typical structural relaxation, subdiffusive behavior, and a breakdown of the Stokes-Einstein diffusion relation as in glass-forming liquids. In addition, the simulations show that the translational dynamics, reorientation dynamics, and structural relaxation dynamics are well described by the Vogel-Fulcher-Tammann equation like fragile glass forming liquids. Building on previous work that employed ion cage models, it was found that the diffusion dynamics of the cations and anions were well described by a hopping process random walk where the step time is the ion cage lifetime obtained from the cage correlation function. Detailed analysis of the ion cage structures indicated that the electrostatic potential energy of the ion cage dominates the diffusion dynamics of the caged ion. The ion orientational relaxation dynamics showed that ion reorientation is a necessary step for ion cage restructuring. The dynamic ion cage model description of ion diffusion presented here may have implications for designing ILs to control their transport behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Transport Upscaling under Flow Heterogeneity and Matrix-Diffusion in Three-Dimensional Discrete Fracture Networks

For this work, we investigate the combined effects of network scale flow variability and retention due to matrix-diffusion on the scaling behavior of transport through fractured media. Two of the principal mechanisms controlling the transport of solutes through fractured low-permeability media are broad distributions of flow velocities and retention times in the solid matrix. We study the relative impact of these two processes under different initial conditions using a set of three-dimensional discrete fracture network simulations. We use these simulations to develop and calibrate an upscaled continuous time random walk (CTRW) approach for advective transport based on an Ornstein-Uhlenbeck model for the particle velocities that accounts for the fracture-matrix coupling using a compound Poisson process. This CTRW model can be conditioned on the initial solute distribution and allows to observe late-time scaling behavior at distances beyond what is feasible to observe using high-fidelity direct numerical simulations. We determine that the initial distribution of particles leads to marked differences in the persistent long-term scale behavior in the solute travel time distributions, even those undergoing retention due to matrix diffusion through implementation and analysis of the model.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Staircases of passive and active scalar concentration in cellular flow

This paper develops a unified model for staircase formation in both passive and active scalar systems, building upon prior numerical studies by offering new heuristic and physical insights. While prior studies primarily reported numerical results, they did not explore the underlying unifying physics that governs both types of scalar transport; this work addresses that gap by identifying shared mechanisms across both cases. Results of studies of passive and active scalar staircase formation in cellular flows are presented. Staircase formation in cellular flows occurs due to the interplay of fast mixing within cells and slow transport across the inter-cell boundary. The cell boundary emerges as a de facto transport barrier. Special attention is focused on the effects of cellular fluctuations and noise upon staircase structure. A forced, fluctuating vortex array model is used to drive the underlying flow structure. Cellular Peclet number and staircase profile curvature are identified as figures-of-merit to quantify the resiliency of layering. These are related to simple, multi-scatterer scalar random walk models. Results for Peclet number and curvature scaling with flow excitation are presented. We also study staircases of magnetic potential evolving in two-dimensional magnetohydrodynamics as examples of layering of active scalar concentration. Formation of magnetic potential staircases is indeed observed. Flux expulsion inhibits the intercellular transport of magnetic potential and strengthens staircase barriers. Magnetic staircases can be supported against resistive decay by magnetic potential noise forcing. Implications for staircase formation in magnetic confinement experiments are discussed.

Control theory↗
DOE OSTI2020

Markovian transport processes in a heterogeneous, variably saturated watershed: A multi-domain spatial Markov model

Transport processes in watersheds remain challenging to describe accurately with upscaled models because of their intricate complexities at multiple scales that can lead to multi-modal, non-Fickian breakthrough curves. One of the recent advances in solute transport modeling has been the Lagrangian spatial Markov model (SMM), which describes transport using velocity increments and models the correlated transitions between them. The method has been applied previously to saturated porous media and fractured media, but these successful applications suggest it may be useful for describing watershed scale processes as well. An existing, 3-D, heterogeneous variably saturated model of a headwaters catchment and classical random walk particle tracking were used to generate a forward simulation and the resulting transport was analyzed with a spatial Markov framework. Existing SMM methods were used to analyze the transport paths through the complex, unconfined domain and we found that the vadose zone plays an important role in the correlation structure. As such, a new kind of SMM is proposed, termed a multi-domain SMM (MD-SMM), that uses separate correlation models for auto-transitions in the vadose zone and the saturated zone, and for the cross-transitions between them. The MD-SMM improves representation of transport relative to the single domain model and can naturally delineate times when a particle is in the vadose zone or the saturated zone, adding an additional degree of realism to the upscaled model. Lastly, we show that the transition probability matrix for transport in the saturated portion of this unconfined system is approximately Markovian, which is an important validation of the proposed framework.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Geometric Delocalization in Two Dimensions

We demonstrate the existence of transient two-dimensional surfaces where a random-walking particle escapes to infinity in contrast to localization in standard flat two-dimensional space. We first prove that any rotationally symmetric two-dimensional membrane embedded in flat three-dimensional space cannot be transient. Then we formulate a criterion for the transience of a general asymmetric two-dimensional membrane. We use it to explicitly construct a class of transient two-dimensional manifolds with a nontrivial metric and height function but “zero average curvature,” which we dub “tablecloth manifolds.” The absence of the logarithmic infrared divergence of the Laplace-Beltrami operator in turn implies the absence of weak localization, nonexistence of bound states in shallow potentials, and breakdown of the Mermin-Wagner theorem and Kosterlitz-Thouless transition on the tablecloth manifolds, which may be realizable in both quantum simulators and corrugated two-dimensional materials.

Anderson localization↗
DOE OSTI2021

Impact of confined geometries on hopping and trapping of motile bacteria in porous media

We use a random walk particle-tracking (RWPT) approach to elucidate the impact of porous media confinement and cell-cell interactions on bacterial transport. The model employs stochastic alternating motility states consisting of hopping movement and trapping reorientation. The stochastic motility patterns are defined based on direct visualization of individual trajectory data. We validate our model against experimental data, at single-cell resolution, of bacterial E. coli motion in three-dimensional confined porous media. Results show that the model is able to efficiently simulate the spreading dynamics of motile bacteria as it captures the impact of cell-cell interaction and pore confinement, which marks the transition to a late-time subdiffusive regime. Furthermore, the model is able to qualitatively reproduce the observed directional persistence. Lastly, our RWPT model constitutes a meshless simple method which is easy to implement and does not invoke ad hoc assumptions but represents the basis for a multiscale approach to the study of bacterial dispersal in porous systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Metropolis–Hastings Monte Carlo Method for Neutron Emissivity Tomographic Inversion in Tokamak Plasma

The article presents a new approach to plasma neutron emissivity reconstruction based on the Metropolis–Hastings Monte Carlo algorithm. The algorithm is based on a biased random walk. A dedicated computer code generates pseudo-random samples within the domain that contains the solution. The properly chosen objective function ensures convergence to the desired solution. The method allows to overcome the problem of trapping in local minima by accepting, with some finite probability, changes in the direction opposite to that indicated by the objective function. The developed method was tested using a synthetic data set. The results show that the method provides reliable reconstructions and thus can be used as an alternative or complementary method to the commonly applied ones. The method capabilities were also tested against noise present in the line-integrated projections.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Diffusion in Lamellae, Cylinders, and Double Gyroid Block Copolymer Nanostructures

We study transport of penetrants through nanoscale morphologies motivated by common block copolymer morphologies, using confined random walks and coarse-grained simulations. Diffusion through randomly oriented grains is 1/3 for cylinder and 2/3 for lamellar morphologies versus an unconstrained (homopolymer) system, as previously understood. Diffusion in the double gyroid structure depends on volume fraction and is 0.47–0.55 through the minority phase at 30–50 vol.% and 0.73–0.80 through the majority at 50–70 vol.%. In conclusion, among randomly oriented standard minority phase structures with no grain boundary effects, lamellae is preferable for transport.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Efficient Sampling of Complex Interdependent and Multiplex Networks

Efficient sampling of interdependent and multiplex infrastructure networks is critical for effectively applying failure and recovery algorithms in real-world settings, as well as to generate property-preserving reduced-order graph-based ensembles that address topological uncertainties. In this paper, we first explore the performance, i.e. the success in preserving graph properties, of graph sampling algorithms for interdependent and multiplex networks with synthetic and real-world graphs. We simulate sampling algorithms under different parameter settings. These settings include probabilistic graph generators, coupling patterns, and various performance metrics. Our results show that while Random Node and Random Walk sampling algorithms perform best for interdependent networks, Random Edge and Forest Fire sampling algorithms perform best for multiplex networks. Second, we propose and implement a novel similarity-based sampling algorithm for multiplex networks that samples only log(N) number of layers of an N-layer multiplex network while yielding computational savings with performance guarantees. Experimental results show that similarity sampling outperforms complete sampling of all layers while decreasing performance costs from a linear scale to a logarithmic one. Our results also indicate that similarity-based sampling outperforms complete sampling and random selection in nearly all scenarios when tested with real-world data.

Subasi, Omer↗
NASA NTRS1972

Model for radon diffusion through the lunar regolith.

Description of a model for radon diffusion through the lunar regolith in which the atom migrates by random walk. The regolith is represented by a system of randomly oriented baffles in which the mean distance which the atom travels between two collisions takes on the role of a mean free path. The effective mean time between two collisions depends on two entities: the actual mean time-of-flight and the mean sticking time on grain surfaces for one collision. The latter depends strongly on the temperature and the heat of adsorption of radon on regolith materials. Both the mean free path as well as the heat of adsorption are either poorly known or unknown for the lunar regolith; hence these quantities are treated as free parameters. Because of the greatly different mean lifetimes against radioactive decay of Rn219, Rn220, and Rn222, the regolith acts as a powerful 'filter' for these species. Rn222 escape is significant (32%) even for a mean free path of 1 micron, a heat of adsorption of 7.0 kcal/mole and a regolith depth of 4 m. Calculations of radon escape from a 4 m thick regolith, using mean free paths of 1, 10, and 80 microns and heats of adsorption of 4.0, 5.2, and 7.0 kcal/mole show that the Rn222/Rn220 escape ratio can be as small as 7.7 and as large as, or larger than 47. The small value of 7.7 is of particular interest, because it is nearly equal to the escape ratio inferred by Turkevich et al. (1970) from their Surveyor 5 results.

Friesen, L. J.↗