Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Reveal the fast and charge-insensitive lattice diffusion of silver in cubic silicon carbide via first-principles calculations

The lattice diffusion of silver with various charge states in cubic silicon carbide has been investigated by means of high-fidelity spin-polarized density functional theory calculations. The migration energy barrier of the Ag interstitial diffusion is 1.09 eV and 1.11 eV for neutral and q = +1e charge state, respectively, close to the activation energy of Ag diffusion measured in the German HTR fuel program. A general trend is that the migration energy barrier reduces with respect to the increase of charge on Ag, which is much less than the increment in the oxidation energy, suggesting that the lattice diffusion of silver prefers constant neutral state without redox in transition state. Furthermore, our results indicate a scenario that once Ag is deposited to interstitials via the kickout mechanism, it will highly likely perform lattice diffusion across the SiC layer leading to fast release of Ag in Tristructural-Isotropic fuel induced by irradiation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

A finite difference informed random walker (FDiRW) solver for strongly inhomogeneous diffusion problems

In nature, many complex multi-physics coupling problems exhibit strong diffusivity inhomogeneity. For instance, in the context of radionuclide absorption by porous wasteform materials within a flowing waste stream, the difference of species’ diffusivity in solid and liquid phases spans by 3~8 orders of magnitude. To solve the diffusion equations with strongly inhomogeneous diffusivity, traditional discretization-based methods, such as the Finite Difference Method (FDM), require infinitesimally small time steps (<10 -10 ) as high spatial resolutions are employed in most microstructure evolution processes, leading to prohibitively high computational costs. Here, this work developed an integrated numerical approach (FDiRW: Finite Difference informed Random Walk) to tackle this challenge. The idea is that utilizing the Random Walk concept, the fast diffusion is modeled as a superposition of point source’s solution for a concentration distribution while FDM is used to obtain the point source’s solution at each node. A mesh-coarsening algorithm is developed to generate an exclusive coarse mesh for FDiRW approach to maximize its efficiency. The effectiveness of the coarse mesh-based FDiRW approach is validated by benchmarking Finite Difference solutions. Numerical results demonstrated that FDiRW achieves a remarkable 1000x computational efficiency improvement over FDM while preserving desired accuracy for a medium-sized model of 192 × 192 × 192 grids. Finally, as models scale up, a floating-point operations (PLOPs) analysis of the FDiRW algorithm reveals that its computational complexity grows quadratically in terms of the number of nodes employed in computation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A provably stable numerical method for the anisotropic diffusion equation in confined magnetic fields

We present a novel numerical method for solving the anisotropic diffusion equation in magnetic fields confined to a periodic box which is accurate and provably stable. We derive energy estimates of the solution of the continuous initial boundary value problem. A discrete formulation is presented using operator splitting in time with the summation by parts finite difference approximation of spatial derivatives for the perpendicular diffusion operator. Weak penalty procedures are derived for implementing both boundary conditions and parallel diffusion operator obtained by field line tracing. We prove that the fully-discrete approximation is unconditionally stable. Discrete energy estimates are shown to match the continuous energy estimate given the correct choice of penalty parameters. A nonlinear penalty parameter is shown to provide an effective method for tuning the parallel diffusion penalty and significantly minimises rounding errors. Several numerical experiments, using manufactured solutions, the “NIMROD benchmark” problem and a single island problem, are presented to verify numerical accuracy, convergence, and asymptotic preserving properties of the method. Finally, we present a magnetic field with chaotic regions and islands and show the contours of the anisotropic diffusion equation reproduce key features in the field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Data on the comprehensive first-principles diffusion study of the aluminum-magnesium system

First-principles calculation of diffusion coefficients between Mg and Al is investigated comprehensively using density functional theory (DFT). The effect of different uncertainty sources arising from first principles calculations has been investigated systematically. These sources include the diffusion model, energetic, entropic and attempt frequency calculations. Variation in self and impurity diffusion coefficients of Mg and Al in stable phases are quantified using different DFT settings and compared with the experiments. Using the optimal DFT settings, diffusion coefficients in metastable phases of Al and Mg are predicted. The dataset refers to “An integrated experimental and computational study of diffusion and atomic mobility of the aluminum-magnesium system” [1].

