Search NASA⌕ Search

SEARCH · Search NASA

Results for “diffusion processes”

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 1,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Multi-scale model based design of membrane reactor/separator processes for intensified hydrogen production through the water gas shift reaction

This paper aims to first quantify the impact of various diffusion models (Maxwell-Stefan, Wilke, Dusty-Gas) on the predictions of a multi-scale membrane reactor/separator mathematical model, and to then demonstrate this model's use for the design and process intensification of membrane reactor/separator systems for hydrogen production. This multi-scale model captures velocity, temperature and species' concentration profiles along the catalyst pellet's radial direction, and along the reactor's axial direction, by solving the momentum, energy, and species transport equations, accounting for convection, conduction, reaction, and diffusion mechanisms. In the first part of work, the effect of pellet-scale design parameters (mean pore diameter, volumetric porosity, tortuosity factor, etc.) and various species' flux models on the model predictions is studied. In the second part, the study focuses on the comparison, in terms of their process intensification characteristics, of various hydrogen production processes. These include a conventional high-temperature shift reactor (HTSR)/low-temperature shift reactor (LTSR) sequence, a novel HTSR/membrane separator (MS)/LTSR/MS sequence, and a process that involves low-temperature shift membrane reactors-LTSMR in a series.

08 HYDROGEN↗
DOE OSTI2020

Influence of Coating Method on the Performance of Roll-to-Roll Coated PEM Fuel Cell Catalyst Layers

As production of fuel cell vehicles and power systems increases there is a need to understand the process science of fuel cell manufacturing. Industrial production of polymer electrolyte membrane fuel cell (PEMFC) electrodes will rely on continuous roll-to-roll (R2R) coating methods to meet the rates needed for mass production. There are a variety of coating methods that could be used to coat catalyst layers. The different physics of each method may result in different morphologies and performance. Here we present the results of a study comparing fuel cell catalyst layers coated using two R2R coating methods: slot-die and gravure. These two methods were selected because the physics of applying the liquid to the substrate is very different between the two cases. Pt/C catalysts layers were coated directly onto carbon-fiber gas diffusion media to create gas diffusion electrodes (GDEs), without the use of a decal transfer process. These GDEs were assembled into membrane electrode assemblies (MEAs) by hot pressing to Nafion membranes. In situ performance and other electrochemical diagnostics were used to determine the influence of coating method on catalyst layer electrochemical properties. Electron microscopy was utilized to understand the influence of coating method on electrode morphology and its influence on electrochemical properties. We also explored the influence of solvent ratio to understand the potential coupling between ink formulation and coating methodology.

DIRECT ENERGY CONVERSION,MATERIALS SCIENCE↗
DOE OSTI2015

Evaluation of Technetium Getters to Improve the Performance of Cast Stone

Cast Stone has been selected as the preferred waste form for solidification of aqueous secondary liquid effluents from the Hanford Tank Waste Treatment and Immobilization Plant (WTP) process condensates and low-activity waste (LAW) melter off-gas caustic scrubber effluents. Cast Stone is also being evaluated as a supplemental immobilization technology to provide the necessary LAW treatment capacity to complete the Hanford tank waste cleanup mission in a timely and cost effective manner. One of the major radionuclides that Cast Stone has the potential to immobilize is technetium (Tc). The mechanism for immobilization is through the reduction of the highly mobile Tc(VII) species to the less mobile Tc(IV) species by the blast furnace slag (BFS) used in the Cast Stone formulation. Technetium immobilization through this method would be beneficial because Tc is one of the most difficult contaminants to address at the U.S. Department of Energy (DOE) Hanford Site due to its complex chemical behavior in tank waste, limited incorporation in mid- to high-temperature immobilization processes (vitrification, steam reformation, etc.), and high mobility in subsurface environments. In fact, the Tank Closure and Waste Management Environmental Impact Statement for the Hanford Site, Richland, Washington (TC&WM EIS) identifies technetium-99 ( 99 Tc) as one of the radioactive tank waste components contributing the most to the environmental impact associated with the cleanup of the Hanford Site. The TC&WM EIS, along with an earlier supplemental waste-form risk assessment, used a diffusion-limited release model to estimate the release of different contaminants from the WTP process waste forms. In both of these predictive modeling exercises, where effective diffusivities based on grout performance data available at the time, groundwater at the 100-m down-gradient well exceeded the allowable maximum permissible concentrations for 99 Tc. (900 pCi/L). Recent relatively short-term (63 day) leach tests conducted on both LAW and secondary waste Cast Stone monoliths indicated that 99 Tc diffusivities were at or near diffusivities where the groundwater at the 100-m down-gradient well would exceed the allowable maximum permissible 99 Tc concentrations. There is, therefore, a need and an opportunity to improve the retention of Tc in the Cast Stone waste form. One method to improve the performance of the Cast Stone waste form is through the addition of “getters” that selectively sequester Tc inside Cast Stone.

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

