Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Random walk model”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Non-Markovian Hole Excess Noise in Avalanche Amorphous Selenium Thin Films

Enhancing the signal-to-noise ratio in avalanche photodiodes by utilizing impact ionization gain requires materials exhibiting low excess noise factors. Amorphous selenium (a-Se) as a wide bandgap at ~2.1 eV, a solid-state avalanche layer, demonstrates single-carrier hole impact ionization gain and manifests ultralow thermal generation rates. A comprehensive study of the history dependent and non-Markovian nature of hot hole transport in a-Se was modeled using a Monte Carlo (MC) random walk of single hole free flights, interrupted by instantaneous phonon, disorder, hole–dipole, and impact-ionization scattering interactions. The hole excess noise factors were simulated for 0.1–15 μm a-Se thin-films as a function of mean avalanche gain. The hole excess noise factors in a-Se decreases with an increase in electric field, impact ionization gain, and device thickness. The history dependent nature of branching of holes is explained using a Gaussian avalanche threshold distance distribution and the dead space distance, which increases determinism in the stochastic impact ionization process. An ultralow non-Markovian excess noise factor of ~1 was simulated for 100 nm a-Se thin films corresponding to avalanche gains of 1000. Future detector designs can utilize the nonlocal/non-Markovian nature of the hole avalanche in a-Se, to enable a true solid-state photomultiplier with noiseless gain.

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

Examining AGN UV/Optical Variability beyond the Simple Damped Random Walk

We present damped harmonic oscillator (DHO) light-curve modeling for a sample of 12,714 spectroscopically confirmed quasars in the Sloan Digital Sky Survey Stripe 82 region. DHO is a second-order continuous-time autoregressive moving-average process, which can be fully described using four independent parameters: a natural oscillation frequency (ω 0 ), a damping ratio (ξ), a characteristic perturbation timescale (τ perturb ), and an amplitude for the perturbing white noise (σ ϵ ). The asymptotic variability amplitude of a DHO process is quantified by σ DHO —a function of ω 0 , ξ, τ perturb , and σ ϵ . We find that both τ perturb and σ ϵ follow different dependencies with rest-frame wavelength (λ RF ) on either side of 2500 Å, whereas σ DHO follows a single power-law relation with λ RF . After correcting for wavelength dependence, σ DHO exhibits anticorrelations with both the Eddington ratio and the black hole mass, while τ perturb —with a typical value of days in the rest frame—shows an anticorrelation with the bolometric luminosity. Modeling active galactic nuclei (AGN) variability as a DHO offers more insight into the workings of accretion disks close to the supermassive black holes at the center of AGN. The newly discovered short-term variability (characterized by τ perturb and σ ϵ ) and its correlation with bolometric luminosity pave the way for new algorithms that will derive fundamental properties (e.g., Eddington ratio) of AGN using photometric data alone.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1992

Monitoring global geodynamical parameter variations with GPS and low earth satellites

A low earth satellite enhances global geodynamical parameters determination with GPS in two ways. First, it improves the GPS orbits, which in turn improve the estimates of other parameters. Secondly, a low earth satellite completes an orbit cycle in far shorter time (90 to 120 minutes) than do GPS satellites (12 hours); it observes more GPS satellites than a ground receiver does in shorter time and increases the correlation between GPS orbit errors. This reduces the error in the determination of nonrotational coordinate parameters, i.e., geocentric offset of ground tracking sites. Covariance results with the global geodynamical parameters modeled as constants, and as random-walk parameters to closer reflect the actual variations, are compared. The effects of using different GPS data quality and different ground tracking network are studied. The case of using two earth satellites in orthogonal orbital planes is also investigated.

Wu, S. C.↗
NASA NTRS1993

Generalized Levy-walk model for DNA nucleotide sequences

We propose a generalized Levy walk to model fractal landscapes observed in noncoding DNA sequences. We find that this model provides a very close approximation to the empirical data and explains a number of statistical properties of genomic DNA sequences such as the distribution of strand-biased regions (those with an excess of one type of nucleotide) as well as local changes in the slope of the correlation exponent alpha. The generalized Levy-walk model simultaneously accounts for the long-range correlations in noncoding DNA sequences and for the apparently paradoxical finding of long subregions of biased random walks (length lj) within these correlated sequences. In the generalized Levy-walk model, the lj are chosen from a power-law distribution P(lj) varies as lj(-mu). The correlation exponent alpha is related to mu through alpha = 2-mu/2 if 2 < mu < 3. The model is consistent with the finding of "repetitive elements" of variable length interspersed within noncoding DNA.

