Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Conditional diffusion models”

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 793 records · Page 44

Ion and Water Dynamics in the Transition from Dry to Wet Conditions in Salt-Doped PEG

The influence of the water content on ion and water transport mechanisms in polymer membranes under low to moderate hydration conditions remains poorly understood. In this study, we combine ion and water diffusivity (PFG-NMR) measurements with atomistic molecular dynamics simulations to better understand transport processes in hydrated salt-doped poly(ethylene glycol). Above the water percolation threshold, the experimental and simulated diffusivities are in good agreement with the free volume transport models. At low hydration levels, unlike dry systems, ion diffusion cannot be described by polymer segmental dynamics alone. We rationalize such observations by the interplay between ion–water and ion–polymer solvation of cations and between ion–water and cation–anion interactions for anions. Further, we demonstrate that a two-state model combining ion–water solvation and free volume transport can describe water dynamics across the entire hydration range of interest. Here, our findings provide a more encompassing analysis of ion and water transport in hydrated polyelectrolytes, specifically in the low hydration regime.

36 MATERIALS SCIENCE↗

Kinetic Model of Gas Transport in Carbon Nanotube Channels

Carbon nanotube channels represent a rare experimental realization of a nanofluidic channel with molecularly smooth walls and nanometer scale inner diameters. This unique combination of properties gives this channel an ability to support enhanced transport of water and gases with fluxes often exceeding those of conventional channels by several orders of magnitude. Surprisingly, most of these transport enhancement phenomena can be explained using very simple mechanisms that hardly go beyond classical physics concepts. In this paper, we present a simplified analytical model that uses classic kinetic theory formalism to describe gas transport in carbon nanotube channels and highlight the role of surface defects and adsorbates in determining transport efficiency. We also extend this description to include the possibility of gas molecule diffusion along the nanotube walls. Our results show that in all cases the conditions at the nanotube channel walls play a critical role in determining the transport efficiency, and that in some cases obtaining efficient transport has to involve optimization of fluxes from diffusion through the gas phase and along the nanotube surface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dynamic Triple-Mode Sorption and Outgassing in Materials

Moisture uptake and outgassing can be detrimental to a system by altering the chemical and mechanical properties of materials within the system over time. In this work, we conducted isotherm experiments to investigate dynamic moisture sorption and desorption in markedly different materials, i.e., a polymeric material, Sylgard-184 and a ceramic aluminosilicate material, Zircar RS-1200, at different temperatures (30 °C–70 °C) by varying the water activity (0.0–0.90). Sylgard-184 showed a linear sorption and outgassing behavior with no-hysteresis over the entire temperature and water activity range considered here. Whereas, the sorption and outgassing of Zircar RS-1200 was highly non-linear with significant hysteresis, especially at higher water activities, at all temperatures considered here. The type of hysteresis suggested the presence of mesopores in Zircar RS-1200, whereas the lack of hysteresis in Sylgard-184 indicates that it has a nonporous structure. A diffusion model coupled with a dynamic, triple-mode sorption (Langmuir, Henry, and pooling modes) model employed in this study matched our experimental data very well and provides mechanistic insight into the processes. Our triple-mode sorption model was adaptive enough to (1) model these distinctly different materials and (2) predict sorption and outgassing under conditions that are distinctly different from the parameterization experiments.

36 MATERIALS SCIENCE↗

Defects in Na-, K-, and Cd-Doped CuInSe$_2$ : Canonical Thermodynamics Based on Ab Initio Calculations

