Publications

Abstract not provided.

This report summarizes the international collaborations conducted by Sandia funded by the US Department of Energy Office (DOE) of Nuclear Energy Spent Fuel and Waste Science & Technology (SFWST) as part of the Sandia National Laboratories Salt R&D and Salt International work packages. This report satisfies the level-three milestone M3SF-22SN010303063. Several stand-alone sections make up this summary report, each completed by the participants. The sections discuss international collaborations on geomechanical benchmarking exercises (WEIMOS), granular salt reconsolidation (KOMPASS), engineered barriers (RANGERS), numerical model comparison (DECOVALEX) and an NEA Salt Club working group on the development of scenarios as part of the performance assessment development process. Finally, we summarize events related to the US/German Workshop on Repository Research, Design and Operations. The work summarized in this annual update has occurred during the COVID-19 pandemic, and little international or domestic travel has occurred. Most of the collaborations have been conducted via email or as virtual meetings, but a slow return to travel and in-person meetings has begun.

Estimation of two-phase fluid flow properties is important to understand and predict water and gas movement through the vadose zone for agricultural, hydrogeological, and engineering applications, such as containment transport and/or containment of gases in the subsurface. To estimate rock fluid flow properties and subsequently predict physically realistic processes such as patterns and timing of water, gas, and energy (e.g., heat) movement in the subsurface, laboratory spontaneous water imbibition with simultaneous temperature measurement and numerical modeling methods are presented in the FY22 progress report. A multiple-overlapping-continua conceptual model is used to explain and predict observed complex multi-phenomenological laboratory test behavior during spontaneous imbibition experiments. This report primarily addresses two complexities that arise during the experiments: 1) capturing the late-time behavior of spontaneous imbibition tests with dual porosity; and 2) understanding the thermal perturbation observed at or ahead of the imbibing wetting front, which are associated with adsorption of water in initially dry samples. We use numerical approaches to explore some of these issues, but also lay out a plan for further laboratory experimentation and modeling to best understand and leverage these unique observations.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The capability of a 1-D PFLOTRAN model to simulate the S1-3 bentonite saturation experiment has been demonstrated and validated against experimental data. Work remains to be done to refine 1-D PFLOTRAN simulations of the experiment S1-4 which include evaluation of parameter sensitivities on the prediction of material saturation and relative permeabilities. This and further testing of PFLOTRAN capabilities will be done as part of DECOVALEX 2023 Task D contributions by the SNL team in the coming months.

Two-phase fluid flow properties underlie quantitative prediction of water and gas movement, but constraining these properties typically requires multiple time-consuming laboratory methods. The estimation of two-phase flow properties (van Genuchten parameters, porosity, and intrinsic permeability) is illustrated in cores of vitric nonwelded volcanic tuff using Bayesian parameter estimation that fits numerical models to observations from spontaneous imbibition experiments. The uniqueness and correlation of the estimated parameters is explored using different modeling assumptions and subsets of the observed data. The resulting estimation process is sensitive to both moisture retention and relative permeability functions, thereby offering a comprehensive method for constraining both functions. The data collected during this relatively simple laboratory experiment, used in conjunction with a numerical model and a global optimizer, result in a viable approach for augmenting more traditional capillary pressure data obtained from hanging water column, membrane plate extractor, or mercury intrusion methods. This method may be useful when imbibition rather than drainage parameters are sought, when larger samples (e.g., including heterogeneity or fractures) need to be tested that cannot be accommodated in more traditional methods, or when in educational laboratory settings.

Abstract not provided.

Abstract not provided.

Abstract not provided.

A crucial component of field testing is the utilization of numerical models to better understand the system and the experimental data being collected. Meshing and modeling field tests is a complex and computationally demanding problem. Hexahedral elements cannot always reproduce experimental dimensions leading to grid orientation or geometric errors. Voronoi meshes can match complex geometries without sacrificing orthogonality. As a result, here we present a high-resolution 3D numerical study for the BATS heater test at the WIPP that compares both a standard non-deformed cartesian mesh along with a Voronoi mesh to match field data collected during a salt heater experiment.

This paper summarizes the development of post-closure safety assessment for radioactive waste disposal from the point of view of scenarios, which occupy the key point in the process between FEPs and assessment using conceptual, mathematical, and numerical models. Scenarios are used in other fields for similar purposes, but they have a central role in safety assessment for radioactive waste disposal, given the large uncertainties in natural and engineered systems over long time periods. Repository design and assessments are built around a base scenario, which is usually built up from FEPs in a deductive bottom-up fashion. The alternative scenarios are often a perturbation of the base scenario, constructed in a top-down fashion around individual safety functions of key repository features. Despite differences between nations in how they implement scenarios, largely from regulatory differences, the concept of scenarios is beneficial and is used universally in development of deep geological repositories. The methodology has also seen some use outside the field radioactive waste disposal, but its wider adoption might be warranted.

Abstract not provided.

The construction of deep geological repositories (DGR) in salt formations requires penetrating through naturally sealing geosphere layers. While the emplaced nuclear waste is primarily protected by the containment-providing rock zone (CRZ), technical barriers are required, for example during handling. For closure geotechnical barriers seal the repository along the accesses against water or solutions from outside and the possible emission paths for radionuclides contained inside. As these barriers must ensure maintenance-free function on a long-term basis, they typically comprise a set of specialized elements with diversified functions that may be used redundantly. The effects of the individual elements are coordinated so that they are collectively referred to as the Engineered Barrier System (EBS).

Elevated temperature and pressure in the earth's subsurface alters the permeability of salt formations, due to changing properties of the salt-brine interface. Molecular dynamics (MD) simulations are used to investigate the mechanisms of temperature and pressure dependence of liquid-solid interfacial tensions of NaCl, KCl, and NaCl-KCl brines in contact with (100) salt surfaces. Salt-brine dihedral angles vary between 55 and 76° across the temperature (300-450 K) and pressure range (0-150 MPa) evaluated. Temperature-dependent brine composition results in elevated dihedral angles of 65-80°, which falls above the reported salt percolation threshold of 60°. Mixed NaCl-KCl brine compositions increased this effect. Elevated temperatures excluded dissolved Na+ ions from the interface, causing the strong temperature dependence of the liquid-solid interfacial tension and the resulting dihedral angle. Therefore, at higher temperature, pressure, and brine concentrations Na-Cl systems may underpredict the dihedral angle. Higher dihedral angles in more realistic mixed brine systems maintain low permeability of salt formations due to changes in the structure and energetics of the salt-brine interface.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.