Engineering Papers⌕ Search

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.

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

70 records · Page 4
All sourcesNASA NTRSDOE OSTI
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.↗
DOE OSTI2024

Unified understanding of the impact of semiflexibility, concentration, and molecular weight on macromolecular-scale ring diffusion

Conformationally fluctuating, globally compact macromolecules such as polymeric rings, single-chain nanoparticles, microgels, and many-arm stars display complex dynamic behaviors due to their rich topological structure and intermolecular organization. Synthetic rings are hybrid objects with conformations that display both ideal random walk and compact globular features, which can serve as models of genomic DNA. To date, emphasis has been placed on the effect of ring molecular weight on their unusual behaviors. Here, we combine simulations and a microscopic force-level theory to build a unified understanding for how key aspects of ring dynamics depend on different tunable molecular properties including backbone rigidity, monomer concentration, degree of traditional entanglement, and molecular weight. Our large-scale molecular dynamics simulations of ring melts with very different backbone stiffnesses reveal unanticipated behaviors which agree well with our generalized theory. This includes a universal master curve for center-of-mass diffusion constants as a function of molecular weight scaled by a chemistry and thermodynamic state-dependent critical molecular weight that generalizes the concept of an entanglement cross-over for linear chains. The key physics is how backbone rigidity and monomer concentration induced changes of the entanglement length, interring packing, degree of interpenetration, and liquid compressibility slow down space-time dynamic-force correlations on macromolecular scales. A power law decay of the center-of-mass diffusion constant with inverse molecular weight squared is the first consequence, followed by an ultraslow activated hopping transport regime. Our results set the stage to address slow dynamics and kinetic arrest in different families of compact synthetic and biological polymeric systems.

Science & Technology - Other Topics↗
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.↗
DOE OSTI2024

The persistence of memory in ionic conduction probed by nonlinear optics

Predicting practical rates of transport in condensed phases enables the rational design of materials, devices and processes. This is especially critical to developing low-carbon energy technologies such as rechargeable batteries. For ionic conduction, the collective mechanisms, variation of conductivity with timescales and confinement, and ambiguity in the phononic origin of translation, call for a direct probe of the fundamental steps of ionic diffusion: ion hops. However, such hops are rare-event large-amplitude translations, and are challenging to excite and detect. Here we use single-cycle terahertz pumps to impulsively trigger ionic hopping in battery solid electrolytes. This is visualized by an induced transient birefringence, enabling direct probing of anisotropy in ionic hopping on the picosecond timescale. The relaxation of the transient signal measures the decay of orientational memory, and the production of entropy in diffusion. We extend experimental results using in silico transient birefringence to identify vibrational attempt frequencies for ion hopping. Using nonlinear optical methods, we probe ion transport at its fastest limit, distinguish correlated conduction mechanisms from a true random walk at the atomic scale, and demonstrate the connection between activated transport and the thermodynamics of information.

25 ENERGY STORAGE↗
NASA NTRS2017

Tortuosity Computations of Porous Materials using the Direct Simulation Monte Carlo

Low-density carbon fiber preforms, used as thermal protection systems (TPS) materials for planetary entry systems, have permeable, highly porous microstructures consisting of interlaced fibers. Internal gas transport in TPS is important in modeling the penetration of hot boundary-layer gases and the in-depth transport of pyrolysis and ablation products. The gas effective diffusion coefficient of a porous material must be known before the gas transport can be modeled in material response solvers; however, there are very little available data for rigid fibrous insulators used in heritage TPS.The tortuosity factor, which reflects the efficiency of the percolation paths, can be computed from the effective diffusion coefficient of a gas inside a porous material and is based on the micro-structure of the material. It is well known, that the tortuosity factor is a strong function of the Knudsen number. Due to the small characteristic scales of porous media used in TPS applications (typical pore size of the order of 50 micron), the transport of gases can occur in the rarefied and transitional regimes, at Knudsen numbers above 1. A proper way to model the gas dynamics at these conditions consists in solving the Boltzmann equation using particle-based methods that account for movement and collisions of atoms and molecules.In this work we adopt, for the first time, the Direct Simulation Monte Carlo (DSMC) method to compute the tortuosity factor of fibrous media in the rarefied regime. To enable realistic simulations of the actual transport of gases in the porous medium, digitized computational grids are obtained from X-ray micro-tomography imaging of real TPS materials. The SPARTA DSMC solver is used for simulations. Effective diffusion coefficients and tortuosity factors are obtained by computing the mean-square displacement of diffusing particles.We first apply the method to compute the tortuosity factors as a function of the Knudsen number for computationally designed materials such as random cylindrical fibers and packed bed of spheres with prescribed porosity. Results are compared to literature values obtained using random walk methods in the rarefied and transitional regime and a finite-volume method for the continuum regime. We then compute tortuosity factors for a real carbon fiber material with a transverse isotropic structure (FiberForm), quantifying differences between through-thickness and in-plain tortuosities at various Knudsen regimes.

