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

TU-AB-BRA-07: Distortion-Free 3D Diffusion MRI On An MRI-Guided Radiotherapy System for Longitudinal Tumor Response Assessment

Purpose: To develop a reliable, 3D distortion-free diffusion MRI technique for longitudinal tumor response assessment and MRI-guided adaptive radiotherapy(RT). Methods: A diffusion prepared 3D turbo spin echo readout (DP-TSE) sequence was developed and compared with the conventional diffusion-weighted single-shot echo-planar-imaging (DW-ssEPI) sequence in a commercially available diffusion phantom, and one head-and-neck and one brain cancer patient on an MRI-guided RT system (ViewRay). In phantom study, the geometric fidelity was quantified as the ratio between the left-right (RL) and anterior-posterior (AP) dimension. Ten slices were measured on DP-TSE, DW-ssEPI and standard TSE images where the later was used as the geometric reference. ADC accuracy was verified at both 0°C (reference ADC available) and room temperature with a range of diffusivity between 0.35 and 2.0*10{sup −3}mm{sup 2}/s. The ADC reproducibility was assessed based on 8 room-temperature measurements on 6 different days. In the pilot single-slice in-vivo study, CT images were used as the geometric reference, and ADC maps from both diffusion sequences were compared. Results: Distortion and susceptive-related artifact were severe in DW-ssEPI, with significantly lower RL/AP ratio (0.9579±0.0163) than DP-TSE (0.9990±0.0031) and TSE (0.9995±0.0031). ADCs from the two diffusion sequences both matched well with the vendor-provided values at 0°C; however DW-ssEPI fails to provide accurate ADC for high diffusivity vials at room temperature due to high noise level (10 times higher than DP-TSE). The DP-TSE sequence had excellent ADC reproducibility with <4% ADC variation among 8 separate measurements. In patient study, DP-TSE exhibited substantially improved geometric reliability. ROI analysis in ADC maps generated from DP-TSE and DW-ssEPI showed <5% difference where high b-value images were excluded from the latter approach due to excessive noise level. Conclusion: A diffusion MRI sequence with excellent geometric fidelity, accurate and highly reproducible ADC measurements was proposed for longitudinal tumor response assessment using an MRI-guided RT system. Yu Gao acknowledges research support from ViewRay.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

Ion Diffusivity through the Solid Electrolyte Interphase in Lithium-Ion Batteries

Understanding the transport properties of the solid electrolyte interface (SEI) is a critical piece in the development of lithium ion batteries (LIB) with better performance. We studied the lithium ion diffusivity in the main components of the SEI found in LIB with silicon anodes and performed classical molecular dynamics (MD) simulations on lithium fluoride (LiF), lithium oxide (Li 2 O) and lithium carbonate (Li 2 CO 3 ) in order to provide insights and to calculate the diffusion coefficients of Li-ions at temperatures in the range of 250 K to 400 K, which is within the LIB operating temperature range. We find a slight increase in the diffusivity as the temperature increases and since diffusion is noticeable at high temperatures, Li-ion diffusion in the range of 130 to 1800 K was also studied and the diffusion mechanisms involved in each SEI compound were analyzed. We observed that the predominant mechanisms of Li-ion diffusion included vacancy assisted and knock-off diffusion in LiF, direct exchange in Li 2 O, and vacancy and knock-off in Li 2 CO 3 . Moreover, we also evaluated the effect of applied electric fields in the diffusion of Li-ions at room temperature.

25 ENERGY STORAGE↗
DOE OSTI2016

Report on simulation of fission gas and fission product diffusion in UO 2

