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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Bulk Diffusion of Cl through O Vacancies in α-Cr 2 O 3 : A Density Functional Theory Study

Many metals form passive oxides which are broken or weakened by anionic-induced degradation leading to material failure and anionic diffusion is an important step of this degradation process. The diffusion of anionic species through oxides involves a combination of diffusion along grain boundaries, cracks, and channels, as well as bulk diffusion via point defects which is the focus of this study. Using density functional theory, we study bulk diffusion of Cl through O vacancies in α -Cr 2 O 3 as a model system for passive metal oxides. Little is known about Cl diffusion and bonding characteristics, so we benchmark our work through comparison to numerous studies on O diffusion in α -Cr 2 O 3 to analyze similarities and differences between the O and Cl diffusion in the passivation and degradation process of α -Cr 2 O 3 respectively. Unlike O diffusion, the lowest diffusion barrier for Cl is cross-plane diffusion between two (0001) planes through a vacant cation site but the much shorter in-plane diffusion path within the same coordination polyhedron has 35% higher barrier. This work provides the basis for considering the contributions of Cl bulk-diffusion in the overall diffusion kinetics of Cl through α -Cr 2 O 3 .

Oware Sarfo, Kofi↗
DOE OSTI2024

Vacancy diffusion barrier spectrum and diffusion correlation in multicomponent alloys

Vacancy diffusion serves a crucial role in many important kinetic behaviors and properties of multicomponent alloys. Essential questions, however, persist regarding how chemical complexity affects diffusion and what unique characteristics, if any, set these alloys apart from traditional metals. Using neural network kinetics model, we study vacancy diffusion in NbMoTa alloy across a broad temperature range (2600 to 800 K). Unlike pure metals, the two key diffusion parameters—diffusion correlation factor f and activation energy ΔG m —are not constant in alloys, but instead substantially decrease with decreasing temperature. This temperature dependence arises from a reduced number of active vacancy jump pathways at lower temperatures, leading to more correlated diffusion. Upon examining vacancy diffusion throughout the entire compositional space of the Nb-Mo-Ta system, we discover that the slowest vacancy diffusion surprisingly occurs in the non-equimolar region, rather than the equimolar concentration where the configurational entropy is highest. The diffusion barrier spectrum, characterizing the diffusion energy landscape, is an intrinsic material characteristic, which controls both f and ΔG m and, thereby, the diffusivity. Lastly, we find that the vacancy diffusion rate drops noticeably in the presence of local chemical order in the NbMoTa system, particularly for MoTa alloys with long-range B2 order.

Defects↗
DOE OSTI2023

Atomic Diffusion, Segregation, and Grain Boundary Migration in Nickel-Based Alloys from Molecular Dynamics Simulations

Grain boundary diffusion and metal mobility in alloys control material performance in many applications and yet remain poorly understood at a mechanistic level. With advances in accessible time and length scales for computational molecular simulations, and recent force field developments, we now possess tools to help unravel those mechanisms. Using large-scale molecular dynamics simulations, here we examined vacancy-mediated diffusion processes in Ni-5Cr alloy with low and high-energy grain boundaries. We show that atomic diffusion inside the grain boundary plane is about four times higher than bulk diffusion, at any temperature, and exhibits a typical Arrhenius behavior with a very small energy barrier (0~.8 eV for Cr and 0.7 eV for Ni within 1300-1600 K). Additionally, the fastest diffusing species inverts; Cr diffusion was faster than Ni in the bulk but slower in the grain boundaries. This is attributed to the creation of high cohesive energy clusters of Cr at the grain boundary. Grain boundary migration was also observed to be temperature dependent and appears to be two times higher in the 5% Cr alloy than in pure Ni, highlighting the important role of the alloying element on grain boundary motion.

Simonnin, Pauline GN↗
DOE OSTI2021

Mutual Diffusion Coefficients of Na 2 SO 4 at C = 0.02989 mol·dm -3 and 298.15 K by Using Rayleigh Interferometry with Free Diffusion Boundary Conditions: Experimental Test of the Effect of the C 3/2 Concentration Dependence of Refractive Index at Low Concentration

