Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusivity”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Neural network kinetics for exploring diffusion multiplicity and chemical ordering in compositionally complex materials

Diffusion involving atom transport from one location to another governs many important processes and behaviors such as precipitation and phase nucleation. The inherent chemical complexity in compositionally complex materials poses challenges for modeling atomic diffusion and the resulting formation of chemically ordered structures. Here, we introduce a neural network kinetics (NNK) scheme that predicts and simulates diffusion-induced chemical and structural evolution in complex concentrated chemical environments. The framework is grounded on efficient on-lattice structure and chemistry representation combined with artificial neural networks, enabling precise prediction of all path-dependent migration barriers and individual atom jumps. To demonstrate the method, we study the temperature-dependent local chemical ordering in a refractory NbMoTa alloy and reveal a critical temperature at which the B2 order reaches a maximum. The atomic jump randomness map exhibits the highest diffusion heterogeneity (multiplicity) in the vicinity of this characteristic temperature, which is closely related to chemical ordering and B2 structure formation. The scalable NNK framework provides a promising new avenue to exploring diffusion-related properties in the vast compositional space within which extraordinary properties are hidden.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Role of Sr doping and external strain on relieving bottleneck of oxygen diffusion in La 2- x Sr x CuO 4-δ

In many complex oxides, the oxygen vacancy formation is a promising route to modify the material properties such as a superconductivity and an oxygen diffusivity. Cation substitutions and external strain have been utilized to control the concentration and diffusion of oxygen vacancies, but the mechanisms behind the controls are not fully understood. Using first-principles calculations, we find how Sr doping and external strain greatly enhances the diffusivity of oxygen vacancies in La 2-x Sr x CuO 4-δ (LSCO) in the atomic level. In hole-doped case (2x > δ), the formation energy of an apical vacancy in the LaO layer is larger than its equatorial counterpart by 0.2 eV that the bottleneck of diffusion process is for oxygen vacancies to escape equatorial sites. Such an energy difference can be reduced and even reversed by either small strain (< 1.5%) or short-range attraction between Sr and oxygen vacancy, and in turn, the oxygen diffusivity is greatly enhanced. For fully compensated hole case (2x ≦ δ), the formation energy of an apical vacancy becomes too high that most oxygen vacancies cannot move but would be trapped at equatorial sites. From our electronic structure analysis, we found that the contrasting change in the formation energy by Sr doping and external strain is originated from the different localization natures of electron carrier from both types of oxygen vacancies.

25 ENERGY STORAGE↗
DOE OSTI2021

Fast Na diffusion and anharmonic phonon dynamics in superionic Na 3 PS 4

The design of new solid electrolytes (SEs) hinges on identifying and tuning relevant descriptors. Phonons describe the atomic dynamics in crystalline materials and provide a basis to encode possible minimum energy pathways for ion migration, but anharmonic effects can be large in SEs. Identifying and controlling the pertinent phonon modes coupled most strongly with ionic conductivity, and assessing the role of anharmonicity, could therefore pave the way for discovering and designing new SEs via phonon engineering. In this work, we investigate phonons in Na 3 PS 4 and their coupling to fast Na diffusion, using a combination of neutron scattering, ab initio molecular dynamics (AIMD), and extended molecular dynamics based on machine-learned potentials. We identify that anharmonic soft-modes at the Brillouin zone boundary of the anharmonically stabilized cubic phase constitute key phonon modes that control the Na diffusion process in Na 3 PS 4 . We demonstrate how these strongly anharmonic phonon modes enable Na-ions to hop along the minimum energy pathways. Further, the quasi-elastic neutron scattering (QENS) measurements, supplemented with ab initio and machine-learned molecular dynamics simulations, probe the Na diffusivity and diffusion mechanism. These results offer detailed microscopic insights into the dynamic mechanism of fast Na diffusion and provide an avenue to search for further Na solid electrolytes.

42 ENGINEERING↗
DOE OSTI2021

Stoichiometric tuning of lattice flexibility and Na diffusion in NaAlSiO 4 : quasielastic neutron scattering experiment and ab initio molecular dynamics simulations