In this work. we extend a canonical-thermodynamic method for computing intrinsic point defect concentrations as a function of chemical stoichiometry and temperature to include extrinsic defect contributions, while applying this method to study defects in CuInSe 2 . This technique relies on a large set of defect formation energies calculated from first principles, which require corrections for known errors arising from spurious interactions. Guided by recent experimental work exhibiting the complex interplay between Na, K, and Cd incorporation, we examine the behavior of dominant defects as the material composition varies under experimentally relevant conditions. In addition to identifying the regions of composition-parameter space relevant to the incorporation of impurities Na, K, and Cd, and stable against the formation of secondary compounds, we further study defect kinetics susceptible to the presence of these impurities. From this analysis, we propose a simple model for the enhanced diffusion of Cu-vacancies mediated by K, which could lead to the enhanced incorporation of Cd observed experimentally.

42 ENGINEERING↗

Conditional Latent Diffusion for High-Resolution Prediction of Electrochemical Surface Morphology

A conditionally guided generative latent diffusion process that is trained on a set of experimental processing parameters and their associated resulting electron microscope images of the electrodeposition process is able to interpolate between processing parameters in a physically consistent way. Electrodeposition of rhenium with pulse and pulse-reverse waveforms is used as a model system, and the process is adaptable to other electrodeposition, electropolishing, or corrosion processes. The method is able to extrapolate, predicting estimates of material morphologies for experimental setups unseen in the training data. The results are demonstrated with experimental data.

36 MATERIALS SCIENCE↗

M4FT-16LL080303052-State of Knowledge for Colloid Facilitated Radionuclide Transport and Update on Actinide Diffusion in Bentonite Backfill

This progress report (Level 4 Milestone Number M4FT-16LL080303052) summarizes research conducted at Lawrence Livermore National Laboratory (LLNL) within the Crystalline Disposal R&D Activity Number FT-16LL080303051 and Crystalline International Collaborations Activity Number FT-16LL080303061. The focus of this research is the interaction of radionuclides with Engineered Barrier System (EBS) and host rock materials at various physico-chemical conditions relevant to subsurface repository environments. They include both chemical and physical processes such as solubility, sorption, and diffusion. The colloid-facilitated transport effort focused on preparation of a draft manuscript summarizing the state of knowledge and parameterization of colloid-facilitated transport mechanisms in support of reactive transport and performance assessment models for generic crystalline repositories. This draft manuscript is being submitted as a level 3 milestone with LANL as the primary author. LLNL’s contribution to that effort is summarized only briefly in the present report. A manuscript summarizing long-term U(VI) diffusion experiments through bentonite backfill material was recently accepted for publication; the contents of that manuscript are summarized in the present report. The Np(IV) diffusion experiments were started mid-year and are ongoing. The completion of these experiments is planned for early FY17. Our progress in quantifying Np(IV) diffusion in bentonite backfill is summarized in the present report. Our involvement with the NEA TDB project was summarized in a recent Argillite Disposal activity report. It is not included in this report.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Shallow Cumulus in WRF Parameterizations Evaluated against LASSO Large-Eddy Simulations

Abstract Representation of shallow cumulus is a challenge for mesoscale numerical weather prediction models. These cloud fields have important effects on temperature, solar irradiance, convective initiation, and pollutant transport, among other processes. Recent improvements to physics schemes available in the Weather Research and Forecasting (WRF) Model aim to improve representation of shallow cumulus, in particular over land. The DOE LES ARM Symbiotic Simulation and Observation Workflow (LASSO) project provides several cases that we use here to test the new physics improvements. The LASSO cases use multiple large-scale forcings to drive large-eddy simulations (LES), and the LES output is easily compared to output from WRF single-column simulations driven with the same initial conditions and forcings. The new Mellor–Yamada–Nakanishi–Niino (MYNN) eddy diffusivity mass-flux (EDMF) boundary layer and shallow cloud scheme produces clouds with timing, liquid water path (LWP), and cloud fraction that agree well with LES over a wide range of those variables. Here we examine those variables and test the scheme’s sensitivity to perturbations of a few key parameters. We also discuss the challenges and uncertainties of single-column tests. The older, simpler total energy mass-flux (TEMF) scheme is included for comparison, and its tuning is improved. This is the first published use of the LASSO cases for parameterization development, and the first published study to use such a large number of cases with varying cloud amount. This is also the first study to use a more precise combined infrared and microwave retrieval of LWP to evaluate modeled clouds.

