Search NASA⌕ Search

SEARCH · Search NASA

Results for “Contaminant 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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Defining a Metric to Assess Leachability of Contaminants from Cementitious Waste Forms in Lab Testing Related to Hanford - 20478

This paper reports on the determination of performance metrics for the release of four contaminants (technetium, iodine, nitrate, and chromium) from a theoretical supplemental low-activity waste (LAW) inventory of grout in the Hanford Integrated Disposal Facility (IDF). The performance metric was defined as a retardation factor, R, that can be used to quantify a waste form's ability to retain one contaminant compared with an un-retarded species (e.g., sodium). By setting the target R value at the magnitude required to maintain concentrations of the contaminants below the drinking water standard at the point of compliance for the IDF, a bar for grout leach performance can be established for screening laboratory leach test data. The R can be converted to a corresponding laboratory-observed diffusivity value against which to assess leach test data. The individual performance metrics were determined through an iterative set of simulations using a variable saturation and transport flow model of a grout waste form inventory in the IDF to calculate the resulting concentration-based requirements at the point of compliance. (authors)

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

Mixing Cell Model: A One-Dimensional Numerical Model for Assessment of Water Flow and Contaminant Transport in the Unsaturated Zone

This report describes the Mixing Cell Model code, a one dimensional model for water flow and solute transport in the unsaturated zone under steady state or transient flow conditions. The model is based on the principles and assumptions underlying mixing-cell model formulations. The unsaturated zone is discretized into a series of independent mixing cells. Each cell may have unique hydrologic, lithologic, and sorptive properties. Ordinary differential equations describe the material (water and solute) balance within each cell. Water-flow equations are derived from the continuity equation, assuming that unit gradient conditions exist at all times in each cell. Pressure gradients are considered implicitly through model discretization. Unsaturated hydraulic conductivity and moisture contents are determined by the material specific moisture-characteristic curves. Solute-transport processes include explicit treatment of advective processes, first order chain decay, and linear sorption reactions. Dispersion is addressed through implicit and explicit dispersion. Implicit dispersion is an inherent feature of all mixing-cell models and originates from the formulation of the problem in terms of mass balance around fully mixed volume elements. Expressions are provided that relate implicit dispersion to the physical dispersion of the system. Two FORTRAN codes were developed to solve the water flow and solute-transport equations: (1) the Mixing Cell Model for Flow (MCMF) solves transient water-flow problems and (2) the Mixing Cell Model for Transport (MCMT) solves the solute-transport problem. The transient water-flow problem is typically solved first by estimating the water flux through each cell in the model domain as a function of time using the MCMF code. These data are stored in either ASCII or binary files that are later read by the solute transport code (MCMT). Code output includes solute pore water concentrations, water and solute inventories in each cell and at each specified output time, and water and solute fluxes through each cell and specified output time. Computer run times for coupled transient water flow and solute transport were typically several seconds on a 2 GHz Intel Pentium IV desktop computer. The model was benchmarked against analytical solutions and finite element approximations to the partial differential equations (PDE) describing unsaturated flow and transport. Differences between the maximum solute flux estimated by the mixing cell model and the PDE models were typically less than two percent. This revision includes an option for a fixed concentration lower boundary condition for diffusive fluxes for versions 020321 and later.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Sorption kinetics and stability of conventional adsorbents for mercury remediation

In-situ remediation of mercury at numerous contaminated sites worldwide is a challenging and costly endeavor due to the persistency of this contaminant. In this study, we evaluated eight commercially available sorbent media ranging from carbon-, clays- and silica-based materials (PBC– Biochar, eSorb – Sorbster, nsPAC – Powdered Activated Carbon, fsPAC – Powdered Activated Carbon with Mackinawite, F300 – Filtrasorb 300, Si-SH – Silica Thiol, eBind – RemBind, Q-Clay – Organoclay PM-199), for their effectiveness in sorbing mercury (Hg 2+ ) and mercury complexed with dissolved organic matter (Hg-DOM). Under the chosen experimental conditions of this study, results showed that in the absence of DOM, the kinetic rates of Hg 2+ sorption onto the evaluated sorbents were in the order of 0.31 min-1 (Si-SH) to 2.98 min -1 (nsPAC), whereas in the presence of DOM, the rates varied from 0.16 min-1 (F300) to 0.95 min -1 (nsPAC). The measured sorption capacity for Hg 2+ in the absence of DOM varied from 3.02 mg/g (Q-Clay) to 35.15 mg/g (Si-SH), whereas in the presence of DOM, calculated partition coefficient (KD) ranged from 69.7 mL/g (Q-Clay) to 41,510 mL/g (Si-SH). Furthermore, kinetic data suggest liquid film diffusion was the rate-limiting steps governing mercury sorption onto the studied media. Overall, the obtained study parameters (kinetics/isotherm) are particularly important in informing robust engineering designs for deployment of the vast majority of evaluated sorbents. Thus, sorbent-based strategies offer viable solutions for cost-effective cleanup of mercury at industrially contaminated sites.

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