We have performed quasielastic neutron scattering (QENS) experiments up to 1243 K and ab initio molecular dynamics (AIMD) simulations to investigate the Na diffusion in various phases of NaAlSiO4 (NASO), namely, low-carnegieite (L-NASO; trigonal), high-carnegieite (H-NASO; cubic) and nepheline (N-NASO; hexagonal) phases. The QENS measurements reveal Na ions localized diffusion behavior in L-NASO and N-NASO, but long-range diffusion behavior in H-NASO. Furthermore, the AIMD simulation supplemented the QENS measurements and showed that excess Na ions in H-NASO enhance the host network flexibility and activate the AlO 4 /SiO 4 tetrahedra rotational modes. These framework modes enable the long-range diffusion of Na across a pathway of interstitial sites. The simulations also show Na diffusion in Na-deficient N-NASO through vacant Na sites along the hexagonal c-axis.

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

Optical/X-ray/radio view of Abell 1213: A galaxy cluster with anomalous diffuse radio emission

Context. Abell 1213, a low-richness galaxy system, is known to host an anomalous radio halo detected in data of the Very Large Array (VLA). It is an outlier with regard to the relation between the radio halo power and the X-ray luminosity of the parent clusters. Aims. Our aim is to analyze the cluster in the optical, X-ray, and radio bands to characterize the environment of its diffuse radio emission and to shed new light on its nature. Methods. We used optical data from the Sloan Digital Sky Survey to study the internal dynamics of the cluster. We also analyzed archival XMM-Newton X-ray data to unveil the properties of its hot intracluster medium. Finally, we used recent data from the LOw Frequency ARray (LOFAR) at 144 MHz, together with VLA data at 1.4 GHz, to study the spectral behavior of the diffuse radio source. Results. Both our optical and X-ray analysis reveal that this low-mass cluster exhibits disturbed dynamics. In fact, it is composed of several galaxy groups in the peripheral regions and, in particular, in the core, where we find evidence of substructures oriented in the NE–SW direction, with hints of a merger nearly along the line of sight. The analysis of the X-ray emission adds further evidence that the cluster is in an unrelaxed dynamical state. At radio wavelengths, the LOFAR data show that the diffuse emission is ~510 kpc in size. Moreover, there are hints of low-surface-brightness emission permeating the cluster center. Conclusions. The environment of the diffuse radio emission is not what we would expect for a classical halo. The spectral index map of the radio source is compatible with a relic interpretation, possibly due to a merger in the N–S or NE–SW directions, in agreement with the substructures detected through the optical analysis. The fragmented, diffuse radio emissions at the cluster center could be attributed to the surface brightness peaks of a faint central radio halo.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

A horseshoe-shaped ring of diffuse emission detected at 1.4 GHz

We identify a horseshoe-shaped ring (HSR) of diffuse emission in J1407+0453 from the Faint Images of Radio Sky at Twenty-cm (FIRST) survey using the Very Large Array telescope at 1.4 GHz. An optical galaxy, SDSS J140709.01+045302.1, is present near the limb of the HSR of J1407+0453, with a spectroscopic redshift ofz = 0.13360. The total extent of the source, including the diffuse emission of J1407+0453, is 65 arcsec (with a physical extent of 160 kpc), whereas the diameter of the HSR is approximately 10 arcsec (25 kpc). The flux density of the HSR is ∼47 mJy at 1400 MHz, whereas the flux densities of the whole diffuse emission of J1407+0453 at 1400 MHz and 150 MHz are 172 mJy and 763 mJy, respectively. We measure the radio luminosity of HSR J1407+0453 as 1.94 × 10 24 W Hz −1 , with a spectral index,α 150 1400 = −0.67. The black hole mass of J1407+0453 is 5.8 × 10 8 M ⊙ . We compare the HSR of diffuse emission of J1407+0453 with other discovered diffuse circular sources. The possible formation scenarios for J1407+0453 are discussed, so as to understand the nature of the source. We present a spectral index map of source J147+0453 to study the spectral properties of the source.

Astronomy & Astrophysics↗
DOE OSTI2020

Accelerated kinetic Monte Carlo: A case study; vacancy and dumbbell interstitial diffusion traps in concentrated solid solution alloys

Vacancy and self-interstitial atomic diffusion coefficients in concentrated solid solution alloys can have a non-monotonic concentration dependence. Here, the kinetics of monovacancies and $\langle$100$\rangle$ dumbbell interstitials in Ni–Fe alloys are assessed using lattice kinetic Monte Carlo (kMC). The non-monotonicity is associated with superbasins, which impels using accelerated kMC methods. Detailed implementation prescriptions for first passage time analysis kMC (FPTA-kMC), mean rate method kMC (MRM-kMC), and accelerated superbasin kMC (AS-kMC) are given. The accelerated methods are benchmarked in the context of diffusion coefficient calculations. The benchmarks indicate that MRM-kMC underestimates diffusion coefficients, while AS-kMC overestimates them. In this application, MRM-kMC and AS-kMC are computationally more efficient than the more accurate FPTA-kMC. Our calculations indicate that composition dependence of migration energies is at the origin of the vacancy’s non-monotonic behavior. In contrast, the difference between formation energies of Ni–Ni, Ni–Fe, and Fe–Fe dumbbell interstitials is at the origin of their non-monotonic diffusion behavior. Additionally, the migration barrier crossover composition—based on the situation where Ni or Fe atom jumps have lower energy barrier than the other one—is introduced. KMC simulations indicate that the interplay between composition dependent crossover of migration energy and geometrical site percolation explains the non-monotonic concentration-dependence of atomic diffusion coefficients.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Robust estimation of charge carrier diffusivity using transient photoluminescence microscopy