Tortuosity↗
NASA NTRS1989

Computations of turbulent evaporating sprays

A computational study of turbulent evaporating sprays is reported. The major focus is to examine the sensitivity of the vaporization behavior of turbulent sprays to the transient liquid-phase processes. Three models considered to represent these processes are the thin skin, infinite diffusion, and diffusion limit models. Favre-averaged equations with k-epsilon-g turbulence model are employed for the gas phase. The Lagrangian approach with a stochastic separated flow method is used for the liquid phase where the effects of gas turbulence on droplet trajectories and interphase transport rates are considered using random-walk computations. Also the variable-property effects are considered in detail. Results indicate that, depending upon the boiling temperature and heat of vaporization of the fuel considered, the vaporization behavior of turbulent sprays may be quite sensitive to the modeling of transient liquid-phase processes. Thus, it is important that for most hydrocarbon fuels these processes be adequately represented in any comprehensive spray computations. The present results also provide further support to the conclusions of earlier studies which have been based on simplified spray configurations.

Aggarwal, S. K.↗
NASA NTRS1995

On the upstream boundary of electron foreshocks in the solar wind

The upstream boundary of electron foreshocks is defined as the path of the fastest electrons reflected by collisionless shocks and moving along the magnetic field in the solar wind. Considerable levels of magnetic fluctuations are found in these regions of the solar wind, and their effect is to create both a broadening and a fine structure of the electron foreshock boundary. The magnetic structure is studied by means of a 3-D numerical simulation of a turbulent magnetic field. Enhanced, anomalous diffusion is found, (Delta x(exp 2)) varies as s(sup alpha), where alpha is greater than 1 for typical values of the parameters (here, Delta x(exp 2) is the mean square width of the tangent magnetic surface and s is the field line length). This corresponds to a Levy flight regime for the magnetic field line random walk, and allows very efficient electron propagation perpendicular to the magnetic field. Implications on the observations of planetary foreshocks and of the termination shock foreshock are considered.

Zimbardo, G.↗
DOE OSTI2023

Parallelized domain decomposition for multi-dimensional Lagrangian random walk mass-transfer particle tracking schemes

Lagrangian particle tracking schemes allow a wide range of flow and transport processes to be simulated accurately, but a major challenge is numerically implementing the inter-particle interactions in an efficient manner. This article develops a multi-dimensional, parallelized domain decomposition (DDC) strategy for mass-transfer particle tracking (MTPT) methods in which particles exchange mass dynamically. We show that this can be efficiently parallelized by employing large numbers of CPU cores to accelerate run times. In order to validate the approach and our theoretical predictions we focus our efforts on a well-known benchmark problem with pure diffusion, where analytical solutions in any number of dimensions are well established. In this work, we investigate different procedures for “tiling” the domain in two and three dimensions (2-D and 3-D), as this type of formal DDC construction is currently limited to 1-D. An optimal tiling is prescribed based on physical problem parameters and the number of available CPU cores, as each tiling provides distinct results in both accuracy and run time. We further extend the most efficient technique to 3-D for comparison, leading to an analytical discussion of the effect of dimensionality on strategies for implementing DDC schemes. Increasing computational resources (cores) within the DDC method produces a trade-off between inter-node communication and on-node work. For an optimally subdivided diffusion problem, the 2-D parallelized algorithm achieves nearly perfect linear speedup in comparison with the serial run-up to around 2700 cores, reducing a 5 h simulation to 8 s, while the 3-D algorithm maintains appreciable speedup up to 1700 cores.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1985

