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 271 records · Page 15
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.↗
DOE OSTI2011

TIME-DEPENDENT PERPENDICULAR TRANSPORT OF FAST CHARGED PARTICLES IN A TURBULENT MAGNETIC FIELD

We present an analytic derivation of the temporal dependence of the perpendicular transport coefficient of charged particles in magnetostatic turbulence, for times smaller than the time needed for charged particles to travel the turbulence correlation length. This time window is left unexplored in most transport models. In our analysis all magnetic scales are taken to be much larger than the particle gyroradius, so that perpendicular transport is assumed to be dominated by the guiding center motion. Particle drift from the local magnetic field lines (MFLs) and magnetic field line random walk are evaluated separately for slab and three-dimensional (3D) isotropic turbulence. Contributions of wavelength scales shorter and longer than the turbulence coherence length are compared. In contrast to the slab case, particles in 3D isotropic turbulence unexpectedly diffuse from local MFLs; this result questions the common assumption that particle magnetization is independent of turbulence geometry. Extensions of this model will allow for a study of solar wind anisotropies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Universality Classes of Interaction Structures for NK Fitness Landscapes

Kauffman’s NK-model is a paradigmatic example of a class of stochastic models of genotypic fitness landscapes that aim to capture generic features of epistatic interactions in multilocus systems. Genotypes are represented as sequences of L binary loci. The fitness assigned to a genotype is a sum of contributions, each of which is a random function defined on a subset of k≤L loci. These subsets or neighborhoods determine the genetic interactions of the model. Whereas earlier work on the NK model suggested that most of its properties are robust with regard to the choice of neighborhoods, recent work has revealed an important and sometimes counter-intuitive influence of the interaction structure on the properties of NK fitness landscapes. Here we review these developments and present new results concerning the number of local fitness maxima and the statistics of selectively accessible (that is, fitness-monotonic) mutational pathways. In particular, we develop a unified framework for computing the exponential growth rate of the expected number of local fitness maxima as a function of L, and identify two different universality classes of interaction structures that display different asymptotics of this quantity for large k. Moreover, we show that the probability that the fitness landscape can be traversed along an accessible path decreases exponentially in L for a large class of interaction structures that we characterize as locally bounded. Finally, we discuss the impact of the NK interaction structures on the dynamics of evolution using adaptive walk models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Latent Stochastic Differential Equations for Modeling Quasar Variability and Inferring Black Hole Properties

Quasars are bright and unobscured active galactic nuclei (AGN) thought to be powered by the accretion of matter around supermassive black holes at the centers of galaxies. The temporal variability of a quasar’s brightness contains valuable information about its physical properties. The UV/optical variability is thought to be a stochastic process, often represented as a damped random walk described by a stochastic differential equation (SDE). Upcoming wide-field telescopes such as the Rubin Observatory Legacy Survey of Space and Time (LSST) are expected to observe tens of millions of AGN in multiple filters over a ten year period, so there is a need for efficient and automated modeling techniques that can handle the large volume of data. Latent SDEs are machine learning models well suited for modeling quasar variability, as they can explicitly capture the underlying stochastic dynamics. In this work, we adapt latent SDEs to jointly reconstruct multivariate quasar light curves and infer their physical properties such as the black hole mass, inclination angle, and temperature slope. Our model is trained on realistic simulations of LSST ten year quasar light curves, and we demonstrate its ability to reconstruct quasar light curves even in the presence of long seasonal gaps and irregular sampling across different bands, outperforming a multioutput Gaussian process regression baseline. Our method has the potential to provide a deeper understanding of the physical properties of quasars and is applicable to a wide range of other multivariate time series with missing data and irregular sampling.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Precipitation and Release of Solar Energetic Particles from the Solar Coronal Magnetic Field

