Search NASA⌕ Search

SEARCH · Search NASA

Results for “Diffusion & random walks”

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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

Grid-free simulation of diffusion using random walk methods

The simulation of the diffusion of a continuum field by the random walk (RW) displacement of a set of particles is considered. Elements of the gradients of the diffusive concentration are transported by computational particles. It is demonstrated that, by the use of concentration gradients in the RW process, statistical errors are reduced and each realization of the numerical solution is a representation of the exact solution. The algorithm is grid-free, and the computational elements move to follow the gradients; hence, the algorithm is self-adaptive, and uniform resolution is achieved for all times.

Ghoniem, A. F.↗
NASA NTRS1973

Monte-Carlo analysis of rarefied-gas diffusion including variance reduction using the theory of Markov random walks

Molecular diffusion through a rarefied gas is analyzed by using the theory of Markov random walks. The Markov walk is simulated on the computer by using random numbers to find the new states from the appropriate transition probabilities. As the sample molecule during its random walk passes a scoring position, which is a location at which the macroscopic diffusing flow variables such as molecular flux and molecular density are desired, an appropriate payoff is scored. The payoff is a function of the sample molecule velocity. For example, in obtaining the molecular flux across a scoring position, the random walk payoff is the net number of times the scoring position has been crossed in the positive direction. Similarly, when the molecular density is required, the payoff is the sum of the inverse velocity of the sample molecule passing the scoring position. The macroscopic diffusing flow variables are then found from the expected payoff of the random walks.

Perlmutter, M.↗
NASA NTRS1976

Low pressure radon diffusion - A laboratory study and its implications for lunar venting

Results of a study of radon migration through columns of fine particulate materials, at total pressures of 0.02-0.2 torr, are reported. Materials studied were: NBS Glass Spheres (SRM 1003), Emerson & Cuming Eccospheres (IG-101), activated coconut charcoal, Lipaci obsidian, and W-1 Standard Diabase. Rates of diffusion were used to derive heats of adsorption for radon on the materials tested. The most reliable values found clustered around 8-9 kcal/mole. These high heats of adsorption, if typical for most materials, combined with low percentages of radon emanation by lunar soils found by other researchers, imply that random walk diffusion will not be an important mechanism for redistributing the radon and the radon daughters produced in the lunar regolith. In particular, since random walk migration is not a sufficient mechanism to account for localized high concentrations of radon-222 and its daughter polonium-210 observed by the Apollo 15 and 16 command modules, an alternative mechanism is proposed, in which radon would be swept to the surface by other gases during intermittent venting events.

Friesen, L. J.↗
NASA NTRS1991

Basic research in solar physics

This grant, dating back more than 20 years has supported a variety of investigations of the chromospheres and coronae of the Sun and related cool stars by the Principal Investigator, his postdocs and graduate students, and colleagues at other institutions. This work involved studies of radiative transfer and spectral line formation theory, and the application of these techniques to the analysis of spectra obtained from space and ground-based observatories in the optical, ultraviolet, x-ray and radio portions of the spectrum. Space observations have included the analysis of spectra from OSO-7, Skylab, SMM, and the HRTS rocket experiments. Recent work has concentrated on the interaction of magnetic fields, plasma and radiation in the outer atmospheres of the Sun and other magnetically active stars with different fundamental parameters. Our study of phenomena common to the Sun and stars, the 'solar-stellar connection', can elucidate the fundamental physics, because spatially-resolved observations of the Sun provide us with the 'groundtruth,' while interpretation of stellar data permit us to isolate those parameters critical to stellar activity. Recently, we have studied the differences in physical properties between solar regions of high magnetic flux density and the surrounding plasma. High-resolution CN and CO spectroheliograms have been used to model the thermal inhomogeneities driven by unstable CO cooling, and we have analyzed spatially resolved UV spectra from HRTS to model the thermal structure and energy balance of small-scale structures. The study of nonlinear relations between atmospheric radiative losses and the photospheric magnetic flux density has been continued. We have also proposed a new model for the decay of plages by random walk diffusion of magnetic flux. Our analysis of phenomena common to the Sun and stars included the application of available spectroscopic diagnostics, establishing evidence that the atmospheres of the least active stars are heated at a 'basal' rate that is also found in the centers of solar supergranules, and using the Doppler-imaging technique to measure the position, size, and brightness of stellar active regions. We are computing multi-component models for solar and stellar atmospheres, and models for coronal loops and for the transition-region down flows. The study of solar and stellar flares permits us to assess the role of turbulent energy transport, to pinpoint the mechanism behind Type I radio bursts, to determine whether plasma radiation or cyclotron maser is responsible for microwave flares on M dwarfs, and to extend our knowledge of the basic physics pertinent to cyclotron-maser processes operating on the Sun.