NASA Discipline Number 14-10↗
DOE OSTI2013

Interplay between inhibited transport and reaction in nanoporous materials

This work presents a detailed formulation of reaction and diffusion dynamics of molecules in confined pores such as mesoporous silica and zeolites. A general reaction-diffusion model and discrete Monte Carlo simulations are presented. Both transient and steady state behavior is covered. Failure of previous mean-field models for these systems is explained and discussed. A coarse-grained, generalized hydrodynamic model is developed that accurately captures the interplay between reaction and restricted transport in these systems. This method incorporates the non-uniform chemical diffusion behavior present in finite pores with multi-component diffusion. Two methods of calculating these diffusion values are developed: a random walk based approach and a driven diffusion model based on an extension of Fick's law. The effects of reaction, diffusion, pore length, and catalytic site distribution are investigated. In addition to strictly single file motion, quasi-single file diffusion is incorporated into the model to match a range of experimental systems. The connection between these experimental systems and model parameters is made through Langevin dynamics modeling of particles in confined pores.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Stochastic agent-based model for predicting turbine-scale raptor movements during updraft-subsidized directional flights

Rapid expansion of wind energy development across the world has highlighted the need to better understand turbine-caused avian mortality. The risk to golden eagles (Aquila chrysaetos) is of particular concern due to their small population size and conservation status. Golden eagles subsidize their flight in part by soaring in orographic updrafts, which can place them in conflict with wind turbines utilizing the same low-altitude wind resource. Understanding the behavior of soaring raptors in varying atmospheric conditions can therefore be relevant to predicting and mitigating their risk of collision. We present a predictive movement model that simulates individual paths of golden eagles during directional flight (such as migration) that is subsidized by orographic updraft. We modeled eagles in a 50 km by 50 km study area in Wyoming containing three wind power plants with documented golden eagle collisions with turbines. The movement model is applicable to any region where ground elevation is known at turbine scale (50 m) and wind conditions are known at facility scale (3 km). For a given set of atmospheric conditions, the model simulates movements of thousands of orographic soaring eagles to produce a density map quantifying the relative probability of eagle presence. We validated the simulated tracks with GPS telemetry data showing four directional tracks made by golden eagles transiting through the area in 2019 and 2020. For each eagle track, validation was performed using the ratio of the model-simulated eagle presence likelihood with uniform eagle presence and the presence computed using directed random-walk movements. We found that the predictive performance of the model was significantly better (likelihood ratio 1) for low-altitude movements than high-altitude movements that can involve thermal-soaring. We employed the model to produce seasonal presence maps for migrating golden eagles. We found significant turbine-level variations in eagle presence between northerly and southerly migration routes through the study area. Overall, the proposed model offers a generalizable, probabilistic, and predictive tool to assist wind energy developers, ecologists, wildlife managers, and industry consultants in estimating the potential for conflict between soaring birds and wind turbines, thereby reducing the need for site-specific data on golden eagle movements.

17 WIND ENERGY↗
NASA NTRS1985

Structure of ducted particle-laden turbulent jets

New and existing measurements have been used to develop and evaluate models of particle-laden jets and sprays, yielding encouraging results for a stochastic separated flow (SSF) model which treats effects of interphase slip and turbulence on the basis of random walk computations of particle or drop motion. An evaluation of the SSF model using the measurements of Modarress et al. (1982, 1983) yielded results comparable to earlier work, suggesting that data from coflowing jets in ducts should be accompanied by static pressure measurements, since even small pressure gradients appreciably influence jet properties.

Zhang, Q.-F.↗
NASA NTRS2008

Multiscale Modeling of Ablation and Pyrolysis in PICA-Like materials

During atmospheric entry of planetary probes, the thermal protection system (TIPS) of the probe is exposed to high temperatures under low pressures. In these conditions, carbonous fibrous TIPS materials may undergo oxidation leading to mass loss and wall recession called ablation. This work aims to improve the understanding of material/environment interactions through a study of the coupling between oxygen transport in the Knudsen regime, heterogeneous oxidation of carbon, and surface recession. A 3D Random Walk Monte Carlo simulation tool is used for this study. The fibrous architecture of a model material, consisting of high porosity random array of carbon fibers, is numerically represented on a 3D Cartesian grid. Mass transport in the Knudsen regime from the boundary layer to the surface, and inside this porous material is simulated by random walk. A reaction probability is used to simulate the heterogeneous oxidation reaction. The surface recession of the fibers is followed by front tracking using a simplified marching cube approach. The output data of the simulations are ablation velocity and dynamic evolution of the material porosity. A parametric study is carried out to analyze the material behavior as a function of Knudsen number for the porous media (length of the mean free path compared to the mean pore diameter) and the intrinsic reactivity of the carbon fibers. The model is applied to Stardust mission reentry conditions and explains the unexpected behavior of the TIPS material that underwent mass loss in volume.