A comparative study of machine learning models for predicting the state of reactive mixing

Mixing phenomena are important mechanisms controlling flow, species transport, and reaction processes in fluids and porous media. Accurate predictions of reactive mixing are critical for many Earth and environmental science problems such as contaminant fate and remediation, macroalgae growth, and plankton biomass evolution. Here, to investigate the evolution of mixing dynamics under different scenarios (e.g., anisotropy, fluctuating velocity fields), a finite-element-based numerical model was built to solve the fast, irreversible bimolecular reaction-diffusion equations to simulate a range of reactive-mixing scenarios. A total of 2,315 simulations were performed using different sets of model input parameters comprising various spatial scales of vortex structures in the velocity field, time-scales associated with velocity oscillations, the perturbation parameter for the vortex-based velocity, anisotropic dispersion contrast (i.e., ratio of longitudinal-to-transverse dispersion), and molecular diffusion. The outputs comprised concentration profiles of reactants and products. The inputs to and outputs from these simulations were concatenated into feature and label matrices, respectively, to train 20 different machine learning (ML) models intended to emulate system behavior. These 20 ML emulators, based on linear methods, Bayesian methods, ensemble learning methods, and multilayer perceptrons (MLPs), were trained to classify the state of mixing and predict three quantities of interest (QoIs) characterizing species production, decay (i.e., average concentration, square of average concentration), and degree of mixing (i.e., variances of species concentration). Unsurprisingly, linear classifiers and regressors failed to reproduce the QoIs; however, ensemble methods (classifiers and regressors) and the MLP model accurately classified the state of reactive mixing and the QoIs. Among ensemble methods, random forest and decision-tree-based AdaBoost faithfully predicted the QoIs. At run time, trained ML emulators produced results times faster than the finite-element simulations. Due to their low computational expense and high accuracy, ensemble and MLP models are excellent emulators for these numerical simulations and great utilities in uncertainty quantification exercises, which can require 1,000s of forward model runs.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

User’s Manual for RESRAD-BUILD Code V.4: Vol. 1 – Methodology and Models Used in RESRAD-BUILD Code

The RESRAD-BUILD computer code models radionuclide release and transport in indoor environments and performs pathway analyses to evaluate the potential radiological dose and risk incurred by an individual who works or lives in a building contaminated with radioactive material or housing radioactively contaminated furniture or equipment. The code provides four geometries to characterize a radiation source: point, line, area, and volume, in which radionuclides are homogeneously distributed. Radionuclides contained in a source are considered to be released to the indoor air due to various processes including erosion (mechanically or weathering), diffusion (for tritium and radon in a volume source), or emanation (radon in a point, line, or area source). The release can proceed through different time phases with different rates. In RESRAD-BUILD Version 4.0, a dynamic ventilation model is implemented to simulate the fate and transport of source material particles and radionuclides after their releases. This dynamic ventilation model considers (1) air exchange between rooms in the building and between the rooms and the outdoor environment, (2) deposition from air to floor, (3) resuspension from the floor to the air, and (4) periodical vacuuming that reduces the floor deposition. The fate and transport modeling provides estimates of radionuclide concentrations in the source, in the air, and on the floor at different times, which are then integrated over the exposure duration for the calculation of radiation doses and cancer risks. A single run of the RESRAD-BUILD code can model a building with up to 9 rooms, 10 sources, and 10 receptors. The potential radiation dose and cancer risk incurred by each receptor are calculated for seven exposure pathways: (1) external radiation directly from the sources (accounting for shielding), (2) external radiation from radioactive particles deposited on the floors, (3) external radiation from airborne radionuclides, (4) inhalation of airborne radionuclides, (5) inhalation of radon and radon progenies, (6) inadvertent ingestion of radioactive particles directly from the source, and (7) ingestion of radioactive particles deposited on the floors. Various exposure scenarios can be modeled with RESRAD-BUILD, including but are not limited to, office worker, renovation worker, decontamination worker, building visitor, and resident. Both deterministic and probabilistic analyses can be performed to obtain results in both text reports and graphic displays.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2025

