ANS IHLRWM 2017 - 16th International High-Level Radioactive Waste Management Conference: Creating a Safe and Secure Energy Future for Generations to Come - Driving Toward Long-Term Storage and Disposal
Numerical simulation of a repository for heatgenerating nuclear waste in fractured crystalline rock requires a method for simulating coupled heat and fluid flow and reactive radionuclide transport in both porous media (bentonite buffer, surface sediments) and fractured rock (the repository host rock). Discrete fracture networks (DFNs), networks of two-dimensional planes distributed in a three-dimensional domain, are commonly used to simulate isothermal fluid flow and particle transport in fractures, but unless coupled to a continuum, are incapable of simulating heat conduction through the rock matrix, and therefore incapable of capturing the effects of thermally driven fluid fluxes or of coupling chemical processes to thermal processes. We present a method for mapping a stochastically generated DFN to a porous medium domain that allows representation of porous and fractured media in the same domain, captures the behavior of radionuclide transport in fractured rock, and allows simulation of coupled heat and fluid flow including heat conduction through the matrix of the fractured rock. We apply the method within Sandia's Geologic Disposal Safety Assessment (GDSA) framework to conduct a post-closure performance assessment (PA) of a generic repository for commercial spent nuclear fuel in crystalline rock. The three-dimensional, kilometer-scale model domain contains approximately 4.5 million grid cells; grid refinement captures the detail of 3, 360 individual waste packages in 42 disposal drifts. Coupled heat and fluid flow and reactive transport are solved numerically with PFLOTRAN, a massively parallel multiphase flow and reactive transport code. Simulations of multiple fracture realizations were run to 1 million years, and indicate that, because of the channeled nature of fracture flow, thermally-driven fluid fluxes associated with peak repository temperatures may be a primary means of radionuclide transport out of the saturated repository. The channeled nature of fracture flow gives rise to unique challenges in uncertainty and sensitivity quantification, as radionuclide concentrations at any given location outside the repository depend heavily on the distribution of fractures in the domain.
Extreme-scale computational science increasingly demands multiscale and multiphysics formulations. Combining software developed by independent groups is imperative: no single team has resources for all predictive science and decision support capabilities. Scientific libraries provide high-quality, reusable software components for constructing applications with improved robustness and portability. However, without coordination, many libraries cannot be easily composed. Namespace collisions, inconsistent arguments, lack of third-party software versioning, and additional difficulties make composition costly. The Extreme-scale Scientific Software Development Kit (xSDK) defines community policies to improve code quality and compatibility across independently developed packages (hypre, PETSc, SuperLU, Trilinos, and Alquimia) and provides a foundation for addressing broader issues in software interoperability, performance portability, and sustainability. The xSDK provides turnkey installation of member software and seamless combination of aggregate capabilities, and it marks first steps toward extreme-scale scientific software ecosystems from which future applications can be composed rapidly with assured quality and scalability.
The geological disposal of nuclear waste is based on the multi-barrier concept, comprising various engineered and natural barriers, to confine the radioactive waste and isolate it from the biosphere. Some of the planned repositories for high-level nuclear waste will be hosted in fractured crystalline rock formations. The potential of these formations to act as natural transport barriers is related to two coupled processes: diffusion into the rock matrix and sorption onto the mineral surfaces available in the rock matrix. Different in situ and laboratory experiments have pointed out the ubiquitous heterogeneous nature of the rock matrix: mineral surfaces and pore space are distributed in complex microstructures and their distribution is far from being homogeneous (as typically assumed by Darcy-scale coarse reactive transport models). In this work, we use a synthetically generated fracture–matrix system to assess the implications of grain-scale physical and mineralogical heterogeneity on cesium transport and retention. The resulting grain-scale reactive transport model is solved using high-performance computing technologies, and the results are compared with those derived from two alternative models, denoted as upscaled models, where mineral abundance is averaged over the matrix volume. In the grain-scale model, the penetration of cesium into the matrix is faster and the penetration front is uneven and finger-shaped. The analysis of the cesium breakthrough curves computed at two different points in the fracture shows that the upscaled models provide later first-arrival time estimates compared to the grain-scale model. The breakthrough curves computed with the three models converge at late times. These results suggest that spatially averaged upscaled parameters of sorption site distribution can be used to predict the late-time behavior of breakthrough curves but could be inadequate to simulate the early behavior.
The Spent Fuel and Waste Science and Technology (SFWST) Campaign of the U.S. Department of Energy (DOE) Office of Nuclear Energy (NE) is conducting research and development (R&D) on generic deep geologic disposal systems (i.e., repositories). This report describes specific activities in FY 2016 associated with the development of a Defense Waste Repository (DWR)a for the permanent disposal of a portion of the HLW and SNF derived from national defense and research and development (R&D) activities of the DOE.
The Used Fuel Disposition Campaign (UFDC) of the U.S. Department of Energy (DOE) Office of Nuclear Energy (NE), Office of Fuel Cycle Technology (OFCT) is conducting research and development (R&D) on geologic disposal of used nuclear fuel (UNF) and high-level nuclear waste (HLW). Two of the high priorities for UFDC disposal R&D are design concept development and disposal system modeling (DOE 2011). These priorities are directly addressed in the UFDC Generic Disposal Systems Analysis (GDSA) work package, which is charged with developing a disposal system modeling and analysis capability for evaluating disposal system performance for nuclear waste in geologic media (e.g., salt, granite, clay, and deep borehole disposal). This report describes specific GDSA activities in fiscal year 2016 (FY 2016) toward the development of the enhanced disposal system modeling and analysis capability for geologic disposal of nuclear waste. The GDSA framework employs the PFLOTRAN thermal-hydrologic-chemical multi-physics code and the Dakota uncertainty sampling and propagation code. Each code is designed for massively-parallel processing in a high-performance computing (HPC) environment. Multi-physics representations in PFLOTRAN are used to simulate various coupled processes including heat flow, fluid flow, waste dissolution, radionuclide release, radionuclide decay and ingrowth, precipitation and dissolution of secondary phases, and radionuclide transport through engineered barriers and natural geologic barriers to the biosphere. Dakota is used to generate sets of representative realizations and to analyze parameter sensitivity.
The R&D program from the DOE Used Fuel Disposition Campaign (UFDC) has documented key advances in coupled Thermal-Hydrological-Mechanical-Chemical (THMC) modeling of clay to simulate its complex dynamic behavior in response to thermal and hydrochemical feedbacks. These efforts have been harnessed to assess the isolation performance of heat-generating nuclear waste in a deep geological repository in clay/shale/argillaceous rock formations. This report describes the ongoing disposal R&D efforts on the advancement and refinement of coupled THMC process models, hydrothermal experiments on barrier clay interactions, used fuel and canister material degradation, thermodynamic database development, and reactive transport modeling of the near-field under non-isothermal conditions. These play an important role to the evaluation of sacrificial zones as part of the EBS exposure to thermally-driven chemical and transport processes. Thermal inducement of chemical interactions at EBS domains enhances mineral dissolution/precipitation but also generates mineralogical changes that result in mineral H2O uptake/removal (hydration/dehydration reactions). These processes can result in volume changes that can affect the interface / bulk phase porosities and the mechanical (stress) state of the bentonite barrier. Characterization studies on bentonite barrier samples from the FEBEX-DP international activity have provided important insight on clay barrier microstructures (e.g., microcracks) and interactions at EBS interfaces. Enhancements to the used fuel degradation model outlines the need to include the effects of canister corrosion due the strong influence of H2 generation on the source term.