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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

EGR Dilution and Fuel Property Effects on High-Efficiency Spark-Ignition Flames

Modern spark ignition internal combustion engines rely on fast combustion rates and high dilution to achieve high brake thermal efficiencies. To accomplish this, new engine designs have moved towards increased tumble ratios and stroke-to-bore ratios. Increased tumble ratios correlate positively with increases in turbulent kinetic energy and improved fuel and residual gas mixing, all of which favor faster and more efficient combustion. Longer stroke-to-bore ratios allow higher geometric compression ratios and use of late intake valve closing to control peak compression pressures and temperatures. The addition of dilution to improve efficiency is limited by the resulting increase in combustion instabilities manifested by cycle-to-cycle variability. A number of effects - preferential diffusion, turbulence-combustion interactions, stochastic flow patterns, laminar-turbulent flame kernel transitions, and relative length and velocity scales between flame and turbulence - are believed to be responsible for the increase in cycle-to-cycle variations, where their contributions are likely interlinked. Several studies have shown the influence of stochastic flow characteristics on the nature of combustion instabilities, such as velocity patterns on flame kernel formation and cycle-to-cycle variations in residual gas. However, few have focused on the specific effects of fuel properties. The objective of this work is to contrast the effects of dilution on propane stoichiometric combustion against gasoline. Dilution tolerance experiments were conducted in a purpose-built high stroke-to-bore ratio single cylinder engine with both gasoline and LPG. Three-dimensional full cycle computational fluid dynamics (CFD) simulations employing a level-set combustion approach and Reynolds averaged Navier-Stokes (RANS) turbulence modeling was used to qualitatively assess the changes in length and velocity scales for turbulence and the flame. The experimental results showed that LPG can tolerate higher exhaust gas recirculation (EGR) dilution under a variety of conditions. Analysis of CFD simulations showed that propane flames are likely less sensitive to influences from the flow field due to less thickening of the flame and higher effective flame speeds.

33 ADVANCED PROPULSION SYSTEMS↗
NASA NTRS1985

Cosmic ray modulation in a random anisotropic magnetic field

Inhomogeneities of the interplanetary magnetic field can be divided into small scale and large scale ones as may be required by the character of the problem of cosmic ray (CR) propagation. CR propagation in stochastic magnetic fields is of diffusion character. The main contribution into the scattering of CR particles is made by their interaction with inhomogeneities of the magnetic field H which have characteristic dimensions 1 of the order of Larmor radius R=cp/eH of particle (p is the absolute value of particle momentum, e is particle charge, c is velocity of light). Scattering of particles on such inhomogeneities leads to their diffusion mostly along a magnetic field with characteristic dimensions of variation in space exceeding the mean free path.

Dorman, L. I.↗
NASA NTRS2005

Numerical Model for Cosmic Rays Species Production and Propagation in the Galaxy

In recent years, considerable progress has been made in studying the propagation and origin of cosmic rays, as new and more accurate data have become available. Many models have been developed to study cosmic ray interactions and propagation showed flexibility in resembling various astrophysical conditions and good agreement with observational data. However, some astrophysical problems cannot be addressed using these models, such as the stochastic nature of the cosmic rays source, small-scale structures and inhomogeneities in the interstellar gas that can affect radioactive secondary abundance in cosmic rays. We have developed a new model and a corresponding computer code that can address some of these limitations. The model depends on the expansion of the backward stochastic solution of the general diffusion transport equation (Zhang 1999) starting from an observer position to solve a group of diffusion transport equations each of which represents a particular element or isotope of cosmic ray nuclei. In this paper we are focusing on key abundance ratios such as B/C, sub-Fe/Fe, (10)Be/(9)Be, (26)Al/(27)Al, (36)Cl/(37)Cl and (54)Mn/(55)Mn, which all have well established cross sections, to evaluate our model. The effect of inhomogeneity in the interstellar medium is investigated. The contribution of certain cosmic ray nuclei to the production of other nuclei is addressed. The contribution of various galactic locations to the production of cosmic ray nuclei observed at solar system is also investigated.

Farahat, Ashraf↗
NASA NTRS1998

Ocean Turbulence V: Mesoscale Modeling in Level Coordinates. The Effect of Random Nature of Density