Transient microscopy has emerged as a powerful tool for imaging the diffusion of excitons and free charge carriers in optoelectronic materials. In many excitonic materials, extraction of diffusion coefficients can be simplified because of the linear relationship between signal intensity and local excited state population. However, in materials where transport is dominated by free charge carriers, extracting diffusivities accurately from multidimensional data is complicated by the nonlinear dependence of the measured signal on the local charge carrier density. To obtain accurate estimates of charge carrier diffusivity from transient microscopy data, statistically robust fitting algorithms coupled to efficient 3D numerical solvers that faithfully relate local carrier dynamics to raw experimental measurables are sometimes needed. Here, we provide a detailed numerical framework for modeling the spatiotemporal dynamics of free charge carriers in bulk semiconductors with significant solving speed reduction and for simulating the corresponding transient photoluminescence microscopy data. To demonstrate the utility of this approach, we apply a fitting algorithm using a Markov chain Monte Carlo sampler to experimental data on bulk CdS and methylammonium lead bromide (MAPbBr 3 ) crystals. Parameter analyses reveal that transient photoluminescence microscopy can be used to obtain robust estimates of charge carrier diffusivities in optoelectronic materials of interest, but that other experimental approaches should be used for obtaining carrier recombination constants. Additionally, simplifications can be made to the fitting model depending on the experimental conditions and material systems studied. Our open-source simulation code and fitting algorithm are made freely available to the scientific community.

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

Diffusion–convection model of runaway electrons due to large magnetohydrodynamic perturbations in post-thermal quench plasmas

Systematic test particle tracing simulations for runaway electrons (REs) are performed for six post-thermal quench equilibria from DIII-D and ITER, where large scale, kink-like n = 1 (n is the toroidal mode number) magnetohydrodynamic (MHD) instabilities are found. The modeled particle guiding center orbits allow extraction of the effective diffusion–convection coefficients of REs in the presence of large three-dimensional (3D) perturbations up to 10% of the equilibrium toroidal field. With a fixed spatial distribution of the field perturbation, the RE transport coefficients along the plasma radial coordinate track reasonably well with the surface-averaged perturbation level. A substantial variation in the value of the transport coefficients—by three orders of magnitude in most cases, however, occurs with varying launching location of REs along the plasma radius. Large 3D perturbations almost always lead to comparable diffusion and convection processes, meaning that diffusion alone is insufficient to describe the particle motion. At lower (but still high) level of perturbation, the RE convection is found to be dominant over diffusion. A similar observation is made when the perturbation is too strong. In the presence of large perturbation, the dependence of the RE transport on the particle energy is sensitive to the spatial distribution of the perturbation. Based on numerically obtained RE transport coefficients, an analytic fitting model is proposed to quantify the particle diffusion and convection processes due to large MHD events in post-thermal quench plasmas. The model is shown to reasonably well reproduce the direct test particle tracing results for the RE loss fraction and can, thus, be useful for incorporating into other kinetic RE codes in order to simulate the RE beam evolution in the presence of large 3D perturbations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Onset of kinetic effects on Rayleigh–Taylor instability: Advective–diffusive asymmetry

In nature and engineering applications, the Rayleigh–Taylor instability (RTI) occurs over a wide range of Atwood, Reynolds, Mach, and Knudsen numbers. At low Atwood, Mach, and Knudsen numbers, the classic advective instability causes quasi-symmetric bubble and spike growth on the two sides of the interface. However, recent findings suggest that at high degrees of rarefaction, advective effects are suppressed and molecular diffusion leads to planar growth of the density fronts on either side of the interface. This study aims to investigate the flow physics of the transition from advective to diffusive behavior, focusing on the onset of the kinetic effects. Using the gas kinetic methodology, RTI is simulated over a range of Knudsen and Mach numbers in the transition regime. The simulation results reveal the various stages of transformation from advective instability to diffusive transport. For the first time, the study demonstrates the existence of a Knudsen–Mach parameter regime where the bubble side exhibits advective instability, while the other side shows a planar density front due to molecular diffusion, rather than the canonical advective spike shape. The dominance of different mechanisms on the two sides of the interface leads to the advective–diffusive asymmetry. In conclusion, the findings of this study can lead to a more comprehensive understanding of RTI over a wide range of Mach and Knudsen numbers.

