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 739 records · Page 41

Relationship between Surface Self-Diffusion and Bubble Mobility in Solids: Theory and Atomistic Simulation

Diffusion processes in solids are characterized by complex mechanisms occurring at the atomistic level. The relevant theoretical concepts are sufficiently developed in the case of the diffusion of point defects. At the same time, the diffusion of objects such as clusters of vacancies or nanometer-size cavities in the crystal lattice is described only in terms of continuous-medium approximation. This paper attempts to link the existing theory of diffusion of bubbles in a crystalline matrix with the atomistic description of the process on the basis of the molecular dynamics method. For simplicity, the diffusion of cavities in a two-dimensional Lennard-Jones lattice is considered. A controlled molecular dynamics method is proposed to speed up the direct computation of the diffusion rate of bubbles. Based on the results of simulation, a practical interpretation is given to a parameter of the continuum theory such as the rate of surface self-diffusion. The applicability range of the continuum theory is demonstrated, and principles for the refinement of the theory in the case of minimum-size bubbles are formulated.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Atomic-level mechanisms of short-circuit diffusion in materials

This paper reviews the recent progress in understanding the atomic mechanisms of short-circuit diffusion along materials interfaces, such as grain and interphase boundaries, as well as lattice and interfacial dislocations/disconnections. Recent atomistic computer simulations have shown that short-circuit diffusion is dominated by collective atomic rearrangements in the form of strings and rings of mobile atoms. The process is dynamically heterogeneous in space and time and has many features in common with atomic dynamics in supercooled glass-forming liquids. We discuss examples of grain boundary, interphase boundary, and dislocation diffusion in metals and alloys, including the solute effect on the diffusion rates and mechanisms. Interphase boundaries are exemplified by Al–Si interfaces with diverse orientation relationships and atomic structures. The hierarchy of short-circuit diffusion paths in materials is reviewed by comparing the rates of grain boundary, interphase boundary, and dislocation diffusion. Future directions in the field of short-circuit diffusion in defect core regions are discussed.

36 MATERIALS SCIENCE↗

Diffusion and Leaching Behavior of Radionuclides in Category 3 Waste Encasement Concrete and Soil Fill Material – Summary Report

