This document summarizes research performed under the SNL LDRD entitled - Computational Mechanics for Geosystems Management to Support the Energy and Natural Resources Mission. The main accomplishment was development of a foundational SNL capability for computational thermal, chemical, fluid, and solid mechanics analysis of geosystems. The code was developed within the SNL Sierra software system. This report summarizes the capabilities of the simulation code and the supporting research and development conducted under this LDRD. The main goal of this project was the development of a foundational capability for coupled thermal, hydrological, mechanical, chemical (THMC) simulation of heterogeneous geosystems utilizing massively parallel processing. To solve these complex issues, this project integrated research in numerical mathematics and algorithms for chemically reactive multiphase systems with computer science research in adaptive coupled solution control and framework architecture. This report summarizes and demonstrates the capabilities that were developed together with the supporting research underlying the models. Key accomplishments are: (1) General capability for modeling nonisothermal, multiphase, multicomponent flow in heterogeneous porous geologic materials; (2) General capability to model multiphase reactive transport of species in heterogeneous porous media; (3) Constitutive models for describing real, general geomaterials under multiphase conditions utilizing laboratory data; (4) General capability to couple nonisothermal reactive flow with geomechanics (THMC); (5) Phase behavior thermodynamics for the CO2-H2O-NaCl system. General implementation enables modeling of other fluid mixtures. Adaptive look-up tables enable thermodynamic capability to other simulators; (6) Capability for statistical modeling of heterogeneity in geologic materials; and (7) Simulator utilizes unstructured grids on parallel processing computers.
The PUREX process is one of the most widely adopted approaches of industrial-scale liquid-liquid extraction (LLE) methods for the separation and recovery of U and Pu from dissolved used nuclear fuel. This paper documents the application of chemical equilibrium approaches to the modeling of LLE for the PUREX process at 25°C. In this work, we focus on modeling the extraction of HNO3 and UO2(NO3)2 by the non-electrolyte organic liquid phase Tri-n-butyl Phosphate (TBP) plus a diluent (Amsco 125-82). Six chemical reactions representing the equilibria between the concentrated HNO3-UO2(NO3)2 electrolyte and TBP-diluent organic phase is sufficient to accurately generate the extraction isotherms for the extracted components for a wide range of TBP and acid concentrations in the non-electrolyte and electrolyte phases, respectively. The Pitzer approach is used to compute the activity coefficients of the HNO3 electrolyte phase. No activity coefficient model is adopted for TBP organic complexes where these were assigned unity values. Even with this assumption, the predicted isotherms are in very good agreement with the experimental data and such result serves as a preamble to extend this modeling approach to more compositionally complex systems relevant to LLE of radionuclides.
The U.S. is currently re-evaluating the policy on high-level radioactive waste (HLW) management and has been studying generic disposal system environment (GDSE) concepts to support the development of a long-term strategy for geologic disposal of HLW. The GDSE study focuses on the analysis of different GDSE options, and a salt repository is one of the options currently under study. The immediate goal of the generic salt repository study is to develop the necessary modeling tools to evaluate and improve understanding on the repository system response and processes relevant to long-term HLW disposal in salt. An initial version of the salt GDSE performance assessment model and the preliminary analysis results are discussed, emphasizing key attributes of a salt repository that are potentially important to the long-term safe disposal of HLW. Also discussed are the preliminary results on the repository response to the effects of different waste types (commercial UNF, existing DOE HLW, and hypothetical reprocessing HLW), and radionuclide release scenarios (undisturbed and human intrusion). Soluble, non- to weakly sorbing fission products, particularly 129I, 79Se, and 26Ra are the major dose contributors. However, the conservative assumptions made about their geochemical behaviors contribute to their calculated dose. The paper elaborates on the identified knowledge gaps and path forwards for future R&D efforts to advance understanding of salt repository system performance for HLW disposal.
Key components of the nuclear fuel cycle are short-term storage and long-term disposal of nuclear waste. The latter encompasses the immobilization of used nuclear fuel (UNF) and radioactive waste streams generated by various phases of the nuclear fuel cycle, and the safe and permanent disposition of these waste forms in geological repository environments. The engineered barrier system (EBS) plays a very important role in the long-term isolation of nuclear waste in geological repository environments. EBS concepts and their interactions with the natural barrier are inherently important to the long-term performance assessment of the safety case where nuclear waste disposition needs to be evaluated for time periods of up to one million years. Making the safety case needed in the decision-making process for the recommendation and the eventual embracement of a disposal system concept requires a multi-faceted integration of knowledge and evidence-gathering to demonstrate the required confidence level in a deep geological disposal site and to evaluate long-term repository performance. The focus of this report is the following: (1) Evaluation of EBS in long-term disposal systems in deep geologic environments with emphasis on the multi-barrier concept; (2) Evaluation of key parameters in the characterization of EBS performance; (3) Identification of key knowledge gaps and uncertainties; and (4) Evaluation of tools and modeling approaches for EBS processes and performance. The above topics will be evaluated through the analysis of the following: (1) Overview of EBS concepts for various NW disposal systems; (2) Natural and man-made analogs, room chemistry, hydrochemistry of deep subsurface environments, and EBS material stability in near-field environments; (3) Reactive Transport and Coupled Thermal-Hydrological-Mechanical-Chemical (THMC) processes in EBS; and (4) Thermal analysis toolkit, metallic barrier degradation mode survey, and development of a Disposal Systems Evaluation Framework (DSEF). This report will focus on the multi-barrier concept of EBS and variants of this type which in essence is the most adopted concept by various repository programs. Empasis is given mainly to the evaluation of EBS materials and processes through the analysis of published studies in the scientific literature of past and existing repository research programs. Tool evaluations are also emphasized, particularly on THCM processes and chemical equilibria. Although being an increasingly important aspect of NW disposition, short-term or interim storage of NW will be briefly discussed but not to the extent of the EBS issues relevant to disposal systems in deep geologic environments. Interim storage will be discussed in the report Evaluation of Storage Concepts FY10 Final Report (Weiner et al. 2010).
Imported oil exacerabates our trade deficit and funds anti-American regimes. Nuclear Energy (NE) is a demonstrated technology with high efficiency. NE's two biggest political detriments are possible accidents and nuclear waste disposal. For NE policy, proliferation is the biggest obstacle. Nuclear waste can be reduced through reprocessing, where fuel rods are separated into various streams, some of which can be reused in reactors. Current process developed in the 1950s is dirty and expensive, U/Pu separation is the most critical. Fuel rods are sheared and dissolved in acid to extract fissile material in a centrifugal contactor. Plants have many contacts in series with other separations. We have taken a science and simulation-based approach to develop a modern reprocessing plant. Models of reprocessing plants are needed to support nuclear materials accountancy, nonproliferation, plant design, and plant scale-up.
Carbon capture and sequestration (CCS) is an option to mitigate impacts of atmospheric carbon emission. Numerous factors are important in determining the overall effectiveness of long-term geologic storage of carbon, including leakage rates, volume of storage available, and system costs. Recent efforts have been made to apply an existing probabilistic performance assessment (PA) methodology developed for deep nuclear waste geologic repositories to evaluate the effectiveness of subsurface carbon storage (Viswanathan et al., 2008; Stauffer et al., 2009). However, to address the most pressing management, regulatory, and scientific concerns with subsurface carbon storage (CS), the existing PA methodology and tools must be enhanced and upgraded. For example, in the evaluation of a nuclear waste repository, a PA model is essentially a forward model that samples input parameters and runs multiple realizations to estimate future consequences and determine important parameters driving the system performance. In the CS evaluation, however, a PA model must be able to run both forward and inverse calculations to support optimization of CO{sub 2} injection and real-time site monitoring as an integral part of the system design and operation. The monitoring data must be continually fused into the PA model through model inversion and parameter estimation. Model calculations will in turn guide the design of optimal monitoring and carbon-injection strategies (e.g., in terms of monitoring techniques, locations, and time intervals). Under the support of Laboratory-Directed Research & Development (LDRD), a late-start LDRD project was initiated in June of Fiscal Year 2010 to explore the concept of an enhanced performance assessment system (EPAS) for carbon sequestration and storage. In spite of the tight time constraints, significant progress has been made on the project: (1) Following the general PA methodology, a preliminary Feature, Event, and Process (FEP) analysis was performed for a hypothetical CS system. Through this FEP analysis, relevant scenarios for CO{sub 2} release were defined. (2) A prototype of EPAS was developed by wrapping an existing multi-phase, multi-component reservoir simulator (TOUGH2) with an uncertainty quantification and optimization code (DAKOTA). (3) For demonstration, a probabilistic PA analysis was successfully performed for a hypothetical CS system based on an existing project in a brine-bearing sandstone. The work lays the foundation for the development of a new generation of PA tools for effective management of CS activities. At a top-level, the work supports energy security and climate change/adaptation by furthering the capability to effectively manage proposed carbon capture and sequestration activities (both research and development as well as operational), and it greatly enhances the technical capability to address this national problem. The next phase of the work will include (1) full capability demonstration of the EPAS, especially for data fusion, carbon storage system optimization, and process optimization of CO{sub 2} injection, and (2) application of the EPAS to actual carbon storage systems.
In this report, we summarize our work on developing a production level capability for modeling brine thermodynamic properties using the open-source code Cantera. This implementation into Cantera allows for the application of chemical thermodynamics to describe the interactions between a solid and an electrolyte solution at chemical equilibrium. The formulations to evaluate the thermodynamic properties of electrolytes are based on Pitzer's model to calculate molality-based activity coefficients using a real equation-of-state (EoS) for water. In addition, the thermodynamic properties of solutes at elevated temperature and pressures are computed using the revised Helgeson-Kirkham-Flowers (HKF) EoS for ionic and neutral aqueous species. The thermodynamic data parameters for the Pitzer formulation and HKF EoS are from the thermodynamic database compilation developed for the Yucca Mountain Project (YMP) used with the computer code EQ3/6. We describe the adopted equations and their implementation within Cantera and also provide several validated examples relevant to the calculations of extensive properties of electrolyte solutions.
We present a model to evaluate the water mass balance inside a breached waste package in Yucca Mountain (YM) repository environments. The amount of water as liquid or vapor that can accumulate inside or percolate through the package in the emplacement drift is modeled as a function of the temperature and relative humidity (RH) near the waste package, the dripping rate of water from seepage, the area of failure patches on the waste package, and the extent of waste degradation. The water activity inside the waste package is assumed to be determined by both matric and osmotic potentials in the porous waste degradation products that also includes hygroscopic salts. We implemented the model and conducted a set of Monte Carlo simulations to gain insight into the variability and uncertainty associated with model predictions. The model shows that water vapor diffusion can be as important as the advective seepage flow. In addition, chemical reactions during waste degradation can consume a significant fraction of water accumulated in the waste package.