Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion coefficient”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Multiphysics and Multiscale Modeling of Coupled Transport of Chloride Ions in Concrete

Chloride ions (Cl−)-induced corrosion is one of the main degradation mechanisms in reinforced concrete (RC) structures. In most situations, the degradation initiates with the transport of Cl− from the surface of the concrete towards the reinforcing steel. The accumulation of Cl− at the steel-concrete interface could initiate reinforcement corrosion once a threshold Cl− concentration is achieved. An accurate numerical model of the Cl− transport in concrete is required to predict the corrosion initiation in RC structures. However, existing numerical models lack a representation of the heterogenous concrete microstructure resulting from the varying environmental conditions and the indirect effect of time dependent temperature and relative humidity (RH) on the water adsorption and Cl− binding isotherms. In this study, a numerical model is developed to study the coupled transport of Cl− with heat, RH and oxygen (O2) into the concrete. The modeling of the concrete microstructure is performed using the Virtual Cement and Concrete Testing Laboratory (VCCTL) code developed by the U.S. National Institute of Standards and Technology (NIST). The concept of equivalent maturation time is utilized to eliminate the limitation of simulating concrete microstructure using VCCTL in specific environmental conditions such as adiabatic. Thus, a time-dependent concrete microstructure, which depends on the hydration reactions coupled with the temperature and RH of the environment, is achieved to study the Cl− transport. Additionally, Cl− binding isotherms, which are a function of the pH of the concrete pore solution, Cl− concentration, and weight fraction of mono-sulfate aluminate (AFm) and calcium-silicate-hydrate (C-S-H), obtained from an experimental study by the same authors are utilized to account for the Cl− binding of cement hydration products. The temperature dependent RH diffusion was considered to account for the transport of Cl− with moisture transport. The temperature and RH diffusion in the concrete domain, composite theory, and Cl− binding and water adsorption isotherms are used in combination, to estimate the ensuing Cl− diffusion field within the concrete. The coupled transport process of heat, RH, Cl−, and O2 is implemented in the Multiphysics Object-Oriented Simulation Environment (MOOSE) developed by the U.S. Idaho National Laboratory (INL). The model was verified and validated using data from multiple experimental studies with different concrete mixture proportions, curing durations, and environmental conditions. Additionally, a sensitivity analysis was performed to identify that the water-to-cement (w/c) ratio, the exposure duration, the boundary conditions: temperature, RH, surface Cl− concentration, Cl− diffusion coefficient in the capillary water, and the critical RH are the important parameters that govern the Cl− transport in RC structures. In a case study, the capabilities of the developed numerical model are demonstrated by studying the complex 2D diffusion of Cl− in a RC beam located in two different climatic regions: warm and humid weather in Galveston, Texas, and cold and dry weather in North Minnesota, Minnesota, subjected to time varying temperature, RH, and surface Cl− concentrations.

composite theory↗
DOE OSTI2023

The stability of FeH x and hydrogen transport at Earth’s core mantle boundary

Iron hydride in Earth’s interior can be formed by the reaction between hydrous minerals (water) and iron. Studying iron hydride improves our understanding of hydrogen transportation in Earth’s interior. Our high-pressure experiments found that face-centered cubic (fcc) FeH x (x≤1) is stable up to 165 GPa, and our ab initio molecular dynamics simulations predicted that fcc FeH x transforms to a superionic state under lower mantle conditions. In the superionic state, H-ions in fcc FeH become highly diffusive-like fluids with a high diffusion coefficient of ~3.7 × 10 -4 cm 2 s -1 , which is comparable to that in the liquid Fe-H phase. The densities and melting temperatures of fcc FeH x were systematically calculated. Similar to superionic ice, the extra entropy of diffusive H-ions increases the melting temperature of fcc FeH. Further, the wide stability field of fcc FeH enables hydrogen transport into the outer core to create a potential hydrogen reservoir in Earth’s interior, leaving oxygen-rich patches (ORP) above the core mantle boundary (CMB).

58 GEOSCIENCES↗
DOE OSTI2024

Adsorption and Solvation Modulate Rhodamine B Diffusion in Ethanol/Water-Filled Anodic Alumina Nanopores