Most solar energetic particles (SEPs) are produced in the corona. They propagate through complex coronal magnetic fields subject to scattering and diffusion across the averaged field lines by turbulence. We examine the behaviors of particle transport using a stochastic 3D focused transport simulation in a potential field source surface model of coronal magnetic field. The model is applied to an SEP event on 2010 February 7. We study three scenarios of particle injection at (i) the compact solar flare site, (ii) the coronal mass ejection (CME) shock, and (iii) the EUV wave near the surface. The majority of particles injected on open field lines are able to escape the corona. We found that none of our models can explain the observations of wide longitudinal SEP spread without perpendicular diffusion. If the perpendicular diffusion is about 10% of what is derived from the random walk of field lines at the rate of supergranular diffusion, particles injected at the compact solar flare site can spread to a wide range of longitude and latitude, very similar to the behavior of particles injected at a large CME shock. Stronger pitch-angle scattering results in a little more lateral spread by holding the particles in the corona for longer periods of time. Some injected particles eventually end up precipitating onto the solar surface. Even with a very small perpendicular diffusion, the pattern of the particle precipitation can be quite complicated depending on the detailed small-scale coronal magnetic field structures, which could be seen with future sensitive gamma-ray telescopes.

79 ASTRONOMY AND ASTROPHYSICS↗
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.↗
DOE OSTI2023

Dark matter halo mass functions and density profiles from mass and energy cascade

Abstract Halo abundance and structure play a central role for modeling structure formation and evolution. Without relying on a spherical or ellipsoidal collapse model, we analytically derive the halo mass function and cuspy halo density (inner slope of −4/3) based on the mass and energy cascade theory in dark matter flow. The hierarchical halo structure formation leads to halo or particle random walk with a position-dependent waiting time $$\tau _g$$ τ g . First, the inverse mass cascade from small to large scales leads to the halo random walk in mass space with $$\tau _g\propto m_h^{-\lambda }$$ τ g ∝ m h - λ , where $$m_h$$ m h is the halo mass and $$\lambda$$ λ is a halo geometry parameter with predicted value of 2/3. The corresponding Fokker-Planck solution for halo random walk in mass space gives rise to the halo mass function with a power-law behavior on small scale and exponential decay on large scale. This can be further improved by considering two different $$\lambda$$ λ for haloes below and above a critical mass scale $$m_h^*$$ m h ∗ , i.e. a double- $$\lambda$$ λ halo mass function. Second, a double- $$\gamma$$ γ density profile can be derived based on the particle random walk in 3D space with a position-dependent waiting time $$\tau _g \propto \Phi (r)^{-1} \propto r^{-\gamma }$$ τ g ∝ Φ ( r ) - 1 ∝ r - γ , where $$\Phi$$ Φ is the gravitational potential and r is the particle distance to halo center. Theory predicts $$\gamma =2/3$$ γ = 2 / 3 that leads to a cuspy density profile with an inner slope of −4/3, consistent with the predicted scaling laws from energy cascade. The Press-Schechter mass function and Einasto density profile are just special cases of proposed models. The small scale permanence can be identified due to the scale-independent rate of mass and energy cascade, where density profiles of different halo masses and redshifts converge to the $$-4/3$$ - 4 / 3 scaling law ( $$\rho _h \propto r^{-4/3}$$ ρ h ∝ r - 4 / 3 ) on small scales. Theory predicts the halo number density scales with halo mass as $$\propto m_h^{-1.9}$$ ∝ m h - 1.9 , while the halo mass density scales as $$\propto m_h^{4/9}$$ ∝ m h 4 / 9 . Results were compared against the Illustris simulations. This new perspective provides a theory for nearly universal halo mass functions and density profiles.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

On collisional diffusion in a stochastic magnetic field

The effect of particle collisions on the transport in a stochastic magnetic field in tokamaks is investigated. The model of resonant magnetic perturbations generated by external coils at the plasma edge is used for the stochastic magnetic field. The particle collisions are simulated by a random walk process along the magnetic field lines and the jumps across the field lines at the collision instants. The dependencies of the local diffusion coefficients on the mean free path λ{sub mfp}, the diffusion coefficients of field lines D{sub FL}, and the collisional diffusion coefficients, χ{sub ⊥} are studied. Based on these numerical data and the heuristic arguments, the empirical formula, D{sub r}=χ{sub ⊥}+v{sub ||}D{sub FL}/(1+L{sub c}/λ{sub mfp}), for the local diffusion coefficient is proposed, where L{sub c} is the characteristic length of order of the connection length l{sub c}=πqR{sub 0}, q is the safety factor, R{sub 0} is the major radius. The formula quite well describes the results of numerical simulations. In the limiting cases, the formula describes the Rechester-Rosenbluth and Laval scalings.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Characterizing system dynamics with a weighted and directed network constructed from time series data

