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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Grid-based diffusion Monte Carlo for fermions without the fixed-node approximation

A diffusion Monte Carlo algorithm is introduced that can determine the correct nodal structure of the wave function of a few-fermion system and its ground-state energy without an uncontrolled bias. This is achieved by confining signed random walkers to the points of a uniform infinite spatial grid, allowing them to meet and annihilate one another to establish the nodal structure without the fixed-node approximation. An imaginary-time propagator is derived rigorously from a discretized Hamiltonian, governing a non-Gaussian, sign-flipping, branching, and mutually annihilating random walk of particles. The accuracy of the resulting stochastic representations of a fermion wave function is limited only by the grid and imaginary-time resolutions and can be improved in a controlled manner. Here, the method is tested for a series of model problems including fermions in a harmonic trap as well as the He atom in its singlet or triplet ground state. For the latter case, the energies approach from above with increasing grid resolution and converge within 0.015 E h of the exact basis-set-limit value for the grid spacing of 0.08 a.u. with a statistical uncertainty of 10 –5 E h without an importance sampling or Jastrow factor.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Optimizing temperature distributions for training neural quantum states using parallel tempering

Parametrized artificial neural networks (ANNs) can be very expressive ansatzes for variational algorithms, reaching state-of-the-art energies on many quantum many-body Hamiltonians. Nevertheless, the training of the ANN can be slow and stymied by the presence of local minima in the parameter landscape. One approach to mitigate this issue is to use parallel tempering methods, and in this work, we focus on the role played by the temperature distribution of the parallel tempering replicas. Using an adaptive method that adjusts the temperatures in order to equate the exchange probability between neighboring replicas, we show that this temperature optimization can significantly increase the success rate of the variational algorithm with negligible computational cost by eliminating bottlenecks in the replicas' random walk. Furthermore, we demonstrate this using two different neural networks, a restricted Boltzmann machine and a feedforward network, which we use to study a toy problem based on a permutation invariant Hamiltonian with a pernicious local minimum and the 𝐽 1 −𝐽 2 model on a rectangular lattice.

Neural network simulations↗
NASA NTRS1992

Lifetime of binary asteroids versus gravitational encounters and collisions

We investigate the effect on the dynamics of a binary asteroid in the case of a near encounter with a third body. The dynamics of the binary is modeled as a two-body problem perturbed by an approaching body in the following ways: near encounters and collisions with a component of the system. In each case, the typical value of the two-body energy variation is estimated, and a random walk for the cumulative effect is assumed. Results are applied to some binary asteroid candidates. The main conclusion is that the collisional disruption is the dominant effect, giving lifetimes comparable to or larger than the age of the solar system.

Chauvineau, Bertrand↗
DOE OSTI2015

Study of electron transport in a Hall thruster by axial–radial fully kinetic particle simulation

Electron transport across a magnetic field in a magnetic-layer-type Hall thruster was numerically investigated for the future predictive modeling of Hall thrusters. The discharge of a 1-kW-class magnetic-layer-type Hall thruster designed for high-specific-impulse operation was modeled using an r-z two-dimensional fully kinetic particle code with and without artificial electron-diffusion models. The thruster performance results showed that both electron transport models captured the experimental result within discrepancies less than 20% in thrust and discharge current for all the simulated operation conditions. The electron cross-field transport mechanism of the so-called anomalous diffusion was self-consistently observed in the simulation without artificial diffusion models; the effective electron mobility was two orders of magnitude higher than the value obtained using the classical diffusion theory. To account for the self-consistently observed anomalous transport, the oscillation of plasma properties was speculated. It was suggested that the enhanced random-walk diffusion due to the velocity oscillation of low-frequency electron flow could explain the observed anomalous diffusion within an order of magnitude. The dominant oscillation mode of the electron flow velocity was found to be 20 kHz, which was coupled to electrostatic oscillation excited by global ionization instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

The nanocaterpillar's random walk: diffusion with ligand–receptor contacts

Particles with ligand–receptor contacts bind and unbind fluctuating “legs” to surfaces, whose fluctuations cause the particle to diffuse. Quantifying the diffusion of such “nanoscale caterpillars” is a challenge, since binding events often occur on very short time and length scales. Here, in this work, we derive an analytical formula, validated by simulations, for the long time translational diffusion coefficient of an overdamped nanocaterpillar, under a range of modeling assumptions. We demonstrate that the effective diffusion coefficient, which depends on the microscopic parameters governing the legs, can be orders of magnitude smaller than the background diffusion coefficient. Furthermore it varies rapidly with temperature, and reproduces the striking variations seen in existing data and our own measurements of the diffusion of DNA-coated colloids. Our model gives insight into the mechanism of motion, and allows us to ask: when does a nanocaterpillar prefer to move by sliding, where one leg is always linked to the surface, and when does it prefer to move by hopping, which requires all legs to unbind simultaneously? We compare a range of systems (viruses, molecular motors, white blood cells, protein cargos in the nuclear pore complex, bacteria such as Escherichia coli, and DNA-coated colloids) and present guidelines to control the mode of motion for materials design.

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