One of the methods being considered for safely disposing of Category 3 low-level radioactive wastes is to encase the waste in concrete. Such concrete encasement would contain and isolate the waste packages from the hydrologic environment and would act as an intrusion barrier. The current plan for waste isolation consists of stacking low-level waste packages on a trench floor, surrounding the stacks with reinforced steel, and encasing these packages in concrete. These concrete-encased waste stacks are expected to vary in size with maximum dimensions of 6.4 m long, 2.7 m wide, and 4 m high. The waste stacks are expected to have a surrounding minimum thickness of 15 cm of concrete encasement. These concrete-encased waste packages are expected to withstand environmental exposure (solar radiation, temperature variations, and precipitation) until an interim soil cover or permanent closure cover is installed, and to remain largely intact thereafter. Any failure of concrete encasement may result in water intrusion and consequent mobilization of radionuclides from the waste packages. The mobilized radionuclides may escape from the encased concrete by mass flow and/or diffusion and move into the surrounding subsurface environment. Therefore, it is necessary to assess the performance of the concrete encasement structure and the ability of the surrounding soil to retard radionuclide migration. The retardation factors for radionuclides contained in the waste packages can be determined from measurements of diffusion coefficients for these contaminants through concrete and fill material. Some of the mobilization scenarios include (1) potential leaching of waste form before permanent closure cover is installed; (2) after the cover installation, long-term diffusion of radionuclides from concrete waste form into surrounding fill material; (3) diffusion of radionuclides from contaminated soils into adjoining concrete encasement and clean fill material. Additionally, the rate of diffusion of radionuclides may be affected by the formation of structural cracks in concrete, the carbonation of the buried waste form, and any potential effect of metallic iron (in the form of rebars) on the mobility of radionuclides. The radionuclides iodine-129 ({sup 129}I), technetium-99 ({sup 99}Tc), and uranium-238 ({sup 238}U) are identified as long-term dose contributors in Category 3 waste (Mann et al. 2001; Wood et al. 1995). Because of their anionic nature in aqueous solutions, {sup 129}I, {sup 99}Tc, and carbonate-complexed {sup 238}U may readily leach into the subsurface environment (Serne et al. 1989, 1992a, b, 1993, and 1995). The leachability and/or diffusion of radionuclide species must be measured to assess the long-term performance of waste grouts when contacted with vadose-zone pore water or groundwater. Although significant research has been conducted on the design and performance of cementitious waste forms, the current protocol conducted to assess radionuclide stability within these waste forms has been limited to the Toxicity Characteristic Leaching Procedure, Method 1311 Federal Registry (EPA 1992) and ANSI/ANS-16.1 leach test (ANSI 1986). These tests evaluate the performance under water-saturated conditions and do not evaluate the performance of cementitious waste forms within the context of waste repositories which are located within water-deficient vadose zones. Moreover, these tests assess only the diffusion of radionuclides from concrete waste forms and neglect evaluating the mechanisms of retention, stability of the waste form, and formation of secondary phases during weathering, which may serve as long-term secondary hosts for immobilization of radionuclides. The results of recent investigations conducted under arid and semi-arid conditions (Al-Khayat et al. 2002; Garrabrants et al. 2002; Garrabrants and Kosson 2003; Garrabrants et al. 2004; Gervais et al. 2004; Sanchez et al. 2002; Sanchez et al. 2003) provide valuable information suggesting structural and chemical changes to concrete waste forms which may affect contaminant containment and waste form performance. However, continued research is necessitated by the need to understand: the mechanism of contaminant release; the significance of contaminant release pathways; how waste form performance is affected by the full range of environmental conditions within the disposal facility; the process of waste form aging under conditions that are representative of processes occurring in response to changing environmental conditions within the disposal facility; the effect of waste form aging on chemical, physical, and radiological properties, and the associated impact on contaminant release. Recent reviews conducted by the National Academies of Science recognized the efficacy of cementitious materials for waste isolation, but further noted the significant shortcomings in our current understanding and testing protocol for evaluating the performance of various formulations.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Deep Learning Approaches to Surrogates for Solving the Diffusion Equation for Mechanistic Real-World Simulations

In many mechanistic medical, biological, physical, and engineered spatiotemporal dynamic models the numerical solution of partial differential equations (PDEs), especially for diffusion, fluid flow and mechanical relaxation, can make simulations impractically slow. Biological models of tissues and organs often require the simultaneous calculation of the spatial variation of concentration of dozens of diffusing chemical species. One clinical example where rapid calculation of a diffusing field is of use is the estimation of oxygen gradients in the retina, based on imaging of the retinal vasculature, to guide surgical interventions in diabetic retinopathy. Furthermore, the ability to predict blood perfusion and oxygenation may one day guide clinical interventions in diverse settings, i.e., from stent placement in treating heart disease to BOLD fMRI interpretation in evaluating cognitive function (Xie et al., 2019; Lee et al., 2020). Since the quasi-steady-state solutions required for fast-diffusing chemical species like oxygen are particularly computationally costly, we consider the use of a neural network to provide an approximate solution to the steady-state diffusion equation. Machine learning surrogates, neural networks trained to provide approximate solutions to such complicated numerical problems, can often provide speed-ups of several orders of magnitude compared to direct calculation. Surrogates of PDEs could enable use of larger and more detailed models than are possible with direct calculation and can make including such simulations in real-time or near-real time workflows practical. Creating a surrogate requires running the direct calculation tens of thousands of times to generate training data and then training the neural network, both of which are computationally expensive. Often the practical applications of such models require thousands to millions of replica simulations, for example for parameter identification and uncertainty quantification, each of which gains speed from surrogate use and rapidly recovers the up-front costs of surrogate generation. We use a Convolutional Neural Network to approximate the stationary solution to the diffusion equation in the case of two equal-diameter, circular, constant-value sources located at random positions in a two-dimensional square domain with absorbing boundary conditions. Such a configuration caricatures the chemical concentration field of a fast-diffusing species like oxygen in a tissue with two parallel blood vessels in a cross section perpendicular to the two blood vessels. To improve convergence during training, we apply a training approach that uses roll-back to reject stochastic changes to the network that increase the loss function. The trained neural network approximation is about 1000 times faster than the direct calculation for individual replicas. Because different applications will have different criteria for acceptable approximation accuracy, we discuss a variety of loss functions and accuracy estimators that can help select the best network for a particular application. We briefly discuss some of the issues we encountered with overfitting, mismapping of the field values and the geometrical conditions that lead to large absolute and relative errors in the approximate solution.

