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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS1973

Composition and physics of the lunar atmosphere.

The existence in the lunar atmosphere of helium, neon, argon, and possibly molecular hydrogen has been confirmed by the Apollo 17 mass spectrometer. The observed helium concentrations and distribution agree closely with model predictions for a non-condensable gas based on a solar wind source, thermal escape and a Monte Carlo random walk calculated longitudinal distribution. Heavier gases are lost by photoionization and subsequent sweeping away by the solar wind electric field. The observed nighttime neon concentration of 80,000 molecules per cu cm is consistent with expected amounts. Argon, however, is adsorbed on the lunar surface late at night when the surface temperature is lowest. It shows the expected predawn enhancement exhibited by condensable gases released into the atmosphere at the sunrise terminator. Hydrogen appears to exist in the molecular rather than atomic state. Its observed concentration is less than a factor of 3 higher than that predicted by a model similar to that used for helium.

Hoffman, J. H.↗
DOE OSTI2013

The role of the Kubo number in two-component turbulence

We explore the random walk of magnetic field lines in two-component turbulence by using computer simulations. It is often assumed that the two-component model provides a good approximation for solar wind turbulence. We explore the dependence of the field line diffusion coefficient on the Kubo number which is a fundamental and characteristic quantity in the theory of turbulence. We show that there are two transport regimes. One is the well-known quasilinear regime in which the diffusion coefficient is proportional to the Kubo number squared, and the second one is a nonlinear regime in which the diffusion coefficient is directly proportional to the Kubo number. The so-called percolative transport regime which is often discussed in the literature cannot be found. The numerical results obtained in the present paper confirm analytical theories for random walking field lines developed in the past.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Impact of confined geometries on hopping and trapping of motile bacteria in porous media

We use a random walk particle-tracking (RWPT) approach to elucidate the impact of porous media confinement and cell-cell interactions on bacterial transport. The model employs stochastic alternating motility states consisting of hopping movement and trapping reorientation. The stochastic motility patterns are defined based on direct visualization of individual trajectory data. We validate our model against experimental data, at single-cell resolution, of bacterial E. coli motion in three-dimensional confined porous media. Results show that the model is able to efficiently simulate the spreading dynamics of motile bacteria as it captures the impact of cell-cell interaction and pore confinement, which marks the transition to a late-time subdiffusive regime. Furthermore, the model is able to qualitatively reproduce the observed directional persistence. Lastly, our RWPT model constitutes a meshless simple method which is easy to implement and does not invoke ad hoc assumptions but represents the basis for a multiscale approach to the study of bacterial dispersal in porous systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Direct simulation Monte Carlo modeling of relaxation processes in polyatomic gases

Relaxation processes of polyatomic molecules are modeled and implemented in an in-house Direct Simulation Monte Carlo code in order to enable the simulation of atmospheric entry maneuvers at Mars and Saturn’s Titan. The description of rotational and vibrational relaxation processes is derived from basic quantum-mechanics using a rigid rotator and a simple harmonic oscillator, respectively. Strategies regarding the vibrational relaxation process are investigated, where good agreement for the relaxation time according to the Landau-Teller expression is found for both methods, the established prohibiting double relaxation method and the new proposed multi-mode relaxation. Differences and applications areas of these two methods are discussed. Consequently, two numerical methods used for sampling of energy values from multi-dimensional distribution functions are compared. The proposed random-walk Metropolis algorithm enables the efficient treatment of multiple vibrational modes within a time step with reasonable computational effort. The implemented model is verified and validated by means of simple reservoir simulations and the comparison to experimental measurements of a hypersonic, carbon-dioxide flow around a flat-faced cylinder.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Improved approximations of non-Gaussian probability, force, and energy of a single polymer chain

