Search NASA⌕ Search

SEARCH · Search NASA

Results for “matrix 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.

826 records · Page 46

Phase Field Modeling of Chemical Reaction Related Damage Evolution in Environmental Barrier Coatings

The advent of next-generation engines necessitates materials capable of withstanding temperatures beyond the reach of current superalloys. SiC-based ceramic matrix composites, augmented with environmental barrier coatings (EBCs), present a promising materials solution. Given the active search for effective and durable EBCs, there is a pressing need for modeling tools to understand and predict damage evolution in these materials to help accelerate their development. This study introduces a phase-field model (PFM) designed to simulate the thermally grown oxides (TGO) and phase transformation in the degradation and failure of EBCs. The model accounts for the severe volume expansion due to oxidation, alongside phase transformations and porosity evolution during thermal cycling, offering a comprehensive view of the damage processes. Simulation results are validated against experimental findings reported in the literature, establishing the model's potential as a significant tool for understanding and improving the resilience of EBCs in cyclic oxidative environments.

fast-diffusion path↗

Improvement of small to large grain A15 ratio in Nb 3 Sn PIT wires by inverted multistage heat treatments

The next generation of superconducting accelerator magnets for the Large Hadron Collider at CERN will require large amounts of Nb 3 Sn superconducting wires and the Powder-In-Tube (PIT) process, which utilizes a NbSn 2 -rich powder core within tubes of Nb(7.5wt%Ta) contained in a stabilizing Cu matrix, is a potential candidate. But, the critical current density, J c , is limited by the formation of a large grain (LG) A15 layer which does not contribute to transport current, but occupies 25-30% of the total A15 area. Thus it is important to understand how this layer forms, and if it can be minimized in favor of the beneficial small grain (SG) A15 morphology which carries the supercurrent. The ratio of SG/LG A15 is our metric here, where an increase signals improvement in the wires A15 morphology distribution. We have made a critical new observation that the initiation of the LG A15 formation can be controlled at a wide range of temperatures relative to the formation of the small grain (SG) A15. The LG A15 can be uniquely identified as a decomposition product of the Nb6Sn5(Cu x ), surrounded by a layer of rejected Cu, thus the LG A15 is not only of low pin density, but is not continuous grain to grain. We have found that in single stage reactions limited to 630 °C - 690 °C, the maximum SG A15 layer thickness prior to LG A15 formation is very sensitive to temperature, with a maximum around 670 °C. This result led to the design of four novel heat treatments which all included a short, high temperature stage early in the reaction, followed by a slow cooling to a more typical reaction temperature of 630 °C. We also found that this heat treatment (HT) modification increased the SG A15 layer thickness while simultaneously suppressing LG A15 morphology, with no additional consumption of the diffusion barrier. In the best heat treatment the SG/LG A15 ratio improved by 30%. Unfortunately, J c values suffered slightly, however further exploration of this high temperature reaction region is required to understand the limits to A15 formation in Nb3Sn PIT conductors.

36 MATERIALS SCIENCE↗

Experimental studies on the formation of D{sub 2}O and D{sub 2}O{sub 2} by implantation of energetic D{sup +} ions into oxygen ices

The formation of water (H{sub 2}O) in the interstellar medium is intrinsically linked to grain-surface chemistry; thought to involve reactions between atomic (or molecular) hydrogen with atomic oxygen (O), molecular oxygen (O{sub 2}), and ozone (O{sub 3}). Laboratory precedent suggests that H{sub 2}O is produced efficiently when O{sub 2} ices are exposed to H atoms (∼100 K). This leads to the sequential generation of the hydroxyperoxyl radical (HO{sub 2}), then hydrogen peroxide (H{sub 2}O{sub 2}), and finally H{sub 2}O and a hydroxyl radical (OH); despite a barrier of ∼2300 K for the last step. Recent detection of the four involved species toward ρ Oph A supports this general scenario; however, the precise formation mechanism remains undetermined. Here, solid O{sub 2} ice held at 12 K is exposed to a monoenergetic beam of 5 keV D{sup +} ions. Products formed during the irradiation period are monitored through FTIR spectroscopy. O{sub 3} is observed through seven archetypal absorptions. Three additional bands found at 2583, 2707, and 1195 cm {sup –1} correspond to matrix isolated DO{sub 2} (ν{sub 1}) and D{sub 2}O{sub 2} (ν{sub 1}, ν{sub 5}), and D{sub 2}O (ν{sub 2}), respectively. During subsequent warming, the O{sub 2} ice sublimates, revealing a broad band at 2472 cm{sup –1} characteristic of amorphous D{sub 2}O (ν{sub 1}, ν{sub 3}). Sublimating D{sub 2}, D{sub 2}O, D{sub 2}O{sub 2}, and O{sub 3} products were confirmed through their subsequent detection via quadrupole mass spectrometry. Reaction schemes based on both thermally accessible and suprathermally induced chemistries were developed to fit the observed temporal profiles are used to elucidate possible reaction pathways for the formation of D{sub 2}-water. Several alternative schemes to the hydrogenation pathway (O{sub 2}→HO{sub 2}→H{sub 2}O{sub 2}→H{sub 2}O) were identified; their astrophysical implications are briefly discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Characterization and Modeling of Grain Boundary Chemistry Evolution in Ferritic Steels under Irradiation

