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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

Laser Velocimeter Measurements of the Flow Within a Vane-Island Diffuser

A series of detailed non-intrusive laser velocimeter flow field surveys were performed to validate computational fluid dynamic (CFD) codes for rotating pump components. Three component laser-2-focus (L2F) laser velocimetry was used to acquire data within a flow passage of a diffuser at design flow conditions. The thirteen vane-island type diffuser was designed to allow diffusion in the axial direction. Multiple laser velocimeter surveys were performed throughout several diffuser passages and at multiple diffuser depths (hub to shroud). The axial (3D) component of the velocity at all locations was small. The radial and meridional velocity components nearest the diffuser hub increased in magnitude relative to the flow nearest the shroud as the flow progressed through the diffuser. A continuity check across the diffuser throat and based on the meridional velocity component yielded a match with a facility flow meter of 98.0%.

Ferguson, Thomas V.↗
NASA NTRS2001

Suppression of Soot Formation and Shapes of Laminar Jet Diffusion Flames

Laminar nonpremixed (diffusion) flames are of interest because they provide model flame systems that are far more tractable for analysis and experiments than practical turbulent flames. In addition, many properties of laminar diffusion flames are directly relevant to turbulent diffusion flames using laminar flamelet concepts. Finally, laminar diffusion flame shapes have been of interest since the classical study of Burke and Schumann because they involve a simple nonintrusive measurement that is convenient for evaluating flame shape predictions. Motivated by these observations, the shapes of round hydrocarbon-fueled laminar jet diffusion flames were considered, emphasizing conditions where effects of buoyancy are small because most practical flames are not buoyant. Earlier studies of shapes of hydrocarbon-fueled nonbuoyant laminar jet diffusion flames considered combustion in still air and have shown that flames at the laminar smoke point are roughly twice as long as corresponding soot-free (blue) flames and have developed simple ways to estimate their shapes. Corresponding studies of hydrocarbon-fueled weakly-buoyant laminar jet diffusion flames in coflowing air have also been reported. These studies were limited to soot-containing flames at laminar smoke point conditions and also developed simple ways to estimate their shapes but the behavior of corresponding soot-free flames has not been addressed. This is unfortunate because ways of selecting flame flow properties to reduce soot concentrations are of great interest; in addition, soot-free flames are fundamentally important because they are much more computationally tractable than corresponding soot-containing flames. Thus, the objectives of the present investigation were to observe the shapes of weakly-buoyant laminar jet diffusion flames at both soot-free and smoke point conditions and to use the results to evaluate simplified flame shape models. The present discussion is brief.

Xu, F.↗
NASA NTRS1998

Diffusion Characteristics of Upwind Schemes on Unstructured Triangulations

The diffusive characteristics of two upwind schemes, multi-dimensional fluctuation splitting and dimensionally-split finite volume, are compared for scalar advection-diffusion problems. Algorithms for the two schemes are developed for node-based data representation on median-dual meshes associated with unstructured triangulations in two spatial dimensions. Four model equations are considered: linear advection, non-linear advection, diffusion, and advection-diffusion. Modular coding is employed to isolate the effects of the two approaches for upwind flux evaluation, allowing for head-to-head accuracy and efficiency comparisons. Both the stability of compressive limiters and the amount of artificial diffusion generated by the schemes is found to be grid-orientation dependent, with the fluctuation splitting scheme producing less artificial diffusion than the dimensionally-split finite volume scheme. Convergence rates are compared for the combined advection-diffusion problem, with a speedup of 2-3 seen for fluctuation splitting versus finite volume when solved on the same mesh. However, accurate solutions to problems with small diffusion coefficients can be achieved on coarser meshes using fluctuation splitting rather than finite volume, so that when comparing convergence rates to reach a given accuracy, fluctuation splitting shows a 20-25 speedup over finite volume.

Wood, William A.↗
NASA NTRS1996

Microgravity Turbulent Gas-Jet Diffusion Flames