Phase field-volumetric lattice Boltzmann model of ion uptake in porous nuclear waste form materials under continuous flow

The flow field within the mesopores of sorbent particles plays a crucial role in radionuclide diffusion and ion uptake kinetics, thus, impacting the overall performance of porous nuclear waste form materials. To fundamentally understand the influence of microstructures and material properties on the radionuclide absorption and retention processes requires a coupled multi-physics model that considers the advection and diffusion within the flow field, the reaction at liquid-solid interfaces, and finally, the solid-state diffusion within a complex nanoporous medium. Here, this study employs the volumetric lattice Boltzmann method (VLBM) to accurately and efficiently calculate the steady state velocity field inside the mesopores of sorbent particles. The obtained velocity field is then utilized to calculate the advection of ions in the steady flow. A phase field (PF) model of ion uptake is used to describe the reaction occurring at the solid-liquid interface and diffusion inside the porous medium. The integrated PF-VLBM model is verified in terms of the mass conservation and numerical efficiency and validated qualitatively with experimental observation data. Then, it is applied to study the influence of thermodynamic and kinetic properties, as well as flow field conditions on the ion uptake kinetics. The numerical results demonstrate that the ion uptake kinetics in porous particles has three distinct stages, which is in agreement with the observations in continuous flow experiments. In the first stage, the kinetics is predominantly controlled by the flow field and ion diffusivity in the liquid phase. The kinetics in the second stage is primarily governed by ion diffusivity in the solid phase. In the third stage the system reaches a dynamic equilibrium with a net zero uptake flux at the interface. It is also found that porous structures significantly affect the efficiency and capacity of ion uptake. The simulation results can help to understand the physics behind the observed ion uptake kinetics in experiments and to facilitate the development of constitutive equations that can account for heterogeneous microstructures in engineering performance codes.

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

Energetic and dynamic analysis of transport of Na + and K + through a cyclic peptide nanotube in water and in lipid bilayers

Potential of mean force (PMF) profiles and position-dependent diffusion coefficients of Na + and K + are calculated to elucidate the translocation of ions through a cyclic peptide nanotube, composed of 8 × cyclo[-(D-Leu-Trp) 4- ] rings, in water and in hydrated DMPC bilayers. The PMF profiles and PMF decomposition analysis for the monovalent cations show that favorable interactions of the cations with the CPN as well as the lipid bilayer and dehydration free energy penalties are two major competing factors which determine the free energy surface for ion transport through CPNs both in water and lipid bilayers, and that the selectivity of CPNs to cations mainly arises from favorable interaction energies of cations with CPNs and lipid bilayers that are more dominant than the dehydration penalties. Calculations of the position-dependent diffusion coefficients and dynamic friction kernels of the cations indicate that the dehydration process along with the molecular rearrangements occurring outside the channel and the coupling of the ion motions with the chain-structured water movements inside the channel lead to decrease of the diffusion coefficients far away from the channel entrance and also reduced coefficients inside the channel. Here the PMF and diffusivity profiles for Na + and K + reveal that the energetics of ion transport through the CPN are governed by global interactions of ions with all the components in the system while the diffusivity of ions through the channel is mostly determined by local interactions of ions with the confined water molecules inside the channel. Comparison of Na + and K + ion distributions based on overdamped Brownian dynamics simulations based on the PMF and diffusivity profiles with the corresponding results from molecular dynamics shows good agreement, indicating accuracy of the Bayesian inference method for determining diffusion coefficients in this application. In addition this work shows that position-dependent diffusion coefficients of ions are required to explain the dynamics and conductance of ions through the CPN properly.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Computational Studies of Solid Stoppers for Generating Intense Exotic Radioactive Ion Beams at FRIB