The main result of this paper is the derivation of a new expression for the tracer subgrid term in level coordinates S(l) to be employed in O-GCM. The novel feature is the proper account of the random nature of the density field which strongly affects the transformation from isopycnal to level coordinates of the variables of interest, velocity and tracer fields, their correlation functions and ultimately the subgrid terms. In deriving our result we made use of measured properties of vertical ocean turbulence. The major new results are: 1) the new subgrid expression is different from that of the heuristic GM model, 2) u++(tracer)=1/2u+(thickness), where u++ and u+ are the tracer and thickness bolus velocities. In previous models, u++ = u+, 2) the subgrid for a tracer tau is not the same as that for the density rho even when one accounts for the obvious absence of a diffusion term in the latter. The difference stems from a new treatment of the stochastic nature of the density, 3) the mesoscale diffusivity enters both locally and non-locally, as the integral over all z's from the bottom of the ocean to the level z.

Canuto, V. M.↗
DOE OSTI2019

Droplet size distributions in turbulent clouds: experimental evaluation of theoretical distributions

Precipitation efficiency and optical properties of clouds, both central to determining Earth's weather and climate, depend on the size distribution of cloud particles. Here in this paper, theoretical expressions for cloud droplet size distribution shape are evaluated using measurements from controlled experiments in a convective-cloud chamber. The experiments are a unique opportunity to constrain theory because they are in steady-state and because the initial and boundary conditions are well characterized compared to typical atmospheric measurements. Three theoretical distributions obtained from a Langevin drift-diffusion approach to cloud formation via stochastic condensation are tested: (a) stochastic condensation with a constant removal time-scale; (b) stochastic condensation with a size-dependent removal time-scale; (c) droplet growth in a fixed supersaturation condition and with size-dependent removal. In addition, a similar Weibull distribution that can be obtained from the drift-diffusion approach, as well as from mechanism-independent probabilistic arguments (e.g., maximum entropy), is tested as a fourth hypothesis. Statistical techniques such as the χ2 test, sum of squared errors of prediction, and residual analysis are employed to judge relative success or failure of the theoretical distributions to describe the experimental data. An extensive set of cloud droplet size distributions are measured under different aerosol injection rates. Five different aerosol injection rates are run both for size-selected aerosol particles, and six aerosol injection rates are run for broad-distribution, polydisperse aerosol particles. In relative comparison, the most favourable comparison to the measurements is the expression for stochastic condensation with size-dependent droplet removal rate. However, even this optimal distribution breaks down for broad aerosol size distributions, primarily due to deviations from the measured large-droplet tail. Lastly, a possible explanation for the deviation is the Ostwald ripening effect coupled with deactivation/activation in polluted cloud conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Anharmonic lattice dynamics and superionic transition in AgCrSe 2

Intrinsically low lattice thermal conductivity ( κ l a t ) in superionic conductors is of great interest for energy conversion applications in thermoelectrics. Yet, the complex atomic dynamics leading to superionicity and ultralow thermal conductivity remain poorly understood. Here, we report a comprehensive study of the lattice dynamics and superionic diffusion in A g C r S e 2 from energy- and momentum-resolved neutron and X-ray scattering techniques, combined with first-principles calculations. Our results settle unresolved questions about the lattice dynamics and thermal conduction mechanism in A g C r S e 2 . We find that the heat-carrying long-wavelength transverse acoustic (TA) phonons coexist with the ultrafast diffusion of Ag ions in the superionic phase, while the short-wavelength nondispersive TA phonons break down. Strong scattering of phonon quasiparticles by anharmonicity and Ag disorder are the origin of intrinsically low κ l a t . The breakdown of short-wavelength TA phonons is directly related to the Ag diffusion, with the vibrational spectral weight associated to Ag oscillations evolving into stochastic decaying fluctuations. Furthermore, the origin of fast ionic diffusion is shown to arise from extended flat basins in the energy landscape and collective hopping behavior facilitated by strong repulsion between Ag ions. These results provide fundamental insights into the complex atomic dynamics of superionic conductors.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Novel type of chimera spiral waves arising from decoupling of a diffusible component