42 ENGINEERING↗
DOE OSTI2022

Evidence for superionic H 2 O and diffusive He–H 2 O at high temperature and high pressure

In this work, we present the evidence of superionic phase formed in H 2 O and, for the first time, diffusive H 2 O–He phase, based on time-resolved x-ray diffraction experiments performed on ramp-laser-heated samples in diamond anvil cells. The diffraction results signify a similar bcc-like structure of superionic H 2 O and diffusive He–H 2 O, while following different transition dynamics. Based on time and temperature evolution of the lattice parameter, the superionic H 2 O phase forms gradually in pure H 2 O over the temperature range of 1350–1400 K at 23 GPa, but the diffusive He–H 2 O phase forms abruptly at 1300 K at 26 GPa. We suggest that the faster dynamics and lower transition temperature in He–H 2 O are due to a larger diffusion coefficient of interstitial-filled He than that of more strongly bound H atoms. This conjecture is then consistent with He disordered diffusive phase predicted at lower temperatures, rather than H-disordered superionic phase in He–H 2 O.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Investigations on a fuzzy process: effect of diffusion on calibration and particle identification in Liquid Argon Time Projection Chambers

Ionization electron diffusion in Liquid Argon Time Projection Chambers (LArTPCs) has typically been considered at the detector design stage, but little attention has been given to its effects on calibration and particle identification. We use a GEANT4-based simulation to study how diffusion impacts these techniques, and give consideration to how this effect is simulated. We find that diffusion can cause a drift-dependent bias to both the median and Most Probable Value (MPV) of dQ/dx distributions. The bias is estimated to be ~2.5% (median) and ~5.0% (MPV) for typical maximum drift times in currently running LArTPCs before adding detector specific considerations such as electric field non-uniformities. This information indicates that these metrics should not be used for calibration without care, contrary to the conventional wisdom. The impact of diffusion on the ability of LArTPCs to separate muons and protons is small, and not expected to pose any problems in future detectors. Diffusion may however be a significant source of systematic uncertainty when separating particles of more similar masses (muons and pions, kaons and protons). Separation of such populations may be improved by implementation of a drift-time dependent particle identification.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
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 OSTI2025

Techniques for improved statistical convergence in quantification of eddy diffusivity moments

While recent approaches, such as the macroscopic forcing method (MFM) or Green's function-based approaches, can be used to compute Reynolds-averaged Navier-Stokes closure operators using forced direct numerical simulations, MFM can also be used to directly compute moments of the effective nonlocal and anisotropic eddy diffusivities. The low-order spatial and temporal moments contain limited information about the eddy diffusivity but are often sufficient for quantification and modeling of nonlocal and anisotropic effects. However, when using MFM to compute eddy diffusivity moments, the statistical convergence can be slow for higher-order moments. In this work, we demonstrate that using the same direct numerical simulation (DNS) for all forced MFM simulations improves statistical convergence of the eddy diffusivity moments. We present its implementation in conjunction with a decomposition method that handles the MFM forcing semianalytically and allows for consistent boundary condition treatment, which we develop for both scalar and momentum transport. We demonstrate that for a two-dimensional Rayleigh-Taylor instability case study, using the same DNS for all forced MFM simulations results in convergence with 𝒪⁡(100) simulations rather than 𝒪⁡(1000) simulations. In conclusion, we then demonstrate the impacts of improved convergence on the quantification of the eddy diffusivity.

general physics↗
DOE OSTI2025

Random Walks With Tweedie: A Unified View of Score-Based Diffusion Models [In the Spotlight]

We present a concise derivation for several influential score-based diffusion models that relies on only a few textbook results. Diffusion models have recently emerged as powerful tools for generating realistic, synthetic signals—particularly natural images—and often play a role in state-of-the-art algorithms for inverse problems in image processing. While these algorithms are often surprisingly simple, the theory behind them is not, and multiple complex theoretical justifications exist in the literature. Here, in this study, we provide a simple and largely self-contained theoretical justification for score-based diffusion models that is targeted towards the signal processing community. This approach leads to generic algorithmic templates for training and generating samples with diffusion models. We show that several influential diffusion models correspond to particular choices within these templates and demonstrate that alternative, more straightforward algorithmic choices can provide comparable results. This approach has the added benefit of enabling conditional sampling without any likelihood approximation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Matrix Diffusion Controls Mountain Hillslope Groundwater Ages and Inferred Storage Dynamics

