Search NASA⌕ Search

SEARCH · Search NASA

Results for “Stochastic diffusion”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Stochastic frequency fluctuation super-resolution imaging

The inherent non-linearity of intensity correlation functions can be used to spatially distinguish identical emitters beyond the diffraction limit, as achieved, for example, in super-resolution optical fluctuation imaging (SOFI). Here, we propose a complementary concept based on spectral correlation functions, termed spectral fluctuation super-resolution (SFSR) imaging. Through theoretical and computational analysis, we show that spatially resolving time-frequency correlation functions in the image plane can improve the imaging resolution by a factor of $\sqrt2$ in most cases and up to twofold for strictly two emitters. This improvement is achieved by quantifying the degree of correlation in spectral fluctuations across the spatial domain. Experimentally, SFSR can be implemented using a combination of interferometry and photon-correlation measurements. The method works for non-blinking emitters and stochastic spectral fluctuations with arbitrary temporal statistics. This suggests its utility in super-resolution microscopy of quantum emitters at low temperatures, where spectral diffusion is often more pronounced than emitter blinking.

47 OTHER INSTRUMENTATION↗
DOE OSTI2018

Summer term 2018 work with IMC radiative transport in binary stochastic media

In the quest for new and improved methods to solve thermal radiation transport in stochastic media, this summer project was to characterize radation transport in a deterministic lattice problem that represents a real, binary stochastic media. A representative geometry for a binary stochastic media consisting of a cube of a high volume fraction media containing constant-diameter non-overlappying spheres randomly distributed within was run previously. Previous thermal radiation wave problems in binary stochastic media displayed three characteristic thermal energy deposition wavefronts in the following order, front of the spheres, bulk material, and the backs of the spheres. Several physical aspects of the representative geometry for a stochastic binary media using an optically thin background with optically thick spheres were analyzed. This analysis includes wave speed propagation using a one dimensional diffusion approximation, energy fluence transmission, and material energy deposition. The approaches to analysis included representing one sphere of the high opacity material in a diffusive background of the low opacity material, as well as repeating spheres in a column con guration with reflective boundary conditions to approximate a lattice. Additionally, both atomic mix and Levermore-Pomraning homogenization schemes were compared against the representative cube geometry as well as the sphere column lattice geometry. The wave speed analysis with a one dimensional Marshak diffusion wave was largely ineffective at determining an effective opacity value to which a problem could be homogenized, likely due to low mean free path resolution within the problem. However, the sphere column lattice geometry showed promising results for approximating the energy fluence transmitted through the problem domain. One concern with the sphere column lattice geometry indicated the high forward-bias effect of fluence reaching the problem boundary due to streaming paths around the spheres in the case of small mean free path resolution of the background material.

36 MATERIALS SCIENCE↗
NASA NTRS1976

A kinetic-theory approach to turbulent chemically reacting flows

The paper examines the mathematical and physical foundations for the kinetic theory of reactive turbulent flows, discussing the differences and relation between the kinetic and averaged equations, and comparing some solutions of the kinetic equations obtained by the Green's function method with those obtained by the approximate bimodal method. The kinetic method described consists essentially in constructing the probability density functions of the chemical species on the basis of solutions of the Langevin stochastic equation for the influence of eddies on the behavior of fluid elements. When the kinetic equations are solved for the structure of the diffusion flame established in a shear layer by the bimodal method, discontinuities in gradients of the mean concentrations at the two flame edges appear. This is a consequence of the bimodal approximation of all distribution functions by two dissimilar half-Maxwellian functions, which is a very crude approximation. These discontinuities do not appear when the solutions are constructed by the Green's function method described here.

Chung, P. M.↗
DOE OSTI2016

Finite element formulation of fluctuating hydrodynamics for fluids filled with rigid particles using boundary fitted meshes