Rayleigh interferometry is one of the most precise and accurate methods for the experimental determination of the mutual diffusion coefficients of both electrolytes and non-electrolytes in liquid solutions. For binary solutions at moderate and high concentrations where the concentration dependences of both diffusion coefficient and refractive index can be assumed to be linear or almost linear and the concentration difference ΔC between a pair of solutions undergoing free diffusion is not large compared to their average concentration, and by using proper combinations of the interference fringe positions as they vary with time, under these conditions diffusion coefficients are obtained that are independent of the size of ΔC. However, this situation becomes more complicated at low concentrations for electrolyte solutions where the concentration dependences of both the diffusion coefficient and refractive index are non-linear, with C 1/2 dependence for the diffusion coefficient and C 3/2 dependence for the refractive index. As a test of these effects on calculated mutual diffusion of electrolyte solutions at low concentrations as measured by Rayleigh interferometry, apparent mutual diffusion coefficients D v,a (volume-fixed reference frame) were measured for fixed average concentrations C=(0.02989 ± 0.00001) mol·dm –3 of Na 2 SO 4 (aq) at 298.15 K while varying ΔC/C from 0.66673 to 2, a threefold variation. The dependence of the diffusion coefficient on ΔC/C was found to be linear, yielding D v = 1.041 6 × 10 –9 m 2 ·s –1 for an infinitely small concentration difference between the diffusing solutions.

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

Effect of Confinement of SIA Cluster Diffusion by Impurities on Radiation Defect Accumulation due to 14MeV neutrons in Tungsten

Impurities acting as traps, reduce the self-interstitial atom (SIA) cluster mobility, thereby increasing recombination and in turn enhancing radiation-resistance of polycrystalline materials. However, at elevated temperatures, depending on SIA-trap (impurity) binding energy, SIA clusters may trap and detrap (detach) numerous times. In tungsten, SIA clusters of all sizes glide one-dimensionally (1D) along their Burgers vector direction. Hence, detrapped SIAs would retrace their path, confining their 1D-glide between traps. However, small SIA clusters can change the direction of their 1D-glide by overcoming a rotation barrier. In addition, impurities or solutes are known to reduce this barrier. A lower rotation barrier can result in detrapped SIA clusters that diffuse in a random 1D direction, effectively leading to 3D (net-3D) diffusion. Here we investigate the effect of two cases of SIA diffusion, namely confined-1D diffusion and net-3D diffusion, on the damage accumulation in polycrystalline tungsten under 14 MeV neutron irradiation using the object kinetic Monte Carlo method. Simulations were performed using cascade debris obtained from molecular dynamics simulations with primary knock-on atom (PKA) energies corresponding to 14 MeV neutrons. In the simulations with net-3D diffusion, SIA clusters of all sizes diffuse in a random 1D direction after the clusters detrap. While for the confined 1D diffusion case, SIA clusters larger than size 5 retain their original direction of 1D-glide. In both types of simulations, trapped vacancies are assumed to be permanently immobilized. Unexpectedly, the damage accumulation was lower with confined 1D diffusion than that with net-3D. We present a systematic comparison of the influence of the two cases of SIA diffusion on the radiation damage accumulation as a function of dose rate, detrapping barrier, and impurity concentration at 1025 K.

Accumulation and recovery of lattice damage, tungs↗
DOE OSTI2023

Code for the manuscript titled "Blackout Diffusion: Generative Diffusion Models in Discrete-State Spaces"

We would like to disclose two scripts, written in Jupyter notebook, in which we implement the "Blackout Diffusion Process" described in the manuscript "Blackout Diffusion: Generative Diffusion Models in Discrete-State Spaces" (LA-UR-23-20509), to be submitted to the International Conference of Machine Learning (ICML). The abstract of the manuscript is append below. == Typical generative diffusion models rely on a Gaussian diffusion process for training the backward transformations, which can then be used generate samples from Gaussian noise. However, real world data often takes place in discrete-state spaces, which occur in many scientific applications. Here, we develop a theoretical formulation for arbitrary discrete-state Markov processes in the forward diffusion process. We relate the theory to the existing continuous-state Gaussian diffusion and identify the corresponding reverse-time stochastic process and score function in the continuous-time setting, and the reverse-time mapping in the discrete-time setting. As an example of this framework, we introduce "Blackout Diffusion", which learns to produce samples from an empty image instead of from noise. Numerical experiments on the CIFAR-10 dataset confirm the feasibility of generative diffusion modeling in a discrete space. Generalizing from specific (Gaussian) forward processes to a more general framework also sheds light on how to interpret generative diffusion models and their mathematical structure, which we comment on.

Lin, Yen Ting↗
DOE OSTI2021

Modifying Li + and Anion Diffusivities in Polyacetal Electrolytes: A Pulsed-Field-Gradient NMR Study of Ion Self-Diffusion

