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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Modified Closures in Monte Carlo Algorithms for Diffusive Binary Stochastic Media Transport Problems

In a stochastic medium, the material properties at a given spatial location are known only statistically. The most common approach to solving particle transport problems involving binary stochastic media (BSM) is to use the atomic mix approximation in which the transport problem is solved using homogenized (volume-averaged) material properties. A common deterministic model developed for solving BSM particle transport problems is the Levermore-Pomraning (LP) model. Zimmerman and Adams proposed a Monte Carlo algorithm (Algorithm A) that solves the LP equations and another Monte Carlo algorithm (Algorithm B) that locally preserves the sampled material realization; we refer to these Monte Carlo algorithms as the LP and LPLRP (LP local realization preserving) algorithms, respectively. One-dimensional (1D) planar geometry benchmark studies have shown that the LPLRP algorithm is often significantly more accurate than the LP algorithm for problems with an incident angular flux as well as for problems with an interior source. The LPLRP algorithm implementation in these benchmark comparisons made explicit use of the one-dimensional nature of the problem. The LP model is derived by assuming an 'upwind' closure in the coupling term relating the two materials in a binary stochastic medium. Su and Pomraning developed a modified form of this closure by considering the small correlation length limit and requiring the modified closure to produce the correct exponential decay for a source-free half-line albedo problem in rod geometry. They concluded that the modified closure is generally not inferior to the LP closure and in some cases is significantly better. Brantley investigated the use of the Su-Pomraning (SP) closure in the Monte Carlo LP algorithm (LP-SP) for the suite of benchmark problems and concluded that 1) the LP-SP algorithm was somewhat more accurate overall than the LP algorithm and somewhat less accurate overall than the LPLRP algorithm for an incident angular flux benchmark suite, and 2) the LP-SP algorithm was generally the least accurate of the algorithms for an interior source benchmark suite. Larsen et al. performed an asymptotic analysis of the transport equation in 1D planar geometry for the situation in which the physical system is 1) a random binary stochastic medium with material macroscopic total cross sections and mean slab width values of O(1) or smaller and 2) optically thick with weak absorption and sources at each spatial point and therefore globally diffusive. The asymptotic analysis demonstrates that, under these assumptions, the transport equation limits to the conventional diffusion equation with atomically-mixed (volume-averaged) material properties. Larsen et al. also demonstrate that, under these same assumptions, the LP equations asymptotically limit to an atomically-mixed diffusion equation with a diffusion coefficient that is too large, leading to an unphysical flattening of the ensemble-averaged scalar flux distribution. Vasques and Yadav recently performed an asymptotic analysis of an adjusted Levermore-Pomraning closure in which the Markovian transition functions are rescaled such that the equations asymptotically limit to the correct atomically-mixed diffusion equation; we refer to this model as LP-VY. In this paper, we demonstrate that the adjusted LP closure proposed by Vasques and Yadav is a special case of the SP closure obtained by assuming that the absorption cross sections in the two materials are equal. We further demonstrate that the Su-Pomraning closure has the correct asymptotic limit for the diffusive physical system and is therefore an appropriate closure for these problems. We also describe how to incorporate these various closures into the Monte Carlo LP and LPLRP particle transport algorithms, and we present numerical results comparing the accuracy of the Monte Carlo LP, LP-VY, LP-SP, LPLRP, and LPLRP-SP algorithms for the set of diffusive benchmark problems. We demonstrated that the SP closure has the correct asymptotic limit for the diffusive physical system under investigation in this paper and that the VY closure is a special case of the SP closure obtained by assuming equal absorption in the two materials. Both of these closures are readily implemented in the Monte Carlo LP and LPLRP algorithms for particle transport in BSM. Through numerical comparisons to a diffusive benchmark suite, we confirmed that 1) the LP and LPLRP algorithms do not limit to the appropriate diffusion solution and 2) the LP-VY, LP-SP, and LPLRP-SP algorithms limit to the appropriate diffusion solution and are of comparable accuracy. The LPLRP algorithm has been shown to be generally more accurate than the LP algorithm for non-diffusive problems. However, the numerical results in this paper indicate that the LPLRP algorithm does not limit to the correct diffusion equation for the asymptotic limit under consideration. The work in this paper demonstrates that the SP algorithm does limit to the correct diffusion equation for the asymptotic limit considered. Future work will investigate the possibility of whether the LPLRP-SP algorithm is accurate for both diffusive and non-diffusive problems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Computational determination of a primary diffusion mode in γ U-10Mo under irradiation