The structure of dilute combusting sprays

An experimental and theoretical study of drop processes in a turbulent flame is described. The experiments involved a monodisperse (105 and 180 micro m initial diameter) stream of methanol drops injected at the base of a turbulent methane-fueled diffusion flame burning in still air. The following measurements were made: mean and fluctuating phase velocities, mean drop number flux, drop-size distributions and mean gas-phase temperatures. Measurements were compared with predictions of two separated flow models: (1) deterministic separated flow, where drop-turbulence interactions are ignored; and (2) stochastic separated flow, where drop-turbulence interactions are considered using random-walk computations. The stochastic separated flow analysis yielded best agreement with measurements, since it provides for turbulent dispersion of drops which was important for present test conditions (and probably for most combusting sprays as well). Distinguishing the presence or absence of envelope flames around the drops, however, was relatively unimportant for present test conditions, since the drops spent most of their lifetime in fuel-rich regions of the flow where this distinction is irrelevant.

Shuen, J. S.↗
NASA NTRS1986

Analytic method for calculating properties of random walks on networks

A method for calculating the properties of discrete random walks on networks is presented. The method divides complex networks into simpler units whose contribution to the mean first-passage time is calculated. The simplified network is then further iterated. The method is demonstrated by calculating mean first-passage times on a segment, a segment with a single dangling bond, a segment with many dangling bonds, and a looplike structure. The results are analyzed and related to the applicability of the Einstein relation between conductance and diffusion.

Goldhirsch, I.↗
DOE OSTI2026

Substrate-Mediated Evaporation and Stochastic Evolution of Supported Au Nanoparticles

Here, we use in situ transmission electron microscopy with automated tracking to study supported gold nanoparticles (NPs) during high-temperature vacuum annealing. The average mass loss per NP is governed by a flat, nearly size-independent substrate-mediated evaporation profile. On top of this mean shrinkage, individual NPs show significant fluctuations in apparent growth or shrinkage, and NP volume follows a random-walk-like trajectory. To rationalize both the ensemble-mean behavior and the particle-resolved variability, we develop a self-consistent theory that couples substrate-mediated evaporation to collective 2D Ostwald-type mass exchange through a shared adatom field, described in terms of a renormalized screening length and background concentration. In the experimentally relevant regime, the theory predicts an approximately size-independent mean shrinkage rate and clarifies how net mass loss suppresses classical coarsening. Superimposed on this deterministic drift, we quantify stochastic volume trajectories and capture their fluctuation spectrum with a minimal Langevin description consistent with intermittent adatom attachment and detachment events. In addition, we characterize the lateral diffusive motion of NPs, which is responsible for their coalescence. Altogether, our results highlight that stochasticity is intrinsic at the nanoscale and that predicting the evolution of supported NPs at early and intermediate times requires a unified framework combining substrate-mediated evaporation, collective mass exchange, and stochastic fluctuations.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2022

Rotational Dynamics and Transition Mechanisms of Surface-Adsorbed Proteins

Assembly of biomolecules at solid-water interfaces requires molecules to traverse complex orientation-dependent energy landscapes through processes that are poorly understood, largely due to the dearth of in-situ single molecule measurements and statistical analyses of the rotational dynamics that define directional selection. Emerging capabilities in high-speed atomic force microscopy and machine learning have allowed us to directly determine the orientational energy landscape and observe and quantify the rotational dynamics for protein nanorods on the surface of muscovite mica under a variety of conditions. Comparisons with kinetic Monte Carlo simulations show that the transition rates between adjacent orientation-specific energetic minima can largely be understood through traditional models of in-plane Brownian rotation across a biased energy landscape, with resulting transition rates that are exponential in the energy-barriers between states. However, transitions between more distant angular states are decoupled from barrier height, with jump-size distributions showing a power-law decay that is characteristic of a non-classical Levy-flight random walk, indicating that large jumps are enabled by alternative modes of motion via activated states. The findings provide new insights into the dynamics of biomolecules at solid-liquid interfaces that lead to self-assembly, epitaxial matching and other orientationally anisotropic outcomes and define a general procedure for exploring such dynamics with implications for hybrid biomolecular-inorganic materials design.

Orientational energy landscapes, Rotational dynami↗