60 APPLIED LIFE SCIENCES↗

The Effect of Acetylation on Iron Uptake and Diffusion in Water Saturated Wood Cell Walls and Implications for Decay

Acetylation is widely used as a wood modification process that protects wood from fungal decay. The mechanisms by which acetylation protects wood are not fully understood. With these experiments, we expand upon the literature and test whether previously observed differences in iron uptake by wood were a result of decreased iron binding capacity or slower diffusion. We measured the concentration of iron in 2 mm thick wood sections at 0, 10, and 20% acetylation as a function of time after exposure to iron solutions. The iron was introduced either strongly chelated with oxalate or weakly chelated with acetate. The concentrations of iron and oxalate in solution were chosen to be similar to those found during brown rot decay, while the concentration of iron and acetate matched previous work. The iron content of oxalate-exposed wood increased only slightly and was complete within an hour, suggesting little absorption and fast diffusion, or only slight surface adsorption. The increase in iron concentration from acetate solutions with time was consistent with Fickian diffusion, with a diffusion coefficient on the order of 10 -16 m 2 s -1 . The rather slow diffusion rate was likely due to significant binding of iron within the wood cell wall. The diffusion coefficient did not depend on the acetylation level; however, the capacity for iron absorption from acetate solution was greatly reduced in the acetylated wood, likely due to the loss of OH groups. We explored several hypotheses that might explain why the diffusion rate appears to be independent of the acetylation level and found none of them convincing. Implications for brown rot decay mechanisms and future research are discussed.

54 ENVIRONMENTAL SCIENCES↗

Theory of plane, symmetrical intake diffusers

The present report ties in with the investigations on the inlet diffusers by P. Ruden. The theory developed by Ruden had produced results which found excellent confirmation in wind-tunnel tests and in spite of certain still-existing defects, are technically very promising. The reasons for the new theory of the diffuser forms indicated by Ruden are twofold: first, the arguments adduced in Ruden's theory deal only with one specific operating condition, that is, a certain ratio of mean velocity within the diffuser to flying speed, while in the present report any desired velocity ratios are involved; second, a different choice of parameters and the increased possibilities of variation result in diffuser forms which cannot be reconciled at once with Ruden's theory. The first enables a theoretical check of the measurements made with Ruden's diffusers at variable velocity ratio, the second permits the calculation of diffuser types which in many respects are superior to Ruden's diffusers.

FLOW - CONFORMAL TRANSFORMATION - HODOGRAPH METHOD↗

Tortuosity measurement and the effects of finite pulse widths on xenon gas diffusion NMR studies of porous media