The mechanical behavior of polymers such as their probability density function (PDF), strain energy, and entropic force, has long been described by the non-Gaussian statistical model. Non-Gaussian models are often approximated by the Kuhn-Grün (KG) distribution function, which is derived from the first-order approximation of the complex Rayleigh's exact Fourier integral distribution. Furthermore, the KG function is widely accepted in polymer physics, where the non-Gaussian theory is often used to describe energy of the chains with various flexibility ratios. However, the KG function is shown to be relevant only for long chains and becomes extremely inaccurate for chains with fewer than 40 segments. In comparison to the KG model, other approximations of the non-Gaussian statistical model are often less accurate, and those with higher accuracy are usually too complex to be implemented in large-scale simulations. Here a new accurate approximation of the non-Gaussian PDF, entropic force, and strain energy of a single chain subsequently is developed to describe the mechanics of a polymer chain. With a similar level of complexity, the presented approximations of non-Gaussian PDF, strain energy, and entropic force are at least 10 times more accurate than KG approximations and thus are an excellent alternative option to be used in micromechanical constitutive models.

42 ENGINEERING↗
DOE OSTI2022

A Differentiable Hydrology Approach for Modeling With Time-Varying Parameters

In the analysis of observational data with mathematical models, it is often desirable to let one or more model parameters vary with time to account for changing environmental conditions or to offer more flexibility. A standard approach for dynamic models is to allow for the number of parameters to grow as a function of time, leading to a high-dimensional inference problem for long data records. We propose using Hamiltonian Monte Carlo, a gradient-based Markov chain Monte Carlo (MCMC) method, to perform Bayesian inference for time-varying parameters (TVP) in hydrology models. As derivatives of model error functions with regard to parameters are not available in closed form, we implemented the GR4J rainfall-runoff models in Theano/PyMC3, and Jax/NumPyro, allowing for the application of automatic differentiation techniques. A simulation experiment assessing the viability of these methods for recovering underlying temporal variation in parameters indicates that a discrete Gaussian random walk-based prior appears to be best suited across a range of scenarios while the continuous-time Gaussian process compares relatively poorly. In an analysis of 20 years of daily real-world streamflow records from the Model Parameter Estimation Experiment (MOPEX) hydrology data set with increasing data sparsity and at multiple temporal resolutions, we find that TVP inference with gradient-based MCMC is a flexible and reliable approach for analyzing hydrology models with dynamic parameter sets.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2013

THE ANTICORRELATED NATURE OF THE PRIMARY AND SECONDARY ECLIPSE TIMING VARIATIONS FOR THE KEPLER CONTACT BINARIES

We report a study of the eclipse timing variations in contact binary systems, using long-cadence lightcurves from the Kepler archive. As a first step, observed minus calculated (O - C) curves were produced for both the primary and secondary eclipses of some 2000 Kepler binaries. We find {approx}390 short-period binaries with O - C curves that exhibit (1) random walk-like variations or quasi-periodicities, with typical amplitudes of {+-}200-300 s, and (2) anticorrelations between the primary and secondary eclipse timing variations. We present a detailed analysis and results for 32 of these binaries with orbital periods in the range of 0.35 {+-} 0.05 days. The anticorrelations observed in their O - C curves cannot be explained by a model involving mass transfer, which, among other things, requires implausibly high rates of {approx}0.01 M{sub Sun} yr{sup -1}. We show that the anticorrelated behavior, the amplitude of the O - C delays, and the overall random walk-like behavior can be explained by the presence of a starspot that is continuously visible around the orbit and slowly changes its longitude on timescales of weeks to months. The quasi-periods of {approx}50-200 days observed in the O - C curves suggest values for k, the coefficient of the latitude dependence of the stellar differential rotation, of {approx}0.003-0.013.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2011

River meander modeling and confronting uncertainty.

This study examines the meandering phenomenon as it occurs in media throughout terrestrial, glacial, atmospheric, and aquatic environments. Analysis of the minimum energy principle, along with theories of Coriolis forces (and random walks to explain the meandering phenomenon) found that these theories apply at different temporal and spatial scales. Coriolis forces might induce topological changes resulting in meandering planforms. The minimum energy principle might explain how these forces combine to limit the sinuosity to depth and width ratios that are common throughout various media. The study then compares the first order analytical solutions for flow field by Ikeda, et al. (1981) and Johannesson and Parker (1989b). Ikeda's et al. linear bank erosion model was implemented to predict the rate of bank erosion in which the bank erosion coefficient is treated as a stochastic variable that varies with physical properties of the bank (e.g., cohesiveness, stratigraphy, or vegetation density). The developed model was used to predict the evolution of meandering planforms. Then, the modeling results were analyzed and compared to the observed data. Since the migration of a meandering channel consists of downstream translation, lateral expansion, and downstream or upstream rotations several measures are formulated in order to determine which of the resulting planforms is closest to the experimental measured one. Results from the deterministic model highly depend on the calibrated erosion coefficient. Since field measurements are always limited, the stochastic model yielded more realistic predictions of meandering planform evolutions. Due to the random nature of bank erosion coefficient, the meandering planform evolution is a stochastic process that can only be accurately predicted by a stochastic model.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

