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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Ab-initio informed cluster dynamics simulation of self- and Xe diffusivity in uranium mononitride under irradiation

Uranium mononitride (UN) is one of the ceramic nuclear fuel alternatives to oxide fuels considered for light water reactor and advanced reactor designs, as it presents significant advantages such as high uranium density (better economics) and high thermal conductivity and melting point (increased safety). Self- and fission gas diffusivities need to be better understood, given that they influence key fuel performance phenomena like swelling and fission gas release. Recently, radiation enhanced diffusivity was investigated in UN by means of cluster dynamics simulations relying on empirical potential-based parameterizations, the reliability of which highly depends on the interatomic potential accuracy. Here, in this work, we refine this approach by determining, using ab-initio calculations, the properties of defect clusters containing vacancies, self-interstitials and Xe impurities. We also consider larger clusters than previous studies. The obtained dataset (formation enthalpies, entropies, and migration barriers) is used to parameterize a cluster dynamics model of mobile clusters, and to calculate the defect cluster concentrations under irradiation. This gives us access to the radiation enhanced self- and fission gas diffusivities. Although the resulting diffusivities are close to the values reported in the literature, we find important qualitative differences in the diffusion mechanisms. Capturing the correct mechanisms is crucial to properly describe the chemistry and fission rate dependence of the model.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Investigation of species-mass diffusion in binary-species boundary layers at high pressure using direct numerical simulations

Here, direct numerical simulations of single-species and binary-species temporal boundary layers at high pressure are performed with special attention to species-mass diffusion. The working fluids are nitrogen or a mixture of nitrogen and methane. Mean profiles and turbulent fluctuations of mass fraction show that their qualitative characteristics are different from those of streamwise velocity and temperature, due to the different boundary conditions. In a wall-parallel plane near the wall, the streamwise velocity and temperature have streaky patterns and the fields are similar. However, the mass fraction field at the same location is different from the streamwise velocity and temperature fields indicating that species-mass diffusion is not similar to the momentum and thermal diffusion. In contrast, at the centre and near the edge of the boundary layer, the mass fraction and temperature fields have almost the same pattern, indicating that the similarity between thermal and species-mass diffusion holds away from the wall. The lack of similarity near the wall is traced to the Soret effect that induces a temperature-gradient-dependent species-mass flux. As a result, a new phenomenon has been identified for a non-isothermal binary-species system – uphill diffusion, which in its classical isothermal definition can only occur for three or more species. A quadrant analysis for the turbulent mass flux reveals that near the wall the Soret effect enhances the negative contributions of the quadrants. Due to the enhancement of the negative contributions, small species-concentration fluid tends to be trapped near the wall.

42 ENGINEERING↗
DOE OSTI2025

Bio-Oil Impact on Water Diffusion and Durability of Bitumen: Influence of Aging and Salinity

Bio-oils derived from inexpensive biomass offer environmentally friendly options for modifying bitumen for improved durability or rejuvenating aged material, although their impact on bitumen moisture resistance can be mixed. A better understanding of water diffusion into bitumen and its interactions with bio-oil compounds would aid the development of effective bio-oils. In this study, the effect of water exposure on bio-oil-modified bitumen blends was examined by de-wetting contact angle measurements and differential FTIR spectroscopy. The bio-oils improved the anti-stripping behavior of the bitumen but also increased water absorption which could weaken cohesive strength. Short-term thermal aging increased water diffusion for all bitumen blends including the control, probably due to the presence of more oxidized compounds. Negative peaks of alkanes and polar groups in the differential FTIR data suggest the co-diffusion of surfactant molecules towards the bitumen-water interface. Basic pH increased de-wetting of some bitumen blends from silica, possibly by attacking the silica surface itself. The presence of salt in solution altered the bitumen surface composition through the formation of salt complexes with bitumen or bio-oil compounds. In particular, the formation of calcium-carboxylate complexes appeared to greatly improve anti-stripping effects of bio-oils. Water diffusion into most of the bitumen blends was insensitive to pH or salt concentration except for a few notable outliers. Identifying the compounds responsible for increasing or decreasing water diffusion in these outlier cases would be valuable developing future bio-oil formations that avoid or promote those compounds.

Hung, Albert M.↗
DOE OSTI2024

Effects of Nonequilibrium Atomic Structure on Ionic Diffusivity in LLZO: A Classical and Machine Learning Molecular Dynamics Study