Low enriched uranium (< 20 % 235 U)-molybdenum (U-Mo) monolithic fuel is the primary candidate for high performance research and test reactors, and is in the process of being qualified to replace highly enriched uranium (≥ 20 % 23 5U) fuel. As part of the qualification process, it is critical to understand and predict the behavior of fission gas bubbles under irradiation, which affects fuel swelling and fuel failure. Mechanistic fuel models are being developed that can both reproduce the existing experimental data for fuel swelling, and be further applied to irradiation conditions beyond the experimental scope. Diffusion of species under irradiation conditions is an important parameter in the mechanistic fuel models; however, no temperature-relevant experimental diffusion data exists. In the present work, radiation-enhanced diffusion coefficients of U, Mo, and Xe in λU-10wt.%Mo were calculated in the temperature range between 300 K and 1400 K via rate-theory models and molecular dynamics simulations with an embedded-atom method interatomic potential for the U-Mo-Xe system. Accordingly, total diffusion coefficients under relevant irradiation conditions are determined using previously obtained intrinsic thermal diffusion and radiation-driven diffusion coefficients, as well as the newly calculated radiation-enhanced diffusion coefficients presented herein. Radiation-enhanced diffusion of U and Mo was dominant in the intermediate temperature range, whereas radiation-enhanced diffusion of Xe did not significantly contribute to total diffusion of Xe at the relevant fission rate densities. Radiation-enhanced diffusion of Xe became faster than both intrinsic thermal diffusion and radiation-driven diffusion at a fission rate density of 5 x 10 22 fissions/m 3 /s, which is higher than the typical fission rate density range in research reactors. The temperature regime that radiation-enhanced diffusion of each element dominated was dependent on the fission rate density. Finally, the total diffusion coefficients of U, Mo, and Xe, updated in this work, will be utilized as parameters in the mechanistic fuel models to help predict the behavior of fission gas bubbles under irradiation more accurately.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Investigation of Molecular Diffusion at Block Copolymer Thin Films Using Maximum Entropy Method-Based Fluorescence Correlation Spectroscopy and Single Molecule Tracking

Fluorescence correlation spectroscopy (FCS) has been widely used to investigate molecular diffusion behavior in various samples. The use of the maximum entropy method (MEM) for FCS data analysis provides a unique means to determine multiple distinct diffusion coefficients without a priori assumption of their number. Comparison of the MEM-based FCS method (MEM-FCS) with another method will reveal its utility and advantage as an analytical tool to investigate diffusion dynamics. Herein, we measured diffusion of fluorescent probes doped into nanostructured thin films using MEM-FCS, and validated the results with single molecule tracking (SMT) data. The efficacy of the MEM code employed was first demonstrated by analyzing simulated FCS data for systems incorporating one and two diffusion modes with broadly distributed diffusion coefficients. The MEM analysis accurately afforded the number of distinct diffusion modes and their mean diffusion coefficients. These results contrasted with those obtained by fitting the simulated data to conventional two-component and anomalous diffusion models, which yielded inaccurate estimates of the diffusion coefficients. Subsequently, the MEM analysis was applied to FCS data acquired from hydrophilic dye molecules incorporated into microphase-separated polystyrene-block-poly(ethylene oxide) (PS-b-PEO) thin films characterized under a water-saturated nitrogen atmosphere. The MEM analysis revealed distinct fast and slow diffusion components attributable to molecules diffusing on the film surface and inside the film, respectively. SMT studies of the same materials yielded trajectories for mobile molecules that appear to follow the curved PEO microdomains. Diffusion coefficients obtained from the SMT data were consistent with those obtained for the slow diffusion component detected by MEM-FCS. Furthermore, these results highlight the utility of MEM-FCS and SMT for gaining complementary information on molecular diffusion processes in heterogeneous material systems.

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