Lachaud, Jean↗
NASA NTRS1984

Structure of nonevaporating sprays - Measurements and predictions

Structure measurements were completed within the dilute portion of axisymmetric nonevaporating sprays (SMD of 30 and 87 microns) injected into a still air environment, including: mean and fluctuating gas velocities and Reynolds stress using laser-Doppler anemometry; mean liquid fluxes using isokinetic sampling; drop sizes using slide impaction; and drop sizes and velocities using multiflash photography. The new measurements were used to evaluate three representative models of sprays: (1) a locally homogeneous flow (LHF) model, where slip between the phases was neglected; (2) a deterministic separated flow (DSF) model, where slip was considered but effects of drop interaction with turbulent fluctuations were ignored; and (3) a stochastic separated flow (SSF) model, where effects of both interphase slip and turbulent fluctuations were considered using random sampling for turbulence properties in conjunction with random-walk computations for drop motion. The LHF and DSF models were unsatisfactory for present test conditions-both underestimating flow widths and the rate of spread of drops. In contrast, the SSF model provided reasonably accurate predictions, including effects of enhanced spreading rates of sprays due to drop dispersion by turbulence, with all empirical parameters fixed from earlier work.

Solomon, A. S. P.↗
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.↗
DOE OSTI2012

Equations for the distributions of functionals of a random-walk trajectory in an inhomogeneous medium

Based on the random-trap model and using the mean-field approximation, we derive an equation that allows the distribution of a functional of the trajectory of a particle making random walks over inhomogeneous-lattice site to be calculated. The derived equation is a generalization of the Feynman-Kac equation to an inhomogeneous medium. We also derive a backward equation in which not the final position of the particle but its position at the initial time is used as an independent variable. As an example of applying the derived equations, we consider the one-dimensional problem of calculating the first-passage time distribution. We show that the average first-passage times for homogeneous and inhomogeneous media with identical diffusion coefficients coincide, but the variance of the distribution for an inhomogeneous medium can be many times larger than that for a homogeneous one.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

Nonlocal operators, parabolic-type equations, and ultrametric random walks

In this article, we introduce a new type of nonlocal operators and study the Cauchy problem for certain parabolic-type pseudodifferential equations naturally associated to these operators. Some of these equations are the p-adic master equations of certain models of complex systems introduced by Avetisov, V. A. and Bikulov, A. Kh., “On the ultrametricity of the fluctuation dynamicmobility of protein molecules,” Proc. Steklov Inst. Math. 265(1), 75–81 (2009) [Tr. Mat. Inst. Steklova 265, 82–89 (2009) (Izbrannye Voprosy Matematicheskoy Fiziki i p-adicheskogo Analiza) (in Russian)]; Avetisov, V. A., Bikulov, A. Kh., and Zubarev, A. P., “First passage time distribution and the number of returns for ultrametric random walks,” J. Phys. A 42(8), 085003 (2009); Avetisov, V. A., Bikulov, A. Kh., and Osipov, V. A., “p-adic models of ultrametric diffusion in the conformational dynamics of macromolecules,” Proc. Steklov Inst. Math. 245(2), 48–57 (2004) [Tr. Mat. Inst. Steklova 245, 55–64 (2004) (Izbrannye Voprosy Matematicheskoy Fiziki i p-adicheskogo Analiza) (in Russian)]; Avetisov, V. A., Bikulov, A. Kh., and Osipov, V. A., “p-adic description of characteristic relaxation in complex systems,” J. Phys. A 36(15), 4239–4246 (2003); Avetisov, V. A., Bikulov, A. H., Kozyrev, S. V., and Osipov, V. A., “p-adic models of ultrametric diffusion constrained by hierarchical energy landscapes,” J. Phys. A 35(2), 177–189 (2002); Avetisov, V. A., Bikulov, A. Kh., and Kozyrev, S. V., “Description of logarithmic relaxation by a model of a hierarchical random walk,” Dokl. Akad. Nauk 368(2), 164–167 (1999) (in Russian). The fundamental solutions of these parabolic-type equations are transition functions of random walks on the n-dimensional vector space over the field of p-adic numbers. We study some properties of these random walks, including the first passage time.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Communication: Distinguishing between short-time non-Fickian diffusion and long-time Fickian diffusion for a random walk on a crowded lattice