36 MATERIALS SCIENCE↗
DOE OSTI2026

Multi-Objective design of interlocking metasurfaces using conditional diffusion models

Unit cell design remains a major challenge for interlocking metasurfaces, a promising joining technology for dissimilar materials, due to the complex, competing, multivariate design space and the need for rapid adaptation to varying performance requirements. This study explores Conditional Diffusion Models as a design optimization tool for interlocking metasurfaces. Given the complex, competing, multivariate design space for interlocking metasurfaces, unit cell design remains a major challenge for this joining technology. We trained a conditional diffusion model on 25,000 finite element analysis-simulated interlocking metasurface unit cells to generate designs with tailored thermo-mechanical properties (tensile strength, shear strength, and thermal conductivity) based on specified performance criteria. The model demonstrated a success rate of approximately 72 % in producing designs that met specified property bounds. The conditional diffusion model generated both thermally resistive and conductive designs, revealing clear trends in design characteristics: taller, dendritic structures were advantageous for tensile loads, while shorter, robust designs excelled in shear applications. Our findings indicate that the model's performance is more influenced by the breadth of the design space than by the quantity of training data, highlighting the importance of expansive design domains for generating innovative solutions. This work establishes conditional diffusion models as a highly efficient and adaptable tool for rapid interlocking metasurface unit cell design, paving the way for advancements in multi-material joining technologies, as well as highlighting the justification to leverage conditional diffusion models as design tools across complex design domains.

Conditional diffusion models↗
DOE OSTI2019

Concentration dependence of hydrogen diffusion in α-iron from atomistic perspectives

Evaluation of hydrogen diffusion in structural materials is essential to predict the leakage and embrittlement of hydrogen storage applications. Here, we investigate the atomic-scale diffusion of interstitial hydrogen (H) in α-iron (Fe) over a temperature range from 350 to 900 K with different H concentrations (0.01–5%), employing classical molecular dynamics (MD) simulations. The self-diffusivity of H atoms increases with increasing temperature and decreasing concentration. With low concentrations, the calculated diffusion properties agree well with prior experiments. However, with a higher concentration (≥1%), the H diffusivity at low temperatures deviates from a high-temperature Arrhenius behavior. Through the energetic and structural analysis, we suggest that this deviation is attributed to a reduced mobility due to increased energy barrier by other H interstitials. This work contributes to the effective design of H storage applications by identifying temperature and concentration effects on permeability and addressing possible microstructural transformation.

08 HYDROGEN↗
DOE OSTI2022

Assessment of diffuse-interface methods for compressible multiphase fluid flows and elastic-plastic deformation in solids

This work describes three diffuse-interface methods for the simulation of immiscible, compressible multiphase fluid flows and elastic-plastic deformation in solids. The first method is the localized-artificial-diffusivity approach of Cook, Subramaniam et al., and Adler and Lele, in which artificial diffusion terms are added to the individual phase mass fraction transport equations and are coupled with the other conservation equations. The second method is the gradient-form approach that is based on the quasi-conservative method of Shukla et al., in which the diffusion and sharpening terms (together called regularization terms) are added to the individual phase volume fraction transport equations and are coupled with the other conservation equations. The third approach is the divergence-form approach that is based on the fully conservative method of Jain et al., in which the regularization terms are added to the individual phase volume fraction transport equations and are coupled with the other conservation equations. In the present study, all three diffuse-interface methods are used in conjunction with a four-equation, multicomponent mixture model, in which pressure and temperature equilibria are assumed among the various phases. The primary objective of this work is to compare these three methods in terms of their ability to: maintain constant interface thickness throughout the simulation; conserve mass, momentum, and energy; and maintain accurate interface shape for long-time integration. The second objective of this work is to consistently extend these methods to model interfaces between solid materials with strength. To assess and compare the methods, they are used to simulate a wide variety of problems, including (1) advection of an air bubble in water, (2) shock interaction with a helium bubble in air, (3) shock interaction and the collapse of an air bubble in water, and (4) Richtmyer–Meshkov instability of a copper–aluminum interface. The current work focuses on comparing these methods in the limit of relatively coarse grid resolution, which illustrates the true performance of these methods. In conclusion, this is because it is rarely practical to use hundreds of grid points to resolve a single bubble or drop in large-scale simulations of engineering interest.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Enabling probabilistic learning on manifolds through double diffusion maps

