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 1,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Grain Boundary Effects on Radiation Damage and Tritium Diffusion in Li–Al–O Ceramics from Molecular Dynamics and Experiments

Understanding irradiation damage and tritium transport in LiAlO2 ceramics is essential for their deployment in tritium-producing burnable absorber rods (TPBARs). Grain boundaries (GBs) play an important role in governing radiation response and tritium transport in LiAlO2 and its secondary phase LiAl5O8. Using molecular dynamics simulations, we investigated defect evolution and tritium diffusion during displacement cascades in single-crystal and bicrystal LiAlO2 and LiAl5O8 at 600 K. The results reveal that GBs suppress the damage in the LiAlO2 bicrystals as compared to single crystals, via interstitial emission. The effect is more pronounced for Li defects, than for Al and O. In LiAl5O8 the damage reduction due to introduction of GB is less pronounced overall. Tritium diffusion coefficients at GBs are higher by factors of 2–10 relative to diffusion inside grain bulk. The effect is particularly pronounced in LiAlO2, where S5, S17, and S25 GBs promote rapid tritium transport, whereas LiAl5O8 exhibits slower diffusion due to reduced free volume at its GBs. These findings evince a trade-off that while GBs mitigate radiation damage by absorbing interstitials, they simultaneously provide fast pathways for tritium migration which is not desirable for TPBAR application. The mechanistic insights gained here establish a foundation for microstructural design strategies to balance radiation tolerance and tritium retention in ceramic breeder materials.

Roy, Ankit↗
DOE OSTI2018

Fluid-enhanced surface diffusion controls intraparticle phase transformations

Phase transformations driven by compositional change require mass flux across a phase boundary. In some anisotropic solids, however, the phase boundary moves along a non-conductive crystallographic direction. One such material is Li X FePO 4 , an electrode for lithium-ion batteries. With poor bulk ionic transport along the direction of phase separation, it is unclear how lithium migrates during phase transformations. Here, we show that lithium migrates along the solid/liquid interface without leaving the particle, whereby charge carriers do not cross the double layer. X-ray diffraction and microscopy experiments as well as ab initio molecular dynamics simulations show that organic solvent and water molecules promote this surface ion diffusion, effectively rendering Li X FePO 4 a three-dimensional lithium-ion conductor. Phase-field simulations capture the effects of surface diffusion on phase transformation. Lowering surface diffusivity is crucial towards supressing phase separation. In conclusion, this work establishes fluid-enhanced surface diffusion as a key dial for tuning phase transformation in anisotropic solids.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Hydrogen self-diffusion in single crystal olivine and electrical conductivity of the Earth’s mantle

Nominally anhydrous minerals formed deep in the mantle and transported to the Earth’s surface contain tens to hundreds of ppm wt H 2 O, providing evidence for the presence of dissolved water in the Earth’s interior. Even at these low concentrations, H 2 O greatly affects the physico-chemical properties of mantle materials, governing planetary dynamics and evolution. The diffusion of hydrogen (H) controls the transport of H 2 O in the Earth’s upper mantle, but is not fully understood for olivine ((Mg, Fe) 2 SiO 4 ) the most abundant mineral in this region. Here we present new hydrogen self-diffusion coefficients in natural olivine single crystals that were determined at upper mantle conditions (2 GPa and 750–900 °C). Hydrogen self-diffusion is highly anisotropic, with values at 900 °C of 10 -10.9 , 10 -12.8 and 10 -11.9 m 2 /s along [100], [010] and [001] directions, respectively. Combined with the Nernst-Einstein relation, these diffusion results constrain the contribution of H to the electrical conductivity of olivine to be σH = 10 2.12 S/m·C H2O ·exp -187kJ/mol/(RT) . Comparisons between the model presented in this study and magnetotelluric measurements suggest that plausible H 2 O concentrations in the upper mantle (≤250 ppm wt) can account for high electrical conductivity values (10 -2 –10 -1 S/m) observed in the asthenosphere.

58 GEOSCIENCES↗
DOE OSTI2018

Microscopic diffusion processes measured in living planarians

Living planarian flatworms were probed using quasielastic neutron scattering to measure, on the pico-to-nanosecond time scale and nanometer length scale, microscopic diffusion of water and cell constituents in the planarians. Measurable microscopic diffusivities were surprisingly well defined in such a complex system as living animals. The overall variation in the microscopic diffusivity of cell constituents was found to be far lower than the variation in the microscopic diffusivity of water in planarians in a temperature range of 284.5 to 304.1K.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

