skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Global Sensitivity and Data-Worth Analyses in iTOUGH2: User's Guide


This manual explains the use of local sensitivity analysis, the global Morris OAT and Sobol’ methods, and a related data-worth analysis as implemented in iTOUGH2. In addition to input specification and output formats, it includes some examples to show how to interpret results.

 [1];  [1]
  1. Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States). Earth Sciences Division; Univ. of California, Berkeley, CA (United States)
Publication Date:
Research Org.:
Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States)
Sponsoring Org.:
OSTI Identifier:
Report Number(s):
DOE Contract Number:
Resource Type:
Technical Report
Country of Publication:
United States

Citation Formats

Wainwright, Haruko Murakami, and Finsterle, Stefan. Global Sensitivity and Data-Worth Analyses in iTOUGH2: User's Guide. United States: N. p., 2016. Web. doi:10.2172/1274412.
Wainwright, Haruko Murakami, & Finsterle, Stefan. Global Sensitivity and Data-Worth Analyses in iTOUGH2: User's Guide. United States. doi:10.2172/1274412.
Wainwright, Haruko Murakami, and Finsterle, Stefan. Fri . "Global Sensitivity and Data-Worth Analyses in iTOUGH2: User's Guide". United States. doi:10.2172/1274412.
title = {Global Sensitivity and Data-Worth Analyses in iTOUGH2: User's Guide},
author = {Wainwright, Haruko Murakami and Finsterle, Stefan},
abstractNote = {This manual explains the use of local sensitivity analysis, the global Morris OAT and Sobol’ methods, and a related data-worth analysis as implemented in iTOUGH2. In addition to input specification and output formats, it includes some examples to show how to interpret results.},
doi = {10.2172/1274412},
journal = {},
number = ,
volume = ,
place = {United States},
year = {Fri Jul 15 00:00:00 EDT 2016},
month = {Fri Jul 15 00:00:00 EDT 2016}

Technical Report:

Save / Share:
  • This document describes the development and use of the Integrated Flow Code (IFC), a numerical code and related model to be used for the simulation of time-dependent, two-phase flow in the near field and geosphere of a gas-generating nuclear waste repository system located in an initially fully water-saturated claystone (Opalinus Clay) in Switzerland. The development of the code and model was supported by the Swiss National Cooperative for the Disposal of Radioactive Waste (Nagra), Wettingen, Switzerland. Gas generation (mainly H{sub 2}, but also CH{sub 4} and CO{sub 2}) may affect repository performance by (1) compromising the engineered barriers through excessivemore » pressure build-up, (2) displacing potentially contaminated pore water, (3) releasing radioactive gases (e.g., those containing {sup 14}C and {sup 3}H), (4) changing hydrogeologic properties of the engineered barrier system and the host rock, and (5) altering the groundwater flow field and thus radionuclide migration paths. The IFC aims at providing water and gas flow fields as the basis for the subsequent radionuclide transport simulations, which are performed by the radionuclide transport code (RTC). The IFC, RTC and a waste-dissolution and near-field transport model (STMAN) are part of the Integrated Radionuclide Release Code (IRRC), which integrates all safety-relevant features, events, and processes (FEPs). The IRRC is embedded into a Probabilistic Safety Assessment (PSA) computational tool that (1) evaluates alternative conceptual models, scenarios, and disruptive events, and (2) performs Monte-Carlo sampling to account for parametric uncertainties. The preliminary probabilistic safety assessment concept and the role of the IFC are visualized in Figure 1. The IFC was developed based on Nagra's PSA concept. Specifically, as many phenomena as possible are to be directly simulated using a (simplified) process model, which is at the core of the IRRC model. Uncertainty evaluation (scenario uncertainty, conceptualization uncertainty, parametric uncertainty) is handled by the outer shell of the PSA model; it is not further discussed in this report. Moreover, justifications for the inclusion or exclusion of FEPs as well as for certain simplifying assumptions are available or can be obtained using detailed process models and other supporting information. The IFC is both a numerical code and a model of a repository system. The numerical code is a modification of the multiphase, multicomponent simulator TOUGH2 (Pruess et al., 1999), as implemented within the iTOUGH2 (Finsterle, 2007abc) framework. The code modifications are mainly concerned with the implementation of relevant FEPs as outlined in Nagra (2007a, AN 07-115), as well as removal of processes and features that are not needed within the IFC; the modifications are summarized in Appendix A1. In addition, the IFC includes a model, i.e., a simplified representation of the repository system. Specifically, a computational grid was generated, which includes the emplacement tunnels for spent fuel, high-level wastes, as well as long-lived intermediate-level wastes. Moreover, the model represents engineered barriers (backfill, seals, plugs, etc.), various tunnels and other underground facilities, and includes a simplified representation of the geological structure, i.e., the host rock (including the excavation disturbed zone (EDZ) around the underground openings), confining units, local aquifers, and a highly transmissive zone. The IFC model was designed in close collaboration with Nagra. This report describes all functional requirements of the IFC and how they are implemented in the IFC. The input formats needed to invoke added modeling capabilities are documented. Finally, the IFC model grid is described, and results from a test simulation are presented.« less
  • The precursor of TOUGH2, TOUGH, was originally developed with non-hysteretic characteristic curves. Hysteretic capillary pressure functions were implemented in TOUGH in the late 1980s by Niemi and Bodvarsson (1988), and hysteretic capillary pressure and relative permeability functions were added to iTOUGH2 about ten years later by Finsterle et al. (1998). Recently, modifications were made to the iTOUGH2 hysteretic formulation to make it more robust and efficient (Doughty, 2007). Code development is still underway, with the ultimate goal being a hysteretic module that fits into the standard TOUGH2 (Pruess et al., 1991) framework. This document provides a user's guide for themore » most recent version of the hysteretic code, which runs within iTOUGH2 (Finsterle, 1999a,b,c). The current code differs only slightly from what was presented in Doughty (2007), hence that document provides the basic information on the processes being modeled and how they are conceptualized. This document focuses on a description of the user-specified parameters required to run hysteretic iTOUGH2. In the few instances where the conceptualization differs from that of Doughty (2007), the features described here are the current ones. Sample problems presented in this user's guide use the equation-of-state module ECO2N (Pruess, 2005). The components present in ECO2N are H{sub 2}O, NaCl, and CO{sub 2}. Two fluid phases and one solid phase are considered: an aqueous phase, which primarily consists of liquid H2O and may contain dissolved NaCl and CO{sub 2}; a supercritical phase which primarily consists of CO{sub 2}, but also includes a small amount of gaseous H{sub 2}O; and a solid phase consisting of precipitated NaCl. Details of the ECO2N formulation may be found in Pruess (2005). The aqueous phase is the wetting phase and is denoted ''liquid'', whereas the supercritical phase is the non-wetting phase and is denoted ''gas''. The hysteretic formalism may be applied to other iTOUGH2 equation-of-state modules, as long as the liquid phase is the wetting phase and the gas phase is the non-wetting phase.« less
  • Supercritical fluids exist near magmatic heat sources in geothermal reservoirs, and the high enthalpy fluid is becoming more desirable for energy production with advancing technology. In geothermal modeling, the roots of the geothermal systems are normally avoided but in order to accurately predict the thermal behavior when wells are drilled close to magmatic intrusions, it is necessary to incorporate the heat sources into the modeling scheme. Modeling supercritical conditions poses a variety of challenges due to the large gradients in fluid properties near the critical zone. This work focused on using the iTOUGH2 simulator to model the extreme temperature andmore » pressure conditions in magmatic geothermal systems.« less
  • This User's Guide for the LWR Assemblies data base system is part of the Characteristics Data Base being developed under the Waste Systems Data Development Program. The objective of the LWR Assemblies data base is to provide access at the personal computer level to information about fuel assemblies used in light-water reactors. The information available is physical descriptions of intact fuel assemblies and radiological descriptions of spent fuel disassembly hardware. The LWR Assemblies data base is a user-oriented menu driven system. Each menu is instructive about its use. Section 5 of this guide provides a sample session with the datamore » base to assist the user.« less
  • The user's guide gives a complete description of the overall project, including management, research, development, and review procedures that support the integration and improvement of this database, details of the database structure, organization of the CD-ROM, information on use with geographic information systems, as well as information on links between the database and global change characterization and modeling.