In UO 2 nuclear fuel, the retention and release of fission gas atoms such as xenon (Xe) are important for nuclear fuel performance by, for example, reducing the fuel thermal conductivity, causing fuel swelling that leads to mechanical interaction with the clad, increasing the plenum pressure and reducing the fuel–clad gap thermal conductivity. We use multi-­scale simulations to determine fission gas diffusion mechanisms as well as the corresponding rates in UO 2 under both intrinsic and irradiation conditions. In addition to Xe and Kr, the fission products Zr, Ru, Ce, Y, La, Sr and Ba have been investigated. Density functional theory (DFT) calculations are used to study formation, binding and migration energies of small clusters of Xe atoms and vacancies. Empirical potential calculations enable us to determine the corresponding entropies and attempt frequencies for migration as well as investigate the properties of large clusters or small fission gas bubbles. A continuum reaction-­diffusion model is developed for Xe and point defects based on the mechanisms and rates obtained from atomistic simulations. Effective fission gas diffusivities are then obtained by solving this set of equations for different chemical and irradiation conditions using the MARMOT phase field code. The predictions are compared to available experimental data. The importance of the large Xe U3O cluster (a Xe atom in a uranium + oxygen vacancy trap site with two bound uranium vacancies) is emphasized, which is a consequence of its high mobility and high binding energy. We find that the Xe U3O cluster gives Xe diffusion coefficients that are higher for intrinsic conditions than under irradiation over a wide range of temperatures. Under irradiation the fast-­moving Xe U3O cluster recombines quickly with irradiation-induced interstitial U ions, while this mechanism is less important for intrinsic conditions. The net result is higher concentration of the Xe U3O cluster for intrinsic conditions than under irradiation. We speculate that differences in the irradiation conditions and their impact on the Xe U3O cluster can explain the wide range of diffusivities reported in experimental studies. However, all vacancy-­mediated mechanisms underestimate the Xe diffusivity compared to the empirical radiation-­enhanced rate used in most fission gas release models. We investigate the possibility that diffusion of small fission gas bubbles or extended Xe-­vacancy clusters may give rise to the observed radiation-­enhanced diffusion coefficient. These studies highlight the importance of U divacancies and an octahedron coordination of uranium vacancies encompassing a Xe fission gas atom. The latter cluster can migrate via a multistep mechanism with a rather low effective barrier, which together with irradiation-induced clusters of uranium vacancies, gives rise to the irradiation-enhanced diffusion coefficient observed in experiments.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Theoretical Analysis of Experimental Data of Sodium Diffusion in Oxidized Molybdenum Thin Films

In this work, the diffusion process of sodium (Na) in molybdenum (Mo) thin films while it was deposited on soda lime glass (SLG) was studied. A small amount of oxygen was present in the chamber while the direct-current (DC) magnetron sputtering was used for the deposition. The substrate temperatures were varied to observe its effect. Such molybdenum films, with or without oxidations, are often used in thin film solar cells, either as back contact or as hole transport layers. Secondary ion mass spectrometry (SIMS) was used to quantify the concentration of the species. A grain diffusion mechanistic model incorporating the effect of grain and grain boundary geometrical shape and size was developed. The model was used to provide an in-depth theoretical analysis of the sodium diffusion in molybdenum thin films that lead to the measured SIMS data. It was observed that not only diffusion coefficients should be considered when analyzing diffusion processes in thin films but also the ratio of grain boundary size to grain size. Both depend on substrate temperature and directly affect the amount of diffused species in the film. The data were analyzed under the light of the film growth speed versus diffusion front speed, the effect of oxygen content, and the effect of substrate temperature on the overall diffusion process. The temperature inversely affects the ratio of grain boundary size and grain size and directly affects the diffusion coefficient, which leads to a preferable temperature at which the highest amount of alkali can be found in the film.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Perpendicular Diffusion of Solar Energetic Particles: Model Results and Implications for Electrons

The processes responsible for the effective longitudinal transport of solar energetic particles (SEPs) are still not completely understood. We address this issue by simulating SEP electron propagation using a spatially 2D transport model that includes perpendicular diffusion. By implementing, as far as possible, the most reasonable estimates of the transport (diffusion) coefficients, we compare our results, in a qualitative manner, to recent observations at energies of 55–105 keV, focusing on the longitudinal distribution of the peak intensity, the maximum anisotropy, and the onset time. By using transport coefficients that are derived from first principles, we limit the number of free parameters in the model to (i) the probability of SEPs following diffusing magnetic field lines, quantified by a∈[0,1], and (ii) the broadness of the Gaussian injection function. It is found that the model solutions are extremely sensitive to the magnitude of the perpendicular diffusion coefficient and relatively insensitive to the form of the injection function as long as a reasonable value of a = 0.2 is used. We illustrate the effects of perpendicular diffusion on the model solutions and discuss the viability of this process as a dominant mechanism by which SEPs are transported in longitude. Lastly, we try to quantity the effectiveness of perpendicular diffusion as an interplay between the magnitude of the relevant diffusion coefficient and the SEP intensity gradient driving the diffusion process. It follows that perpendicular diffusion is extremely effective early in an SEP event when large intensity gradients are present, while the effectiveness quickly decreases with time thereafter.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Improved diffusion coefficients generated from Monte Carlo codes