Spiral waves composed of coherent traveling waves surrounding a core containing stochastically distributed stationary areas are found in numerical simulations of a three-variable reaction-diffusion system with one diffusible species. In the spiral core, diffusion of this component (w) mediates transitions between dynamic states of the subsystem formed by the other two components, whose dynamics is more rapid than that of w. Diffusive coupling between adjacent sites can be “on” or “off” depending on the subsystem state. The incoherent structures in the spiral core are produced by this decoupling of the slow diffusive component from the fast non-diffusing subsystem. The phase diagram reveals that the region of incoherent behavior in chimera spirals grows drastically, leading to modulation and breakup of the spirals, in the transition zones between 1{sup n-1} and 1{sup n} local mixed-mode oscillations.

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

Path optimization with limited sensing ability

We propose a computational strategy to find the optimal path for a mobile sensor with limited coverage to traverse a cluttered region. The goal is to find one of the shortest feasible paths to achieve the complete scan of the environment. We pose the problem in the level set framework, and first consider a related question of placing multiple stationary sensors to obtain the full surveillance of the environment. By connecting the stationary locations using the nearest neighbor strategy, we form the initial guess for the path planning problem of the mobile sensor. Then the path is optimized by reducing its length, via solving a system of ordinary differential equations (ODEs), while maintaining the complete scan of the environment. Furthermore, we use intermittent diffusion, which converts the ODEs into stochastic differential equations (SDEs), to find an optimal path whose length is globally minimal. To improve the computation efficiency, we introduce two techniques, one to remove redundant connecting points to reduce the dimension of the system, and the other to deal with the entangled path so the solution can escape the local traps. Numerical examples are shown to illustrate the effectiveness of the proposed method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

An alternative realization of the exact non-Markovian stochastic Schrödinger equation

Based on the path integral approach, we derive a new realization of the exact non-Markovian stochastic Schrödinger equation (SSE). The main difference from the previous non-Markovian quantum state diffusion (NMQSD) method is that the complex Gaussian stochastic process used for the forward propagation of the wave function is correlated, which may be used to reduce the amplitude of the non-Markovian memory term at high temperatures. The new SSE is then written into the recently developed hierarchy of pure states scheme, in a form that is more closely related to the hierarchical equation of motion approach. Numerical simulations are then performed to demonstrate the efficiency of the new method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1974

A relation between cosmic-ray fluctuations, gradient, and diffusion coefficient

The motion of charged particles in a stochastic magnetic field is considered via a generalized quasi-linear expansion of Liouville's equation. The result is an equation relating cosmic-ray scintillations to particle gradients and to magnetic-field fluctuations (or diffusion coefficient). The resulting theory may be regarded as an example of a fluctuation-dissipation phenomenon, in which the diffusion coefficient plays the role of the dissipative parameter. The resonant interaction between particles and the random interplanetary magnetic field is considered explicitly, and it is shown that observed scintillations of high-energy (about 1 GeV) cosmic rays may be reasonably explained by the model.

Owens, A. J.↗
DOE OSTI2023

Constraining the Reaction Rate of Criegee Intermediates with Carboxylic Acids during the Multiphase Ozonolysis of Aerosolized Alkenes

Criegee intermediates (CIs) are central to the ozonolysis of unsaturated organic compounds, controlling the formation of small fragmentation reaction products and higher molecular weight oligomeric compounds through competing unimolecular and bimolecular pathways. In particular, the reaction of thermalized CI with carboxylic acids (RCOOH) to produce oligomeric α-acyloxyalkyl hydroperoxides (AAHPs) is of interest to the chemical transformation of organic compounds in the atmosphere. In the gas phase, the CI + RCOOH reaction proceeds at supercollisional rates, consistent with a barrierless reaction mechanism. However, the rate of this reaction in condensed phases remains uncertain. Here, we report the results of aerosol flow tube studies of the ozonolysis of submicron organic aerosol containing mixtures of an internal n-alkene (Z-9-tricosene, tri) and a saturated acid (2-hexyldecanoic acid, HDA). The decay of the reactants and formation of reaction products are monitored using atmospheric pressure chemical ionization mass spectrometry (APCI–MS). A 6-fold decrease in the reactive uptake of ozone is observed with increasing acid concentration, consistent with a previous study using alcohol additives. The decay of HDA with ozone exposure shows that HDA scavenges CIs at rates comparable to competing unimolecular pathways involving CIs. The decay kinetics of HDA are analyzed using a simple kinetic model in order to constrain the ratio of the CI + RCOOH and unimolecular CI rate constants. Empirical fitting to this model, combined with theoretical estimates of the unimolecular loss rate, yields a CI + RCOOH rate constant of 1.85 ± 0.27 × 10 –19 cm 3 molec –1 s –1 , which is six orders of magnitude slower than the expected diffusion limit in a liquid organic matrix. Stochastic kinetic simulations are conducted with different values of the CI + RCOOH rate constant to validate this result and show that the diffusion-limited rate is indeed much too rapid to reproduce the experimental results.

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