A First-Principles Density Function Theory Study of Tritium Diffusion in Li 2 ZrO 3 : Application for Producing Tritium

The ceramic material Li 2 ZrO 3 has superior thermo-physical and thermo-chemical properties and is highly compatible with other blanket materials used in nuclear reactors. Like LiAlO 2 , it could be used in the form of an annular pellet in tritium-producing burnable absorber rods (TPBARs) to produce tritium( 3 H) upon thermal neutron irradiation of lithium isotope ( 6 Li). The radiation damaged pellet crystal contains vacancies, defects of its constituent elements, and several other trapping sites which hinder the 3 H diffusion process and releasing behavior. In this study, we investigate the diffusion mechanisms of 3 H and O 3 H species in Li 2 ZrO 3 ceramic pellet in order to understand the effects on diffusion barriers and diffusion coefficients due to the presence of interstitial and substitutional Li defects, hydroxide (O– 3 H) vacancy defect, and of the interactions of 3 H with O-vacancies in the radiation damaged Li 2 ZrO 3 crystal. We consider several possible diffusion pathways of interstitial and substitutional 3 H and its species in a defective supercell and calculate the activation energy barrier profiles. We find that the smallest activation energy barrier is 0.3 eV and corresponding diffusion coefficient at 600 K is 1.93 × 10 –9 m 2 /s for 3 H substitutional diffusion. The smallest 3 H interstitial diffusion barrier and diffusion coefficient are found to be 0.34 eV and 3.25 × 10 –9 m 2 /s. By examining several channels for diffusion, our results show the most likely diffusion mechanism of 3 H migration is occurring along the c-direction by exchange of 3 H within and between different Li sites. Here, our calculated results reveal that the smallest energy barrier for O 3 H diffusion is 0.75 eV which is when O 3 H is diffusing to the O–Li vacancy pair and the corresponding diffusion coefficient is 6.31 × 10 –13 m 2 /s. In terms of 3 H diffusion, the obtained results indicate that the performance of Li 2 ZrO 3 could be better than γ-LiAlO 2 , a widely used ceramic material in TPBAR for producing 3 H.

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

Predicting the influence of long-range molecular interactions on macroscopic-scale diffusion by homogenization of the Smoluchowski equation

The macroscopic diffusion constant for a charged diffuser is in part dependent on (1) the volume excluded by solute “obstacles” and (2) long-range interactions between those obstacles and the diffuser. Increasing excluded volume reduces transport of the diffuser, while long-range interactions can either increase or decrease diffusivity, depending on the nature of the potential. We previously demonstrated [P. M. Kekenes-Huskey et al., Biophys. J. 105, 2130 (2013)] using homogenization theory that the configuration of molecular-scale obstacles can both hinder diffusion and induce diffusional anisotropy for small ions. As the density of molecular obstacles increases, van der Waals (vdW) and electrostatic interactions between obstacle and a diffuser become significant and can strongly influence the latter's diffusivity, which was neglected in our original model. Here, we extend this methodology to include a fixed (time-independent) potential of mean force, through homogenization of the Smoluchowski equation. We consider the diffusion of ions in crowded, hydrophilic environments at physiological ionic strengths and find that electrostatic and vdW interactions can enhance or depress effective diffusion rates for attractive or repulsive forces, respectively. Additionally, we show that the observed diffusion rate may be reduced independent of non-specific electrostatic and vdW interactions by treating obstacles that exhibit specific binding interactions as “buffers” that absorb free diffusers. Finally, we demonstrate that effective diffusion rates are sensitive to distribution of surface charge on a globular protein, Troponin C, suggesting that the use of molecular structures with atomistic-scale resolution can account for electrostatic influences on substrate transport. This approach offers new insight into the influence of molecular-scale, long-range interactions on transport of charged species, particularly for diffusion-influenced signaling events occurring in crowded cellular environments.

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