Confinement of solvents and solutes within nanoporous materials frequently leads to the emergence of unique mass transport behaviors that, once fully understood, may lead to improved chemical separations. Here, the diffusion of Rhodamine B (RhB) dye within 10 and 20 nm diameter anodic aluminum oxide (AAO) nanopores filled with binary ethanol/water mixtures is investigated. Mixture compositions spanning from pure ethanol to pure water are employed. The results of confocal fluorescence correlation spectroscopy studies reveal that RhB diffusion occurs by a two-component mechanism comprising composition-dependent fast and slow motions, characterized by diffusion coefficients D f and D s . The results are consistent with those of previous studies performed under more limited conditions [J. Phys. Chem. C, 2023, 127, 411-420]. The fast component scales with mixture viscosity and is assigned to hindered bulk-like diffusion in central pore regions. Slow diffusion likely involves adsorption of RhB to the pore surface and may be described by a desorption mediated mechanism. The occurrence of RhB adsorption to the AAO surface is verified at the single-molecule level by wide-field fluorescence imaging of membrane cross-sectional surfaces. Unique composition dependent trends in the autocorrelation amplitude and in D s that mimic bulk RhB solubility are revealed. D s is found to be smallest in pure ethanol and pure water and largest in intermediate mixtures. These results suggest that RhB surface adsorption is strongest in the pure liquids and weakest in mixtures of intermediate composition, where the dye is least soluble, and most soluble, respectively. As a result, molecular dynamics simulations reveal that a water layer appears on the pore surface under most conditions, while RhB is solvated primarily by ethanol. The composition dependence of RhB diffusion is concluded to reflect its solvation dependent interactions with the pore walls.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Machine Learning-Based Upscaling of Finite-Size Molecular Dynamics Diffusion Simulations for Binary Fluids

Molecular diffusion coefficients calculated using molecular dynamics (MD) simulations suffer from finite-size (i.e., finite box size and finite particle number) effects. Results from finite-sized MD simulations can be upscaled to infinite simulation size by applying a correction factor. For self-diffusion of single-component fluids, this correction has been well-studied by many researchers including Yeh and Hummer (YH); for binary fluid mixtures, a modified YH correction was recently proposed for correcting MD-predicted Maxwell–Stephan (MS) diffusion rates. In this study we use both empirical and machine learning methods to identify improvements to the finite-size correction factors for both self-diffusion and MS diffusion of binary Lennard-Jones (LJ) fluid mixtures. Using artificial neural networks (ANNs), the error in the corrected LJ fluid diffusion is reduced by an order of magnitude versus existing YH corrections, and the ANN models perform well for mixtures with large dissimilarities in size and interaction energies where the YH correction proves insufficient.

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

Fast and Universal Kohn-Sham Density Functional Theory Algorithm for Warm Dense Matter to Hot Dense Plasma

Understanding many processes, e.g., fusion experiments, planetary interiors, and dwarf stars, depends strongly on microscopic physics modeling of warm dense matter and hot dense plasma. This complex state of matter consists of a transient mixture of degenerate and nearly free electrons, molecules, and ions. This regime challenges both experiment and analytical modeling, necessitating predictive ab initio atomistic computation, typically based on quantum mechanical Kohn-Sham density functional theory (KS-DFT). However, cubic computational scaling with temperature and system size prohibits the use of DFT through much of the warm dense matter regime. A recently developed stochastic approach to KS-DFT can be used at high temperatures, with the exact same accuracy as the deterministic approach, but the stochastic error can converge slowly and it remains expensive for intermediate temperatures (< 50 eV). Here we have developed a universal mixed stochastic-deterministic algorithm for DFT at any temperature. This approach leverages the physics of KS-DFT to seamlessly integrate the best aspects of these different approaches. We demonstrate that this method significantly accelerated self-consistent field calculations for temperatures from 3 to 50 eV, while producing stable molecular dynamics and accurate diffusion coefficients.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Simulations of the churning mode: Toroidally symmetric plasma convection and turbulence around the X-points in a snowflake divertor

