Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Langevin 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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Predictive Integrated Modeling of Heated Tokamak Plasmas (Final Technical Report)

The research carried out by the Lehigh University fusion group during this reporting period: (1) Self-consistent predictive core-pedestal ITER scenario modeling; (2) Microtearing instabilities and electron thermal transport in standard and spherical tokamak discharges; (3) Neural network model of the Multi-Mode anomalous transport module; (4) Anomalous transport in KSTAR, JET, EAST, Alcator C-Mod and DIII-D discharges; (5) Theory for transport in magnetized plasmas; (6) Control-oriented modeling of the q-profile dynamic response in preparation for advanced scenario control in KSTAR and in EAST; (7) Modeling of edge harmonic oscillation in DIII-D QH mode discharges; (8) Elucidating plasma dynamics in Hasegawa-Wakatani turbulence by information geometry; (9) Non-linear Langevin and Fractional Fokker–Planck equations for anomalous diffusion by Levy stable processes; and (10) Investigation of the plasma shaping effects on the H-mode pedestal structure.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1983

Termolecular ionic recombination at high ambient gas density

It is shown that the importance of diffusion relative to drift in ion-ion recombination is a decreasing function of the separation between the ions (contrary to a common supposition). Langevin's formula for the recombination coefficient is obtained very simply using a physical model which allows for both diffusion and drift.

Bates, D. R.↗
DOE OSTI2024

Shot noise and universal Fano factor as a characterization of strongly correlated metals

Shot noise measures out-of-equilibrium current fluctuations and is a powerful tool to probe the nature of current-carrying excitations in quantum systems. Recent shot-noise measurements in the heavy-fermion strange metal YbRh 2 Si 2 exhibit a strong suppression of the Fano factor ( F )—the ratio of the current noise to the average current in the dc limit. This system is representative of metals in which electron correlations are extremely strong. Here we carry out the first theoretical study on the shot noise of diffusive metals in the regime of strong correlations. A Boltzmann-Langevin equation formulation is constructed in a quasiparticle description in the presence of strong correlations. We find that F = 3 / 4 in such a correlation regime. Thus, we establish the aforementioned Fano factor as universal to Fermi liquids, and we show that the Fano factor suppression observed in experiments on YbRh 2 Si 2 necessitates a loss of the quasiparticles. Our work opens the door to systematic theoretical studies of shot noise as a means of characterizing strongly correlated metallic phases and materials. Published by the American Physical Society 2024

Wang, Yiming (ORCID:0000000277582921)↗
DOE OSTI2020

A collective elastic fluctuation mechanism for decoupling and stretched relaxation in glassy colloidal and molecular liquids

We propose a microscopic theory for the decoupling of self-diffusion and structural relaxation in glass-forming liquids within the Elastically Collective Nonlinear Langevin Equation (ECNLE) activated dynamics framework. Our central hypothesis is that the heterogeneity relevant to this problem is static fluctuations of local density on the scale of 3–4 particle diameters and how this changes local packing correlations. These fluctuations modify the degree of dynamical cage expansion that mechanistically couples intracage large amplitude hopping and longer range collective elasticity in ECNLE theory. Decoupling only emerges in the deeply supercooled regime where the strongly temperature dependent elastic barrier becomes non-negligible relative to its noncooperative local analog. The theory makes predictions for various aspects of the decoupling phenomenon, including apparent fractional power law Stokes-Einstein behavior, that appear to be consistent with experiments and simulations on hard sphere fluids and molecular liquids. Of central importance is a microscopic connection between the barrier fluctuation variance and most probable barrier height. Sensible results are also obtained for the nonexponential stretching of a generic relaxation time correlation function and its temperature evolution. Nonuniversality can arise from the relative importance of the local and collective barriers (related to fragility) and the precise magnitude of the length scale that defines the transition from local cage to elastic physics. Comparison is made with a traplike model based on a Gaussian distribution of barriers.

Xie, Shi-Jie (ORCID:0000000333174720)↗
DOE OSTI2026

Substrate-Mediated Evaporation and Stochastic Evolution of Supported Au Nanoparticles