A Systematic Search for MeV–GeV Pulsar Wind Nebulae without Gamma-Ray Detected Pulsars

An increasing number of pulsar wind nebulae (PWNe) are being identified in the TeV band by ground-based Imaging Air Cherenkov Telescopes such that they constitute the dominant source class of Galactic TeV emitters. However, MeV–GeV PWN counterparts are still largely lacking. To date, only a dozen PWNe are identified by the Fermi–Large Area Telescope (LAT) in the MeV–GeV band. Most PWNe are located along the Galactic plane embedded within the prominent, diffuse Galactic γ-ray emission, which makes these sources difficult to disentangle from the bright diffuse background. We present a systematic search for γ-ray counterparts to known PWNe in the 300 MeV–2 TeV energy band using the Fermi–LAT. We target the locations of previously identified PWNe that lack detected Fermi–LAT pulsars to minimize associated pulsar contamination. The sample includes six previously identified Fermi PWNe and eight Fermi–LAT sources associated with PWNe. We report the analysis of 58 regions of interest and classify Fermi–LAT detected sources as either a likely PWN or a candidate PWN counterpart based on their morphological and spectral characteristics across the broadband spectrum. There are nine unidentified Fermi–LAT sources that we consider as likely PWN counterparts, which, if confirmed to be PWNe, would greatly increase the PWN population detected by the Fermi–LAT from 12 to 21. The remaining Fermi–LAT detected sources are considered weaker PWN candidates. A second approach in the systematic search for γ-ray emitting PWNe will involve studying the off-pulse phases of Fermi–LAT detected pulsars for the presence of an obscured PWN and will be reported in a subsequent paper.

