TOUGH2, T2VOC, iTOUGH2, TOUGHREACTV1.2, TMVOCV1.0, ECO2N.
Last updated versions are available at http://esd.lbl.gov/research/projects/tough/
|Program name||Package id||Status||Status date|
|Package ID||Orig. computer||Test computer|
|ESTS0219/06||Many Computers||PC Pentium II 450,IBM RISC6000 WS,SUN W.S.,DEC ALPHA W.S.|
|ESTS0219/09||IBM PC,UNIX W.S.||Linux-based PC,PC Windows|
|ESTS0219/10||MacOS,Linux-based PC,PC Windows|
TOUGH2V2.0 is a new and improved version of TOUGH2 for simulating fluid flow and heat transfer in porous media. It is upwardly compatible with and includes all of the capabilities of the earlier version, including the flexibility to handle different fluid mixtures (water, water with tracer; water, CO2; water, air; water, air, with vapor pressure lowering; and water, hydrogen), facilities for processing of geometric data (computational grids), and an internal version control system to ensure referenceability of code applications. Improvements in TOUGH2V2.0 include (1) several new EOS modules for different fluid mixtures, including brines and water-soluble and volatile tracers, (2) enhanced capabilities for previously released fluid property modules, (3) description of diffusion and dispersion in multiphase systems, (4) strongly coupled flow and transport processes, (5) coupling between flow in geothermal reservoirs and wellbores in two-phase conditions, (6) tracer transport with sorption and radioactive decay, (7) flow in media with strong heterogeneity, and (8) a new package of preconditioned conjugate gradient routines for more robust solution of numerically difficult problems. In addition, numerous enhancements were made to facilitate applications to more diverse and demanding flow problems. The TZVOC code for three phase flow of water, air, and a non-aqueous phase liquid (NAPL) has been fully integrated into and is now part of the TOUGH2V2.0 program package. TOUGH2V2.0 comes with a new self-contained users guide that includes technical specifications, a complete reference of input data formats and a collection of sample problems.
TOUGH2V2.0 is a multi-dimensional numerical model for simulating the coupled transport of water, vapor, air, and other fluids, and heat in porous and fractured media. The program provides options for specifying injection or withdrawal of heat and fluids. (more in OTHER PROG/OPER SYS INFO).
iTOUGH2 (inverse TOUGH2) is a computer program that provides inverse modeling capabilities for the TOUGH2 code, a simulator for multiphase, multicomponent, non-isothermal flows in fractured-porous media.
TOUGHREACT is a numerical simulation program for chemically reactive non-isothermal flows of multiphase fluids in porous and fractured media. The program is written in Fortran 77 and was developed by introducing reactive chemistry into the multiphase flow code TOUGH2. Interactions between mineral assemblages and fluids can occur under local equilibrium or kinetic rates. The gas phase can be chemically active. Precipitation and dissolution reactions can change formation porosity and permeability. The code is distributed with a comprehensive user's guide that includes worked sample problems addressing geothermal reservoirs and hydrothermal systems, nuclear waste isolation, groundwater quality, sequestration of carbon dioxide in saline acquifers, and supergene copper enrichment.
Additional information is available on the TOUGH2 homepage, at http://www-esd.lbl.gov/TOUGHREACT/.
TMVOC is a simulator for multicomponent, multiphase, nonisothermal flows of water, soil gas, and several volatile organic chemicals (VOCs). It is based on the M2NOTS-code developed by Adeyinka Adenekan. TMVOC operates within the freamework of TOUGH2, Version 2.0 (T2V2). The code consists of two modules (groups of routines) named *t2fm.f* and *emvoc.f* that must be linked to several standard T2V2 modules.
TMVOC is designed for applications to contamination problems that involve hydrocarbon fuel or organic solvent spills in saturated and unsaturated zones. It can model contaminant behavior under "natural" environmental conditions, as well as for engineered systems, such as siol vapor extraction, groundwater pumping or stem-assisted source remediation.
ECO2N is a fluid property module for the TOUGH2 simulator (Version 2.0) that was designed for applications to geologic sequestration of CO2 in saline aquifers. It includes a comprehensive description of the thermodynamics and thermophysical properties of H2O - NaCl - CO2 mixtures, that reproduces fluid properties largely within experimental error for the temperature, pressure and salinity conditions of interest (10 deg.C <= T <= 110 deg.C; P <= 600 bar; salinity up to full halite saturation). Flow processes can be modeled isothermally or non-isothermally, and phase conditions represented may include a single (aqueous or CO2-rich) phase, as well as two-phase mixtures. Fluid phases may appear or disappear in the course of a simulation, and solid salt may precipitate or dissolve.
TOUGH2V2.0 performs a fully-coupled, simultaneous solution of mass and energy balance for multi-component, multi-phase fluids in permeable media coupled with a multiphase extension of Darcy's Law, and Fickian diffusion and dispersion. The governing equations are discretized using integral finite difference for space and fully-implicit first-order finite differences in time. The resulting non-linear algebraic equations are solved by Newton-Raphson iteration. Iteration is continued until all residuals are reduced below user-specified convergence tolerances. The linear equations arising at each iteration step can be solved with a variety of techniques, including direct matrix methods, and iterative approaches using preconditioned conjugate gradients.
iTOUGH2 solves the inverse problem by automatically calibrating a TOUGH2 model against observed data. Any TOUGH2 input parameter can be estimated based on any observation for which a corresponding TOUGH2 output can be calculated. An objective function measures the difference between the model calculation and the observed data, and a minimization algorithm proposes new parameter sets that iteratively improve the match. Once the best estimate parameter set is identified, iTOUGH2 performs an extensive error analysis which provides statistical information about residuals, estimation uncertainties, and the ability to discriminate among model alternatives. Furthermore, an uncertainty propagation analysis allows one to quantify prediction errors.
iTOUGH2 runs in three application modes:
Uncertainty Propagation Analysis
TMVOC input files are very similar to T2VOC input files. Slight revisions in T2VOC input files are required in data blocks CHEMP for VOC parameters. User options selectable by means of the parameter NKIN are available to initialize a TMVOC run from T2VOC-stuly initial conditions, with automatic internal conversion input data prior to execution.
ECO2N: The partitioning of H2O and CO2 between liquid and gas phases is modeled as a function of temperature, pressure, and salinity, using the recently developed correlations of Spycher and Pruess (2005). Dissolution and precipitation of salt is treated by means of local equilibrium solubility. Associated changes in fluid porosity and permeability may also be modeled. All phases - gas, liquid, solid - may appear or disappear in any grid block during the course of a simulation. Thermodynamic conditions covered include a temperature range from ambient to 100 deg.C (approximately), pressures up to 600 bar, and salinity from zero to fully saturated. These parameter ranges should be adequate for most conditions encountered during disposal of CO2 into deep saline aquifers.
Highly variable depending on problem size and simulation time. Many simple flow problems with a few hundred grid blocks can be run in a matter of minutes on workstations. Large, complex 3-D problems with 50,000 - 100,000 grid blocks in heterogeneous media may run for days.
During the testing at the NEA DB, running time was case dependent, in a PC-Windows/XP Intel Pentium 3.2GHz:
Running Time (sec)
To demonstrate different initialization options
Radial flow from a CO2 injection well
CO2 discharge along a fault zone
CO2 injection into a 2-D layered brine formation
TOUGH2V2.0 has a flexible modular structure that allows modeling of different fluid mixtures. The integral finite difference method used is equivalent to conventional finite differences for regular grid systems, but offers additional flexibility by allowing irregular grids to represent geologic or geometric features. One, two, and three-dimensional systems and single or multiple continuum media are all treated on the same footing. There is an extensive record of code verifications and worked problems.
TOUGH2V2.0 supersedes the previous versions of TOUGH2, which in turn superseded the TOUGH Code. Utility programs to facilitate TOUGH2 applications, such as preprocessing of calculational grids, specification of initial and boundary conditions, and interfacing with graphics programs, are freely available for downloading from the TOUGH2 homepage on the web at http://www-esd.lbl.gov/TOUGH2/.
TOUGHREACT: An intimate knowledge of the TOUGH2V2.0 code is required. It is recommended that TOUGH2V2.0 be licensed along with TOUGHREACT.
TOUGH2 V2.0, T2VOC. TMVOCV1.0 runs interrelated with TOUGH2 V2.0.
TOUGH2, Version 2.0.
CO2TAB3B is an auxiliary program for generating a table of thermophysical properties (CO2TAB file) to be used with TOUGH2/ECO2N.
Updated source code can be found in: http://esd.lbl.gov/research/projects/tough/
|Package ID||Status date||Status|
|ESTS0219/06||20-JUN-2000||Tested at NEADB|
|ESTS0219/09||03-JUN-2010||Tested at NEADB|
Background references for TMVOCV1.0:
K. Pruess, C. Oldenburg and G. Moridis:
TOUGH2 User's guide, Version 2.0, Lawrence Berkeley National Laboratory Report LBNL-43134, Nov 1999.
Numerical Modeling of Multiphase Transport of Multicomponent Organic Contaminants and Heat in the Subsurface, PhD thesis, University of California at Berkeley, Berkeley, CA 94720, 1992.
Adenekan, A.E., T.W. Patzek and K. Pruess:
Modeling of Multiphase Transport of Multicomponent Organic Contaminants and Heat in the Subsurface: Numerical Model Formulation, Water Resour. Res., Vol. 29, No. 11, pp. 3727-3740, 1993.
Background references for TOUGHREACTV1.2:
Xu, T., and K. Pruess:
Modeling multiphase non-isothermal fluid flow and reactive geochemical transport in variably saturated fractured rocks: 1. Methodology, American Journal of Science, 301, 16-33, 2001
A flow problem with up to 800 grid blocks, 2400 connections (interfaces) between them, and 3 equations per grid block requires approximately 4 Megabytes of RAM. More detailed information is available in the user's guide.
Tested at the NEA DB on PC-Windows/XP using Compaq Visual Fortran 6.6A and PC-linux/Fedora12 using Intel Fortran Compiler V10.1.
No executables are included in the package.
|Package ID||Computer language|
Operating System-Independent (Requires 64-bit compiler to provide double precision arithmetic).
ECO2N is 'plug-compatible' with TOUGH2, Version 2.0. The code is intrinsically single-precision, but requires 64-bit arithmetic, which on 32-bit processors (PCs and workstations) may be realized by using compiler options for generating 64-bit arithmetic.
ECO2N may be linked with standard TOUGH2 modules like any of the fluid property modules included in the TOUGH2 V 2.0 package.
Although primarily designed for geothermal reservoir studies, and isolation of high level nuclear waste in partially saturated geological media, TOUGH2V2.0 should also be useful for a wider range of problems in heat and moisture transfer, and in the drying of porous materials. The simulator was primarily developed for problems involving strongly heat driven flow. To describe these phenomena a multiphase approach to fluid and heat flow is used, which fully accounts for the movement of gaseous and liquid phases, their transport of latent and sensible heat, and phase transitions between liquid and vapor. TOUGH2V2.0 takes account of fluid flow in both liquid and gaseous phases occurring under pressure, viscous, gravity forces according to Darcy's Law. Interference between the phases is represented by means of relative permeability functions. The code handles multiphase diffusion effects coupled with phase partitioning, and represents capillary and phase adsorption effects for the liquid phase. Heat transport occurs by convection and diffusion, accounting for sensible and latent heat, and by conduction with thermal conductivity dependent on water saturation. TOUGH2V2.0 can also efficiently model isothermal flow problems, such as unsaturated flow described by Richard's equation, or water and tracer transport in groundwater aquifers.
TOUGH2V2.0 includes the proprietary software package MA28 from the Harwell Subroutine Library.
TMVOC - Storage requirements increase with number of grid blocks (mesh size). Maximum problem dimensions can be specified through PARAMETER statements in the main program, or through an INCLUDE file.
While testing the code at the NEA DB, a bug was found while testing ECO2N in TEST4. For high pressures (P & P(NPK), with P(NPK) is the maximum value in CO2TAB file, the EWASG- extrapolation is used. However, this extrapolation was meant for the "low-pressure end". It was never intended for high pressures. The code stops execution. This bug was reported to the author (K. Pruess) and the author confirms and suggested a modification in the program to be posted in the weblink "bugs and fixes"(http://esd.lbl.gov/research/projects/tough/).
In order to avoid this EWASG- extrapolation users could extent the P(NPK) in CO2TAB to higher values (this is the solution adopted by NEA in TEST4).
S. Finsterle, K. Pruess, C. Oldenburg, G. Moridis, Tianfu Xu, E. Sonnenthal, N. Spycher
Lawrence Berkeley National Lab.,
Mail Stop 90-1116
Berkeley, CA 94720
Keywords: environment, equations of state, fluid flow, gaseous wastes, geologic structures, geothermal, heat flow, hydrology, multiphase flow, porous materials, radioactive waste storage, two-phase flow, water pollution.