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

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

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

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

Diffuse Reflectance Spectroscopy and Principal Component Analysis to Retrospectively Determine Production History of Plutonium Dioxide

Diffuse reflectance spectroscopy measurements in the shortwave infrared (930–1600 nm) spectral region were acquired for Pu 2 (C 2 O 4 ) 3 •9H 2 O and its thermal decomposition product, PuO 2 . We analyzed a total of eight PuO 2 samples that were produced at different calcination temperatures (300, 350, 450, 525, 600, 675, 750, and 900 °C). Our goal was to identify spectroscopic fingerprints that could be used to gain retrospective information regarding the production parameters of these important nuclear compounds. The diffuse reflectance spectrum of Pu 2 (C 2 O 4 ) 3 •9H 2 O features several broad bands that currently preclude detailed analysis. However, all PuO 2 samples produced relatively sharp spectral features that got sharper and more intense for samples that were produced at higher calcination temperatures. The electronic band observed at 1433 nm in the diffuse reflectance spectra of PuO 2 was found to be a sensitive indicator of crystallinity; a result that is corroborated by ancillary Raman spectroscopy measurements. Principal component analysis of diffuse reflectance spectra was able to clearly rank and categorize PuO 2 samples based on the calcination temperature that was employed during their production. Thus, we show herein that important retrospective information pertaining to the process history of PuO 2 can be gained through the relatively simplistic combination of diffuse reflectance spectroscopy and principal component analysis. This discovery presents a new method for determining the provenance and process history of PuO 2 and should have an impact in the fields of nuclear forensics and nuclear nonproliferation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2018

Modeling Radionuclide Releases from TRISO Particles by Simultaneous Diffusion Through and Corrosion of the Silicon Carbide Barrier Layer

TRISO nuclear fuel particles that are less than 1 mm in diameter are designed with multiple barrier layers to retain fission products both during reactor operations and for long-term geological disposal. The primary barrier is a 35 μm thick silicon carbide (SiC-a highly impermeable semi-metal) layer for which data are available on the diffusion of short-lived fission products at high temperatures (> 1000 °C). However, for a geological repository, this layer may contact brine and hence corrode even at ambient temperatures. As an initial approach to assess the effectiveness of the SiC barrier for geological repositories, ranges of fission product diffusivities and corrosion rates for SiC are modeled concurrently with the simultaneous effect of radioactive decay. Using measured corrosion rates of SiC, if the diffusivity is more than about 10 -20 m 2 /s, fission product releases may occur before the SiC barrier has corroded to the point of breach. For diffusivities less than about 10 -21 m 2 /s there may not be significant diffusional releases prior to SiC barrier removal/breach by corrosion. This work shows the importance of estimating diffusivities in SiC at geological repository temperatures, and highlights the relevance of evaluating the porosity/permeability evolution of the SiC layer in a geologic environment.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2019

Neptunium(IV) Diffusion through Bentonite Clay

DOE Fellow, Silvina A. Di Pietro, completed a 10-week internship at Lawrence Livermore National Laboratory (LLNL) in Livermore, California. The objective of the study was to determine the diffusion profile of neptunium (Np) in bentonite clay (potential geo-engineered barrier for nuclear waste repositories). Under the mentorship of Drs. Claudia Joseph (Karlsruhe Institute of Technology, Germany) and Mavrik Zavarin, she prepared a series of core abrasion peeling samples. This procedure was a continuation of a Np diffusion experiment conducted by Dr. Joseph during her post-doc residency at LLNL in the year 2016. The second part of the experiment (from now on will be referred as abrasive peeling) focuses on understanding (1) the Np movement through the clay core and (2) the influence of temperature (25°C, 65°C, and 85°C) for the Np transport across the clay core under anaerobic conditions. End of 2016, the diffusion core samples were embedded in epoxy resin to stabilize them for long-term. The abraded core powder samples vary in material (epoxy vs. clay), Np activity and concentration, temperature exposed during the diffusion experiment, and mass, or amount of clay peeled off. Throughout the course of the 10-week internship program at LLNL, the experiments implemented and discussed in this report aimed to answer the questions: (1) How is the Np(IV) diffusion profile through bentonite influenced by the temperature? and (2) Will Np(IV) migrate through the geoengineered barrier of a nuclear waste repository? In addition, Ph.D. candidate Di Pietro received invaluable guidance in experimental design and development that will benefit her ongoing research endeavors.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