Polyacetal electrolytes have been demonstrated as promising alternatives to liquid electrolytes and poly(ethylene oxide) (PEO) for rechargeable lithium-ion batteries; however, the relationship between polymer structure and ion motion is difficult to characterize. Here, we study structure-property trends in ion diffusion with respect to polymer composition for a systematic series of five polyacetals with varying ratios of ethylene oxide (EO) to methylene oxide (MO) units, denoted as P(xEO-yMO), and PEO. We first use 7 Li and 19 F pulsed-field-gradient NMR spectroscopy to measure cation and anion self-diffusion, respectively, in polymer/lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) salt mixtures. At 90 °C, we observe modest changes in Li + diffusivity across all polymer compositions, while anion (TFSI - ) self-diffusion coefficients decrease significantly with increasing MO content. At a given reduced temperature (T - T g ), all polyacetal electrolytes exhibit faster Li + self-diffusion than PEO. Intriguingly, P(EO-MO) and P(EO-2MO) also show slower TFSI - anion self-diffusion than PEO at a given reduced temperature. Molecular dynamics simulations reveal that shorter distances between acetal oxygen atoms (O-CH 2 -O) compared to ether oxygens (O-CH 2 -CH 2 -O) promote more diverse, often asymmetric, Li + coordination environments. Finally, Raman spectra reveal that anion-rich ion clusters in P(EO-MO) and P(EO-2MO) lead to decreased anion diffusivity, which along with increased cation diffusivity, support the viability of polyacetals as high-performance polymer electrolytes.

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

On the Importance of Using Event-Specific Wave Diffusion Rates in Modeling Diffuse Electron Precipitation

A few to tens of keV electron precipitation that carries substantial energy source down to the upper atmosphere to create aurora is manifested as an important magnetosphere-ionosphere coupling process. The precipitation is usually caused by scattering processes associated with plasma waves in the magnetosphere. The scattering process is often quantified by wave diffusion rates that indicate how fast an electron is scattered. Global models commonly use diffusion coefficients that are derived from statistical wave models. However, due to the statistical nature, many localized, transient features could be smeared out. In this study, we investigate electron precipitation using event-specific diffusion coefficients that are obtained based on simultaneous in-situ measured/inferred, rather than statistical, chorus wave dynamics. We find that the application of the event-specific diffusion coefficients associated with a more dynamic and intense chorus wave model leads more electrons, particularly at several to tens of keV in the dawn-to-noon sector at L > 3, to precipitate than using statistical coefficients. Here, the new simulation roughly captures both the intensity and variability of the precipitating flux as detected by the NOAA/POES satellites. Ionospheric electron density in the lower E region (100–120 km) observed by the mid-latitude Millstone Hill radar is also much better reproduced, while the case using statistical diffusion coefficients underestimates the ionization rate. This study implies the importance of using event-specific diffusion rates in simulating the diffuse electron precipitation and understanding the magnetosphere-ionosphere coupling.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Blackout Diffusion: Generative Diffusion Models in Discrete-State Spaces

Typical generative diffusion models rely on a Gaussian diffusion process for training the backward transformations, which can then be used to generate samples from Gaussian noise. However, real world data often takes place in discrete-state spaces, including many scientific applications. Here, we develop a theoretical formulation for arbitrary discrete-state Markov processes in the forward diffusion process using exact (as opposed to variational) analysis. We relate the theory to the existing continuous-state Gaussian diffusion as well as other approaches to discrete diffusion, and identify the corresponding reverse-time stochastic process and score function in the continuous-time setting, and the reverse-time mapping in the discrete-time setting. As an example of this framework, we introduce “Blackout Diffusion”, which learns to produce samples from an empty image instead of from noise. Numerical experiments on the CIFAR-10, Binarized MNIST, and CelebA datasets confirm the feasibility of our approach. Generalizing from specific (Gaussian) forward processes to discrete-state processes without a variational approximation sheds light on how to interpret diffusion models, which we discuss.

Santos, Javier E.↗
DOE OSTI2021

Role of a Multivalent Ion–Solvent Interaction on Restricted Mg 2+ Diffusion in Dimethoxyethane Electrolytes

The diffusion behavior of Mg 2+ in electrolytes is not as readily accessible as that from Li + or Na + utilizing PFG NMR, due to the low sensitivity, poor resolution, and rapid relaxation encountered when attempting 25 Mg NMR. In MgTFSI 2 /DME solutions, “bound” DME (coordinating to Mg 2+ ) and “free” DME (bulk) are distinguishable from 1 H NMR. With the exchange rates between them obtained from 2D 1 H EXSY NMR, we can extract the self-diffusivities of free DME and bound DME (which are equal to that of Mg 2+ ) before the exchange occurs using PFG diffusion NMR measurements coupled with analytical formulas describing diffusion under two-site exchange. Further, the high activation enthalpy for exhange (65–70 kJ/mol) can be explained by the structural change of bound DME as evidenced by its reduced C–H bond length. Comparison of the diffusion behaviors of Mg 2+ , TFSI – , DME, and Li + reveals a relative restriction to Mg 2+ diffusion that is caused by the long-range interaction between Mg 2+ and solvent molecules, especially those with suppressed motions at high concentrations and low temperatures.

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