We have extended the utility of NMR as a technique to probe porous media structure over length scales of approximately 100-2000 microm by using the spin 1/2 noble gas 129Xe imbibed into the system's pore space. Such length scales are much greater than can be probed with NMR diffusion studies of water-saturated porous media. We utilized Pulsed Gradient Spin Echo NMR measurements of the time-dependent diffusion coefficient, D(t), of the xenon gas filling the pore space to study further the measurements of both the pore surface-area-to-volume ratio, S/V(p), and the tortuosity (pore connectivity) of the medium. In uniform-size glass bead packs, we observed D(t) decreasing with increasing t, reaching an observed asymptote of approximately 0.62-0.65D(0), that could be measured over diffusion distances extending over multiple bead diameters. Measurements of D(t)/D(0) at differing gas pressures showed this tortuosity limit was not affected by changing the characteristic diffusion length of the spins during the diffusion encoding gradient pulse. This was not the case at the short time limit, where D(t)/D(0) was noticeably affected by the gas pressure in the sample. Increasing the gas pressure, and hence reducing D(0) and the diffusion during the gradient pulse served to reduce the previously observed deviation of D(t)/D(0) from the S/V(p) relation. The Pade approximation is used to interpolate between the long and short time limits in D(t). While the short time D(t) points lay above the interpolation line in the case of small beads, due to diffusion during the gradient pulse on the order of the pore size, it was also noted that the experimental D(t) data fell below the Pade line in the case of large beads, most likely due to finite size effects.

Non-NASA Center↗

Fission Product Diffusion in Silicon Carbide: As Revealed by Computer Simulations and Experimental Measurements

The diffusion and release of fission products through silicon carbide in the tri-structural isotropic (TRISO) fuel particles remain unsolved for decades. The underlying mechanism is quite challenging to be determined. To help unveil the mysterious story, the current work applies molecular dynamics method to show the stability of silver, palladium, ruthenium and iodine as an interstitial and their atomic diffusion along coincident site lattice (CSL) boundary, especially ?3 grain boundary (GB). The major finding presents a much faster diffusion along GB than in bulk for all elements considered. The reasonably close estimate to experiments and simulations where available has confirmed the important role of grain boundary diffusion of Ag and Pd in SiC. However, the discrepancy addressed in Ag with measurements from fuel studies suggest a more complicated mechanism, which might be in correlation with high energy grain boundaries or the presence of crack. The subsequent characterization of Ru and I distribution in SiC-PyC-SiC diffusion couples, which have been ion irradiated at 900?C to 10 dpa and 20 dpa, has performed by secondary ion mass spectrometry (SIMS) analysis. The experiment measurements correlate well with the grain boundary diffusion by simulation, which provide further evidences that the grain boundary diffusion cannot be neglected once the fission products are accessible at the grain boundary.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Bounce- and MLT-averaged diffusion coefficients in a physics-based magnetic field geometry obtained from RAM-SCB for the March 17 2013 storm

Local acceleration via whistler wave and particle interaction plays a significant role in particle dynamics in the radiation belt. In this work we explore gyro-resonant wave-particle interaction and quasi-linear diffusion in different magnetic field configurations related to the March 17 2013 storm. We consider the Earth's magnetic dipole field as a reference and compare the results against non-dipole field configurations corresponding to quiet and stormy conditions. The latter are obtained with the ring current-atmosphere interactions model with a self-consistent magnetic field RAM-SCB, a code that models the Earth's ring current and provides a realistic modeling of the Earth's magnetic field. By applying quasi-linear theory, the bounce- and MLT-averaged electron pitch angle, mixed term, and energy diffusion coefficients are calculated for each magnetic field configuration. For radiation belt (~1 MeV) and ring current (~100 keV) electrons, it is shown that at some MLTs the bounce-averaged diffusion coefficients become rather insensitive to the details of the magnetic field configuration, while at other MLTs storm conditions can expand the range of equatorial pitch angles where gyro-resonant diffusion occurs and significantly enhance the diffusion rates. When MLT average is performed at drift shell L = 4.25 (a good approximation to drift average), the diffusion coefficients become quite independent of the magnetic field configuration for relativistic electrons, while the opposite is true for lower energy electrons. These results suggest that, at least for the March 17 2013 storm and for L ≲ 4.25, the commonly adopted dipole approximation of the Earth's magnetic field can be safely used for radiation belt electrons, while a realistic modeling of the magnetic field configuration is necessary to describe adequately the diffusion rates of ring current electrons.

79 ASTRONOMY AND ASTROPHYSICS↗

Au and Pt Diffusion in Electrodeposited Amorphous Sb 2 Te 3 Thin Films