Ferritic/martensitic (FM) steels such as HT-9, T-91 and NF12 with chromium concentrations in the range of 9-12 at.% Cr and high Cr ferritic steels (oxide dispersion strengthened steels with 12-18% Cr) are receiving increasing attention for advanced nuclear applications, e.g. cladding and duct materials for sodium fast reactors, pressure vessels in Generation IV reactors and first wall structures in fusion reactors, thanks to their advantages over austenitic alloys. Predicting the behavior of these alloys under radiation is an essential step towards the use of these alloys. Several radiation-induced phenomena need to be taken into account, including phase separation, solute clustering, and radiation-induced segregation or depletion (RIS) to point defect sinks. RIS at grain boundaries has raised significant interest because of its role in irradiation assisted stress corrosion cracking (IASCC) and corrosion of structural materials. Numerous observations of RIS have been reported on austenitic stainless steels where it is generally found that Cr depletes at grain boundaries, consistently with Cr atoms being oversized in the fcc Fe matrix. While FM and ferritic steels are also subject to RIS at grain boundaries, unlike austenitic steels, the behavior of Cr is less clear with significant scatter and no clear dependency on irradiation condition or alloy type. In addition to the lack of conclusive experimental evidence regarding RIS in F-M alloys, there have been relatively few efforts at modeling RIS behavior in these alloys. The need for predictability of materials behavior and mitigation routes for IASCC requires elucidating the origin of the variable Cr behavior. A systematic detailed high-resolution structural and chemical characterization approach was applied to ion-implanted and neutron-irradiated model Fe-Cr alloys containing from 3 to 18 at.% Cr. Atom probe tomography analyses of the microstructures revealed slight Cr clustering and segregation to dislocations and grain boundaries in the ion-irradiated alloys. More significant segregation was observed in the neutron irradiated alloys. For the more concentrated alloys, irradiation did not affect existing Cr segregation to grain boundaries and segregation to dislocation loops was not observed perhaps due to a change in the dislocation loop structure with increasing Cr concentration. Precipitation of α’ was observed in the neutron irradiated alloys containing over 9 at.% Cr. However ion irradiation appears to suppress the precipitation process. Initial low dose ion irradiation experiments strongly suggest a cascade recoil effect. The systematic analysis of grain boundary orientation on Cr segregation was significantly challenged by carbon contamination during ion irradiation or by existing carbon and therefore carbide formation at grain boundaries (sensitization). The combination of the proposed systematic experimental approach with atomistic modeling of diffusion processes directly addresses the programmatic need for developing and benchmarking predictive models for material degradation taking into account atomistic kinetics parameters

36 MATERIALS SCIENCE↗

Effectiveness of an in-situ permeable reactive barrier in removing strontium-90 from groundwater - 15418

An in-situ permeable reactive barrier (PRB) is being used to remove strontium-90 from groundwater along the 100-N Area shoreline at the U.S. Department of Energy Hanford Site in southeastern Washington. The PRB was formed by emplacing apatite, a calcium phosphate mineral, in the subsurface sediments via injection of calcium-citrate and sodium-phosphate solution into the aquifer through a network of vertical wells. The relatively slow biodegradation of citrate allowed sufficient time for the injected solution to disperse in the PRB treatment zone as the free calcium and phosphate combined to form apatite. Substitution of strontium for calcium in the apatite crystal structure is thermodynamically favorable. As groundwater flows through the PRB toward the Columbia River, strontium-90 is sequestered into the apatite crystal matrix. The sequestered strontium-90 undergoes radioactive decay in place within the PRB. Groundwater monitoring wells installed between the PRB and the Columbia River are used to monitor the PRB effectiveness. The concentration of strontium-90 in the monitoring wells following apatite emplacement are compared to the pre-injection concentrations. The performance goal of the PRB is to reduce the concentration of strontium-90 in groundwater down gradient of the PRB by 90 percent, compared to the pre-injection concentrations. Groundwater samples in Fall 2013 show the percent reduction in strontium-90 concentrations ranged from 71 to 98 percent in eleven of the twelve monitoring wells. The data indicate that the PRB is effective in reducing the flux of strontium-90 toward the Columbia River and is thereby instrumental in protection of the Columbia River. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Effects of dose rate on the microstructure and deuterium retention in γ-LiAlO 2 pellets