Halide Substitution Effects on Lithium-Ion Diffusion in Protonated Antiperovskites

Solid-state electrolytes (SSEs) for all-solid-state lithium-ion batteries are generating intense interest because these batteries can improve the safety and performance compared with devices fabricated with conventional, flammable liquid electrolytes. In these SSEs, it has been suggested that Li + ion diffusion in the grain boundaries is hindered and is a critical determinant of the overall ionic conductivity (σ). However, Li + ion diffusivities in the grain (D G ) and the grain boundary (D GB ) are difficult to determine experimentally, with few techniques capable of distinguishing the individual contributions. Here, we distinguished the D G and D GB for the protonated lithium antiperovskites (pLiAPs) SSEs: Li 2 OHCl, Li 2 OHBr, Li 2 OHF 0.1 Cl 0.9 , Li 2 OHF 0.1 Br 0.9 , and Li 2 OHCl 0.3 7Br 0.63 . The measurements were obtained directly from 7 Li pulsed-field gradient nuclear magnetic resonance (PFG-NMR) at 353 K. The 7 Li PFG-NMR echo profiles were composed of two primary components with additional secondary oscillatory components – the so-called NMR diffraction phenomenon. The length scale separating the two main components corresponds to a diffusion length of ~1.7 µm, which is thought to be the average grain size (by diameter). The short-range (≤ 1.7 µm) diffusion component associated with D G (≈10 -11 m 2 /s) varied minimally with halide substitution, while the long-range (≥ 1.7 µm) component D GB (≈10 -12 to 10 -15 m 2 /s) was highly sensitive to the substitution of halides and closely correlated with s. In addition, from the comparison of the ratio D GB /D G to D t (the Li + ion diffusion coefficient estimated from the rotational correlation time, t c ), it was determined that the contribution of D G to σ is negligible; 0.01 ~ 0.04 in the pLiAPs studied here. Finally, these insights provide fundamental understanding of the halide substitution effects on Li + ion grain versus grain boundary diffusion, and suggest that careful engineering of the grain boundaries at the microscopic level is necessary to achieve high-performance pLiAP SSEs.

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

Graph-theoretical KMC Framework for Calculating Effective Diffusivity in TPBAR Components: Effective Diffusivity of Tritium in α-Zr Grain Boundaries

We report the development of a computer simulation tool based on the Graph-theoretical kinetic Monte Carlo (GT-KMC) approach that can be used to simulate diffusion and derive effective diffusivity of species in complex structures of TPBAR components. The GT-KMC framework was implemented in AKSOME, an on-lattice self-learning kinetic Monte Carlo tool developed at PNNL and was previously used to study solute diffusion in metal alloys. This newly developed tool will be helpful in simulating atom diffusion and extracting their effective diffusivities in grain boundaries and interfaces in various TPBAR components. Once benchmarked, the tool will be used to simulate the diffusion of tritium (hydrogen) along grain boundaries in α-Zr using the activation energy barrier data previously obtained from Density Functional Theory calculations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Production and diffusion of H 2 O 2 during the interaction of a direct current pulsed atmospheric pressure plasma jet on a hydrogel

The interaction of cold atmospheric pressure plasma jets with hydrogels has been used as a model system to study the interaction of plasmas with tissues. In this study, we analyze the diffusion of reactive oxygen species (in particular H 2 O 2 ) and quantify the amount of plasma-produced H 2 O 2 species that penetrates into a gelatin hydrogel. We show that the diffusion constant of H 2 O 2 in 10% gelatin hydrogel is similar to its diffusion constant in water and that the production of H 2 O 2 in the hydrogel is significantly less than the production of H 2 O 2 in distilled water for the same plasma operation conditions suggesting that the scavenging of OH radicals at the plasma-gel interface significantly reduces the H 2 O 2 production.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Cation–Ligand Interactions Dictate Salt Partitioning and Diffusivity in Ligand-Functionalized Polymer Membranes