Mineral reaction kinetics constrain the length scale of rock matrix diffusion

Mass transport by aqueous fluids is a dynamic process in shallow crustal systems, redistributing nutrients as well as contaminants. Rock matrix diffusion into fractures (void space) within crystalline rock has been postulated to play an important role in the transient storage of solutes. The reacted volume of host rock involved, however, will be controlled by fluid-rock reactions. Here we present the results of a study which focusses on defining the length scale over which rock matrix diffusion operates within crystalline rock over timescales that are relevant to safety assessment of radioactive and other long-lived wastes. Through detailed chemical and structural analysis of natural specimens sampled at depth from an active system (Toki Granite, Japan), we show that, contrary to commonly proposed models, the length scale of rock matrix diffusion may be extremely small, on the order of centimetres, even over timescales of millions of years. This implies that in many cases the importance of rock matrix diffusion will be minimal. Additional analyses of a contrasting crystalline rock system (Carnmenellis Granite, UK) corroborate these results.

58 GEOSCIENCES↗
DOE OSTI2020

Carbon nanotube porin diffusion in mixed composition supported lipid bilayers

Carbon nanotube porins (CNTPs), short pieces of carbon nanotubes capable of self-inserting into a lipid bilayer, represent a simplified model of biological membrane channels. We have used high-speed atomic force microscopy (HS-AFM) and all-atom molecular dynamics (MD) simulations to study the behavior of CNTPs in a mixed lipid membrane consisting of DOPC lipid with a variable percentage of DMPC lipid added to it. HS-AFM data reveal that the CNTPs undergo diffusive motion in the bilayer plane. Motion trajectories extracted from the HS-AFM movies indicate that CNTPs exhibit diffusion coefficient values broadly similar to values reported for membrane proteins in supported lipid bilayers. The data also indicate that increasing the percentage of DMPC leads to a marked slowing of CNTP diffusion. MD simulations reveal a CNTP-lipid assembly that diffuses in the membrane and show trends that are consistent with the experimental observations.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

The role of geography in the complex diffusion of innovations

The urban–rural divide is increasing in modern societies calling for geographical extensions of social influence modelling. Improved understanding of innovation diffusion across locations and through social connections can provide us with new insights into the spread of information, technological progress and economic development. In this work, we analyze the spatial adoption dynamics of iWiW, an Online Social Network (OSN) in Hungary and uncover empirical features about the spatial adoption in social networks. During its entire life cycle from 2002 to 2012, iWiW reached up to 300 million friendship ties of 3 million users. We find that the number of adopters as a function of town population follows a scaling law that reveals a strongly concentrated early adoption in large towns and a less concentrated late adoption. We also discover a strengthening distance decay of spread over the life-cycle indicating high fraction of distant diffusion in early stages but the dominance of local diffusion in late stages. The spreading process is modelled within the Bass diffusion framework that enables us to compare the differential equation version with an agent-based version of the model run on the empirical network. Although both model versions can capture the macro trend of adoption, they have limited capacity to describe the observed trends of urban scaling and distance decay. We find, however that incorporating adoption thresholds, defined by the fraction of social connections that adopt a technology before the individual adopts, improves the network model fit to the urban scaling of early adopters. Controlling for the threshold distribution enables us to eliminate the bias induced by local network structure on predicting local adoption peaks. Finally, we show that geographical features such as distance from the innovation origin and town size influence prediction of adoption peak at local scales in all model specifications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Chemical instability leads to unusual chemical-potential-independent defect formation and diffusion in perovskite solar cell material CH 3 NH 3 PbI 3