In this work, we propose a novel method to transform a time series into a weighted and directed network. For a given time series, we first generate a set of segments via a sliding window, and then use a doubly symbolic scheme to characterize every windowed segment by combining absolute amplitude information with an ordinal pattern characterization. Based on this construction, a network can be directly constructed from the given time series: segments corresponding to different symbol-pairs are mapped to network nodes and the temporal succession between nodes is represented by directed links. With this conversion, dynamics underlying the time series has been encoded into the network structure. We illustrate the potential of our networks with a well-studied dynamical model as a benchmark example. Results show that network measures for characterizing global properties can detect the dynamical transitions in the underlying system. Moreover, we employ a random walk algorithm to sample loops in our networks, and find that time series with different dynamics exhibits distinct cycle structure. That is, the relative prevalence of loops with different lengths can be used to identify the underlying dynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1999

Approximation for Horizontal Photon Transport in Cloud Remote Sensing Problems

The effect of horizontal photon transport within real-world clouds can be of consequence to remote sensing problems based on plane-parallel cloud models. An analytic approximation for the root-mean-square horizontal displacement of reflected and transmitted photons relative to the incident cloud-top location is derived from random walk theory. The resulting formula is a function of the average number of photon scatterings, and particle asymmetry parameter and single scattering albedo. In turn, the average number of scatterings can be determined from efficient adding/doubling radiative transfer procedures. The approximation is applied to liquid water clouds for typical remote sensing solar spectral bands, involving both conservative and non-conservative scattering. Results compare well with Monte Carlo calculations. Though the emphasis is on horizontal photon transport in terrestrial clouds, the derived approximation is applicable to any multiple scattering plane-parallel radiative transfer problem. The complete horizontal transport probability distribution can be described with an analytic distribution specified by the root-mean-square and average displacement values. However, it is shown empirically that the average displacement can be reasonably inferred from the root-mean-square value. An estimate for the horizontal transport distribution can then be made from the root-mean-square photon displacement alone.

Plantnick, Steven↗
DOE OSTI2014

History dependent quantum random walks as quantum lattice gas automata

Quantum Random Walks (QRW) were first defined as one-particle sectors of Quantum Lattice Gas Automata (QLGA). Recently, they have been generalized to include history dependence, either on previous coin (internal, i.e., spin or velocity) states or on previous position states. These models have the goal of studying the transition to classicality, or more generally, changes in the performance of quantum walks in algorithmic applications. We show that several history dependent QRW can be identified as one-particle sectors of QLGA. This provides a unifying conceptual framework for these models in which the extra degrees of freedom required to store the history information arise naturally as geometrical degrees of freedom on the lattice.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Hydrodynamics of operator spreading and quasiparticle diffusion in interacting integrable systems

In this work, we address the hydrodynamics of operator spreading in interacting integrable lattice models. In these models, operators spread through the ballistic propagation of quasiparticles, with an operator front whose velocity is locally set by the fastest quasiparticle velocity. In interacting integrable systems, this velocity depends on the density of the other quasiparticles, so equilibrium density fluctuations cause the front to follow a biased random walk, and therefore to broaden diffusively. Ballistic front propagation and diffusive front broadening are also generically present in nonintegrable systems in one dimension; thus, although the mechanisms for operator spreading are distinct in the two cases, these coarse-grained measures of the operator front do not distinguish between the two cases. We present an expression for the front-broadening rate; we explicitly derive this for a particular integrable model (the “Floquet-Fredrickson-Andersen” model), and argue on kinetic grounds that it should apply generally. In conclusion, our results elucidate the microscopic mechanism for diffusive corrections to ballistic transport in interacting integrable models.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Turbulent-like flows in quasi two-dimensional dense suspensions of motile colloids