Using a reduced MHD model, extended to include field-aligned thermal conduction, we present numerical simulations of the churning mode (CM): a toroidally symmetric, non-linear plasma vortex in the vicinity of the null points in a snowflake (SF) divertor [D. D. Ryutov et al., Phys. Scr. 89, 088002 (2014)]. Simulations are carried out across a range of inter-null separations, d xx ⁠, and inter-null orientations, θ⁠, primarily in conditions relevant to the MAST-U tokamak. We find that, when d xx is small, the CM induces additional transport across the X-points when β pm ≳8%, where β pm is the ratio of the plasma pressure in the null region to poloidal magnetic pressure at the midplane. This transport also increases approximately linearly as d xx is reduced. A diffusive model of this transport is shown to predict the total transport across the null points, where diffusion coefficients of up to ∼ 10 2 m 2 s −1 centered on a small region around the X-points are used. However, the CM also results in significant changes to the flux surfaces in the null region, which is not captured by this diffusive model. The changes in magnetic geometry mean the fractional exhaust power delivered to each divertor leg is highly sensitive to β pm ⁠, d xx ⁠, and θ⁠. For small values of θ⁠, the CM can induce a change in topology, redirecting exhaust power from a secondary divertor leg on the high field side to one on the low field side. Similar behavior is found in the fraction of exhaust power going to the inner and outer divertor. Such changes in the flux surfaces may not be captured by Grad–Shafranov solvers and so may be a source of error in the magnetic reconstruction of SF experiments. We consistently find that the fractional exhaust power going to a secondary divertor leg on the high field side is small, consistent with SF experiments.

Computer simulation↗
DOE OSTI2026

Langevin and Fokker-Planck analyses for diffusion-mediated passing of circular and discorectangular species in two-dimensional channels

The propensity for pairs of diffusing species to pass each other within narrow channels or pores is of basic interest as a first-passage-type problem. It is also of relevance for solution-phase transport in nanoporous materials, and in particular for catalytic conversion reactions where high yield requires that product species can efficiently pass reactant species to exit the pores. Here, we analyze a two-dimensional model with nonoverlapping circular and discorectangular species confined to a rectangular channel, and where passing is mediated by Brownian dynamics in an implicit solvent. For narrower channels where passing is still possible, the discorectangle must align with the channel to pass the circular species. Behavior of the passing propensity, 𝑃, can be assessed by strongly damped Langevin simulations, or within an equivalent Fokker-Planck equation (FPE) formalism. The latter corresponds to a diffusion problem in a “higher-dimensional channel” with a constriction. We assess the variation of the passing propensity, 𝑃, for a broad range of channel width including its scaling just above the threshold where passing is sterically blocked. Analysis of 𝑃 versus the rotational diffusion coefficient 𝐷 𝑟 of the discorectangle reveals a significant decrease in 𝑃 for lower 𝐷 𝑟 for moderate channel width. This prompts a direct analysis of the regime where 𝐷 𝑟 → 0 , for which the FPE can be reduced to a three-dimensional diffusion problem, precise analysis of which is facilitated by adaptive-mesh finite element methods. The dependence of 𝑃 on the aspect ratio of the discorectangle is also assessed.