Low-emittance beams of rare isotopes at low energies (less than 3 MeV/u) are required for many high precision measurements in nuclear science research. Such beams will be produced at the Facility for Rare Isotope Beams (FRIB) by first slowing down and thermalizing or stopping fast fragmentation beams, and subsequently reaccelerating them to the required energies for experiments. The use of a solid stopper has the potential to provide high-intensity reaccelerated beams of many species. We performed simulations of the processes critical for generating beams in solid stoppers, including beam heating, external heating, diffusive release, effusive transport, and surface sticking. Using a variety of software tools, we investigated the viability of numerous solid stopper system designs for FRIB. We present our computational approach, the physical processes we modeled, the codes we used, and our results. We conclude with a proof-of-concept solid stopper system that releases ions sufficiently quickly for experiments with unstable reaccelerated beams at FRIB.

47 OTHER INSTRUMENTATION↗
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 OSTI2025

Understanding Oxide–Metal Interactions During Hot Isostatic Pressing to Diffusion Bond Aluminum Alloy 6061 Plates

The interaction between Mg, Si, and Al 2 O 3 during hot isostatic pressing diffusion bonding of aluminum alloy 6061 (AA6061) plates was investigated through thermodynamic calculations and experimental microstructural characterization. Thermodynamic calculations as functions of temperature, pressure, and composition revealed that the interaction among Mg, Si, and Al 2 O 3 yields Mg 2 Si and either MgO + Al or MgAl 2 O 4 + Al, facilitating the reduction of Al 2 O 3 and allowing Al/Al metallic bonds to form. Total pressure variation had a negligible influence on the oxygen partial pressure, and consequently, the reaction product formation. Oxygen partial pressure variation as a function of temperature and initial amount of Al 2 O 3 determined the formation of either MgO or MgAl 2 O 4 . Experimental Hot Isostatic Pressure (HIP) bonding at 723 K and 833 K under a constant pressure of 1017 atm documented the cooling rate-dependent formation of β–Mg 2 Si precipitates. High-resolution transmission electron microscopy imaging and selected area electron diffraction patterns verified the formation of β–Mg 2 Si and MgO at the interface but did not detect MgAl 2 O 4 . In conclusion, findings from this study clarify the role of thermochemical interactions in oxide disruption and bonding mechanisms during HIP diffusion bonding of AA6061 and provide guidance for optimizing joining processes for monolithic nuclear fuel assemblies.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2015

Radiation defect dynamics in Si at room temperature studied by pulsed ion beams

The evolution of radiation defects after the thermalization of collision cascades often plays the dominant role in the formation of stable radiation disorder in crystalline solids of interest to electronics and nuclear materials applications. Here, we explore a pulsed-ion-beam method to study defect interaction dynamics in Si crystals bombarded at room temperature with 500 keV Ne, Ar, Kr, and Xe ions. The effective time constant of defect interaction is measured directly by studying the dependence of lattice disorder, monitored by ion channeling, on the passive part of the beam duty cycle. The effective defect diffusion length is revealed by the dependence of damage on the active part of the beam duty cycle. Results show that the defect relaxation behavior obeys a second order kinetic process for all the cases studied, with a time constant in the range of ∼4–13 ms and a diffusion length of ∼15–50 nm. Both radiation dynamics parameters (the time constant and diffusion length) are essentially independent of the maximum instantaneous dose rate, total ion dose, and dopant concentration within the ranges studied. However, both the time constant and diffusion length increase with increasing ion mass. This demonstrates that the density of collision cascades influences not only defect production and annealing efficiencies but also the defect interaction dynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Interfacial Void Formation and Self-Healing in Oxide Scales on Al-containing High-Entropy Alloy