Acharyya, A. [Univ. of Southern Denmark, Odense (D↗
DOE OSTI2021

ForSE: A GAN-based Algorithm for Extending CMB Foreground Models to Subdegree Angular Scales

We present ForSE (Foreground Scale Extender), a novel Python package that aims to overcome the current limitations in the simulation of diffuse Galactic radiation, in the context of cosmic microwave background (CMB) experiments. ForSE exploits the ability of generative adversarial neural networks (GANs) to learn and reproduce complex features present in a set of images, with the goal of simulating realistic and non-Gaussian foreground radiation at subdegree angular scales. This is of great importance in order to estimate the foreground contamination to lensing reconstruction, delensing, and primordial B-modes for future CMB experiments. We applied this algorithm to Galactic thermal dust emission in both total intensity and polarization. Our results show how ForSE is able to generate small-scale features (at 12') having as input the large-scale ones (80'). The injected structures have statistical properties, evaluated by means of the Minkowski functionals, in good agreement with those of the real sky and which show the correct amplitude scaling as a function of the angular dimension. Furthermore, the obtained thermal dust Stokes Q and U full-sky maps as well as the ForSE package are publicly available for download.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Nature of contaminants introduced in silicon wafers during molecular beam epitaxy chamber annealing

Epitaxial monolithic III–V/Si tandem solar cells are one of the most promising technologies to be adopted by the industry after the efficiency of the current market dominating single junction silicon solar cell saturates at its fundamental limit. One of the key limitations of this technology is the degradation of silicon wafers during in situ annealing in the molecular beam epitaxy chamber. Determining the nature of contaminants in this process is key to improve the efficiency of epitaxial tandem solar cells. However, to date, the nature of contaminants from molecular beam epitaxy chambers remains unknown. In this work, we use photoluminescence imaging, lifetime spectroscopy, and deep level transient spectroscopy to measure the electronic properties of extrinsic impurities incorporated during annealing in the molecular beam epitaxy chamber. Photoluminescence images reveal that at least two impurities diffuse into silicon wafers during molecular beam epitaxy annealing. One is highly localized, while the other one is distributed uniformly across the whole wafer. Phosphorus diffusion is found to confine the localized impurity within the diffused layer but is ineffective at preventing the indiffusion of other impurities. Lifetime spectroscopy shows that metastable impurities with characteristic similar to Cr and CrB in our molecular beam epitaxy annealed silicon wafers. No evidence of Fe or FeB was observed. The emission rates and concentrations of the electrically active defects were measured with deep-level transient spectroscopy: The emission rates of detected defects do not match that of known Cr-related defects.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A potential emerging issue concerning repair welding of out-of-core PWR components involving tritium exposure and 3 He retention

This paper identifies a potential but previously unrecognized risk of helium-induced embrittlement and cracking during repair welding of out-of-core pressurized water reactor (PWR) components exposed to tritium-contaminated coolant. While previous weldability concerns centered on 4 He accumulation in neutron-irradiated alloys located within the in-core or near-core regions, new measurements show that 3 He generated by tritium decay can accumulate in out-of-core components. Because hydrogen isotopes readily diffuse to grain boundaries and become trapped there, significant 3 He generation at grain boundaries may lead to cracking during weld repairs. Initial data from far-below-core PWR flux thimble tubes confirm the presence of 3 He levels above known cracking thresholds for repair welds. In conclusion, these findings indicate that out-of-core regions should be considered when defining safe weld repair windows in reactors operating for 60–100 years.

Extended PWR lifetimes↗
DOE OSTI2020

The Influence of Aqueous Dispersivity on Transport of Radiological Material at LANL MDA G - 20450

Aqueous dispersivity is a phenomenon by which a solute moving in groundwater mixes with water that does not contain the solute and becomes diluted. Solutes take paths of varying lengths through pore networks or experience slower velocities due to interactions with pore walls. Aqueous dispersivity leads to breakthrough curves with earlier arrival than expected in purely advective systems and persistent 'tails,' or late-arriving, dilute concentrations of the solute. Dispersivity can be longitudinal (in the direction of flow) or transverse (orthogonal to flow). LANL has disposed of radioactive material at Material Disposal Area (MDA) G since the late 1950's. Most of the disposed materials at MDA G are low-level waste (LLW), but other types of radioactive wastes are disposed or stored at the site. For example, there is a collection of remote-handled transuranic waste (RH-TRU) waste in 33 shafts, which were intended to be stored retrievably. Aqueous-phase radionuclides included in models of contaminant transport at MDA G are subject to dispersive processes in porous media of the vadose zone and saturated zone. To model dispersivity at MDA G probabilistically, we created site-specific distributions of aqueous dispersivity. First, we performed a comprehensive review of peer-reviewed publications, reports from other performance assessments of MDA G, and documentation from sites that were considered similar based on expert opinion, such as Yucca Mountain. Data were included from a wide range of sources, so that distributions developed based on the data reflect the current state of knowledge in properties of the subsurface environment. Data were grouped based on value type (e.g., determined via experiment or used in a different model), rock type, zone type (saturated or unsaturated), and site, and were weighted in the distributions according to their relevance to the site and zone type. We combined Monte Carlo simulation and bootstrapping techniques to produce approximate distributions of mean aqueous dispersivity. Distributions of dispersivity for saturated and unsaturated zone transport are differentiated by data weight and the length scales over which transport occurs. We find that underestimating or disregarding transverse dispersivity may lead to overstating simulated well capture in the saturated zone at MDA G. Neglecting transverse dispersivity is a common assumption due in part to use of one-dimensional transport models that do not consider solute spreading orthogonal to flow. Also, we find that variability in longitudinal dispersivity is associated with substantial travel time alterations in the MDA G unsaturated zone and in simulated contaminant breakthrough times at a downstream observation plane in the MDA G saturated zone. Finally, at large scales (>1000 m), a linear trend between the dispersivity coefficient and length over which transport occurs may overestimate the dispersivity coefficient. Future modeling work with application to solute spreading could include interacting advective and diffusive zones in the highly fractured basalt below MDA G. (authors)

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

Point defects and impurities in fluorite PuO 2

The native surface oxide of plutonium plays a critical role in ensuring the stability and safe storage of the underlying metal; consequently, understanding the role of defects and impurities in determining the properties of the oxide layer is critical. Here, in this study, we use hybrid density-functional theory calculations to evaluate the electronic structure and defect chemistry of PuO 2 , the most stable of the native oxide phases, including both native and extrinsic defects. We find that oxygen vacancies (𝑉 O ) form readily in PuO 2 , as do polarons. Electron polarons (𝜂 − ) are the lowest-energy acceptor species in PuO 2 , while the charge compensating donor species will shift from 𝑉 O under O-poor conditions to hole polarons (𝜂 + ) under O-rich conditions. Nitrogen and fluorine can substitute readily for oxygen atoms under O-poor conditions, while fluorine can also incorporate in an interstitial configuration (F$^−_i$) under more O-rich conditions. Carbon and chlorine incorporation in PuO 2 will be very limited. We also evaluate the kinetic barriers for oxygen-related defects, which we find to diffuse readily when present. Our results provide valuable insights into the critical role and variable chemistry of point defects and impurities in PuO 2 , which in turn have important implications for the safe storage of the underlying metal layer. In short, exposure of freshly prepared plutonium to reactive nitrogen- and fluorine-containing contaminants should be avoided, while carbon- or chlorine-containing contaminants are less likely to incorporate readily into the oxide.

Materials science↗
DOE OSTI2019

Contaminant Cation Effect on Oxygen Transport through the Ionomers of Polymer Electrolyte Membrane Fuel Cells

We characterized the effects of cobalt (Co 2+ ) and other cation contaminants on the oxygen (O 2 ) transport properties of the PFSA ionomer used in polymer electrolyte membrane fuel cells (PEMFCs) and gained insight into the mechanisms by which contaminant cations inhibit O 2 transport. Such cations can be released by alloy catalysts and environmental conditions and pose a significant challenge to maintaining high current density performance with low platinum (Pt) loadings. We used a test cell capable of isolating the ionomer from the membrane electrode assembly (MEA), allowing for O 2 transport resistance (R O2 ) measurements using a limiting current technique. We contaminated ionomer membranes with Li + , Na + , Ni 2+ , Co 2+ , and Ce 3+ and found a general increase in R O2 for increased contamination levels and decreased water activity. In addition, our Co 2+ results indicated distinct concentration-dependent regimes. The other cation-form ionomers allowed us to separate the impacts of ion pair strength, multivalency, and reduced water uptake. We believe that these factors cause a more compressed hydrophilic domain and tortuous O 2 diffusion path, and a commensurate increase in R O2 . Lastly, we studied the impact of Co 2+ on an operating PEMFC and found an increase in R O2 consistent with the results of our isolated membrane tests.

25 ENERGY STORAGE↗
DOE OSTI2023

Evaluation of Zero-Net-Rate Pumping Tests

Accurately estimating the distribution of aquifer properties is key to understanding contaminant movement in the subsurface. The distribution of aquifer properties is typically addressed using slug or constant-rate well tests, and the pros and cons of these tests are well known. Slug tests are appealing because they avoid removing contaminated water, but their results are affected by well skin and the small volume of displaced water limits the volume of aquifer that can be evaluated. Constant-rate well tests have the disadvantage of requiring disposal of potentially contaminated water, but they can generate properties that are more representative than slug tests, and they can be used to estimate well efficiency and storativity, which are difficult to characterize using slug tests. Periodic pumping tests are appealing because they have many of the advantages and few of the disadvantages of slug and constant-rate well tests. Periodic pumping tests involve cycling the pumping rate with a regular period and measuring the resulting response in monitoring wells. Some of these tests involve imposing a periodic rate on a constant mean rate. However, other tests involve moving water out and back into the well at rates that are balanced so the net rate after each period is zero. Some zero-net-rate (ZNR) well tests use rates that follow a sinusoidal pattern, whereas others use a square wave pattern that switches between a constant rate of pumping and a constant rate of injection to achieve a zero-net rate. ZNR well tests appear to be a useful compromise between slug and constant rate tests, but methods for conducting and analyzing the results from these tests have received limited evaluation. Periodic pumping tests are typically evaluated through recorded pressure responses to a disturbance source. Recent work has demonstrated that the strain field in the vadose zone shows a response to disturbances in the underlying aquifer. Various instrument designs are available that can record vertical and horizontal strain to very precise resolutions (10-9 ε). Measuring the strain in the vadose zone can present lower costs than measuring pressure in a monitoring well, so using strain could improve the resolution of ZNR tests, as well as make these tests potentially cheaper. A ZNR periodic pumping system was constructed and used to perform tests in the Clemson, South Carolina area. The system was designed to pump and inject with a periodic square-wave, so pumping occurs at a constant rate for half the period and is followed by injection at the same rate for half the period resulting in no net gain or loss of water from the aquifer. The system is designed to generate flow rates from 1 to 3.5 gallons per minute (gpm) with a capacity of 875-gallon per half period. Pressure and strain monitoring points are located in the vicinity of the pumping well. A series of 10 ZNR periodic pumping tests were conducted at the field site, with periods ranging from 2 to 540 minutes. Pressure data were collected and recorded in four monitoring wells in various locations around the pumping well. In addition, vertical and horizontal strain and tilt data were recorded in various locations around the pumping well. Traditional aquifer tests, including slug tests and constant-rate pumping tests were conducted at the field site to provide a baseline of aquifer parameter estimates that are compared the results of ZNR tests. Pressure data measured at monitoring wells during ZNR tests were analyzed utilizing an analytical solution (Streltsova, 1988) that assumes a confined aquifer and estimates hydraulic diffusivity and transmissivity using the time lag and amplitude of the pressure signal. The storativity can be separated out from the hydraulic diffusivity using the transmissivity estimate. The time-lag and amplitude of the pressure signals were estimated using a Fourier transform. The time-lag of the primary period of the pressure increases as a linear function of distance from the pumping well, and it increases as the square root of the period of the ZNR test over a range of periods spanning two orders of magnitudes. The first few harmonics follow similar trends, but the higher frequency harmonics diverge from this trend. Strain and tilt data from various instruments in the vadose and saturated zones are periodic during ZNR tests. The time lags of the strain components increase roughly linearly with distance from the well. Many of the time lags of the strain data from the vadose zone are similar to the time lags observed in the pressure data from a similar distance and pumping period. Theoretical experiments were performed to understand how pressure and strain responds to ZNR periodic pumping. The time-lags of the simulated pressure responses match the field data trends. The results were used to validate the Streltsova (1988) solution for periodic pumping to estimate hydraulic diffusivity. Hydraulic diffusivity was estimated to be 1.8 x 10-2 m2s-1 < Dh < 9.0 x 10-2 m2s-1 using pressure data from all the wells during ZNR tests. Transmissivity was estimated to be 0.8 x 10-4 m2s-1 < T < 4.3 x 10-4 m2s-1. Assuming Storativity = T/Dh gives 1.1 x 10-3 < S < 19 x10-3. Transmissivity is used to estimate the hydraulic conductivity, K, by dividing by the assumed aquifer thickness. Two constant-rate pumping tests were conducted and analyzed using two conceptual models: The Hantush (1961) solution was used to analyze data assuming confined conditions and the Neuman (1974) solution was used to analyze data assuming unconfined conditions. The Hantush (1961) solution estimated transmissivity to be 1.6 x 10-4 m2s-1 < T < 1.8 x 10-4 m2s-1 when using data from all the wells. Storativity was estimated to be 4.4 x 10-3 < S < 8.2 x 10-3. Assuming Dh = T/S gives 0.22 x 10-2 m2s-1 < Dh < 3.6 x 10-2 m2s-1. The Neuman (1974) solution estimated transmissivity to be 0.80 x 10-4 m2s-1 < T < 0.84 x 10-4 m2s-1. Total storativity (S + Sy) was estimated to be 22 x 10-3 < S < 110 x 10-3. Storativity excluding the Sy term was estimated to be 2.2 x 10-3 < S < 2.5 x 10-3. Assuming Dh = T/S gives 0.073 x 10-2 m2s-1 < S < 0.36 x 10-2 m2s-1. Two slug tests were performed on the pumping well, PW-2. The Bouwer and Rice (1976) solution for slug tests was used to analyze the data. This solution works for both confined and unconfined conditions, so both conceptual models were covered through the analysis. Hydraulic conductivity was estimated to be 1.95 x 10-6 ms-1 < K < 2.49 x 10-6 ms-1. Statistical comparison show generally no significant difference in the estimates of Dh, K, or S made using ZNR and conventional pumping tests.. However, the total storage estimated by the Neuman (1974) solution for unconfined settings is larger than that estimate using the Streltsove analysis applied to the ZNR tests. This is likely because the Neuman solution considers delayed yield from storage that the water table, whereas the Streltsova (1988) solution assumes confined conditions. The analysis also indicates that there is a statistically significant difference between parameters estimated with conventional slug tests and those measured with either constant-rate pumping tests or ZNR tests. Hydraulic diffusivity is estimated from the time-lag and the distance from the well. The time-lags of the strain data from the vadose zone are similar to the lags of the pressure data from a similar distance. This indicates that hydraulic diffusivities estimated from the strain data measured in the vadose zone would be similar to the estimates from the pressure data measured in the aquifer. There appear to be sign reversals in some of the strain data that were corrected to estimate the time lag. The dissolved oxygen (CO) concentration was measured during several ZNR tests to evaluate the feasibility of using the procedure to increase contaminant degradation kinetics that are related to DO. The results indicate that DO can be increased during the cycling of ZNR tests in some cases.

Smith-Jones, Austin↗
DOE OSTI2023

Spectral induced polarization monitoring of induced calcite precipitation in subsurface sediments

Co-precipitation of contaminants within the crystalline structure of calcite is a promising natural attenuation or remedial technology being considered at contaminated sites. We explore the sensitivity of the spectral induced polarization (SIP) method to induced calcite precipitation in natural sediments as a path forward to non-invasively monitor these sites. We performed time-lapse column experiments using phased (I-IV) injections over 40 days on natural sediments from the Hanford Site (Washington State, USA). In the phased injections, abiotic calcite precipitation was induced and confirmed to have occurred. Previous work in glass beads and homogeneous sand was limited to high frequency detection of calcite, however in this work we observed the development of two polarization mechanisms, one at high frequency (>100 Hz) and one at low frequency (< 100 Hz). Based on the characteristic frequencies from the SIP high and low frequency regimes, characteristic length scales (L) were computed where the adsorption mode of Na+ versus Ca2+ was compared by using diffusion coefficients corresponding to Na+ versus an arithmetically averaged value for Na+ and Ca2+. Using the diffusion coefficient of Na+, the high frequency L was found to correlate well with the size of the calcite crystals. The low frequency L correlated well with the individual natural sediment grain sizes within the columns. During late experimental times (day 36 and 40), the characteristic low frequency in two of the experimental columns shifted to lower frequencies (<0.001 Hz) which may signify SIP sensitivity of the formed calcite with the sediment grains. In field applications, the development of a low frequency polarization length scale to monitor calcite precipitation is promising for field monitoring applications, however further laboratory work needs to be performed to examine the SIP sensitivity of calcite formation in the presence of to natural sediments

calcium carbonate, calcite, precipitation, spectra↗
DOE OSTI2025

Minimizing particle aggregation in Sm 2 Fe 17 N 3 powders: A CaO-assisted reduction-diffusion approach

We report a novel synthesis of Sm 2 Fe 17 N 3 powders using a CaO-assisted reduction-diffusion (RD) approach. CaO plays a crucial role during mechanochemical processing – acting both as a dispersant and a surface coating agent, which helps to prevent agglomeration due to sintering of Sm 2 Fe 17 particles during the RD step. Together with the added dispersant, the CaO by-product formed during RD can be easily removed in the washing step resulting in fewer aggregated Sm 2 Fe 17 N 3 particles and lower oxygen contamination in the final product. The impact of varying CaO amounts was examined, and synthesis conditions were optimized to achieve phase-pure Sm 2 Fe 17 N 3 powders with less aggregation of magnetic particles. The powders synthesized with addition of 1 wt% CaO as dispersant exhibited the highest hard-magnetic properties: a coercivity (H c ) of 10.7 kOe and a maximum energy product ((BH) max ) of 17.3 MGOe. By densifying the Sm 2 Fe 17 N 3 powders using high-pressure spark plasma sintering, a bulk magnet with a (BH) max of 21.1 MGOe with 88 % of theoretical density was produced. In conclusion, reducing aggregation of the Sm 2 Fe 17 N 3 increases coercivity and remanence of these magnets.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Mitigation of Molten Salt Corrosion

The desire to increase the peak operational temperature of Concentrated Solar Power (CSP) plants to lower the levelized cost of CSP-derived electricity, has led to consideration of an earth-abundant, low-cost, thermally-stable, molten MgCl 2 -KCl-NaCl (40-40-20 mole%) salt for high-temperature thermal energy storage and heat transfer. The reliable, long-term containment of such a molten salt is challenging, however, due to the corrosion of cost-effective, high-temperature structural metal alloys in this salt (particularly if this salt is contaminated with oxygen-bearing species). This project has been focused on: i) development of equipment for evaluating corrosion kinetics in molten chlorides under well-controlled thermal, fluid dynamic, and atmospheric conditions at CSP-relevant temperatures, ii) use of such equipment to evaluate the influence of forced convection on corrosion kinetics (to determine the role of liquid phase diffusion on such corrosion), and iii) identifying a corrosion-mitigation strategy that can be industrially applied to achieve a projected annual corrosion ≤30 μm at CSP-relevant temperatures. The overall goals of this project have been: i) to understand the rate-limiting step(s) and associated kinetic laws governing corrosion in molten chlorides at 750°C, ii) to determine key operational parameters affecting the corrosion kinetics in molten salts, and iii) to demonstrate the tailoring of such parameters to appreciably reduce such corrosion.