The expected anisotropy in solid inflation

Solid inflation is an effective field theory of inflation in which isotropy and homogeneity are accomplished via a specific combination of anisotropic sources (three scalar fields that individually break isotropy). This results in specific observational signatures that are not found in standard models of inflation: a non-trivial angular dependence for the squeezed bispectrum, and a possibly long period of anisotropic inflation (to drive inflation, the ''solid'' must be very insensitive to any deformation, and thus background anisotropies are very slowly erased). In this paper we compute the expected level of statistical anisotropy in the power spectrum of the curvature perturbations of this model. To do so, we account for the classical background values of the three scalar fields that are generated on large (superhorizon) scales during inflation via a random walk sum, as the perturbation modes leave the horizon. Such an anisotropy is unavoidably generated, even starting from perfectly isotropic classical initial conditions. The expected level of anisotropy is related to the duration of inflation and to the amplitude of the squeezed bispectrum. If this amplitude is close to its current observational limit (so that one of the most interesting predictions of the model can be observed in the near future), we find that a level of statistical anisotropy F{sup 2} gives frozen and scale invariant vector perturbations on superhorizon scales.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

TURBULENCE IN THE SOLAR WIND MEASURED WITH COMET TAIL TEST PARTICLES

By analyzing the motions of test particles observed remotely in the tail of Comet Encke, we demonstrate that the solar wind undergoes turbulent processing enroute from the Sun to the Earth and that the kinetic energy entrained in the large-scale turbulence is sufficient to explain the well-known anomalous heating of the solar wind. Using the heliospheric imaging (HI-1) camera on board NASA's STEREO-A spacecraft, we have observed an ensemble of compact features in the comet tail as they became entrained in the solar wind near 0.4 AU. We find that the features are useful as test particles, via mean-motion analysis and a forward model of pickup dynamics. Using population analysis of the ensemble's relative motion, we find a regime of random-walk diffusion in the solar wind, followed, on larger scales, by a surprising regime of semiconfinement that we attribute to turbulent eddies in the solar wind. The entrained kinetic energy of the turbulent motions represents a sufficient energy reservoir to heat the solar wind to observed temperatures at 1 AU. We determine the Lagrangian-frame diffusion coefficient in the diffusive regime, derive upper limits for the small scale coherence length of solar wind turbulence, compare our results to existing Eulerian-frame measurements, and compare the turbulent velocity with the size of the observed eddies extrapolated to 1 AU. We conclude that the slow solar wind is fully mixed by turbulence on scales corresponding to a 1–2 hr crossing time at Earth; and that solar wind variability on timescales shorter than 1–2 hr is therefore dominated by turbulent processing rather than by direct solar effects.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1981

The stochastic evolution of asteroidal regoliths and the origin of brecciated and gas-rich meteorites

A model is constructed which views regolith evolution on asteroids as a stochastic process. Average values are shown to be poor descriptors of regolith depth. The utility of the average depth is not significantly increased by avoiding large craters or thick ejecta deposits, a procedure adopted in previous regolith studies. The statistical uncertainty associated with regolith depth severely limits the power of regolith models in predicting parent-body size for brecciated meteorites. A Monte Carlo algorithm was used to simulate the random walks and corresponding charged-particle irradiation histories of grains in regoliths. On rocky asteroids, only about 20 percent of the grains was exposed to solar cosmic ray ions. Results based on present-day conditions in the asteroid belt agree well with irradiation features observed in gas-rich meteorites. An origin during epochs of early solar system evolution is not required.

Housen, K. R.↗
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.↗
DOE OSTI2018

Construction of distinct discrete time scattering quantum walk formulations on the honeycomb lattice