Rahman, Md Khaledur [Ames Laboratory (AMES), Ames,↗
DOE OSTI2022

Effects of Hydrogen Bonding on Nuclear Data Development of Liquid Anhydrous HF

Anhydrous Hydrogen Fluoride (HF) at high temperatures and pressures is used to process and manufacture nuclear fuel. As HF is often used directly with uranium, correct neutron thermal scattering cross sections are crucial to criticality safety applications. Classical molecular dynamics (CMD) simulation of the flexible HF system was used to create the thermal scattering law (TSL) and cross sections. The initial 2-site model is used in LAMMPS, and it can not capture the H-bond. To correctly represent the H-bond effects, a second, 3-site model was constructed in GROMACS. The 3-site model handled H-bonds by connecting a massless charge to the molecule. Key model parameters were compared to experimental data to verify the approach and models. To get the normalized VACF, the model was compared using hydrogen and fluorine bond length, density, potential energy, and diffusion coefficient. The phonon DOSs for both models were derived from the normalized VACF. DOSs were used to estimate the TSL ( S ( α, β )) and neutron thermal scattering cross sections for hydrogen in HF. The TSLs were evaluated using the FLASSH code with the Schofield diffusion model. It was observed that the representation of the hydrogen bonding changes the TSL's diffusional contributions. This is represented in the low energy scattering cross section, where intermolecular binding effects shift the cross section.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Effective Fragment Potential-Based Molecular Dynamics Studies of Diffusion in Acetone and Hexane

To facilitate more reliable descriptions of transport properties in liquids, molecular dynamics (MD) simulations are performed based on the effective fragment potential (EFP) method derived from first-principles quantum mechanics (in contrast to MD based upon empirically fitted potentials). The EFP method describes molecular interactions in terms of Coulomb, polarization/induction, dispersion, exchange-repulsion, and charge-transfer interactions. The EFP MD simulations described in this paper, performed on hexane and acetone, are able to track the mean-square displacement of molecules for sufficient time to reliably extract translational diffusion coefficients. Finally, the results reported here are in reasonable agreement with experiment.

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

Simultaneously Enhancing Exciton/Charge Transport in Organic Solar Cells by an Organoboron Additive

Efficient exciton diffusion and charge transport play a vital role in advancing the power conversion efficiency (PCE) of organic solar cells (OSCs). Here, a facile strategy is presented to simultaneously enhance exciton/charge transport of the widely studied PM6:Y6-based OSCs by employing highly emissive trans-bis(dimesitylboron)stilbene (BBS) as a solid additive. BBS transforms the emissive sites from a more H-type aggregate into a more J-type aggregate, which benefits the resonance energy transfer for PM6 exciton diffusion and energy transfer from PM6 to Y6. Transient gated photoluminescence spectroscopy measurements indicate that addition of BBS improves the exciton diffusion coefficient of PM6 and the dissociation of PM6 excitons in the PM6:Y6:BBS film. Transient absorption spectroscopy measurements confirm faster charge generation in PM6:Y6:BBS. Moreover, BBS helps improve Y6 crystallization, and current-sensing atomic force microscopy characterization reveals an improved charge-carrier diffusion length in PM6:Y6:BBS. Finally, owing to the enhanced exciton diffusion, exciton dissociation, charge generation, and charge transport, as well as reduced charge recombination and energy loss, a higher PCE of 17.6% with simultaneously improved open-circuit voltage, short-circuit current density, and fill factor is achieved for the PM6:Y6:BBS devices compared to the devices without BBS (16.2%).

14 SOLAR ENERGY↗
DOE OSTI2023

Halide Substitution Effects on Lithium-Ion Diffusion in Protonated Antiperovskites

Solid-state electrolytes (SSEs) for all-solid-state lithium-ion batteries are generating intense interest because these batteries can improve the safety and performance compared with devices fabricated with conventional, flammable liquid electrolytes. In these SSEs, it has been suggested that Li + ion diffusion in the grain boundaries is hindered and is a critical determinant of the overall ionic conductivity (σ). However, Li + ion diffusivities in the grain (D G ) and the grain boundary (D GB ) are difficult to determine experimentally, with few techniques capable of distinguishing the individual contributions. Here, we distinguished the D G and D GB for the protonated lithium antiperovskites (pLiAPs) SSEs: Li 2 OHCl, Li 2 OHBr, Li 2 OHF 0.1 Cl 0.9 , Li 2 OHF 0.1 Br 0.9 , and Li 2 OHCl 0.3 7Br 0.63 . The measurements were obtained directly from 7 Li pulsed-field gradient nuclear magnetic resonance (PFG-NMR) at 353 K. The 7 Li PFG-NMR echo profiles were composed of two primary components with additional secondary oscillatory components – the so-called NMR diffraction phenomenon. The length scale separating the two main components corresponds to a diffusion length of ~1.7 µm, which is thought to be the average grain size (by diameter). The short-range (≤ 1.7 µm) diffusion component associated with D G (≈10 -11 m 2 /s) varied minimally with halide substitution, while the long-range (≥ 1.7 µm) component D GB (≈10 -12 to 10 -15 m 2 /s) was highly sensitive to the substitution of halides and closely correlated with s. In addition, from the comparison of the ratio D GB /D G to D t (the Li + ion diffusion coefficient estimated from the rotational correlation time, t c ), it was determined that the contribution of D G to σ is negligible; 0.01 ~ 0.04 in the pLiAPs studied here. Finally, these insights provide fundamental understanding of the halide substitution effects on Li + ion grain versus grain boundary diffusion, and suggest that careful engineering of the grain boundaries at the microscopic level is necessary to achieve high-performance pLiAP SSEs.

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

Investigation of desorption of hydrogen gas from polymer matrix using thermal desorption analysis and finite element modeling

Polymer materials are used widely as sealing materials and liners in hydrogen storage and transportation infrastructure. During operating conditions, polymers are exposed to high pressure hydrogen gas periodically. When in contact, the hydrogen gas diffuses through the polymer material with ease due to its smaller sized molecules compared to the polymer molecules. In the event of rapid decompression of outside hydrogen gas, the hydrogen molecules within the polymer diffuse out. This process is known as the desorption of polymer and can cause blisters or cracking causing permanent damage to the polymer. Here, we used a combined experimental and modeling approach to investigate the desorption of the hydrogen gas in the polymer. Thermal desorption analysis (TDA) experiments were performed on the polymer sample exposed to the high-pressure hydrogen gas to measure the diffusion coefficient and equilibrium hydrogen concentration. Experiments also provided the quantity of hydrogen within the polymer during the desorption process which followed Fick’s law. However, for higher pressurization values, some anomaly was observed in the quantity of hydrogen remaining in the polymer sample. A continuum mechanics based coupled diffusion–deformation-damage model was deployed to simulate the gas desorption process in the material using the finite element method to understand its effect on the material. It was found that the anomaly in TDA results were due to the blister formation on the polymer surface which resulted in longer hydrogen retention. Simulations also showed that blister formation can be affected by multiple parameters such as the addition of filler particles, and pre-existing cavity size.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Use of Hydrogel Electrolyte in Zn-MnO2 Rechargeable Batteries: Characterization of Safety, Performance, and Cu2+ Ion Diffusion

Achieving commercially acceptable Zn-MnO2 rechargeable batteries depends on the reversibility of active zinc and manganese materials, and avoiding side reactions during the second electron reaction of MnO2. Typically, liquid electrolytes such as potassium hydroxide (KOH) are used for Zn-MnO2 rechargeable batteries. However, it is known that using liquid electrolytes causes the formation of electrochemically inactive materials, such as precipitation Mn3O4 or ZnMn2O4 resulting from the uncontrollable reaction of Mn3+ dissolved species with zincate ions. In this paper, hydrogel electrolytes are tested for MnO2 electrodes undergoing two-electron cycling. Improved cell safety is achieved because the hydrogel electrolyte is non-spillable, according to standards from the US Department of Transportation (DOT). The cycling of “half cells” with advanced-formulation MnO2 cathodes paired with commercial NiOOH electrodes is tested with hydrogel and a normal electrolyte, to detect changes to the zincate crossover and reaction from anode to cathode. These half cells achieved ≥700 cycles with 99% coulombic efficiency and 63% energy efficiency at C/3 rates based on the second electron capacity of MnO2. Other cycling tests with “full cells” of Zn anodes with the same MnO2 cathodes achieved ~300 cycles until reaching 50% capacity fade, a comparable performance to cells using liquid electrolyte. Electrodes dissected after cycling showed that the liquid electrolyte allowed Cu ions to migrate more than the hydrogel electrolyte. However, measurements of the Cu diffusion coefficient showed no difference between liquid and gel electrolytes; thus, it was hypothesized that the gel electrolytes reduced the occurrence of Cu short circuits by either (a) reducing electrode physical contact to the separator or (b) reducing electro-convective electrolyte transport that may be as important as diffusive transport.

25 ENERGY STORAGE↗
DOE OSTI2020

Hydrogen Diffusivity Measurements of YH 1.87 Moderator Material with Incoherent Quasielastic Neutron Scattering

Yttrium hydride is an optimal choice of moderator material for thermal neutron spectrum reactors requiring small core volumes. The Transformation Challenge Reactor (TCR) program at Oak Ridge National Laboratory (ORNL) aims to develop an additively manufactured microreactor core by leveraging recent advances in materials, manufacturing, data analytics, and high-fidelity modeling and simulation. Yttrium hydride was selected as the moderator material for TCR due to its high moderating power, which enables the spatially efficient reactor design and thermal stability that is superior when compared with zirconium hydride. Hydrogen desorption from yttrium hydride is expected at elevated temperatures, posing safety and performance concerns. As a lifetime component in advanced reactors, appropriate strategies to mitigate hydrogen release from yttrium hydride over long-term reactor operation are required. Understanding hydrogen mobility in yttrium hydride provides the scientific basis for developing a hydrogen barrier for hydride moderator and predicting the hydrogen redistribution within the material during various operational conditions. This study investigated the hydrogen motion in YH 1.87 as a function of temperature ranging from 200 to 1,173 K using incoherent quasielastic neutron scattering (IQNS) at the Spallation Neutron Source (SNS). The results are presented and discussed in terms of hydrogen self-diffusion coefficients, activation energy for translational hydrogen motion, and hydrogen jump distances. Based on the IQNS data, YH 1.87 is found to be stable at high temperatures with no translational H motion below 1,023 K within the IQNS detection limits. In the temperature range of 1,073–1,173 K, hydrogen diffusivity in YH 1.87 is $D = 4.57 X 10^{⁻3} [m^2/s] exp (⁻\frac{1.73[eV]}{RT}$). The hydrogen concentration’s dependence on hydrogen diffusivity is also discussed in this report.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Developing the Science Basis for Understanding Polymer Encapsulant Degradation Mechanisms: DuraMAT 2.0 Final Project Report

Polymeric encapsulants are essential materials in photovoltaic modules, protecting sensitive electronics from the environment while providing mechanical integrity to the multilayered assembly. However, these polymeric materials are susceptible to degradation processes driven by the ingress of environmental species, ultraviolet radiation, thermal stresses, and mechanical loading. In this study, we employ a combined atomistic simulation and accelerated aging experimental approach to study the molecular-scale mechanisms of encapsulant degradation. Classical molecular dynamics simulations quantify the diffusion of environmental and degradation species through the polymer matrix, producing composition-specific diffusion coefficients. Reactive simulations characterize activation energy barriers and reaction rate constants for key chemical pathways. In parallel, thermal-desorption analyses coupled with mass spectrometry monitor the emergence and concentration profiles of degradation products under controlled stressor conditions. By integrating simulation and experiment, we establish quantitative correlations between polymer composition, species diffusivity, and chemical reactivity. We anticipate that these relations and quantitative values could serve as high-fidelity inputs to reaction-diffusion models, enabling physics-informed lifetime predictions and guiding the design of more durable encapsulant materials for solar energy applications.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Assessment the leaching characteristics and long-term leaching behavior of some radionuclides from synthesized zeolite cement matrix

Highlights: • Synthesis and characterization of zeolite cement synthesized from fly ash • Leaching tests were used to assess the release rate of {sup 85}Sr and {sup 60}Co radionuclides. • Simple mathematical models were used to detect leaching mechanisms. • Long-term leaching behavior of {sup 85}Sr and {sup 60}Co from cementitious matrix was assessed. The development of cementitious materials remains a vital goal to produce valuable products with good mechanical, physical and chemical properties suitable for the safe disposal of concentrated radionuclides resulted from the treatment of contaminated solutions. In this paper, zeolite cement with good properties was laboratory synthesized from industrial by-products fly ash and characterized using various characterization techniques. The leaching characteristics of strontium and cobalt radionuclides from the synthesized zeolite cement were scrutinized conforming with the standard leaching methodology of International Atomic Energy Authority (IAEA). A mechanical strength assessment was executed to characterize the extent of immobilization process of the solidified matrices. The cumulative leaching fraction (φ) of the two studied radionuclides was found to be less than 5% in all examined conditions, which implies the applicability of the IAEA-recommended methodology for estimating the diffusion coefficient. The experimental leaching data were regressed nonlinearly to various mathematical kinetic models to assess the controlling leaching mechanism and to determine the leaching parameters. The regression results indicated that strontium ({sup 85}Sr) and cobalt ({sup 60}Co) leaching resulted from two succeeding mechanisms: as first order kinetic reaction, and then diffusion. The calculated values of leachability indices signify that the performance of each of the studied matrices is within an acceptable range. A simplified mathematical model, rooted in the first order reaction and diffusion mechanisms, was simulated to predict the radionuclides leaching rates from zeolite cement matrices. By comparing the synthesized zeolite cement with other cementitious materials, it can be concluded that the synthesized material can be classified as an efficient material suitable to immobilize {sup 85}Sr and {sup 60}Co from radioactive wastes. The acquired findings demonstrated that the studied immobilized waste matrices have acceptable mechanical effectiveness.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Data-driven methods for diffusivity prediction in nuclear fuels

The growth rate of structural defects in nuclear fuels under irradiation is intrinsically related to the diffusion rates of the defects in the fuel lattice. The generation and growth of atomistic structural defects can significantly alter the performance characteristics of the fuel. This alteration of functionality must be accurately captured to qualify a nuclear fuel for use in reactors. Predicting the diffusion coefficients of defects and how they impact macroscale properties such as swelling, gas release, and creep is therefore of significant importance in both the design of new nuclear fuels and the assessment of current fuel types. In this article, we apply data-driven methods focusing on machine learning (ML) to determine various diffusion properties of two nuclear fuels—uranium oxide and uranium nitride. We show that using ML can increase, often significantly, the accuracy of predicting diffusivity in nuclear fuels in comparison to current analytical models. We also illustrate how ML can be used to quickly develop fuel models with parameter dependencies that are more complex and robust than what is currently available in the literature. In conclusion, these results suggest there is potential for ML to accelerate the design, qualification, and implementation of nuclear fuels.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Oxygen diffusion in garnet: Experimental calibration and implications for timescales of metamorphic processes and retention of primary O isotopic signatures

Knowledge of oxygen diffusion in garnet is crucial for a correct interpretation of oxygen isotope signatures in natural samples. A series of experiments was undertaken to determine the diffusivity of oxygen in garnet, which remains poorly constrained. The first suite included high-pressure (HP), nominally dry experiments performed in piston-cylinder apparatus at: (1) T = 1050–1600 °C and P = 1.5 GPa and (2) T = 1500 °C and P = 2.5 GPa using yttrium aluminum garnet (YAG; Y3Al5O12) cubes. Second, HP H2O-saturated experiments were conducted at T = 900 °C and P = 1.0–1.5 GPa, wherein YAG crystals were packed into a YAG + Corundum powder, along with 18O-enriched H2O. Third, 1 atm experiments with YAG cubes were performed in a gas-mixing furnace at T = 1500–1600 °C under Ar flux. Finally, an experiment at T = 900 °C and P = 1.0 GPa was done using a pyrope cube embedded into pyrope powder and 18O-enriched H2O. Experiments using grossular were not successful. Profiles of 18O/(18O+16O) in the experimental charges were analyzed with three different secondary ion mass spectrometers (SIMS): sensitive high-resolution ion microprobe (SHRIMP II and SI), CAMECA IMS-1280, and NanoSIMS. Considering only the measured length of 18O diffusion profiles, similar results were obtained for YAG and pyrope annealed at 900 °C, suggesting limited effects of chemical composition on oxygen diffusivity. However, in both garnet types, several profiles deviate from the error function geometry, suggesting that the behavior of O in garnet cannot be fully described as simple concentration-independent diffusion, certainly in YAG and likely in natural pyrope as well. The experimental results are better described by invoking O diffusion via two distinct pathways with an inter-site reaction allowing O to move between these pathways. Modeling this process yields two diffusion coefficients (D values) for O, one of which is approximately two orders of magnitude higher than the other. Taken together, Arrhenius relationships are:logDm2s-1=-7.2(±1.3)+(-321(±32)kJmol-12.303RT) for the slow pathway, andlogDm2s-1=-5.4(±0.7)+(-321(±20)kJmol-12.303RT) for the fast pathway. We interpret the two pathways as representing diffusion following vacancy and inter-stitial mechanisms, respectively. Regardless, our new data suggest that the slow mechanism is prevalent in garnet with natural compositions, and thus is likely to control the retentivity of oxygen isotopic signatures in natural samples. The diffusivity of oxygen is similar to Fe-Mn diffusivity in garnet at 1000–1100 °C and Ca diffusivity at 850 °C. However, the activation energy for O diffusion is larger, leading to lower diffusivities at P-T conditions characterizing crustal metamorphism. Therefore, original O isotopic signatures can be retained in garnets showing major element zoning partially re-equilibrated by diffusion, with the uncertainty caveat of extrapolating the experimental data to lower temperature conditions.

Geochemistry & Geophysics↗