Here, we use in situ transmission electron microscopy with automated tracking to study supported gold nanoparticles (NPs) during high-temperature vacuum annealing. The average mass loss per NP is governed by a flat, nearly size-independent substrate-mediated evaporation profile. On top of this mean shrinkage, individual NPs show significant fluctuations in apparent growth or shrinkage, and NP volume follows a random-walk-like trajectory. To rationalize both the ensemble-mean behavior and the particle-resolved variability, we develop a self-consistent theory that couples substrate-mediated evaporation to collective 2D Ostwald-type mass exchange through a shared adatom field, described in terms of a renormalized screening length and background concentration. In the experimentally relevant regime, the theory predicts an approximately size-independent mean shrinkage rate and clarifies how net mass loss suppresses classical coarsening. Superimposed on this deterministic drift, we quantify stochastic volume trajectories and capture their fluctuation spectrum with a minimal Langevin description consistent with intermittent adatom attachment and detachment events. In addition, we characterize the lateral diffusive motion of NPs, which is responsible for their coalescence. Altogether, our results highlight that stochasticity is intrinsic at the nanoscale and that predicting the evolution of supported NPs at early and intermediate times requires a unified framework combining substrate-mediated evaporation, collective mass exchange, and stochastic fluctuations.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

Anharmonic phonons, thermal expansion, and nuclear quantum effects in Zn

The anharmonic behavior of phonons and thermal expansion of hexagonal zinc were studied from 15 to 690 K by inelastic neutron scattering (INS) and ab initio simulations. Phonon spectra were measured for Q-points covering the full Brillouin zone, giving the phonon density of states (DOS) and dispersions along high-symmetry directions. The dispersions were sharp at 15 K, but diffuse intensity was observed at energies above them. The dispersions broadened with temperature T and the diffuse intensity became stronger. This diffuse intensity appeared in all INS measurements and simulations, except for classical molecular dynamics at 15 K. The temperature-dependent effective potential (TDEP) method, which included the nuclear quantum effect from zero-point vibrational dynamics, was used to calculate the free energy and thermal expansion. For T < 100 K nuclear quantum effects were important for obtaining the correct negative thermal expansion, and path integral molecular dynamics (PIMD) was particularly effective for obtaining the negative thermal expansion in the basal plane. A Heisenberg-Langevin model for interacting phonons coupled to a thermal bath was able to reproduce the shape and intensity of the diffuse spectral features.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Bringing discrete-time Langevin splitting methods into agreement with thermodynamics

In light of the recently published complete set of statistically correct Grønbech–Jensen (GJ) methods for discrete-time thermodynamics, we revise a differential operator splitting method for the Langevin equation in order to comply with the basic GJ thermodynamic sampling features, namely, the Boltzmann distribution and Einstein diffusion, in linear systems. This revision, which is based on the introduction of time scaling along with flexibility of a discrete-time velocity attenuation parameter, provides a direct link between the ABO splitting formalism and the GJ methods. This link brings about the conclusion that any GJ method has at least weak second order accuracy in the applied time step. It further helps identify a novel half-step velocity, which simultaneously produces both correct kinetic statistics and correct transport measures for any of the statistically sound GJ methods. Explicit algorithmic expressions are given for the integration of the new half-step velocity into the GJ set of methods. Finally, numerical simulations, including quantum-based molecular dynamics (QMD) using the QMD suite Los Alamos Transferable Tight-Binding for Energetics, highlight the discussed properties of the algorithms as well as exhibit the direct application of robust, time-step-independent stochastic integrators to QMD.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1983

Stochastic differential equations and turbulent dispersion

Aspects of the theory of continuous stochastic processes that seem to contribute to an understanding of turbulent dispersion are introduced and the theory and philosophy of modelling turbulent transport is emphasized. Examples of eddy diffusion examined include shear dispersion, the surface layer, and channel flow. Modeling dispersion with finite-time scale is considered including the Langevin model for homogeneous turbulence, dispersion in nonhomogeneous turbulence, and the asymptotic behavior of the Langevin model for nonhomogeneous turbulence.

Durbin, P. A.↗
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 OSTI2020

Stratification of polymer mixtures in drying droplets: Hydrodynamics and diffusion

In this work, we study the evaporation-induced stratification of a mixture of short and long polymer chains in a drying droplet using molecular simulations. We systematically investigate the effects of hydrodynamic interactions (HI) on this process by comparing hybrid simulations accounting for HI between polymers through the multiparticle collision dynamics technique with free-draining Langevin dynamics simulations neglecting the same. We find that the dried supraparticle morphologies are homogeneous when HI are included but are stratified in core–shell structures (with the short polymers forming the shell) when HI are neglected. The simulation methodology unambiguously attributes this difference to the treatment of the solvent in the two models. We rationalize the presence (or absence) of stratification by measuring phenomenological multicomponent diffusion coefficients for the polymer mixtures. The diffusion coefficients show the importance of not only solvent backflow but also HI between polymers in controlling the dried supraparticle morphology.

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