Monte Carlo codes are becoming more widely used for reactor analysis. Some of these applications involve the generation of diffusion theory parameters including macroscopic cross sections and diffusion coefficients. Two approximations used to generate diffusion coefficients are assessed using the Monte Carlo code MC21. The first is the method of homogenization; whether to weight either fine-group transport cross sections or fine-group diffusion coefficients when collapsing to few-group diffusion coefficients. The second is a fundamental approximation made to the energy-dependent P1 equations to derive the energy-dependent diffusion equations. Standard Monte Carlo codes usually generate a flux-weighted transport cross section with no correction to the diffusion approximation. Results indicate that this causes noticeable tilting in reconstructed pin powers in simple test lattices with L2 norm error of 3.6%. This error is reduced significantly to 0.27% when weighting fine-group diffusion coefficients by the flux and applying a correction to the diffusion approximation. Noticeable tilting in reconstructed fluxes and pin powers was reduced when applying these corrections. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2018

Effect of Weld Consumable Conditioning on the Diffusible Hydrogen and Subsequent Residual Stress and Flexural Strength of Multipass Welded P91 Steels

P91 steel weld joint was prepared using the shielded metal arc welding process and four different conditions of weld consumable that provide the different levels of diffusible hydrogen in deposited metal. In the present research, the effects of diffusible hydrogen content on the flexural strength, lower critical stress, and tensile strength of P91 steel welds were also determined with respect to different electrode conditions. To investigate the effect of diffusible hydrogen on multipass welding, top and root side flexural tests were performed. The residual stresses (axial stress and transverse stress) were also measured using the blind hole drilling method for different conditions of welding consumable. The peak value of residual stresses was measured at the center of the weld fusion zone. The maximum value of transverse stress was measured to be 355 MPa for case II (6.21 mL/100 g of diffusible hydrogen), while the maximum axial stress was about 218 MPa for case IV (12.43 mL/100 g of diffusible hydrogen). A three-dimensional finite element simulation was also performed to predict the residual stress distribution and thermal profile along the welded joint. The experimentally determined residual stresses correlated well with the numerically estimated residual stresses. The diffusible hydrogen content was not observed to have any significant effect on the residual stresses. The corrected residual stress values were also predicted by considering the plasticity-induced error. However, the flexural performance of the welded joint was affected by the diffusible hydrogen content. The top and root flexural strength was measured to be optimum for the low level of diffusible hydrogen content, and the values decreased with an increase in diffusible hydrogen content.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Study on impact of Cr and Mo on diffusion of H in 2.25Cr1Mo steel using first-principle calculations

Chrome-molybdenum steel (2.25Cr1Mo steel) is one of the main materials of a steam generator (SG) in high-temperature reactor-pebblebed modules (HTR-PM). It is essential to analyze the source term of tritium in this material, because the behavior of tritium in a SG is important for performing source term analysis in normal and accident conditions. In this article, the diffusion behavior of H atom in 2.25Cr1Mo steel was calculated to estimate the diffusivity of tritium using first-principles density functional theory. To develop and simplify the model of hydrogen diffusion in 2.25Cr1Mo steel, the impact of Cr and Mo on the diffusion of hydrogen in bcc-Fe were first calculated, all the possible diffusion paths were considered, and the minimum energy path was obtained. The diffusion activation energy and pre-exponential factor of the diffusion coefficient were obtained from Vienna Ab initio Simulation Package combined with the climbing image-nudged elastic band method. Here the results indicate that the minimum energy path for the impurity H atom is from one tetrahedral interstitial site to an adjacent tetrahedral interstitial site. The function of the diffusion coefficient of H in 2.25Cr1Mo steel with temperature T can be expressed as D = 1.486 X 10 -7 x ($\frac{-16.350kJ/mol}{RT}$) (M2/s). The diffusion coefficient of our calculation and some of the previous experiments have an excellent quantitative agreement, which indicates the reliability of our crystalline model and the practicability of the present theoretical approach. More importantly, the computational results in this work can be treated as a good screening method to collect reasonable experimental data, which will provide a good reference for tritium source term evaluation in the SG of the HTR-PM.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Long-term diffusion of U(VI) in bentonite: Dependence on density