The exceptional high-temperature oxidation resistance of Al-containing high-entropy alloys (HEAs) is often attributed to the formation of a protective α-Al 2 O 3 scale. However, the dynamic, atomic-scale mechanisms governing the stability of this scale—including interfacial void formation and the often-postulated but rarely visualized “self-healing” capacity—remain poorly understood. Herein, we reveal the complex evolution of the triple-layer oxide scale on an Al 10 CoCrFeNi HEA through combined electron microscopy and diffraction study. We show that interfacial voids are an inherent consequence of the scaling process, originating from two distinct mechanisms: the Kirkendall effect at the interface between the γ-Al 2 O 3 /α-Al 2 O 3 and alloy driven by cationic diffusion imbalance and volumetric contraction due to phase transformations at the spinel/Cr 2 O 3 interface. Crucially, we provide microstructural evidence consistent with an intrinsic self-healing response. This process is driven by coupled inward diffusion of oxygen and outward diffusion of metal cations, leading to the in-situ formation of transient θ-Al 2 O 3 and spinel phases that partially fill and seal the voids. Here, these results provide atomic-scale insights into the phase evolution, defect formation, and self-repair of oxide scales in HEAs—highlighting pathways to enhance their oxidation resistance in extreme environments.

High-entropy alloy↗
DOE OSTI2024

Dynamics of High Molecular Weight Cylindrical and Lamellar Block Copolymers with X‐ray Photon Correlation Spectroscopy

Here, the structure and dynamics of polystyrene (PS)-b-poly(ethylene oxide) block copolymers (BCPs) are studied. The BCPs exhibit microphase-separated cylindrical and lamellar morphologies. Structural dynamics are measured with X-ray photon correlation spectroscopy in the small-angle regime. Morphologies and domain sizes are evaluated using small-angle X-ray scattering (SAXS), scanning electron microscopy, and atomic force microscopy. Different solvent processing conditions are investigated. Grain sizes evaluated using SAXS are found to depend on processing only for the rubbery majority BCP. The structural relaxation times are examined as a function of PS volume fraction, temperature, morphology, and structural sizes. Well above the glass transition temperature (T g ) of PS, all samples exhibit stretched autocorrelation decays and diffusive dynamics. Near T g of PS, the dynamics of all samples are anomalous with compressed autocorrelation decays and hyperdiffusive dynamics. This transition occurs at 153 °C or 1.13 T g of PS. In the diffusive regime (at high temperature), structural relaxation times are dependent on the processing method. Near PS T g (at low temperature), structural relaxation times scale with the PS volume fraction. Structural relaxation times do not correlate with grain size, indicating that the out-of-equilibrium state of PS dominates the structural dynamics of these strongly phase-segregated BCPs.

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

Nature of Moisture-Induced fogging defects in scintillator plastic

Polyvinyl toluene (PVT) scintillator plastic may degrade in the field due to inward water diffusion at elevated temperatures, which exceeds the saturation limit at lower temperatures, and after long periods of time (e.g., many years), leads to the formation of disk-like defects that attenuate scintillation light, leading to detector degradation. In this work, using fractography and high-magnification optical and electron microscopy to characterize the water-induced defects, a model of the fogging process is hypothesized as follows: excess water present at low temperatures diffuses to spheroids to minimize contact with the hydrophobic polymer. Through the hydrophobic effect, the entropy of water increases by forming nanoclusters which minimizes the contact between the water and the PVT. As the water nanoclusters grow, they break and fold the polymer into densely-packed crystalline regions creating more space for water within the spheroids. The polymer outside the spheroid resists the shrinkage which builds up tension within the spheroid. Once the tensile stress exceeds the yield strength of the plastic, the spheroid is torn in half resulting in a defect. Excess water then drains into the cavities along the disk thereby further increasing its entropy. Slower cooling (over 1 day) leads to larger spheroids and hence, larger “permanent” defects. Freezing causes some defects to further grow due to the expansion of water to ice. These findings imply that the remaining lifetime of scintillator plastic in the field could be predicted using temperature and humidity data thereby mitigating security risks of degraded radiation portal monitors.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Ni-base superalloy powder-processed porous layer for gas cooling in extreme environments