This report presents experimental findings obtained from November 2024 to September 2025. Defect accumulation and microstructural evolution during ion irradiation at elevated temperatures are governed by two competing processes: defect production, driven by dose rate, and defect recovery, controlled by defect diffusion, interaction, and annihilation. At a given dose, the resulting microstructural evolution depends on both the dose rate and irradiation temperature. As a continuation of our FY24 tritium science project, which investigated temperature effects at a fixed dose rate, this study focuses on dose-rate effects at a fixed temperature to provide deeper insights into the irradiated microstructure and compositional changes in γ-LiAlO 2 pellets. The study aims to elucidate the impact of dose rate on microstructure, precipitate morphology, deuterium retention, and lithium volatility in γ-LiAlO 2 pellets under sequential 120 keV He + and 80 keV D 2 + ion irradiation. Three dose rates of 7.3×10 4 , 2.9×10 4 and 6.8×10 5 dpa/s were applied to achieve the same total ion fluence of 2×10 17 (He + +D + )/cm 2 at 500 °C, corresponding to a maximum combined dose of 7.55 dpa at ~255 nm. The irradiated pellets were subsequently characterized using scanning transmission electron microscopy (STEM) and time-of-flight secondary ion mass spectrometry (ToF-SIMS). The microstructural response of γ-LiAlO 2 to ion irradiation was found to be strongly dose-rate dependent. At medium and high dose rates, irradiation produced a surface amorphized layer and an underlying crystalline layer containing LiAl 5 O 8 precipitates, with implanted gases accumulating and forming blisters at the crystalline-amorphous interface. It remains to be investigated whether the amorphized layer was produced by He + ion irradiation prior to D 2 + ion irradiation. In contrast, at low dose rates, the material remained crystalline, with cavities, likely gas-filled, distributed around precipitates, within the γ-LiAlO 2 matrix, and along grain boundaries. While precipitate morphology exhibited anisotropy, their size showed little sensitivity to dose rate in the applied range of this study. This result, however, does not rule out the possibility that further reductions in dose rate could influence precipitate size. High dose-rate irradiation enhanced protonium–deuterium isotopic exchange and lithium depletion in the amorphized region. Collectively, the results show that dose rate governs amorphization, gas redistribution, isotopic exchange, and lithium depletion, providing important insights into the mechanisms underlying structural evolution in γ-LiAlO 2 pellets under reactor-relevant irradiation conditions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DECOVALEX-2023: Task F1 Final Report