As a contribution to the safety assessment of nuclear waste repositories, U(VI) diffusion through the potential buffer material MX-80 bentonite was investigated at three clay dry densities over six years. Synthetic MX-80 model pore water was used as background electrolyte. Speciation calculations showed that Ca 2 UO 2 (CO 3 ) 3 (aq) was the main U(VI) species. The in- and out-diffusion of U(VI) was investigated separately. U(VI) diffused about 3 mm, 1.5 mm, and 1 mm into the clay plug at ρ = 1.3, 1.6, and 1.9 g/cm 3 , respectively. No through-diffusion of the U(VI) tracer was observed. However, leaching of natural uranium contained in the clay occurred and uranium was detected in all receiving reservoirs. As expected, the effective and apparent diffusion coefficients, D e and D a , decreased with increasing dry density. The D a values for the out-diffusion of natural U(VI) were in good agreement with previously determined values. Surprisingly, D a values for the in-diffusion of U(VI) were about two orders of magnitude lower than values obtained in short-term in-diffusion experiments reported in the literature. Some potential reasons for this behavior that were evaluated are changes of the U(VI) speciation within the clay (precipitation, reduction) or changes of the clay porosity and pore connectivity with time. By applying Archie's law and the extended Archie's law, it was estimated that a significantly smaller effective porosity must be present for the long-term in-diffusion of U(VI). Finally, the results suggest that long-term studies of key transport phenomena may reveal additional processes that can directly impact long-term repository safety assessments.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2021

Artificial neural network prediction of self-diffusion in pure compounds over multiple phase regimes

Artificial neural networks (ANNs) were developed to accurately predict the self-diffusion constants for pure components in liquid, gas and super critical phases. The ANNs were tested on an experimental database of 6625 self-diffusion constants for 118 different chemical compounds. The presence of multiple phases results in a heavy skew in the distribution of diffusion constants and multiple approaches were used to address this challenge. First, an ANN was developed with the raw diffusion values to assess what the main drawbacks of this direct method were. The first approach for improving the predictions involved taking the log 10 of diffusion to provide a more uniform distribution and reduce the range of target output values used to develop the ANN. The second approach involved developing individual ANNs for each phase using the raw diffusion values. Results show that the log transformation leads to a model with the best self-diffusion constant predictions and an overall average absolute deviation (AAD) of 6.56%. The resultant ANN is a generalized model that can be used to predict diffusion across all three phases and over a diverse group of compounds. The importance of each input feature was ranked using a feature addition method revealing that the density of the compound has the largest impact on the ANN prediction of self-diffusion constants in pure compounds.

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

Anisotropic Diffusion of a Charged Tritium Interstitial in Li 2 Ti O 3 from First-Principles Calculations

Tritium diffusion in the nuclear fusion reactor breeder material Li 2 TiO 3 is a fundamentally important process for the tritium release kinetics. The energy barrier of tritium diffusion in Li 2 TiO 3 was reported with a considerable uncertainty in previous experiments. Here, we perform the systematic density-functional-theory (DFT) studies for the diffusion processes of a positively-charged tritium, which is the most preferable charge state of the tritium interstitial in a single Li 2 TiO 3 crystal. By calculating various local-diffusion minimum-energy paths, we find that the diffusion of the tritium is strongly anisotropic along different crystalline directions. The most favorable diffusion paths appear within a Li 6 atomic single layer of Li 2 TiO 3 and the corresponding DFT diffusion barrier is 0.281 eV, while the diffusion barrier for most favorable diffusion paths crossing a Li 2 Ti 4 atomic layer is 0.868 eV.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Impact of Alloy-Disorder-Induced Localization on Hole Diffusion in Highly Excited c -Plane and m -Plane (In, Ga) N Quantum Wells