Extreme high temperature conditions demand novel solutions for hot gas filters and coolant access architectures, i.e., porous layers on exposed components. These high temperatures, for example in current turbine engines, are at or exceeding current material limits for high temperature oxidation/corrosion, creep resistance, and, even, melting temperature. Thus novel blade designs allowing greater heat removal are required to maintain airfoil temperatures below melting and/ or rapid creep deformation limits. Gas atomized Ni-base superalloy powders were partially sintered into porous layers to allow full-surface, transpirational cooling of the surface of airfoils. Furthermore, these powder-processed porous layers were fully characterized for surface, morphology, cross-sectional microstructure, and mechanical strength characteristics. A sintering model based on pure Ni surface diffusion correlated well with the experimental results and allowed reasonable control over the partial sintering process to obtain a specified level of porosity within the porous layer.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Interfacial reaction in SiC{sub f}/Ti-6Al-4V composite by using transmission electron microscopy

The interfacial reactions of continuous SiC fiber reinforced Ti-6Al-4V matrix composite (SiC{sub f}/Ti-6Al-4V composite) and continuous SiC fiber coated by C reinforced Ti-6Al-4V matrix composite (SiC{sub f}/C/Ti-6Al-4V composite) were investigated by using micro-beam electron diffraction (MBED) and energy disperse spectroscopy (EDS) on transmission electron microscopy (TEM). The sequence of the interfacial reactions in the as-processed and exposed at 900°C for 50h SiC{sub f}/Ti-6Al-4V composites can be described as SiC||TiC||Ti{sub 5}Si{sub 3} + TiC||Ti-6Al-4V and SiC||TiC||Ti{sub 5}Si{sub 3}||TiC||Ti{sub 5}Si{sub 3}||TiC||Ti{sub 5}Si{sub 3}||Ti-6Al-4V, respectively. Additionally, both in as-processed and exposed composites, Ti{sub 3}SiC{sub 2} and Ti{sub 3}Si are absent at the interfaces. For the SiC{sub f}/C/Ti-6Al-4V composite exposed at 900 °C for 50 h, the sequence of the interfacial reaction can be described as SiC||C||TiC{sub F}||TiC{sub C}||Ti-6Al-4V before C coating is completely consumed by interfacial reaction. When interfacial reaction consumes C coating completely, the sequence of the interfacial reaction can be described as SiC||TiC||Ti{sub 5}Si{sub 3}||TiC||Ti-6Al-4V. Furthermore, in SiC{sub f}/C/Ti-6Al-4V composite, C coating can absolutely prevent Si diffusion from SiC fiber to matrix. Basing on these results, the model of formation process of the interfacial reaction products in the composites was proposed. - Highlights: • We obtained the sequence of the interfacial reactions in the as-processed and exposed at 900 °C for 50 h SiC{sub f}/Ti-6Al-4 V composites as well as in the SiC{sub f}/C/Ti-6Al-4 V composite exposed at 900 °C for 50 h. • We verified that both in as-processed and exposed SiC{sub f}/Ti-6Al-4 V composites, Ti{sub 3}SiC{sub 2} and Ti{sub 3}Si are absent at the interfaces. • Carbon coating can absolutely prevent silicon diffusion from SiC fiber to matrix. • Basing on these results, the model of formation process of the interfacial reaction products in the composites was proposed.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Bio‐Imaging quorum sensing signal molecules in a soil‐mimic gel

