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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Brownian dynamics of confined rigid bodies

We introduce numerical methods for simulating the diffusive motion of rigid bodies of arbitrary shape immersed in a viscous fluid. We parameterize the orientation of the bodies using normalized quaternions, which are numerically robust, space efficient, and easy to accumulate. We construct a system of overdamped Langevin equations in the quaternion representation that accounts for hydrodynamic effects, preserves the unit-norm constraint on the quaternion, and is time reversible with respect to the Gibbs-Boltzmann distribution at equilibrium. We introduce two schemes for temporal integration of the overdamped Langevin equations of motion, one based on the Fixman midpoint method and the other based on a random finite difference approach, both of which ensure that the correct stochastic drift term is captured in a computationally efficient way. We study several examples of rigid colloidal particles diffusing near a no-slip boundary and demonstrate the importance of the choice of tracking point on the measured translational mean square displacement (MSD). We examine the average short-time as well as the long-time quasi-two-dimensional diffusion coefficient of a rigid particle sedimented near a bottom wall due to gravity. For several particle shapes, we find a choice of tracking point that makes the MSD essentially linear with time, allowing us to estimate the long-time diffusion coefficient efficiently using a Monte Carlo method. However, in general, such a special choice of tracking point does not exist, and numerical techniques for simulating long trajectories, such as the ones we introduce here, are necessary to study diffusion on long time scales.

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

Modeling electrokinetic flows with the discrete ion stochastic continuum overdamped solvent algorithm

In this article we develop an algorithm for the efficient simulation of electrolytes in the presence of physical boundaries. In previous work the discrete ion stochastic continuum overdamped solvent (DISCOS) algorithm was derived for triply periodic domains, and was validated through ion-ion pair correlation functions and Debye-Hückel-Onsager theory for conductivity, including the Wien effect for strong electric fields. In extending this approach to include an accurate treatment of physical boundaries we must address several important issues. First, the modifications to the spreading and interpolation operators necessary to incorporate interactions of the ions with the boundary are described. Next we discuss the modifications to the electrostatic solver to handle the influence of charges near either a fixed potential or dielectric boundary. An additional short-ranged potential is also introduced to represent interaction of the ions with a solid wall. Finally, the dry diffusion term is modified to account for the reduced mobility of ions near a boundary, which introduces an additional stochastic drift correction. Herein, several validation tests are presented confirming the correct equilibrium distribution of ions in a channel. Additionally, the methodology is demonstrated using electro-osmosis and induced-charge electro-osmosis, with comparison made to theory and other numerical methods. Notably, the DISCOS approach achieves greater accuracy than a continuum electrostatic simulation method. We also examine the effect of under-resolving hydrodynamic effects using a “dry diffusion” approach, and find that considerable computational speedup can be achieved with a negligible impact on accuracy.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Characterizing acid diffusion lengths in chemically amplified resists from measurements of deprotection kinetics

The acid-catalyzed deprotection of glassy poly(4-hydroxystyrene-co-tertbutyl acrylate) films was studied with infrared absorbance spectroscopy and stochastic simulations. Experimental data were interpreted with a simple description of subdiffusive acid transport coupled to second-order acid loss. This model predicts key attributes of observed deprotection rates, such as fast reaction at short times, slow reaction at long times, and a nonlinear dependence on acid loading. Fickian diffusion is approached by increasing the post-exposure bake temperature or adding plasticizing agents to the polymer resin. These findings demonstrate that acid mobility and overall deprotection kinetics are coupled to glassy matrix dynamics. To complement the analysis of bulk kinetics, acid diffusion lengths were calculated from the anomalous transport model and compared with nanopattern line widths. The consistent scaling between experiments and simulations suggests that the anomalous diffusion model could be further developed into a predictive lithography tool.

Patil, Abhijit A.↗
DOE OSTI2020

Fundamentals of Au(111) Surface Dynamics: Coarsening of Two-Dimensional Au Islands