On-the-fly generation of differential resonance scattering probability distribution functions for Monte Carlo codes

Current Monte Carlo codes use one of three models to model neutron scattering in the epithermal energy range: (1) the asymptotic scattering model, (2) the free gas scattering model, or (3) the S({alpha},{beta}) model, depending on the neutron energy and the specific Monte Carlo code. The free gas scattering 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 using the free gas scattering model in the vicinity of the resonances in the lower epithermal range can under-predict resonance absorption due to the up-scattering phenomenon. Existing methods all involve performing the collision analysis in the center-of-mass frame, followed by a conversion back to the laboratory frame. In this paper, we will present a new sampling methodology that (1) accounts for the energy-dependent scattering cross sections in the collision analysis and (2) acts 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 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 these methods showed very close comparison to results using the reference Doppler-broadened rejection correction (DBRC) scheme. (authors)

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2013

Computational Mechanics for Heterogeneous Materials

The subject of this work is the development of models for the numerical simulation of matter, momentum, and energy balance in heterogeneous materials. These are materials that consist of multiple phases or species or that are structured on some (perhaps many) scale(s). By computational mechanics we mean to refer generally to the standard type of modeling that is done at the level of macroscopic balance laws (mass, momentum, energy). We will refer to the flow or flux of these quantities in a generalized sense as transport. At issue here are the forms of the governing equations in these complex materials which are potentially strongly inhomogeneous below some correlation length scale and are yet homogeneous on larger length scales. The question then becomes one of how to model this behavior and what are the proper multi-scale equations to capture the transport mechanisms across scales. To address this we look to the area of generalized stochastic process that underlie the transport processes in homogeneous materials. The archetypal example being the relationship between a random walk or Brownian motion stochastic processes and the associated Fokker-Planck or diffusion equation. Here we are interested in how this classical setting changes when inhomogeneities or correlations in structure are introduced into the problem. Aspects of non-classical behavior need to be addressed, such as non-Fickian behavior of the mean-squared-displacement (MSD) and non-Gaussian behavior of the underlying probability distribution of jumps. We present an experimental technique and apparatus built to investigate some of these issues. We also discuss diffusive processes in inhomogeneous systems, and the role of the chemical potential in diffusion of hard spheres is considered. Also, the relevance to liquid metal solutions is considered. Finally we present an example of how inhomogeneities in material microstructure introduce fluctuations at the meso-scale for a thermal conduction problem. These fluctuations due to random microstructures also provide a means of characterizing the aleatory uncertainty in material properties at the mesoscale.

36 MATERIALS SCIENCE↗
NASA NTRS1976

The interplanetary magnetic field and energetic electrons from Jupiter

The reported four-month 'season' for observation of Jovian electron events at AU is shown to be a plausible consequence of field-line random walk in the solar wind. Values of the field-line random-walk coefficient deduced earlier in connection with observations of solar protons, photospheric supergranulation, and the interplanetary magnetic-field power spectrum, lead to an approximate value for the spread of Jovian electrons at 1 AU, which corresponds to an approximate four-month duration of events. This model is consistent with a point source at Jupiter.

Jokipii, J. R.↗
DOE OSTI2014

A new Green's function Monte Carlo algorithm for the solution of the two-dimensional nonlinear Poisson–Boltzmann equation: Application to the modeling of the communication breakdown problem in space vehicles during re-entry

The objective of this paper is the exposition of a recently-developed, novel Green's function Monte Carlo (GFMC) algorithm for the solution of nonlinear partial differential equations and its application to the modeling of the plasma sheath region around a cylindrical conducting object, carrying a potential and moving at low speeds through an otherwise neutral medium. The plasma sheath is modeled in equilibrium through the GFMC solution of the nonlinear Poisson–Boltzmann (NPB) equation. The traditional Monte Carlo based approaches for the solution of nonlinear equations are iterative in nature, involving branching stochastic processes which are used to calculate linear functionals of the solution of nonlinear integral equations. Over the last several years, one of the authors of this paper, K. Chatterjee has been developing a philosophically-different approach, where the linearization of the equation of interest is not required and hence there is no need for iteration and the simulation of branching processes. Instead, an approximate expression for the Green's function is obtained using perturbation theory, which is used to formulate the random walk equations within the problem sub-domains where the random walker makes its walks. However, as a trade-off, the dimensions of these sub-domains have to be restricted by the limitations imposed by perturbation theory. The greatest advantage of this approach is the ease and simplicity of parallelization stemming from the lack of the need for iteration, as a result of which the parallelization procedure is identical to the parallelization procedure for the GFMC solution of a linear problem. The application area of interest is in the modeling of the communication breakdown problem during a space vehicle's re-entry into the atmosphere. However, additional application areas are being explored in the modeling of electromagnetic propagation through the atmosphere/ionosphere in UHF/GPS applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2012