In this paper, we present a finite element implementation of fluctuating hydrodynamics with a moving boundary fitted mesh for treating the suspended particles. The thermal fluctuations are incorporated into the continuum equations using the Landau and Lifshitz approach [1]. The proposed implementation fulfills the fluctuation–dissipation theorem exactly at the discrete level. Since we restrict the equations to the creeping flow case, this takes the form of a relation between the diffusion coefficient matrix and friction matrix both at the particle and nodal level of the finite elements. Brownian motion of arbitrarily shaped particles in complex confinements can be considered within the present formulation. A multi-step time integration scheme is developed to correctly capture the drift term required in the stochastic differential equation (SDE) describing the evolution of the positions of the particles. The proposed approach is validated by simulating the Brownian motion of a sphere between two parallel plates and the motion of a spherical particle in a cylindrical cavity. The time integration algorithm and the fluctuating hydrodynamics implementation are then applied to study the diffusion and the equilibrium probability distribution of a confined circle under an external harmonic potential.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

The large deviation principle and steady-state fluctuation theorem for the entropy production rate of a stochastic process in magnetic fields

Fluctuation theorem is one of the major achievements in the field of nonequilibrium statistical mechanics during the past two decades. There exist very few results for steady-state fluctuation theorem of sample entropy production rate in terms of large deviation principle for diffusion processes due to the technical difficulties. Here we give a proof for the steady-state fluctuation theorem of a diffusion process in magnetic fields, with explicit expressions of the free energy function and rate function. The proof is based on the Karhunen-Loève expansion of complex-valued Ornstein-Uhlenbeck process.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2014

Tuning Monotonic Basin Hopping: Improving the Efficiency of Stochastic Search as Applied to Low-Thrust Trajectory Optimization