A gas-jet diffusion flame is similar to the flame on a Bunsen burner, where a gaseous fuel (e.g., propane) flows from a nozzle into an oxygen-containing atmosphere (e.g., air). The difference is that a Bunsen burner allows for (partial) premixing of the fuel and the air, whereas a diffusion flame is not premixed and gets its oxygen (principally) by diffusion from the atmosphere around the flame. Simple gas-jet diffusion flames are often used for combustion studies because they embody the mechanisms operating in accidental fires and in practical combustion systems. However, most practical combustion is turbulent (i.e., with random flow vortices), which enhances the fuel/air mixing. These turbulent flames are not well understood because their random and transient nature complicates analysis. Normal gravity studies of turbulence in gas-jet diffusion flames can be impeded by buoyancy-induced instabilities. These gravitycaused instabilities, which are evident in the flickering of a candle flame in normal gravity, interfere with the study of turbulent gas-jet diffusion flames. By conducting experiments in microgravity, where buoyant instabilities are avoided, we at the NASA Lewis Research Center hope to improve our understanding of turbulent combustion. Ultimately, this could lead to improvements in combustor design, yielding higher efficiency and lower pollutant emissions. Gas-jet diffusion flames are often researched as model flames, because they embody mechanisms operating in both accidental fires and practical combustion systems (see the first figure). In normal gravity laboratory research, buoyant air flows, which are often negligible in practical situations, dominate the heat and mass transfer processes. Microgravity research studies, however, are not constrained by buoyant air flows, and new, unique information on the behavior of gas-jet diffusion flames has been obtained.

Source record↗
NASA NTRS1999

Effects of g-Jitter on Diffusion in Binary Liquids

The microgravity environment offers the potential to measure the binary diffusion coefficients in liquids without the masking effects introduced by buoyancy-induced flows due to Earth s gravity. However, the background g-jitter (vibrations from the shuttle, onboard machinery, and crew) normally encountered in many shuttle experiments may alter the benefits of the microgravity environment and introduce vibrations that could offset its intrinsic advantages. An experiment during STS-85 (August 1997) used the Microgravity Vibration Isolation Mount (MIM) to isolate and introduce controlled vibrations to two miscible liquids inside a cavity to study the effects of g-jitter on liquid diffusion. Diffusion in a nonhomogeneous liquid system is caused by a nonequilibrium condition that results in the transport of mass (dispersion of the different kinds of liquid molecules) to approach equilibrium. The dynamic state of the system tends toward equilibrium such that the system becomes homogeneous. An everyday example is the mixing of cream and coffee (a nonhomogeneous system) via stirring. The cream diffuses into the coffee, thus forming a homogeneous system. At equilibrium the system is said to be mixed. However, during stirring, simple observations show complex flow field dynamics-stretching and folding of material interfaces, thinning of striation thickness, self-similar patterns, and so on. This example illustrates that, even though mixing occurs via mass diffusion, stirring to enhance transport plays a major role. Stirring can be induced either by mechanical means (spoon or plastic stirrer) or via buoyancy-induced forces caused by Earth s gravity. Accurate measurements of binary diffusion coefficients are often inhibited by buoyancy-induced flows. The microgravity environment minimizes the effect of buoyancy-induced flows and allows the true diffusion limit to be achieved. One goal of this experiment was to show that the microgravity environment suppresses buoyancy-induced convection, thereby mass diffusion becomes the dominant mechanism for transport. Since g-jitter transmitted by the shuttle to the experiment can potentially excite buoyancy-induced flows, we also studied the effects of controlled vibrations on the system.

Duval, Walter M. B.↗
NASA NTRS2013

Generating a Simulated Fluid Flow Over an Aircraft Surface Using Anisotropic Diffusion

A fluid-flow simulation over a computer-generated aircraft surface is generated using a diffusion technique. The surface is comprised of a surface mesh of polygons. A boundary-layer fluid property is obtained for a subset of the polygons of the surface mesh. A pressure-gradient vector is determined for a selected polygon, the selected polygon belonging to the surface mesh but not one of the subset of polygons. A maximum and minimum diffusion rate is determined along directions determined using a pressure gradient vector corresponding to the selected polygon. A diffusion-path vector is defined between a point in the selected polygon and a neighboring point in a neighboring polygon. An updated fluid property is determined for the selected polygon using a variable diffusion rate, the variable diffusion rate based on the minimum diffusion rate, maximum diffusion rate, and angular difference between the diffusion-path vector and the pressure-gradient vector.