DRIFT: Diffusivity Regulation of Isotopes by nanoconFinementT

Spontaneous isotope fractionation has been reported under nanoconfinement conditions in naturally occurring systems, but the origin of this phenomena is currently unknown. Two existing hypotheses have been proposed, one based on changes in the solvation environment of the isotopes that reduces the non-mass dependent hydrodynamics contribution to diffusion. The other is that isotopes have mass-dependent surface adsorption, varying their total diffusion through nanoconfined channels. To investigate these hypotheses, benchtop experiments, nuclear magnetic resonance (NMR) spectroscopy, and molecule scale modeling were applied. Classical molecular dynamics simulations identified that the Na + and Cl - hydration shells across the three different salt solutions ( 22 Na 35 Cl, 23 Na 35 Cl, 24 Na 35 Cl) did not vary as a function of the Na + isotope, but that there was a significant pore size effect, with larger hydration shells at larger pore sizes. Additionally, while total adsorption times did not vary as a function of the Na + isotope or pore size, the free ion concentration, or those adsorbed on the surface for <5% of the simulation time did exhibit isotope dependence. Experimentally, challenges occurred developing a repeatable experiment, but NMR characterization of water diffusion rates through ordered alumina membranes was able to identify the existence of two distinct water environments associated with water inside and outside the pore. Further NMR studies could be used to confirm variation in hydration shells and diffusion rates of dissolved ions in water. Ultimately, mass-dependence adsorption is a primary driver of variations in isotope diffusion rates, rather than variation in hydration shells that occur under nanoconfinement.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Model Parameter Development for Complex Materials: Species-Specific Diffusion Barriers in 316 Stainless Steel from Systematic DFT Calculations

Vacancy-mediated diffusion barriers in 316 stainless steel have been systematically calculated using density functional theory to provide essential parameters for mesoscale microstructure evolution models. A statistical sampling approach employing 210 nudged elastic band calculations across multiple special quasi-random structures captures the effects of local chemical environments in this concentrated alloy. The computational methodology addresses challenges specific to chemically disordered systems, including proper magnetic treatment throughout multi-step calculations and validation against experimental structural properties. The calculated activation barriers reveal clear species-dependent diffusion behavior with the hierarchy Ni >> Fe ˜ Cr >> Mo. Nickel exhibits the highest barriers (0.74–1.31 eV, mean 1.045 eV), confirming its role as the slowest-diffusing major component. Iron and chromium show similar moderate barriers averaging 0.587 eV and 0.522 eV, respectively. Remarkably, molybdenum demonstrates exceptionally low barriers (0.12–0.28 eV, mean 0.194 eV), suggesting much higher mobility than previously recognized and potentially significant implications for precipitation kinetics and microstructure evolution. The barrier ranges remain consistent across different 316 SS compositions, supporting parameter transferability for modeling applications. The overall mean barrier of 0.64 eV provides a practical approximation for phase field simulations, while species-specific values enable detailed treatments of diffusion-controlled processes. This systematic approach establishes a validated framework for generating diffusion parameters in other concentrated alloys where experimental data are limited, while providing the first systematic set of species-specific barriers for predictive modeling of 316 stainless steel microstructure evolution.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Atomic Diffusivities of Yttrium, Titanium and Oxygen Calculated by Ab Initio Molecular Dynamics in Molten 316L Oxide-Dispersion-Strengthened Steel Fabricated via Additive Manufacturing

Oxide-dispersion-strengthened (ODS) steels have long been viewed as a prime solution for harsh environments. However, conventional manufacturing of ODS steels limits the final product geometry, is difficult to scale up to large components, and is expensive due to multiple highly involved, solid-state processing steps required. Additive manufacturing (AM) can directly incorporate dispersion elements (e.g., Y, Ti and O) during component fabrication, thus bypassing the need for an ODS steel supply chain, the scale-up challenges of powder processing routes, the buoyancy challenges associated with casting ODS steels, and the joining issues for net-shape component fabrication. In the AM process, the diffusion of the dispersion elements in the molten steel plays a key role in the precipitation of the oxide particles, thereby influencing the microstructure, thermal stability and high-temperature mechanical properties of the resulting ODS steels. In this work, the atomic diffusivities of Y, Ti, and O in molten 316L stainless steel (SS) as functions of temperature are determined by ab initio molecular dynamics simulations. The latest Vienna Ab initio Simulation Package (VASP) package that incorporates an on-the-fly machine learning force field for accelerated computation is used. At a constant temperature, the time-dependent coordinates of the target atoms in the molten 316L SS were analyzed in the form of mean square displacement in order to obtain diffusivity. The values of the diffusivity at multiple temperatures are then fitted to the Arrhenius form to determine the activation energy and the pre-exponential factor. Given the challenges in experimental measurement of atomic diffusivity at such high temperatures and correspondingly the lack of experimental data, this study provides important physical parameters for future modeling of the oxide precipitation kinetics during AM process.