Charge Carrier Diffusion Length Determination in 4H-SiC Schottky Alpha Detectors

When developing semiconductor radiation detectors detailed characterization of the detector behavior is required. The ability of a semiconductor radiation detector to collect the charge liberated by incident ionizing radiation is a critical characteristic. Several factors influence the charge collection efficiency of a semiconductor radiation detector, including electron-hole trapping, applied electric field strength, and dopant concentrations and dimensions of low doped regions. Applying a reverse bias electric potential to a low doped semiconductor region creates a free charge carrier depleted region and establishes an electric field which can remove subsequently freed charge in the form of an electrical signal. When ionizing radiation is incident upon a semiconductor electron-hole pairs are generated. These charge carriers are created along the path of the ionizing radiation regardless of the doping concentration. However, only the charge freed in the low doped epitaxial region of a detector will be collected as a result of the applied electric field. When the electric field is insufficient to develop a depleted region equal to or larger than the range of the ionizing radiation in the detector material some charge freed by the ionizing radiation outside of the sensitive region will be lost, thus potentially resulting in less than complete charge collection. The incomplete collection of charge complicates the analysis of the incident particle energy spectrum. Charge freed in the highly doped bulk region beyond the low doped epitaxial region may diffuse back into the low doped depleted region and be subsequently included in the integral pulse. The amount of charge collected as a result of this diffusion mechanism is dependent upon the diffusion length of the charge carriers in the material. The diffusion length is subsequently dependent on the material properties of the detector including the number of defects, trapping centers, and temperature. In order to fully characterize the properties of 4H-SiC as a radiation detector the diffusion length was calculated by adjusting the theoretical diffusion length value used in the calculation of the charge collection efficiency to agree with empirically observed charge collection efficiencies. The discrepancy between the observed and theoretical charge collection efficiency at 50 to 150 volts applied bias indicates an unexplained charge collection phenomenon which is contributing to the larger than predicted charge collection efficiency of the detectors. This voltage range corresponds to an intermediate applied bias, indicating that either the theoretical depletion width is too short, the diffusion length varies with the composition of the region that the charge carriers are diffusing through, or some other phenomenon is involved. If the charge carrier diffusion length varies depending on whether or not the diffusing charge is freed in unbiased epitaxy or bulk layer, then a function describing the variation of the diffusion layer could be determined by measuring the diffusion length at the applied voltages used previously, but over a range of epitaxial layer thicknesses. Future work includes the characterization of the effects of high temperature operation on the diffusion length of 4H-SiC in addition to measuring the effects of varying energy alpha emitter electrodeposition on the observed diffusion length. Additional Schottky barrier contacts, namely those with a nickel-platinum contact, will be tested in various sizes which have been developed in order to eliminate certain failure mechanisms and explore the potential effects of Schottky contact size on charge collection efficiency. (authors)

36 MATERIALS SCIENCE↗
DOE OSTI2017

Comparison of Experimental Methods for Estimating Matrix Diffusion Coefficients for Contaminant Transport Modeling

Diffusion cell and diffusion wafer experiments were conducted to compare methods for estimating matrix diffusion coefficients in rock core samples from Pahute Mesa at the Nevada Nuclear Security Site (NNSS). A diffusion wafer method, in which a solute diffuses out of a rock matrix that is pre-saturated with water containing the solute, is presented as a simpler alternative to the traditional through-diffusion (diffusion cell) method. Both methods yielded estimates of matrix diffusion coefficients that were within the range of values previously reported for NNSS volcanic rocks. The difference between the estimates of the two methods ranged from 14 to 30%, and there was no systematic high or low bias of one method relative to the other. From a transport modeling perspective, these differences are relatively minor when one considers that other variables (e.g., fracture apertures, fracture spacings) influence matrix diffusion to a greater degree and tend to have greater uncertainty than diffusion coefficients. For the same relative random errors in concentration measurements, the diffusion cell method yields diffusion coefficient estimates that have less uncertainty than the wafer method. However, the wafer method is easier and less costly to implement and yields estimates more quickly, thus allowing a greater number of samples to be analyzed for the same cost and time. Given the relatively good agreement between the methods, and the lack of any apparent bias between the methods, the diffusion wafer method appears to offer advantages over the diffusion cell method if better statistical representation of a given set of rock samples is desired.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2018