Modeling nanoparticle charge distribution in the afterglow of non-thermal plasmas and comparison with measurements

Particle charging in the afterglows of non-thermal plasmas typically take place in a non-neutral space charge environment. We model the same by incorporating particle-ion collision rate constant models, developed in prior work by analyzing particle-ion trajectories calculated using Langevin Dynamics simulations, into species transport equations for ions, electrons and charged particles in the afterglow. A scaling analysis of particle charging and additional Langevin Dynamics calculations of the particle-ion collision rate constant are presented to extend the range of applicability to ion electrostatic to thermal energy ratios of 300 and diffusive Knudsen number (that scales inversely with gas pressure) up to 2000. The developed collision rate constant models are first validated by comparing predictions of particle charge against measured values in a stationary, non-thermal DC plasma from past PK-4 campaigns published in Phys. Rev. Lett. 93(8): 085001 and Phys. Rev. E 72(1): 016406). The comparisons reveal excellent agreement within ±35% for particles of radius 0.6,1.0,1.3 μm in the gas pressure range of ~20-150 Pa. The experiments to probe particle charge distributions by Sharma et al. (J. Physics D: Appl. Phys. 53(24): 245204) are modeled using the validated particle-ion collision rate constant models and the calculated charge fractions are compared with measurements. The comparisons reveal that the ion/electron concentration and gas temperature in the afterglow critically influence the particle charge and the predictions are generally in qualitative agreement with the measurements. Along with critical assessment of the modeling assumptions, several recommendations are presented for future experimental design to probe charging in afterglows.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Combining machine-learned and empirical force fields with the parareal algorithm: application to the diffusion of atomistic defects

We numerically investigate an adaptive version of the parareal algorithm in the context of molecular dynamics. This adaptive variant has been originally introduced in [1]. We focus here on test cases of physical interest where the dynamics of the system is modelled by the Langevin equation and is simulated using the molecular dynamics software LAMMPS. In this work, the parareal algorithm uses a family of machine-learning spectral neighbor analysis potentials (SNAP) as fine, reference, potentials and embedded-atom method potentials (EAM) as coarse potentials. We consider a self-interstitial atom in a tungsten lattice and compute the average residence time of the system in metastable states. Our numerical results demonstrate significant computational gains using the adaptive parareal algorithm in comparison to a sequential integration of the Langevin dynamics. We also identify a large regime of numerical parameters for which statistical accuracy is reached without being a consequence of trajectorial accuracy.

36 MATERIALS SCIENCE↗
NASA NTRS1975

Langevin equation with stochastic damping - Possible application to critical binary fluid

We solve the familiar Langevin equation with stochastic damping to represent the motion of a Brownian particle in a fluctuating medium. A connection between the damping and the random driving forces is proposed which preserves quite generally the Einstein relation between the diffusion and mobility coefficients. We present an application to the case of a Brownian particle in a critical binary mixture.

Jasnow, D.↗
DOE OSTI2022

A Stochastic Reduced-Order Model for Statistical Microstructure Descriptors Evolution

Integrated computational materials engineering (ICME) models have been a crucial building block for modern materials development, relieving heavy reliance on experiments and significantly accelerating the materials design process. However, ICME models are also computationally expensive, particularly with respect to time integration for dynamics, which hinders the ability to study statistical ensembles and thermodynamic properties of large systems for long time scales. To alleviate the computational bottleneck, we propose to model the evolution of statistical microstructure descriptors as a continuous-time stochastic process using a non-linear Langevin equation, where the probability density function (PDF) of the statistical microstructure descriptors, which are also the quantities of interests (QoIs), is modeled by the Fokker–Planck equation. In this work, we discuss how to calibrate the drift and diffusion terms of the Fokker–Planck equation from the theoretical and computational perspectives. The calibrated Fokker–Planck equation can be used as a stochastic reduced-order model to simulate the microstructure evolution of statistical microstructure descriptors PDF. Considering statistical microstructure descriptors in the microstructure evolution as QoIs, we demonstrate our proposed methodology in three integrated computational materials engineering (ICME) models: kinetic Monte Carlo, phase field, and molecular dynamics simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Prediction and observation of intermodulation sidebands from anharmonic phonons in NaBr