Rodriguez, David L.↗
DOE OSTI2014

Atomistic modeling of intrinsic and radiation-enhanced fission gas (Xe) diffusion in UO2 +/- x: Implications for nuclear fuel performance modeling

Based on density functional theory (DFT) and empirical potential calculations, the diffusivity of fission gas atoms (Xe) in UO2 nuclear fuel has been calculated for a range of non-stoichiometry (i.e. UO2x), under both out-of-pile (no irradiation) and in-pile (irradiation) conditions. This was achieved by first deriving expressions for the activation energy that account for the type of trap site that the fission gas atoms occupy, which includes the corresponding type of mobile cluster, the charge state of these defects and the chemistry acting as boundary condition. In the next step DFT calculations were used to estimate migration barriers and internal energy contributions to the thermodynamic properties and calculations based on empirical potentials were used to estimate defect formation and migration entropies (i.e. pre-exponentials). The diffusivities calculated for out-of-pile conditions as function of the UO2x nonstoichiometrywere used to validate the accuracy of the diffusion models and the DFT calculations against available experimental data. The Xe diffusivity is predicted to depend strongly on the UO2x non-stoichiometry due to a combination of changes in the preferred Xe trap site and in the concentration of uranium vacancies enabling Xe diffusion, which is consistent with experiments. After establishing the validity of the modeling approach, it was used for studying Xe diffusion under in-pile conditions, for which experimental data is very scarce. The radiation-enhanced Xe diffusivity is compared to existing empirical models. Finally, the predicted fission gas diffusion rates were implemented in the BISON fuel performance code and fission gas release from a Risø fuel rod irradiation experiment was simulated. 2014 Elsevier B.V. All rights

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2019

Simulating Uranium(VI) Diffusion in Sodium-Montmorillonite as a Function of 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 in these systems. However, their prediction is complicated by the complex mineralogical structure of montmorillonite clay, leading to two types of surface sites and clay porosities, and by the dependence of U(VI) solution speciation on chemical 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, at the high degrees of clay compaction expected in future barrier systems, the ability of U(VI) species to access small, highly negatively charged interlayer spaces may ultimately control diffusive fluxes. In this modeling study, we investigate how the charge of U(VI) solution species affects U(VI) sorption and diffusion behavior. For this purpose, we simplified complex chemical solution speciation diagrams to charge speciation diagrams. Furthermore, we coupled an existing U(VI) sorption model (in PHREEQC) with U(VI) diffusion calculations (in GLE) to simulate U(VI) mobility in lab-scale through-diffusion experiments over various chemical conditions. Our results suggest that there are two important parameters that govern U(VI) diffusion in these systems. Uranium(VI) surface complexation on clay edge surfaces leads to contaminant retardation, while the access of U(VI) solution species to diffuse layers in clay interlayer spaces controls diffusive fluxes. Hence on the field scale, these parameters will determine how early a plume may cross site boundaries and pose risks to drinking water resources, and the levels of U(VI) concentrations to be expected in contaminated waters.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2016

Simulation of energy-dependent electron diffusion processes in the Earth's outer radiation belt

The radial and local diffusion processes induced by various plasma waves govern the highly energetic electron dynamics in the Earth's radiation belts, causing distinct characteristics in electron distributions at various energies. In this study, we present our simulation results of the energetic electron evolution during a geomagnetic storm using the University of California, Los Angeles 3-D diffusion code. Following the plasma sheet electron injections, the electrons at different energy bands detected by the Magnetic Electron Ion Spectrometer (MagEIS) and Relativistic Electron Proton Telescope (REPT) instruments on board the Van Allen Probes exhibit a rapid enhancement followed by a slow diffusive movement in differential energy fluxes, and the radial extent to which electrons can penetrate into depends on energy with closer penetration toward the Earth at lower energies than higher energies. We incorporate radial diffusion, local acceleration, and loss processes due to whistler mode wave observations to perform a 3-D diffusion simulation. Here, our simulation results demonstrate that chorus waves cause electron flux increase by more than 1 order of magnitude during the first 18 h, and the subsequent radial extents of the energetic electrons during the storm recovery phase are determined by the coupled radial diffusion and the pitch angle scattering by EMIC waves and plasmaspheric hiss. The radial diffusion caused by ULF waves and local plasma wave scattering are energy dependent, which lead to the observed electron flux variations with energy dependences. Lastly, this study suggests that plasma wave distributions in the inner magnetosphere are crucial for the energy-dependent intrusions of several hundred keV to several MeV electrons.