Au(111) surfaces play a central role in many applications, yet studies of fundamental aspects of their dynamics are limited. Thus, using scanning tunneling microscopy (STM) at 300 K, we analyze the coarsening of first-layer 2D Au islands directly on the Au(111) substrate and also of second-layer 2D Au islands. Specifically, we monitor the decay of Au first-layer islands with areas of about 100–500 nm 2 in the vicinity of larger islands or extended step edges over a period of approximately 40 h, the relevant time scale for this process. Experimentally observed behavior is captured by analytical theory for terrace-diffusion-limited decay incorporating DFT results for the Au terrace diffusion barrier and the adatom formation energy.Lastly, experimental observations of second-layer island decay are also compared with appropriate analytical theory and stochastic simulations, thereby determining the effective Ehrlich–Schwoebel barrier for Au on Au(111).

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

An exact and efficient first passage time algorithm for reaction–diffusion processes on a 2D-lattice

We present an exact and efficient algorithm for reaction–diffusion–nucleation processes on a 2D-lattice. The algorithm makes use of first passage time (FPT) to replace the computationally intensive simulation of diffusion hops in KMC by larger jumps when particles are far away from step-edges or other particles. Our approach computes exact probability distributions of jump times and target locations in a closed-form formula, based on the eigenvectors and eigenvalues of the corresponding 1D transition matrix, maintaining atomic-scale resolution of resulting shapes of deposit islands. We have applied our method to three different test cases of electrodeposition: pure diffusional aggregation for large ranges of diffusivity rates and for simulation domain sizes of up to 4096×4096 sites, the effect of diffusivity on island shapes and sizes in combination with a KMC edge diffusion, and the calculation of an exclusion zone in front of a step-edge, confirming statistical equivalence to standard KMC simulations. The algorithm achieves significant speedup compared to standard KMC for cases where particles diffuse over long distances before nucleating with other particles or being captured by larger islands.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Coupled Interfacial and Bulk Kinetics Govern the Timescales of Multiphase Ozonolysis Reactions

Chemical transformations in aerosols impact the lifetime of particle phase species, the fate of atmospheric pollutants, and both climate- and health-relevant aerosol properties. Timescales for multiphase reactions of ozone in atmospheric aqueous phases are governed by coupled kinetic processes between the gas phase, the particle interface, and its bulk, which respond dynamically to reactive consumption of O 3 . However, models of atmospheric aerosol reactivity often do not account for the coupled nature of multiphase processes. Herein, to examine these dynamics, we use new and prior experimental observations of aqueous droplet reaction kinetics, including three systems with a range of surface affinities and ozonolysis rate coefficients ( trans -aconitic acid (C 6 H 6 O 6 ), maleic acid (C 4 H 4 O 4 ), and sodium nitrite (NaNO 2 )). Using literature rate coefficients and thermodynamic properties, we constrain a simple two-compartment stochastic kinetic model which resolves the interface from the particle bulk and represents O 3 partitioning, diffusion, and reaction as a coupled kinetic system. Our kinetic model accurately predicts decay kinetics across all three systems, demonstrating that both the thermodynamic properties of O 3 and the coupled kinetic and diffusion processes are key to making accurate predictions. An enhanced concentration of adsorbed O 3 , compared to gas and bulk phases is rapidly maintained and remains constant even as O 3 is consumed by reaction. Multiphase systems dynamically seek to achieve equilibrium in response to reactive O 3 loss, but this is hampered at solute concentrations relevant to aqueous aerosol by the rate of O 3 arrival in the bulk by diffusion. As a result, bulk-phase O 3 becomes depleted from its Henry's law solubility. This bulk-phase O 3 depletion limits reaction timescales for relatively slow-reacting organic solutes with low interfacial affinity (i.e., trans-aconitic and maleic acids, with k rxn ≈ 10 3 -10 4 M -1 s -1 ), which is in contrast to fast-reacting solutes with higher surface affinity (i.e., nitrite, with k rxn ≈ 10 5 M -1 s -1 ) where surface reactions strongly impact the observed decay kinetics.

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

Irreversible reactions and diffusive escape: Stationary properties

We study three basic diffusion-controlled reaction processes—annihilation, coalescence, and aggregation. We examine the evolution starting with the most natural inhomogeneous initial configuration where a half-line is uniformly filled by particles, while the complementary half-line is empty. We show that the total number of particles that infiltrate the initially empty half-line is finite and has a stationary distribution. We determine the evolution of the average density from which we derive the average total number N of particles in the initially empty half-line; e.g. for annihilation $\langle N\rangle = \frac{3}{16}+\frac{1}{4\π}$ . For the coalescence process, we devise a procedure that in principle allows one to compute P(N), the probability to find exactly N particles in the initially empty half-line; we complete the calculations in the first non-trivial case (N = 1). As a by-product we derive the distance distribution between the two leading particles.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Quasielastic neutron scattering measurements and ab initio MD-simulations on single ion motions in molten NaF