DECOVALEX-2023 Task F is a comparison of models and methods for post-closure performance assessment (PA) of a deep geologic repository for radioactive waste. The general aims of Task F are to build confidence in the models, methods, and software used for PA and to stimulate additional research and development in PA methodologies. The task objectives are to motivate development of PA modelling skills and capabilities, to examine the influence of model choices on calculated repository performance, and to compare the uncertainties introduced by model choices to other sources of uncertainty. Task F involves no actual experiment or site. It is a PA modelling exercise that requires the conceptual development of hypothetical repository designs and geologic settings. Because three of the teams were interested in salt and the rest of the teams were interested in crystalline rock, Task F was split into two branches: Task F1 for crystalline rock and Task F2 for salt. This report is for Task F1, crystalline rock. Teams from seven countries (Canada, Czech Republic, Germany, Korea, Sweden, Taiwan, and United States) participated in Task F1. The teams worked together to define the features, events, and processes of the reference case repository and established a set of performance measures. In addition, they defined a set of benchmark problems designed to test and compare modelling capabilities for fracture flow and transport at different scales. The repository design and benchmark problems are documented in a Task Specification that evolved over time as the group honed the specifications. The benchmark problems verified that each team can aptly model flow and transport in fractured media in 1-, 2-, and 3-dimensions. Two general approaches were used for the 3-dimensional benchmarks: discrete fracture network (DFN) and equivalent continuous porous medium (ECPM). DFN modelling involves explicit meshing of each fracture while ECPM modelling aims to capture the effective porosity and directional permeability of each cell in a space-filling mesh as affected by intersecting fractures. In some models, a combination of the two is used, i.e., DFN for large known fractures and ECPM for the rest of the domain. Transport is solved by using either the advection-dispersion equation or particle tracking. Although some variation is observed among model breakthrough curves in the benchmark problems, there is strong agreement in breakthrough behaviour up to at least the 75 th percentile for all benchmarks. At the 90 th percentile, breakthrough results show larger differences, suggesting several models retain substantially higher fractions of tracer in regions of slower moving water. In addition to the flow and transport benchmarks, several teams completed the source term benchmark, verifying capabilities for modelling radionuclide decay and ingrowth, waste package breach, instant release fractions, fuel matrix degradation rates, and radionuclide solubility limitations. The reference case is conceptualized as a generic spent fuel repository at a depth of 450 m in fractured crystalline rock. The repository has 50 parallel backfilled drifts, each with 50 deposition holes 6 m apart. Each deposition hole contains a 4-PWR waste package and bentonite buffer. The rock domain is 5 km in length, 2 km in width, and 1 km in depth. It has 6 deterministic fractured deformation zones and a multitude of stochastic fractures. Teams generally used the ECPM approach for the entire rock or a hybrid approach in which the deterministic fracture zones are modelled with a DFN and the rest of the rock is modelled by ECPM. Of the reference case problems specified, only the results of the initial reference case problem are compared in this report. The initial problem focuses on transport from the deposition holes to the surface, i.e., it neglects waste package performance. Tracers are released at all waste package locations at time zero and tracked for their releases to the near field and ground surface. The water fluxes calculated at the ground surface entry and exit regions of the domain are similar for all models except for two that have considerably lower fluxes. For tracer transport, large differences are observed among models in the magnitude of tracer transported. Much of the difference appears to be due to how the repository is implemented and hence the different degrees of repository simplification. Models that exclude the drifts, buffer, and backfill from the domain tend to show greater release of tracers and radionuclides from the repository. The initial study presented here indicates that major differences in modelling important processes within the repository (e.g., diffusion through buffer and backfill) can produce broadly different release and transport results, especially when those processes are excluded. Even for the models that included all specified features, events, and processes, the results show significant differences and demonstrate the importance of examining multiple modelling approaches in performance assessment. The differences in results observed in this study are expected to motivate teams to either increase complexity in future versions of the reference case models or to improve methods to account for the effects of simplified features and processes. Either way, future improvements in these models are expected to produce results that more closely agree.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Anisotropy in thermal conductivity of graphite flakes–SiC{sub p}/matrix composites: Implications in heat sinking design for thermal management applications

Within the frame of heat dissipation for electronics, a very interesting family of anisotropic composite materials, fabricated by liquid infiltration of a matrix into preforms of oriented graphite flakes and SiC particles, has been recently proposed. Aiming to investigate the implications of the inherent anisotropy of these composites on their thermal conductivity, and hence on their potential applications, materials with matrices of Al–12 wt.% Si alloy and epoxy polymer have been fabricated. Samples have been cut at a variable angle with respect to the flakes plane and thermal conductivity has been measured by means of two standard techniques, namely, steady state technique and laser flash method. Experimental results are presented and discussed in terms of current models, from which important technological implications for heat sinking design can be derived. - Highlights: • Anisotropy in thermal conductivity of graphite flakes-based composites is evaluated. • Samples are cut in a direction forming a variable angle with the oriented flakes. • For angles 0° and 90°, thermal conductivity does not depend on sample geometry. • For intermediate angles, thermal conductivity strongly depends on sample geometry. • “Thin” samples must be thicker than 600 μm, “thick” samples must be encapsulated.

36 MATERIALS SCIENCE↗

A Versatile Integrated Ambient Ionization Source Platform

The pursuit of high-throughput sample analysis from complex matrix demands development of multiple ionization techniques with complementary specialties. A versatile integrated ambient ionization source (iAmIS) platform is proposed in this work, based on the idea of integrating multiple functions, enhancing the efficiency of current ionization techniques, extending the applications, and decreasing the cost of the instrument. The design of the iAmIS platform combines flowing atmospheric pressure afterglow (FAPA) source/direct analysis in real time (DART), dielectric barrier discharge ionization (DBDI)/low-temperature plasma (LTP), desorption electrospray ionization (DESI), and laser desorption (LD) technique. All individual and combined ionization modes can be easily attained by modulating parameters. In particular, the FAPA/DART&DESI mode can realize the detection of polar and nonpolar compounds at the same time with two different ionization mechanisms: proton transfer and charge transfer. The introduction of LD contributes to the mass spectrometry imaging and the surface-assisted laser desorption (SALDI) under ambient condition. Compared with other individual or multi-mode ion source, the iAmIS platform provides the flexibility of choosing different ionization modes, broadens the scope of the analyte detection, and facilitates the analysis of complex samples. .

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

QMCPACK: Advances in the development, efficiency, and application of auxiliary field and real-space variational and diffusion quantum Monte Carlo