58 GEOSCIENCES↗
DOE OSTI2012

Modeling Intragranular Diffusion in Low-Connectivity Granular Media

Diffusive exchange of solutes between bulk water in an aquifer and water in the intragranular pores of the solid phase remains confusing after decades of study. In a previous paper, we reviewed some of the explanations, and suggested that the disparities between observation and theory were largely due to low connectivity of the intragranular pores. Low connectivity indicates that a useful conceptual framework is percolation theory, which guided our analysis. The present study was initiated to improve the finite difference (FD) model presented in the previous paper, and to test that new model rigorously against new random walk (RW) simulations of diffusion in low-connectivity porous spheres starting from non-equilibrium. The new FD model calculates diffusion separately in the infinite cluster and the finite clusters, and closely matches the new, more complex RW results. The percolation-theory based description of the new model is fairly simple, and can readily be incorporated into existing FD models. The simulations showed that the combination of low intragranular pore connectivity, and out-diffusion initiated at diffusive non-equilibrium, can produce diffusive behavior that appears as if the solute had undergone slow sorption, even in the absence of any sorption process. This mechanism may help explain some hitherto confusing aspects of intragranular diffusion.

diffusion↗
DOE OSTI2014

Diffusion length in nanoporous TiO{sub 2} films under above-band-gap illumination

We determined the carrier diffusion lengths in TiO{sub 2} nanoporous layers of dye-sensitized solar cells by using scanning photocurrent microscopy using an ultraviolet laser. Here, we excited the carrier directly in the nanoporous layers where the diffusion lengths were found to 140 μm as compared to that of visible illumination measured at 90 μm. The diffusion length decreased with increasing laser modulation frequency, in which we determined the electron lifetimes and the diffusion coefficients for both visible and UV illuminations. The diffusion lengths have been studied in terms of the sintering temperatures for both cells with and without binding molecules. We found a strong correlation between the diffusion length and the overall light-to-current conversion efficiency, proving that improving the diffusion length and hence the interparticle connections, is key to improving cell efficiency.

36 MATERIALS SCIENCE↗
DOE OSTI2011

Progress Report for Diffusion Welding of the NGNP Process Application Heat Exchangers (Revision 1)

The U.S. Department of Energy selected the high temperature gas-cooled reactor as the basis for the Next Generation Nuclear Plant (NGNP). The NGNP will demonstrate the use of nuclear power for electricity, hydrogen production, and process heat applications. The NGNP Project is currently investigating the use of metallic, diffusion welded, compact heat exchangers to transfer heat from the primary (reactor side) heat transport system to the secondary heat transport system. An intermediate heat exchanger will transfer this heat to downstream applications such as hydrogen production, process heat, and electricity generation. The channeled plates that make up the heat transfer surfaces of the intermediate heat exchanger will have to be assembled into an array by diffusion welding. This report describes the preliminary results of a scoping study that evaluated the diffusion welding process parameters and the resultant mechanical properties of diffusion welded joints using Alloy 800H. The long-term goal of the program is to progress towards demonstration of small heat exchanger unit cells fabricated with diffusion welds. Demonstration through mechanical testing of the unit cells will support American Society of Mechanical Engineers rules and standards development, reduce technical risk, and provide proof of concept for heat exchanger fabrication methods needed to deploy heat exchangers in several potential NGNP configurations.1 Researchers also evaluated the usefulness of modern thermodynamic and diffusion computational tools (Thermo-Calc and Dictra) in optimizing the parameters for diffusion welding of Alloy 800H. The modeling efforts suggested a temperature of 1150°C for 1 hour with an applied pressure of 5 MPa using 15 μm nickel foil as joint filler to reduce chromium oxidation on the welded surfaces. Good agreement between modeled and experimentally determined concentration gradients was achieved

08 HYDROGEN↗
DOE OSTI2021

Atomistic and mesoscale simulations to determine effective diffusion coefficient of fission products in SiC