Membranes are an attractive alternative to current thermal separations due to their scalability and energy efficiency in desalinating water. Unfortunately, many of the conventional membrane materials available today are unable to differentiate between ionic solutes, especially alkali cations, compromising their use in ion–ion separations. Inspired by the ion-specific interactions exhibited by biological ion channels, recent research efforts have focused on synthesizing and characterizing new polymeric materials that incorporate ligands into polymer networks to bias solubility and/or diffusivity of one cationic species over another. Despite these efforts, little is known about the influence of incorporating ligands into polymer membranes on solubility and diffusivity of the complexing species. In this study, we first build a qualitative model of salt partitioning, diffusivity, and permeability in generic cation-complexing ligand-functionalized polymer membranes. Next, to validate our model and hypotheses, we perform atomistic molecular dynamics simulations of a 12-crown-4-functionalized membrane in the presence of alkali halide salts at low concentration. Generally, cation complexation enhances cation solubility but decreases diffusivity. Interestingly, the reduction in diffusivity is predicted to be larger than the enhancement in solubility for materials which operate by the mechanisms proposed in our physical picture, ultimately resulting in a reduction in the permeability of the selectively complexing ion.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Role of surface diffusion in formation of unique reactivity for graphite oxidation: Time-resolved measurements in a pulsed diffusion reactor

Quantification of oxidation kinetics is essential to develop graphitic materials for diverse applications: from refractories found in gas-cooled nuclear reactors to catalysts needed for chemical manufacturing. In this work, using well-defined highly oriented pyrolytic graphite, low-pressure isotopic transient experiments combined with controlled annealing periods, we resolve the role of surface diffusion and quantify oxidation kinetics with nanomole-precision. We observe an unexpected increase in reactivity following annealing which is explained by the role of surface diffusion increasing the probability for trapping mobile oxygen at more reactive edge sites. Here, the locus of adsorption and spillover to the basal plane is distinct from the trapping location creating a more active oxygen species. Isotopic products reflect the population dynamics of oxygen added at the edge and surface diffusion that relocates basal plane oxygen to more reactive edge sites. Since this process proceeds in parallel with direct oxidation reactions, it is not likely to be observed using steady-state or conventional ‘bulk’ characterization techniques. Our unique time-resolved non-equilibrium measurement in a well-defined transport regime, enables observation of three distinct behaviors: short-term deactivation due to the balance of rates in oxygen supply/product formation, reactivity increases due to surface diffusion and longer-term reactivity increase with oxygen accumulation.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Influence of Hydrogen Isotopes on Vacancy Formation and Antisite Defect Diffusion in Palladium and Vanadium Metals

Density functional theory calculations are performed to study the interaction between hydrogen isotopes, vacancy, and antisite defects in Pd and V. Various defect configurations and defect-defect distances are systematically explored. Binding energies and migration energy barriers are calculated and tabulated. The study provides atomistic data for subsequent mesoscale simulations of tritium, vacancy, and antisite diffusion. In Pd, a strong binding is found between a vacancy and tritium (0.16 eV). This tritium binding increases the thermal concentration of vacancies by a factor of ~10 at 500 ?C. The increase in vacancy concentration enhances V antisite diffusion in Pd by a factor of ~6 at 500 ?C. The influence of tritium is even stronger in V, with a tritium and vacancy binding energy of 0.38 eV. Such a strong binding increases the thermal concentration of vacancies by a factor of ~300 at 500 ?C. The increase in vacancy concentration enhances Pd antisite diffusion in V by a factor of ~640 at 500 ?C. Vanadium and Pd exhibit a strong driving force to intermix with a formation energy of -1.57 eV for V antisite in Pd and -1.05 eV for Pd antisite in V. The results suggest vanadium diffusion into Pd is energetically stronger than the reverse. Zero-point-energy corrections are taken into account and calculations for hydrogen and tritium are presented.

Setyawan, Wahyu↗
DOE OSTI2021

Water dynamics in C–S–H and M-S-H cement pastes: A revised jump-diffusion and rotation-diffusion model

The Quasi-Elastic Neutron Scattering (QENS) spectra from four cement pastes are re-analyzed by a new revised jump-diffusion and rotation-diffusion model (rJRM). From the QENS fit, it can be seen that the rJRM is an improved model to fit QENS spectra within the whole detected neutron energy transfer and scattering vector. By the rJRM fitting, the structure parameters extracted from QENS spectra show that both the additives aluminum-silicate nanotubes (ASN) and carboxyl group functionalized ASN (ASN-COOH) can improve magnesium-silicate-hydrate (M-S-H) toward calcium-silicate-hydrate (C–S–H) direction in mechanical properties, but the improvement is weakened with decreasing temperatures. The extracted dynamical parameters show that there is a dynamic anomaly near 230 K in not only translational but also rotational diffusion of water confined in all the investigated samples. In conclusion, the anomaly in rotational diffusion is new compared with those results obtained by the QENS fit using other models.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