SEARCH · Engineering Papers

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.

SourceResearch connectionFrom yearThrough year

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

At least 37 records · Page 2

All sources NASA NTRS DOE OSTI

NASA NTRS1990

Growth morphology with anisotropic surface kinetics

The morphological evolution of crystals growing from an incongruent vapor phase is studied using a Monte Carlo model, and the full range of growth morphologies is recovered. The diffusion in the bulk nutrient and the anisotropy in the interface kinetics are morphologically destabilizing and stabilizing, respectively. For a given set of simulation parameters and lattice symmetries there is a critical size, which scales linearly with the mean free path in the vapor, beyond which a crystal cannot retain its stable, macroscopically faceted growth shape. Surface diffusion stabilizes faceted growth on the shorter scale of the mean surface diffusion length. In simulations with a uniform drift superimposed on the random walk nutrient transport, crystal faces oriented toward the drift show enhanced morphological stability compared to the purely diffusive situation. Rotational drifts with periodic reversal of direction are morphologically stabilizing for all crystal facets.

Xiao, Rong-Fu↗

NASA NTRS1999

Dynamical Signatures of Living Systems

One of the main challenges in modeling living systems is to distinguish a random walk of physical origin (for instance, Brownian motions) from those of biological origin and that will constitute the starting point of the proposed approach. As conjectured, the biological random walk must be nonlinear. Indeed, any stochastic Markov process can be described by linear Fokker-Planck equation (or its discretized version), only that type of process has been observed in the inanimate world. However, all such processes always converge to a stable (ergodic or periodic) state, i.e., to the states of a lower complexity and high entropy. At the same time, the evolution of living systems directed toward a higher level of complexity if complexity is associated with a number of structural variations. The simplest way to mimic such a tendency is to incorporate a nonlinearity into the random walk; then the probability evolution will attain the features of diffusion equation: the formation and dissipation of shock waves initiated by small shallow wave disturbances. As a result, the evolution never "dies:" it produces new different configurations which are accompanied by an increase or decrease of entropy (the decrease takes place during formation of shock waves, the increase-during their dissipation). In other words, the evolution can be directed "against the second law of thermodynamics" by forming patterns outside of equilibrium in the probability space. Due to that, a specie is not locked up in a certain pattern of behavior: it still can perform a variety of motions, and only the statistics of these motions is constrained by this pattern. It should be emphasized that such a "twist" is based upon the concept of reflection, i.e., the existence of the self-image (adopted from psychology). The model consists of a generator of stochastic processes which represents the motor dynamics in the form of nonlinear random walks, and a simulator of the nonlinear version of the diffusion equation which represents the mental dynamics. It has been demonstrated that coupled mental-motor dynamics can simulate emerging self-organization, prey-predator games, collaboration and competition, "collective brain," etc.

Zak, M.↗

NASA NTRS1985

Simulations of the mean solar magnetic field during sunspot cycle 21

Regarding new bipolar magnetic regions as sources of flux, the evolution of the photospheric magnetic field during 1976-1984 was computed and the corresponding evolution of the mean line-of-sight field as seen from earth was derived. A good, but imperfect, agreement was obtained between the observed mean field and the field computed for a nominal choice of flux transport parameters. The response of the computed mean field to variations in the transport parameters and the source properties was determined. The results suggest that the mean-field evolution is a random-walk process with dissipation. New eruptions of flux produce the random walk, and together differential rotation, meridional flow, and diffusion provide the dissipation. The net effect of each new source depends on its strength and orientation and on the time elapsed before the next eruption.

Sheeley, N. R., Jr.↗

NASA NTRS1995

Large-scale magnetic variances near the South Solar Pole

We summarize recent Ulysses observations of the variances over large temporal scales in the interplanetary magnetic field components and their increase as Ulysses approached the South Solar Pole. A model of these fluctuations is shown to provide a very good fit to the observed amplitude and temporal variation of the fluctuations. In addition, the model predicts that the transport of cosmic rays in the heliosphere will be significantly altered by this level of fluctuations. In addition to altering the inward diffusion and drift access of cosmic rays over the solar poles, we find that the magnetic fluctuations also imply a large latitudinal diffusion, caused primarily by the associated field-line random walk.

Jokipii, J. R.↗

NASA NTRS1995

Field lines and magnetic surfaces in a two-component slab/2D model of interplanetary magnetic fluctuations

A two-component model for the spectrum of interplanetary magnetic fluctuations was proposed on the basis of ISEE observations, and has found an intriguing level of application in other solar wind studies. The model fluctuations consist of a fraction of 'slab' fluctuations, varying only in the direction parallel to the locally uniform mean magnetic field B(0) and a complement of 2D (two-dimensional) fluctuations that vary in the directions transverse to B(0). We have developed an spectral method computational algorithm for computing the magnetic flux surfaces (flux tubes) associated with the composite model, based upon a precise analogy with equations for ideal transport of a passive scalar in planar two dimensional geometry. Visualization of various composite models will be presented, including the 80 percent 2D/ 20 percent slab model with delta B/B(0) approximately equals 1 and a minus 5/3 spectral law, that is thought to approximately represent a snapshot of solar wind turbulence. Characteristically, the visualizations show that flux tubes, even when defined as regular on some plane, shred and disperse rapidly as they are viewed along the parallel direction. This diffusive process, which generalizes the standard picture of field line random walk, will be discussed in detail. Evidently, the traditional picture that flux tubes randomize like strands of spaghetti with a uniform tangle along the axial direction is in need of modification.