Chemistry↗
DOE OSTI2020

Diffusion coefficients in the envelopes of white dwarfs

The diffusion of elements is a key process in understanding the unusual surface composition of white dwarfs and their spectral evolution. The diffusion coefficients of Paquette et al. have been widely used to model diffusion in white dwarfs. In this work, we perform new calculations of the coefficients of interdiffusion and ionic thermal diffusion with (1) a more advanced model that uses a recent modification of the calculation of the collision integrals that is more suitable for the partially ionized, partially degenerate, and moderately coupled plasma and (2) classical molecular dynamics. The coefficients are evaluated for silicon and calcium in white dwarf envelopes of hydrogen and helium. A comparison of our results with Paquette et al. shows that the latter systematically underestimates the coefficient of interdiffusion yet provides reliable estimates for the relatively weakly coupled plasmas found in nearly all types of stars, as well as in white dwarfs with hydrogen envelopes. In white dwarfs with cool helium envelopes (${T}_{\mathrm{eff}}$ < 15,000 K), the difference grows to more than a factor of two. We also explored the effect of the ionization model used to determine the charges of the ions and found that it can be a substantial source of discrepancy between different calculations. Finally, we consider the relative diffusion timescales of Si and Ca in the context of the pollution of white dwarf photospheres by accreted planetesimals and find factor of ≳3 differences between calculations based on Paquette et al. and our model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Probabilistic estimation of depth-resolved profiles of soil thermal diffusivity from temperature time series

Abstract. Improving the quantification of soil thermal and physical properties is key to achieving a better understanding and prediction of soil hydro-biogeochemical processes and their responses to changes in atmospheric forcing. Obtaining such information at numerous locations and/or over time with conventional soil sampling is challenging. The increasing availability of low-cost, vertically resolved temperature sensor arrays offers promise for improving the estimation of soil thermal properties from temperature time series, and the possible indirect estimation of physical properties. Still, the reliability and limitations of such an approach need to be assessed. In the present study, we develop a parameter estimation approach based on a combination of thermal modeling, sliding time windows, Bayesian inference, and Markov chain Monte Carlo simulation to estimate thermal diffusivity and its uncertainty over time, at numerous locations and at an unprecedented vertical spatial resolution (i.e., down to 5 to 10 cm vertical resolution) from soil temperature time series. We provide the necessary framework to assess under which environmental conditions (soil temperature gradient, fluctuations, and trend), temperature sensor characteristics (bias and level of noise), and deployment geometries (sensor number and position) soil thermal diffusivity can be reliably inferred. We validate the method with synthetic experiments and field studies. The synthetic experiments show that in the presence of median diurnal fluctuations ≥ 1.5 ∘C at 5 cm below the ground surface, temperature gradients > 2 ∘C m−1, and a sliding time window of at least 4 d the proposed method provides reliable depth-resolved thermal diffusivity estimates with percentage errors ≤ 10 % and posterior relative standard deviations ≤ 5 % up to 1 m depth. Reliable thermal diffusivity under such environmental conditions also requires temperature sensors to be spaced precisely (with accuracy to a few millimeters), with a level of noise ≤ 0.02 ∘C, and with a bias defined by a standard deviation ≤ 0.01 ∘C. Finally, the application of the developed approach to field data indicates significant repeatability in results and similarity with independent measurements, as well as promise in using a sliding time window to estimate temporal changes in soil thermal diffusivity, as needed to potentially capture changes in bulk density or water content.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Probabilistic estimation of depth-resolved profiles of soil thermal diffusivity from temperature time series: Supporting Data