The diffusion coefficient of holes can provide knowledge about carrier localization in (In, Ga)N, where the carrier dynamics are altered by randomly fluctuating potential landscape. In group-III nitrides, the diffusivity of holes is difficult to measure by electrical methods but it can be studied using optical techniques. Here, we investigate the dependence of the hole diffusion coefficient on direction and carrier density in c-plane and m-plane (In, Ga)N structures by employing the light-induced transient-grating technique. We show that the hole diffusion coefficient is anisotropic in the m-plane structure, where it is several times larger along the a crystallographic direction than along the c direction. Such anisotropy is observed within the broad range of carrier densities from 10 18 to 10 20 cm –3 . The diffusivity changes nonmonotonously with increasing photoexcitation, this dependence being different in thick and thin layers. We argue that an unexpectedly high diffusion coefficient at low carrier densities in thick quantum wells can be a signature of efficient hole transport via percolative paths occurring due to compositional disorder. In turn, a decrease of diffusivity with the excitation can reflect the effect of Coulomb blockade of these paths. Lastly, we demonstrate that disorder impacts carrier diffusivity even at carrier densities above 10 19 cm –3 , where the overflow of localized states must be included to explain the observed increase of the diffusion coefficient with the carrier density.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Diffusion in Jammed Particle Packs

Using random walk simulations we explore diffusive transport through monodisperse sphere packings over a range of packing fractions Φ in the vicinity of the jamming transition at Φ c . Various diffusion properties are computed over several orders of magnitude in both time and packing pressure. Two well-separated regimes of normal “Fickian” diffusion, where the mean squared displacement is linear in time, are observed. The first corresponds to diffusion inside individual spheres, while the latter is the long-time bulk diffusion. The intermediate anomalous diffusion regime and the long-time value of the diffusion coefficient are both shown to be controlled by particle contacts, which in turn depend on proximity to Φ c . The time required to recover normal diffusion t * scales as (Φ–Φ c ) –0.5 and the long-time diffusivity D∞~(Φ–Φ c ) 0.5 , or D ∞ ~1/ t *. Here, it is shown that the distribution of mean first passage times associated with the escape of random walkers between neighboring particles controls both t * and D∞ in the limit Φ→Φ c .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

Characterization of Grain‐Size Distribution, Thermal Conductivity, and Gas Diffusivity in Variably Saturated Binary Sand Mixtures

Core Ideas Parameterized models are relevant to applications where different porous media mixtures are used. Models include a grain‐size distribution function to describe bimodal behavior for binary mixtures. Improved model describes observed thermal conductivity–saturation relations for binary mixtures. Combined Buckingham–Penman model used to describe observed gas diffusivity–air content relations. This work is relevant for proper simulation of mixed porous media. Characterization of differently textured porous materials, as well as different volumetric porous media mixtures, in relation to mass and heat transport is vital for many engineering and research applications. Functional relations describing physical (e.g., grain‐size distribution, total porosity), thermal, and gas diffusion properties of porous media and mixtures are necessary to optimize the design of porous systems that involve heat and gas transport processes. However, only a limited number of studies provide characterization of soil physical, thermal, and gas diffusion properties and the functional relationships of these properties under varying soil water contents, especially for soil mixtures, complicating optimization efforts. To better understand how mixing controls the physical, thermal, and gas diffusion properties of porous media, a set of laboratory experiments was performed using five volumetric mixtures of coarse‐ and fine‐grained sand particles. For each mixture, the grain‐size distribution (GSD), thermal conductivity, and gas diffusivity were obtained and parameterized using existing and suggested parametric models. Results show that the extended, two‐region Rosin–Rammler particle‐size distribution model proposed in this study could successfully describe the bimodal behavior of the GSD of binary mixtures. Further, the modified Côté and Konrad thermal conductivity model adequately described the thermal conductivity–water saturation relations observed in different mixtures. The proposed simple soil‐gas diffusivity descriptive model parameterized the upper limit, average, and lower limit behavior in gas diffusivity–air content relations in apparently texture‐invariant gas diffusivity data. Results further show a close analogy between gas diffusivity and thermal conductivity and their variation with saturation across different binary mixtures. Overall, the results of the study provide useful numerical insight into the physical, thermal, and gas transport characteristics of binary mixtures, with wide implications for future engineering and research applications that involve multicomponent porous systems.