54 ENVIRONMENTAL SCIENCES↗

Modeling Extreme Deformations in Lithium Ion Batteries

A simultaneously coupled modeling approach to study the electrochemical and thermal behavior of lithium-ion batteries under large mechanical deformation has been developed. The thermo-electrochemical pseudo-2D (P2D) battery model is coupled with a mechanical material model. Mechanical, thermal, and electrochemical models are implemented as user-defined sub-routines in the commercial multi-physics code LS-DYNA. The mechanical strain experienced by anode, cathode and separator results in thickness and porosity changes in each layer which in turn influences electrochemical behavior. The evolution of concentration profiles and cell potential are studied under different mechanical loading conditions. Internal short-circuits caused by mechanical deformation and corresponding physical behaviors are also elucidated. We discuss the competing effects of improved transport at higher temperature due to the internal short-circuit versus a drop in the effective ionic conductivity and electrolyte diffusivity due to mechanical deformation.

47 OTHER INSTRUMENTATION↗

Metal oxide candidates for thermochemical water splitting obtained with a generative diffusion model

Generative diffusion models (DMs) for inorganic crystalline materials are being actively investigated for their potential to expand the chemical and structural design spaces for known functional materials. Generative candidates are particularly useful for applications where few functional, let alone commercially viable, materials currently exist, such as metal oxides for thermochemical water-splitting, which have strict requirements for defect thermodynamics and host stability. Here, we critically examine generated metal oxides from the M ATTER G EN DM conditioned on select chemical systems for thermochemical water splitting applications. Perhaps most notably, we find that M ATTER G EN predicts a novel, thermodynamically stable, quinary metal oxide, Ba 2 SrInFeO 6 , although this compound represents an ordered and layered substitution within the same A 3 B 2 O 6 structural prototype as its two ternary end members. Detailed density functional theory calculations and spin configuration sampling for this material and its possible decomposition products—beyond what existed in M ATTER G EN training data—are required to quantitatively validate hull energy predictions and conclusions of stability. Furthermore, the material exhibits oxygen defect formation energies appropriate for thermochemical water splitting, warranting targeted investigation in an experimental validation campaign, along with other future M ATTER G EN candidates in this application space.

36 MATERIALS SCIENCE↗

Comparison of mercury (Hg) bioaccumulation with mono- and mixed Lemna minor and Spirodela polyrhiza cultures

Mercury (Hg) is a prevalent and harmful contaminant that persists in the environment. For phytoremediation, it is important to discover which plants can bioaccumulate meaningful amounts of Hg while also tolerating its toxicity. Additionally, increasing biodiversity could create a more resilient and self-sustaining system for remediation. This study explores whether mixed populations of Lemna minor and Spirodela polyrhiza can better bioaccumulate and tolerate Hg than monocultures. Mono- and mixed cultures of L. minor and S. polyrhiza were grown in mesocosms of 0.5 μg/L or 100 μg/L Hg ( HgCl2) spiked water for 96 h. Change in weight of duckweed was used to assess Hg tolerance. Diffusive gradients in thin-films (DGTs) were used as surrogate monitoring devices for bioavailable levels of Hg. For biomass growth, the mixed culture of the L. minor was greater than the monoculture at the high dose. The L. minor accumulated more Hg in the mixed culture at the low dose while the S. polyrhiza was higher in the mixed at the high dose. Hg speciation in water was modeled using Windermere Humic Aqueous Model 7 (WHAM7) to compare the bioavailable species indicated by the DGTs. Potentially due to the controlled conditions, the WHAM7 output of bioavailable Hg was almost 1:1 to that estimated by the DGTs, indicating good predictive capability of geochemical modeling and passive sampler DGT on metal bioavailability. Altogether, the mixed cultures statistically performed as well as or better than the monocultures when tolerating and bioaccumulating Hg. However, there needs to be further work to see if the significant differences translate into practical differences worth the extra resources to maintain multiple species.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗

Extension of the DIF3D Code for Molten Salt Reactor Analysis

The molten salt reactor (MSR) with the fuel dissolved into the liquid salt has been selected as one of the six Generation-IV reactors due to its excellent characteristics in the aspects of sustainability, economy, passive safety, and resource utilization. In a MSR, the fuel salt is circulated throughout the whole system and thus some of delayed neutron precursors generated in the core decay in position of low importance or even out of the core. As a result, the multiplication factor (k-eff) of MSR is smaller than that of a reactor of stationary fuel when all the other conditions are same. Therefore, the neutronics code developed for stationary fuel cannot be applied directly to MSR analysis. In this paper, we describe the extension and application of the DIF3D code to MSR neutronics analysis. The DIF3D code system is a multi-group steady-state neutron diffusion and transport theory solver. It provides three flux solution options: finite-difference diffusion theory, nodal diffusion theory, and variational nodal transport theory methods. Eigenvalue, adjoint, fixed source and criticality (concentration) search problems are permitted. The DIF3D code has recently been modified to model the flowing fuel. To verify the new capability of DIF3D, the MSFR [4] core was analyzed with a wide range of velocity field of fuel salt and time spent outside of the core.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Phase-field modeling of diffusional phase behaviors of solid surfaces: A case study of phase-separating Li X FePO 4 electrode particles

In this paper, we present a comprehensive phase-field model for simulating diffusion-mediated kinetic phase behaviors near the surface of a solid particle. The model incorporates elastic inhomogeneity and anisotropy, diffusion mobility anisotropy, interfacial energy anisotropy, and Cahn–Hilliard diffusion kinetics. The free energy density function is formulated based on the regular solution model taking into account the possible solute-surface interaction near the surface. The coherency strain energy is computed using the Fourier-spectral iterative-perturbation method due to the strong elastic inhomogeneity with a zero surface traction boundary condition. Employing a phase-separating Li X FePO 4 electrode particle for Li-ion batteries as a model system, we perform parametric three-dimensional computer simulations. The model permits the observation of surface phase behaviors that are different from the bulk counterpart. For instance, it reproduces the theoretically well-established surface modes of spinodal decomposition of an unstable solid solution: the surface mode of coherent spinodal decomposition and the surface-directed spinodal decomposition mode. We systematically investigate the influences of major factors on the kinetic surface phase behaviors during the diffusional process. Finally, our simulation study provides insights for tailoring the internal phase microstructure of a particle by controlling the surface phase morphology.

25 ENERGY STORAGE↗

Comparing the structures and photophysical properties of two charge transfer co-crystals

Organic co-crystals have emerged as a promising class of semiconductors for next-generation optoelectronic devices due to their unique photophysical properties. This paper presents a joint experimental-theoretical study comparing the crystal structure, spectroscopy, and electronic structure of two charge transfer co-crystals. Reported herein is a novel co-crystal Npe:TCNQ, formed from 4-(1-naphthylvinyl)pyridine (Npe) and 7,7,8,8-tetracyanoquinodimethane (TCNQ) via molecular self-assembly. This work also presents a revised study of the co-crystal composed of Npe and 1,2,4,5-tetracyanobenzene (TCNB) molecules, Npe:TCNB, herein reported with a higher-symmetry (monoclinic) crystal structure than previously published. Npe:TCNB and Npe:TCNQ dimer clusters are used as theoretical model systems for the co-crystals; the geometries of the dimers are compared to geometries of the extended solids, which are computed with periodic boundary conditions density functional theory. UV-Vis absorption spectra of the dimers are computed with time-dependent density functional theory and compared to experimental UV-Vis diffuse reflectance spectra. Both Npe:TCNB and Npe:TCNQ are found to exhibit neutral character in the S 0 state and ionic character in the S 1 state. In conclusion, the high degree of charge transfer in the S 1 state of both Npe:TCNB and Npe:TCNQ is rationalized by analyzing the changes in orbital localization associated with the S 1 transitions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Liquid chloroform structure from computer simulation with a full ab initio intermolecular interaction potential