Exits in order: How crowding affects particle lifetimes

Diffusive processes are often represented using stochastic random walk frameworks. The amount of time taken for an individual in a random walk to intersect with an absorbing boundary is a fundamental property that is often referred to as the particle lifetime, or the first passage time. The mean lifetime of particles in a random walk model of diffusion is related to the amount of time required for the diffusive process to reach a steady state. Mathematical analysis describing the mean lifetime of particles in a standard model of diffusion without crowding is well known. However, the lifetime of agents in a random walk with crowding has received much less attention. Since many applications of diffusion in biology and biophysics include crowding effects, here we study a discrete model of diffusion that incorporates crowding. Using simulations, we show that crowding has a dramatic effect on agent lifetimes, and we derive an approximate expression for the mean agent lifetime that includes crowding effects. Our expression matches simulation results very well, and highlights the importance of crowding effects that are sometimes overlooked.

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

Development of morphogen gradient: The role of dimension and discreteness

The fundamental processes of biological development are governed by multiple signaling molecules that create non-uniform concentration profiles known as morphogen gradients. It is widely believed that the establishment of morphogen gradients is a result of complex processes that involve diffusion and degradation of locally produced signaling molecules. We developed a multi-dimensional discrete-state stochastic approach for investigating the corresponding reaction-diffusion models. It provided a full analytical description for stationary profiles and for important dynamic properties such as local accumulation times, variances, and mean first-passage times. The role of discreteness in developing of morphogen gradients is analyzed by comparing with available continuum descriptions. It is found that the continuum models prediction about multiple time scales near the source region in two-dimensional and three-dimensional systems is not supported in our analysis. Using ideas that view the degradation process as an effective potential, the effect of dimensionality on establishment of morphogen gradients is also discussed. In addition, we investigated how these reaction-diffusion processes are modified with changing the size of the source region.

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

PARTICLE ACCELERATION AT THE HELIOSPHERIC TERMINATION SHOCK WITH A STOCHASTIC SHOCK OBLIQUITY APPROACH

Observations by the plasma and magnetic field instruments on board the Voyager 2 spacecraft suggest that the termination shock is weak with a compression ratio of {approx}2. However, this is contrary to the observations of accelerated particle spectra at the termination shock, where standard diffusive shock acceleration theory predicts a compression ratio closer to {approx}2.9. Using our focused transport model, we investigate pickup proton acceleration at a stationary spherical termination shock with a moderately strong compression ratio of 2.8 to include both the subshock and precursor. We show that for the particle energies observed by the Voyager 2 Low Energy Charged Particle (LECP) instrument, pickup protons will have effective length scales of diffusion that are larger than the combined subshock and precursor termination shock structure observed. As a result, the particles will experience a total effective termination shock compression ratio that is larger than values inferred by the plasma and magnetic field instruments for the subshock and similar to the value predicted by diffusive shock acceleration theory. Furthermore, using a stochastically varying magnetic field angle, we are able to qualitatively reproduce the multiple power-law structure observed for the LECP spectra downstream of the termination shock.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Protein drift-diffusion dynamics and phase separation in curved cell membranes and dendritic spines: Hybrid discrete-continuum methods