We review recent advances in the capabilities of the open source ab initio Quantum Monte Carlo (QMC) package QMCPACK and the workflow tool Nexus used for greater efficiency and reproducibility. The auxiliary field QMC (AFQMC) implementation has been greatly expanded to include k-point symmetries, tensor-hypercontraction, and accelerated graphical processing unit (GPU) support. These scaling and memory reductions greatly increase the number of orbitals that can practically be included in AFQMC calculations, increasing the accuracy. Advances in real space methods include techniques for accurate computation of bandgaps and for systematically improving the nodal surface of ground state wavefunctions. Results of these calculations can be used to validate application of more approximate electronic structure methods, including GW and density functional based techniques. To provide an improved foundation for these calculations, we utilize a new set of correlation-consistent effective core potentials (pseudopotentials) that are more accurate than previous sets; these can also be applied in quantum-chemical and other many-body applications, not only QMC. These advances increase the efficiency, accuracy, and range of properties that can be studied in both molecules and materials with QMC and QMCPACK.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Bimolecular reaction of CH{sub 3} + CO in solid p-H{sub 2}: Infrared absorption of acetyl radical (CH{sub 3}CO) and CH{sub 3}-CO complex

We have recorded infrared spectra of acetyl radical (CH{sub 3}CO) and CH{sub 3}-CO complex in solid para-hydrogen (p-H{sub 2}). Upon irradiation at 248 nm of CH{sub 3}C(O)Cl/p-H{sub 2} matrices, CH{sub 3}CO was identified as the major product; characteristic intense IR absorption features at 2990.3 (ν{sub 9}), 2989.1 (ν{sub 1}), 2915.6 (ν{sub 2}), 1880.5 (ν{sub 3}), 1419.9 (ν{sub 10}), 1323.2 (ν{sub 5}), 836.6 (ν{sub 7}), and 468.1 (ν{sub 8}) cm{sup −1} were observed. When CD{sub 3}C(O)Cl was used, lines of CD{sub 3}CO at 2246.2 (ν{sub 9}), 2244.0 (ν{sub 1}), 1866.1 (ν{sub 3}), 1046.7 (ν{sub 5}), 1029.7 (ν{sub 4}), 1027.5 (ν{sub 10}), 889.1 (ν{sub 6}), and 723.8 (ν{sub 7}) cm{sup −1} appeared. Previous studies characterized only three vibrational modes of CH{sub 3}CO and one mode of CD{sub 3}CO in solid Ar. In contrast, upon photolysis of a CH{sub 3}I/CO/p-H{sub 2} matrix with light at 248 nm and subsequent annealing at 5.1 K before re-cooling to 3.2 K, the CH{sub 3}-CO complex was observed with characteristic IR features at 3165.7, 3164.5, 2150.1, 1397.6, 1396.4, and 613.0 cm{sup −1}. The assignments are based on photolytic behavior, observed deuterium isotopic shifts, and a comparison of observed vibrational wavenumbers and relative IR intensities with those predicted with quantum-chemical calculations. This work clearly indicates that CH{sub 3}CO can be readily produced from photolysis of CH{sub 3}C(O)Cl because of the diminished cage effect in solid p-H{sub 2} but not from the reaction of CH{sub 3} + CO because of the reaction barrier. Even though CH{sub 3} has nascent kinetic energy greater than 87 kJ mol{sup −1} and internal energy ∼42 kJ mol{sup −1} upon photodissociation of CH{sub 3}I at 248 nm, its energy was rapidly quenched so that it was unable to overcome the barrier height of ∼27 kJ mol{sup −1} for the formation of CH{sub 3}CO from the CH{sub 3} + CO reaction; a barrierless channel for formation of a CH{sub 3}-CO complex was observed instead. This rapid quenching poses a limitation in production of free radicals via bimolecular reactions in p-H{sub 2}.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Noise from undetected sources in Dark Energy Survey images

ABSTRACT For ground-based optical imaging with current CCD technology, the Poisson fluctuations in source and sky background photon arrivals dominate the noise budget and are readily estimated. Another component of noise, however, is the signal from the undetected population of stars and galaxies. Using injection of artifical galaxies into images, we demonstrate that the measured variance of galaxy moments (used for weak gravitational lensing measurements) in Dark Energy Survey (DES) images is significantly in excess of the Poisson predictions, by up to 30 per cent, and that the background sky levels are overestimated by current software. By cross-correlating distinct images of ‘empty’ sky regions, we establish that there is a significant image noise contribution from undetected static sources (US), which, on average, are mildly resolved at DES resolution. Treating these US as a stationary noise source, we compute a correction to the moment covariance matrix expected from Poisson noise. The corrected covariance matrix matches the moment variances measured on the injected DES images to within 5 per cent. Thus, we have an empirical method to statistically account for US in weak lensing measurements, rather than requiring extremely deep sky simulations. We also find that local sky determinations can remove most of the bias in flux measurements, at a small penalty in additional, but quantifiable, noise.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Quantitative estimation of localization errors of 3d transition metal pseudopotentials in diffusion Monte Carlo