The ionic stochastic motions in the molten alkali halide NaF are investigated by quasielastic neutron scattering and first principles molecular dynamics simulation. Quasielastic neutron scattering was employed to extract the diffusion behavior of the sodium ions in the melt. An extensive first principles based simulation on a box of up to 512 particles has been performed to complement the experimental data. From that large box, a smaller 64-particle box has then been simulated over a runtime of 60 ps. A good agreement between calculated and neutron data on the level of spectral shape has been obtained. The obtained sodium diffusion coefficients agree very well. The simulation predicts a fluorine diffusion coefficient similar to the sodium one. Applying the Nernst-Einstein equation, a remarkable large cross correlation between both ions can be deduced. The velocity cross correlations demonstrate a positive correlation between the ions over a period of 0.1 ps. That strong correlation is evidence that the unlike ions do not move completely statistically independent and have a strong association over a short period of time.

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

Ultrafast vibrational spectroscopy (2D-IR) of CO{sub 2} in ionic liquids: Carbon capture from carbon dioxide’s point of view

The CO{sub 2}ν{sub 3} asymmetric stretching mode is established as a vibrational chromophore for ultrafast two-dimensional infrared (2D-IR) spectroscopic studies of local structure and dynamics in ionic liquids, which are of interest for carbon capture applications. CO{sub 2} is dissolved in a series of 1-butyl-3-methylimidazolium-based ionic liquids ([C{sub 4}C{sub 1}im][X], where [X]{sup −} is the anion from the series hexafluorophosphate (PF{sub 6}{sup −}), tetrafluoroborate (BF{sub 4}{sup −}), bis-(trifluoromethyl)sulfonylimide (Tf{sub 2}N{sup −}), triflate (TfO{sup −}), trifluoroacetate (TFA{sup −}), dicyanamide (DCA{sup −}), and thiocyanate (SCN{sup −})). In the ionic liquids studied, the ν{sub 3} center frequency is sensitive to the local solvation environment and reports on the timescales for local structural relaxation. Density functional theory calculations predict charge transfer from the anion to the CO{sub 2} and from CO{sub 2} to the cation. The charge transfer drives geometrical distortion of CO{sub 2}, which in turn changes the ν{sub 3} frequency. The observed structural relaxation timescales vary by up to an order of magnitude between ionic liquids. Shoulders in the 2D-IR spectra arise from anharmonic coupling of the ν{sub 2} and ν{sub 3} normal modes of CO{sub 2}. Thermal fluctuations in the ν{sub 2} population stochastically modulate the ν{sub 3} frequency and generate dynamic cross-peaks. These timescales are attributed to the breakup of ion cages that create a well-defined local environment for CO{sub 2}. The results suggest that the picosecond dynamics of CO{sub 2} are gated by local diffusion of anions and cations.

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

Microtearing turbulence: Magnetic braiding and disruption limit

A realistic reduced model involving a large poloidal spectrum of microtearing modes is used to probe the existence of some stochasticity of magnetic field lines. Stochasticity is shown to occur even for the low values of the magnetic perturbation δB/B devoted to magnetic turbulence that have been experimentally measured. Because the diffusion coefficient may strongly depend on the radial (or magnetic-flux) coordinate, being very low near some resonant surfaces, and because its evaluation implicitly makes a normal diffusion hypothesis, one turns to another indicator appropriate to diagnose the confinement: the mean residence time of magnetic field lines. Their computation in the microturbulence frame points to the existence of a disruption limit, namely of a critical order of magnitude of δB/B above which stochasticity is no longer benign yet, leads to a macroscopic loss of confinement in some tens to hundred of electron toroidal excursions. Since the level of magnetic turbulence δB/B has been measured to grow with the plasma electron density, this would also be a density limit.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Assimilating partial observation to enhance feedback control of stochastic dynamical systems