The motion of cells and molecules through biological environments is often hindered by the presence of other cells and molecules. A common approach to modeling this kind of hindered transport is to examine the mean squared displacement (MSD) of a motile tracer particle in a lattice-based stochastic random walk in which some lattice sites are occupied by obstacles. Unfortunately, stochastic models can be computationally expensive to analyze because we must average over a large ensemble of identically prepared realizations to obtain meaningful results. To overcome this limitation we describe an exact method for analyzing a lattice-based model of the motion of an agent moving through a crowded environment. Using our approach we calculate the exact MSD of the motile agent. Our analysis confirms the existence of a transition period where, at first, the MSD does not follow a power law with time. However, after a sufficiently long period of time, the MSD increases in proportion to time. This latter phase corresponds to Fickian diffusion with a reduced diffusivity owing to the presence of the obstacles. Our main result is to provide a mathematically motivated, reproducible, and objective estimate of the amount of time required for the transport to become Fickian. Our new method to calculate this crossover time does not rely on stochastic simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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.↗
DOE OSTI2016

A NON-PARAMETRIC APPROACH TO CONSTRAIN THE TRANSFER FUNCTION IN REVERBERATION MAPPING

Broad emission lines of active galactic nuclei stem from a spatially extended region (broad-line region, BLR) that is composed of discrete clouds and photoionized by the central ionizing continuum. The temporal behaviors of these emission lines are blurred echoes of continuum variations (i.e., reverberation mapping, RM) and directly reflect the structures and kinematic information of BLRs through the so-called transfer function (also known as the velocity-delay map). Based on the previous works of Rybicki and Press and Zu et al., we develop an extended, non-parametric approach to determine the transfer function for RM data, in which the transfer function is expressed as a sum of a family of relatively displaced Gaussian response functions. Therefore, arbitrary shapes of transfer functions associated with complicated BLR geometry can be seamlessly included, enabling us to relax the presumption of a specified transfer function frequently adopted in previous studies and to let it be determined by observation data. We formulate our approach in a previously well-established framework that incorporates the statistical modeling of continuum variations as a damped random walk process and takes into account long-term secular variations which are irrelevant to RM signals. The application to RM data shows the fidelity of our approach.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

The Role of Diffusion in the Transport of Energetic Electrons during Solar Flares

The transport of the energy contained in suprathermal electrons in solar flares plays a key role in our understanding of many aspects of flare physics, from the spatial distributions of hard X-ray emission and energy deposition in the ambient atmosphere to global energetics. Historically the transport of these particles has been largely treated through a deterministic approach, in which first-order secular energy loss to electrons in the ambient target is treated as the dominant effect, with second-order diffusive terms (in both energy and angle) generally being either treated as a small correction or even neglected. Here, we critically analyze this approach, and we show that spatial diffusion through pitch-angle scattering necessarily plays a very significant role in the transport of electrons. We further show that a satisfactory treatment of the diffusion process requires consideration of non-local effects, so that the electron flux depends not just on the local gradient of the electron distribution function but on the value of this gradient within an extended region encompassing a significant fraction of a mean free path. Our analysis applies generally to pitch-angle scattering by a variety of mechanisms, from Coulomb collisions to turbulent scattering. We further show that the spatial transport of electrons along the magnetic field of a flaring loop can be modeled rather effectively as a Continuous Time Random Walk with velocity-dependent probability distribution functions of jump sizes and occurrences, both of which can be expressed in terms of the scattering mean free path.

79 ASTRONOMY AND ASTROPHYSICS↗
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 OSTI2018

On a Local Version of the Bak–Sneppen Model

A major difficulty in studying the Bak–Sneppen model is in effectively comparing it with well-understood models. This stems from the use of two geometries: complete graph geometry to locate the global fitness minimizer, and graph geometry to replace the species in the neighborhood of the minimizer. Over the years a number of models inspired by Bak–Sneppen were studied, usually by introducing different or new features (e.g. discretizing fitness, randomized neighbors or population size). We present a variant that only uses features present in Bak–Sneppen, and whose difference from the Bak–Sneppen is that only the graph geometry is used for the evolution. This allows to obtain the stationary distribution through random walk dynamics while preserving the geometric nature of the model. We use this to show that for constant-degree graphs, the stationary fitness distribution converges to an IID law as the number of vertices tends to infinity. We also discuss exponential ergodicity through coupling, and avalanches for the model.

97 MATHEMATICS AND COMPUTING↗