This dataset consists of soil temperature time series that were used to estimate soil thermal diffusivity and its uncertainty trough the probabilistic modelling approach developed and presented in the article named "Probabilistic estimation of depth-resolved profiles of soil thermal diffusivity from temperature time series" and published in Earth Surface Dynamics. There are two compressed (.zip) files that contains synthetic (Synthetic_soiltemp_Teller.zip) and field (Field_soiltemp_Teller.zip) data. There is one MATLAB file that requires MATLAB to execute but any text editor can open it. The Synthetic_soiltemp_Teller.zip file includes 5 comma-delimited data files (.csv) each of which contains soil temperature time series generated through forward modeling (i.e., heat-conduction process in a heterogeneous medium using an explicit finite difference method) to mimic various types of temperature gradients, trend and fluctuations. These synthetic soil temperatures were used to investigate the impact of different environmental conditions on the uncertainty of thermal diffusivity estimates. The Field_soiltemp_Teller.zip file contains 28 comma-delimited data files (.csv) out of which (a) 27 files includes soil temperature time series recorded from 27 temperature probes located in a site along Teller Road about 40 km northwest of Nome, Alaska (64.72°N, 165.94°W), (b) one includes the name and coordinates of the 27 probes. These field soil temperatures were used to infer soil thermal diffusivity at numerous locations and depths in a discontinuous permafrost environment, and to evaluate the links between the estimated soil thermal diffusivity values and soil physical properties. The comma-delimited data files (.csv) of the synthetic and field soil temperature time series includes date and time (UTC) in the first column and soil temperature from 5 cm below the ground surface to 1.05 m depth (with 5 or 10 cm spacing between sensors) in the other columns. The measurements were acquired every 15 minutes. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1971

Development of lithium diffused radiation resistant solar cells, part 2

The work performed to investigate the effect of various process parameters on the performance of lithium doped P/N solar cells is described. Effort was concentrated in four main areas: (1) the starting material, (2) the boron diffusion, (3) the lithium diffusion, and (4) the contact system. Investigation of starting material primarily involved comparison of crucible grown silicon (high oxygen content) and Lopex silicon (low oxygen content). In addition, the effect of varying growing parameters of crucible grown silicon on lithium cell output was also examined. The objective of the boron diffusion studies was to obtain a diffusion process which produced high efficiency cells with minimal silicon stressing and could be scaled up to process 100 or more cells per diffusion. Contact studies included investigating sintering of the TiAg contacts and evaluation of the contact integrity.

Payne, P. R.↗
NASA NTRS1971

Electron pitch-angle diffusion driven by oblique whistler-mode turbulence.

A general description of cyclotron harmonic resonant pitch-angle scattering is presented. Quasi-linear diffusion coefficients are prescribed in terms of the wave normal distribution of plasma wave energy. Numerical computations are performed for the specific case of relativistic electrons interacting with a band of low frequency whistler-mode turbulence. A parametric treatment of the wave energy distribution permits normalized diffusion coefficients to be presented graphically solely as a function of the electron pitch-angle. The diffusion coefficients generally decrease with increasing cyclotron harmonic number. Higher harmonic diffusion is insignificant at very small electron pitch-angles, but becomes increasingly important as the pitch-angle increases. One thus expected the rate of pitch-angle scattering to decrease with increasing electron energy, since the resonant value of the latter varies proportionately with harmonic number. This indicates that, in mirror-type magnet field geometrics, such as the earth's radiation belts, the diffusion losses of high energy electrons are likely to be appreciably slower than those at low energy.

Lyons, L. R.↗
NASA NTRS1972

Inner-zone energetic-electron repopulation by radial diffusion.

Quantitative study of the intrusion of natural electrons into the inner radiation zone during and after the geomagnetic storm of Sept. 2, 1966. It is shown that the transport is consistent with a radial-diffusion mechanism in which the first two invariants are conserved. Except for the three-day period of the storm main phase when data were missing, the radial-diffusion coefficient is found to have a value which could be produced by a variation of a large-scale electric field across the magnetosphere having an amplitude of 0.28 mV/m and a period of 1600 sec. Electric fields having approximately these characteristics have been inferred from previous observations of the motion of whistler ducts within the plasmapause. If fields of this amplitude and period exist throughout the magnetosphere, the radial diffusion of all geomagnetically trapped particles except the high-energy inner-zone protons is strongly influenced by electric-field variations. A comprehensive review of previously reported radial-diffusion coefficients shows reasonable agreement for L less than about 3.0, but serious discrepancies among reported values exist for determinations made in the outer zone. These discrepancies cannot be explained by the simple theory of radial diffusion due to variation of large-scale electric or magnetic fields.