Deepagoda, T. K. K. Chamindu↗
DOE OSTI2017

An efficient approach for treating composition-dependent diffusion within organic particles

Mounting evidence demonstrates that under certain conditions the rate of component partitioning between the gas and particle phase in atmospheric organic aerosol is limited by particle-phase diffusion. To date, however, particle-phase diffusion has not been incorporated into regional atmospheric models. An analytical rather than numerical solution to diffusion through organic particulate matter is desirable because of its comparatively small computational expense in regional models. Current analytical models assume diffusion to be independent of composition and therefore use a constant diffusion coefficient. To realistically model diffusion, however, it should be composition-dependent (e.g. due to the partitioning of components that plasticise, vitrify or solidify). This study assesses the modelling capability of an analytical solution to diffusion corrected to account for composition dependence against a numerical solution. Results show reasonable agreement when the gas-phase saturation ratio of a partitioning component is constant and particle-phase diffusion limits partitioning rate (<10% discrepancy in estimated radius change). However, when the saturation ratio of the partitioning component varies, a generally applicable correction cannot be found, indicating that existing methodologies are incapable of deriving a general solution. Until such time as a general solution is found, caution should be given to sensitivity studies that assume constant diffusivity. Furthermore, the correction was implemented in the polydisperse, multi-process Model for Simulating Aerosol Interactions and Chemistry (MOSAIC) and is used to illustrate how the evolution of number size distribution may be accelerated by condensation of a plasticising component onto viscous organic particles.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Viscosities, diffusion coefficients, and mixing times of intrinsic fluorescent organic molecules in brown limonene secondary organic aerosol and tests of the Stokes–Einstein equation

Viscosities and diffusion rates of organics within secondary organic aerosol (SOA) remain uncertain. Using the bead-mobility technique, we measured viscosities as a function of water activity ( a w ) of SOA generated by the ozonolysis of limonene followed by browning by exposure to NH 3 (referred to as brown limonene SOA or brown LSOA). These measurements together with viscosity measurements reported in the literature show that the viscosity of brown LSOA increases by 3–5 orders of magnitude as the a w decreases from 0.9 to approximately 0.05. In addition, we measured diffusion coefficients of intrinsic fluorescent organic molecules within brown LSOA matrices using rectangular area fluorescence recovery after photobleaching. Based on the diffusion measurements, as the a w decreases from 0.9 to 0.33, the average diffusion coefficient of the intrinsic fluorescent organic molecules decreases from 5.5 × 10 - 9 to 7.1 × 10 - 13 cm 2 s -1 and the mixing times of intrinsic fluorescent organic molecules within 200 nm brown LSOA particles increases from 0.002 to 14 s. These results suggest that the mixing times of large organics in the brown LSOA studied here are short ( < 1 h) for a w and temperatures often found in the planetary boundary layer (PBL). Since the diffusion coefficients and mixing times reported here correspond to SOA generated using a high mass loading ( ~1000 µ g m -3 ), biogenic SOA particles found in the atmosphere with mass loadings ≤10 µ g m -3 are likely to have higher viscosities and longer mixing times (possibly 3 orders of magnitude longer). These new measurements of viscosity and diffusion were used to test the accuracy of the Stokes–Einstein relation for predicting diffusion rates of organics within brown LSOA matrices. The results show that the Stokes–Einstein equation gives accurate predictions of diffusion coefficients of large organics within brown LSOA matrices when the viscosity of the matrix is as high as 10 2 to 10 4 Pa s. These results have important implications for predicting diffusion and mixing within SOA particles in the atmosphere.

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