We have calculated the intermolecular interaction energies of the chloroform dimer in 12 orientations using the second-order Møller-Plesset perturbation theory. Single point energies of important geometries were calibrated by the coupled cluster with single and double and perturbative triple excitation method. Dunning's correlation consistent basis sets up to aug-cc-pVQZ have been employed in extrapolating the interaction energies to the complete basis set limit values. With the ab initio potential data we constructed a 5-site force field model for molecular dynamics simulations. We compared the simulation results with recent experiments and obtained quantitative agreements for the detailed atomwise radial distribution functions. Our results were also consistent with previous results using empirical force fields with polarization effects. Moreover, the calculated diffusion coefficients reproduced the experimental data over a wide range of thermodynamic conditions. To the best of our knowledge, this is the first ab initio force field which is capable of competing with existing empirical force fields for liquid chloroform.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

THE CIRCUMGALACTIC MEDIUM OF MASSIVE GALAXIES AT z {approx} 3: A TEST FOR STELLAR FEEDBACK, GALACTIC OUTFLOWS, AND COLD STREAMS

We present new results on the kinematics, thermal and ionization state, and spatial distribution of metal-enriched gas in the circumgalactic medium (CGM) of massive galaxies at redshift {approx}3, using the Eris suite of cosmological hydrodynamic ''zoom-in'' simulations. The reference run adopts a blastwave scheme for supernova feedback that produces large-scale galactic outflows, a star formation recipe based on a high gas density threshold, metal-dependent radiative cooling, and a model for the diffusion of metals and thermal energy. The effect of the local UV radiation field is added in post-processing. The CGM (defined as all gas at R > 0.2 R{sub vir} = 10 kpc, where R{sub vir} is the virial radius) contains multiple phases having a wide range of physical conditions, with more than half of its heavy elements locked in a warm-hot component at T > 10{sup 5} K. Synthetic spectra, generated by drawing sightlines through the CGM, produce interstellar absorption-line strengths of Ly{alpha}, C II, C IV, Si II, and Si IV as a function of the galactocentric impact parameter (scaled to the virial radius) that are in broad agreement with those observed at high redshift by Steidel et al. The covering factor of absorbing material declines less rapidly with impact parameter for Ly{alpha} and C IV compared to C II, Si IV, and Si II, with Ly{alpha} remaining strong (W{sub Ly{alpha}} > 300 mA) to {approx}> 5 R{sub vir} = 250 kpc. Only about one third of all the gas within R{sub vir} is outflowing. The fraction of sightlines within one virial radius that intercept optically thick, N{sub H{sub I}}>10{sup 17.2} cm{sup -2} material is 27%, in agreement with recent observations by Rudie et al. Such optically thick absorption is shown to trace inflowing ''cold'' streams that penetrate deep inside the virial radius. The streams, enriched to metallicities above 0.01 solar by previous episodes of star formation in the main host and in nearby dwarfs, are the origin of strong (N{sub C{sub II}}>10{sup 13} cm{sup -2}) C II absorption with a covering factor of 22% within R{sub vir} and 10% within 2 R{sub vir}. Galactic outflows do not cause any substantial suppression of the cold accretion mode. The central galaxy is surrounded by a large O VI halo, with a typical column density N{sub O{sub VI}} {approx}> 10{sup 14} cm{sup -2} and a near unity covering factor maintained all the way out to 150 kpc. This matches the trends recently observed in star-forming galaxies at low redshift by Tumlinson et al. Our zoom-in simulations of this single system appear then to reproduce quantitatively the complex baryonic processes that determine the exchange of matter, energy, and metals between galaxies and their surroundings.