Here we address the scattering construction of discrete time quantum walks on the honeycomb lattice. We write the system general (unitary) one step time evolution operator in terms of: (i) complete arbitrary scattering matrices Γ{sup (j,k)}, defined on the sites (j,k) of the lattice and; and (ii) topological directional functions Φ, which represent the distinct ways we can keep track of the propagation directions along the three bonds attached to each (j,k). By imposing the Φ’s to comply with all the honeycomb translational and point group symmetries (in the case of Γ{sup (j,k)}=Γ), we obtain in total ten independent model versions, presenting different dynamical features. To study some of their traits, we consider the idea of ‘characteristic paths’ CPs (closely related to classical random walks), determining the CPs for each one of the ten formulations. We then discuss many numerical examples of time evolution dynamics for our scattering quantum walks. For the calculations we use few Γ’s, as the Discrete Fourier Transform (DFT) and Grover matrices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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.↗
DOE OSTI2017

On-the-Fly Generation of Differential Resonance Scattering Probability Distribution Functions for Monte Carlo Codes

Current Monte Carlo codes use one of three models: (1) the asymptotic scattering model, (2) the free gas scattering model, or (3) the S(α,β) model, depending on the neutron energy and the specific Monte Carlo code. This thesis addresses the consequences of using the free gas scattering model, which assumes that the neutron interacts with atoms in thermal motion in a monatomic gas in thermal equilibrium at material temperature, T. Most importantly, the free gas model assumes the scattering cross section is constant over the neutron energy range, which is usually a good approximation for light nuclei, but not for heavy nuclei where the scattering cross section may have several resonances in the epithermal region. Several researchers in the field have shown that the exact resonance scattering model is temperaturedependent, and neglecting the resonances in the lower epithermal range can under-predict resonance absorption due to the upscattering phenomenon mentioned above, leading to an over-prediction of keff by several hundred pcm. Existing methods to address this issue involve changing the neutron weights or implementing an extra rejection scheme in the free gas sampling scheme, and these all involve performing the collision analysis in the center-of-mass frame, followed by a conversion back to the laboratory frame to continue the random walk of the neutron. The goal of this paper was to develop a sampling methodology that (1) accounted for the energydependent scattering cross sections in the collision analysis and (2) was performed in the laboratory frame,avoiding the conversion to the center-of-mass frame. The energy dependence of the scattering cross section was modeled with even-ordered polynomials (2nd and 4th order) to approximate the scattering cross section in Blackshaw’s equations for the moments of the differential scattering PDFs. These moments were used to sample the outgoing neutron speed and angle in the laboratory frame on-the-fly during the random walk of the neutron. Results for criticality studies on fuel pin and fuel assembly calculations using methods developed in this dissertation showed very close comparison to results using the reference Dopplerbroadened rejection correction (DBRC) scheme.

42 ENGINEERING↗
DOE OSTI2024

Cylinders’ percolation: Decoupling and applications

In this paper we establish a strong decoupling inequality for the cylinder’s percolation process introduced by Tykesson and Windisch (Probab. Theory Related Fields 154 (2012) 165–191). This model features a very strong dependency structure, making it difficult to study, and this is why such decoupling inequalities are desirable. It is important to notice that the type of dependencies featured by cylinder’s percolation is particularly intricate, given that the cylinders have infinite range (unlike some models like Boolean percolation) while at the same time being rigid bodies (unlike processes such as random interlacements). Here our work introduces a new notion of fast decoupling, proves that it holds for the model in question and finishes with an application. More precisely, we prove that for a small enough density of cylinders, a random walk on a connected component of the vacant set is transient for all dimensions d≥3.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Photon Next-Event Estimator Implementation in MCATK

A photon next-event fluence estimator at a point has been implemented in the Monte Carlo Application Toolkit (MCATK). The next-event estimator provides an expected value estimator for the flux at a point due to all source and collision events. An advantage of the next-event estimator over track-length estimators, which are normally employed in MCATK, is that flux estimates can be made in locations that have no random walk particle tracks. The next-event estimator allows users to calculate radiographs and estimate response for detectors outside of the modeled geometry. The next-event estimator is not yet accessable through the MCATK FlatAPI for C and Fortran. The next-event estimator in MCATK has been tested against MCNP6 using 5 suites of test problems. No issues were found in the MCATK implementation. One issue was found in the exclusion radius approximation in MCNP6. The theory, implementation, and testing are described in this document.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