The silicon carbide (SiC) layer in tristructural isotropic (TRISO) particles serves as the barrier to prevent escape of fission products produced in the fuel kernel. Knowing the diffusion coefficient of fission products through SiC is critical to determining whether fission gas can escape from the particle. It has been observed in experiments that Ag accumulated in grain boundaries and triple junctions in SiC. It is hypothesized that grain boundary diffusion is the primary pathway by which fission products penetrate the SiC layer. In this report, the effective diffusion coefficient of the fission product Ag through the grain boundary network is calculated using a combination of atomistic and phase-field methods. The grain boundary diffusion coefficient is calculated using molecular dynamics simulations. The bulk diffusion coefficient is determined using a combination of density functional theory and nudged elastic band methods. An effective diffusion coefficient is calculated, accounting for the grain structure using a phase-field method. The effective diffusion coefficient will be incorporated into Bison and fission product release calculations are compared to available experimental data.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2012

Comparison of homogenized and enhanced diffusion solutions of model PWR problems

Model problem comparisons in slab geometry are made between two forms of homogenized diffusion theory and enhanced diffusion theory. The pin-cell discontinuity factors for homogenized diffusion calculations are derived from homogenized variational nodal P1 response matrices and from standard finite differencing. Enhanced diffusion theory consists of applying quasi-reflected interface conditions to reduce variational nodal Pn response matrices to one degree of freedom per interface, without homogenization within the cell. As expected both homogenized diffusion methods preserve reaction rates exactly if the discontinuity factors are derived from the P 11 reference solutions. If no reference lattice solution is available, discontinuity factors may be approximated from single cells with reflected boundary conditions; the computational effort is then comparable to calculating the enhanced diffusion response matrices. In this situation enhanced diffusion theory gives the most accurate results and finite difference discontinuity factors the least accurate. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

The Sluggish Diffusion of Cations in CeO 2 Probed through Molecular Dynamics and Metadynamics Simulations

Cation diffusion in fluorite‐structured CeO 2 , though far slower than anion diffusion, is an important, high‐temperature process because it governs diverse fabrication and degradation phenomena. Herein, cation diffusion is studied by means of classical molecular dynamics and metadynamics simulations. Three different mechanisms are examined: migration involving an isolated cerium vacancy, migration involving a cerium vacancy in a defect associate with an oxygen vacancy, and migration involving a cation divacancy. For each mechanism, defect diffusion coefficients are calculated as a function of temperature, from which the respective activation enthalpy of defect migration is obtained. Through comparisons with experimental cation diffusion data (specifically, of the absolute magnitude of the cation diffusivity as well as its activation enthalpy), it is concluded that cation diffusion takes place predominantly neither by isolated vacancies nor by cation vacancy–oxygen vacancy associates but by cation divacancies.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Self-diffusion Measurements of Liquid Sn Using the Shear Cell Technique and Stable Density Layering

By utilizing the shear cell technique and achieving stable density layering with the addition of an alloying element Bi, the self-diffusion coefficients of liquid Sn were measured on the ground at 573 K, 773 K, and 973 K (300 °C, 500 °C, and 700 °C). Moreover, the impurity diffusion coefficients of Bi in the liquid Sn were simultaneously measured, to confirm the suppression of natural convection. From the experimental results, natural convection was confirmed to be suppressed, given that the impurity diffusion coefficients of Bi were in good agreement with the microgravity reference data. Upon changing the amounts of added Bi within 5 at. pct Bi, the self-diffusion coefficients of liquid Sn did not vary significantly. Moreover, the SnBi system within 5 at. pct Bi can be regarded as a dilute solution by calculating the activity of Sn in the SnBi system beforehand. The self-diffusion coefficients of the liquid Sn were in good agreement with the power law of temperature dependence from the microgravity reference data. After confirming the suppression of natural convection and that the SnBi system is a dilute solution, the self-diffusion coefficient of liquid Sn was determined as 5.00 ± 0.16 × 10{sup −9} m{sup 2} s{sup −1} at 773 K (500 °C). The effectiveness of the experimental method used in this study for the measurement of the self-diffusion coefficients was confirmed, under the condition that the liquid system could be regarded as a dilute solution.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Self-diffusion of Ti interstitial based point defects and complexes in TiC