Here, we present a generative learning framework for probabilistic sampling that extends Probabilistic Learning on Manifolds (PLoM), which is designed to generate statistically consistent realizations of a random vector in a finite-dimensional Euclidean space, informed by a (representative) set of observations. In its original form, PLoM constructs a reduced-order probabilistic model by combining three main components: (a) kernel density estimation to approximate the underlying probability measure, (b) Diffusion Maps to characterize the manifold of the data, and (c) a reduced-order Itô Stochastic Differential Equation (ISDE) to sample from the learned distribution. However, its sampling dynamics are posed in the ambient space and the retained number of reduced coordinates is chosen by projection-reconstruction error. In practice, this often (i) requires more coordinates than the data’s intrinsic dimension to achieve stable sampling and (ii) lacks a smooth, basis-independent lifting back to the data domain; moreover, standard Diffusion Maps emphasize harmonic eigenfunctions and can miss non-harmonic latent structure. We address these limitations by decoupling geometry learning from sampling: a first Diffusion Maps pass identifies non-harmonic coordinates on which we formulate a full-order ISDE directly in the latent space, while Double Diffusion Maps captures multiscale geometric features and Geometric Harmonics (GH) learns a smooth lifting map to the ambient variables that is independent of the particular diffusion basis. This hybrid design preserves the system’s dynamical richness with a compact geometric representation and enables principled out-of-sample inference. The effectiveness and robustness of the proposed method are illustrated through two numerical studies: one based on data generated from two-dimensional Hermite polynomial functions and another based on high-fidelity simulations of a detonation wave in a reactive flow.

Double diffusion maps↗
DOE OSTI2023

Pore connectivity influences mass transport in natural rocks: Pore structure, gas diffusion and batch sorption studies

For this work, six rocks (one granodiorite, one limestone, two chalks, one mudstone, and one dolostone) with different extents of heterogeneity at six different particle sizes (from 75 to 8000 μm) were studied to describe the effects of pore connectivity on mass transport. The methods applied were (i) porosity measurement of granular rocks, (ii) analyses of gas-phase diffusive transport in a bed of packed particles, along with a solid quartz method at these six particle sizes being developed to identify the contribution of intraparticle diffusion, and (iii) batch sorption tests of multiple ions (anions and cations) with subsequent analyses of inductively coupled plasma-mass spectrometry. Granular porosity measurement results reveal that with decreasing particle sizes, the effective porosities for the “heterogenous” group of rocks (Grimsel granodiorite and Edwards limestone) increase, whereas the porosities of another “homogeneous” group (two Israel chalk samples, Japan mudstone, and Wyoming dolostone) remain constant. Gas diffusion results show that the intraparticle gas diffusion coefficient among these two sample groups, varying in the magnitude of 10 -8 to 10 -6 m 2 /s, are not directly correlated to the porosity differences. Moreover, the batch sorption work displays a different affinity of rocks for various tracers. For Grimsel granodiorite, Japan mudstone, and Wyoming dolostone, the adsorption capacity of Sm 3+ and Eu 3+ increases as the particle size decreases. In general, this integrated research of grain size distribution, granular rock porosity, intraparticle diffusivity, and ionic sorption capacity gives insights into the pore connectivity effect on both physical and chemical transport behaviors for different lithologies and/or different particle sizes.

58 GEOSCIENCES↗
DOE OSTI2021

Molecular dynamics studies of lattice defect effects on tritium diffusion in zirconium