The necessarily approximate evaluation of non-local pseudopotentials in diffusion Monte Carlo (DMC) introduces localization errors. In this paper, we estimate these errors for two families of non-local pseudopotentials for the first-row transition metal atoms Sc–Zn using an extrapolation scheme and multideterminant wavefunctions. Sensitivities of the error in the DMC energies to the Jastrow factor are used to estimate the quality of two sets of pseudopotentials with respect to locality error reduction. The locality approximation and T-moves scheme are also compared for accuracy of total energies. After estimating the removal of the locality and T-moves errors, we present the range of fixed-node energies between a single determinant description and a full valence multideterminant complete active space expansion. The results for these pseudopotentials agree with previous findings that the locality approximation is less sensitive to changes in the Jastrow than T-moves yielding more accurate total energies, however not necessarily more accurate energy differences. For both the locality approximation and T-moves, we find decreasing Jastrow sensitivity moving left to right across the series Sc–Zn. The recently generated pseudopotentials of Krogel et al. reduce the magnitude of the locality error compared with the pseudopotentials of Burkatzki et al. by an average estimated 40% using the locality approximation. The estimated locality error is equivalent for both sets of pseudopotentials when T-moves is used. Finally, for the Sc–Zn atomic series with these pseudopotentials, and using up to three-body Jastrow factors, our results suggest that the fixed-node error is dominant over the locality error when a single determinant is used.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Transport properties of damaged materials. Cementitious barriers partnership

The objective of the Cementitious Barriers Partnership (CBP) project is to develop tools to improve understanding and prediction of the long-term structural, hydraulic, and chemical performance of cementitious barriers used in low-level waste storage applications. One key concern for the long-term durability of concrete is the degradation of the cementitious matrix, which occurs as a result of aggressive chemical species entering the material or leaching out in the environment, depending on the exposure conditions. The objective of the experimental study described in this report is to provide experimental data relating damage in cementitious materials to changes in transport properties, which can eventually be used to support predictive model development. In order to get results within a reasonable timeframe and to induce as much as possible uniform damage level in materials, concrete samples were exposed to freezing and thawing (F/T) cycles. The methodology consisted in exposing samples to F/T cycles and monitoring damage level with ultrasonic pulse velocity measurements. Upon reaching pre-selected damage levels, samples were tested to evaluate changes in transport properties. Material selection for the study was motivated by the need to get results rapidly, in order to assess the relevance of the methodology. Consequently, samples already available at SIMCO from past studies were used. They consisted in three different concrete mixtures cured for five years in wet conditions. The mixtures had water-to-cement ratios of 0.5, 0.65 and 0.75 and were prepared with ASTM Type I cement only. The results showed that porosity is not a good indicator for damage caused by the formation of microcracks. Some materials exhibited little variations in porosity even for high damage levels. On the other hand, significant variations in tortuosity were measured in all materials. This implies that damage caused by internal pressure does not necessarily create additional pore space in the microstructure, but likely creates new thin pathways between existing pore space for species to travel. These results have a significant impact on modeling efforts. Models relating porosity to tortuosity and permeability are unlikely to provide the correct basis for predicting long-term durability of concrete sustaining internal pressures and microcrack formation. Other avenues like the modeling of internal crystallization pressure need to be explored.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Fuel Matrix Degradation Model Development Update: Alloy Corrosion Rates and Hydrogen Generation