Amorphous Sb 2 Te 3 thin films are synthesized using electrodeposition in different thicknesses to explore their solid electrolytic nature with respect to some noble metals. Au and Pt are used as the substrate materials and thermally diffused into the Sb 2 Te 3 films under passivated and nonpassivated conditions. Rutherford backscattering spectrometry is used to study the Au and Pt diffusion into the films as a function of depth. It is found that Au diffuses into the Sb 2 Te 3 thin films even at room temperature. In contrast, there is no observable diffusion of Pt at room temperature. At a higher temperature (i.e., 200 °C), both Au and Pt diffuse into the Sb 2 Te 3 thin films and Au diffusion is more significant than Pt. Overall, the findings suggest that amorphous Sb 2 Te 3 can be used as a good solid electrolyte to permit the diffusion of some noble metals such as Au and Pt.

36 MATERIALS SCIENCE↗

Thermal diffusivity recovery and defect annealing kinetics of self-ion implanted tungsten probed by insitu transient grating spectroscopy

Tungsten is a promising candidate material for plasma-facing armour components in future fusion reactors. A key concern is irradiation-induced degradation of its normally excellent thermal transport properties. In this comprehensive study, thermal diffusivity degradation in ion-implanted tungsten and its evolution from room temperature (RT) to 1073 K is considered. Five samples were exposed to 20 MeV self-ions at RT to achieve damage levels ranging from 3.2 x 10 -4 to 3.2 displacements per atom (dpa). Transient grating spectroscopy with in-situ heating was then used to study thermal diffusivity evolution as a function of temperature. Using a kinetic theory model, an equivalent point defect density is estimated from the measured thermal diffusivity. The results showed a prominent recovery of thermal diffusivity between 450 K and 650 K, which coincides with the onset of mono-vacancy mobility. After 1073 K annealing samples with initial damage of 3.2 x 10 -3 dpa or less recover close to the pristine value of thermal diffusivity. For doses of 3.2 x 10 -2 dpa or higher, on the other hand, a residual reduction in thermal diffusivity remains even after 1073 K annealing. Transmission electron microscopy reveals that this is associated with extended, irradiation-induced dislocation structures that are retained after annealing. A sensitivity analysis shows that thermal diffusivity provides an efficient tool for the assessing total defect content in tungsten up to 1000 K.

36 MATERIALS SCIENCE↗

Application of PDF mixing models to premixed flames with differential diffusion

Differential diffusion alters the balance of reaction and diffusion in turbulent premixed combustion, affecting the performance and emissions of combustion devices. Modelling combustion devices with Probability or Filtered Density Function (PDF or FDF) methods provides an exact treatment for the change in composition due to chemical reaction, while molecular mixing has to be modelled. Previous PDF molecular mixing models do not account for differential diffusion in a manner which satisfies realizability requirements. A new approach for treating differential diffusion, which ensures realizability, is proposed for pairwise-exchange mixing models in general, and applied in the Interaction by Exchange with the Mean (IEM) model of Dopazo [26], and in the Euclidean Minimum Spanning Tree (EMST) model of Subramaniam and Pope [5]. The new differential diffusion models are referred to as IEM-DD and EMST-DD respectively. Results from two and three-dimensional DNS of turbulent premixed methane–air combustion show that mixing rates and conditional statistics of species mass fractions depend on species diffusivities and the combustion regime. Zero-dimensional PDF model results obtained for the two-dimensional DNS case show that the EMST-DD model best reproduces the features that characterize differential diffusion in the DNS. The essential feature of the EMST-DD model, which accounts for its success in turbulent premixed combustion, is that differential mixing rates are imposed within a model which mixes locally in composition space.

Energy & Fuels↗

IoT-based retrofit information diffusion in future smart communities