Tritium diffusion in α-Zr containing point defects such as vacancies or self-interstitial atoms (SIAs) is simulated using molecular dynamics. Point defects rapidly aggregate to form extended defects, such as 3D nanoclusters and Frank loops. The geometry of extended defects is affected by the presence of tritium. At low temperature and in the absence of tritium, vacancies aggregate to form stacking fault pyramids. Addition of tritium at these temperatures promotes aggregation of vacancies to form 3D nanoclusters, within which the tritium concentration can be sufficiently high to suggest that these defects may serve as nucleation sites for hydride precipitation. Trapping of tritium in vacancy nanocluster reduces the calculated bulk diffusivity by an amount proportional to the vacancy concentration. At high temperature, vacancy clusters change shape to form planar basal dislocation loops, which bind tritium less strongly, leading to a sharp reduction in the fraction of trapped tritium and a corresponding increase in tritium diffusivity at high temperature. In contrast, SIAs increase tritium diffusion through α-Zr. Analysis of atomic trajectories shows that tritium does not interact directly with SIAs. In conclusion, diffusion enhancement is instead related to expansion of the lattice.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Diffusion and migration in polymer electrolytes

Mixtures of neutral polymers and lithium salts have the potential to serve as electrolytes in next-generation rechargeable Li-ion batteries. The purpose of this review is to expose the delicate interplay between polymer-salt interactions at the segmental level and macroscopic ion transport at the battery level. Since complete characterization of this interplay has only been completed in one system: mixtures of poly(ethylene oxide) and lithium bis(trifluoromethanesulfonyl)imide (PEO/LiTFSI), we focus on data obtained from this system. We begin with a discussion of the activity coefficient, followed by a discussion of six different diffusion coefficients: the Rouse motion of polymer segments is quantified by D seg , the self-diffusion of cations and anions is quantified by D self,+ and D self,- , and the build-up of concentration gradients in electrolytes under an applied potential is quantified by Stefan-Maxwell diffusion coefficients, D 0+ , D 0- , and D +- . The Stefan-Maxwell diffusion coefficients can be used to predict the velocities of the ions at very early times after an electric field is applied across the electrolyte. The surprising result is that D 0- is negative in certain concentration windows. A consequence of this finding is that at these concentrations, both cations and anions are predicted to migrate toward the positive electrode at early times. We describe the controversies that surround this result. Knowledge of the Stefan-Maxwell diffusion coefficients enable prediction of the limiting current. We argue that the limiting current is the most important characteristic of an electrolyte. Excellent agreement between theoretical and experimental limiting current is seen in PEO/LiTFSI mixtures. What sequence of monomers that, when polymerized, will lead to the highest limiting current remains an important unanswered question. It is our hope that the approach presented in this review will guide the development of such polymers.

25 ENERGY STORAGE↗
DOE OSTI2018

TiN diffusion barrier for stable W/SiC(0001) interfaces in inert ambient at high temperature

The effect of high-temperature annealing on tungsten (W) films deposited on silicon carbide (SiC) with and without a titanium nitride (TiN) diffusion barrier was examined as a function of time. Evolutions in phase composition, surface morphology, and roughness from annealing at 1273 K were investigated for up to 24 h. Without a TiN diffusion barrier, solid state reactions between the W film and SiC substrate led to the formation of W 5 Si 3 , W 2 C, and WC species and the rise of an inhomogeneous surface structure that was initially web-like and later discontinuous. Severe roughening on the order of the initial film thickness was observed. Incorporation of a 100 nm TiN diffusion barrier suppressed the formation of W 5 Si 3 and W 2 C and only trace WC could be detected due to species diffusion through TiN grain boundaries. Changes in the surface structure and roughness were minimal. Our results warrant consideration of TiN as an effective diffusion barrier for W on SiC systems where structural stability at high temperatures is highly desired.

36 MATERIALS SCIENCE↗
DOE OSTI2019

A Quasielastic Neutron Scattering Study of Water Diffusion in Model Anion Exchange Membranes over Localized and Extended Volume Increments