Here, in this paper, we present a novel methodology to tackle feedback optimal control problems in scenarios where the exact state of the controlled process is unknown. It integrates data assimilation techniques and optimal control solvers to manage partial observation of the state process, a common occurrence in practical scenarios. Traditional stochastic optimal control methods assume full state observation, which is often not feasible in real-world fluid dynamics control problems. Our approach underscores the significance of utilizing observational data to inform control policy design. Specifically, we introduce a kernel learning backward stochastic differential equation (SDE) filter to enhance data assimilation efficiency and propose a sample-wise stochastic optimization method within the stochastic maximum principle framework. We demonstrate the efficacy and accuracy of our method in the control of advection-diffusion-reaction flow problem and the Dubins airplane maneuvering problem with model uncertainty.

data driven↗
DOE OSTI2016

A stochastic model of solid state thin film deposition: Application to chalcopyrite growth

Developing high fidelity quantitative models of solid state reaction systems can be challenging, especially in deposition systems where, in addition to the multiple competing processes occurring simultaneously, the solid interacts with its atmosphere. In this work, we develop a model for the growth of a thin solid film where species from the atmosphere adsorb, diffuse, and react with the film. The model is mesoscale and describes an entire film with thickness on the order of microns. Because it is stochastic, the model allows us to examine inhomogeneities and agglomerations that would be impossible to characterize with deterministic methods. We demonstrate the modeling approach with the example of chalcopyrite Cu(InGa)(SeS){sub 2} thin film growth via precursor reaction, which is a common industrial method for fabricating thin film photovoltaic modules. The model is used to understand how and why through-film variation in the composition of Cu(InGa)(SeS){sub 2} thin films arises and persists. We believe that the model will be valuable as an effective quantitative description of many other materials systems used in semiconductors, energy storage, and other fast-growing industries.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2013

PROBING DYNAMICS OF ELECTRON ACCELERATION WITH RADIO AND X-RAY SPECTROSCOPY, IMAGING, AND TIMING IN THE 2002 APRIL 11 SOLAR FLARE

Based on detailed analysis of radio and X-ray observations of a flare on 2002 April 11 augmented by realistic three-dimensional modeling, we have identified a radio emission component produced directly at the flare acceleration region. This acceleration region radio component has distinctly different (1) spectrum, (2) light curves, (3) spatial location, and, thus, (4) physical parameters from those of the separately identified trapped or precipitating electron components. To derive evolution of physical parameters of the radio sources we apply forward fitting of the radio spectrum time sequence with the gyrosynchrotron source function with five to six free parameters. At the stage when the contribution from the acceleration region dominates the radio spectrum, the X-ray- and radio-derived electron energy spectral indices agree well with each other. During this time the maximum energy of the accelerated electron spectrum displays a monotonic increase with time from {approx}300 keV to {approx}2 MeV over roughly one minute duration indicative of an acceleration process in the form of growth of the power-law tail; the fast electron residence time in the acceleration region is about 2-4 s, which is much longer than the time of flight and so requires a strong diffusion mode there to inhibit free-streaming propagation. The acceleration region has a relatively strong magnetic field, B {approx} 120 G, and a low thermal density, n{sub e} {approx}< 2 Multiplication-Sign 10{sup 9} cm{sup -3}. These acceleration region properties are consistent with a stochastic acceleration mechanism.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

On the anisotropic advection-diffusion equation with time dependent coefficients

The advection-diffusion equation with time dependent velocity and anisotropic time dependent diffusion tensor is examined in regard to its non-classical transport features and to the use of a non-orthogonal coordinate system. Although this equation appears in diverse physical problems, particularly in particle transport in stochastic velocity fields and in underground porous media, a detailed analysis of its solutions is lacking. In order to study the effects of the time-dependent coefficients and the anisotropic diffusion on transport, we solve analytically the equation for an initial Dirac delta pulse. Here, we discuss the solutions to three cases: one based on power-law correlation functions where the pulse diffuses faster than the classical rate ~t, a second case specically designed to display slower rate of diffusion than the classical one, and a third case to describe hydrodynamic dispersion in porous media

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Generative diffusion model surrogates for mechanistic agent-based biological models