Linsky, Jeffrey L.↗
NASA NTRS2014

Tuning Monotonic Basin Hopping: Improving the Efficiency of Stochastic Search as Applied to Low-Thrust Trajectory Optimization

Trajectory optimization methods using monotonic basin hopping (MBH) have become well developed during the past decade [1, 2, 3, 4, 5, 6]. An essential component of MBH is a controlled random search through the multi-dimensional space of possible solutions. Historically, the randomness has been generated by drawing random variable (RV)s from a uniform probability distribution. Here, we investigate the generating the randomness by drawing the RVs from Cauchy and Pareto distributions, chosen because of their characteristic long tails. We demonstrate that using Cauchy distributions (as first suggested by J. Englander [3, 6]) significantly improves monotonic basin hopping (MBH) performance, and that Pareto distributions provide even greater improvements. Improved performance is defined in terms of efficiency and robustness. Efficiency is finding better solutions in less time. Robustness is efficiency that is undiminished by (a) the boundary conditions and internal constraints of the optimization problem being solved, and (b) by variations in the parameters of the probability distribution. Robustness is important for achieving performance improvements that are not problem specific. In this work we show that the performance improvements are the result of how these long-tailed distributions enable MBH to search the solution space faster and more thoroughly. In developing this explanation, we use the concepts of sub-diffusive, normally-diffusive, and super-diffusive random walks (RWs) originally developed in the field of statistical physics.

autonomous↗
NASA NTRS1994

On the effects of a weak magnetic field on turbulent transport

We discuss the effects of a weak large-scale magnetic field on turbulent transport. We show by means of a series of two-dimensional numerical experiments that turbulent diffusion can be effectively suppressed by a (large scale) magnetic field whose energy is small compared to equipartition. The suppression mechanism is associated with a subtle modification of the Lagrangian energy spectrum, and it does not require any substantial reduction of the turbulent amplitude. We exploit the relation between diffusion and random walking to emphasize that the effect of a large-scale magnetic field is to induce a long-term memory in the field of turbulence. The implications of the general case of three-dimensional transport are briefly discussed.

Cattaneo, F.↗
NASA NTRS2023

Does the way we do science foster discovery?

Freedom to explore the unknown is key to scientific discovery. Maximizing modern individualistic measures of scientific productivity like citations and number of publications may impede the progress of science as a whole.

discovery↗
NASA NTRS2000

Inertial Effects in Suspension Dynamics