THE BUILDUP OF A SCALE-FREE PHOTOSPHERIC MAGNETIC NETWORK

We use a global Monte Carlo simulation of the formation of the solar photospheric magnetic network to investigate the origin of the scale invariance characterizing magnetic flux concentrations visible on high-resolution magnetograms. The simulations include spatially and temporally homogeneous injection of small-scale magnetic elements over the whole photosphere, as well as localized episodic injection associated with the emergence and decay of active regions. Network elements form in response to cumulative pairwise aggregation or cancellation of magnetic elements, undergoing a random walk on the sphere and advected on large spatial scales by differential rotation and a poleward meridional flow. The resulting size distribution of simulated network elements is in very good agreement with observational inferences. We find that the fractal index and size distribution of network elements are determined primarily by these post-emergence surface mechanisms, and carry little or no memory of the scales at which magnetic flux is injected in the simulation. Implications for models of dynamo action in the Sun are briefly discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1980

The random walk of a drilling laser beam

The disregistry of holes drilled with a pulse laser beam in 330-micron-thick single-crystal silicon-on-sapphire wafers is examined. The exit positions of the holes were displaced from the hole entrance positions on the opposing face of the wafer, and this random displacement increased with the number of laser pulses required. A model in which the bottom of the drill hole experiences small random displacements during each laser pulse is used to describe the experimental observations. It is shown that the average random displacement caused by each pulse is only a few percent of the hole diameter and can be reduced by using as few laser pulses as necessary while avoiding the cracking and spalling of the wafer that occur with a hole drilled with a single pulse.

Anthony, T. R.↗
DOE OSTI2012

TRANSIT AND ECLIPSE ANALYSES OF THE EXOPLANET HD 149026b USING BLISS MAPPING

The dayside of HD 149026b is near the edge of detectability by the Spitzer Space Telescope. We report on 11 secondary-eclipse events at 3.6, 4.5, 3 Multiplication-Sign 5.8, 4 Multiplication-Sign 8.0, and 2 Multiplication-Sign 16 {mu}m plus three primary-transit events at 8.0 {mu}m. The eclipse depths from jointly fit models at each wavelength are 0.040% {+-} 0.003% at 3.6 {mu}m, 0.034% {+-} 0.006% at 4.5 {mu}m, 0.044% {+-} 0.010% at 5.8 {mu}m, 0.052% {+-} 0.006% at 8.0 {mu}m, and 0.085% {+-} 0.032% at 16 {mu}m. Multiple observations at the longer wavelengths improved eclipse-depth signal-to-noise ratios by up to a factor of two and improved estimates of the planet-to-star radius ratio (R{sub p} /R{sub *} = 0.0518 {+-} 0.0006). We also identify no significant deviations from a circular orbit and, using this model, report an improved period of 2.8758916 {+-} 0.0000014 days. Chemical-equilibrium models find no indication of a temperature inversion in the dayside atmosphere of HD 149026b. Our best-fit model favors large amounts of CO and CO{sub 2}, moderate heat redistribution (f = 0.5), and a strongly enhanced metallicity. These analyses use BiLinearly-Interpolated Subpixel Sensitivity (BLISS) mapping, a new technique to model two position-dependent systematics (intrapixel variability and pixelation) by mapping the pixel surface at high resolution. BLISS mapping outperforms previous methods in both speed and goodness of fit. We also present an orthogonalization technique for linearly correlated parameters that accelerates the convergence of Markov chains that employ the Metropolis random walk sampler. The electronic supplement contains light-curve files.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1991

A walk through the planned CS building

Using the architectural plan views of our future computer science building as test objects, we have completed the first stage of a Building walkthrough system. The inputs to our system are AutoCAD files. An AutoCAD converter translates the geometrical information in these files into a format suitable for 3D rendering. Major model errors, such as incorrect polygon intersections and random face orientations, are detected and fixed automatically. Interactive viewing and editing tools are provided to view the results, to modify and clean the model and to change surface attributes. Our display system provides a simple-to-use user interface for interactive exploration of buildings. Using only the mouse buttons, the user can move inside and outside the building and change floors. Several viewing and rendering options are provided, such as restricting the viewing frustum, avoiding wall collisions, and selecting different rendering algorithms. A plan view of the current floor, with the position of the eye point and viewing direction on it, is displayed at all times. The scene illumination can be manipulated, by interactively controlling intensity values for 5 light sources.