Community-scale building retrofits are not merely scaled-up versions of single-building retrofits. They involve complex challenges, such as reconciling individual interests with collective goals and managing the dynamic interplay between buildings through mechanisms like power grids and social connections. Internet of Things (IoT) connectivity holds the potential to leverage these interplays to balance individual and collective interests effectively in smart communities. One critical aspect of this interplay is information diffusion, which shapes how retrofit decisions spread among neighbors, influencing individual choices and ultimately impacting community-level retrofit outcomes. In other words, IoT-based smart devices automatically push tailored retrofit notifications to homeowners, which completely changes the format of information diffusion in the future. To investigate this influence by such information diffusion, the study used CityBES to simulate energy performance for different retrofits and applied an information diffusion model to analyze how decisions spread in a networked community of 192 buildings. The diffusion process was modeled on a weighted, directed network, capturing the dynamics of information flow and decision-making across 16 scenarios. Individual retrofit benefits were evaluated through payback years, while community-level retrofit outcomes were assessed using greenhouse gas (GHG) emission reductions. The results demonstrate that easier information diffusion among neighbors encourages households to prioritize retrofit measures that align with the majority’s optimal choices, even at the expense of individual financial benefits. In this case, such collective prioritization enhanced community-level retrofit performance, increasing GHG emission reductions by up to 29.4 %. However, this improvement came with trade-offs, as the average payback period for households extended by approximately 1.74 years. These findings highlight the potential of IoT-based information diffusion in future smart communities to coordinate individual interests with collective goals, ultimately accelerating community-level building retrofits.

Shu, Lei↗

Hydrogen diffusion induced dislocation transformations in a nickel superalloy

The diffusion of hydrogen in metals and alloys induces embrittlement that can adversely affect the structural properties. We examine the adsorption and diffusion of hydrogen in Inconel-718 (IN-718), and scrutinize the ensuing effects on the dislocation behavior in the alloy to elucidate the fundamental mechanisms of hydrogen-microstructure interactions from classical molecular simulations. Hydrogen adsorption increases with time until the surface saturates, while hydrogen diffusion exhibits strong temperature dependence, with diffusion coefficients converging above 1300 K regardless of the initial hydrogen concentration in the alloy. The diffusion in IN-718 is significantly sluggish than in pure Ni, Fe, or Cr, and is strongly impacted by hydrogen concentrations, resulting in an order of magnitude higher diffusion coefficient for hydrogen (10-14 m2/s relative to 10-15 m2/s) at high concentrations, especially below 600 K. Hydrogen diffusion coefficient varies from 10-12 to 10-15 m2/s in IN-718 depending on temperature (500–1400 K). More critically, our results reveal that increasing hydrogen concentration induces microstructural changes in the alloy, transforming perfect dislocations into stair-rods and Shockley partials, with higher temperatures favoring the latter. The results are significant for hydrogen fuel applications to gain insights into the materials chemistry for designing safer and more efficient propulsion systems, particularly in high-performance environments related to controlled hydrogen combustion applications.

Banerjee, Tanumoy↗

Diffusion of Zr, Ru, Ce, Y, La, Sr and Ba fission products in UO 2

The diffusivity of the solid fission products (FP) Zr (Zr 4+ ), Ru (Ru 4+ , Ru 3+ ), Ce (Ce 4+ ), Y (Y 3+ ), La (La 3+ ), Sr (Sr 2+ ) and Ba (Ba 2+ ) by a vacancy mechanism has been calculated, using a combination of density functional theory (DFT) and empirical potential (EP) calculations. The activation energies for the solid fission products are compared to the activation energy for Xe fission gas atoms calculated previously. Apart from Ru, the solid fission products all exhibit higher activation energy than Xe. Furthermore, for all solid FPs except Y 3+ , the migration of the FP has lower barrier than the migration of a neighboring U atom, making the latter the rate limiting step for direct migration. An indirect mechanism, consisting of two successive migrations around the FP, is also investigated. The calculated diffusivities show that most solid fission products diffuse with rates similar to U self-diffusion. But, Ru, Ba and Sr exhibit faster diffusion than the other solid FPs, with Ru 3+ and Ru 4+ diffusing even faster than Xe for T < 1200 K. The diffusivities correlate with the observed fission product solubility in UO 2 , and the tendency to form metallic and oxide second phase inclusions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Simulation of radiation driven fission gas diffusion in UO 2 , ThO 2 and PuO 2