A quantum Langevin model, like models from optomechanics, was developed for phonons. It predicts intermodulation phonon sidebands (IPSs) in anharmonic crystals. Ab initio calculations of anharmonic phonons in rock-salt NaBr showed these spectral features as many-body effects. Modern inelastic neutron scattering measurements on a crystal of NaBr at 300 K revealed diffuse intensity at high phonon energy from a predicted upper IPS. The transverse optical (TO) part of the new features originates from phonon intermodulation between the transverse acoustic (TA) and TO phonons. The longitudinal optical spectral features originate from three-phonon coupling between the TA modes and the TO lattice modes. The partner lower IPS proves to be an intrinsic localized mode. Interactions with the thermal bath broaden and redistribute the spectral weight of the IPS pair. These sidebands are a probe of the anharmonicity and quantum noise of phonons in NaBr and suggest novel interactions between photons and phonons.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Formation of transfermium elements in reactions with Pb 208

Within the Langevin framework, we investigate the dynamics of the fusion process for production of transfermium elements in reactions of Ca 48 , Ti 50 , Cr 54 , and Fe 58 with Pb 208 . After the reacting nuclei have made contact, the early dynamical stage is dominated by the dissipation of the initial radial kinetic energy, while the subsequent shape evolution is diffusive. The probability for surmounting the inner barrier and forming a compound system is obtained by simulating the evolution as a Metropolis random walk in a five-dimensional potential-energy landscape. Good agreement with the available data is obtained, especially for the maximal formation probability. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Heavy flavor quenching and flow: the roles of initial condition, pre-equilibrium evolution, and in-medium interaction

Within an advanced Langevin-hydrodynamics framework coupled to a hybrid fragmentation-coalescence hadronization model, we study heavy flavor quenching and flow in relativistic heavy-ion collisions. We investigate how the initial heavy quark spectrum, the in-medium energy loss and hadronization mechanisms of heavy quarks, the evolution profile of the pre-equilibrium stage, the medium flow, and the temperature dependence of heavy quark diffusion coefficients influence the suppression and elliptic flow of heavy mesons at the RHIC and the LHC. Our results show that the different modeling of initial conditions, pre-equilibrium evolution, and in-medium interactions can individually yield uncertainties of approximately 10-40% in D meson suppression and flow at a low transverse momentum. Furthermore, we also find that proper combinations of collisional versus radiative energy loss, coalescence versus fragmentation in hadronization, and the inclusion of medium flow are the most important factors for describing the suppression and elliptic flow of heavy mesons.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Coarse-grained explicit-solvent molecular dynamics simulations of semidilute unentangled polyelectrolyte solutions

In this study, we present results from explicit-solvent coarse-grained molecular dynamics (MD) simulations of fully charged, salt-free, and unentangled polyelectrolytes in semidilute solutions. The inclusion of a polar solvent in the model allows for a more physical representation of these solutions at concentrations, where the assumptions of a continuum dielectric medium and screened hydrodynamics break down. The collective dynamic structure factor of polyelectrolytes, S(q, t), showed that at q > q*, where q* = 2π/ξ is the polyelectrolyte peak in the structure factor S(q) and ξ is the correlation length, the relaxation time obtained from fits to stretched exponential was $\tau$ KWW ~ q -3 , which describes unscreened Zimm-like dynamics. This is in contrast to implicit-solvent simulations using a Langevin thermostat where $\tau$ KWW ~ q -2 . At q < q*, a crossover region was observed that eventually transitions to another inflection point $\tau$ KWW ~ q -2 at length scales larger than ξ for both implicit- and explicit-solvent simulations. The simulation results were also compared to scaling predictions for correlation length, ξ ~ c$-½\atop{p}$, specific viscosity, η sp ~ c$½\atop{p}$, and diffusion coefficient, D ~ c$0\atop{p}$, where c p is the polyelectrolyte concentration. The scaling prediction for ξ holds; however, deviations from the predictions for η sp and D were observed for systems at higher c p , which are in qualitative agreements with recent experimental results. This study highlights the importance of explicit-solvent effects in molecular dynamics simulations, particularly in semidilute solutions, for a better understanding of polyelectrolyte solution behavior.

36 MATERIALS SCIENCE↗