Groundwater age distributions provide fundamental insights on coupled water and biogeochemical processes in mountain watersheds. Field-based studies have found mixtures of young and old-aged groundwater in mountain catchments underlain by bedrock; yet, the processes that dictate these groundwater age distributions are poorly understood. In this work, we use the coupled ParFlow-CLM integrated hydrologic and EcoSLIM particle tracking models to simulate groundwater age distributions on a lower montane hillslope in the East River Watershed, Colorado (USA). We develop a convolution-based approach to propagate fracture-matrix diffusion processes to the EcoSLIM advection-dominated age distributions. We compare observed 3 H and 4 He concentrations from two groundwater wells against model predictions that have varying advective transport times and matrix diffusion magnitudes. Based on a Monte Carlo analysis that considers uncertain matrix and fracture parameters, we find that matrix diffusion is needed to jointly predict 3 H and 4 He observations at both wells. The advection-dominated age distributions lack adequate mixing of young and old-aged water to capture the observed co-occurrence of 3 H and 4 He. The model scenario that best matches the 3 H, 4 He, and water level observations when considering both advective flowpath and matrix diffusion mixing processes has a dynamic bedrock groundwater reservoir that is susceptible to considerable storage losses during low-snow periods. This dynamic groundwater system amplifies the need to assimilate deeper bedrock groundwater into watershed hydro-biogeochemical predictions. This work further highlights the importance of considering matrix diffusion when interpreting environmental tracers in bedrock groundwater systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

A Study of Ferrocene Diffusion in Toluene/Tween 20/1-Butanol/Water Microemulsions for Redox Flow Battery Applications

Redox flow batteries (RFBs) possess multiple advantages as a flexible energy storage solution. However, RFB researchers are still facing many challenges in finding an appropriate electrolyte. Microemulsions have recently been proposed as a promising alternative RFB electrolyte because of their ability to accommodate organic redox species with fast electron transfer rates that are not soluble in aqueous phase, while still offering the high conductivity of an aqueous salt electrolyte. In this work, we focused on understanding the transport of ferrocene (Fc) in a toluene/Tween 20/1-butanol/water model microemulsion and studied the compositional influence on Fc diffusion. The results show that Fc redistributes among the oil, surfactant, and water microenvironments, and the corresponding diffusion and partition coefficients are quantified. Thus, a tortuous path diffusion model is proposed to describe the mass transport of Fc to an electrode surface. Diffusion coefficients are also obtained by pulsed-field gradient nuclear magnetic resonance (PFG NMR), while the values for Fc diffusion are substantially higher than those from electrochemical measurements, suggesting that they measure samples in different ways. The current contributions from each microenvironment indicate that the Fc permeability is much higher in the oil, even though the electron transfer reaction is likely occurring in the surfactant.

25 ENERGY STORAGE↗
DOE OSTI2025

Ab Initio Insights to Metal Passivation: Diffusion and Defect Formation in Amorphous Zirconia and Alumina

Metal passivation refers to the formation of protective oxide films on metals, which shield them from further corrosion and oxidation, playing a crucial role in maintaining their stability. The mesoscopic Point Defect Model has successfully predicted passivity as a steady state process where oxide growth from oxygen vacancies at the metal/film interface competes with oxide dissolution at the film/environment interface. In this work, informed by the Point Defect Model parameters, we use first-principles calculations to calculate defect formation and atomic diffusion in amorphous materials and correlate these descriptors with the behavior and growth of the oxide film. Focusing on amorphous zirconia and alumina, we demonstrate that defect formation energies exhibit significant variability in amorphous systems. In alumina, vacancies dominate, with cation and anion vacancies occurring at comparable concentrations. Diffusion calculations for stoichiometric amorphous alumina and zirconia, as well as oxygen-deficient zirconia, reveal faster diffusion in the oxygen-deficient case, highlighting the impact of defects on transport. Comparison of calculated self-diffusion coefficients for the dominant defect species with experimentally measured oxide thicknesses shows a clear correlation, suggesting that first-principles-derived diffusivity information can serve as a key descriptor for surface passivation film growth.

Farnell, Mackinzie S. [University of California, B↗