Mechanistic, multicellular, agent-based models are commonly used to investigate tissue, organ, and organism-scale biology at single-cell resolution. The Cellular-Potts Model (CPM) is a powerful and popular framework for developing and interrogating these models. CPMs become computationally expensive at large space- and time- scales making application and investigation of developed models difficult. Surrogate models may allow for the accelerated evaluation of CPMs of complex biological systems. However, the stochastic nature of these models means each set of parameters may give rise to different model configurations, complicating surrogate model development. In this work, we leverage denoising diffusion probabilistic models (DDPMs) to train a generative AI surrogate of a CPM used to investigate in vitro vasculogenesis. We describe the use of an image classifier to learn the characteristics that define unique areas of a 2-dimensional parameter space. We then apply this classifier to aid in surrogate model selection and verification. Our CPM model surrogate generates model configurations 20,000 timesteps ahead of a reference configuration and demonstrates approximately a 22x reduction in computational time as compared to native code execution. Our work represents a step towards the implementation of DDPMs to develop digital twins of stochastic biological systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Cosmic-Ray Transport in Heliospheric Magnetic Structures. II. Modeling Particle Transport through Corotating Interaction Regions

The transport of cosmic rays (CRs) in the heliosphere is determined by the properties of the solar wind plasma. The heliospheric plasma environment has been probed by spacecraft for decades and provides a unique opportunity for testing transport theories. Of particular interest for the three-dimensional (3D) heliospheric CR transport are structures such as corotating interaction regions (CIRs), which, due to the enhancement of the magnetic field strength and magnetic fluctuations within and due to the associated shocks as well as stream interfaces, do influence the CR diffusion and drift. In a three-fold series of papers, we investigate these effects by modeling inner-heliospheric solar wind conditions with the numerical magnetohydrodynamic (MHD) framework Cronos (Wiengarten et al., referred as Paper I), and the results serve as input to a transport code employing a stochastic differential equation approach (this paper). While, in Paper I, we presented results from 3D simulations with Cronos, the MHD output is now taken as an input to the CR transport modeling. We discuss the diffusion and drift behavior of Galactic cosmic rays using the example of different theories, and study the effects of CIRs on these transport processes. In particular, we point out the wide range of possible particle fluxes at a given point in space resulting from these different theories. The restriction of this variety by fitting the numerical results to spacecraft data will be the subject of the third paper of this series.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Large Deviations and Importance Sampling for Systems of Slow-Fast Motion

In this paper we develop the large deviations principle and a rigorous mathematical framework for asymptotically efficient importance sampling schemes for general, fully dependent systems of stochastic differential equations of slow and fast motion with small noise in the slow component. We assume periodicity with respect to the fast component. Depending on the interaction of the fast scale with the smallness of the noise, we get different behavior. We examine how one range of interaction differs from the other one both for the large deviations and for the importance sampling. We use the large deviations results to identify asymptotically optimal importance sampling schemes in each case. Standard Monte Carlo schemes perform poorly in the small noise limit. In the presence of multiscale aspects one faces additional difficulties and straightforward adaptation of importance sampling schemes for standard small noise diffusions will not produce efficient schemes. It turns out that one has to consider the so called cell problem from the homogenization theory for Hamilton-Jacobi-Bellman equations in order to guarantee asymptotic optimality. We use stochastic control arguments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Communication: Distinguishing between short-time non-Fickian diffusion and long-time Fickian diffusion for a random walk on a crowded lattice

The motion of cells and molecules through biological environments is often hindered by the presence of other cells and molecules. A common approach to modeling this kind of hindered transport is to examine the mean squared displacement (MSD) of a motile tracer particle in a lattice-based stochastic random walk in which some lattice sites are occupied by obstacles. Unfortunately, stochastic models can be computationally expensive to analyze because we must average over a large ensemble of identically prepared realizations to obtain meaningful results. To overcome this limitation we describe an exact method for analyzing a lattice-based model of the motion of an agent moving through a crowded environment. Using our approach we calculate the exact MSD of the motile agent. Our analysis confirms the existence of a transition period where, at first, the MSD does not follow a power law with time. However, after a sufficiently long period of time, the MSD increases in proportion to time. This latter phase corresponds to Fickian diffusion with a reduced diffusivity owing to the presence of the obstacles. Our main result is to provide a mathematically motivated, reproducible, and objective estimate of the amount of time required for the transport to become Fickian. Our new method to calculate this crossover time does not rely on stochastic simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