To improve the performance of electrochemical devices, it is essential to understand the effects of nonequilibrium motifs in solids, such as grain boundaries, amorphous phases, and highly strained regions, on atomic-scale transport and stability. Molecular dynamics simulations are used to explore the combined effect of far-from-equilibrium atomic structures and the choice of interatomic potential on ionic diffusivity predictions for Li 7 La 3 Zr 2 O 12 (LLZO), a promising solid electrolyte for all-solid-state batteries. Amorphization and high strain are considered using both classical Buckingham interatomic potentials and machine learning force fields. Here we find that both crystalline expansion and amorphization tend to slow diffusion, although the different physical encodings in the two potentials impact the properties in different ways. We trace these variations to a combination of structural and transport factors, the contributions of which are deconvoluted computationally. Graph-based analysis reveals that the variations for amorphous LLZO arise from the connectivity of diffusion pathways within the predicted structures, which generally correlates with diffusivity and is notably higher for structures generated by the machine learning force fields. Our study provides additional insight into the relationship between atomic structure and diffusivity in LLZO, while also highlighting the need for care in choosing and validating potentials to simulate far from equilibrium structures.

25 ENERGY STORAGE↗
DOE OSTI2026

Diffusion of NiH on Ni(111) and of Ni with Adsorbed H on Ni(111) in the Context of Ni Coarsening in Solid Oxide Cells

In the Ni-based hydrogen electrode of a solid oxide cell (SOC), Ni coarsening is an important degradation mechanism and could be enhanced by NiH diffusing on the Ni particle surfaces. Here, in this work, the average lifetime and diffusion distance of NiH on Ni(111) are computed using density-functional theory and kinetic Monte Carlo methods. It is found that NiH is extremely short-lived and, thus, cannot promote coarsening in the SOC. Also, the diffusion of Ni on Ni(111) is shown to be at most slightly accelerated by H along NiH dissociation paths involving the movement of Ni with a nearby H. Based on this result, coarsening is not likely to be dramatically accelerated by the diffusion of Ni with H on Ni(111). However, support is provided for the experimental procedure of measuring the product of surface coverage and single-species diffusivity of Ni on Ni(111) in a hydrogen-rich atmosphere.

Ni coarsening↗
DOE OSTI2020

A Particle-Based Conditional Sampling Scheme for the Simulation of Transport in Fractured Rock With Diffusion Into Stagnant Water and Rock Matrix

In situ experiments and field-scale characterization studies have pointed out that, in fractured crystalline media, groundwater flow is highly channelized. This implies that, at the scale of a single fracture, only part of the fracture surface area is in contact with flowing water, while the rest of in-plane water is essentially stagnant and can be accessed by solutes via molecular diffusion. Despite their importance for contaminant retention, to date, there are no numerical or analytical approaches that could be used to assess the implication of stagnant water zones on solute transport in realistic large-scale Discrete Fracture Network-based models. In this study, we present an efficient and flexible algorithm for the simulation of transport in fractured rock with diffusion into stagnant water and rock matrix. The algorithm is a generalization of a previously developed numerical framework for time domain particle tracking in sparsely fractured rock. The key of the generalization is that total time in fracture (τ f ) is first evaluated using a Monte Carlo sampling and then a second sampling is performed conditioned on τ f . The algorithm has been successfully validated against existing independent solutions and the implication of diffusion into stagnant water and secondary diffusion into the matrix has been assessed for a realistic modeling scenario. The results show that, due to diffusion into stagnant water, contaminants are more strongly retarded. This increased retention is more significant for sorbing species, as a larger number of sorption sites is accessible. A high sensitivity to the flowing channel/stagnant water zone geometry has also been observed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Diffuse auroral precipitation driven by lower-band chorus second harmonics

Diffuse aurora at the Earth’s high latitude regions is mainly caused by the low-energy (0.1–30 keV) electron precipitation which carries the major energy flux into the nightside upper atmosphere. Previous studies have demonstrated that combined scattering by the upper- and lower- band chorus waves acts as the dominant cause of diffuse auroral precipitation, but that is not necessarily the case as these two types of waves do not always occur simultaneously, with the lower-band more often. Here we report that the lower-band chorus satisfying the preferred condition can generate their second harmonics so as to trigger the diffuse auroral electron precipitation. We find that the lower-band chorus alone can only cause the precipitation of electrons greater than 4 keV, while the self-consistently generated second harmonic is weak but still able to result in the electron precipitation below 4 keV. The combined effect of those modes results in the observed pancake electron distributions and the diffuse aurora. Our results clearly demonstrate an alternative but universal mechanism of chorus-driven diffuse aurora in the Earth, which may also apply to the auroral formation in other planetary magnetospheres.

79 ASTRONOMY AND ASTROPHYSICS↗
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↗