Anion exchange membranes provide the basis for potentially more cost-effective fuel cells, electrochemical energy conversion, and energy storage. The membranes consist of a polymer domain with tethered cationic side groups and a pervading aqueous domain that conducts water and anions. Water diffusion is investigated through temperature-dependent quasielastic neutron scattering (QENS) at application relevant water contents, employing the BASIS spectrometer at ORNL. The spectrometer resolves two distinct quasielastic components corresponding to localized water diffusion (on length scales lower than the water domain thickness) and to extended water diffusion (on length scales approaching the water domain thickness). Water diffusion coefficients in membranes are time and length scale dependent and water exchanges between diffusion on a localized and extended volume increment. An exchange model for QENS is developed to describe the exchange of water between the two volume increments. In conclusion, the populations and lifetimes of water in each volume are discussed in terms of their connection to long-range water transport and membrane nanomorphology.

Diffusion↗
DOE OSTI2020

Microscopic Theory of Long-Time Center-of-Mass Self-Diffusion and Anomalous Transport in Ring Polymer Liquids

We construct a microscopic theory at the level of segment-scale correlated space–time intermolecular forces for the long-time center-of-mass (CM) diffusion constant and intermediate-time non-Fickian transport in dense solutions and melts of ring polymers. The approach combines ideas of polymer, colloid, and liquid-state statistical mechanics to quantify how the multifractal intra-ring conformational structure and inter-ring packing correlations determine dynamic caging constraints and time-dependent friction. Breakdown of Rouse theory is predicted to occur due to length scale-dependent temporal correlation of forces exerted on pairs of tagged ring segments from surrounding polymers. At large enough degrees of polymerization (N), a stronger scaling of the CM diffusion constant (D ∝ N –2 ) is predicted, with a crossover N D proportional to the product of the system-specific macromolecular volume fraction and dimensionless compressibility. In analogy with caging effects in glass-forming fluids, the theory does appear to begin to fail at sufficiently high N/N D for the center-of-mass diffusivity, likely due to another crossover to an even slower activated transport regime. However, use of N/N D with the theoretically predicted N D , in conjunction with dynamic blob scaling ideas for local Rouse friction, collapses semidilute and concentrated solution simulation data onto a master curve. Based on the same physical ideas employed to predict the diffusion constant, a generalized Langevin equation description is formulated for intermediate-time CM transport. It predicts two subdiffusive regimes that emerge due to the self-similar nature of internal ring structure. Analytic and numerical predictions for the apparent non-Fickian exponent as a function of time, N/N D and dimensionless compressibility, the evolution of the maximum degree of subdiffusive motion with N/N D , and the time scale for recovering Fickian diffusion are made, all of which are in good accordance with melt simulations. Furthermore, the present work sets the stage to address activated dynamics and glass formation on the macromolecular scale.

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

Precise Linker Length and Dynamic Bond Exchange Control Penetrant Diffusion in Dense Vitrimers

Polymer networks with dynamic covalent bonds have been investigated for their self-healing ability, recyclability, and potential as more sustainable materials. Recent results have indicated that in some cases, bond exchange can enhance the transport of penetrants in dense networks, pointing to their potential for separations of membranes. Here, imine dynamic bonds in ethylene oxide (EO) networks with precise linker lengths were synthesized to investigate the transport of N,N′-bis(2,5-di-tert-butylphenyl)-3,4,9,10-perylenedicarboximide (BTBP), a large, anisotropic dye molecule. Networks with mesh sizes smaller than, comparable to, and greater than the size of the penetrant axes were investigated to probe the effects of bond exchange and network confinement on transport. Mesh sizes, which ranged from 0.5 to 1.62 nm, were determined from shear rheology, glass transitions by calorimetry, and probe diffusion coefficients by fluorescence recovery after photobleaching. Permanent networks with identical EO chain lengths were prepared as control samples, and up to a 3 orders of magnitude increase in diffusion coefficient is observed in the dynamic systems for short linkers containing 13 backbone atoms. The longest linkers with 71 backbone atoms show no difference between the permanent and dynamic networks. Linkers shorter than 11 backbone atoms, corresponding to a mesh size smaller than the penetrant small axis, diffusion is no longer observable on the experimental time scale, indicating a sharp cutoff attributed to the precise linkers and narrow mesh size distribution. The dynamic imine exchange time scales were compared to the diffusive hopping times of penetrants and indicate that exchange can occur during a diffusive displacement. Furthermore, these findings provide insights into the factors affecting penetrant transport in dense polymers and inspire the development of next-generation selective polymer membranes.