Matthaeus, W. H.↗

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.↗

NASA NTRS1985

Calculation of charge-state ratios for satellite Tor I

The diffusion of ions in a satellite plasma torus is presently modeled in terms of a one-dimensional random walk in which the particle source is at 0, the particle sink is at an N value that is an integer greater than 2, and the scale size of the diffusion cell is unity. The probability distribution function of the number of steps to exit for an ion is obtained and used in a model which incorporates ionization by electron impact to derive steady state expressions for the ratio of doubly to singly ionized ions, as well as the total number of ions in the torus. The results thus obtained are applied to the torus of the Jovian satellite Io, in order to predict mean residence times for sulfur and oxygen ions.

Summers, D.↗

NASA NTRS1983

A second-order effect of stratospheric vertical motions

A byproduct of stratospheric motion distinct from the normal downgradient transport represented by vertical eddy diffusion is studied. It is suggested that an air parcel executes a random walk over a range that could be several kilometers, and similar amplitudes could be present in wavelike oscillations. The photochemical loss rate should therefore be averaged over this range, which give an effect that is proportional to the curvature of the profile of loss rate, and is particularly large for fluorocarbon-11, CFCl3. Test computations are presented which indicate that the known discrepancy between observations and predictions for this compound can be at least partly reconciled for plausible values of the amplitude.

Hunten, D. M.↗

NASA NTRS1973

Magnetic fields and the distribution of cosmic rays in the Galaxy

If a galactic origin for the bulk of cosmic-ray particles is assumed, it is the magnetic field which regulates the propagation and escape of the particles from the Galaxy. Questions of the anisotropy of the cosmic-ray flux are considered together with magnetic-field models which can confine the particles and produce low anisotropies. In one standard approach it is postulated that the particles random walk through the presumably irregular magnetic field. In view of the difficulties associated with a diffusion hypothesis other confinement mechanisms have been proposed.

Jokipii, J. R.↗

NASA NTRS1986

Can stochastic, dissipative wave fields be treated as random walk generators

A suggestion by Meek et al. (1985) that the gravity wave field be viewed as stochastic, with significant nonlinearities, is applied to calculate diffusivities. The purpose here is to calculate the diffusivity for stochastic wave model and compare it with previous diffusivity estimates. The researchers do this for an idealized case in which the wind velocity changes but slowly, and for which saturation is the principal mechanism by which wave energy is lost. A related calculation was given in a very brief way (Weinstock, 1976), but the approximations were not fully justified, nor were the physical pre-suppositions clearly explained. The observations of Meek et al. (1985) have clarified the pre-suppositions for the researchers and provided a rationalization and improvement of the approximations employed.

Weinstock, J.↗

NASA NTRS1986

The diffusion of individual molecules within a gas

The Direct Simulation Monte Carlo method is used to study the positional history of the individual molecules in a gas that is homogeneous at the macroscopic level and is in Maxwellian equilibrium at the microscopic level. The behavior at small times is characterized by 'persistence of velocity' effects, and a 'random walk' type of dispersal occurs over a longer timescale. It is shown that the rate of dispersal can be directly related to the self-diffusion coefficient. In addition, the diffusion coefficients are obtained directly from one-dimensional calculations, and the local Knudsen number at which the Chapman-Enskog theory breaks down is determined. Results are presented for both simple gases and gas mixtures.

Bird, G. A.↗

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↗

NASA NTRS1987

Numerical methods for one-dimensional reaction-diffusion equations arising in combustion theory

A review of numerical methods for one-dimensional reaction-diffusion equations arising in combustion theory is presented. The methods reviewed include explicit, implicit, quasi-linearization, time linearization, operator-splitting, random walk and finite-element techniques and methods of lines. Adaptive and nonadaptive procedures are also reviewed. These techniques are applied first to solve two model problems which have exact traveling wave solutions with which the numerical results can be compared. This comparison is performed in terms of both the wave profile and computed wave speed. It is shown that the computed wave speed is not a good indicator of the accuracy of a particular method. A fourth-order time-linearized, Hermitian compact operator technique is found to be the most accurate method for a variety of time and space sizes.

Ramos, J. I.↗

NASA NTRS1987

Gyroresonant pitch angle scattering by coherent and incoherent whistler mode waves in the magnetosphere