Tomassian, A. D.↗
NASA NTRS1972

Resonant diffusion in the presence of strong plasma turbulence.

The diffusion equation which describes the evolution of the average one-particle distribution function for an ensemble of strongly turbulent plasmas is derived. The diffusion tensor is a time integral of the autocorrelation tensor of the fluctuations as observed by particles moving along statistically distributed orbits. These orbits contain the effects of fluctuations and thus differ from those encountered in weak turbulence theory. Two statistical orbit effects quadratic in the strength of the fluctuations affect the magnitude of the diffusion: (a) modification of the ensemble average orbits by the fluctuations, and (b) statistical dispersion in particle orbits about the average. The plasma trajectory equations are used to relate each to the diffusion tensor itself when the turbulence is electrostatic. The diffusion tensor is explicitly evaluated for a strongly turbulent unmagnetized plasma.

Birmingham, T. J.↗
NASA NTRS1973

Diffusion across the modified polyethylene separator GX in the heat-sterilizable AgO-Zn battery

Models of diffusion across an inert membrane have been studied using the computer program CINDA. The models were constructed to simulate various conditions obtained in the consideration of the diffusion of Ag (OH)2 ions in the AgO-Zn battery. The effects on concentrations across the membrane at the steady state and on the fluxout as a function of time were used to examine the consequences of stepwise reducing the number of sources of ions, of stepwise blocking the source and sink surfaces, of varying the magnitude of the diffusion coefficient for a uniform membrane, of varying the diffusion coefficient across the membrane, and of excluding volumes to diffusion.

Lutwack, R.↗
NASA NTRS1973

Effects of diffusion on aluminum depletion and degradation of NiAl coatings

Experiments were performed to critically demonstrate the effects of diffusion on the aluminum depletion and degradation of NiAl coatings on superalloys. Pack aluminized IN 100 and Mar-M200 were diffusion annealed in 0.0005 torr vacuum at 1100 C for 300 hours. Aluminum losses due to oxidation and vaporization were minimal. Metallographic and electron microprobe analyses showed considerable interdiffusion of the coating with the substrate, which caused a large decrease in the original aluminum level of the coating. Subsequent cyclic furnace oxidation tests were performed at 1100 C using 1 hour cycles on pre-diffused and as-coated specimens. The pre-diffusion treatment decreased the oxidation protection for both alloys, but more dramatically for IN 100. Identical oxidation tests of bulk NiAl, where such diffusion effects are precluded, showed no signs of degradation at twice the time needed to degrade the coated superalloys.

Smialek, J. L.↗
NASA NTRS1974

The diffusive idealization of charged particle transport in random magnetic fields

The transport of charged particles diffusing in a random magnetic field parallel to a relatively large guiding field is presented. The same coefficient of diffusion is obtained by three methods. Two corrections must be added to the expression in which the diffusive flux is proportional to the gradient of the density. Explicit expressions are given for a characteristic time and a characteristic length which describe the corrections. The well known divergence of the coefficient of diffusion, which is implied by the quasilinear analysis of pitch angle scattering, does not occur if the scattering rate is finite at 90 deg pitch angle. This effect is illustrated by formulas which give the coefficient of diffusion when the quasilinear expression is perturbed by a variable amount of isotropic scattering.

Earl, J. A.↗
NASA NTRS1974

The electron diffusion coefficient in Jupiter's magnetosphere

A steady-state model of Jupiter's electron radiation belt is developed. The model includes injection from the solar wind, radial diffusion, energy degradation by synchrotron radiation, and absorption at Jupiter's surface. A diffusion coefficient of the form D sub RR/R sub J squared = k times R to the m-th power is assumed, and then observed data on synchrotron radiation are used to fit the model. The free parameters determined from this fit are m = 1.95 plus or minus 0.5, k = 1.7 plus or minus 0.5 x 10 to the 9th power per sec, and the magnetic moment of injected particles equals 770 plus or minus 300 MeV/G. The value of m shows quite clearly that the diffusion is not caused by magnetic pumping by a variable solar wind or by a fluctuating convection electric field. The process might be field line exchange driven by atmospheric-ionospheric winds; our diffusion coefficient has roughly the same radial dependence but is considerably smaller in magnitude than the upper bound diffusion coefficients recently suggested for this process by Brice and McDonough (1973) and Jacques and Davis (1972).

Birmingham, T.↗