Diffusion↗
DOE OSTI2023

Model for Humidity-Mediated Diffusion on Aluminum Surfaces and Its Role in Accelerating Atmospheric Aluminum Corrosion

Bare aluminum metal surfaces are highly reactive, which leads to the spontaneous formation of a protective oxide surface layer. Because many subsequent corrosive processes are mediated by water, the structure and dynamics of water at the oxide interface are anticipated to influence corrosion kinetics. Here, using molecular dynamics simulations with a reactive force field, we model the behavior of aqueous aluminum metal ions in water adsorbed onto aluminum oxide surfaces across a range of ion concentrations and water film thicknesses corresponding to increasing relative humidity. We find that the structure and diffusivity of both the water and the metal ions depend strongly on the humidity of the environment and the relative height within the adsorbed water film. Aqueous aluminum ion diffusion rates in water films corresponding to a typical indoor relative humidity of 30% are found to be more than 2 orders of magnitude slower than self-diffusion of water in the bulk limit. Connections between metal ion diffusivity and corrosion reaction kinetics are assessed parametrically with a reductionist model based on a 1D continuum reaction–diffusion equation. Our results highlight the importance of incorporating the properties specific to interfacial water in predictive models of aluminum corrosion.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Influence of Rigidity–Hydration Coupling on Size-Dependent Diffusion in Hydrated Polymer Membranes

Selective ion transport in polymer membranes depends critically on how penetrant motion couples to polymer dynamics and hydration. Yet, the mechanistic interplay between polymer rigidity, water content, and penetrant size remains poorly understood, especially in the regime where the penetrant diameter, polymer Kuhn length, and correlation length are comparable. Here, we employ coarse-grained molecular dynamics simulations to systematically investigate penetrant diffusion in hydrated polymer networks across a broad range of water volume fractions, chain rigidities, and penetrant sizes. The results reveal a transition from a decoupled regime, where small penetrants diffuse nearly independently of polymer relaxation, to a coupled regime in which large penetrants require cooperative polymer motion for transport. Increasing polymer rigidity amplifies the sensitivity of diffusivity to hydration, particularly at low water content, leading to pronounced deviations from Stokes−Einstein scaling. Comparison with scaling theories and free-volume models shows that classical nanoparticle-based frameworks fail to capture this intermediate regime. To address this gap, we extend the Yasuda model to incorporate polymer rigidity through a single parameter that quantifies the dynamic contribution of chain stiffness to free-volume fluctuations. The resulting model collapses diffusivity data across all sizes, water contents, and rigidities, providing a unified description of penetrant transport in hydrated polymer matrices. Furthermore, these findings establish polymer rigidity as a key, tunable determinant of diffusion and offer a framework for interpreting size-dependent transport in ion-selective membranes.

diffusion↗
DOE OSTI2024

Surface Self-Diffusion Induced Sintering of Nanoparticles

Despite the critical role of sintering phenomena in constraining the long-term durability of nanosized particles, a clear understanding of nanoparticle sintering has remained elusive due to the challenges in atomically tracking the neck initiation and discerning different mechanisms. Through the integration of in situ transmission electron microscopy and atomistic modeling, this study uncovers the atomic dynamics governing the neck initiation of Pt–Fe nanoparticles via a surface self-diffusion process, allowing for coalescence without significant particle movement. Real-time imaging reveals that thermally activated surface morphology changes in individual nanoparticles induce significant surface self-diffusion. Further, the kinetic entrapment of self-diffusing atoms in the gaps between closely spaced nanoparticles leads to the nucleation and growth of atomic layers for neck formation. This surface self-diffusion-driven sintering process is activated at a relatively lower temperature compared to the classic Ostwald ripening and particle migration and coalescence processes. The fundamental insights have practical implications for manipulating the morphology, size distribution, and stability of nanostructures by leveraging surface self-diffusion processes.

36 MATERIALS SCIENCE↗