This work is being performed as part of the DOE NE Spent Fuel and Waste Science and Technology Campaign, Argillite and Crystalline Rock R&D work packages: SF-19AN01030101 and SF-19AN01030201. This document meets the milestone M4SF-19AN010301013 for Argillite R&D and the milestone M4SF-19AN010302013 for Crystalline R&D. The fuel matrix degradation (FMD) model calculates the degradation rate of spent UO2 fuel based on fundamental electrochemical theory and thermodynamics. The model provides the quantitative basis for radionuclide source term predictions in the generic disposal system analysis (GDSA) performance assessment (PA) model. The FMD model has been implemented in a manner that facilitates its integration with the GDSA-PA code and a preliminary, integrated FMD-GDSA model has been successfully tested (Jerden et al., 2017). The specific focus of on-going work is to accurately and quantitatively represent the generation of H2 in a breached waste package and its effect on the spent fuel degradation rate so that this key process can be accurately represented in PA models. It has been shown experimentally that millimolar concentrations of dissolved H2 in contact with spent fuel will inhibit oxidative dissolution and decrease the fuel degradation rates by as much as three orders of magnitude (e.g., Röllin et al., 2001, Ollila, 2008). This will lead to a significant decrease in the radionuclide source term values used in repository PA models. Sensitivity studies using the FMD model have shown that the dissolved H2 concentration is the dominant environmental variable affecting the UO2 spent fuel dissolution rate as long as a source of H2 persists in the breached waste package (Jerden et al., 2015). Therefore, on-going experimental and modeling efforts focus on H2 production, accumulation and interfacial reactions within a breached waste package. The anoxic corrosion of metallic engineering materials within a breached waste package will be the main source of H2 in crystalline and argillite rock repository systems. The most important H2-producing alloys present within a waste package include stainless steel, carbon steel, and aluminum alloys and the Zircaloy fuel cladding. A steel corrosion module was added to the FMD model to account for these H2 sources and couple the H2 generation rate with fuel degradation processes (Jerden et al., 2017). As discussed below, there is a need for experimental data from electrochemical corrosion tests to quantify the H2 generation rates due to the corrosion of these waste package alloys as a function of key environmental variables (Eh, pH, groundwater chemistry). This report presents results from series of electrochemical corrosion tests with materials present in waste packages that provide information needed to parameterize and validate the FMD model. The alloys tested include 316L stainless steel (316SS), AISI 4320 carbon steel (C-steel) and Zircaloy-4 fuel cladding (unirradiated). Potentiodynamic polarization and potentiostatic polarization tests were performed in buffer solutions at pH 4, 7 and 10 with NaCl concentrations of 0, 4 and 100 millimolar (mM). The potentiostatic tests were performed at voltages varying from 0.5 volts down to -0.66 volts, vs the standard hydrogen electrode (SHE), and all tests were performed at laboratory ambient temperatures (~22 °C). The corrosion rates measured during the electrochemical tests were used as direct input to a new FMD/in-package chemistry model that is being developed to provide accurate long-term spent fuel degradation rates as chemical conditions (Eh, pH, speciation) evolve within a breached waste package. This new model was developed by coupling the FMD model with the reactive transport code X1t, which is a module within the Geochemist’s Workbench (GWB) reaction path modeling code (Jerden et al., 2019). The reactive transport model consists of a 1D domain discretized with 21 reaction diffusion cells and includes a single waste packag

Jerden, James↗

Validation Data Acquisition in HTTF during PCC Events