Dense bacterial suspensions exhibit turbulent-like flows at low Reynolds numbers, driven by the activity of the microswimmers. In this study, we develop a model system to examine these dynamics using motile colloids that mimic bacterial locomotion. The colloids are powered by the Quincke instability, which causes them to spontaneously roll in a random-walk pattern when exposed to a square-wave electric field. We experimentally investigate the flow dynamics in dense suspensions of these Quincke random walkers under quasi two-dimensional conditions, where the particle size is comparable to the gap between the electrodes. The results suggest a scaling regime in the energy spectrum ∼k −4 at high wavenumbers, observed consistently across activity levels and particle concentrations. We observe that velocity time correlations decay within a single period of the square-wave field, yet an anti-correlation appears between successive field applications, indicative of a dynamic structural memory of the ensemble.

Luo, Rui [Northwestern Univ., Evanston, IL (United↗
DOE OSTI2018

Near Critical Preferential Attachment Networks have Small Giant Components

Preferential attachment networks with power law exponent τ>3 are known to exhibit a phase transition. There is a value ρ{sub c}>0 such that, for small edge densities ρ≤ρ{sub c} every component of the graph comprises an asymptotically vanishing proportion of vertices, while for large edge densities ρ>ρ{sub c} there is a unique giant component comprising an asymptotically positive proportion of vertices. In this paper we study the decay in the size of the giant component as the critical edge density is approached from above. We show that the size decays very rapidly, like exp (−c/√(ρ−ρ{sub c})) for an explicit constant c>0 depending on the model implementation. This result is in contrast to the behaviour of the class of rank-one models of scale-free networks, including the configuration model, where the decay is polynomial. Our proofs rely on the local neighbourhood approximations of Dereich and Mörters (Ann Probab 41(1):329–384, 2013) and recent progress in the theory of branching random walks (Gantert et al. in Ann Inst Henri Poincaré Probab Stat 47(1):111–129, 2011).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1976

Energy distributions and radiation transport in uranium plasmas

An approximate analytic model, based on continuous electron slowing, has been used for survey calculations. Where more accuracy is required, a Monte Carlo technique is used which combines an analytic representation of Coulombic collisions with a random walk treatment of inelastic collisions. The calculated electron distributions have been incorporated into another code that evaluates both the excited atomic state densities within the plasma and the radiative flux emitted from the plasma.

Miley, G. H.↗
DOE OSTI2011

Brownian Shape Motion on Five-Dimensional Potential-Energy Surfaces:Nuclear Fission-Fragment Mass Distributions

Although nuclear fission can be understood qualitatively as an evolution of the nuclear shape, a quantitative description has proven to be very elusive. In particular, until now, there existed no model with demonstrated predictive power for the fission-fragment mass yields. Exploiting the expected strongly damped character of nuclear dynamics, we treat the nuclear shape evolution in analogy with Brownian motion and perform random walks on five-dimensional fission potential-energy surfaces which were calculated previously and are the most comprehensive available. Test applications give good reproduction of highly variable experimental mass yields. This novel general approach requires only a single new global parameter, namely, the critical neck size at which the mass split is frozen in, and the results are remarkably insensitive to its specific value.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Curvature perturbations from stochastic particle production during inflation

We calculate the curvature power spectrum sourced by spectator fields that are excited repeatedly and non-adiabatically during inflation. In the absence of detailed information of the nature of spectator field interactions, we consider an ensemble of models with intervals between the repeated interactions and interaction strengths drawn from simple probabilistic distributions. We show that the curvature power spectra of each member of the ensemble shows rich structure with many features, and there is a large variability between different realizations of the same ensemble. Such features can be probed by the cosmic microwave background (CMB) and large scale structure observations. They can also have implications for primordial black hole formation and CMB spectral distortions. The geometric random walk behavior of the spectator field allows us to calculate the ensembleaveraged power spectrum of curvature perturbations semi-analytically. For sufficiently large stochastic sourcing, the ensemble averaged power spectrum shows a scale dependence arising from the time spent by modes outside the horizon during the period of particle production, in spite of there being no preferred scale in the underlying model. We find that the magnitude of the ensemble-averaged power spectrum overestimates the typical power spectra in the ensemble because the ensemble distribution of the power spectra is highly non-Gaussian with fat tails.

79 ASTRONOMY AND ASTROPHYSICS↗