Khorramabadi, Delnaz↗
DOE OSTI2021

Transient cycling of nitrogen, organic carbon and oxygen within the free-flowing Columbia River corridor: Linking exposure time dependent biogeochemical reactions to river stage fluctuations (Final project report)

The objective of this project was to develop next-generation techniques for representing the transformations of complex reaction networks like those impacting transient river-corridors, then apply them to develop upscaling tools. The project contributed significantly to the development of Lagrangian “mass transfer particle tracking” (MTPT) tools that are the first numerical framework that explicitly separates mixing and spreading processes; a crucial distinction for accurate representations of reaction rates across scales. The MTPT approach uses an operator splitting scheme where physical transport processes are simulated using classical random walk methods and each particle is treated as a “container” that carries with it any number of chemical species. Mixing between particles is simulated using a colocation-probability based mass transfer kernel and reactions are evaluated on each particle after mixing. The method is stable, accurate, and also allows the explicit modeling of residence time distributions of the individual containers of mass. The MTPT scheme allowed us to explore several research questions related to reactive transport in river corridors and complex reaction networks. We found that it is possible to quantify the exposure time distributions (ETDs) of reactants and that this does suggest some pathways to upscaling. However, comparable approximations can be obtained using simplified (reduced dimensionality) MTPT simulations in less time with similar uncertainty, so development of extensive ETD-based methods was not productive. The major findings are that i) the MTPT schemes are robust and highly accurate across a range of arbitrarily complex reactions, ii) the methods can be efficiently parallelized and the parallel performance characteristics of MTPT are predictable, and iii) the MTPT tools allow simultaneous tracking of residence time. The development of these research tools into efficient software packages continues. The methods are already available to the community because open-source, working examples have been included with all publications.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2011

Simulation of thermal ionization in a dense helium plasma by the Feynman path integral method

The region of equilibrium states is studied where the quantum nature of the electron component and a strong nonideality of a plasma play a key role. The problem of negative signs in the calculation of equilibrium averages a system of indistinguishable quantum particles with a spin is solved in the macroscopic limit. It is demonstrated that the calculation can be conducted up to a numerical result. The complete set of symmetrized basis wave functions is constructed based on the Young symmetry operators. The combinatorial weight coefficients of the states corresponding to different graphs of connected Feynman paths in multiparticle systems are calculated by the method of random walk over permutation classes. The kinetic energy is calculated using a viral estimator at a finite pressure in a statistical ensemble with flexible boundaries. Based on the methods developed in the paper, the computer simulation is performed for a dense helium plasma in the temperature range from 30000 to 40000 K. The equation of state, internal energy, ionization degree, and structural characteristic of the plasma are calculated in terms of spatial correlation functions. The parameters of a pseudopotential plasma model are estimated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Stochastic field-line wandering in magnetic turbulence with shear. I. Quasi-linear theory

We investigate the random walk of magnetic field lines in magnetic turbulence with shear. In the first part of the series, we develop a quasi-linear theory in order to compute the diffusion coefficient of magnetic field lines. We derive general formulas for the diffusion coefficients in the different directions of space. We like to emphasize that we expect that quasi-linear theory is only valid if the so-called Kubo number is small. We consider two turbulence models as examples, namely, a noisy slab model as well as a Gaussian decorrelation model. For both models we compute the field line diffusion coefficients and we show how they depend on the aforementioned Kubo number as well as a shear parameter. It is demonstrated that the shear effect reduces all field line diffusion coefficients.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Primary fission fragment mass yields across the chart of nuclides

In this work, we have calculated a complete set of primary fission fragment mass yields, Y(A), for heavy nuclei across the chart of nuclides, including those of particular relevance to the rapid neutron capture process (r process) of nucleosynthesis. We assume that the nuclear shape dynamics are strongly damped, which allows for a description of the fission process via Brownian shape motion across nuclear potential-energy surfaces. The macroscopic energy of the potential was obtained with the Finite-Range Liquid-Drop Model (FRLDM), while the microscopic terms were extracted from the single-particle level spectra in the fissioning system by the Strutinsky procedure for the shell energies and the BCS treatment for the pairing energies. For each nucleus considered, the fission fragment mass yield, Y(A), is obtained from 50 000 to 500 000 random walks on the appropriate potential-energy surface. The full mass and charge yield, Y(Z,A), is then calculated by invoking the Wahl systematics. With this method, we have calculated a comprehensive set of fission-fragment yields from over 3800 nuclides bounded by 80 ≤ Z ≤ 130 and A ≤ 330; these yields are provided as an ASCII formatted database in the Supplemental Material. We compare our yields to known data and discuss general trends that emerge in low-energy fission yields across the chart of nuclides.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