14 SOLAR ENERGY↗
DOE OSTI2020

Galaxy Cluster Contribution to the Diffuse Extragalactic Ultraviolet Background

The diffuse ultraviolet background radiation has been mapped over most of the sky with 2′ resolution using data from the Galaxy Evolution Explorer survey. We utilize this map to study the correlation between the UV background and clusters of galaxies discovered via the Sunyaev–Zeldovich effect in the Planck survey. We use only high Galactic latitude (|b|>60{sup ∘}) galaxy clusters to avoid contamination by Galactic foregrounds, and we only analyze clusters with a measured redshift. This leaves us with a sample of 142 clusters over the redshift range of 0.02 ≤ z ≤ 0.72, which we further subdivide into four redshift bins. In analyzing our stacked samples binned by redshift, we find evidence for a central excess of UV background light compared to local backgrounds for clusters with z < 0.3. We then stacked these z < 0.3 clusters to find a statistically significant excess of 12 ± 2.3 photon cm{sup −2} s{sup −1} sr{sup −1} Å{sup −1} over the median of ∼380 photon cm{sup −2} s{sup −1} sr{sup −1} Å{sup −1} measured around random blank fields. We measure the stacked radial profile of these clusters, and find that the excess UV radiation decays to the level of the background at a radius of ∼1 Mpc, roughly consistent with the maximum radial extent of the clusters. Analysis of possible physical processes contributing to the excess UV brightness indicates that non-thermal emission from relativistic electrons in the intracluster medium and faint, unresolved UV emission from cluster member galaxies and intracluster light are likely the dominant contributors.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Choosing the Best Modeling Platform for Radiological Risk Assessment Models - 20468