Rhizosphere, the narrow region surrounding plant roots directly influenced by root exudates and the root associated microbiome, plays an important role in plant productivity and the rhizosphere is well known for stimulating microbial metabolic activities. How microbial communities interact to form stable, metabolically interconnected functional communities is an area of intense interest. The question remains on how microbially produced secreted molecules that function as intercellular communication signals shape the structure and function of microbial communities. Our goal is to detect and quantify signal molecules produced by microbes or plants in the root‐soil environment, spatially and temporarily. As a proof‐of‐concept, we are imaging diffusible extracellular microbial metabolites involved in a bacteria‐bacteria communication process called quorum‐sensing. Quorum‐sensing relies on the accumulation of high concentrations of signal molecules in the environment to control bacterial gene expression, influencing rhizosphere colonization and plant health. Local concentration of quorum sensing molecules are determined using aptamer based sensors. Aptamers that can specifically bind the desired signal molecules are selected through SELEX and immobilized on the surface of nano‐porous membranes. Binding of the signal molecules and aptamer covered surface results in changes in surface charge distribution and steric hinderance and thus modifying the transmembrane ionic transport. Changes in transmembrane impedance can be measured through electrochemical impedance spectroscopy methods to monitor local concentrations of signal molecules. Sensor responses were determined for different concentrations of signal molecules and results showed detection of C4‐HSL in a soil‐mimic solution with K D of 10 nM. We demonstrate that the quorum‐sensing signal molecule C4‐homoserine lactone and (C4‐HSL) can rapidly diffuse in a soil‐mimic gel, providing evidence for our use of a soil‐mimic for further aptasensor development and validation. The sensors were then inserted into soil mimic gel for monitoring C4‐HSL diffusion and the resulted impedance changes were used to determine the C4‐HSL concentration at different positions in the gel. Measurements of local C4‐HSL concentrations variation and numerical solution of diffusion equation were used to create 4D images of C4‐HSL molecule diffusion in the soil‐mimic gel.

Jiang, Nianyu↗
DOE OSTI2023

Impact of oxidation and carbonation on the release rates of iodine, selenium, technetium, and nitrogen from a cementitious waste form

Evaluation of the long-term retention mechanisms and potential release rates for the primary constituents of potential concern (COPCs) (i.e., Tc, I, Se, and nitrate) is necessary to determine if Cast Stone, a radioactive waste form, can meet performance objectives under near-surface disposal scenarios. In this work, a mineral and parameter set accounting for the solubility of I and Se in Cast Stone was developed based on pH-dependent and monolithic diffusion leaching test results, to extend a geochemical speciation model previously developed. The impact of oxidation and carbonation as environmental aging processes on the retention properties of Cast Stone for primary COPCs was systematically estimated. Physically, the effective diffusion coefficients of 4 COPCs in Cast Stone were increased after carbonation and/or oxidation, reflecting an increase in permeability to diffusion. Chemically, i) pH & pe conditions in the original Cast Stone were favorable for the stabilization of Tc, but not for I, Se, and N; ii) oxidation (with/without carbonation) of Cast Stone changed the pe & pH conditions to be detrimental for Tc stabilization; and iii) carbonation (with/without oxidation) of Cast Stone modified the pH & pe conditions to be beneficial for the stabilization of I (in system with Ag added) and Se.

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

Final Report: Dominant Mechanisms of Uranium-Phosphate Reactions in Subsurface Sediments