A test particle approach is used to compare gyroresonant pitch angle scattering of energetic electrons by coherent versus incoherent whistler mode waves, for the case in which the coherent wave amplitude is below the nonlinear phase trapping threshold. Wave packets of 400 ms duration propagating along the magnetic field at L = 4 within the plasmasphere are considered, and the wave-induced pitch angle scattering along the propagation path from one hemisphere to the other and the resulting precipitation flux are computed. An incoherent wave spectrum is simulated by random modulation of the wave frequency at intervals of 1 ms, thereby generating signals with nearly constant power spectral density over a bandwidth of 2 kHz centered at 5.5 kHz. The associated pitch angle scattering is compared with that of a monochromatic 5.5-kHz signal of 400 ms duration. Results of the test particle analysis are compared with those expected on the basis of a classical diffusion treatment, and an expression is derived for an effective “diffusion” coefficient for pitch angle scattering by coherent waves. The trajectory followed by a particle when interacting with incoherent waves essentially represents a random walk in velocity space, while for coherent waves the pitch angle of the particle varies in a well-defined manner. In spite of the fact that individual particle scatterings are typically larger for coherent waves, the peak precipitation fluxes induced by incoherent waves are found to be approximately the same as those for coherent waves having the same total power. This results from the fact that incoherent waves interact with particles over a wider range of energies. As a consequence, the energy spectrum and the temporal extent of transient precipitation pulses due to incoherent wave packets are broader than those for equivalent coherent ones.

Umran S Inan↗

NASA NTRS1993

A diffusion-limited aggregation model for the evolution of drainage networks

We propose a modified diffusion-limited aggregation (DLA) model for the evolution of fluvial drainage networks. Random walkers are introduced randomly on a grid, and each two-dimensional random walk proceeds until the walker finds a drainage network on which to accrete. This model for headward growth of drainage networks generates drainage patterns remarkably similar to actual drainages. The model also predicts statistical features which agree with actual networks, including the frequency-order (bifurcation) ratio (R(sub b) = 3.98) and the stream length-order (R(sub r) = 2.09). Using the definition of network fractal dimension D = log R(sub b)/log R(sub r), we find that our DLA model gives D = 1.87, near the observed range of D approximately equal to 1.80 - 1.85.

Masek, Jeffrey G.↗

NASA NTRS1982

Random element method for numerical modeling of diffusional processes

The random element method is a generalization of the random vortex method that was developed for the numerical modeling of momentum transport processes as expressed in terms of the Navier-Stokes equations. The method is based on the concept that random walk, as exemplified by Brownian motion, is the stochastic manifestation of diffusional processes. The algorithm based on this method is grid-free and does not require the diffusion equation to be discritized over a mesh, it is thus devoid of numerical diffusion associated with finite difference methods. Moreover, the algorithm is self-adaptive in space and explicit in time, resulting in an improved numerical resolution of gradients as well as a simple and efficient computational procedure. The method is applied here to an assortment of problems of diffusion of momentum and energy in one-dimension as well as heat conduction in two-dimensions in order to assess its validity and accuracy. The numerical solutions obtained are found to be in good agreement with exact solution except for a statistical error introduced by using a finite number of elements, the error can be reduced by increasing the number of elements or by using ensemble averaging over a number of solutions.

Ghoniem, A. F.↗

NASA NTRS2011

Computing Temperatures in Optically Thick Protoplanetary Disks

We worked with a Monte Carlo radiative transfer code to simulate the transfer of energy through protoplanetary disks, where planet formation occurs. The code tracks photons from the star into the disk, through scattering, absorption and re-emission, until they escape to infinity. High optical depths in the disk interior dominate the computation time because it takes the photon packet many interactions to get out of the region. High optical depths also receive few photons and therefore do not have well-estimated temperatures. We applied a modified random walk (MRW) approximation for treating high optical depths and to speed up the Monte Carlo calculations. The MRW is implemented by calculating the average number of interactions the photon packet will undergo in diffusing within a single cell of the spatial grid and then updating the packet position, packet frequencies, and local radiation absorption rate appropriately. The MRW approximation was then tested for accuracy and speed compared to the original code. We determined that MRW provides accurate answers to Monte Carlo Radiative transfer simulations. The speed gained from using MRW is shown to be proportional to the disk mass.

radiative transfer↗

NASA NTRS1989

The effect of an interaction of magnetic flux and supergranulation on the decay of magnetic plages

This paper studies how the properties of large-scale convection affect the decay of plages. The plage decay, caused by the random-walk dispersion of flux tubes, is suggested to be severely affected by differences between the mean size of cellular openings within and around plages. The smaller cell size within a plage largely explains the smaller diffusion coefficient within plages as compared to that of the surrounding regions. The semipermeability of the plage periphery, together with the dependence of the diffusion coefficient on the flux-tube density, can explain the observed slow decay of plages (predicting a typical life time of about a month for a medium-sized plage), the existence of a well-defined plage periphery, and the observed characteristic mean magnetic flux density of about 100 G. One effect of the slowed decay of the plage by the semipermeability of the plage periphery is the increase of the fraction of the magnetic flux that can cancel with flux of the opposite polarity along the neutral line to as much as 80 percent as compared to at most 50 percent in the case of nonuniform diffusion. This may explain why only a small fraction of the magnetic flux is observed to escape from the plage into the surrounding network.

Schrijver, C. J.↗