79 ASTRONOMY AND ASTROPHYSICS↗

Fundamental Study of Disposition and Release of Methane in a Shale Gas Reservoir

The recent boom in shale gas production through hydrofracturing has reshaped the energy production landscape in the United States. Wellbore production rates vary greatly among the wells within a single field and decline rapidly with time, thus bring up a serious concern with the sustainability of shale gas production. Shale gas production starts with creating a fracture network by injecting a pressurized fluid in a wellbore. The induced fractures are then held open by proppant particles. During production, gas releases from the mudstone matrix, migrates to nearby fractures, and ultimately reaches a production wellbore. Given the relatively high permeability of the induced fractures, gas release and migration in low-permeability shale matrix is likely to be a limiting step for long-term wellbore production. Therefore, a clear understanding of the underlying mechanisms of methane disposition and release in shale matrix is crucial for the development of new technologies to maximize gas production and recovery. Shale is a natural nanocomposite material with distinct characteristics of nanometer-scale pore sizes, extremely low permeability, high clay contents, significant amounts of organic carbon, and large spatial heterogeneities. Our work has shown that nanopore confinement plays an important role in methane disposition and release in shale matrix. Using molecular simulations, we show that methane release in nanoporous kerogen matrix is characterized by fast release of pressurized free gas (accounting for ~ 30 - 47% recovery) followed by slow release of adsorbed gas as the gas pressure decreases. The first stage is driven by the gas pressure gradient while the second stage is controlled by gas desorption and diffusion. The long-term production decline appears controlled by the second stage of gas release. We further show that diffusion of all methane in nanoporous kerogen behaves differently from the bulk phase, with much smaller diffusion coefficients. The MD simulations also indicate that a significant fraction (3 - 35%) of methane deposited in kerogen can potentially become trapped in isolated nanopores and thus not recoverable. We have successfully established experimental capabilities for measuring gas sorption and desorption on shale and model materials under a wide range of physical and chemical conditions. Both low and high pressure measurements show significant sorption of CH 4 and CO 2 onto clays, implying that methane adsorbed on clay minerals could contribute a significant portion of gas-in-place in an unconventional reservoir. We have also studied the potential impact of the interaction of shale with hydrofracking fluid on gas sorption. We have found that the CH 4 -CO 2 sorption capacity for the reacted sample is systematically lower (by a factor of ~2) than that for the unreacted (raw) sample. This difference in sorption capacity may result from a mineralogical or surface chemistry change of the shale sample induced by fluid-rock interaction. Our results shed a new light on mechanistic understanding gas release and production decline in unconventional reservoirs.

04 OIL SHALES AND TAR SANDS↗

A Phase-Field Model for In-Space Manufacturing of Binary Alloys

The integrity of the final printed components is mostly dictated by the adhesion between the particles and phases that form upon solidification, which is a major problem in printing metallic parts using available In-Space Manufacturing (ISM) technologies based on the Fused Deposition Modeling (FDM) methodology. Understanding the melting/solidification process helps increase particle adherence and allows to produce components with greater mechanical integrity. We developed a phase-field model of solidification for binary alloys. The phase-field approach is unique in capturing the microstructure with computationally tractable costs. The developed phase-field model of solidification of binary alloys satisfies the stability conditions at all temperatures. The suggested model is tuned for Ni-Cu alloy feedstocks. We derived the Ginzburg-Landau equations governing the phase transformation kinetics and solved them analytically for the dilute solution. We calculated the concentration profile as a function of interface velocity for a one-dimensional steady-state diffuse interface neglecting elasticity and obtained the partition coefficient, k, as a function of interface velocity. Numerical simulations for the diluted solution are used to study the interface velocity as a function of undercooling for the classic sharp interface model, partitionless solidification, and thin interface.

36 MATERIALS SCIENCE↗

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↗