Comparison of experimental methods for estimating matrix diffusion coefficients for contaminant transport modeling

Here, diffusion cell and diffusion wafer experiments were conducted to compare methods for estimating effective matrix diffusion coefficients in rock core samples from Pahute Mesa at the Nevada Nuclear Security Site (NNSS). A diffusion wafer method, in which a solute diffuses out of a rock matrix that is pre-saturated with water containing the solute, is presented as a simpler alternative to the traditional through-diffusion (diffusion cell) method. Both methods yielded estimates of effective matrix diffusion coefficients that were within the range of values previously reported for NNSS volcanic rocks. The difference between the estimates of the two methods ranged from 14 to 30%, and there was no systematic high or low bias of one method relative to the other. From a transport modeling perspective, these differences are relatively minor when one considers that other variables (e.g., fracture apertures, fracture spacings) influence matrix diffusion to a greater degree and tend to have greater uncertainty than effective matrix diffusion coefficients. For the same relative random errors in concentration measurements, the diffusion cell method yields effective matrix diffusion coefficient estimates that have less uncertainty than the wafer method. However, the wafer method is easier and less costly to implement and yields estimates more quickly, thus allowing a greater number of samples to be analyzed for the same cost and time. Given the relatively good agreement between the methods, and the lack of any apparent bias between the methods, the diffusion wafer method appears to offer advantages over the diffusion cell method if better statistical representation of a given set of rock samples is desired.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Density functional theory calculations of self- and Xe diffusion in U 3 Si 2

Uranium silicide, U 3 Si 2 , has been proposed as an advanced nuclear fuel to be used in light water reactors (LWRs). Development of this alternative to the predominant current fuel, UO 2 , is motivated by enhanced accident tolerance as a result of higher thermal conductivity as well as improved fuel cycle economics through increased uranium density. In order to accurately model the fuel performance of U 3 Si 2 , the diffusion rate of point defects, which is related to self-diffusion, and of fission gas atoms must be determined. DFT calculations are used to predict the U and Si point defect concentrations, the corresponding self-diffusivities, the preferred Xe trap site and the Xe diffusivity. Effects of irradiation are not considered. A low defect formation energy and a high entropy for Si interstitials give rise to Si-rich non-stoichiometry at elevated temperatures. Both U and Si self-diffusion and Xe diffusion are anisotropic as a consequence of the tetragonal crystal structure of U 3 Si 2 . Si diffusion occurs by interstitial mechanisms in both the a-b plane and along the c axis, while the U c axis diffusion rate is controlled by a vacancy mechanism. Interstitial diffusion of U is very fast in the a-b plane of the U 3 Si 2 crystal structure. Xe atoms prefer to occupy U vacancy trap sites. The highest Xe diffusion rate occurs by a vacancy mechanism in both the a-b plane and along the c axis. The diffusion rate is similar in the a-b plane and along the c axis. U and Si self-diffusion and Xe diffusion are all faster in U3Si2 than intrinsic U and Xe diffusion in conventional UO 2 nuclear fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Seismic Diffusivity and the Influence of Heterogeneity on Injection‐Induced Seismicity