To date, the mechanism of Ti atom self-diffusion is unproven. Prior theoretical work mostly focused on Ti vacancy based mediators, but these do not reproduce the experimental activation energy or entropy. In this work, in density functional theory calculations, Ti interstitials and related defect complexes are systematically considered as possible mediators of Ti self-diffusion. Among these defects, the defect complex of two C vacancies tightly bound to a Ti dumbbell is found to have the lowest formation energy. A sustainable migration of the complex, in a translational or rotational fashion, is enabled in the presence of another (free) carbon vacancy nearby the complex, and thus the rate of Ti self-diffusion by this mechanism is dependent on the concentration of carbon vacancies. The calculated activation energy of the complex agrees well with the experimental value in TiC 0.97 . Similar analyses of the Ti self-diffusion mechanisms mediated by Ti interstitials or dumbbells yield much higher activation energies, but the corresponding migration energies are evaluated to be less than 1 eV, which suggests they can be possible mediators of the radiation-enhanced Ti self-diffusion in TiC. To fully enable the comparison with experiments that are typically conducted at temperatures as high as 2500 K, we also consider the temperature dependent vibrational contribution to the activation energy of the defect complex. The vibrational contribution imposes an additive effect on the defect formation energy, while the migration energies are lowered due to the thermal expansion of the lattice. When combined, these factors give an excellent agreement with the experiments. Furthermore, this work gives strong support to the concept that Ti interstitial based defect complexes are likely diffusion mediators for Ti atom self-diffusion in TiC, further establishes a solid basis for large-scale modeling, and may eventually pave the way to accurately predicting defect-controlled diffusional processes.

36 MATERIALS SCIENCE↗
DOE OSTI2024

HTO and selenate diffusion through compacted Na-, Na–Ca-, and Ca-montmorillonite

Radionuclide transport in smectite clay barrier systems used for nuclear waste disposal is controlled by diffusion, with adsorption significantly retarding transport rates. While a relatively minor component of spent nuclear fuel, 79 Se is a major driver of the safety case for spent fuel disposal due to its long half-life (3.3×10 5 yr) and its low adsorption to clay (K D < 10 L/kg), thus a thorough understanding of Se diffusion through clay is critical for understanding the long-term safety of spent fuel disposal systems. Through-diffusion experiments with tritiated water (HTO, conservative tracer) and Se(VI) were conducted with a well-characterized, purified montmorillonite source clay (SWy-2) under a constant ionic strength (0.1 M) and three different electrolyte compositions: Na + , Ca 2+ , and a Na + -Ca 2+ mixture at pH 6.5 in order to probe the effects of electrolyte composition and interlayer cation composition on clay microstructure, Se(VI) aqueous speciation, and ultimately diffusion. Further, the results were modeled using a reactive transport modeling approach to determine values of porosity (ε), D e (effective diffusion coefficient), and K D (distribution coefficient for adsorption). HTO diffusive flux was higher in Ca-montmorillonite (D e =1.68×10 -10 m 2 s -1 ) compared to Na-montmorillonite (De=7.83×10 -11 m 2 s -1 ). This increase in flux is likely due to a greater degree of clay layer stacking in the presence of Ca 2+ compared to Na + , which leads to larger inter-particle pores. Overall, the Se(VI) flux was much lower than the HTO flux due to anion exclusion, with Se(VI) flux following the order Ca (D e = 1.03×10 -11 m 2 s -1 ) > Na–Ca (D e = 2.12×10 -12 m 2 s -1 ) > Na (D e = 1.28×10 -12 m 2 s -1 ). These differences in Se(VI) flux are due to a combination of factors, including (1) larger accessible porosity in Ca-montmorillonite due to clay layer stacking and smaller electrostatic effects compared to Na-montmorillonite, (2) larger accessible porosity for neutral-charge CaSeO4 species which makes up 32% of aqueous Se(VI) in the pure Ca system, and (3) possibly higher Se(VI) adsorption for Ca-montmorillonite. Through a combination of experimental and modeling work, this study highlights the compounding effects that electrolyte and counterion compositions can have on radionuclide transport through clay. Diffusion models that neglect these effects are not transferable from laboratory experimental conditions to in situ repository conditions.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