Trajectory optimization methods using MBH have become well developed during the past decade. An essential component of MBH is a controlled random search through the multi-dimensional space of possible solutions. Historically, the randomness has been generated by drawing RVs from a uniform probability distribution. Here, we investigate the generating the randomness by drawing the RVs from Cauchy and Pareto distributions, chosen because of their characteristic long tails. We demonstrate that using Cauchy distributions (as first suggested by Englander significantly improves MBH performance, and that Pareto distributions provide even greater improvements. Improved performance is defined in terms of efficiency and robustness, where efficiency is finding better solutions in less time, and robustness is efficiency that is undiminished by (a) the boundary conditions and internal constraints of the optimization problem being solved, and (b) by variations in the parameters of the probability distribution. Robustness is important for achieving performance improvements that are not problem specific. In this work we show that the performance improvements are the result of how these long-tailed distributions enable MBH to search the solution space faster and more thoroughly. In developing this explanation, we use the concepts of sub-diffusive, normally-diffusive, and super-diffusive RWs originally developed in the field of statistical physics.

mission design↗
DOE OSTI2014

Communication: Dominance of extreme statistics in a prototype many-body Brownian ratchet

Many forms of cell motility rely on Brownian ratchet mechanisms that involve multiple stochastic processes. We present a computational and theoretical study of the nonequilibrium statistical dynamics of such a many-body ratchet, in the specific form of a growing polymer gel that pushes a diffusing obstacle. We find that oft-neglected correlations among constituent filaments impact steady-state kinetics and significantly deplete the gel's density within molecular distances of its leading edge. These behaviors are captured quantitatively by a self-consistent theory for extreme fluctuations in filaments' spatial distribution.

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

Application of Scalar Filtered Density Function to Turbulent Flows Under Supercritical Condition

Abstract The scalar filtered density function (FDF) methodology is extended and employed for large eddy simulation (LES) of turbulent flows under supercritical condition. To describe real fluid behavior, the extended methodology incorporates the generalized heat and mass diffusion models along with real fluid thermodynamic relations which are derived using the cubic Peng–Robinson equation of state. These models are implemented within the stochastic differential equations comprising the scalar FDF transport. Simulations are conducted of a temporally developing mixing layer under supercritical condition and the results are assessed by comparing with data generated by direct numerical simulation (DNS) of the same layer. The consistency of the proposed FDF methodology is assessed. The LES-FDF predictions are shown to agree favorably with the DNS data and exhibit several key features pertaining to supercritical turbulent flows.

Energy & Fuels↗
NASA NTRS1987

Kinetic aspects of tail dynamics - Theory and simulation

Kinetic theories relevant to the geomagnetic tail are reviewed. The topics discussed include kinetic instabilities, simulations, and current-sheet particle acceleration. Tearing mode and reconnection theories are emphasized. Kinetic treatment is appropriate for these topics since the tail plasma is collisionless. Fluid calculations are appropriate when stochastic processes dominate and for studies where long wavelengths are important. However, fluid treatments of tearing modes and reconnection require a finite resistivity in the diffusion region. Thus, although 'anomalous resistivity' can be guessed or in some cases calculated, ideally the kinetic treatment is often to be preferred. Particle motion and acceleration in the current sheet can give rise to beam-like distributions in the plasma-sheet boundary layer. Studies of current-sheet particle motion have also been used as the basis for 'kinetic' tail equilibrium models. Furthermore, quite recently current-sheet particle motion is used directly in Coroniti's explosive tail reconnection model. The 'inertial conductivity' from the equilibrium models provides the 'dissipation' necessary for reconnection.

Speiser, T. W.↗
NASA NTRS1989

Computations of turbulent evaporating sprays

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

Aggarwal, S. K.↗
NASA NTRS1999

Ocean Turbulence: Mesoscale Modeling in Isopycnal Coordinates the role of the Spectrum of Vertical Shear - Part 4

We study the tracer subgrid term in isopycnal coordinates, S(sub I). We employ two ingredients: the experimental data on vertical spectra of ocean turbulence measured by Gargett et al.(1981) and the stochastic approach recently developed by Dukowicz and Smith (1997). Our result confirms that S(sub I) is made of two parts: an advection and a diffusion term. However, the tracer bolus velocity u** consists of two terms u** = u(sub 1) + u(sub 2) while in the GM model there is only a term related to u(sub 1) which is shown to be: u(sub 1) = k(bar-q)(sup -1)(delta)(sub rho) where bar-q is the thickness weighted average potential vorticity, a result in agreement with the recent suggestions by Treguier et al. (1997), Lee et al. (1997) and Greatbatch (1998). The second component u(sub 2) IS new. We compute it in the geostrophic approximation using the Gargett et al. data (1981) on ocean vertical turbulence. We find that u(sub 2) much greater than u(sub 1) and that u(sub 2) is orthogonal to u(sub 1).

Canuto, V. M.↗
DOE OSTI2023

Statistical Learning for Nonlinear Model Reduction from Local Simulations of Stochastic and Particle- and Agent-Based Systems

Stochastic physical systems across the sciences that have very high-dimensional state spaces, with a large number of fast degrees of freedom that force direct simulators to proceed by integration steps that are orders of magnitude smaller than events of interests (e.g., particle collisions). Examples range from molecular motion to dynamics of large populations of cells. A grand challenge in the simulation and understanding of such systems is the systematic construction of accurate, interpretable, reduced models, enabling faster simulations, revealing fundamental properties of the dynamics, and predicting phenomena of interest that the original simulator could not reached with sufficient accuracy or within a given computational budget. In this projected we developed novel statistical estimation/machine learning techniques for analyzing and building empirical reduced models for important families of high-dimensional stochastic systems, in particular: - we developed techniques for estimating interaction kernels in interacting particle- and agent-based systems, which are ubiquitous in Physics, Biology and many other sciences, given observed trajectories of the system; - we developed techniques for nonlinear model reduction for high-dimensional stochastic systems that have a small number of unknown, nonlinear slow variables, and a large number of fast modes, that are possibly of large magnitude, given observed short trajectories of the system in the form of bursts of trajectories from different initial conditions; - we developed novel techniques for estimating linear dynamical systems on graphs when both the dynamics and the underlying graph are unknown, and we have a sparse set of space-time observations; - we considered the problem of estimating an unknown nonlinear observation function of a standard process (e.g. Brownian motion), so that we can recognized if an observed dynamics is "just" a nonlinear version of a known dynamics; we also developed benchmarks for learning algorithms aimed at learning and classifying diffusion processes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

STOCHASTIC PARTICLE ACCELERATION IN TURBULENCE GENERATED BY MAGNETOROTATIONAL INSTABILITY

We investigate stochastic particle acceleration in accretion flows. It is believed that magnetorotational instability (MRI) generates turbulence inside accretion flows and that cosmic rays (CRs) are accelerated by the turbulence. We calculate equations of motion for CRs in the turbulent fields generated by MRI with the shearing box approximation and without back reaction to the field. Our results show that the CRs randomly gain or lose their energy through interaction with the turbulent fields. The CRs diffuse in the configuration space anisotropically: the diffusion coefficient in the direction of the unperturbed flow is about 20 times higher than the Bohm coefficient, while those in the other directions are only a few times higher than the Bohm. The momentum distribution is isotropic and its evolution can be described by the diffusion equation in momentum space where the diffusion coefficient is a power-law function of the CR momentum. We show that the shear acceleration works efficiently for energetic particles. We also cautiously note that in the shearing box approximation, particles that cross the simulation box many times along the radial direction undergo unphysical runaway acceleration by the Lorentz transformation, which needs to be taken into account with special care.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Collisional dependence of Alfven mode saturation in tokamaks

Saturation of \alfven modes driven unstable by a distribution of high energy particles as a function of collisionality is investigated with a guiding center code, using numerical eigenfunctions produced by linear theory and numerical high energy particle distributions. The most important resonance is found and it is shown that when the resonance domain is bounded, not allowing particles to collisionlessly escape, the saturation amplitude is given by the balance of the resonance mixing time with the time for nearby particles to collisionally diffuse across the resonance width. Saturation amplitudes are in agreement with theoretical predictions as long as the mode amplitude is not so large that it produces stochastic loss from the resonance domain.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

An agent-based approach to study the diffusion rate and the effect of policies on joint placement of photovoltaic panels and green roof under climate change uncertainty

As two of the highest trending green technologies, photovoltaic panels and green roofs are proven to be effective practices for energy generation and energy saving. The achievable impact from the widespread installation of such technologies is, however, not clearly established. This is mainly because the degree of this impact highly depends on the inherently uncertain environmental and climate factors, as well as the unknown adoption rates of these technologies, which in turn depend on different characteristics of decision makers and interactions among them. To that end, this study aims to investigate the diffusion rate of these green technologies under uncertainties caused by climate change, characteristics of adopters, and their interactions. An integrated framework is developed to capture the interplay between financial and attitudinal aspects, as well as the uncertainties due to both the stochastic nature of system parameters and the interactions among agents involving human beings. Specifically, this framework consists of a integer programming model to optimize the green roof and/or photovoltaic panel installation settings for a given building under climate change uncertainty, and an agent-based model to factor in the role of human behavior and interactions. A case study for the city of Knoxville, TN, is presented to evaluate the effects of different policies on the diffusion rate of the green technologies of interest. The results show that the affordability of green technologies and public awareness are the key drivers of the adoption of these technologies, which highlight the important role of the decision makers in impacting the diffusion rate.

14 SOLAR ENERGY↗
DOE OSTI2014

The formation of kappa-distribution accelerated electron populations in solar flares

Driven by recent RHESSI observations of confined loop-top hard X-ray sources in solar flares, we consider stochastic acceleration of electrons in the presence of Coulomb collisions. If electron escape from the acceleration region can be neglected, the electron distribution function is determined by a balance between diffusive acceleration and collisions. Such a scenario admits a stationary solution for the electron distribution function that takes the form of a kappa distribution. We show that the evolution toward this kappa distribution involves a 'wave front' propagating forward in velocity space, so that electrons of higher energy are accelerated later; the acceleration timescales with energy according to τ{sub acc} ∼ E {sup 3/2}. At sufficiently high energies escape from the finite-length acceleration region will eventually dominate. For such energies, the electron velocity distribution function is obtained by solving a time-dependent Fokker-Planck equation in the 'leaky-box' approximation. Solutions are obtained in the limit of a small escape rate from an acceleration region that can effectively be considered a thick target.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Collisional dependence of Alfvén mode saturation in tokamaks

Saturation of Alfvén modes driven unstable by a distribution of high energy particles as a function of collisionality is investigated with a guiding center code, using numerical eigenfunctions produced by linear theory and numerical high energy particle distributions. The most important resonance is found and it is shown that when the resonance domain is bounded, not allowing particles to collisionlessly escape, the saturation amplitude is given by the balance of the resonance mixing time with the time for nearby particles to collisionally diffuse across the resonance width. Finally, saturation amplitudes are in agreement with theoretical predictions as long as the mode amplitude is not so large that it produces stochastic loss from the resonance domain.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