Abstract The spatiotemporal patterns of injection‐induced seismicity (IIS) are commonly interpreted with the concept of a triggering front, which propagates in a diffusion‐like manner with an associated diffusivity parameter. Here, we refer to this diffusivity as the “seismic diffusivity.” Several previous studies implicitly assume that seismic diffusivity is equivalent to the effective hydraulic diffusivity of the subsurface, which describes the behavior of the mean pressure field in heterogeneous porous media. Seismicity‐based approaches for hydraulic characterization or simulations of IIS using domains of homogeneous equivalent porous media are implicitly based on this assumed equivalence. However, seismicity is expected to propagate with the threshold triggering pressure, and thus not be controlled by the evolution of the mean pressure field. We present numerical simulations of fluid injection to compare the seismic and effective hydraulic diffusivities in heterogeneous formations (including fractured rock). The numerical model combines uncoupled, linear pressure diffusion with the Mohr‐Coulomb failure criterion to simulate IIS. We demonstrate that connected pathways of relatively high hydraulic diffusivity in heterogeneous media (particularly in fractured rock domains) allow the threshold triggering pressure to propagate more rapidly than predicted by the effective hydraulic diffusivity. As a result, the seismic diffusivity is greater than the effective hydraulic diffusivity in heterogeneous porous media, possibly by an order of magnitude or more. Additionally, we present a case study of IIS near Soultz‐sous‐Forêts where seismic diffusivity is found to be at least one order of magnitude larger than the effective hydraulic diffusivity.

Haagenson, Ryan↗
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↗
DOE OSTI2022

Effect of charge and anisotropic diffusivity on ion exchange kinetics in nuclear waste form materials

Motifs of radionuclide absorption particles (absorbers or fillers) for the nuclear waste treatment usually have crystal structures with nano-sized tunnels for fast ion diffusion. On the other hand, the accumulation of ions generates an electric field which affects the ion diffusion. Therefore, the strong anisotropic kinetic properties and the electric field affect ion exchange kinetics, especially in an assembly of these particles which may block the fast diffusion tunnels for each other. In this work, we developed a mesoscale model to describe the effect of anisotropic kinetic properties and electric field on ion exchange kinetics. With the model we simulated the effect of anisotropic diffusivity, electrochemical potential, particle morphology, size and aggregation on ion exchange kinetics during batch experiments. Here, tt is found that 1) ion exchange results in charge accumulation inside the particle because different ions have different diffusivities and different diffusion driving forces. The electric neutrality is not preserved; 2) the electric field speeds up the slower diffusion ions (either uptake or release) and slows down the faster diffusion ions to reduce charge accumulation and maintains charge neutrality; 3) the ion exchange becomes very slow at the later stages because a) diffusion driving force (electrochemical gradient, Cs + concentration in the solution and Na + inside the particle) continuously decreases, as observed in batch experiments; and b) the diffusivity decreases due to the Cs + concentration dependence of ion diffusivity. This can explain the observed experimental phenomena that the radionuclide capacity usually cannot be reached; 4) for needle-like particle and their clusters, both surface area and diffusion distance along the direction with the fastest diffusivities have significant impact on the ion exchange kinetics.

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

Perturbative thermal diffusivity from partial sawtooth crashes in Alcator C-Mod

Perturbative thermal diffusivity has been measured on Alcator C-Mod via the use of the extended-time-to-peak method on heat pulses generated by partial sawtooth crashes. Perturbative thermal diffusivity governs the propagation of heat pulses through a plasma. It differs from power balance thermal diffusivity, which governs steady state thermal transport. Heat pulses generated by sawtooth crashes have been used extensively in the past to study heat pulse thermal diffusivity, but the details of the sawtooth event typically lead to non-diffusive 'ballistic' transport, making them an unreliable measure of perturbative diffusivity on many tokamaks. Partial sawteeth are common on numerous tokamaks, and generate a heat pulse without the 'ballistic' transport that often accompanies full sawteeth. This is the first application of the extended-time-to-peak method of diffusivity calculation to partial sawtooth crashes. Furthermore, this analysis was applied to over 50 C-Mod shots containing both L- and I-Mode. Results indicate correlations between perturbative diffusivity and confinement regime (L- versus I-mode), as well as correlations with local temperature, density, the associated gradients, and gradient scale lengths (a/L Te and a/L n ). In addition, diffusivities calculated from partial sawteeth are compared to perturbative diffusivities calculated with the nonlinear gyrokinetic code GYRO. We find that standard ion-scale simulations (ITG/TEM turbulence) under-predict the perturbative thermal diffusivity, but new multi-scale (ITG/TEM coupled with ETG) simulations can match the experimental perturbative diffusivity within error bars for an Alcator C-Mod L-mode plasma. Perturbative diffusivities extracted from heat pulses due to partial sawteeth provide a new constraint that can be used to validate gyrokinetic simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Diffusion of aromatic hydrocarbons in hierarchical mesoporous H-ZSM-5 zeolite