A validation experimental campaign was conducted in an Integral Effect Test (IET) facility of High Temperature Gas Reactors (HTGR), the High-Temperature Test Facility (HTTF) at Oregon State University (OSU). The HTTF simulates Depressurized and Pressurized Conduction Cooldowns (DCC and PCC). This campaign required the development of a new laser spectroscopic diagnostic to measure velocity in the challenging conditions of the HTTF: low speed (~1 m/s) gas flows at HTGR prototypical temperature and 1/10th pressure. This was a collaborative effort between co-PIs at The George Washington University (GW), Oregon State University (OSU), and NASA Langley Research Center. The main accomplishments of this project include the record for dynamic range for velocimetry, highest accuracy obtained with this technique, successful deployment in an IET leading to new validation matrix for CFD. These are detailed below and in manuscript appended to this executive summary. For this project, we introduced a new class of laser spectroscopic diagnostics to Thermal-Hydraulics to measure velocity of gases; furthermore, the diagnostic was demonstrated in-situ in an IET during DCC events. In such scenarios, particles used in mainstream techniques, like Particle Image Velocimetry (PIV) are not appropriate as they settle down too rapidly and also contaminate the experimental facility. Molecular tracers stay mixed within the working gas and can seed the flow in a technique called Molecular Tagging Velocimetry (MTV). In MTV a molecular tracer is photo-dissociated by a first (write) laser into radicals or molecules. The pattern created by the write laser is then interrogated with planar laser-induced fluorescence (PLIF), the read pulse(s), which are recorded with a camera. The pattern is probed and matched at two times (interval or probe time, dt), resulting in a time-of-flight velocimetry technique. This project demonstrated a new application range for MTV in gases. MTV has been extensively used for high-speed gas dynamics: molecules do not suffer from inertia associated with particles in regions of high accelerations (such as shock waves). For high-speed flows, the probe time is on the order of 0.1-10 µs. In time-of-flight techniques, the velocity is reconstructed from the recorded displacement, dx, between time interval dt: U = dx/dt. To first order, the resolution of the pattern matching is typically a (fixed) fraction of a pixel; therefore the resolution of the velocity measurement increases with dx and, for a given flow velocity, dt. For the low-speed flows of interest here, this requires long probe time, and therefore finding tracers that are stable for a long time. After an extensive literature survey we identified 6 tracers that would be suitable. We also worked with Spectra Physics to design a custom laser system that is flexible enough to probe these various tracers. The final instrument is made of 4 lasers that are integrated into a compact and movable cart to facilitate the deployment of the instrument off-site. One unique aspect of the laser system is the tunable dye laser (needed to do PLIF of the radicals). By being pumped by two Nd:YAG lasers, this laser is able to provide two pulses with any probe time (dt) for wavelengths between 200 and 900 nm with a linewidth of 3 pm. Accommodating the dual-pulse capability required customizing the pump lasers and the dye laser, making this system unique. Additionally, a custom UV-intensified CCD camera capable of frame-straddling and exposure time as short as 10 ns was acquired. More than two-third of the instrument cost (~$200k) was covered by a cost sharing from the PI’s startup funds. Of the 6 potential tracers identified, two were tested extensively (H2O and N2O) at high-temperature and elevated pressure at GW prior to conducting the tests at OSU. For both tracers, an Ar-F laser (193 nm, 10 ns, 10 mJ/pulse) photo-dissociates trace amounts (0.1 – 5% molar fraction) of the gases into OH and NO. OH is probed at 281.905 nm and NO at 226.186 nm (4-30mJ, 10 ns pulse). NO is very stable in inert environment and we could reconstruct velocities precisely with probe delay times as long at 40 ms, compared to 2 ms for OH, resulting in increased precision. Therefore, for the tests in the HTTF, N2O was selected as the final tracer and its optimal concentration was determined to be 0.5% molar. In the GW laboratory, we obtained a velocity measurement precision of 4 mm/s at conditions mimicking the HTTF tests. We also demonstrated the benefit of molecular tracers over particles and the flexibility of our instrument by continuously measuring velocity from 1000 to 0.1 m/s in a blow down test. This represents a dynamic range on velocity of 80 dB, which is unprecedented for velocimetry! Once the diagnostic was optimized, two two-week campaigns were conducted in the HTTF. This experimental campaign required numerous practical solutions to successfully and rapidly deploy this very sensitive table-top diagnostic in-situ. This was accomplished thanks to very extensive planning from all parties involved in the project. Due to constraints on availability of the HTTF, only adiabatic DCC tests were conducted. Specifically, the velocity was measured at the exit of the hot leg into the Reactor Cavity Simulation Tank (RCST). We were able to record the gas velocity associated with the air ingress for long times (up to one hour). Depending on the initial conditions (density of the gases in the Reactor Pressure Vessel and the RCST) the lock exchange lasted 30-60 s past the valve openings. Large vortices were identified within the shear layer at the interface of the cross flow. These vortices are an indication of shear-driven instability between the counter streams of different density. After the lock-exchange, a small residual flow was still present and its strength decayed exponentially. It finds its source in a buoyancy difference between the gases in both vessels due to the mixing and diffusion of the gases in the cross-over duct and the RPV lower plenum. We did not expect to capture this part of the transient, which would have been less significant (if not negligible) in the presence of a heated core; the precision of our measurement in the HTTF was 6 mm/s, which exceeded our expectations for such a complex test! Using commercial software (Fluent) we simulated the experiment with Computational Fluid Dynamic (CFD) using Reynolds Averaged Navier-Stokes (RANS) turbulence models. We were not able to accurately reproduce the shear layer instability, mixing, and diffusion captured in our experimental data, even with large 3D meshes (>10 M points). We attributed this to the inadequacy of the RANS models for the flow of interest: our experimental Reynolds number was too low and was outside the applicability range of the RANS models. Instead, we should have employed CFD with less modeling, such as Large Eddy (LES) or even Direct Numerical (DNS) simulations. Unfortunately, we did not have the computational resources to do so. Finally, it should be noted that this experimental technique holds great promises for measuring velocity in thermalhydraulics. MTV has an application domain (highlighted in the narrative) that is complimentary to more established (and simpler) techniques, such as PIV. It has attracted the attention of several thermalhydraulics researchers and sponsors to measure velocity in harsh environments. The spectroscopic nature of this technique also enables to expand it to measure temperature and concentration simultaneously to velocity. In lieu of a narrative, 4 detailed journal publications are appended to this executive summary to describe the various development and accomplishments during this project. They include 2 published manuscripts (Experiments in Fluids and Measurement Science and Technology), and two invited contributions following the NURETH-17 conference (Nuclear Science and Technology). Our current Nuclear Engineering and Design draft is not included at this time.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