Phosphate addition is an in situ remediation approach that may enhance the sequestration of uranium without requiring sustained reducing conditions. However, the geochemical factors that determine the dominant immobilization mechanisms upon phosphate addition are insufficiently understood to design efficient remediation strategies or accurately predict U(VI) transport. The overall objective of our project is to determine the dominant mechanisms of U(VI)-phosphate reactions in subsurface environments. Our research approach seeks to determine the U(VI)-phosphate solid that form in the presence of different groundwater cations, characterize the effects of phosphate on U(VI) adsorption and precipitation on smectite and iron oxide minerals, examples of two major reactive mineral phases in contaminated sediments, and investigate how phosphate affects U(VI) speciation and fate during water flow through sediments from contaminated sites. The research activities conducted for this project have generated a series of major findings. U(VI) phosphate solids from the autunite mineral family are the sole phases to form during precipitation, with uranyl orthophosphate not occurring despite its predicted greater stability. Calcium phosphates may take up substantial quantities of U(VI) through three different removal processes (adsorption, coprecipitation, and precipitation) but the dominance of each process varies with the pathway of reaction. Phosphate co-adsorbs with U(VI) onto smectite mineral surfaces, forming a mixed uranium-phosphate surface complex over a wide range of conditions. However, this molecular-scale association of uranium and phosphate has not effect on the overall extent of uptake. In contrast, phosphate enhanced U(VI) adsorption to iron oxide minerals at acidic pH conditions but suppresses such adsorption at neutral and alkaline pH, despite forming mixed uranium-phosphate surface complexes during adsorption. Nucleation barriers exist that inhibit U(VI) phosphate solids from precipitating in the presence of smectite and iron oxide minerals as well as sediments from contaminated sites. Phosphate addition enhances retention of U(VI) by sediments from the Rifle, CO and Hanford, WA field research sites, areas containing substantial uranium contamination of groundwater. This enhanced retention is through adsorption processes. Both fast and slow uptake and release behavior is observed, indicating that diffusion of uranium between sediment grains has a substantial effect of U(VI) fate in flowing groundwater systems. This project has revealed the complexity of U(VI)-phosphate reactions in subsurface systems. Distinct chemical processes occur in acidic and alkaline groundwater systems. For the latter, calcium phosphate formation, solution complexation, and competition between phosphate and uranium for adsorption sites may serve to either enhance or inhibit U(VI) removal from groundwater. Under the groundwater conditions present at many contaminated sites in the U.S., phosphate appears to general enhance U(VI) retention and limit transport. However, formation of low-solubility uranium phosphate solids does not occur under field-relevant conditions, despite this being the desired product of phosphate-based remediation approaches. In addition, simple equilibrium approaches fail to well-predict uranium fate in contaminated sediments amended with phosphate, with reactive transport models that include reaction rates and mass transport through occluded domains needed to properly describe the system. Phosphate addition faces challenges to being effective as a stand-alone groundwater treatment approach but would prove beneficial as an add-on to other treatment methods that will further limit uranium migration in the subsurface.

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

Predicting Correlation Coefficients for Monte Carlo Eigenvalue Simulations

Monte Carlo methods are most often considered as a reference for neutron transport simulations since very limited approximations are made about nuclear data and system geometry. To report uncertainty of any tally evaluated as generation averages, the sample variance is divided by the number of active generations, which is based on the assumption that the neutron generations are independent. Correlation effects between neutrons in multiplying systems, particularly when performing power iteration to evaluate eigenvalues have been observed in previous work. Neglecting the correlation effect results in an underestimate of uncertainty reported by Monte Carlo calculations. Previous work has also proposed methods to predict the underestimation ratio. Yamamoto et al expanded the fission source distribution with diffusion equation modes, performed numerical simulation of the AR (autoregressive) process of the expansion coefficients and used the correlation of the AR process to predict that of the Monte Carlo eigenvalue simulation. Sutton applied the discretized phase space (DPS) approach to predict the underestimation ratio but the method cannot predict the ratio when one neutron generates offspring in different phase space regions or generates a random number of offspring. This paper presents a method to predict the correlation effect with the model of multi-type branching processes (MBP). The method requires simulations for one generation of neutrons without knowing the source distribution and can predict the underestimation ratio for the cases where the traditional DPS approach does not work. The generation-to-generation correlation determines the convergence rate of active generations, the bias of variance estimator for each generation and the underestimation ratio of variance estimator for tallies averaged over active generations. The generation-to-generation correlation is characterized by the Auto-Correlation Coefficients (ACC) between tallies from different generations. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