The present work analyses the dynamics of a suspension of heavy particles in shear flow. The magnitude of the particle inertia is given by the Stokes number St = m(gamma/6(pi)a, which is the ratio of the viscous relaxation time of a particle tau(sub p) = m=6pi(eta)a to the flow time gamma(sup -1). Here, m is the mass of the particle, a is its size, eta is the viscosity of the suspending fluid and gamma is the shear rate. The ratio of the Stokes number to the Reynolds number, Re = (rho)f(gamma)a(exp 2)/eta, is the density ratio rho(sub p)/rho(sub f). Of interest is to understand the separate roles of particle (St) and fluid (Re) inertia in the dynamics of suspensions. In this study we focus on heavy particles, rho(sub p)/rho(sub f) much greater than 1, for which the Stokes number is finite, but the Reynolds number is sufficiently small for inertial forces in the fluid to be neglected; thus, the fluid motion is governed by the Stokes equations. On the other hand, the probability density governing the statistics of the suspended particles satisfies a Fokker-Planck equation that accounts for both configuration and momentum coordinates, the latter being essential for finite St. The solution of the Fokker-Planck equation is obtained to O(St) via a Chapman-Enskog type-procedure, and the conditional velocity distribution so obtained is used to derive a configuration-space Smoluchowski equation with inertial corrections. The inertial effects are responsible for asymmetry in the relative trajectories of two spheres in shear flow, in contrast to the well known symmetric structure in the absence of inertia. Finite St open trajectories in the plane of shear suffer a downward lateral displacement resulting from the inability of a particle of finite mass to follow the curvature of the zero-Stokes-number pathlines. In addition to the induced asymmetry, the O(St) inertial perturbation dramatically alters the nature of the near-field trajectories. The stable closed orbits (for St = 0) in the plane of shear now spiral in, approaching particle-particle contact in the limit. All trajectories starting from an initial offset of O(St(sup 1/2) or less (which remain open for St = 0) also spiral in. The asymmetry of the trajectories leads to a non-Newtonian rheology and diffusive behavior. The latter because a given particle (moving along a finite St open trajectory) suffers a net displacement in the transverse direction after a single interaction. A sequence of such uncorrelated displacements leads to the particle executing a random walk. The inertial diffusivity tensor is anisotropic on account of differing strengths of interaction in the gradient and vorticity directions. Since the entire region (constituting an in finite area) of closed orbits in the plane of shear spirals onto contact for #finite St, the latter represents a singular surface for the pair-distribution function. The exact form of the pair-distribution function at contact is still, however, indeterminate in the absence of non-hydrodynamic effects. It should also be noted that finite St non-rectilinear flows do not support a spatially uniform number density owing to the cross-streamline inertial migration of particles.

J. F. Brady↗
NASA NTRS1974

The abundances of nuclei in the cosmic radiation

The relative abundances are treated as a consequence of processes in cosmic ray transport occurring during passage of the radiation through interstellar material at high velocity. Some of the subjects mentioned are nuclear fragmentation and the production of secondary nuclei, nuclear reactions, energy loss and nuclear decay, ionization, the range-energy relation and propagation variables, capture and loss of electrons, the propagation of nuclei, the transport equation, equilibrium solutions, energy-dependent path length distribution, exponential path length distributions, discrete spectra, sources, supernovae, and the origin of the abundances. The connection between the space-time features of the sources, the material traversed, and the effects of magnetic fields is established by describing the particle-field interaction as a diffusive or random-walk process.

Reames, D. V.↗
NASA NTRS1983

Numerical solution of the problem of flame propagation by the use of the random element method

A numerical, grid-free algorithm is presented for one-dimensional reaction-diffusion model of laminar flame propagation in premixed gases. It is based on the random element method we developed for the analysis of diffusional processes. The effect of combustion is taken into account by applying the principle of fractional steps to separate the process of diffusion, modeled by the random walk of computational elements, from the exothermic effects of chemical reaction, monitoring their strength. The validity of the algorithm is demonstrated by application to flame propagation problems for which exact solutions exist. The flame speed evaluated by its use oscillates around the exact value at a relatively small amplitude, while the temperature and species concentration profiles are self-correcting in their convergence to the exact solution. A satisfactory resolution is obtained by the use of quite a small number of computational elements which automatically adjust their distribution of fit sharp gradients.

Ghoniem, A. F.↗
NASA NTRS1986

Fluctuation and rotation in diffusion-influenced monomolecular reactions

Models of partially absorbing boundaries for calculating diffusion-influenced chemical reactions are examined analytically. The limitations of diffusion-equation and random-walk models are indicated; the presence of an ambiguous length parameter in a standard radiation-boundary-condition (RBC) model is demonstrated; the RBC model for a boundary with oscillation (at a rate much faster than the rate of thermal collisions) between absorbing and nonabsorbing states is considered; an alternative stochastically switching absorber (SSA) model is introduced and shown to be physically consistent; the equivalence of one-dimensional RBC and SSA models under some conditions is shown; one time passage is calculated for an axially symmetric three-dimensional SSA model with anisotropic absorption; and the analogy between rotation and stochastic switching is explored.

Bashford, Donald↗
NASA NTRS2015

Analysis of Fibrous Felts for Flexible Ablators Using Synchrotron Hard X-Ray Micro-Tomography

We analyzed the material properties of low-density felts that are used as substrates for new-generation flexible and conformal carbon/phenolic ablators, and compared them with those of a rigid carbon fiber preform that is used to manufacture rigid carbon/phenolic ablators. Micro-tomography measurements were obtained using synchrotron X-rays, allowing the characterization of the materials microstructure at the scale of the fibers. Using the tomography voxels as computational grids, we computed tortuosity and room temperature conductivity. In addition we performed micro-scale simulations of the oxidation of carbon fibers using a random walk model for oxygen diffusion and a probability law to model surface reactions.

Ablation↗
NASA NTRS1994

The rotation of photospheric magnetic fields: A random walk transport model

In an earlier study of solar differential rotation, we showed that the transport of magnetic flux across latitudes acts to establish quasi-stationary patterns, therby accounting for the observed rigid rotation of the large-scale photospheric field. In that paper, the effect of supergranular convection was represented by a continuum diffusion, limiting the applicability of the calculations to large spatial scales. Here we extend the model to scales comparable to that of the supergranulation itself by replacing the diffusive transport with a discrete random walk process. Rotation curves are derived by cross-correlating the simulated photospheric field maps for a variety of time lags and spatial resolutions. When the lag between maps is relatively short less than or approximately = 15 days), the midlatitude correlation functions show two distinct components: a broad feature associated with the large-scale unipolar patterns and a narrow feature originating from small magnetic structures encompossing from one to several supergranular cells. By fitting the broad component we obtain the rigid rotation profile of the patterns, whereas by fitting the narrow component, we recover the differential rate of the photospheric plasma itself. For time lags of 1 month or greater, only the broad feature associated with the long-lived patterns remains clearly identifiable in the simulations.

Wang, Y. -M.↗
NASA NTRS1995

Putting quasilinear theory to the test: A comparison of analytical theory with results from a trajectory code

Parallel and perpendicular diffusion coefficients were computed numerically by following particle orbits in a simulated magnetic field. The simulated field was chosen to have delta B/B(sub o) small, so as to provide a test of quasilinear theory in a regime where the theory should be most accurate. The simulation space is large enough to contain many magnetic field correlation lengths, so that effects of field line random walk can be studied. After presenting results for parallel diffusion, we will focus on two controversial issues relating to perpendicular diffusion: (1) Do quasilinear descriptions of perpendicular diffusion retain any validity for particles whose Larmor radius is smaller than a correlation length? (2) Does field line random walk lead to particle diffusion in the usual sense, or does it produce 'compound' diffusion for which particles spread out proportionally to t(exp 1/4) instead of t(exp 1/2)?

Bieber, J. W.↗
NASA NTRS2014

Anomalous Transport of High Energy Cosmic Rays in Galactic Superbubbles

High-energy cosmic rays may exhibit anomalous transport as they traverse and are accelerated by a collection of supernovae explosions in a galactic superbubble. Signatures of this anomalous transport can show up in the particles' evolution and their spectra. In a continuous-time-random- walk (CTRW) model assuming standard diffusive shock acceleration theory (DSA) for each shock encounter, and where the superbubble (an OB stars association) is idealized as a heterogeneous region of particle sources and sinks, acceleration and transport in the superbubble can be shown to be sub-diffusive. While the sub-diffusive transport can be attributed to the stochastic nature of the acceleration time according to DSA theory, the spectral break appears to be an artifact of transport in a finite medium. These CTRW simulations point to a new and intriguing phenomenon associated with the statistical nature of collective acceleration of high energy cosmic rays in galactic superbubbles.

Barghouty, Nasser F.↗
NASA NTRS1972

Variance reduction in Monte Carlo analysis of rarefied gas diffusion

The present analysis uses the Monte Carlo method to solve the problem of rarefied diffusion between parallel walls. The diffusing molecules are evaporated or emitted from one of two parallel walls and diffused through another molecular species. The analysis treats the diffusing molecule as undergoing a Markov random walk and the local macroscopic properties are found as the expected value of the random variable, the random walk payoff. By biasing the transition probabilities and changing the collision payoffs the expected Markov walk payoff is retained but its variance is reduced so that the M. C. result has a much smaller error.

Perlmutter, M.↗
NASA NTRS1972

Variance reduction in Monte Carlo analysis of rarefied gas diffusion.

The problem of rarefied diffusion between parallel walls is solved using the Monte Carlo method. The diffusing molecules are evaporated or emitted from one of the two parallel walls and diffuse through another molecular species. The Monte Carlo analysis treats the diffusing molecule as undergoing a Markov random walk, and the local macroscopic properties are found as the expected value of the random variable, the random walk payoff. By biasing the transition probabilities and changing the collision payoffs, the expected Markov walk payoff is retained but its variance is reduced so that the Monte Carlo result has a much smaller error.

Perlmutter, M.↗