We develop methods for investigating protein drift-diffusion dynamics in heterogeneous cell membranes and the roles played by geometry, diffusion, chemical kinetics, and phase separation. Our hybrid stochastic numerical methods combine discrete particle descriptions with continuum-level models for tracking the individual protein drift-diffusion dynamics when coupled to continuum fields. We show how our approaches can be used to investigate phenomena motivated by protein kinetics within dendritic spines. The spine geometry is hypothesized to play an important biological role regulating synaptic strength, protein kinetics, and self-assembly of clusters. We perform simulation studies for model spine geometries varying the neck size to investigate how phase-separation and protein organization is influenced by different shapes. We also show how our methods can be used to study the roles of geometry in reaction-diffusion systems including Turing instabilities. Furthermore, our methods provide general approaches for investigating protein kinetics and drift-diffusion dynamics within curved membrane structures.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2013

DETERMINATION OF STOCHASTIC ACCELERATION MODEL CHARACTERISTICS IN SOLAR FLARES

Following our recent paper, we have developed an inversion method to determine the basic characteristics of the particle acceleration mechanism directly and non-parametrically from observations under the leaky box framework. Earlier, we demonstrated this method for obtaining the energy dependences of the escape time and pitch angle scattering time. Here, by converting the Fokker-Planck equation to its integral form, we derive the energy dependences of the energy diffusion coefficient and direct acceleration rate for stochastic acceleration in terms of the accelerated and escaping particle spectra. Combining the regularized inversion method of Piana et al. and our procedure, we relate the acceleration characteristics in solar flares directly to the count visibility data from RHESSI. We determine the timescales for electron escape, pitch angle scattering, energy diffusion, and direct acceleration at the loop top acceleration region for two intense solar flares based on the regularized electron flux spectral images. The X3.9 class event shows dramatically different energy dependences for the acceleration and scattering timescales, while the M2.1 class event shows a milder difference. The discrepancy between the M2.1 class event and the stochastic acceleration model could be alleviated by a turbulence spectrum that is much steeper than the Kolmogorov-type spectrum. A likely explanation of the X3.9 class event could be that the escape of electrons from the acceleration region is not governed by a random walk process, but instead is affected by magnetic mirroring, in which the scattering time is proportional to the escape time and has an energy dependence similar to the energy diffusion time.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1994

Multivariate space - time analysis of PRE-STORM precipitation

This paper presents the methodologies and results of the multivariate modeling and two-dimensional spectral and correlation analysis of PRE-STORM rainfall gauge data. Estimated parameters of the models for the specific spatial averages clearly indicate the eastward and southeastward wave propagation of rainfall fluctuations. A relationship between the coefficients of the diffusion equation and the parameters of the stochastic model of rainfall fluctuations is derived that leads directly to the exclusive use of rainfall data to estimate advection speed (about 12 m/s) as well as other coefficients of the diffusion equation of the corresponding fields. The statistical methodology developed here can be used for confirmation of physical models by comparison of the corresponding second-moment statistics of the observed and simulated data, for generating multiple samples of any size, for solving the inverse problem of the hydrodynamic equations, and for application in some other areas of meteorological and climatological data analysis and modeling.

Polyak, Ilya↗
DOE OSTI2016

GALACTIC COSMIC-RAY INTENSITY MODULATION BY COROTATING INTERACTION REGION STREAM INTERFACES AT 1 au

We present a new model that couples galactic cosmic-ray (GCR) propagation with magnetic turbulence transport and the MHD background evolution in the heliosphere. The model is applied to the problem of the formation of corotating interaction regions (CIRs) during the last solar minimum from the period between 2007 and 2009. The numerical model simultaneously calculates the large-scale supersonic solar wind properties and its small-scale turbulent content from 0.3 au to the termination shock. Cosmic rays are then transported through the background, and thus computed, with diffusion coefficients derived from the solar wind turbulent properties, using a stochastic Parker approach. Our results demonstrate that GCR variations depend on the ratio of diffusion coefficients in the fast and slow solar winds. Stream interfaces inside the CIRs always lead to depressions of the GCR intensity. On the other hand, heliospheric current sheet (HCS) crossings do not appreciably affect GCR intensities in the model, which is consistent with the two observations under quiet solar wind conditions. Therefore, variations in diffusion coefficients associated with CIR stream interfaces are more important for GCR propagation than the drift effects of the HCS during a negative solar minimum.

79 ASTRONOMY AND ASTROPHYSICS↗