Below 1000 K it is thought that fission gas diffusion in nuclear fuel during irradiation occurs through atomic mixing due to radiation damage. Here we present a molecular dynamics (MD) study of Xe, Kr, Th, U, Pu and O diffusion due to irradiation. It is concluded that the ballistic phase does not sufficiently account for the experimentally observed diffusion. Thermal spike simulations are used to confirm that electronic stopping remedies the discrepancy with experiment and the predicted diffusivities lie within the scatter of the experimental data. Here, our results predict that the diffusion coefficients are ordered such that D* 0 > D* Kr > D* Xe > D* U . For all species >98.5% of diffusivity is accounted for by electronic stopping. Fission gas diffusivity was not predicted to vary significantly between ThO 2 , UO 2 and PuO 2 , indicating that this process would not change greatly for mixed oxide fuels.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Atomistic modeling of out-of-pile xenon diffusion by vacancy clusters in UO 2

We report that diffusion of Xe atoms in UO 2 fuel is important for nuclear fuel performance, and is enabled by interaction with U and O vacancies. Previous work using atomistic calculations based on density functional theory (DFT) and empirical potentials (EP) focused on the role of small vacancy clusters (Xe U2Oy , y = 0, 1, 2) for Xe transport, but the model predictions failed to reproduce the experimental observations by consistently underestimating Xe diffusivity. In this work, where we focused on out-of-pile conditions (i.e. in the absence of irradiation), we have explored two of the uncertainties associated with the model: the DFT methodology and the types of clusters considered. Additionally, we found that using the energy barriers obtained with GGA + U DFT allows to reconcile simulation results and experimental observations in the intrinsic regime (i.e. out-of-pile conditions). The Xe UO cluster is the preferred configuration at high temperature, while the mobile Xe U2O cluster takes over below 1720 K. The latter cluster is also the main contributor to Xe diffusion and good agreement is obtained with current models empirically fitted to experiments for diffusion at high temperature. A simple expression that captures most of Xe diffusivity is proposed, relying on Xe UO and Xe U2O clusters only. We also determined the out-of-pile concentration and diffusivity of extended clusters Xe U z O y , 3 ≤ x ≤ 9, and determined that these do not contribute significantly to the overall diffusivity under intrinsic conditions. However, we speculate that the Xe U 4 O 3 and Xe U 8 O 9 clusters, which have relatively low migration barrier, may play a role under irradiation conditions, because the radiation-induced enhanced vacancy concentration would rapidly increase the fraction of large Xe-vacancy clusters.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Diffusion behavior of lanthanide-additive compounds (Ce 4 Sb 3 , Ce 2 Sb, and CeTe) against HT9 and Fe

Antimony and tellurium have been identified as promising additives in metallic fuel, which can immobilize free-lanthanide fission products into stable intermetallic compounds in order to mitigate the fuel-cladding chemical interaction. Ce 4 Sb 3 , Ce 2 Sb, and CeTe are the primary compounds formed by Sb or Te with the lanthanide Ce present in the fuel. If these compounds are present at the outer periphery of the fuel, they will come in contact with and react with the cladding after the fuel swells. The present study investigates the reactivity of these compounds with two cladding materials, HT9 and Fe. The diffusion couple tests between these compounds and HT9 or Fe were conducted at 853 K. Scanning electron microscopy and transmission electron microscopy were used to analyze the morphology, microstructure, and phase distribution of the diffusion region. It was observed that the diffusion region thickness formed by the three compounds was significantly reduced compared to free Ce. Moreover, there was no observed diffusion or reaction between Ce 4 Sb 3 or Ce 2 Sb with either HT9 or Fe. CeTe was found to diffuse and react with HT9, forming Cr 3 Te 4 and TeFe at the diffusion region, as well as to penetrate into Fe, mainly by intergranular diffusion.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