Diffusion kinetics of 3 He in pyroxene and plagioclase and applications to cosmogenic exposure dating and paleothermometry in mafic rocks

In this study, we investigate the diffusivity of cosmogenic 3 He in a variety of plagioclase and pyroxene compositions, and its application to paleothermometry and exposure dating in these minerals, through stepwise degassing experiments. While cosmogenic 3 He has been utilized for exposure dating in pyroxene for decades due to its retentivity, plagioclase, often found along with pyroxene in mafic rocks, is generally less retentive of cosmogenic noble gas. However, the diffusivity of 3 He in either plagioclase or pyroxene has not yet been measured quantitatively. A challenge in measuring diffusion kinetics by step-degassing experiments in poorly retentive minerals is the fact that significant amounts of He can be lost prior to the experiment. To address this issue, we apply a forward “multiple diffusion domain” (MDD) inversion model that includes model predictions of initial gas loss during irradiation and storage of the samples to account for this observation and add constraints to the diffusion parameters. We find that 3 He diffusivity in plagioclase appears to be highly variable. This variability can be explained by the MDD inversion models' inability to constrain the diffusion parameters when significant gas has been lost during irradiation and/or prolonged storage prior to experiment analysis, resulting in an overestimation of 3 He retentivity. Plagioclase samples that were kept frozen after irradiation, to limit this initial gas loss, yielded the most reliable estimate of diffusion kinetics. We find that 3 He in plagioclase is diffusively lost at Earth's surface temperatures on a timescale of 100 years and is therefore unsuitable for surface temperature paleothermometry. On the contrary, we find cosmogenic 3 He in pyroxene to be retentive at Earth's surface temperatures on a 1-million-year timescale.

Bergelin, Marie [Berkeley Geochronology Center (BG↗
DOE OSTI2019

A Prediction of Uranium(VI) Diffusion in Montmorillonite at Various Chemical Solution Conditions

Most nuclear waste disposal options include compacted bentonite, predominately consisting of montmorillonite clay, as part of a barrier system to minimize contaminant mobility. Uranium (U) is the primary element in spent nuclear fuel, and a potential contaminant of water resources. Uranium sorption onto clay and its slow diffusive transport away from waste canisters are expected to limit U(VI) mobility, and possibly also control nuclear fuel degradation rates based on mass action considerations. A prediction of U(VI) diffusion in these systems is complicated by two main factors: (1) the complex mineralogical structure of montmorillonite clay, leading to two types of clay porosities and surface sites, and (2) the dependence of U(VI) sorption on U(VI) solution speciation and the specific chemical system conditions. For instance, depending on pH and ionic strength, U(VI) may sorb onto montmorillonite due to weak ion exchange reactions or the formation of stable surface complexes. Furthermore, in contrast to conventional porous media, U(VI) diffusion in montmorillonite is not only limited by contaminant sorption, but also by a (potential) lack of access to the full clay porosity for some U(VI) solution species. For instance, at the high degrees of clay compaction expected in future barrier systems, the partial or full exclusion of U(VI) from negatively-charged clay interlayer spaces may substantially decrease U(VI) fluxes. In this modeling study, we specifically investigated how the charge of U(VI) solution species affects U(VI) sorption and diffusion behavior. For this purpose, we first simplified complex chemical solution speciation diagrams to charge speciation diagrams. Then, we coupled an existing U(VI) surface complexation model (Tournassat et al., Geochim. Cosmochim. Acta, 2018) with an analytical solution of the diffusion equation to simulate U(VI) mobility in lab-scale through-diffusion experiments for various chemical conditions. Our results suggest that there are two important parameters that govern U(VI) diffusion in these systems. First, U(VI) surface complexation on clay edge surfaces leads to contaminant retardation. On the field scale, this will determine how early a plume may cross site boundaries and pose risks to drinking water resources. Second, the degree of access of U(VI) solution species to diffuse layers in clay interlayer spaces largely controls U(VI) diffusive fluxes. In the field, this parameter will govern the level of U(VI) concentrations that can be expected in contaminated waters.

Pistorino, Jonathan C.↗