Radiological risk assessments, in the form of performance or safety assessments, are often required under regulations or guidance for remediation of contaminated land, decommissioning of contaminated buildings or structures, and radioactive waste disposal. These risk assessments are usually supported by fate and transport models that address decay and ingrowth of radionuclides, as well as their movement through engineered systems and the natural environment. These models are often projected thousands, or more, years into the future, largely because the radioactive species change through decay and ingrowth, and hence the magnitude of the radioactive effect changes with time. There are many computer codes that are available to address this type of modeling. They range from addressing specific pathways or processes such as infiltration of water, groundwater, surface water, air, biota, diffusion and advection of water and gases, to those that try to couple all processes together to evaluate the impact of fate and transport through the entire system to places in space and time to which access is assumed. These different types of codes are sometimes separated with the monikers process-level and systems-level codes, although it is often not clear that this separation does justice to the capabilities of the many codes that are available to evaluate fate and transport of radionuclides. The focus of this paper is the latter group of modeling codes. Several systems level modeling codes exist and are used. There are differences between these codes in terms of utility, flexibility, complexity and cost. The purpose of this paper is to compare a few of these codes in the context of work currently being performed by the International Atomic Energy Agency (IAEA) Modeling and Data for Radiological Impact Assessments (MODARIA) II Working Group 1 (WG1). The MODARIA II WG1's main focus is how stakeholder engaged decision analysis can, or should, be applied to radiological contamination problems so that better, longstanding, sustainable, solutions are reached. However, the WG1 also recognizes the potential impact of the modeling tools that are chosen to address radiological risk, which is often a primary objective of decision making for radiological problems. Other objectives might also be important, such as constraining costs, obtaining financing, minimizing impact on ecosystems, saving cultural resources, saving jobs, farmland, environmental justice, etc., in a full decision analysis for a given radiological contamination problem, but none of these other objectives have the same types of complex modeling needs as minimize radiological dose. Consequently, a further focus of WG1 is to evaluate the potential impacts on decision making of the choice of fate and transport, and risk assessment, modeling codes that are used to support decision making. The WG1 will produce a report at the end of 2020 that will focus on an approach to effective decision making and stakeholder engagement. The report will also consider the role that performance assessment modeling should play in the decision-making process, including the impact of the choice of modeling tools or computer codes on risk-informed decision making. Several sites around the World have been made available by Member States for these model comparisons, and several modeling tools have been considered. However, the focus of this paper is on two of the sites, one in Belgium and one in Ukraine, and on three of the tools: NORMALYSA (NORM And Legacy Site Assessment); GoldSim{sup C}, and AMBER{sup C}. The final report from this working group will also cover other modeling tools, including RESRAD, and PC-Cream{sup R}. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