Hierarchical mesoporous zeolites exhibit higher catalytic activities and longer lifetime compared to the traditional microporous zeolites due to improved diffusivity of substrate molecules and their enhanced access to the zeolite active sites. Understanding diffusion of biomass pyrolysis vapors and their upgraded products in such materials is fundamentally important during catalytic fast pyrolysis (CFP) of lignocellulosic biomass, since diffusion makes major contribution to determine shape selectivity and product distribution. However, diffusivities of biomass relevant species in hierarchical mesoporous zeolites are poorly characterized, primarily due to the limitations of the available experimental technology. In this work, molecular dynamics (MD) simulations are utilized to investigate the diffusivities of several selected coke precursor molecules, benzene, naphthalene, and anthracene, in hierarchical mesoporous H-ZSM-5 zeolite. The effects of temperature and size of mesopores on the diffusivity of the chosen model compounds are examined. The simulation results demonstrate that diffusion within the microspores as well as on the external surface of mesoporous H-ZSM-5 dominates only at low temperature. At pyrolysis relevant temperatures, mass transfer is essentially conducted via diffusion along the mesopores. Furthermore, the results illustrate the heuristic diffusion model, such as the extensively used Knudsen diffusion, overestimates the diffusion of bulky molecules in the mesopores, thus making MD simulation a powerful and compulsory approach to explore diffusion in zeolites.

09 BIOMASS FUELS↗
DOE OSTI2014

Molecular dynamics simulation of framework flexibility effects on noble gas diffusion in HKUST-1 and ZIF-8

Molecular dynamics simulations were used to investigate trends in noble gas (Ar, Kr, Xe) diffusion in the metal-organic frameworks HKUST-1 and ZIF-8. Diffusion occurs primarily through inter-cage jump events, with much greater diffusion of guest atoms in HKUST-1 compared to ZIF-8 due to the larger cage and window sizes in the former. We compare diffusion coefficients calculated for both rigid and flexible frameworks. For rigid framework simulations, in which the framework atoms were held at their crystallographic or geometry optimized coordinates, sometimes dramatic differences in guest diffusion were seen depending on the initial framework structure or the choice of framework force field parameters. When framework flexibility effects were included, argon and krypton diffusion increased significantly compared to rigid-framework simulations using general force field parameters. Additionally, for argon and krypton in ZIF-8, guest diffusion increased with loading, demonstrating that guest-guest interactions between cages enhance inter-cage diffusion. No inter-cage jump events were seen for xenon atoms in ZIF-8 regardless of force field or initial structure, and the loading dependence of xenon diffusion in HKUST-1 is different for rigid and flexible frameworks. Diffusion of krypton and xenon in HKUST-1 depends on two competing effects: the steric effect that decreases diffusion as loading increases, and the “small cage effect” that increases diffusion as loading increases. Finally, a detailed analysis of the window size in ZIF-8 reveals that the window increases beyond its normal size to permit passage of a (nominally) larger krypton atom.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Computing eddy-driven effective diffusivity using Lagrangian particles