Methylammonium (MA) lead triiodide (MAPbI 3 ) has recently emerged as a promising solar cell material. But, MAPbI3 is known to have chemical instability, i.e., MAPbI3 is prone to decomposition into MAI and PbI 2 even at moderate temperatures (e.g. 330 K). Here, we show that the chemical instability, as reflected by the calculated negligible enthalpy of formation of MAPbI 3 (with respect to MAI and PbI 2 ), has an unusual and important consequence for defect properties, i.e., defect formation energies in low-carrier-density MAPbI 3 are nearly independent of the chemical potentials of constituent elements and thus can be uniquely determined. This allows straightforward calculations of defect concentrations and the activation energy of ionic conductivity (the sum of the formation energy and the diffusion barrier of the charged mobile defect) in MAPbI 3 . Furthermore, the calculated activation energy for ionic conductivity due to V$+\atop{1}$ diffusion is in excellent agreement with the experimental values, which demonstrates unambiguously that V$+\atop{1}$ is the dominant diffusing defect and is responsible for the observed ion migration and device polarization in MAPbI3 solar cells. The calculated low formation energy of a Frenkel pair (V$+\atop{1}$ -I$-\atop{i}$ and low diffusion barriers of V$+\atop{1}$ and Image I$-\atop{i}$ suggest that the iodine ion migration and the resulting device polarization may occur even in single-crystal devices and grain-boundary-passivated polycrystalline thin film devices (which were previously suggested to be free from ion-migration-induced device polarization), leading to device degradation. Moreover, the device polarization due to the Frenkel pair (which has a relatively low concentration) may take a long time to develop and thus may avoid the appearance of the current–voltage hysteresis at typical scan rates.

14 SOLAR ENERGY↗
DOE OSTI2017

Diffusive confinement of free radical intermediates in the OH radical oxidation of semisolid aerosols

Multiphase chemical reactions (gas + solid/liquid) involve a complex interplay between bulk and interface chemistry, diffusion, evaporation, and condensation. Reactions of atmospheric aerosols are an important example of this type of chemistry: the rich array of particle phase states and multiphase transformation pathways produce diverse but poorly understood interactions between chemistry and transport. Their chemistry is of intrinsic interest because of their role in controlling climate. Their characteristics also make them useful models for the study of principles of reactivity of condensed materials under confined conditions. Previously, we have reported a computational study of the oxidation chemistry of a liquid aliphatic aerosol. In this study, we extend the calculations to investigate nearly the same reactions at a semisolid gas-aerosol interface. A reaction-diffusion model for heterogeneous oxidation of triacontane by hydroxyl radicals (OH) is described, and its predictions are compared to measurements of aerosol size and composition, which evolve continuously during oxidation. Our results are also explicitly compared to those obtained for the corresponding liquid system, squalane, to pinpoint salient elements controlling reactivity. The diffusive confinement of the free radical intermediates at the interface results in enhanced importance of a few specific chemical processes such as the involvement of aldehydes in fragmentation and evaporation, and a significant role of radical-radical reactions in product formation. The simulations show that under typical laboratory conditions semisolid aerosols have highly oxidized nanometer-scale interfaces that encapsulate an unreacted core and may confer distinct optical properties and enhanced hygroscopicity. This highly oxidized layer dynamically evolves with reaction, which we propose to result in plasticization. The validated model is used to predict chemistry under atmospheric conditions, where the OH radical concentration is much lower. The oxidation reactions are more strongly influenced by diffusion in the particle, resulting in a more liquid-like character.

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

Hybrid DFT investigation of the energetics of Mg ion diffusion in α-MoO 3

We report that rechargeable batteries that utilize divalent Mg ions as the charge carrier species can in principle achieve substantially greater volumetric energy densities than conventional Li-ion batteries. One significant impediment to the development of commercially viable Mg-ion batteries is the slow rate of Mg ion diffusion through otherwise promising cathode materials. Accurate prediction of the activation energies associated with this diffusion process using density functional theory (DFT) is especially challenging due to self-interaction errors intrinsic to DFT that lead to over-delocalization of the d-electrons. One effective but highly computationally demanding approach to reducing self-interaction errors is the use of hybrid functionals, which incorporate a fraction of exact Hartree–Fock exchange. In this work, we assess the effects of exact exchange on computed activation energies for ion diffusion in one potential cathode material, α-MoO 3 . In contrast to previous studies that primarily utilize non-hybrid functionals, we perform nudged elastic band calculations in which the nuclear coordinates are fully converged using both hybrid functionals and k-point sampling. Additionally, it is found that while non-hybrid functionals indicate the existence of thermodynamically accessible channels for bulk Mg ion diffusion in all three dimensions, hybrid functionals predict that some of these channels are largely inaccessible under typical charge/discharge conditions. Furthermore, it is demonstrated that certain commonly used approximations for incorporating the effects of Hartree–Fock exchange are inadequate for this system, including DFT+U calculations and the use of single-point hybrid calculations using atomic positions obtained using non-hybrid functionals.

25 ENERGY STORAGE↗
DOE OSTI2019

Diffusive ferromagnetic roller gas

An ensemble of actively rotating ferromagnetic particles is used to realize an active roller gas. Here, we investigate the diffusive properties of such a gas in experiments and simulations. We reveal that ferromagnetic rollers demonstrate a normal (Fickian) diffusion with a characteristic linear growth of the mean-squared displacement, while statistics of displacements stay non-Gaussian. At short times the system has a bimodal distribution of the displacements that transitions with time to a quasi-Gaussian distribution (Gaussian core with overpopulated tails) for a range of studied particle number densities. Inert particles introduced into the active roller gas exhibit similar diffusive behavior. Here, the results provide insights into diffusive properties of active colloidal systems with activity originating from spinning degrees of freedom.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Diffusion characteristics of water molecules in a lamellar structure formed by triblock copolymers

The distribution and diffusion of water molecules are playing important roles in determining self-assembly and transport properties of polymeric systems. Small-angle neutron scattering (SANS) experiments and molecular dynamics (MD) simulation have been applied to understand the distribution of water molecules and their dynamics in the lamellar membrane formed by Pluronic L62 block copolymers. Penetration of water molecules into the polyethylene oxide (PEO) layers of the membranes has been estimated using scattering length density (SLD) profiles obtained from SANS measurements, which agree well with the molecular distribution observed from MD simulations. The water diffusion coefficient at different regions of the lamellar membrane was further investigated using MD simulation. The diffusion characteristic shows a transition from normal to anomalous diffusion as the position of the water molecule changes from the bulk to PEO and to the polypropylene oxide (PPO) layer. Here we find that water molecules within the PEO or PPO layers follow subdiffusive dynamics, which can be interpreted by the model of fractional Brownian motion.

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

Communication: Diverse nanoscale cluster dynamics: Diffusion of 2D epitaxial clusters

The dynamics of nanoscale clusters can be distinct from macroscale behavior described by continuum formalisms. For diffusion of 2D clusters of N atoms in homoepitaxial systems mediated by edge atom hopping, macroscale theory predicts simple monotonic size scaling of the diffusion coefficient, D N ~ N –β , with β = 3/2. However, modeling for nanoclusters on metal(100) surfaces reveals that slow nucleation-mediated diffusion displaying weak size scaling β < 1 occurs for “perfect” sizes N p = L 2 and L(L+1) for integer L = 3,4,… (with unique square or near-square ground state shapes), and also for N p +3, N p +4,…. In contrast, fast facile nucleation-free diffusion displaying strong size scaling β ≈ 2.5 occurs for sizes N p +1 and N p +2. D N versus N oscillates strongly between the slowest branch (for N p +3) and the fastest branch (for N p +1). All branches merge for N = O(10 2 ), but macroscale behavior is only achieved for much larger N = O(10 3 ). Here, this analysis reveals the unprecedented diversity of behavior on the nanoscale.

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

A diffusion-limited reaction model for self-propagating Al/Pt multilayers with quench limits

A diffusion-limited reaction model was calibrated for Al/Pt multilayers ignited on oxidized silicon, sapphire, and tungsten substrates, as well as for some Al/Pt multilayers ignited as free-standing foils. The model was implemented in a finite element analysis code and used to match experimental burn front velocity data collected from several years of testing at Sandia National Laboratories. Moreover, both the simulations and experiments reveal well-defined quench limits in the total Al + Pt layer (i.e., bilayer) thickness. At these limits, the heat generated from atomic diffusion is insufficient to support a self-propagating wave front on top of the substrates. Quench limits for reactive multilayers are seldom reported and are found to depend on the thermal properties of the individual layers. Here, the diffusion-limited reaction model is generalized to allow for temperature- and composition-dependent material properties, phase change, and anisotropic thermal conductivity. Utilizing this increase in model fidelity, excellent overall agreement is shown between the simulations and experimental results with a single calibrated parameter set. However, the burn front velocities of Al/Pt multilayers ignited on tungsten substrates are over-predicted. Finally, possible sources of error are discussed and a higher activation energy (from 41.9 kJ/mol.at. to 47.5 kJ/mol.at.) is shown to bring the simulations into agreement with the velocity data observed on tungsten substrates. Finally, this higher activation energy suggests an inhibited diffusion mechanism present at lower heating rates.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Comprehensive evaluation and parametric sensitivity of interatomic potential models for diffusion kinetics of Cr 2 O 3 in molecular dynamics

While molecular dynamics (MD) has proven to be a promising approach to investigate the diffusion properties, the grand challenge resides in evaluating potential model parameters to accurately replicate experimentally measured properties. The Buckingham potential model with Columbic interaction is widely employed in MD simulations of chromia (Cr 2 O 3 ) systems, as it allows for reasonable computational cost and accuracy. However, considering the well-known limitation of classical potential models in simultaneous reproduction of various physical phenomena, further comprehensive evaluation of the potential is required for calculation of diffusion properties. In this study, we benchmark the performance of three different Buckingham models with the experimental data by calculating structural, thermodynamic, and mechanical properties of defect-free Cr 2 O 3 , and diffusion properties of Cr 2 O 3 with vacancy defects. Available Buckingham models display limited accuracies, consolidating the necessity of retraining the potential parameters for all properties impacting the diffusion dynamics. Oversimplification in parameterization procedures is suggested to impede the universal performance in property reproduction. Furthermore, this research also demonstrates effective guidelines for choosing a proper parameter set of current Buckingham potential for MD simulation with Cr 2 O 3 depending on properties and for potential reparameterization.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Solid-phase reactive inter-diffusion of Mg/B multilayers

A better understanding of the reaction of Mg and B in the solid-phase regime is needed for the low-temperature synthesis of MgB 2 films. Here in this paper, we study the kinetics of reactive inter-diffusion of Mg and B multilayers on glassy carbon substrates in the temperature range of 400-650°C. Results show that, at these temperatures, inter-diffusion is characterized by a single activation energy of ~0.45eV. The formation of the superconducting MgB 2 phase with critical temperatures of 25–31 K occurs at reaction temperatures of 450°C and above, with the rate of inter-diffusion obeying a power law with a kinetic exponent of ~0.3. This suggests that rate-limiting processes are the nucleation and growth of MgB 2 grains rather than diffusion and interfacial reactions. Implications of these results to the low-temperature synthesis of MgB 2 films are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

The impact of nanoparticle softness on its tracer diffusion coefficient in all polymer nanocomposites

The diffusion of nanoparticles in a polymer matrix is an area of current interest. However, a complete understanding is still limited as it is often difficult to quantify the much slower motion of nanoparticles in a polymer matrix. To combat this problem, we have developed a protocol to measure the diffusion coefficient of soft nanoparticles in a linear polymer matrix. Recently developed synthetic control over soft nanoparticle structures combined with this protocol provides a pathway to separately elucidate the effects of the molecular weight and nanoparticle softness on its diffusive behavior. These results indicate that the nanoparticle softness and deformability dictate its motion. Increasing the cross-linking density of the nanoparticle for all molecular weights increases its hardness and suppresses its motion in the linear matrix. Additionally, the nanoparticle molecular weight dependence deviates from the exponential dependence for star polymers suggesting that these nanoparticles benefit from the cooperative motion of the matrix to open pathways for the nanoparticle. Finally, comparison of these experimentally determined values to the Stokes–Einstein theory demonstrates that the nanoparticles diffuse much slower than a hard sphere. This is interpreted to indicate that there exist additional interactions between the nanoparticle and polymer matrix that are not captured by Stokes–Einstein, including threading or entanglement of the linear chain with the nanoparticle.

Rostom, Sahar↗
DOE OSTI2023

Diffusion of Sn donors in β-Ga 2 O 3

Diffusion of the n-type dopant Sn in β-Ga 2 O 3 is studied using secondary-ion mass spectrometry combined with hybrid functional calculations. The diffusion of Sn from a Sn-doped bulk substrate with surface orientation (001) into an epitaxial layer is observed after heat treatments in the temperature range of 1050–1250 °C. Calculated formation energies of Sn-related and intrinsic defects show that the migration of Sn is mediated by Ga vacancies ( V Ga ) through the formation and dissociation of intermittent mobile V Ga Sn Ga complexes. The evolution of the Sn concentration vs depth profiles after heat treatments can be well described by a reaction–diffusion model. Using model parameters guided by the hybrid functional calculations, we extract a V Ga Sn Ga complex migration barrier of 3.0 ± 0.4 eV with a diffusion coefficient of 2 × 10 –1 cm 2 /s. The extracted migration barrier is consistent with our theoretical predictions using the nudged elastic band method, which shows migration barriers of 3.42, 3.15, and 3.37 eV for the [100], [010], and [001] directions, respectively.

36 MATERIALS SCIENCE↗