A novel method to derive effective diffusivity from Lagrangian particle trajectory data sets is developed and then analyzed relative to particle-derived meridional diffusivity for eddy-driven mixing in an idealized circumpolar current. Quantitative standard dispersion- and transport-based mixing diagnostics are defined, compared and contrasted to motivate the computation and use of effective diffusivity derived from Lagrangian particles. We compute the effective diffusivity by first performing scalar transport on Lagrangian control areas using stored trajectories computed from online Lagrangian In-situ Global High-performance particle Tracking (LIGHT) using the Model for Prediction Across Scales Ocean (MPAS-O). Furthermore, the Lagrangian scalar transport scheme is compared against an Eulerian scalar transport scheme. Spatially-variable effective diffusivities are computed from resulting time-varying cumulative concentrations that vary as a function of cumulative area. The transport-based Eulerian and Lagrangian effective diffusivity diagnostics are found to be qualitatively consistent with the dispersion-based diffusivity. All diffusivity estimates show a region of increased subsurface diffusivity within the core of an idealized circumpolar current and results are within a factor of two of each other. The Eulerian and Lagrangian effective diffusivities are most similar; smaller and more spatially diffused values are obtained with the dispersion-based diffusivity computed with particle clusters.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Fast self-diffusion of ions in CH 3 NH 3 PbI 3 : the interstiticaly mechanism versus vacancy-assisted mechanism

The stability of organic-inorganic halide perovskites is a major challenge for their applications and has been extensively studied. Among the possible underlying reasons, ion self-diffusion has been inferred to play important roles. While theoretical studies congruously support that iodine is more mobile, experimental studies only observe the direct diffusion of the MA ion and possible diffusion of iodine. The discrepancy may result from the incomplete understanding of ion diffusion mechanisms. With the help of first-principles calculations, we studied ion diffusion in CH3NH3PbI3 (MAPbI3) through not only the vacancy-assisted mechanisms presumed in previous theoretical studies, but also the neglected interstiticaly mechanisms. We found that compared to the diffusion through the vacancy-assisted mechanism, MA ion diffusion through the interstiticaly mechanism has a much smaller barrier which could explain experimental observations. For iodine diffusion, both mechanisms can yield relatively small barriers. Depending on the growth conditions, defect densities of vacancies and interstitials can vary and so do the diffusion species as well as diffusion mechanisms. Our work thus supports that both MA and iodine ion diffusion could contribute to the performance instability of MAPbI3. While being congruous with experimental results, our work fills the research gap by providing a full understanding of ion diffusion in halide perovskites.

14 SOLAR ENERGY↗
DOE OSTI2014

Lithium ion diffusion in Li β-alumina single crystals measured by pulsed field gradient NMR spectroscopy

The lithium ion diffusion coefficient of a 93% Li β-alumina single crystal was measured for the first time using pulsed field gradient (PFG) NMR spectroscopy with two different crystal orientations. The diffusion coefficient was found to be 1.2 × 10{sup −11} m{sup 2}/s in the direction perpendicular to the c axis at room temperature. The Li ion diffusion coefficient along the c axis direction was found to be very small (6.4 × 10{sup −13} m{sup 2}/s at 333 K), which suggests that the macroscopic diffusion of the Li ion in the β-alumina crystal is mainly two-dimensional. The diffusion coefficient for the same sample was also estimated using NMR line narrowing data and impedance measurements. The impedance data show reasonable agreement with PFG-NMR data, while the line narrowing measurements provided a lower value for the diffusion coefficient. Line narrowing measurements also provided a relatively low value for the activation energy and pre-exponential factor. The temperature dependent diffusion coefficient was obtained in the temperature range 297–333 K by PFG-NMR, from which the activation energy for diffusion of the Li ion was estimated. The activation energy obtained by PFG-NMR was smaller than that obtained by impedance measurements, which suggests that thermally activated defect formation energy exists for 93% Li β-alumina single crystals. The diffusion time dependence of the diffusion coefficient was observed for the Li ion in the 93% Li β-alumina single crystal by means of PFG-NMR experiments. Motion of Li ion in fractal dimension might be a possible explanation for the observed diffusion time dependence of the diffusion coefficient in the 93% Li β–alumina system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